|
|
||||||||




* Beth Israel Deaconess Medical Center, Harvard Medical School, Boston, MA 02215;
Oregon National Primate Research Center, Oregon Health & Science University, Beaverton, OR 97006; and
Dana-Farber Cancer Institute, Harvard Medical School, Boston, MA 02115
| Abstract |
|---|
|
|
|---|
T and G
A replacements together accounting for almost one-third of all mutations. Analysis of replacement and silent mutations in the framework and CDRs suggests that the Igs were subjected to affinity selection. These data demonstrate that the process of Ab maturation is not seriously disrupted in GCs during the early stages of immunodeficiency virus infection, and that Env-specific Igs developing in GCs are subject to extensive somatic mutation and profound selection pressures. | Introduction |
|---|
|
|
|---|
Critical processes occur in lymphoid follicles that lead to the generation of mature Ab responses. An iterative process of selection, expansion, mutation, and further selection culminates in the generation of high-affinity memory B lymphocytes and plasma cells. Specific features of the somatic mutations observed in expressed Igs have been well-characterized (3, 4) and are similar in humans and mice. Some disease states are associated with atypical patterns of somatic mutation in expressed Igs (5, 6, 7, 8, 9). In light of the follicular abnormalities seen in the setting of HIV/SIV infections, the processes that transpire in GCs may also show abnormalities in HIV/SIV-infected individuals.
The Ab response to the envelope glycoproteins (Env) of HIV-1 is of particular importance because of the ability of some anti-Env Abs to neutralize viral infectivity, protect against infection (10, 11, 12, 13, 14), exert selective pressure on HIV-1 and SIV following infection in vivo (15, 16, 17, 18, 19, 20, 21, 22, 23), and contribute to the control of viremia (24, 25, 26). Importantly, neutralizing Abs effective against a broad range of HIV-1 or SIV isolates are produced upon virus exposure only after a substantial delay, months or years following infection (16, 27, 28, 29, 30, 31, 32), suggesting that maturation of the immune response is necessary for the development of such broadly reactive Abs. Changes in the specificity of circulating neutralizing Abs continue to occur throughout at least the first 5 years of HIV-1 infection (33). The molecular mechanisms that underlie this evolution in virus-specific Abs have not been identified, but it has been suggested that GC processes may be critical for development and maintenance of virus-neutralizing Abs (34).
Macaques infected with recombinant viruses that express HIV-1 Env on an SIV backbone (simian HIVs (SHIVs)) have proven valuable for studying humoral immunity to HIV-1 (11, 12, 24, 35). A high degree of conservation of Ig sequences between humans and these monkeys increases their value as animal models for studying human Ab responses. Studies have shown that homologous human and macaque H (36, 37, 38, 39, 40) and L chain (41) V regions typically have >90% nucleotide identity, and even >95% identity for some V region genes.
In previous studies of SHIV-infected rhesus monkeys, we have shown that lymphoid GCs formed during primary and chronic immunodeficiency virus infection include a high frequency of Env-specific B cells (42). In the present report, we extend those studies with an analysis of Ig sequences from Env-binding GCs.
| Materials and Methods |
|---|
|
|
|---|
Recombinant HIV-1 89.6 gp120 (rgp120-89.6) was prepared by Dr. R. Doms (University of Pennsylvania, Philadelphia, PA) as previously described (43, 44) from supernatants of cells infected with recombinant vaccinia viruses expressing HIV envelope glycoproteins, and purified by lectin affinity chromatography. The rgp120-89.6 was biotinylated in 100-µg batches using the FluoReporter Minibiotin-XX protein labeling kit (Molecular Probes) according to the manufacturers instructions. Unreacted biotin was removed by filtration over a ChromaSpin LC-30 column (BD Clontech). Preservation of conformation-dependent epitopes in the biotinylated Env was documented as described (42).
SHIV clones, animal infections, and sample handling
The SHIVs used for this study are described elsewhere (45). In brief, the env, tat, and rev genes of HIV-1 HXBc2 were introduced into an infectious molecular clone of SIVmac239. The env gene of the resulting virus was replaced with the env gene of the HIV-1 clone 89.6 to generate SHIV-89.6 (46). Serial passage of SHIV-89.6 in monkeys generated the highly pathogenic SHIV-89.6P quasispecies, from which SHIV-KB9 was cloned. SHIV-KB9ct and SHIV-KB9(305) are hybrid viruses that combine aspects of the parental 89.6 and more pathogenic KB9 clones. These viruses were selected for the present study because they expressed a gp120 identical with or cross-reactive with that of the rgp120-89.6 used as an inverse immunohistochemical probe and show increased pathogenicity compared with the parental SHIV-89.6 (45, 47).
The rhesus monkeys used in this study were maintained at the Oregon National Primate Research Center in accordance with the guidelines of the committee on Animals for the Harvard Medical School and the Guide for the Care and Use of Laboratory Animals (48). Healthy rhesus monkeys were inoculated i.v. with cell-free SHIV. Subject Mm18284 was infected with SHIV-KB9ct, and subject Mm19941 was infected with SHIV-KB9(305). Some immunologic and virologic data from these monkeys have been previously reported (42, 45, 47). Infection of both animals was confirmed by serologic and virologic assays (45, 47). Before infection, subject Mm18284 had a CD4+ cell count of 1200/µl. Following infection the CD4+ cell count fell to a nadir of 350 on postinfection day (p.i.d.) 17 before stabilizing at
500 by p.i.d. 28. Plasma viral load as measured by ELISA detection of the SHIV Gag protein p27 peaked at 1.4 ng/ml on p.i.d. 14, and declined to below the limit of detection by p.i.d. 21 (47). Mm19941 had a baseline CD4+ cell count of 2000/µl. Following infection the CD4+ count fell to 831 by p.i.d. 21 and rose gradually to stabilize at
950 on p.i.d. 78. Cumulative p27 antigenemia in this animal through p.i.d. 21 was 0.95 ng/ml; viral p27 was below the limit of detection beyond that date (45). The sustained >50% decrease in absolute CD4 counts in the blood of both animals following infection is evidence of the pathogenicity of the SHIV types used in these studies
Monkeys were anesthetized with ketamine-HCl 5 mg/kg, i.m. for all blood sampling, and were euthanized by an overdose of sodium pentobarbital administered i.v. at 3 mo (Mm19941) and 9 mo (Mm18284) after SHIV infection. These time point were chosen to assess GC functions during the early stages of chronic infection, at a time when Env-specific GCs are plentiful (42). Pieces of spleen were frozen in Tissue-Tek cryoembedding medium (Fisher Scientific) immediately upon removal, and stored at 80°C until use. Frozen sections were cut at 10-µm thickness, thaw-mounted onto Superfrost Plus slides (Fisher Scientific), immediately refrozen, and then stored at 80°C.
Inverse immunohistochemistry
Slide-mounted tissue sections of spleen from SHIV-infected monkeys were fixed in cold acetone for 10 min, air dried, and rehydrated in TBS (pH 8). Slides were then processed using the Sequenza immunostaining workstation (Shandon Lipshaw). Inverse immunostaining was enhanced using catalyzed signal amplification with the Renaissance TSA-Indirect kit (PerkinElmer). Each slide was first treated with Tris-NaCl buffer with a proprietary blocking reagent (Renaissance kit) to reduce nonspecific staining, and then incubated for 1 h at room temperature with biotinylated rHIV-1 envelope gp120-89.6 at a concentration of 0.2 mcg/ml (1.67 nM). Slides were rinsed three times with TBS containing 0.05% Tween 20 (TBST) after each subsequent staining step. To catalyze amplification of the biotin signal, the slides were incubated with streptavidin-HRP (1/100 dilution, Renaissance kit) for 30 min at room temperature and subsequently incubated with biotinyl tyramide (1/50 dilution) for 10 min. The amplified biotin signal, indicating bound biotinylated rgp120, was detected histochemically using alkaline phosphatase-conjugated streptavidin-biotin complex (DakoCytomation) and Vector Blue chromogenic substrate (Vector Laboratories) together with levamisole (DakoCytomation) to suppress endogenous alkaline phosphatase activity. After counterstaining with either nuclear fast red (Vector Laboratories) or Gills hematoxylin (Fisher Scientific), slides were covered with Crystal/Mount permanent mounting medium (Fisher Scientific) and air dried.
Microdissection of Env-binding GCs
Slide-mounted splenic frozen sections stored at 80°C were fixed in cold acetone for 10 min and air dried. Sections were rehydrated in nuclease-free water, stained in nuclear fast red (Vector Laboratories), dehydrated through graded alcohols and xylenes, and air dried. The targeted GCs were microdissected using an Arcturus PixCell II laser microdissecting workstation (Arcturus Engineering) and the following method. After Env-binding GCs were identified by inverse immunohistochemical staining, the same GCs were located on adjacent tissue sections using landmarks such as position within the section, orientation relative to other GCs, and distinctive features such as GC contour or the presence of vessels within the structure. Cells from each targeted GC were isolated from the surrounding tissue using laser capture microdissection (49) with CapSure HS membranes (Arcturus Engineering). Multiple overlapping laser pulses of 30-µm diameter spot size were used to capture substantial portions of each GC from the slide. Care was taken to avoid capturing cells from the surrounding mantle zones and other regions outside the targeted zone. Digitized images of immunostained tissue sections and the microdissection process were prepared using the Arcturus Engineering workstation and formatted using Adobe Photoshop version 4 software.
Cloning and sequence analysis of Ig genes from Env-binding GCs
mRNA present in the captured cells was purified using the RNeasy kit (Qiagen) according to the manufacturers instructions. cDNA was prepared using the AMV Reverse Transcription System (Promega). The resulting cDNA (1 µl/reaction) was amplified using nested or seminested PCR catalyzed by a proofreading enzyme, PfuTurbo DNA Polymerase (Stratagene) to provide high copying fidelity and optimal yield from the minute samples obtained using microdissection. Each round of the PCR was performed in 50 µl using oligonucleotide primers designed to selectively amplify rearranged Ig H or L chains (Table I) for 30 cycles using the following parameters: 20 s at 94°C, 30 s at 55°C, 120 s at 72°C. Each PCR was concluded with an extension step at 72°C for 10 min. Products of the second-round PCRs were size selected by electrophoresis in 1.5% agarose, visualized by ethidium bromide staining and cut from the gel. The PCR products were purified using the Qiaex II Gel Extraction System (Qiagen), and cloned into Escherichia coli using the TopoBlunt cloning kit (Invitrogen Life Technologies). Plasmids containing inserts were isolated from expanded cultures using QIAprep Spin Miniprep kit (Qiagen), and sequenced on both strands using an Applied Biosystems Prism 377XL Automated DNA Sequencer and primers that bound outside the insert.
|
Somatic mutations in expressed Ig V regions are usually identified with reference to the sequence of their probable precursor in the germline. Because the macaque Ig loci have not been comprehensively catalogued, macaque Ig cDNA sequences were compared with their nearest homologs among the human germline variable (V), diversity (D), and joining (J) region minigenes, which are closely related to those of macaques (36, 37, 38, 39, 41). The nearest homolog minigenes were identified from the complete dataset of human germline Ig sequences based on sequence similarity (50). To facilitate analysis of the CDR3 regions, the germline segments were trimmed and joined to create hypothetical junctional sequences; these junctions were optimally assembled to minimize the number of nucleotide differences with the macaque clones. When V, D, and J segments were insufficient to generate a CDR3 of the observed length, N and P nucleotides were not added; these positions are shown as blanks in the alignment.
Somatic mutations were tallied twice, under different assumptions. First, any nucleotide difference when compared with the nearest human germline V region and hypothetical junctional sequence was assumed to result from a somatic mutation (assumption 1). Second, it was assumed that any nucleotides present in every member of a lineage were encoded in the macaque germline, so that only intraclonal sequence differences were counted as somatic mutations (assumption 2). Statistical software packages SAS version 8.2 and Splus version 3.4 were used to read in the nucleotide sequence data and determine the number of occurrences for each possible type of substitution (e.g., from nucleotide X to nucleotide Y) for each clone under both assumptions 1 and 2.
Because the genealogical relationships among the clones were inferred but not known with certainty, statistical analyses were based on a population model: the incidence of mutation at each nucleotide position was tallied cumulatively for each occurrence, even when the same substitution occurred in multiple members of the population.
Nucleotide differences between the clones were tallied for the entire V region excluding positions that overlap the PCR primers, which correspond with the initial 22 nt of the IgH framework region 1 (FR1) region or 18 nt of the Ig
and Ig
FR1. FR and CDR in the cloned Igs were identified based on standard Kabat criteria (E. A. Kabat, Kabat database of sequences of proteins of immunological interest, National Institutes of Health,
www.kabatdatabase.com/
), summarized online at
www.rubic.rdg.ac.uk/abeng
.
Nucleotide mutations that result in amino acid substitutions (replacement mutations, R) or do not change the amino acid sequence (silent mutations, S) were counted separately for FR1, 2, and 3 and CDR1, 2, and 3 of each macaque clone. Enumeration of R and S mutations was performed using an automated algorithm that incorporates the methods of Nei (51, 52), implemented online at
http://hiv-web.lanl.gov/content/hiv-db/SNAP/WEBSNAP/SNAP.html
(B. T. Korber, Synonymous Non-Synonymous Analysis Program, Los Alamos National Laboratory, 53). The excess or paucity of replacement mutations in the FRs or CDRs of each sequence was compared with that expected by chance using a multinomial distribution model, implemented online at
www-stat.stanford.edu/Ig
(54).
For analysis of purine or pyrimidine transition or transversion mutations, the test for whether the rates of mutation from T to C, T to A, and T to G were all equal given that T was mutated (and similarly for mutations from C, from A, and from G) was an exact Poisson test (i.e., the multinomial test of ratios after conditioning on the sum of the three mutations). Comparisons between pairs of mutation rates used the exact Poisson test (i.e., the binomial test of ratios after conditioning on the sum). The tests of differences between lineages within chain (H or L) and animal used the exact Zelen test for homogeneity of odds ratios. An adjustment was made for multiple comparisons using Holms method. That is, when there were n relevant comparisons, the smallest of the n p values was compared with 0.05/n, the next smallest p value was compared with 0.05/(n-1), and so on until one of the p values was not smaller than its bound (i.e., 0.05/n), in which case that comparison and all comparisons with larger p values were declared nonsignificant. The multiple comparison adjustment was done separately for each different initial nucleotide.
Within the motif RGYW and its reverse complement WRCY, only the underlined G or C was considered a potential hot spot preferentially targeted for mutation (55, 56). Comparison of mutation frequency for hot spot and non-hot spot sites was made using the Fisher exact text.
| Results |
|---|
|
|
|---|
Splenic tissue frozen sections derived from two rhesus monkeys chronically infected with moderately pathogenic variants of SHIV-89.6 were stained by inverse immunohistochemistry using biotinylated recombinant HIV-1 envelope gp120 (Env) as an Ag probe. Splenic tissues from both monkeys demonstrated numerous GCs containing cells that bound Env. For Mm19941, 35% of >1500 examined GCs contained Env-binding cells, while for Mm18284, 3% of >400 evaluated GCs contained such cells. Both intensely stained and entirely unstained GCs were frequently apparent on a single tissue section in close proximity (data not shown), indicating that the Ag probe was bound to B cells or Abs produced locally rather than to serum Abs nonspecifically bound to GC follicular dendritic cells. One Env-binding GC from each animal, shown in Fig. 1, A and B, was selected for further evaluation.
|
Samples from each tissue section were processed separately. mRNA present in the captured cells was purified, reverse transcribed into cDNA, and amplified using nested PCR with oligonucleotide primers designed to selectively amplify rearranged Ig H or L chains (Table I). Products of the second-round PCRs were size selected, cloned, and both strands were sequenced.
Sequence analysis of Igs cloned from microdissected GCs
The recovered nucleotide sequences from each GC could be grouped on the basis of sequence similarity into multimember clonal lineages (Table II). Within each lineage, all sequences have similar V-D-J or V-J junctions but differ by somatic mutations distributed throughout the variable domain. The V-D-J-C and V-J-C joins were all "in frame" (i.e., productive rearrangements), and no premature termination codons were found in any of the clones. Data from nine clonal lineages were combined to create an aggregate data set that was analyzed statistically.
|
H chain. cDNA from GC2 was amplified using oligonucleotide primers designed to amplify Ig H chain V regions encoding IgG or IgM Abs. The IgG-primed reactions revealed an amplicon of the expected size only for the VH3-primed reaction. Nucleic acid sequencing of VH3-IgG clones yielded a clonal lineage, designated 2S1H3, with nine distinct members (Fig. 2A).
|
L chain.
cDNA from GC 2 was amplified using oligonucleotide primers designed to amplify Ig L chain V regions of either
or
isotype. Electrophoresis of PCR products revealed an amplicon of the expected size only among the VL5-primed reaction. Nucleic acid sequencing of clones from GC2 yielded one multimember clonal lineage, designated 2S2L5, with six distinct members (Fig. 2B).
Ig clones from GC4
H chains. cDNA from GC4 was amplified using oligonucleotide primers designed to amplify Ig H chain V regions encoding IgA, IgD, IgE, IgG, or IgM Abs. Only IgG-specific PCRs produced an amplicon of the expected size. Sequencing of clones from the IgG-specific PCR revealed three multimember clonal IgG lineages, designated 4S1H3, 4S2H3, and 4S3H3, whose sequences are presented in Fig. 2, CE, respectively.
L chains.
cDNA from GC4 was amplified using oligonucleotide primers designed to amplify Ig L chain V regions of either
or
isotype.
L chain PCRs yielded Ig sequences that belonged to three distinct multimember clonal lineages, designated 4S4L5, 4S5L2, and 4S6L8, whose sequences are presented in Fig. 2, FH.
L chain PCRs yielded Ig sequences that could be assigned to a single multimember clonal lineage.4S7K2, presented in Fig. 2I.
Frequency and rate of somatic mutations in V regions of the Ig clones
The somatic mutation rate was calculated based upon two different assumptions, as described in Materials and Methods. Under assumption 1, the 66 Ig sequences cloned from microdissected GCs of SHIV-infected monkeys have an average of 37.1 apparent somatic mutations per sequence, with a range of 24.356.8. This corresponds to an average mutation frequency of 0.11 mutations/nucleotide. No significant difference was found between H and L chain clones in this regard. Under assumption 2, the Ig sequences have an average of 16.5 apparent somatic mutations per sequence, with a range of 2.037.9. The corresponding mutation frequency is 0.05 mutations/nucleotide. Under assumption 2, the mutation frequency of H chains (0.015) was significantly less than that of L chains (0.07; p < 0.02).
Bias in somatic mutation of Ig V regions
Some disease states are associated with atypical patterns of somatic mutation in expressed Igs (5, 6, 7, 8). We therefore analyzed the pattern of mutations observed in GC B cells from these SHIV-infected monkeys to determine whether they are indicative of abnormalities.
The frequency with which mutations were targeted to nucleotides within the previously described RGYW/WRCY "hot spot" motif was statistically analyzed for the 66 Ig clones obtained from the SHIV-infected monkey spleens. Within and across animals, for both H and L chains, there was a higher probability of mutation at hot spots than elsewhere (p < 0.0001). Preferential targeting of hot spots was also evident within each lineage (Table III), despite the smaller sample size. Overall, the percentage of hot spot sites that were mutated was more than twice as large as that seen in non-hot spots.
|
Under assumption 1, there were striking deviations from randomness in the mutational pattern, which was similar to that of normal human B cells in showing a disproportionate (increased) frequency of transition mutations. For example, tests were done evaluating the relative likelihood of pairs of specific mutations (e.g., from T to C vs T to A). Such comparisons were only performed for lineages where, given that a nucleotide of a particular type (T, C, A, or G) mutates, it was not equally likely to mutate to each of the other three nucleotides. After pairwise comparisons, it was found that there were 28 (of a possible 93) significant differences, and 20 of these favored transition mutations. The other eight significant differences occurred within a nucleotide group (e.g., from T to A vs T to G). In all comparisons between transition and transversion mutations showing a statistically significant difference, transition mutations were favored. Overall, 58% of mutations were transitions (Table IV), similar to the normal human pattern but contrasting with the stochastic prediction of 33% transitions.
|
Mutational bias attributable to cytosine deaminase (57) was specifically examined. Under assumption 1, for eight of the nine lineages, the nucleotide C when mutated was most often replaced by T. Mutational bias in mutations of C was statistically significant for seven of nine lineages (range of p values, 0.0008-<0.0001). Similarly, eight of nine lineages showed preferential replacement of G by A, and the deviation from randomness was statistically significant for all of these (range of p values, 0.021-<0.0001). Combining the results for all lineages, the C
T and G
A replacements together constituted almost one-third of all mutations (Table IV). Under assumption 2 also, the nucleotide C when mutated was most often replaced by T, and G was most often replaced by A. Thus, calculations based on either assumption support the conclusion that somatic mutational processes are operating normally in these GCs.
Selection pressures on somatic mutations in Ig V regions
To assess the influence of positive and negative selection on the Ig sequences cloned from the Env-binding GCs of SHIV-infected monkeys, the ratio of replacement to silent mutations was determined for the FRs and CDRs of each V(D)J region by comparison with the nearest human germline homologs (Table II). The excess or paucity of replacement mutations in the FRs or CDRs of each sequence, tallied both under assumptions 1 and 2, was compared with that expected by chance using a multinomial distribution model (54). The fraction of clones from each lineage that showed statistically significant evidence for a selection bias under assumption 1 is shown in Table II. Remarkably, 63 of 66 clones had a statistically significant (p < 0.05) paucity of R mutations in the FRs. This indicates that nearly all the clones had survived substantial negative selection pressures in vivo. The absence of premature termination codons, out-of-frame V(D)J junctions and other inactivating mutations among the cloned sequences provides additional evidence of negative selection. A majority of the clones (37 of 66) also had a statistically significant excess of R mutations in the CDRs, reflecting positive selection pressures in vivo.
We repeated the analysis under assumption 2 (data not shown). Although this assumption drastically reduced the total number of apparent somatic mutations, 35 of 66 clones still showed a statistically significant paucity of R mutations in their FRs, and 27 clones had a statistically significant excess of R mutations in the CDRs. Therefore, calculations based on either assumption support the conclusion that positive and negative selection processes are operating normally in these GCs.
| Discussion |
|---|
|
|
|---|
The abundance of mutations in the Igs indicates that the B cells were subjected to a profound degree of antigenic stimulation in the GC microenvironment. Although it was not formally demonstrated that the studied Igs actually bind Env or another SHIV protein, the abundance of Env-binding cells in the studied GCs provides direct evidence that the viral envelope is a potent immunogen for B cells in fully functional GCs.
Somatic mutation among GC B cells typically occurs at a rate of 103 mutations per base pair per cell per generation (66, 67). The cell cycle of GC B cells has been calculated at between 6 (68) and 12 h (69). Assuming the same rate of mutation applies to the tissues from which the clones were derived, calculations show that the average clone in our dataset underwent between 50 and 113 cycles of replication. Based on the 612 hour cycling time, this number of generations could arise in 28 wk of continuous antigenic stimulation. Because the tissues were obtained at least 13 wk after SHIV infection, it cannot be determined whether diversification of the clones we studied occurred rapidly during acute infection or accumulated more gradually over subsequent weeks.
Using the lower estimate of 50 generations, this number of cell divisions is sufficient to produce nearly 800 trillion cells, if all survived. However, the largest GCs contain a few million cells at most. The difference between these values highlights the stringency of negative selection enforced during clonal expansion in these GCs.
In the present study, we found that GC function was preserved in two monkeys in the steady state stage of SHIV infection who had
50% reductions in CD4+ T cell counts. This corresponds with the protracted lymphadenopathy stage of HIV-1 infection in the early years following infection, although the subset of CD4+ T lymphocytes lost in these monkeys differs from those targeted for destruction in HIV-1-infected humans (70). However, GC function likely declines with more advanced immunodeficiency. HIV-1-infected individuals with very low numbers of CD4+ cells often develop follicular involution, with GCs depleted of lymphocytes and follicular dendritic cells.
The analysis of data in this study was subject to several limitations. First, any characterization of GC B cell lineages will necessarily be incomplete because the iterative process of mutation and selection will delete intermediate ancestors, branches, and entire lineages as an indispensable component of the process leading to affinity maturation. Other B cells present in the GC will be missed through sampling error or because mutations in their Ig sequence disrupt homology with the PCR primers. In any case, because the H and L chain lineages characterized in this study all showed a consistent pattern of nucleotide substitution and selection, it is unlikely that any missing sequences would have altered our conclusions.
A more serious limitation in any analysis of macaque Ig sequence data stems from the substantial gaps in our knowledge of macaque Ig germline sequences. These gaps preclude the use of macaque Ig germline sequences as a basis for comparison with the Ig clones sequenced in this study. Instead, the sequenced macaque Ig cDNAs were compared with the human genomic Ig repertoire. This germline repertoire is known to be very similar but not identical with that of the macaques (36, 37, 38, 39, 41). As a result, where a difference from the human germline sequence was shared by all clones in a macaque lineage, we could not determine whether the differences were introduced by somatic mutation in an ancestral clone or instead were encoded in the macaque germline.
To rule out possible bias introduced by such errors of assignment, we performed most analyses under two different assumptions. First, it was assumed that any nucleotide difference between a macaque sequence and its nearest human germline homolog resulted from a single somatic mutation. This assumption will tend to overestimate the frequency of somatic mutations. Second, it was assumed that any substitutions present in every member of a lineage were encoded in the macaque germline, and were therefore excluded from analysis. This more conservative assumption will undercount somatic mutations, including those retained because they contribute to the specificity of all Igs of a lineage. Importantly, the conclusion that somatic mutation and selection were operating as they should be in the GCs of the SHIV-infected monkeys was supported under either assumption. There also is little published information regarding the normal pattern of somatic mutation in macaques. The lack of a formal baseline or comparator would be of greater consequence had our results suggested that macaque GCs introduce somatic mutations according to some novel pattern. However, because the pattern was the same in normal humans and rodents as in the SHIV-infected macaques we studied, it is likely that the pattern in normal macaques is also the same.
There are few published studies that address the role of GCs during immune responses to viral infection in humans or nonhuman primates. Analyses of somatic mutations in mAbs derived from individuals exposed to viruses suggest that the processes of somatic mutation and Ag-affinity selection can play a major role in shaping the Ab response to viruses, including HIV-1 (58, 59, 71, 72, 73). The present study directly demonstrates that GCs are a major site for somatic mutation and selection of antiviral Abs in primates. Additional studies are needed to determine the contribution of somatic mutation outside of GCs (74, 75, 76) to shaping the Ab response to an infecting virus.
| Acknowledgments |
|---|
| Disclosures |
|---|
|
|
|---|
|
|
| Footnotes |
|---|
1 D.H.M. was supported by Public Health Service (PHS) Grants AI-01587 and AI-45370 from the National Institute of Allergy and Infectious Diseases (NIAID). E.H.S. participated in these studies as a Medical Student Fellow of the Howard Hughes Medical Institute. M.K.A. was supported by Grant AI-054292 from the NIAID, Grant CA-75922 from the National Cancer Institute, Grants RR-016025 and RR018107 from the National Center for Research Resources, and a PHS institutional grant from the National Center for Research Resources supporting the Oregon National Primate Research Center, Grant RR-000163. N.L.L. was supported by PHS Grant AI-067854 from the NIAID Center for HIV/AIDS Vaccine Immunology and Grant AI-20729 from the NIAID. O.G.G. and R.S.G. were supported by PHS Grant AI-28691 to the Dana-Farber Cancer Institute Center for AIDS Research from the NIAID. ![]()
2 Current address: Massachusetts General Hospital, Wang Ambulatory Care Center-835, Boston, MA 02114. ![]()
3 Current address: University of Michigan Hospital Psychiatry 9d9812 University Hospital 0120, 1500 East Medical Center Drive, Ann Arbor, MI 48109-0118. ![]()
4 Current address: Department of Surgery, Mount Sinai School of Medicine, 5 East 98th Street, New York, NY 10029. ![]()
5 Address correspondence and reprint requests to Dr. Norman L. Letvin, Viral Pathogenesis Division, RE-113/Beth Israel Deaconess Medical Center, 330 Brookline Avenue, Boston, MA 02215. E-mail address: nletvin{at}bidmc.harvard.edu ![]()
6 Abbreviations used in this paper: GC, germinal center; SHIV, simian HIV; p.i.d., postinfection day; FR, framework. ![]()
Received for publication July 18, 2005. Accepted for publication April 19, 2006.
| References |
|---|
|
|
|---|
genes in systemic lupus erythematosus. Arthritis Res. 3: 368-374. [Medline]
locus of primates. Genomics 25: 623-629. [Medline]
-deficient mice. Nature 382: 462-466. [Medline]This article has been cited by other articles:
![]() |
U. Hershberg, M. Uduman, M. J. Shlomchik, and S. H. Kleinstein Improved methods for detecting selection by mutation analysis of Ig V region sequences Int. Immunol., May 1, 2008; 20(5): 683 - 694. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |