Previous Article | Next Article ![]()
Journal of Virology, July 2006, p. 7226-7234, Vol. 80, No. 14
0022-538X/06/$08.00+0 doi:10.1128/JVI.02014-05
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
C. T. T. Edwards,1,
N. McCarthy,2
J. Fox,3
H. Brown,1
A. Milicic,1
N. Mackie,3
T. Pillay,1
J. W. Drijfhout,4
S. Dustan,3
J. R. Clarke,3
E. C. Holmes,2
H. T. Zhang,1
K. Pfafferott,1
P. J. Goulder,1
M. O. McClure,3
J. Weber,3
R. E. Phillips,1* and
S. Fidler3
The James Martin 21st Century School, Peter Medawar Building for Pathogen Research, South Parks Road, Oxford OX1 3SY, United Kingdom,1 Department of Zoology, University of Oxford, Oxford OX1 3PS, United Kingdom,2 Division of Medicine, Wright Fleming Institute, Imperial College, St. Mary's Hospital, 4th Floor, Norfolk Place, Paddington, London W2 1PG, United Kingdom,3 Department of Immunohematology and Blood Transfusion, Leiden University Medical Center, P.O. Box 9600, 2300 RC Leiden, The Netherlands4
Received 21 September 2005/ Accepted 5 April 2006
|
|
|---|
|
|
|---|
The selection pressure imposed by HLA class I on HIV-1 results in viral escape mutations that facilitate evasion of the immune response (25, 26) and may compromise vaccine efficacy (4). Recently, it has been shown that certain of these mutations can be transmitted both sexually (2, 12) and vertically (18). At the population level, the overall extent to which different escape mutations are transmitted, and also their impact on the ability of the new host to prime an immune response, is unknown.
Studies of genetic host-pathogen interactions are limited by the requirement for large numbers of patients. In a study of over 400 patients with chronic HIV-1 infection, 64 associations were described between specific amino acids in the HIV-1 RT gene and HLA class I molecules (of which 12 remained statistically significant after correction for multiple comparisons) (21), potentially shaping HIV-1 protein variation and driving HIV-1 viral evolution (19). The extent to which this is true is dependent on the dynamics of transmission and the rates of mutation reversion or accumulation (11, 18). Here, we present a permutation-based statistical approach for testing global host-pathogen genetic associations in populations and, in conjunction with detailed genetic analysis of transmission between pairs of individuals, we show that variants of HIV-1 internal proteins, including those with escape properties, are being transmitted passively with no evidence of limiting selection.
|
|
|---|
Diagnosis of acute HIV infection. Criteria for acute HIV-1 infection were documented seronegative HIV-1 antibody test within the previous 6 months, acute symptomatic seroconversion illness, positive HIV-1 DNA PCR, or positive HIV-1 RNA quantification (Chiron 3.0; Chiron Inc., Emeryville, Calif.) in the absence of an antibody response. In order to confirm primary HIV-1 infection, sequential serum samples were later tested blindly with a "detuned" enzyme-linked immunosorbent assay (Abbott).
PCR amplification and sequencing. For the analysis of the acute and chronic cohorts, viral RNA was extracted from patient plasma and the HIV-1 RT gene was sequenced using the HIV-1 genotyping kit (version 2; ABI) according to the manufacturer's instructions. For analysis of the four transmission pairs, sequences were derived from viral RNA, except for pair 3, for which only proviral DNA was available. HIV-1 virion RNA encoding gag, pol, and nef was reverse transcribed, amplified, cloned, and sequenced. Viral RNA was extracted using the NucleoSpin RNA virus kit (Macherey-Nagel GmbH). Peripheral blood mononuclear cells (PBMC) were obtained from heparinized blood using Ficoll gradient centrifugation, and DNA was extracted from a minimum of 1 x 106 cells using the Purigene DNA extraction kit (Gentra) according to the manufacturer's instructions. Viral RNA was reverse transcribed with the Reverse-iT 1st Strand synthesis blender (ABgene, United Kingdom) and antisense primers, described previously (23). These primers were also used for PCR amplification. Amplification from both cDNA and proviral DNA was achieved with two rounds of nested PCR (95°C for 2 min, followed by 95°C for 30 s, 50°C for 30 s, and 75°C for 1.5 to 3 min for 35 cycles, and finally 75°C for 10 min). PCR products were gel purified (gel purification kit; QIAGEN), ligated into the pCR4-TOPO vector (Invitrogen), and sequenced with M13 forward and reverse primers using ABI BigDye Terminator (v. 3.0; Applied Biosystems). Trace files were checked and assembled using the Staden software package (www.mrc-lmb.cam.ac.uk/pubseq/staden_home.html). Between 10 and 48 clones were sequenced for each time point.
Phylogenetics. To illustrate the relationship between donor and recipient viruses in each of the transmission pairs, gag p24 nucleotide sequences obtained closest to the time of transmission were aligned manually using Se-Al and a neighbor-joining tree was constructed using PAUP*, assuming the HKY85 model of nucleotide substitution. The gag gene was used, as multiple clonal sequences for all eight members of the four transmission pairs were available, thereby facilitating a high-resolution phylogenetic comparison. The relationship between the four transmission pairs was confirmed using HIV-1 RT (data not shown). Subtype A and B reference sequences were obtained from the HIV Sequence Database (www.hiv.lanl.gov).
IFN-
ELISPOT assay.
Wells of sterile 96-well microtiter plates containing polyvinylidene difluoride (Millipore) were coated with mouse anti-human gamma interferon (IFN-
) monoclonal antibody 1-D1K (0.5 µg/ml; Mabtech). Frozen PBMC were added to the wells with peptide at a final concentration 105 to 109 M. Following overnight incubation, 0.5 µg/ml biotinylated anti-human IFN-
monoclonal antibody clone 7-B6-1 (Mabtech) was added. Streptavidine-conjugated alkaline phosphatase was added (0.5 µg/ml) and, following incubation for 40 min, spots were developed using chromogenic alkaline phosphatase substrate (Bio-Rad) and counted using the AID version 2.9 ELISPOT plate reader. The number of spots in the negative control was subtracted, and the results were normalized to give the number of spot-forming cells per 106 PBMC.
Binding assay. Thawed cryopreserved B-cell lines expressing the relevant HLA class I molecules were subjected to acid elution on ice using 0.263 M citric acid and 0.123 M Na2HPO4 in equal volumes. The cells were washed in RPMI-2% fetal calf serum and then incubated with fluorescein isothiocyanate-labeled high-affinity reference peptide (150 µM to 900 µM) and test peptide at threefold dilutions from 184 µM to 0.02 µM. Irrelevant peptides were included as negative controls. Following 24 h of incubation at 4°C, the cells were analyzed by flow cytometry. Percent inhibition of variant peptide against the fluorescent peptide was used as a measure of binding affinity.
Statistics. Two-tailed Spearman's nonparametric correlation coefficients were used to compare the frequencies of variant alleles in the different patient cohorts. In order to test for associations between HLA class I molecules and polymorphic amino acids, the distribution of the number of mutations corresponding to the associated HLA type that would be likely to occur by chance was estimated by permutation. The amino acid sequences of the patients were assigned randomly among the two-digit HLA class I types of the cases, and a score was calculated as the total number of times that a mutation corresponding to the associated HLA type occurred within the population. This was repeated 19,999 times to estimate the distribution of this score expected under the null hypothesis of no association, and the observed score was compared to this distribution. Fisher's exact test was used for the determination of specific HLA-amino acid associations. To investigate the determinants of transmission in the four transmitter pairs, the median frequencies of transmitted and nontransmitted donor group epitope sequences were compared using the nonparametric Wilcoxon matched pairs test.
|
|
|---|
![]() View larger version (15K): [in a new window] |
FIG. 1. Population-level distribution of amino acid variation in RT in two patient cohorts. Variant frequency in the HIV-1 RT gene (amino acids 1 to 250) across a population of 62 drug-naïve chronically infected patients with subtype B HIV-1 recruited at St. Mary's Hospital, London (a) and a population of 101 individuals sampled shortly after infection recruited at St. Mary's Hospital, London (b). (c) Differences in polymorphism frequencies for the HIV-1 RT gene in acute and chronic patients. *, site at which the difference in proportions is statistically significant.
|
![]() View larger version (22K): [in a new window] |
FIG. 2. Analyses of variant frequencies in acute versus chronic patients using Spearman's nonparametric correlation coefficient for all polymorphic sites in RT (a) or only sites at which variation is associated with the presence of the restricting HLA class I allele (b). Correlations are shown between acute seroconverters from St. Mary's Hospital and chronically infected patients from St. Mary's Hospital.
|
Furthermore, one can also model the impact of HLA on the viral population by correcting for attributable mutations, according to reported odds ratios (ORs) (21). We adjusted the prevalence of the 12 sites in RT most strongly correlated with HLA class I to reflect an absence of selection (i.e., if the variability at these sites was random). We found that at the 12 tested sites, 15 out of 129 mutations (13%; 95% confidence interval, 7% to 18%) could be attributed to HLA class I associations. If we assumed an absence of HLA-mediated influence, the removal of these polymorphisms in chronic infection resulted in a significant reduction in the global number of polymorphisms at these sites (P = 0.03 by Fisher's exact test), although not at the individual amino acids. The polymorphism distributions in acute and chronic infections were highly similar prior to this correction and, if the polymorphism distribution in acute infection can be shown to be random (i.e., no associations with HLA class I), then the HIV-1 signature evident in chronic infection is contributed to by HLA-induced selection pressure.
In acute HIV-1 infection, the statistical associations between viral polymorphisms and HLA class I molecules are lost. We, therefore, investigated polymorphisms in the acute seroconverters to see if associations with HLA class I molecules were present. The similarity between the frequency of HLA class I-adapted variants in the acute and chronic cohorts could reflect either the result of de novo adaptation of the virus to HLA in seroconverters or, alternatively, passive transmission of the HLA restriction pattern of the donor population. In the former case, associations between variation at a particular site and expression of a particular HLA class I allele in the recipient would be present in our acute infection cohort. To test this possibility, we again focused on the 12 amino acid-HLA class I associations in RT that have been identified in chronic infection (see Table S1 in the supplemental material) (21). For each acutely infected patient, we counted the number of times a variant was observed in association with expression of the appropriate HLA class I allele. We then summed the number of these observed associations for all individuals in the cohort to provide our test statistic, or score. For example, if each of the 12 associations was identified twice in the entire cohort, this would give a score of 2 x 12, or 24. The expected distribution of scores if there were no association between amino acid variation and HLA was estimated by randomly assigning the HLA allelic profiles of cases to the sequences in our data set and then recalculating the score. This process was repeated 19,999 times, creating a null distribution with a mean of 42.4, a standard deviation of 4.6, and an approximately Gaussian distribution (Fig. 3). The observed score for the seroconverter cohort was 46 (P = 0.38 for a two-tailed test). This analysis provides no evidence for a significant association between amino acid variation at these 12 selected sites and the HLA class I alleles expressed in acutely infected patients. This contrasts with the strong association that has been reported between HLA types and RT variation in chronically infected patients (21). If HLA associations of the strength reported in these chronically infected patients were present in our acute cohort, we would expect a score of 82.7, significantly different from the mean of the null distribution and the observed score of 46 (P < 0.0001). From this result, we conclude that the presence of HIV-1 variants in acute infection is not the result of de novo adaptation to the HLA alleles expressed by the new host.
![]() View larger version (15K): [in a new window] |
FIG. 3. HLA class I-independent variation in acutely infected patients. For each site known to adapt to the class I-restricted CTL response (see Table S1 in the supplemental material), the number of times a variant was observed in association with an HLA allele known to direct variation at that site was counted. The observed score for the seroconverter cohort was 46. Patient sequences were then randomly assigned to different HLA class I alleles expressed in the seroconverter cohort. Repetition of this process generated the null distribution of scores expected if the amino acid variants expressed by an infecting virus were unrelated to the HLA type of the patient. The mean for the null distribution was 42.4. The expected score for a cohort of previously described chronically infected patients (calculated using published odds ratios) (21) is 82.7, which is significantly different from the null distribution (P < 0.0001).
|
To examine this further we looked for variation in the 25 optimal epitopes previously documented in the fragment of the HIV-1 RT protein for which we had data (amino acids 29 to 225) (www.hiv.lanl.gov/content/immunology/maps/maps/html; see also Table S2 in the supplemental material). The second amino acid (position 2) and the carboxy-terminal amino acid are usually key to epitope binding in the HLA class I groove, and mutations at these sites frequently enable immune escape. We concentrated on these positions but also on any other intraepitopic sites at which there was significant variation. This analysis revealed four additional HLA class-amino acid associations (see Table S2). One of these was at the carboxy terminus of the HLA B*44-restricted EELRQHLLR epitope (aa 203 to 211). Interestingly, this association was only present in conjunction with HLA B*4403 and not other B*44 alleles. Three other associations at other amino acid sites were identified in the B*51-restricted epitope EI9 (aa 42 to 50), the B*35-restricted epitope VY10 (aa 118 to 127), and the B58-restricted epitope KL10 (aa 201 to 210) (see Table S2 in the supplemental material). Although these data support the previously published assertion that HLA class I is contributing to HIV-1 adaptation within human hosts, the extent of this selection appears less than that predicted from earlier reports (21).
In summary, the analysis of the population-level data in the acute and chronic cohorts suggests that transmission of HLA class I-restricted T-cell antigens is passive. In addition, statistical associations of HLA class I with specific viral polymorphisms are predominantly lost at acute seroconversion but become apparent in patients with chronic HIV-1 infection.
Viral diversity is decreased following transmission. To examine the passive transmission of HIV-1 internal protein variants at greater resolution, we studied four HIV-1-infected transmitter pairs. In these four pairs, virus samples were obtained from both donor and recipient within 2 to 6 weeks (mean, 3.8 weeks) of the time when acute infection was diagnosed in the index case. HIV-1 virion RNA encoding gag, pol, and nef was reverse transcribed, amplified, cloned, and sequenced. To compare the viruses circulating in each individual, a neighbor-joining tree was constructed using gag p24 nucleotide sequences (Fig. 4a). Reference sequences from the HIV Sequence Database (www.hiv.lanl.gov) were included for comparison. For subtype B, these reference sequences were selected from patients sampled in London between 1993 and 2003 and closely represent the genetic diversity prevalent in our study population. For subtype A, reference sequences were obtained from the Los Alamos HIV-1 Subtype Reference Alignments Database. For pair 1 the transmitted virus clustered with subtype A sequences, whereas for pairs 2, 3, and 4 the sequences conformed to subtype B. The close similarity between the virus circulating in the peripheral blood of each pair, relative to the database sequences and virus in other pairs, confirms the authenticity of the transmission pairs. Compared to each donor, shorter branch lengths leading to the recipient sequences illustrate a substantial loss of diversity (P < 0.0001, paired two-sample t test), consistent with a genetic bottleneck at the time of transmission. The average reduction in nucleotide diversity, determined from the number of substitutions per site, was 81.3% for gag p17, 71.2% for gag p24, and 91.0% for nef (see Table S3 in the supplemental material).
![]() View larger version (27K): [in a new window] |
FIG. 4. (a) Neighbor-joining phylogenetic tree of gag p24 nucleotide sequences from donor-recipient pairs, assuming the HKY85 model of nucleotide substitution. Sequences were 843 bp long. The number of sequences for each patient was as follows: D1, 48; R1, 43; D2, 19; R2, 14; D3, 29; R3, 25; D4, 38; R4, 27. Subtype A and B reference sequences were obtained from the HIV Sequence Database (www.hiv.lanl.gov). Horizontal branch lengths are drawn to scale. Divergence is the number of substitutions per site, expressed as a percentage. (b) Retention of HIV-1 amino acid variants during sexual transmission. gag p24 from four donor-recipient pairs is shown. The variant frequency (proportion of amino acids at a particular site that differ from the subtype consensus sequence) is plotted on the vertical axis. Different-colored bars represent different amino acids at a given site. The horizontal axis represents the linear amino acid sequence from position 1 to 278. Regions of gag p24 restricted by the HLA alleles expressed in the infected patient are marked by red bars.
|
Transmission of viral variants in HIV-1 internal proteins is frequency dependent with no evidence for selection. To investigate more closely the determinants of transmission in HIV-1 internal proteins, we focused on CTL antigens restricted by the HLA class I alleles of each donor. For each polymorphic site within known CTL epitopes in the HIV-1 gag, pol, and nef genes, variant frequency in the donor close to the time of transmission was plotted against whether the variant was transmitted (as judged by detection in the recipient). We found that transmission of a variant was dependent on its frequency in the donor (P < 0.0001, two-tailed Mann-Whitney test; n = 53) with high-frequency variants being more likely to be transmitted (Fig. 5).
![]() View larger version (13K): [in a new window] |
FIG. 5. Frequency-dependent transmission of HLA class I antigenic variants, as shown by donor frequencies of transmission of all intraepitope, polymorphic sites. The median frequencies for the transmitted and nontransmitted variants are shown and were found to be significantly different (P < 0.0001, two-tailed Mann-Whitney test).
|
production in an ELISPOT assay. We found that when only escape variants were analyzed, their frequency in the donor was the main determinant of transmission (P < 0.0001, two-tailed Mann-Whitney test). This analysis might be biased by the inevitable transmission of variants that had reached fixation in the donor. It was, therefore, repeated using a more stringent approach, selecting only escape variants that had not reached fixation (n = 24). In this analysis, transmission remained dependent on donor frequency (P = 0.0029, two-tailed Mann-Whitney test). This finding of frequency-dependent transmission indicates that any coexisting selection is insufficient to impede transmission of the majority of epitope variants from immunogenic proteins, such as HIV-1 Gag, Pol, and Nef. |
|
|---|
In a large survey we identified a pattern of polymorphic variation which was conserved across different patient cohorts and, specifically, at sites associated with CTL selection pressure. This polymorphic signature is more than a reflection of functional constraint. This is supported by the contrasting idiosyncratic distribution of polymorphisms within infected individuals, the loss of similarity between acute and chronic patterns on eliminating known HLA class I-associated sites in the latter, and the previously reported finding of a high proportion of sites in RT under HLA-driven selection pressure (21).
We found that the frequency of variants present in acutely infected patients was closely correlated with that in chronic patients but that the reported statistical associations of variant sites with particular HLA alleles was absent. These data suggest that the majority of variants observed in the acutely infected patients were transmitted by the donor and were unlikely to have arisen de novo as a result of selection pressure in the new host.
Nevertheless, the rate at which new viral adaptation takes place at different loci is unlikely to be uniform but will reflect sites where selection pressures are strongest. The presence, in acute infection, of the strong association between HLA-B51 and RT amino acid 135 suggests that CTL-induced selection pressure is acting on at least one intraepitope site. This amino acid lies at the C-terminal end of the TAFTIPSI epitope, which is a recognized site for escape mutation (27). A limitation of this statistical approach is that this strong signal in acute infection could also be the result of rapid reversion to wild type in non-HLA-B51 patients, as well as reflecting early escape. Such early emergence of immune escape is comparable with the speed with which simian immunodeficiency virus can adapt to major histocompatibility complex class I in macaques, where new escape variants have been detected within 4 weeks of infection (3). Although the differential rate of reversion of transmitted escape mutations due to fitness costs has not yet been fully defined, rapid reversion has been observed (11, 18).
We identified a number of amino acid polymorphisms in known epitopes at which variation was associated with the restricting HLA class I molecule. However, we might have expected to find a greater number of reported associations in our chronic cohort. This might be because our cohort is smaller (n = 128 versus n = 473) and, although we could test previously reported sites, we were not powered to screen the HIV-1 genome for new HLA class I associations. In addition, we assumed that the distribution of HLA class I molecules among HIV-infected individuals in Perthwhere the original study (21) was carried outand London would be similar, and this is not certain.
From our observations of four episodes of sexual transmission in homosexual men, in which samples were available in all parties close to the time of infection, we found the successful transfer of CTL epitope variants was dependent on their frequency in the donor. Cytotoxic T-lymphocyte escape may impair viral replicative fitness (11, 18); however, when a broad spectrum of variants was analyzed we saw no evidence these viruses have diminished capacity for transmission. This contrasts with drug resistance mutations, which may also reduce the replicative capacity of HIV-1 but appear to be less likely to be transmitted (6, 17, 28).
We conclude that during the transmission of HLA class I-restricted antigens derived from immunogenic internal viral proteins, dominant variants present in the donor are able to establish infection in the recipient. This has important implications for viral pathogenesis, as many of the transmitted variants will encode CTL escape properties. This study has only looked at the variability of the transmitted virus. The fate of transmitted variants and the extent to which mutants revert following transmission are preserved or are acquired de novo in response to the new host will determine the significance of the passive transmission of internal protein variants. In HIV-infected populations the characteristic signature pattern of amino acid variation in RT is remarkably stable, but each infected individual selects an idiosyncratic pattern of variants. This variability constitutes the viral "challenge" to which naïve partners are exposed and places an enormous burden on any protective vaccine.
The work was supported by the Medical Research Council (J.F.) and the Wellcome Trust (R.E.P., J.N.W., and P.J.G.), including the program grant SPARTAC. P.J.G. is supported by NIH grant AI046995-06.
We are indebted to the physicians and patients who contributed to the Spanish Swiss Intermittent Therapy Trial (SSITT).
Supplemental material for this article may be found at http://jvi.asm.org/. ![]()
A.J.F. and C.T.T.E. 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»