Previous Article | Next Article ![]()
Journal of Virology, September 2008, p. 9216-9227, Vol. 82, No. 18
0022-538X/08/$08.00+0 doi:10.1128/JVI.01041-08
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

*
Chanson J. Brumme,1,
Jonathan Carlson,2,3,
Hendrik Streeck,1
Mina John,4
Quentin Eichbaum,1
Brian L. Block,1
Brett Baker,1
Carl Kadie,2
Martin Markowitz,5
Heiko Jessen,6
Anthony D. Kelleher,7
Eric Rosenberg,1
John Kaldor,7
Yuko Yuki,8
Mary Carrington,8
Todd M. Allen,1
Simon Mallal,4
Marcus Altfeld,1
David Heckerman,2 and
Bruce D. Walker1,9
Partners AIDS Research Center, Massachusetts General Hospital, Harvard Medical School, Boston, Massachusetts,1 Microsoft Research, Redmond, Washington,2 Department of Computer Science, University of Washington, Seattle, Washington,3 Centre for Clinical Immunology and Biomedical Statistics, Royal Perth Hospital and Murdoch University, Perth, Australia,4 Aaron Diamond AIDS Research Center, New York, New York,5 Jessen Praxis, Berlin, Germany,6 National Centre in HIV Epidemiology and Clinical Research, University of New South Wales, Sydney, Australia,7 Cancer and Inflammation Program, Laboratory of Experimental Immunology, SAIC-Frederick, Inc., NCI-Frederick, Frederick, Maryland,8 Howard Hughes Medical Institute, Chevy Chase, Maryland9
Received 18 May 2008/ Accepted 1 July 2008
|
|
|---|
80% of published CTL epitopes. Five of the 10 most rapidly evolving epitopes were restricted by protective HLA alleles (HLA-B*13/B*51/B*57/B*5801; P = 0.01), supporting the importance of a strong early CTL response in HIV control. Consistent with known fitness costs of escape, B*57-associated mutations in Gag were among the most rapidly reverting positions upon transmission to non-B*57-expressing individuals, whereas many other HLA-associated polymorphisms displayed slow or negligible reversion. Overall, an estimated minimum of 30% of observed substitutions in Gag/Pol and 60% in Nef were attributable to HLA-associated escape and reversion events. Results underscore the dominant role of immune pressures in driving early within-host HIV evolution. Dramatic differences in escape and reversion rates across codons, genes, and HLA restrictions are observed, highlighting the complexity of viral adaptation to the host immune response. |
|
|---|
It is now understood that CTL escape occurs along generally predictable pathways based on the host HLA profile (2, 17, 27), and that escape and reversion represent major forces driving viral evolution and diversity at both the individual and population levels (2, 11, 14, 34, 39, 52, 57-59, 65, 67, 68, 70). However, due in part to the lack of a consistent definition of HLA-associated mutation, as well as the lack of large longitudinal HIV sequence datasets (2, 9, 46, 61), HLA-driven viral adaptation in early HIV infection remains incompletely characterized. In this study, we employ a strict definition of HLA-attributable substitution based on a comprehensive, predefined list of HLA-associated polymorphisms in HIV type 1 (HIV-1) subtype B in order to estimate the proportion of viral evolution attributable to HLA-associated selection pressures in Gag, Pol, and Nef in the first year of HIV infection in a cohort of 98 untreated, subtype B-infected seroconverters. In addition, we systematically compute the rates of escape and estimate the rates of the reversion of these HLA-associated polymorphisms in early HIV infection, revealing marked differences in the rates of escape and reversion across codons, genes, and HLA restrictions.
|
|
|---|
The date of HIV infection was estimated using clinical history (where available) by subtracting 4 weeks from the baseline (for stages 1a or 2a), by subtracting 6 weeks from the baseline (for stages 1b or 2b), by calculating the midpoint between the last negative and the first positive EIA (for stage 3a), or by subtracting 4 months from the baseline (for stage 3b). Subjects were monitored for a median of 425 days after the estimated infection date. The study was approved by the respective institutional review boards and was conducted in accordance with human experimentation guidelines of Massachusetts General Hospital. All subjects provided written informed consent.
HIV RNA genotyping and HLA typing. HIV-1 RNA genotyping of gag, pol, and nef was performed on serial plasma samples (median of 5 samples/patient). In general, patients were monitored at baseline, 1 month, and every 2 or 3 months thereafter, although in some cases more frequent sampling was performed. HIV-1 gag (codons 1 to 500; HXB2 nucleotides [nt] 790 to 2289), a portion of pol (protease codons 1 to 99 and reverse transcriptase [RT] codons 1 to 400; nt 2253 to 3749), and nef (codons 1 to 206; nt 8797 to 9414) were amplified by nested reverse transcription-PCR from extracted plasma HIV RNA using gene-specific primers. Population (bulk) sequencing was performed on an Applied Biosystems 3730 automated DNA sequencer. Data were analyzed using Sequencher (Genecodes). Nucleotide mixtures were assigned if the secondary peak height exceeded 25% of the dominant peak height. Data were aligned to the HIV-1 subtype B reference HXB2 sequence (GenBank accession no. K03455) using a modified NAP algorithm (42). All subjects harbored subtype B infections, as confirmed by comparing gag/pol/nef sequences to all HIV subtype reference sequences (http://hiv-web.lanl.gov/content/hiv-db/mainpage.html). HLA class I typing was performed by standard sequence-specific PCR or sequence-based typing.
Predefined list of HLA-associated polymorphisms.
Previous studies of immune-driven HIV evolution have lacked a consistent definition of HLA-attributable mutation. We address this by defining HLA-associated escape and reversion based on a predefined, comprehensive list of HLA-associated polymorphisms identified in a cross-sectional analysis of >1,200 chronically infected, antiretroviral-naïve individuals from three published cohort studies in Canada (17, 18), the United States (45), and western Australia (11) using phylogenetically corrected methods and a q-value correction for multiple tests (q < 0.2, which translates to a 20% false discovery rate) (60, 75). As described previously (17; D. Heckerman, C. J. Brumme, M. John, J. Carlson, R. Haubrich, S. Riddler, L. Swenson, I. Tao, S. Szeto, D. Chan, C. Kadie, N. Frahm, C. Brander, B. D. Walker, Z. L. Brumme, P. R. Harrigan, and S. Mallal, presented at the 15th International Workshop on HIV Dynamics and Evolution, Santa Fe, NM, 27 to 30 April 2008), HLA-associated polymorphisms are dichotomized based on the presence or absence of selection pressure: escape amino acids represent the residue most likely to emerge under HLA-restricted CTL selection pressure, while reversion amino acids represent the immunologically susceptible residues most likely to reemerge following transmission to an HLA-unmatched individual. In general, the reversion amino acids at any particular HLA-associated position tend to be consensus, while the escaped forms tend to be nonconsensus, but this is not always the case. The list contains over 1,000 unique HLA-associated polymorphisms occurring at
20% of codons in Gag/Pol and
50% of codons in Nef (Heckerman et al., 15th International Workshop on HIV Dynamics and Evolution). Of these,
30% map within optimally defined CTL epitopes (http://www.hiv.lanl.gov/content/immunology/tables/optimal_ctl_summary.html) (Fig. 1).
![]() View larger version (23K): [in a new window] |
FIG. 1. Locations of HLA-associated substitutions within optimally defined CTL epitopes. Published sequences of the 71 optimally defined CTL epitopes harboring at least one HLA-associated polymorphic site, based on analysis of over 1,200 persons with chronic infection, are listed. HLA-associated polymorphic residues are indicated in red. Overlapping and/or variant epitopes restricted by the same HLA allele were removed to avoid double counting of escape mutations. A complete list of all optimally defined epitopes is available at http://www.hiv.lanl.gov/content/immunology/tables/optimal_ctl_summary.html.
|
180 optimally defined epitopes in Gag, Pol, and Nef, specific sites and pathways of HLA-driven substitutions were defined for 74 epitopes. Of these 74, 3 were restricted by HLA alleles not observed in the seroconverter cohort and were excluded from analysis; the remaining 71 epitopes, along with their escape sites, are listed in Fig. 1. Data analysis and definitions. (i) Rates of escape and frequencies of epitope recognition in early infection. The time to escape for each CTL epitope was defined as the number of days that elapsed between the estimated infection date and the first detection of a full or partial amino acid change consistent (for HLA, codon, and the direction of amino acid selection) with the predefined list of HLA-associated polymorphisms (Fig. 1). Note that this definition is independent of whether an amino acid is consensus or nonconsensus: in the presence of the restricting HLA, any observed full or partial amino acid change away from the reversion form and/or toward the escaped form was considered escape. Note that some epitopes (e.g., B*57 TW10) have more than one escapable site; in these cases, we computed the time to the first escape event at any of these sites. In addition, the baseline sequence was analyzed for the presence of HLA-associated polymorphisms to identify cases where escape may have occurred very early (or, indistinguishably, rare cases where an escaped variant may have been acquired at transmission). When such a polymorphism was present at baseline in a person with the restricting allele, it was counted as an escape mutation.
Corresponding epitope recognition frequencies were derived from an independent, partially published, cross-sectional data set of 289 individuals enrolled through the AIEDRP network who were screened for HIV-1-specific CTL responses against optimally defined epitopes using a gamma interferon (IFN-
) enzyme-linked immunospot (ELISpot) assay at a single time point
8 weeks following the initial presentation with primary HIV infection, as described previously (7, 76); this effectively means that individuals were tested anywhere between
3 and 8 months after the estimated infection date. In addition, a subset of seroconverters for whom peripheral blood mononuclear cell samples were available during study follow-up (n = 8) were longitudinally investigated for HLA-restricted CTL responses to optimally described epitopes using IFN-
ELISpot assays.
(ii) Estimating rates and frequencies of reversion, and the total proportion of amino acid substitutions attributable to HLA-associated selection pressures. Rates and frequencies of reversion were conservatively estimated using methods similar to those used for escape, with one important difference. Since the vast majority of reversions are toward the consensus, it was not possible to infer reversion events occurring prior to the baseline. In this analysis, therefore, time zero was defined as the estimated infection date, and the baseline sequence was treated as the transmitted sequence. Reversion was defined as the presence of a specific known HLA-associated polymorphism at the baseline time point in an individual not bearing this HLA allele, followed by full or partial reversion toward the predefined immunologically susceptible (usually the consensus) amino acid during follow-up.
Similarly, the proportion of overall amino acid substitutions attributable to HLA-associated selection pressures (escape and/or reversion) on a gene-wide basis was conservatively calculated as the fraction of the total observed substitutions that achieved an exact match to an HLA-associated polymorphism in the predefined list.
|
|
|---|
180 published optimally defined epitopes in Gag, Pol, and Nef (Heckerman et al., 15th International Workshop on HIV Dynamics and Evolution). Of these, 71 were restricted by HLA alleles observed in the presently described seroconverter cohort (Fig. 1) and, thus, were evaluated for rates of escape in the first year of infection. Epitope escape rates were characterized using Kaplan-Meier methods, where the time to escape was defined as the time elapsed between the estimated infection date and the first observation of a full or partial substitution that was consistent with the predefined list of HLA-associated polymorphisms. Evidence for escape was observed in 57 of 71 (80%) of these epitopes at the 1-year time point (Fig. 2). Escaping epitopes were relatively equally distributed across Gag, Pol, and Nef. Among evolving epitopes, escape was generally rapid: 42 (74%) of evolving epitopes escaped within the first 6 months of infection (an estimate that included cases where the specific HLA-associated polymorphism was already present at baseline in the context of the restricting allele).
![]() View larger version (18K): [in a new window] |
FIG. 2. First-year rates of escape among HLA-restricted, optimally defined CTL epitopes. The rates of escape in optimally defined CTL epitopes in Gag (blue diamonds), Pol (red squares), and Nef (green triangles) in the first year of infection are shown. Vertical bars indicate 95% confidence intervals. Epitopes were restricted to those containing a predefined HLA-associated polymorphism (Fig. 1). Note that this analysis incorporates an estimate of cases where the epitope may have escaped prior to baseline sampling, performed by analyzing each individual's baseline HIV sequence for known HLA-restricted escape mutations. Note that in doing so, we cannot discriminate cases of the very early selection of escape mutations from rare cases where an escaped variant may have been acquired at transmission.
|
![]() View larger version (29K): [in a new window] |
FIG. 3. Rates of CTL escape generally reflect known immunodominance hierarchies of CTL epitope recognition in primary infection. The recognition frequencies of specific HLA-restricted epitopes (assessed by IFN- ELISpot assays measuring CTL responses to optimally described HIV peptides in a cohort of 289 individuals in primary HIV infection) (7, 76) are indicated in the left panels. Kaplan-Meier curves representing the corresponding rates of escape are indicated in the right panels in matching colors. Note that there are two cases where Kaplan-Meier curves for a pair of epitopes are superimposed so that only one curve is visible; these are B*08 EV9 and DL9 (which do not escape) and B*07 FL9 and TM9 (which overlap and share a single escaping site; Fig. 1).
|
In order to characterize the relationship between the frequency of CTL responses and the corresponding frequencies of escape, we analyzed a data set of 289 individuals in primary infection screened for CTL responses against optimally defined epitopes using an IFN-
ELISpot assay (7, 76). As expected, the rates of escape generally mirrored epitope-targeting frequencies (Fig. 3). For example, with the exception of KF11, the escape order for B*57-restricted epitopes generally reflected the established immunodominance hierarchies in primary infection (7), with TW10 representing the most frequently targeted (70%) and most rapidly escaping B*57-restricted epitope. HLA-A*11 and HLA-B*08 alleles also illustrate this relationship (Fig. 3). However, for some alleles, notably B*07, we observed a relatively poor relationship between the frequency of recognition and escape (Fig. 3).
Overall, we observed a robust positive correlation between frequencies of epitope recognition and escape in Gag (Spearman's r = 0.5; P = 0.01) and Pol (r = 0.8; P = 0.002) epitopes, highlighting the ability of HIV to rapidly adapt to CTL-induced selection pressures, even in more conserved regions of the viral proteome (Fig. 4a, b). Of interest, B*27-KK10, the most frequently (90%) targeted epitope in early infection, escaped at a slower rate relative to other frequently targeted epitopes, consistently with previous reports (35, 39, 49). B*27-KK10 ranked as the 20th most rapidly evolving epitope overall, with 5 of 11 B*27-positive individuals selecting the L268M mutation during follow-up and an additional two exhibiting both R264K and L268M at the earliest time point (indicating either the transmission of these mutations or very early escape within this epitope [10]).
![]() View larger version (15K): [in a new window] |
FIG. 4. Correlation between frequencies of recognition and escape by protein. Spearman's rank correlation was used to characterize the relationship between the frequencies of recognition of CTL epitopes in Gag, Pol, and Nef with their corresponding frequencies of escape in the first year of infection. A regression line was drawn to highlight trends. All tested CTL epitopes are represented; those recognized and/or evolving at frequencies of > 40% are labeled with the epitope name and HLA restriction.
|
We observed a number of cases where epitope escape frequencies exceeded recognition frequencies in Gag and Pol (for example, A*24-KW9 [Gag] and A*03-QR9 [RT], the second and sixth most rapidly evolving epitopes overall, as well as a number of B*07-restricted epitopes; Fig. 4). This type of discordance, however, was observed most frequently in Nef and contributed to the poor correlation between recognition and escape frequencies in this protein (r = 0.09; P = 0.7) (Fig. 4c).
Rates of reversion in CTL epitopes. Reversion was conservatively defined as the presence of a known HLA-associated polymorphism at the baseline time point in an individual not bearing this HLA allele, followed by full or partial reversion toward the immunologically susceptible (usually the consensus) amino acid during follow-up. Within the first year of infection, reversions were observed at 29 (6%), 14 (3%), and 26 (13%) Gag, Pol and Nef codons, respectively, while an additional 9, 7, and 10 sites reverted after the first year. The majority of reverting codons (65%) were within published HLA-restricted CTL epitopes. On average, reversions within Gag epitopes restricted by protective HLA alleles (B*13, B*51, B*57, and B*58) occurred more rapidly than reversions outside these regions (P = 0.02); however, no such differences were observed in Pol or Nef. The rates of reversion at known escape sites within optimally described epitopes in the first year of infection are summarized in Fig. 5. The most rapidly reverting mutations within Gag CTL epitopes were B*57-associated escape mutations at codons 147 and 242 in the IW9 and TW10 epitopes, respectively, while the most rapidly reverting epitope-associated mutation within Nef was at codon 135 (residue 2 of A*24-RW8). Although it is important to emphasize that the rates of escape and reversion are not directly comparable due to the differences in calculations (see Materials and Methods and Discussion), it is nevertheless interesting that these three mutations were relatively rapidly escaping as well, implying a substantial fitness cost. However, the factors that shape HIV evolution are complex and, thus, one would not expect rates of escape to correlate with rates of reversion in all cases. RT codon 135 (residue 8 of B*51-TI8), for example, was not observed to revert despite the relatively frequent transmission of the escaped variant (n = 21), implying a negligible fitness cost to this mutation.
![]() View larger version (20K): [in a new window] |
FIG. 5. Locations and first-year rates of commonly observed reversions within published HLA-restricted epitopes in Gag, Pol, and Nef. Conservative estimates of the rates of reversion of known HLA-associated escape mutations in the absence of the restricting HLA (expressed as percent reversions/person-month) at Gag, Pol, and Nef codons within published HLA-restricted CTL epitopes are shown, along with the total number of observations (transmitted mutations) at each codon. A minimum of five observed cases of the transmission of the escaped variant was required for display. Note that in contrast to the escape analysis, this analysis does not take into consideration reversions that may have occurred prior to the baseline sampling. Thus, estimated reversion rates represent considerable underestimates of the true rate and are not directly comparable to epitope escape rates. PR, protease.
|
70% of these observations. |
|
|---|
The current study overcomes many of the limitations of previous investigations in this area, as it represents the first large-scale study of HLA-driven HIV evolution in acute clade B infection that employs a rational list of HLA-associated polymorphisms (defined through the analysis of a large [n > 1,200] independent data set) to systematically characterize the rates of immune-driven evolution on a protein-wide basis in the context of all HLA restrictions. The study focused on Gag, Pol, and Nef, because these proteins have been featured in candidate CTL vaccine design strategies, including the recent failed STEP vaccine trial (74); thus, developing a more in-depth understanding of CTL-driven evolution in these proteins is of key importance. We estimate that a minimum of 30 to 35% of nonsynonymous substitutions in Gag/Pol and a minimum of 60% of substitutions in Nef are attributable to HLA-restricted immune pressures in the first year of infection (of which the majority represent reversions), confirming dramatically different levels of HLA-mediated adaptation in HIV proteins (17).
Using Kaplan-Meier analysis, we systematically computed the rates of escape in published CTL epitopes, revealing markedly different kinetics of escape among them. Although we were not able to evaluate CTL responses in this cohort due to a lack of cryopreserved cells, we were able to integrate the sequence data derived here with epitope recognition frequencies characterized in a cohort of 289 persons with primary HIV infection (7, 76). A robust correlation between the epitope response frequencies and the rates of escape in Gag and Pol was observed, illustrating that in general, if the immune response mounts pressure against an epitope, in most cases HLA-driven mutations will appear shortly thereafter. There were relatively few examples of frequently targeted epitopes that maintained their sequences for extended periods of time without evidence of escape.
Of interest, protective HLA alleles (21, 41, 47) were overrepresented among rapidly escaping epitopes, consistent with the observation that these alleles impose stronger in vivo immune selection pressures than others (33, 72). In particular, B*57, the strongest host genetic factor associated with protection against HIV progression (21), restricted the most rapidly escaping epitope overall (TW10), and all but one of the remaining B*57-restricted epitopes displayed rapid escape as well. The somewhat paradoxical observation that certain alleles remain protective despite rapid and frequent escape may be explained by the fact that a key correlate of protection is the ability to mount a broad and robust CTL response very early in infection (4, 64), which indirectly manifests itself as CTL-driven sequence changes at these key sites. On average, B*57-expressing individuals respond to a mean of 4.7 B*57-restricted epitopes (genome-wide) in early infection, the highest number for all HLA restrictions, compared to a mean of two or three epitopes for the vast majority of other alleles (7, 76). From the perspective of HLA imprinting on HIV sequences at the population level, HLA-B*57 represents one of the alleles with the largest total number of identified HLA-associated polymorphic sites in Gag/Pol/Nef, with associations identified at 29 unique residues in these proteins (compared to 17 and 18 for B*07 and B*08, respectively) (Heckerman et al., 15th International Workshop on HIV Dynamics and Evolution). Furthermore, the earliest CTL responses likely are the most potent and possess the greatest antiviral potential; in no other stage of disease does one observe such a dramatic decline in viremia from the acute peak to set-point levels (19).
A number of additional reasons may explain why certain alleles remain protective despite the rapid selection of HLA-driven substitutions. First, not all of the identified HLA-associated polymorphisms confer an absolute escape phenotype: often, the variant will maintain at least some capacity to bind HLA and/or be cross-recognized by CTL (77), and there is evidence that B*57-restricted CTL have a superior ability to cross-recognize peptide variants compared to that of less protective alleles (77). Furthermore, the ability to mount a de novo CTL response against a selected variant (3, 30, 78) may be greater in earlier disease while immune function is relatively intact. Finally, recent data suggest that the breadth of CTL responses to Gag contributes strongly to HIV control (1, 28, 50, 79), and that the fitness costs of escape mutations in Gag are substantial (16, 23, 62, 73), suggesting that long-term protective effects are due to strong immune selective pressures driving viral evolution toward less-fit forms (6). Indeed, the observed rapid and frequent reversion of escape mutations within B*57-TW10 and IW9 in Gag supports this hypothesis (Fig. 4). Taken together, these results support the importance of a strong early CTL response in the control of HIV viremia, even if a consequence of this is the rapid selection of mutations.
The observation of discordant cases where in vivo escape frequencies exceeded in vitro recognition frequencies merits discussion. The fact that epitope recognition (7, 76) and escape were evaluated in independent cohorts accounts for a portion of these discordances (escape frequencies in the seroconverter cohort represent cumulative longitudinal frequencies calculated at 1 year postinfection, whereas recognition frequencies represent measurements taken at a single time point between 3 and 8 months following the estimated date of infection). The remainder may be attributed to an underestimation of CTL responses and/or an overestimation of escape frequencies, as follows. First, as CTL responses to wild-type peptides decline following escape in the autologous viral sequence (2), the recognition frequencies of rapidly escaping epitopes could have been underestimated if subjects were screened after escape had occurred. Indeed, for the eight seroconverters for whom longitudinal peripheral blood mononuclear cells were available, we observed a decline in responses to B*51 TI8-RT following escape in both B*51-expressing individuals in this subgroup (not shown). Similarly, CTL responses may be underestimated in cases where the autologous viral sequence differs from the tested epitope sequence (5, 62), a fact that would more substantially affect variable proteins such as Nef.
Another contributing factor may be the overestimation of escape. This is particularly an issue in cases where the identified escape variant for a given epitope represents a commonly occurring (or in some cases, the consensus) residue. Both A*24 KW9-Gag and A*03 QR9-RT are examples of such cases: for both epitopes, the immunologically susceptible form (position 1 of KW9 and position 9 of QR9, corresponding to Gag and RT codons 28 and 277, respectively) represents a nonconsensus amino acid, meaning that the subtype B consensus residue is considered the escaped form. The fact that our analysis (Fig. 2 to 4) included subjects exhibiting an escaped residue at baseline could have resulted in an overestimation of escape for these epitopes if these mutations were present at transmission. Finally, it is worth noting that in both cases (as well as for all B*07-restricted epitopes examined here), the published epitope sequence features at least one residue in its HLA-associated escaped form. Work is ongoing to assess whether these cases represent examples where the population HIV consensus has adapted to frequently observed HLA alleles (57); in any event, this observation merits consideration, as it could result in the systematic underestimation of CTL response frequencies to such epitopes.
There are a number of additional limitations that are important to discuss. Most importantly, in the absence of knowing the transmitted sequence, we cannot definitively identify sequence changes occurring prior to baseline sampling. Although the rates of escape analysis take into consideration likely escape events occurring prior to baseline sampling, it is not possible to infer reversion events (or to calculate the total amount of evolution) occurring prior to baseline sampling. Thus, not only does this limitation render the reversion rates to be substantial underestimates of the true rates, it also renders them not directly comparable to the escape rates. Nevertheless, the analysis still yields informative data regarding the locations, frequencies, and estimated rates of reversion, which illuminate potential positions at which mutations likely occur at substantial costs to fitness. For example, it was interesting that two of the most rapidly reverting codons were Gag 147 and 242 within B*57-TW10 and IW9 epitopes, respectively (39, 58), consistent with measurable costs to viral fitness associated with escape in B*57-restricted Gag epitopes (16, 62). Conversely, the rapidly selected B*51-associated escape mutation at RT codon 135 (residue 8 of B*51-TI8) is not observed to revert despite relatively frequent transmission, suggesting a relatively minor fitness cost.
If the rates of reversion and/or escape differ substantially before and after baseline sampling, this difference could contribute to errors in our estimates. Note, however, that no significant difference in the proportion of HLA-attributable evolution was observed among individuals captured within <3 (n = 61) or >3 to 6 months (n = 37) after infection. The inconsistent detection of minority variants below a threshold of
10 to 20% of the circulating species is a known limitation of the bulk PCR and sequencing techniques employed (55, 56); however, comprehensive clonal sequencing was not feasible in a cohort of this size and length of follow-up.
In addition, there are some limitations associated with defining HLA-associated substitutions through the analysis of an independent large clade B data set. As with any method, there will be both false-positive as well as false-negative results within this list (20) (J. Carlson, submitted for publication); thus, this list will not necessarily contain all escape mutations previously reported in the literature. Potential differences in the ethnic composition in the two cohorts also should be acknowledged. Information on ethnicity was unavailable; however, the HLA composition of both cohorts reflected expected allele frequencies among North American populations (not shown). Finally, although many mutational pathways are broadly predictable at the population level (17, 18, 65, 70), unique HLA-driven escape, reversion, or secondary/compensatory changes also will occur in individual patients; thus, employing a population-level definition of HLA-associated substitution at the individual level may underestimate the amount of evolution attributable to HLA. Nevertheless, the use of a predefined list of HLA-associated substitutions allows the comprehensive, systematic classification of amino acid substitutions both within and outside CTL epitopes, as well as the ability to investigate escape and reversion across entire proteins and across all HLA restrictions. Indeed, a similar classification technique was employed in a recent study of transmitted CTL escape mutations in HIV clade C (36). The current study investigated HLA-associated polymorphisms in Gag, Pol, and Nef only; however, as comprehensive lists of HLA-associated polymorphisms become available for additional genes, the analysis of escape and reversion rates in other HIV proteins will become possible. Finally, it is important to acknowledge the incompleteness as well as potential bias toward common and/or protective HLA alleles in the current published CTL epitope lists. As new epitopes are discovered, however, the rates of escape and reversion can easily be calculated using the current data set.
The relatively short follow-up period limits our ability to assess the long-term impact of early escape on HIV disease outcomes. In addition, we were unable to confirm the results of recent studies that demonstrate reduced viral loads in individuals transmitting HLA-associated escape mutations (22, 36); however, the lack of statistical power to detect these effects in the current study must be noted.
Given these limitations, our estimates that
30 to 35% and
60% of overall substitutions in Gag/Pol and Nef, respectively, are attributable to CTL escape and reversion confirm a substantial role of HLA-associated immune pressures in driving early within-host HIV evolution (2, 9, 39, 59). Selective forces responsible for the non-HLA-attributable fraction may include CD4+ T-cell responses (48, 71), adaptations to other host factors (12), reversions of such mutations selected in previous hosts, or random drift. A small number of substitutions in protease/RT (and possibly Gag (26) could be due to the reversion of transmitted resistance mutations. Indeed, a small number of resistance-associated polymorphisms and surveillance mutations (54) were observed in this cohort; note, however, that there was no overlap between major resistance-associated and HLA-associated polymorphic sites in protease/RT.
Substantial variability in the degree to which different viral proteins adapt to HLA-associated immune pressures during early HIV infection underscores the importance of selecting appropriate immunogens and determining how to best incorporate sequence diversity in vaccine design (5, 15, 32, 66, 69). In particular, strongly targeted yet slowly escaping epitopes (such as B*57-KF11 and A*25-EW10) may be of particular interest, as may be strongly targeted epitopes that can escape only at substantial fitness costs to the virus (such as B*57-TW10) (16). Specifically, immunogenic viral regions exhibiting high mutational barriers to escape may be good vaccine candidates, as they may mediate the effective long-term control of HIV CTL (39). On the other hand, epitopes that escape at considerable costs to viral replicative capacity also may be relevant to vaccine design: a vaccine capable of inducing immune responses that drive HIV toward crippled forms also may be effective in reducing viral replication to levels that slow disease progression (6). Indeed, the observation that HLA-B*57-restricted TW10 and IW9 are among the fastest escaping and the fastest reverting epitopes in Gag strongly supports the fitness costs of CTL escape as a key mediator of long-term viremia control in B*57-expressing individuals.
Taken together, our results confirm a substantial role of HLA-associated selection pressures on early within-host HIV evolution and support further research into CTL-based vaccine strategies incorporating information on common escape pathways, despite recent setbacks in the HIV vaccine field (8, 74).
We thank Toshiyuki Miura, Mark Brockman, Nicole Frahm, and Christian Brander for helpful discussions and Philip Goulder for the critical reading of the manuscript.
Z.L.B. is supported by a postdoctoral fellowship from the Canadian Institutes for Health Research (CIHR). H.S. is supported by a fellowship from the Deutscher Akademischer Austauschdienst (DAAD). D.H., J.C., and C.K. are funded by Microsoft Research. This study was supported in part by the Howard Hughes Medical Institute, the Mark and Lisa Schwartz Foundation, the Harvard University Center for AIDS Research (CFAR), National Institutes of Health grant R01AI50429, and Acute Infection Early Disease Research Network (AIEDRP) grant U01AI052403. This project has been funded in part with federal funds from the National Cancer Institute, National Institutes of Health, under contract N01-CO-12400. This research was supported in part by the Intramural Research Program of the NIH, National Cancer Institute, Center for Cancer Research.
Sources of support were not involved in the design and conduct of the study, nor were they involved in the collection, analysis, and interpretation of the data or in the preparation, review, or approval of the manuscript.
The content of this publication does not necessarily reflect the views or policies of the Department of Health and Human Services, nor does the mention of trade names, commercial products, or organizations imply endorsement by the U.S. government.
Published ahead of print on 9 July 2008. ![]()
These authors contributed equally to this work. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»