Previous Article
Journal of Virology, July 2005, p. 9363-9366, Vol. 79, No. 14
0022-538X/05/$08.00+0 doi:10.1128/JVI.79.14.9363-9366.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Intrapatient Escape in the A*0201-Restricted Epitope SLYNTVATL Drives Evolution of Human Immunodeficiency Virus Type 1 at the Population Level
Charles T. T. Edwards,1
Katja J. Pfafferott,1
Philip J. R. Goulder,1
Rodney E. Phillips,1 and
Edward C. Holmes2*
Nuffield Department of Clinical Medicine, University of Oxford, The Peter Medawar Building for Pathogen Research, South Parks Road, Oxford OX1 3SY,1
Department of Zoology, University of Oxford, South Parks Road, Oxford OX1 3PS, United Kingdom2
Received 29 November 2004/
Accepted 26 March 2005

ABSTRACT
The hypothesis that the intrapatient emergence of cytotoxic
T-lymphocyte escape variants contributes to the evolution of
human immunodeficiency virus type 1 at the population (interpatient)
level was tested using the HLA-A*0201-restricted
gag p17 epitope
SLYNTVATL. Using a simple experimental design, we investigated
the evolutionary processes operating within this epitope among
patients while compensating for the confounding influence of
intrapatient natural selection. Using this approach, we revealed
a pattern of A*0201-driven escape within patients, followed
by the sustained transmission of these escape variants among
patients irrespective of their HLA type.

TEXT
During infection with human immunodeficiency virus type 1 (HIV-1),
the cytotoxic T-lymphocyte (CTL) response plays a central role
in controlling viral replication (
3,
16,
17,
27) and has been
shown to drive the emergence of escape variants capable of evading
CTL recognition (
4,
13,
22). If these escape variants are continually
transmitted, it is possible that the CTL response will also
influence HIV-1 genetic variation at the population (interpatient)
level. The accumulation of escape variants across the infected
population is supported by the observation that rare HLA supertypes
are associated with improved control of viral replication and
slower disease progression (
26). However, although the horizontal
transmission of escape variants has been inferred from the detection
of CTL escape mutants early in the acute phase of infection
(
1,
7,
11,
18), whether this process has a major effect on HIV-1
evolution at the population level is unclear. If escape variants
carry a fitness cost in patients unable to mount a response
against the epitope in question (i.e., who lack the appropriate
HLA type), then the persistent transmission of these variants
is unlikely. Indeed, the reversion of CTL escape mutations in
patients with an incompatible genetic background has been demonstrated
(
18), and the need for compensatory mutations suggests that
some CTL escape mutations incur a cost (
15). In the analogous
situation of drug-resistant variants, compensatory mutations
may enhance the fitness of resistant strains (
9) but fail to
prevent reversion to susceptible genotypes upon transmission
to drug-naive recipients (
8).
How intrapatient CTL escape contributes to evolution at the population level is crucial in understanding long-term patterns of change in HIV. To address this question, we examined the HLA-A*0201-restricted gag p17 epitope SLYNTVATL (SL9) in subtype B patients. CTL escape has been intensively studied in this epitope (5, 10, 12, 14), and the A*0201 HLA type is by far the most common within Caucasians (44% phenotype frequency) (19), making it a strong candidate to influence evolution at the population level.
By combining multiple patients within a single analysis, it is possible to analyze evolutionary dynamics at the population level. However, it is first necessary to exclude the confounding effect of selection operating within each of the individual patients. This was achieved by only analyzing those patients unable to mount a CTL response against SL9, which enabled us to determine whether A*0201+-associated escape mutants were transmitted through the population because of their enhanced fitness or emerged within each patient individually.
In our initial analysis, we confirmed that intrapatient positive selection within SL9 is limited to A*0201+ patients. Sequences were obtained from two previously published studies (5, 6) (accession numbers AF017813 to AF017980, AF028563 to AF028587, AF060031 to AF060073, and AF073382 to AF073441), totaling 18 HLA-typed patients, with the number of clones from each ranging from 7 to 43 (average, 15.6) (Table 1). First, sequences were aligned manually using the Se-Al program (23). Maximum-likelihood phylogenetic trees were then estimated using the PAUP* package (25) assuming the general time-reversible (GTR + I +
) model of nucleotide substitution, with all parameters estimated from the data (available on request). For each patient phylogeny, an ancestral reconstruction using maximum likelihood was performed with the DATAMONKEY program (21). The number of synonymous (dS) and nonsynonymous (dN) substitutions per site for each codon was then determined. Under positive selection, we expect dN > dS. A per-site comparison of the difference between dN and dS in the A*0201+ and A*0201 patient groups allowed us to test whether the strength of positive selection was the same in each group. For amino acids 71 to 91 in gag p17, we found a significant difference only at T84 (A*0201 mean dN dS = 0.22; A*0201+ mean dN dS = 0.60; P = 0.010 using a one-tailed t test) (Table 1). Although not an explicit test of intrapatient positive selection, this demonstrates that nonsynonymous changes at T84 are limited largely to A*0201+ patients. Changes at this site have been shown to facilitate escape from A*0201 (10).
We next examined whether transmission of T84 variants influences
evolution at the population level. We combined sequences from
multiple A*0201
patients within a single analysis (accession
numbers
AY656030 to
AY656080) (
18). By excluding A*0201
+ patients,
this ensures that any selection observed is unlikely to be the
consequence of multiple independent selection events operating
within the patient group. Full-length
gag sequences were collected
from a total of 36 subtype B-infected patients visiting clinics
in and around the London area (
18). Of these, 17 were A*0201
and should therefore be incapable of mounting an SL9-restricted
response. We tested for positive selection among the A*0201
patients by estimating d
N/d
S for each codon using the CODEML
program (
28). A comparison of two models of codon evolution
allows investigation of selection pressures based on their fit
to the data. In model 8 (M8), codons are distributed into 11
categories of d
N/d
S, 10 of which have a value of

1, with an
extra category for which d
N/d
S is estimated freely and can therefore
be >1 (implying positive selection). This was compared to
the neutral M8N model in which the 11th d
N/d
S category is fixed
at 1 (implying neutrality), using a likelihood ratio test with
one degree of freedom. Evidence for positive selection is provided
if M8 is a significantly better description of the data than
M8N and d
N/d
S for the 11th category is >1. The posterior
probability that the assignment of sites to this category is
correct is estimated using Bayesian inference. To ensure that
our analysis is not compromised by frequent recombination (
2),
we tested for genetic exchange using the GENECONV program prior
to selection analysis (
24). No evidence for recombination was
found in the A*0201
data set (results available on request).
We found evidence for positive selection at site T84 (d
N/d
S = 1.8; posterior probabilities, >0.95; Table
2; the result
from the analysis of all 36 patients is also shown). Selection
was also observed at Y79, a site which is similarly associated
with escape from A*0201 (
10). This supports the notion that
mutations at these sites are advantageous in A*0201
+ patients
and sustained in the wider population through interpatient transmission.
SL9 is also overlapped by epitopes restricted by A*01, A*1101,
A*3002, and B*0801 (Fig.
1). To exclude their possible influence,
we repeated the analyses using patients negative for all these
alleles, as well as A*0201. To obtain the necessary numbers,
we combined sequences from multiple cohorts in the United Kingdom,
Boston, and Barbados, using 21 patients out of a total of 92
(
18). This analysis again revealed evidence for selection at
sites Y79 and T84 (d
N/d
S = 2.4; posterior probabilities, >0.95;
Table
2), confirming that the selection observed is operating
at a level above the individuals from whom the sequences were
derived. To further validate our hypothesis of population level
evolution, we tested the assumption of sustained transmission
by mapping the selected sites onto a phylogenetic tree. Using
the principle of parsimony, we were able to show that changes
at both the Y79 and T84 sites were located on branches that
precede the divergence between patients, consistent with our
assumption that these variants are transmitted (Fig.
2).
Although we were able to detect positive selection operating
within SL9 at the population level, that A*0201 is the allele
responsible for driving this diversification is only clear for
T84, since nonsynonymous changes at this site are limited to
A*0201
+ patients. Nevertheless, the high frequency of A*0201
compared to other alleles that restrict overlapping epitopes
(Fig.
1) (
19), and the well-documented role of A*0201 in selecting
for intrapatient escape through variation in Y79 and T84 (
10,
14,
20), is supportive of our hypothesis of positive selection
at both sites.
In summary, our analyses demonstrate that escape from the A*0201-restricted CTL response is also having a major effect on evolution within the SL9 epitope at the wider population level. Our study therefore illustrates the potential impact of intrapatient immune selection on the long-term evolution of HIV-1.

ACKNOWLEDGMENTS
This work was supported by the Wellcome Trust.

FOOTNOTES
* Corresponding author. Present address: Department of Biology, The Pennsylvania State University, Mueller Laboratory, University Park, PA 16802. Phone: (814) 863-4689. Fax: (814) 865-9131. E-mail:
ech15{at}psu.edu.


REFERENCES
1 - Allen, T. M., M. Altfeld, X. G. Yu, K. M. O'Sullivan, M. Lichterfeld, S. Le Gall, M. John, B. R. Mothe, P. K. Lee, E. T. Kalife, D. E. Cohen, K. A. Freedberg, D. A. Strick, M. N. Johnston, A. Sette, E. S. Rosenberg, S. A. Mallal, P. J. Goulder, C. Brander, and B. D. Walker. 2004. Selection, transmission, and reversion of an antigen-processing cytotoxic T-lymphocyte escape mutation in human immunodeficiency virus type 1 infection. J. Virol. 78:7069-7078.[Abstract/Free Full Text]
2 - Anisimova, M., R. Nielsen, and Z. Yang. 2003. Effect of recombination on the accuracy of the likelihood method for detecting positive selection at amino acid sites. Genetics 164:1229-1236.[Abstract/Free Full Text]
3 - Borrow, P., H. Lewicki, B. H. Hahn, G. M. Shaw, and M. B. Oldstone. 1994. Virus-specific CD8+ cytotoxic T-lymphocyte activity associated with control of viremia in primary human immunodeficiency virus type 1 infection. J. Virol. 68:6103-6110.[Abstract/Free Full Text]
4 - Borrow, P., H. Lewicki, X. Wei, M. S. Horwitz, N. Peffer, H. Meyers, J. A. Nelson, J. E. Gairin, B. H. Hahn, M. B. Oldstone, and G. M. Shaw. 1997. Antiviral pressure exerted by HIV-1-specific cytotoxic T lymphocytes (CTLs) during primary infection demonstrated by rapid selection of CTL escape virus. Nat. Med. 3:205-211.[CrossRef][Medline]
5 - Brander, C., K. E. Hartman, A. K. Trocha, N. G. Jones, R. P. Johnson, B. Korber, P. Wentworth, S. P. Buchbinder, S. Wolinsky, B. D. Walker, and S. A. Kalams. 1998. Lack of strong immune selection pressure by the immunodominant, HLA-A*0201-restricted cytotoxic T lymphocyte response in chronic human immunodeficiency virus-1 infection. J. Clin. Investig. 101:2559-2566.[Medline]
6 - Brander, C., O. O. Yang, N. G. Jones, Y. Lee, P. Goulder, R. P. Johnson, A. Trocha, D. Colbert, C. Hay, S. Buchbinder, C. C. Bergmann, H. J. Zweerink, S. Wolinsky, W. A. Blattner, S. A. Kalams, and B. D. Walker. 1999. Efficient processing of the immunodominant, HLA-A*0201-restricted human immunodeficiency virus type 1 cytotoxic T-lymphocyte epitope despite multiple variations in the epitope flanking sequences. J. Virol. 73:10191-10198.[Abstract/Free Full Text]
7 - Furutsuki, T., N. Hosoya, A. Kawana-Tachikawa, M. Tomizawa, T. Odawara, M. Goto, Y. Kitamura, T. Nakamura, A. D. Kelleher, D. A. Cooper, and A. Iwamoto. 2004. Frequent transmission of cytotoxic-T-lymphocyte escape mutants of human immunodeficiency virus type 1 in the highly HLA-A24-positive Japanese population. J. Virol. 78:8437-8445.[Abstract/Free Full Text]
8 - Gandhi, R. T., A. Wurcel, E. S. Rosenberg, M. N. Johnston, N. Hellmann, M. Bates, M. S. Hirsch, and B. D. Walker. 2003. Progressive reversion of human immunodeficiency virus type 1 resistance mutations in vivo after transmission of a multiply drug-resistant virus. Clin. Infect. Dis. 37:1693-1698.[CrossRef][Medline]
9 - Garcia-Lerma, J. G., S. Nidtha, K. Blumoff, H. Weinstock, and W. Heneine. 2001. Increased ability for selection of zidovudine resistance in a distinct class of wild-type HIV-1 from drug-naive persons. Proc. Natl. Acad. Sci. USA 98:13907-13912.[Abstract/Free Full Text]
10 - Goulder, P. J., M. A. Altfeld, E. S. Rosenberg, T. Nguyen, Y. Tang, R. L. Eldridge, M. M. Addo, S. He, J. S. Mukherjee, M. N. Phillips, M. Bunce, S. A. Kalams, R. P. Sekaly, B. D. Walker, and C. Brander. 2001. Substantial differences in specificity of HIV-specific cytotoxic T cells in acute and chronic HIV infection. J. Exp. Med. 193:181-194.[Abstract/Free Full Text]
11 - Goulder, P. J., C. Brander, Y. Tang, C. Tremblay, R. A. Colbert, M. M. Addo, E. S. Rosenberg, T. Nguyen, R. Allen, A. Trocha, M. Altfeld, S. He, M. Bunce, R. Funkhouser, S. I. Pelton, S. K. Burchett, K. McIntosh, B. T. Korber, and B. D. Walker. 2001. Evolution and transmission of stable CTL escape mutations in HIV infection. Nature 412:334-338.[CrossRef][Medline]
12 - Goulder, P. J., C. Pasquier, E. C. Holmes, B. Liang, Y. Tang, J. Izopet, K. Saune, E. S. Rosenberg, S. K. Burchett, K. McIntosh, M. Barnardo, M. Bunce, B. D. Walker, C. Brander, and R. E. Phillips. 2001. Mother-to-child transmission of HIV infection and CTL escape through HLA-A2-SLYNTVATL epitope sequence variation. Immunol. Lett. 79:109-116.[CrossRef][Medline]
13 - Goulder, P. J., R. E. Phillips, R. A. Colbert, S. McAdam, G. Ogg, M. A. Nowak, P. Giangrande, G. Luzzi, B. Morgan, A. Edwards, A. J. McMichael, and S. Rowland-Jones. 1997. Late escape from an immunodominant cytotoxic T-lymphocyte response associated with progression to AIDS. Nat. Med. 3:212-217.[CrossRef][Medline]
14 - Goulder, P. J., A. K. Sewell, D. G. Lalloo, D. A. Price, J. A. Whelan, J. Evans, G. P. Taylor, G. Luzzi, P. Giangrande, R. E. Phillips, and A. J. McMichael. 1997. Patterns of immunodominance in HIV-1-specific cytotoxic T lymphocyte responses in two human histocompatibility leukocyte antigens (HLA)-identical siblings with HLA-A*0201 are influenced by epitope mutation. J. Exp. Med. 185:1423-1433.[Abstract/Free Full Text]
15 - Kelleher, A. D., C. Long, E. C. Holmes, R. L. Allen, J. Wilson, C. Conlon, C. Workman, S. Shaunak, K. Olson, P. Goulder, C. Brander, G. Ogg, J. S. Sullivan, W. Dyer, I. Jones, A. J. McMichael, S. Rowland-Jones, and R. E. Phillips. 2001. Clustered mutations in HIV-1 gag are consistently required for escape from HLA-B27-restricted cytotoxic T lymphocyte responses. J. Exp. Med. 193:375-386.[Abstract/Free Full Text]
16 - Klein, M. R., C. A. van Baalen, A. M. Holwerda, S. R. Kerkhof Garde, R. J. Bende, I. P. Keet, J. K. Eeftinck-Schattenkerk, A. D. Osterhaus, H. Schuitemaker, and F. Miedema. 1995. Kinetics of Gag-specific cytotoxic T lymphocyte responses during the clinical course of HIV-1 infection: a longitudinal analysis of rapid progressors and long-term asymptomatics. J. Exp. Med. 181:1365-1372.[Abstract/Free Full Text]
17 - Koup, R. A., J. T. Safrit, Y. Cao, C. A. Andrews, G. McLeod, W. Borkowsky, C. Farthing, and D. D. Ho. 1994. Temporal association of cellular immune responses with the initial control of viremia in primary human immunodeficiency virus type 1 syndrome. J. Virol. 68:4650-4655.[Abstract/Free Full Text]
18 - Leslie, A. J., K. J. Pfafferott, P. Chetty, R. Draenert, M. M. Addo, M. Feeney, Y. Tang, E. C. Holmes, T. Allen, J. G. Prado, M. Altfeld, C. Brander, C. Dixon, D. Ramduth, P. Jeena, S. A. Thomas, A. St. John, T. A. Roach, B. Kupfer, G. Luzzi, A. Edwards, G. Taylor, H. Lyall, G. Tudor-Williams, V. Novelli, J. Martinez-Picado, P. Kiepiela, B. D. Walker, and P. J. Goulder. 2004. HIV evolution: CTL escape mutation and reversion after transmission. Nat. Med. 10:282-289.[CrossRef][Medline]
19 - Marsh, S. G. E., P. Parnham, and L. D. Barber. 2000. The HLA facts book. Academic Press, Inc., New York, N.Y.
20 - Oxenius, A., D. A. Price, A. Trkola, C. Edwards, E. Gostick, H. T. Zhang, P. J. Easterbrook, T. Tun, A. Johnson, A. Waters, E. C. Holmes, and R. E. Phillips. 2004. Loss of viral control in early HIV-1 infection is temporally associated with sequential escape from CD8+ T cell responses and decrease in HIV-1-specific CD4+ 61 and CD8+ T cell frequencies. J. Infect. Dis. 190:713-721.[CrossRef][Medline]
21 - Pond, S. L., S. D. Frost, and S. V. Muse. 2005. HyPhy: hypothesis testing using phylogenies. Bioinformatics 21:676-679.[Abstract/Free Full Text]
22 - Price, D. A., P. J. Goulder, P. Klenerman, A. K. Sewell, P. J. Easterbrook, M. Troop, C. R. Bangham, and R. E. Phillips. 1997. Positive selection of HIV-1 cytotoxic T lymphocyte escape variants during primary infection. Proc. Natl. Acad. Sci. USA 94:1890-1895.[Abstract/Free Full Text]
23 - Rambaut, A. 2002. Se-Al. Sequence alignment editor v2.0a11. [Online.] http://evolve.zoo.ox.ac.uk/.
24 - Sawyer, S. A. 1999. GENECONV: a computer package for the statistical detection of gene conversion. [Online.] http://www.math.wustl.edu/
sawyer.
25 - Swofford, R. 2002. PAUP*. Phylogenetic analysis using parsimony (* and other methods), version 4 ed. Sinauer Associates, Sunderland, Mass.
26 - Trachtenberg, E., B. Korber, C. Sollars, T. B. Kepler, P. T. Hraber, E. Hayes, R. Funkhouser, M. Fugate, J. Theiler, Y. S. Hsu, K. Kunstman, S. Wu, J. Phair, H. Erlich, and S. Wolinsky. 2003. Advantage of rare HLA supertype in HIV disease progression. Nat. Med. 9:928-935.[CrossRef][Medline]
27 - Walker, B. D., S. Chakrabarti, B. Moss, T. J. Paradis, T. Flynn, A. G. Durno, R. S. Blumberg, J. C. Kaplan, M. S. Hirsch, and R. T. Schooley. 1987. HIV-specific cytotoxic T lymphocytes in seropositive individuals. Nature 328:345-348.[CrossRef][Medline]
28 - Yang, Z. 1997. PAML: a program package for phylogenetic analysis by maximum likelihood. Comput. Appl. Biol. Sci. 13:555-556.[Free Full Text]
Journal of Virology, July 2005, p. 9363-9366, Vol. 79, No. 14
0022-538X/05/$08.00+0 doi:10.1128/JVI.79.14.9363-9366.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
This article has been cited by other articles:
-
Brumme, Z. L., Brumme, C. J., Carlson, J., Streeck, H., John, M., Eichbaum, Q., Block, B. L., Baker, B., Kadie, C., Markowitz, M., Jessen, H., Kelleher, A. D., Rosenberg, E., Kaldor, J., Yuki, Y., Carrington, M., Allen, T. M., Mallal, S., Altfeld, M., Heckerman, D., Walker, B. D.
(2008). Marked Epitope- and Allele-Specific Differences in Rates of Mutation in Human Immunodeficiency Type 1 (HIV-1) Gag, Pol, and Nef Cytotoxic T-Lymphocyte Epitopes in Acute/Early HIV-1 Infection. J. Virol.
82: 9216-9227
[Abstract]
[Full Text]
-
Kan-Mitchell, J., Bajcz, M., Schaubert, K. L., Price, D. A., Brenchley, J. M., Asher, T. E., Douek, D. C., Ng, H. L., Yang, O. O., Rinaldo, C. R. Jr., Benito, J. M., Bisikirska, B., Hegde, R., Marincola, F. M., Boggiano, C., Wilson, D., Abrams, J., Blondelle, S. E., Wilson, D. B.
(2006). Degeneracy and Repertoire of the Human HIV-1 Gag p1777-85 CTL Response.. J. Immunol.
176: 6690-6701
[Abstract]
[Full Text]