- Open Access
Unique features of HLA-mediated HIV evolution in a Mexican cohort: a comparative study
- Santiago Avila-Rios1, 2Email author,
- Christopher E Ormsby1,
- Jonathan M Carlson3,
- Humberto Valenzuela-Ponce1,
- Juan Blanco-Heredia1,
- Daniela Garrido-Rodriguez1,
- Claudia Garcia-Morales1,
- David Heckerman3,
- Zabrina L Brumme4, 5,
- Simon Mallal6,
- Mina John6,
- Enrique Espinosa1 and
- Gustavo Reyes-Teran1Email author
© Avila-Rios et al; licensee BioMed Central Ltd. 2009
- Received: 25 April 2009
- Accepted: 10 August 2009
- Published: 10 August 2009
Mounting evidence indicates that HLA-mediated HIV evolution follows highly stereotypic pathways that result in HLA-associated footprints in HIV at the population level. However, it is not known whether characteristic HLA frequency distributions in different populations have resulted in additional unique footprints.
The phylogenetic dependency network model was applied to assess HLA-mediated evolution in datasets of HIV pol sequences from free plasma viruses and peripheral blood mononuclear cell (PBMC)-integrated proviruses in an immunogenetically unique cohort of Mexican individuals. Our data were compared with data from the IHAC cohort, a large multi-center cohort of individuals from Canada, Australia and the USA.
Forty three different HLA-HIV codon associations representing 30 HLA-HIV codon pairs were observed in the Mexican cohort (q < 0.2). Strikingly, 23 (53%) of these associations differed from those observed in the well-powered IHAC cohort, strongly suggesting the existence of unique characteristics in HLA-mediated HIV evolution in the Mexican cohort. Furthermore, 17 of the 23 novel associations involved HLA alleles whose frequencies were not significantly different from those in IHAC, suggesting that their detection was not due to increased statistical power but to differences in patterns of epitope targeting. Interestingly, the consensus differed in four positions between the two cohorts and three of these positions could be explained by HLA-associated selection. Additionally, different HLA-HIV codon associations were seen when comparing HLA-mediated selection in plasma viruses and PBMC archived proviruses at the population level, with a significantly lower number of associations in the proviral dataset.
Our data support universal HLA-mediated HIV evolution at the population level, resulting in detectable HLA-associated footprints in the circulating virus. However, it also strongly suggests that unique genetic backgrounds in different HIV-infected populations may influence HIV evolution in a particular direction as particular HLA-HIV codon associations are determined by specific HLA frequency distributions. Our analysis also suggests a dynamic HLA-associated evolution in HIV with fewer HLA-HIV codon associations observed in the proviral compartment, which is likely enriched in early archived HIV sequences, compared to the plasma virus compartment. These results highlight the importance of comparative HIV evolutionary studies in immunologically different populations worldwide.
- Human Leukocyte Antigen
- Escape Mutation
- Human Leukocyte Antigen Allele
- Codon Pair
- Plasma Virus
The cytotoxic CD8+ T lymphocyte (CTL) response has been identified as an important selective pressure driving HIV evolution within an infected host [1–5]. Strong lines of evidence support the importance of the CTL response in HIV control, including the temporal correlation between the appearance of HIV-specific CTLs in vivo and the decline of viremia in the early stages of HIV infection , as well as the lack of control of virus levels after experimental depletion of CD8+ cells in rhesus macaques prior to simian immunodeficiency virus (SIV) infection . CTLs recognize and destroy infected cells through the binding of their T cell receptor (TCR) to viral peptides (epitopes) presented on the surface of infected cells by highly polymorphic molecules encoded by class I human leukocyte antigen (HLA) genes. Each HLA allele encodes a unique HLA molecule capable of presenting a broad range of possible epitopes derived from various areas of the HIV proteome. CTL recognition of these peptide-HLA complexes may be associated with different functional outcomes in the infection [8–10]. Importantly, as a result of CTL-mediated selective pressure, immune escape mutations are selected that hinder viral peptide binding to HLA molecules, prevent peptide processing before their presentation or lower TCR affinity of specific CTL clones to peptide-HLA complexes [4, 11–13]. Therefore, both the processes of antigen presentation to CTLs and CTL escape are HLA-restricted .
Depending on their costs to viral fitness, some CTL escape mutations can be transmitted and maintained in a new host [15–18], even without the presence of the originally selective HLA allele [11, 19–21]. Additionally, there is evidence supporting the notion that some immune escape mutations can accumulate in a large number of individuals and become fixed in the circulating virus consensus sequence, driving HIV evolutionary changes at the population level [9, 19, 22–24]. As a result, specific HLA epitopes could become extinct in the viral population, allowing HIV adaptation to HLA-associated immune control in a certain region [22, 23]. The relative impact of different factors that could influence the persistence of escape mutations in a large number of individuals remains incompletely understood [22, 25]. The variety of these factors–such as the extent of reversion of immune escape mutations in the absence of the selecting HLA allele [15, 16, 18], selection of compensatory mutations that restore viral fitness , founder effects , conflicting evolutionary forces on clustered epitopes , development of novel CTL responses to escape variants [28, 29], inter-clade differences in the circulating viruses [14, 30], immunodominance hierarchies of CTL responses [25, 31, 32], and HLA allele frequency distributions in different populations –highlight the complexity of viral adaptation to the immune response at the population level [9, 25].
In spite of this complexity, mounting evidence indicating that a large number of CTL escape mutations are reproducibly selected in the context of specific HLA restrictions has led to the hallmark observation that HIV evolution follows generally predictable mutational patterns in response to specific HLA-restricted immune responses (reviewed in ). This "HLA footprint effect" on HIV has been shown at the population level through correlative associations between the presence (or absence) of polymorphisms at specific positions of the viral sequence and the expression of specific HLA alleles [24, 25, 33–35]. Detection of HLA-HIV polymorphism associations is potentially limited by important confounding effects, namely HIV phylogeny, HIV codon covariation, and linkage disequilibrium of HLA alleles [10, 14]. Several studies have accounted for some of these confounding effects explicitly [24, 30, 34]; more recently, a comprehensive evolutionary model considering all these confounding sources was proposed . This phylogenetic dependency network model was shown to be able to reconstruct previously defined escape and compensatory mutation pathways and agrees with emerging data on patterns of epitope targeting. The existence of this kind of comprehensive models represents an opportunity to systematically study HIV evolution in immunogenetically different populations and assess the importance of different HLA backgrounds in HIV evolution at the population level. Due to the extensive polymorphism of HLA genes, allelic and haplotypic frequency distributions in distinct infected populations vary widely . Given the highly consistent effect of HLA-restricted selection on HIV evolution and the distinct HLA allele distributions in differing populations, it is likely that specific HLA-HIV polymorphism associations will be preferentially observed in different populations, determining unique characteristics of HIV evolution in different human groups [9, 36]. To explore this possibility, HLA-mediated HIV evolution at the population level was studied in a cohort of clade B-infected individuals from Central/Southern Mexico, and compared to previously reported studies in a large multicenter cohort of predominantly clade B-infected individuals from British Columbia, Canada; Western Australia; and the USA (the International HIV Adaptation Collaborative [IHAC] cohort) (Brumme ZL, John M, et al, PLoS ONE 2009, in press) [14, 25, 34, 37]. In order to determine to what extent HLA imprinting on HIV is a general phenomenon, it is informative to study HIV evolution in an immunogenetically unique population that possibly reflects a different selective pressure to that observed in other studied populations. The Mexican population is known to have a unique immunogenetic background characterized by the admixture of mainly Amerindian and Caucasian HLA haplotypes [38, 39]. To our knowledge, Latin American cohorts have not been the primary subject of HIV evolutionary studies. Our data suggest that the unique HLA frequency distribution in a previously uncharacterized, immunogenetically unique HIV-infected population is imprinting HIV evolution in a unique way. This fact underlines the importance of systematically expanding our understanding of CTL escape and HIV evolution in immunogenetically distinct populations. This knowledge has important implications for the design of CTL-based vaccines and treatment strategies.
Peripheral blood samples were prospectively obtained from 303 chronically-infected, HIV positive, antiretroviral treatment-naïve individuals from Central/Southern Mexico. Participating individuals were recruited with written informed consent at different health centers in Mexico City and from the states of Puebla, Jalisco, Oaxaca, Guerrero, the State of Mexico and Chiapas. Blood samples were shipped to and processed at the Center for Research in Infectious Diseases of the National Institute of Respiratory Diseases in Mexico City. All ethical issues related to this project were evaluated and approved by the Institutional Bioethics and Science Committee. For each patient, plasma aliquots and peripheral blood mononuclear cells (PBMCs) were obtained and cryopreserved.
HLA frequency and HLA-mediated HIV evolution data obtained from the Mexican cohort were compared with that obtained from a previously described cohort of 1,045 HIV-positive, predominantly Caucasian individuals from British Columbia, Canada (HOMER cohort)  and the large multicenter International HIV Adaptation Combined (IHAC) cohort, including 1,845 predominantly Caucasian individuals from British Columbia, Canada; Western Australia and the USA  (Brumme ZL, John M, et al, PLoS ONE 2009, in press).
Genomic DNA was extracted from at least 6 million PBMCs using QIAmp DNA Blood Mini Kit (QIAGEN, Valencia CA), according to the manufacturer's specifications. Class I HLA A, B and C genes were typed at low/medium resolution for each participating individual by sequence-specific primer polymerase chain reaction (SSP-PCR) using ABC SSP UniTray Kit (Invitrogen, Brown Deer, WI) according to the manufacturer's specifications. Briefly, genomic DNA from each participating individual at 75–125 ng/μL was used as template for 95 PCRs with different sequence-specific primers designed to detect relevant polymorphisms for typing. Reaction products were run on a 2.0% agarose gel (Promega, Madison, WI). Amplification patterns were analyzed with UniMatch v3.2 software using up-to-date data bases to determine HLA groups. All the reactions included an internal amplification control to be validated and each test included a reagent control to detect contamination.
HLA frequency analyses and comparisons
HLA allelic and population frequencies for the Mexican cohort were obtained with the HLA Frequency Analysis tool of the Los Alamos HIV Database http://www.hiv.lanl.gov/content/immunology/tools-links.html. HLA haplotype frequencies were obtained with the Arlequin v3.11 software. Due to the fact that the cohort was composed of non-related individuals with unknown family genetic backgrounds, a gametic phase estimation was carried out for each individual using a pseudo-Bayesian algorithm designed to reconstruct the gametic phase of multi-loci genotypes, included in the Arlequin v3.11 software (Excoffier-Laval-Balding, ELB) . Frequency analysis between the cohort reported here and the Canadian HOMER cohort, the multicenter IHAC cohort and a cohort of HIV-negative individuals from Central/Northern Mexico , was carried out by chi squared test, with post hoc two by two significance determined by Fisher's exact test, corrected for multiple comparisons by q values . Significant values were considered to be q < 0.05. These analyses were carried out with R statistical environment v2.8.1, using the package qvalue v1.1.
HIV polgenotyping from free plasma virus
Viral RNA from free plasma virus was purified from 1 mL of plasma using QIAmp Viral RNA Mini Kit (QIAGEN, Valencia, CA) according to the manufacturer's specifications. A fragment of the viral pol gene including the whole protease (PR) and 335 codons of the reverse transcriptase (RT) was bulk sequenced from plasma viral RNA for each participating individual. Sequences were obtained with a 3100-Avant Genetic Analyzer (Applied Biosystems, Foster City, CA), using ViroSeq HIV-1 Genotyping System (Celera Diagnostics, Alameda, CA) according to the manufacturer's specifications. Briefly, 1.3 Kbp fragments of the pol gene were amplified by RT-PCR from plasma viral RNA. PCR products were purified with ultra filtration columns and quantified in 1.5% agarose gels (Promega, Madison, WI). For each patient, sequencing PCRs were carried out with 7 different primers to assure that the whole genomic region was covered with at least two sequences. Sequences were assembled, aligned to the HXB2 consensus, and manually edited using the ViroSeq v2.7 software provided by the manufacturer.
HIV polgenotyping from PBMC proviral DNA
Genomic DNA was purified as described above. A fragment of approximately 1.5 Kbp covering the whole PR and the first 335 codons of RT was amplified by nested PCR with Platinum Taq DNA Polymerase (Invitrogen, Carlsbad, CA), and primers PR 5' OUTER 5'-CCCTAGGAAAAAGGGCTGTTG-3'/RT 3' OUTER 5'-GTTTTCAGATTTTTAAATGGCTCTTG-3', for the first round of amplification, and PR 5' INNER 5'-TGAAAGATTGTACTGAGAGACAGG-3'/RT 3' INNER 5'-GGCTCTTGATAAATTTGATATGTCC-3' for the second round of amplification. PCR conditions were 1 cycle of 94°C, for 3 min, followed by 35 cycles of 94°C for 30 s, 60°C for 30 s and 72°C for 2 min and a cycle of 72°C for 5 min, with final concentrations of 2 mM Mg++, 0.2 mM dNTPs, 0.4 mM of each primer and 20 ng/μL genomic DNA for both amplification rounds (transferring 10% of the volume of first round PCR product to the second round). In all cases, contamination controls were included. PCR products were purified by QIAquick PCR Purification Kit (QIAGEN, Valencia, CA) according to the manufacturer's specifications, and quantified in 2.0% agarose gels (Promega, Madison, WI). Seven sequencing PCRs were carried out for each patient using seven primer mixes included in the ViroSeq HIV-1 Genotyping System Kit (Celera Diagnostics, Alameda, CA), in order to cover the whole analyzed region with at least two sequences. Bulk proviral pol sequences were obtained with a 3100-Avant Genetic Analyzer (Applied Biosystems, Foster City, CA). Sequences were assembled, aligned to the HXB2 consensus, and edited manually using the ViroSeq v2.7 software.
The phylogenetic dependency network (PDN) model by Carlson, et al , was applied to infer patterns of CTL escape and codon covariation in the plasma and proviral sequence datasets, using the PhyloDv program http://www.codeplex.com/MSCompBio. The PDN model was designed to simultaneously account for HIV codon covariation, linkage disequilibrium among HLA alleles and the confounding effects of HIV phylogeny when attempting to identify HLA-associated polymorphisms in HIV . Briefly, the PDN model is a multivariate model that represents the probabilistic dependencies among a set of target attributes (in this case the presence or absence of amino acids at all codons in an HIV protein) and a set of predictor attributes (in this case the presence or absence of amino acids at all codons other than that for the target attribute in the HIV sequence, as well as the presence or absence of all possible HLA alleles) while correcting for the phylogenetic structure of the viral sequences. A dependency network graphically depicts which HLA and codon attributes predict each target codon attribute, associating a probability distribution for each target codon attribute, conditioned on various HLA and codon attributes. Importantly, each local probability distribution is corrected for the phylogenetic structure of the HIV sequences. To determine the significance of a particular predictor-target pair, the likelihood of a null model that reflects the assertion that the target variable is not under selection pressure from the predictor attribute is compared to an alternative model that reflects the assertion that the target variable is under selection pressure from that predictor attribute. Multiple predictors are added to the model in an iterative fashion using forward selection, in which the most significantly associated attribute is iteratively added to the model until no attribute achieves p < 0.05. The use of a multivariate model minimizes spurious associations caused by the presence of linkage disequilibrium among HLA alleles and HIV codon covariation. For each added predictor attribute, the most significant leaf distribution is recorded (escape, reversion, attraction, or repulsion, see below). The statistical significance of a predictor with respect to a target attribute is computed using false discovery rates (FDRs), which are computed using a likelihood-ratio test in which both the null and the alternative models are conditioned on all significant predictors that were identified in previous iterations of forward selection. For each p-value, we report the corresponding q-value, which is the minimum FDR among rejection regions that include that p-value, as computed using the method of Storey and Tibshirani with the π0 parameter conservatively set to one . Attributes were excluded as possible predictors when the corresponding predictor-target pair had a 2 × 2 contingency table in which any cell of the table had an observed or expected value of three or less.
The precise rules governing the transitions of the target attribute, conditioned on the predictor attributes and the sequence phylogeny, are given by a univariate leaf distribution, which is assumed to be the same for each individual. Four possible leaf distributions are defined: Attraction, having the predictor makes it more likely to have the target; Repulsion, not having the predictor makes it less likely to have the target; Escape, having the predictor makes it less likely to have the target; and Reversion, not having the predictor makes it more likely to have the target. The pair Attraction/Repulsion corresponds to a positive correlation between predictor and target, while the pair Escape/Reversion corresponds to a negative correlation between predictor and target.
General clinical and geographical characteristics of the Mexican cohort
Relevant clinical parameters for a cohort of 303 Mexican individuals.
CD4+ T cell count (cells/μL)*
CD4+ T cell %*
Viral Load (RNA copies/mL)
Log Viral Load
CD4+ T cell count stratification* [n (%)]
201 – 500 célls/μL
50 – 200 célls/μL
Gender [n (%)]
HLA allelic and haplotypic frequencies in a cohort of HIV-positive Mexican individuals
292 HIV-positive individuals from Central/Southern Mexico for whom class I HLA-A, B and C typing was available were used to characterize the immunogenetic background of this cohort. HLA allelic frequencies for the Mexican cohort are shown in the Additional file 1: Figure S1, Table S1. The most frequent alleles at the HLA-A locus were A*02, A*24, A*68 and A*31; the most frequent alleles at the HLA-B locus were B*39, B*35, B*40 and B*15; and the most frequent alleles at the HLA-C locus were Cw*07, Cw*04, Cw*03 and Cw*08 (Additional file 1: Figure S1). Characteristically, more than 60% of the participating individuals expressed A*02, more than 50% expressed Cw*07 and more than a third expressed B*39 and/or B*35 (Additional file 1: Figure S1, Table S1).
Most frequent two-gene class I HLA haplotypes in the cohort of HIV-positive Mexican individuals.†
To our knowledge, this study is the first formal report of class I HLA frequencies in a typical HIV-infected Mexican cohort.
Unique immunogenetic Background in a cohort of HIV-infected, antiretroviral treatment naïve individuals from Central/Southern Mexico
HLA-mediated HIV evolution in a Mexican cohort
HIV evolution mediated by HLA selection at the population level was studied using a 434 amino acid fragment spanning the whole HIV protease and 335 codons of the reverse transcriptase in 280 chronically-infected individuals from this cohort. The phylogenetic dependency network (PDN) model by Carlson et al , currently one of the most comprehensive models to assess HLA-mediated HIV evolution, was applied to infer patterns of CTL escape and codon co-variation in the Mexican cohort. Our results were compared with those previously derived from applying the PDN model to a thoroughly characterized, multi-center, combined cohort of predominantly clade B-infected, antiretroviral treatment-naïve individuals from British Columbia, Canada; Western Australia; and the USA (the IHAC cohort), with a clearly different immunogenetic background compared to the Mexican cohort [14, 34, 37, 46] (Figure 4). The Mexican cohort was also found to be predominantly clade B-infected (99.64%) with only one subtype other than B/recombinant form (CRF_06_cpx) identified (REGA HIV-1 Subtyping Tool 2.0, http://dbpartners.stanford.edu/RegaSubtyping/). A phylogenetic tree for the Mexican pol sequences included in this study is shown in the Additional file 1: Figure S2.
Interestingly, there were only two HIV pol sites previously associated with resistance to antiretroviral drugs that were also predicted to be associated with HLA selective pressure. B*18 was associated with an E to A change in RT position 138. The polymorphism 138A is associated with decreased response to non-nucleoside RT inhibitors (NNRTIs), including etravirine (Stanford University HIV Drug Resistance Database, http://hivdb.stanford.edu/). Similarly, Cw*07 was associated with a lower probability of having a D residue and a tendency for conservation of a V residue in RT position 179. The polymorphism 179D is associated with low level resistance to NNRTIs (Stanford University HIV Drug Resistance Database, http://hivdb.stanford.edu/). These observations show that HLA-mediated evolution can influence antiretroviral drug resistance, both promoting and preventing the presence of drug-resistance-related polymorphisms. This dual pressure phenomenon has been described previously [35, 47]; however, its frequency and population impact in the Mexican cohort will have to be assessed further.
HLA-HIV codon associations found in the Mexican cohort were compared to the ones previously found in a dataset 1845 HIV pol sequences from the combined IHAC cohort  (Brumme ZL, John M, et al, PLoS ONE 2009, in press) (Figure 6, Additional file 1: Table S3). Finding HLA-HIV codon associations in the smaller Mexican cohort that were not observed in the well-powered IHAC cohort could be indicative of unique HLA-driven HIV evolution in immunogenetically distinct cohorts. Not surprisingly, many of the observed HLA-HIV codon associations in the Mexican cohort were also predicted in the IHAC cohort, supporting the observation of highly conserved mutational patterns in HLA-driven HIV evolution (Figure 6, Additional file 1: Table S3). Nevertheless, important differences were also noted between the two cohorts. From the 43 HLA-HIV codon associations observed in the Mexican cohort, 23 were identified as novel associations, not previously observed in the larger IHAC cohort (nor in previous similar studies [25, 30, 34]), representing 18 different HLA-HIV codon pairs. Although several of these are likely to be false positives due to the 20% FDR, the fact that 53% of the present associations were not found in the well-powered IHAC cohort is striking. Furthermore, of these 18 new HLA-HIV codon pairs, 8 involved associations with codons that were not associated with any HLA allele in the IHAC cohort, and 4 involved associations with B*39, which is substantially more frequent in the Mexican cohort than in the IHAC cohort. Remarkably, 17 of the 23 novel associations involved HLA alleles whose frequencies were statistically indistinguishable from those in IHAC, suggesting that their presence is not due to increased statistical power, but rather may be due to differences in patterns of epitope targeting. In addition, 2 of the novel associations involved B*27 or B*08, two alleles that were significantly less frequent in the present cohort than in IHAC, which may reflect differences in epitope targeting between the cohorts or the fixation (and resulting drop in statistical power) of escape mutations in the IHAC cohort . Interestingly, although previously identified as HLA-associated in the IHAC cohort, some HLA-associated HIV codons in the Mexican cohort showed different HLA specificities and/or target amino acids. This was the case in 13 of the 23 novel HLA-HIV codon associations (10 of 18 HLA-HIV codon pairs). For example, Protease 71 V was associated with B*39 in the Mexican cohort, while it was associated with B*15 in the IHAC cohort; RT 245E was associated with B*18 in the Mexican cohort, but with B*57 in the IHAC cohort. Not surprisingly, B*39 was much more frequent in the Mexican cohort (p = 1.80E-44, q = 1.21E-42), while B*15 was more frequent in the IHAC cohort (p = 0.00964, q = 0.0231) (Figure 4). On the other hand B*57 was less frequent in the Mexican cohort (p = 0.000430, q = 0.00180) and no significant difference was found for B*18. These associations can be explained by experimentally confirmed or predicted epitopes both in the Mexican cohort and in the IHAC cohort (Figure 6; Los Alamos HIV Immunology Database, http://www.hiv.lanl.gov/content/immunology/index.html). Both B*15 and B*39 are predicted to have an epitope in position PR 68–76. These observations support the existence of sites in the HIV genome whose sequence variability at the population level reflects active selection pressure by different HLA alleles, and support previous observations that different HLA alleles may drive identical (as well as conflicting) escape mutations .
Point differences in consensus sequences between the Mexican and the HOMER cohorts.
Associated HLA Mexican cohort
Associated HLA HOMER cohort
HLA frequency Mexican cohort
HLA frequency HOMER cohort
As previously observed [22, 49], HLA-B alleles were involved in the majority of the associations (73% of HLA-HIV codon pairs), compared to HLA-A and C alleles in the Mexican cohort (p < 0.001 for both cases). Interestingly, 16.6% of HLA-HIV codon pairs were due to HLA-C alleles. Although in previous studies many of the associations apparently defined by HLA-C alleles represented indirect associations with HLA-B or A alleles due to the HLA linkage disequilibrium phenomenon, the PDN model used in the present study accounts for HLA linkage disequilibrium, minimizing the risk of finding these and other kinds of indirect associations. Nevertheless, it is important to note, that the ability of the PDN model to correct for HLA linkage disequilibrium is positively correlated with sample size and negatively correlated with the strength of linkage disequilibrium. Thus, false positive associations could still be found when strong linkage disequilibrium patterns exist, and random noise makes it difficult to distinguish the true associations. This could be the case for Cw*15-associated position PR 18, which could not be explained by epitope prediction. Cw*15 and B*51 are in strong linkage disequilibrium in the Mexican cohort (Additional file 1: Table S2), and an experimentally confirmed B*51 epitope exists that could explain the association. Although many of the HLA-C associations had high q-values, and could represent false positive associations, some of them were strongly associated and were in consonance with predicted and verified epitope mapping (Additional file 1: Table S3, Figure 6). These observations suggest an important role of HLA-C alleles in shaping HIV evolution at the population level in the Mexican cohort.
On the other hand, B*44 alone was responsible for 16.6% of HLA-HIV codon pairs, followed by B*51 and B*39, each responsible for 10% of the observed HLA-HIV codon pairs (Additional file 1: Table S3, Figure 6). This predominance is not observed in the IHAC cohort, possibly suggesting different patterns of immunodominance and HIV immune escape resulting from different epitope targeting between the two cohorts. These observations further support the existence of differential patterns of HIV selection by HLA alleles in populations worldwide.
Taken together, these results support the existence of highly conserved, universal HLA-mediated mutational patterns or "footprints" on HIV sequences at the population level. However, they also suggest that unique characteristics could exist in HLA-mediated HIV evolution in immunogenetically distinct populations, which can be detected even with cohorts of relatively small number of individuals.
Differences in HIV evolution between free plasma virus and PBMC proviral sequences
Protease and Reverse Transcriptase HLA-HIV codon associations in free plasma virus RNA and PBMC proviral DNA for the Mexican cohort.
Free plasma virus pol sequences†
PBMC proviral pol sequences†
Associated HIV Codon
Associated HIV codon
PR 77Vc, PR 36Md, PR 13Ia
PR 35Ec, RT 293Ia, RT 207Qd, PR 62Id
RT 123Dd, RT 121D d
B*51a, PR 12Td, RT 202Id
B*48b, RT 64Ka
PR 93Ld, PR 39Pa, PR 57Kd, Cw*15b
PR 57Ra, PR 93Ic, Cw*15d, RT 178Mb, PR 39Pa
PR 39Qb, PR 37Db, PR 37Eb
RT 272Pc, Cw*07c
RT 135Ia, RT 200A a
RT 135Ib, RT 200Ab
PR 14Rd, PR 67Cc
PR 64Lb, RT 165Tc, PR 91Tc, RT 135Va
PR 93Ld, PR 39Pa, PR 57Kd
Interestingly, the number of covarying HIV sites common to both compartments was lower than the number of covarying sites observed exclusively in plasma or in proviral sequences (p = 0.0489). Moreover, a large number of covarying sites was seen in proviral HIV sequences possibly reflecting remnants of viral adaptation to previous hosts.
Overall, differences in HLA-mediated selection were observed in the plasma virus and the PBMC provirus compartments, suggesting a highly dynamic HLA-associated evolution in HIV, as many of the HLA-HIV codon associations in the free plasma virus compartment were not evident in the proviral dataset, which likely contains early archived HIV sequences that appear to reflect less adaptation to within host HLA-mediated immune responses.
In this study we have presented evidence suggesting that a unique HLA allele frequency distribution in a cohort of clade B-infected Mexican individuals has left unique footprints on HIV sequences at the population level. We studied HLA-mediated HIV evolution in a clade B-infected Mexican cohort, comparing our data with data from the IHAC cohort, the largest clade B-infected cohort used to assess HLA-mediated evolution so far, which is composed of individuals from Canada, Australia, and the USA  (Brumme ZL, John M, et al, PLoS ONE 2009, in press). The two cohorts were shown to present notably different immunogenetic backgrounds, with an important admixture of Amerindian genes in the population of the Central/Southern part of Mexico (Figure 4). These different immunogenetic backgrounds provided a chance to assess the role of different HLA allele distributions in HLA-mediated selection in two cohorts infected by viruses of the same clade. The present cohort was shown to reflect the typical characteristics of an HIV-infected Mexican cohort, enriched in individuals in relatively advanced stages of HIV disease and presenting a similar HLA allele frequency distribution to the general population (Figure 3), with some specific exceptions (e. g. B*39) that will have to be assessed further in future studies.
Previous studies have suggested that HIV evolution at the population level follows broadly predictable, highly conserved mutational patterns associated with host CTL selective pressure [16, 23, 24, 33–35]. Conclusions obtained from a direct comparison between different studies in different populations have been limited mainly due to the use of different methods and models for assessing HLA-mediated viral evolution in which important sources of confounding are frequently not accounted. We applied the recently described PDN model , which simultaneously accounts for HLA linkage disequilibrium, HIV codon co-variation and viral lineage effects, to clade B pol sequences (Additional file 1: Figure S2) from the present cohort and compared the results to the immunogenetically distinct population of the IHAC cohort. Our data support the observations of highly conserved, universal, HLA-associated footprints in the HIV proteome at the population level, as many of the HLA – HIV codon associations found in the Mexican cohort have consistently been observed in the IHAC cohort, as well as in previous studies with diverse cohorts [24, 33–35, 37].
Interestingly, however, our data also suggest the existence of unique HLA-associated footprints in HIV, which could be influenced by specific HLA frequency distributions in different HIV-infected populations. The unique characteristics of HLA-mediated selection in the Mexican cohort was revealed not only by the presence of unique HLA-HIV codon pairs not detected in the IHAC cohort, but also by the presence of HIV positions previously identified as HLA associated, and with different HLA specificities and/or target amino acids in the two cohorts. The extent to which these unique HLA-associated footprints represent a real biological phenomenon and not a statistical effect will have to be further assessed with experimental data; nevertheless, evidence presented in this study strongly suggests the existence of real differences between the two cohorts. Although the Mexican cohort was much smaller than the IHAC cohort (the power to detect associations increases dramatically with sample size ), resulting in only 20% of the expected associations being confirmed in the present cohort, the fact that 53% of the HLA-HIV codon associations were novel in the Mexican cohort strongly suggests differences in HLA-mediated evolution between the two clade B-infected cohorts. Although these novel associations may represent false negatives from the IHAC cohort, that cohort is large enough that the false negative rate is expected to be quite small and any false negatives are likely to be rare events in it . It is also possible that the novel associations represent false positives in the present cohort; however, with an expected 20% false-positive rate due to the q < 0.2 threshold, the number of novel associations found in the Mexican cohort is striking. In addition, of the 18 novel HLA-HIV codon pairs in the Mexican cohort, at least 11 (61%) can be explained by confirmed or potential CTL epitopes (Figure 6), strongly arguing for the validity of these associations and for the existence of real biological differences in HLA-mediated selection between the two cohorts
The observation of point differences in the population consensus sequences of the two cohorts which were mapped to HLA-associated sites is a piece of evidence that further supports the differential impact of HLA selection in HIV evolution at the population level (Table 3). This was the case of position RT 277, associated with A*03 both in the Mexican cohort and in the IHAC cohort, in which the adapted form 277R has become fixed in the IHAC consensus while the susceptible form 277K has remained in the Mexican consensus. Not surprisingly, the frequency of A*03 was three times higher in the HOMER cohort than in the Mexican cohort (p = 7.08E-10, q = 1.00E-08), supporting an important role of HLA allele frequency in the fixation of HLA escape mutations at the population level (Figure 4, Table 3). Similarly, PR 93 was associated with B*15 in the IHAC cohort with the susceptible form 93I observed in this cohort's consensus, but the adapted form 93L observed in the Mexican consensus. Although no direct HLA association was detected in the Mexican cohort at this site (probably due to statistical power issues), position PR 93 was associated with other HIV sites, such as PR 71, which is HLA associated (Figure 6). Thus, changes in population consensus sequences may be linked with HLA-mediated selection. Also of interest is the observation that B*44 was associated with 5 of the 30 HLA-HIV codon pairs identified in the Mexican cohort. Two of these associations have been described in the IHAC cohort and two have been previously identified as HLA-associated positions with different HLA specificities (Figure 6). The strong influence of B*44 on HLA-mediated HIV evolution in the Mexican cohort could reflect differences in immunodominance hierarchies of CTL responses in the context of different HLA frequency distributions. Whereas strongly immunodominant CTL responses could be masking the effect of other less immunodominant responses in one cohort, these responses could have a greater impact on HLA-mediated HIV evolution in another cohort in which the immunodominant responses are infrequent. It is notable that the frequencies of many strongly immunodominant HLA alleles, such as B*57, B*27, B*08, B*07, A*03, and A*11 , are lower in the Mexican cohort compared to the IHAC cohort (q < 0.05) (Figure 4). It is possible that in the latter cohort, CTL responses restricted by these alleles could be masking the effect of other less immunogenic alleles that are frequently seen in the Mexican population. Indeed, this could be the case for Cw*07, the most frequent HLA-C allele group in the Mexican cohort, which explains 10% of HLA-HIV codon pairs observed in our analysis. These associations are unique to the Mexican cohort, and are supported by predicted epitopes and/or a strong statistical association (q < 0.05) (Figure 6).
The case of B*39 is also noteworthy, being the most frequent HLA-B allele group in the Mexican cohort with a frequency 7 times higher than that observed in the IHAC cohort (p = 1.80E-44, q = 1.21E-42) (Figure 4). B*39 explained another 10% of the HLA-HIV codon pairs identified in the Mexican cohort, suggesting either a strong influence of this allele in HIV evolution in the immunogenetic context of the Mexican cohort or a higher statistical power to detect associations. Interestingly, B*39-restricted associations were primarily escape associations (where possession of B*39 made it less likely to have the target amino acid in question), in which the target amino acid was a residue other than the consensus, suggesting that the consensus residue represents a possible escaped form for B*39 at this position (Additional file 1: Table S3, Figure 6). This could be suggestive of a frequent role of the B*39 allelic group in HIV codon conservation in the Mexican cohort. This HLA-associated conservation of sites has been previously described with highly frequent HLA alleles that promote the accumulation of CTL adapted variants in different populations [22, 23].
Overall, two general key aspects could explain the observation of different associations in cohorts that are infected by viruses of the same clade but which have different HLA frequencies: 1) Different patterns of immunodominance, which argue for real differences in CTL epitope targeting; and 2) Different statistical power to detect associations, which argues for a statistical effect rather than a biological difference. For example, the absence of strong immunodominance patterns in certain populations could potentially facilitate the detection of HIV polymorphisms associated with less immunodominant alleles. Being able to confirm this possibility at the population level strongly relies on low false positive/negative rates. Although the false negative rate on the IHAC data is low, further experimental data is necessary to confirm this point. On the other hand, different HLA frequencies can simply change the statistical power to detect associations, thus supporting the importance of assessing HLA-mediated selection in a diverse set of cohorts. The possibility also exists that a simple statistical power issue could be resolved by combining different cohorts infected by the same viral clade to make a larger reference set, supporting the creation of a universal set of associations that could get updated periodically as new sequences are added. Such a sequence and association database would allow extrapolation from a large reference set to new demographic groups for which collection of cohorts would be difficult. The fact that 17 of the 23 novel HLA-HIV codon associations in the Mexican cohort involved HLA alleles whose frequencies were not significantly different from those in the IHAC cohort strongly suggests that their presence is not due to increased statistical power but rather may be due to differences in patterns of epitope targeting. Furthermore, immunodominance effects as well as statistical power issues depending on HLA frequencies could both exist in the same dataset. Examples of both phenomena have been described above for the Mexican cohort, suggesting that a set of immunogenetically diverse cohorts could greatly enrich HIV evolutionary studies without the need of very large cohorts. It should also be noted that a broad two-digit HLA allele grouping does not reveal all possible divergence in HLA pressure, as a number of HLA subtypes with different peptide-binding motifs can be defined at four-digit level within some allelic groups such as B*35, B*40, B*51, B*58, A*02, all with highly characteristic distributions in different populations [9, 49]. Thus, significant divergence in selection in some cases could be explained by different dominant four-digit subtypes of the broad allele group in the compared cohorts. This fact could have an impact on statistical power to detect associations defined by different subtypes within a broad allele group in different populations and further argues for the unique HLA-associated imprinting of HIV in different populations.
In summary, although important limitations exist for the analysis of HLA-mediated HIV evolution in the Mexican population, including the presence of false positive associations and the low power to detect associations, our analysis yielded strong evidence suggesting that unique characteristics in HLA-mediated HIV evolution in the Mexican cohort indeed exist. These include the striking proportion of unique HLA-HIV codon associations in the Mexican cohort (many of which can be supported by predicted or confirmed CTL epitopes), the presence of HLA-associated differences in the consensus sequence with respect to the HOMER consensus (which reflects differential fixation of CTL escape mutations at the population level with a high dependency on HLA frequency), and the existence of a high proportion of novel associations that involve HLA alleles whose frequencies were similar in the Mexican and the IHAC cohorts (which argues against a statistical power issue in detecting at least some of the significant associations).
To further characterize HLA-mediated HIV evolution, HLA-HIV codon and HIV codon-HIV codon associations were compared in free plasma virus and PBMC proviral DNA in the cohort of Mexican individuals. As shown by graphically depicting the PDNs for the two viral compartments, different mutational patterns and different HLA-HIV codon associations were seen in actively replicating plasma viruses and PBMC-archived proviruses at the population level. A significantly lower number of HLA-HIV codon associations was observed in proviral sequences and there were more distinct than shared HIV codon-HIV codon associations in the two compartments (Figure 7). This could be explained by the observation that proviral sequences frequently represent a stable reservoir of HIV sequences archived early in the course of the infection , whereas plasma viruses represent sequences from later in the course of infection. Thus, the proviral sequences may have been archived before some epitopes were targeted by host CTL responses, or before escape mutations had a chance of being selected at epitopes already being targeted by CTLs, resulting in fewer associations in proviral sequences than in the extant plasma sequences. Indeed, previous studies have shown the presence of HLA-associated escape mutations in plasma viruses that are rare in proviruses within infected individuals . Nevertheless, proviral HLA-HIV codon pairs could not be mapped to known epitopes of early escape  in the present data, although the possibility exists that a larger cohort and analyses in other viral genes could further support this correlation. However, given the differences in escape association that we have observed between the Mexican and IHAC cohorts and the observation that the cohort described  is immunogenetically similar to the IHAC cohort, it may be that the discordance between proviral escape associations reported here and previously reported early-escape epitopes reflects different patterns of CTL epitope targeting and kinetics between the two populations. The proviral associations in the Mexican cohort could thus represent early escape events in a Latin American cohort setting.
Surprisingly, some HLA associations detected for proviral sequences were not seen in the plasma virus dataset. Some of these HLA-HIV codon pairs observed exclusively in proviral sequences have fairly high q-values, possibly suggesting the presence of false positive associations. However, unique proviral associations could also suggest a chronological reshaping of HLA-mediated HIV evolution, reflecting rapidly reverting mutations which are lost soon after transmission to HLA-mismatched individuals. Alternatively, the existence of organ compartmentalization of HIV variants within an infected host and its relation to positive selection has been described . This phenomenon could explain population differences between actively replicating viruses coming from a specific compartment with characteristic selective pressures and archived proviruses, remaining as reservoir(s) originating from different anatomical and/or cellular compartments.
Shared associations between the plasma virus and the PBMC provirus compartments may reflect sites in the viral proteome with continuous CTL targeting throughout the chronic infection, a characteristic that might be of interest in the selection of candidate vaccine targets. On the other hand, these apparently more stable associations could also reflect epitopes with early CTL targeting that has stopped, but for which no reversion has occurred, suggesting low fitness costs for escape. If the latter case were true, some shared associations might be more likely to reach fixation at the population level in the future. This would have implications for our understanding and predictive capabilities of HIV adaptation in human populations.
Similarly, unique coevolving HIV codon pairs were detected in proviral sequences and in plasma virus sequences, perhaps reflecting different patterns of compensatory mutations to the different HLA escape mutations observed in the two compartments. Alternatively, unique proviral HIV codon-HIV codon pairs could be explained as a reorganization of mutational patterns in HIV evolution that reflect escape mutations selected in previous hosts as well as new mutations selected in the current host, while unique plasma virus HIV codon-HIV codon pairs could reflect sequential footprints left by viral adaptation to HLA-restricted responses in chronic infection in the current host. These observations bring up interesting consequences for our understanding of HLA-mediated HIV evolution, suggesting that the appearance and density of the PDNs for a specific population are highly dynamic and could vary in time. The dynamic development of CTL responses over the course of infection within an individual has been previously reported [51, 52]. Further studies in follow-up cohorts or in carefully stratified cross-sectional cohorts might be able to support or refute these observations.
In conclusion, our data derived from analysis of HLA-mediated HIV evolution in a previously uncharacterized, immunogenetically unique cohort from Central/Southern Mexico, support a highly conserved and strongly predictable component of HLA-mediated HIV evolution at the population level, resulting in HLA-associated footprints in the circulating virus worldwide. This effect is evident even after considering important objections to the HLA-HIV population imprinting hypothesis, such as the rapid reversion of a considerable part of the total CTL escape mutations in the absence of the selective HLA allele , the complexity of the CTL response which frequently imposes conflicting selective forces in the same site of the viral sequence , the possibility that an escape variant selected by a specific HLA allele can be targeted by CTL responses restricted by different HLA alleles [18, 30, 57], and the dynamic immunodominance hierarchies observed in HLA-restricted responses [51, 52]. Interestingly, the HLA-mediated evolution analysis in our cohort of Mexican individuals showed additional HLA-HIV codon associations that have not been described in previously studied cohorts, including the large multi-center IHAC cohort with a clearly different immunogenetic background to the Mexican cohort. This fact supports the possibility that these specific associations are not significantly impacting HIV evolution at the population level in other cohorts, but that they are significant in the immunogenetic context of the Mexican population. Comparative HLA-mediated HIV evolution studies, with comparable methods that take into account important confounding factors such as HIV codon covariation, HIV lineage effects and HLA linkage disequilibrium, can thus be useful in identifying these distinct HLA-associated footprints in different populations worldwide. Extending such comparative studies to other immunogenetically distinct cohorts, would allow the reconstruction of a more complete panorama of the impact of HLA selection in HIV evolution worldwide. This knowledge may prove useful for the development of vaccine candidates and the development of therapeutic strategies directed to specific populations. The creation of a universal database of HLA-associated HIV sites applying comprehensive and comparable models to assess HLA-mediated evolution in immunogenetically divergent cohorts from different parts of the world, including cohorts predominantly infected by different viral clades, would greatly improve our understanding of HIV evolution worldwide. Importantly, further experimental evidence will help to understand the limitations imposed by statistical models to detect footprints of HLA-associated evolution in HIV in different populations.
Additionally, a comparison between HLA-mediated evolution in free plasma virus and PBMC proviral sequences suggested a highly dynamic HLA-associated evolution in HIV, as many of the HLA-HIV codon associations observed in the free plasma virus compartment are not evident in the proviral dataset, which is presumably enriched in early HIV sequences and does not reflect the full extent of within-host HLA-driven viral evolution. Moreover, shared HLA-HIV codon associations in both viral compartments could be of interest, reflecting epitopes with continuous CTL targeting throughout the chronic infection or, alternatively, escape mutations with low fitness costs that could reach fixation at the population level in the future. Further studies with larger cohorts and various viral genes could enrich these primary observations and increase our understanding of HIV adaptation to different populations worldwide.
We thank all patients of the Mexican cohort for their participation in this study; the physicians Akio Murakami, María Gomez-Palacio, José L. Sandoval, Daniela de la Rosa, Jorge Ibarra, Ricardo S. Vega and Cristina Sánchez of the Center for Research in Infectious Diseases of the National Institute of Respiratory Diseases in Mexico City, for their help in recruiting patients; Carolina Demeneghi, Mario Preciado, and Silvia del Arenal, for collection of blood samples; Ramón Hernández, and Verónica Quiroz, for viral load, and HIV genotyping assays; Edna Rodríguez for CD4+ T cell count assays; Zeidy Arenas, Sandra Zamora, Rosalinda Hernández, Eduardo López, for their administrative support; Dr. Joel Vázquez, for his technical guidance; Dr. Luis Padilla-Noriega and Dr. Eduardo García-Zepeda for their academic counselling. We thank Dr. Indiana Torres, Dr. Beatriz Ramírez, Dr. Adakatia Armenta, Dr. Jaime Andrade and Dr. Lucero González for providing blood samples of individuals from different states in Mexico. We thank the Doctorate in Biomedical Sciences Program and the National Autonomous University of Mexico, for their support. We thank Dr. Bruce Walker for his mentorship and support, and the BC Centre for Excellence in HIV/AIDS for providing access to data. We thank Dr. Richard Harrigan and the BC Centre for Excellence in HIV/AIDS for providing access to data on the HOMER cohort.
Sponsorship: This work was supported by grants from CONACYT U48159M, Fundación México Vivo, and Comisión de Equidad y Género de la H. Cámara de Diputados. SAR was supported by a scholarship from CONACyT. ZLB is supported by a post-doctoral fellowship from the Canadian Institutes for Health Research (CIHR).
- Allen TM, Altfeld M, Geer SC, Kalife ET, Moore C, O'Sullivan KM, Desouza I, Feeney ME, Eldridge RL, Maier EL, et al: Selective escape from CD8+ T-cell responses represents a major driving force of human immunodeficiency virus type 1 (HIV-1) sequence diversity and reveals constraints on HIV-1 evolution. J Virol. 2005, 79: 13239-13249. 10.1128/JVI.79.21.13239-13249.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Blankson JN, Bailey JR, Siliciano RF: Crosscurrents in HIV-1 evolution. Nat Immunol. 2006, 7: 121-122. 10.1038/ni0206-121.View ArticlePubMedGoogle Scholar
- Borrow P, Lewicki H, Wei X, Horwitz MS, Peffer N, Meyers H, Nelson JA, Gairin JE, Hahn BH, Oldstone MB, Shaw GM: Antiviral pressure exerted by HIV-1-specific cytotoxic T lymphocytes (CTLs) during primary infection demonstrated by rapid selection of CTL escape virus. Nat Med. 1997, 3: 205-211. 10.1038/nm0297-205.View ArticlePubMedGoogle Scholar
- Goulder PJ, Watkins DI: HIV and SIV CTL escape: implications for vaccine design. Nat Rev Immunol. 2004, 4: 630-640. 10.1038/nri1417.View ArticlePubMedGoogle Scholar
- Price DA, Goulder PJ, Klenerman P, Sewell AK, Easterbrook PJ, Troop M, Bangham CR, Phillips RE: Positive selection of HIV-1 cytotoxic T lymphocyte escape variants during primary infection. Proc Natl Acad Sci USA. 1997, 94: 1890-1895. 10.1073/pnas.94.5.1890.PubMed CentralView ArticlePubMedGoogle Scholar
- Koup RA, Safrit JT, Cao Y, Andrews CA, McLeod G, Borkowsky W, Farthing C, Ho DD: Temporal association of cellular immune responses with the initial control of viremia in primary human immunodeficiency virus type 1 syndrome. J Virol. 1994, 68: 4650-4655.PubMed CentralPubMedGoogle Scholar
- Schmitz JE, Kuroda MJ, Santra S, Sasseville VG, Simon MA, Lifton MA, Racz P, Tenner-Racz K, Dalesandro M, Scallon BJ, et al: Control of viremia in simian immunodeficiency virus infection by CD8+ lymphocytes. Science. 1999, 283: 857-860. 10.1126/science.283.5403.857.View ArticlePubMedGoogle Scholar
- Kiepiela P, Ngumbela K, Thobakgale C, Ramduth D, Honeyborne I, Moodley E, Reddy S, de Pierres C, Mncube Z, Mkhwanazi N, et al: CD8+ T-cell responses to different HIV proteins have discordant associations with viral load. Nat Med. 2007, 13: 46-53. 10.1038/nm1520.View ArticlePubMedGoogle Scholar
- Goulder PJ, Watkins DI: Impact of MHC class I diversity on immune control of immunodeficiency virus replication. Nat Rev Immunol. 2008, 8: 619-630. 10.1038/nri2357.PubMed CentralView ArticlePubMedGoogle Scholar
- Carlson JM, Brumme ZL: HIV evolution in response to HLA-restricted CTL selection pressures: a population-based perspective. Microbes Infect. 2008, 10: 455-461. 10.1016/j.micinf.2008.01.013.View ArticlePubMedGoogle Scholar
- Goulder PJ, Brander C, Tang Y, Tremblay C, Colbert RA, Addo MM, Rosenberg ES, Nguyen T, Allen R, Trocha A, et al: Evolution and transmission of stable CTL escape mutations in HIV infection. Nature. 2001, 412: 334-338. 10.1038/35085576.View ArticlePubMedGoogle Scholar
- Phillips RE, Rowland-Jones S, Nixon DF, Gotch FM, Edwards JP, Ogunlesi AO, Elvin JG, Rothbard JA, Bangham CR, Rizza CR, et al: Human immunodeficiency virus genetic variation that can escape cytotoxic T cell recognition. Nature. 1991, 354: 453-459. 10.1038/354453a0.View ArticlePubMedGoogle Scholar
- Draenert R, Le Gall S, Pfafferott KJ, Leslie AJ, Chetty P, Brander C, Holmes EC, Chang SC, Feeney ME, Addo MM, et al: Immune selection for altered antigen processing leads to cytotoxic T lymphocyte escape in chronic HIV-1 infection. J Exp Med. 2004, 199: 905-915. 10.1084/jem.20031982.PubMed CentralView ArticlePubMedGoogle Scholar
- Carlson JM, Brumme ZL, Rousseau CM, Brumme CJ, Matthews P, Kadie C, Mullins JI, Walker BD, Harrigan PR, Goulder PJ, Heckerman D: Phylogenetic dependency networks: inferring patterns of CTL escape and codon covariation in HIV-1 Gag. PLoS Comput Biol. 2008, 4: e1000225-10.1371/journal.pcbi.1000225.PubMed CentralView ArticlePubMedGoogle Scholar
- Friedrich TC, Dodds EJ, Yant LJ, Vojnov L, Rudersdorf R, Cullen C, Evans DT, Desrosiers RC, Mothe BR, Sidney J, et al: Reversion of CTL escape-variant immunodeficiency viruses in vivo. Nat Med. 2004, 10: 275-281. 10.1038/nm998.View ArticlePubMedGoogle Scholar
- Leslie AJ, Pfafferott KJ, Chetty P, Draenert R, Addo MM, Feeney M, Tang Y, Holmes EC, Allen T, Prado JG, et al: HIV evolution: CTL escape mutation and reversion after transmission. Nat Med. 2004, 10: 282-289. 10.1038/nm992.View ArticlePubMedGoogle Scholar
- Peyerl FW, Bazick HS, Newberg MH, Barouch DH, Sodroski J, Letvin NL: Fitness costs limit viral escape from cytotoxic T lymphocytes at a structurally constrained epitope. J Virol. 2004, 78: 13901-13910. 10.1128/JVI.78.24.13901-13910.2004.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu Y, McNevin J, Zhao H, Tebit DM, Troyer RM, McSweyn M, Ghosh AK, Shriner D, Arts EJ, McElrath MJ, Mullins JI: Evolution of human immunodeficiency virus type 1 cytotoxic T-lymphocyte epitopes: fitness-balanced escape. J Virol. 2007, 81: 12179-12188. 10.1128/JVI.01277-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Edwards CT, Pfafferott KJ, Goulder PJ, Phillips RE, Holmes EC: Intrapatient escape in the A*0201-restricted epitope SLYNTVATL drives evolution of human immunodeficiency virus type 1 at the population level. J Virol. 2005, 79: 9363-9366. 10.1128/JVI.79.14.9363-9366.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Furutsuki T, Hosoya N, Kawana-Tachikawa A, Tomizawa M, Odawara T, Goto M, Kitamura Y, Nakamura T, Kelleher AD, Cooper DA, Iwamoto A: Frequent transmission of cytotoxic-T-lymphocyte escape mutants of human immunodeficiency virus type 1 in the highly HLA-A24-positive Japanese population. J Virol. 2004, 78: 8437-8445. 10.1128/JVI.78.16.8437-8445.2004.PubMed CentralView ArticlePubMedGoogle Scholar
- Piontkivska H, Hughes AL: Patterns of sequence evolution at epitopes for host antibodies and cytotoxic T-lymphocytes in human immunodeficiency virus type 1. Virus Res. 2006, 116: 98-105. 10.1016/j.virusres.2005.09.001.View ArticlePubMedGoogle Scholar
- Kawashima Y, Pfafferott K, Frater J, Matthews P, Payne R, Addo M, Gatanaga H, Fujiwara M, Hachiya A, Koizumi H, et al: Adaptation of HIV-1 to human leukocyte antigen class I. Nature. 2009, 458: 641-5. 10.1038/nature07746.PubMed CentralView ArticlePubMedGoogle Scholar
- Leslie A, Kavanagh D, Honeyborne I, Pfafferott K, Edwards C, Pillay T, Hilton L, Thobakgale C, Ramduth D, Draenert R, et al: Transmission and accumulation of CTL escape variants drive negative associations between HIV polymorphisms and HLA. J Exp Med. 2005, 201: 891-902. 10.1084/jem.20041455.PubMed CentralView ArticlePubMedGoogle Scholar
- Moore CB, John M, James IR, Christiansen FT, Witt CS, Mallal SA: Evidence of HIV-1 adaptation to HLA-restricted immune responses at a population level. Science. 2002, 296: 1439-1443. 10.1126/science.1069660.View ArticlePubMedGoogle Scholar
- Brumme ZL, Brumme CJ, Carlson J, Streeck H, John M, Eichbaum Q, Block BL, Baker B, Kadie C, Markowitz M, et al: Marked epitope and allele-specific differences in rates of mutation in HIV-1 Gag, Pol and Nef CTL epitopes in acute/early HIV-1 infection. J Virol. 2008, 82: 9216-27. 10.1128/JVI.01041-08.PubMed CentralView ArticlePubMedGoogle Scholar
- Schneidewind A, Brockman MA, Yang R, Adam RI, Li B, Le Gall S, Rinaldo CR, Craggs SL, Allgaier RL, Power KA, et al: Escape from the dominant HLA-B27-restricted cytotoxic T-lymphocyte response in Gag is associated with a dramatic reduction in human immunodeficiency virus type 1 replication. J Virol. 2007, 81: 12382-12393. 10.1128/JVI.01543-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Iversen AK, Stewart-Jones G, Learn GH, Christie N, Sylvester-Hviid C, Armitage AE, Kaul R, Beattie T, Lee JK, Li Y, et al: Conflicting selective forces affect T cell receptor contacts in an immunodominant human immunodeficiency virus epitope. Nat Immunol. 2006, 7: 179-189. 10.1038/ni1298.View ArticlePubMedGoogle Scholar
- Mason RD, Grant MD: A therapy-related point mutation changes the HLA restriction of an HIV-1 Pol epitope from A2 to B57 and enhances its recognition. Aids. 2005, 19: 981-984. 10.1097/01.aids.0000171415.07034.0d.View ArticlePubMedGoogle Scholar
- Allen TM, Yu XG, Kalife ET, Reyor LL, Lichterfeld M, John M, Cheng M, Allgaier RL, Mui S, Frahm N, et al: De novo generation of escape variant-specific CD8+ T-cell responses following cytotoxic T-lymphocyte escape in chronic human immunodeficiency virus type 1 infection. J Virol. 2005, 79: 12952-12960. 10.1128/JVI.79.20.12952-12960.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Bhattacharya T, Daniels M, Heckerman D, Foley B, Frahm N, Kadie C, Carlson J, Yusim K, McMahon B, Gaschen B, et al: Founder effects in the assessment of HIV polymorphisms and HLA allele associations. Science. 2007, 315: 1583-1586. 10.1126/science.1131528.View ArticlePubMedGoogle Scholar
- Altfeld M, Kalife ET, Qi Y, Streeck H, Lichterfeld M, Johnston MN, Burgett N, Swartz ME, Yang A, Alter G, et al: HLA Alleles Associated with Delayed Progression to AIDS Contribute Strongly to the Initial CD8(+) T Cell Response against HIV-1. PLoS Med. 2006, 3: e403-10.1371/journal.pmed.0030403.PubMed CentralView ArticlePubMedGoogle Scholar
- Streeck H, Lichterfeld M, Alter G, Meier A, Teigen N, Yassine-Diab B, Sidhu HK, Little S, Kelleher A, Routy JP, et al: Recognition of a defined region within p24 gag by CD8+ T cells during primary human immunodeficiency virus type 1 infection in individuals expressing protective HLA class I alleles. J Virol. 2007, 81: 7725-7731. 10.1128/JVI.00708-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Ahlenstiel G, Roomp K, Daumer M, Nattermann J, Vogel M, Rockstroh JK, Beerenwinkel N, Kaiser R, Nischalke HD, Sauerbruch T, et al: Selective pressures of HLA genotypes and antiviral therapy on human immunodeficiency virus type 1 sequence mutation at a population level. Clin Vaccine Immunol. 2007, 14: 1266-1273. 10.1128/CVI.00169-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Brumme ZL, Brumme CJ, Heckerman D, Korber BT, Daniels M, Carlson J, Kadie C, Bhattacharya T, Chui C, Szinger J, et al: Evidence of Differential HLA Class I-Mediated Viral Evolution in Functional and Accessory/Regulatory Genes of HIV-1. PLoS Pathog. 2007, 3: e94-10.1371/journal.ppat.0030094.PubMed CentralView ArticlePubMedGoogle Scholar
- Mueller SM, Schaetz B, Eismann K, Bergmann S, Bauerle M, Schmitt-Haendle M, Walter H, Schmidt B, Korn K, Sticht H, et al: Dual selection pressure by drugs and HLA class I-restricted immune responses on human immunodeficiency virus type 1 protease. J Virol. 2007, 81: 2887-2898. 10.1128/JVI.01547-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Pond SL, Frost SD, Grossman Z, Gravenor MB, Richman DD, Brown AJ: Adaptation to different human populations by HIV-1 revealed by codon-based analyses. PLoS Comput Biol. 2006, 2: e62-10.1371/journal.pcbi.0020062.View ArticlePubMedGoogle Scholar
- Brumme ZL, John M, Brumme CJ, Carlson J, Haubrich R, Riddler S, Swenson L, Tao I, Szeto S, Chan D, et al: Relationship between HLA class I-driven evolution in Gag, Pol and Nef and clinical markers of HIV disease: a multi-center collaborative study. AIDS Vaccine, Abstract P09-01; Cape Town, South Africa. 2008Google Scholar
- Barquera R, Zuniga J, Hernandez-Diaz R, Acuna-Alonzo V, Montoya-Gama K, Moscoso J, Torres-Garcia D, Garcia-Salas C, Silva B, Cruz-Robles D, et al: HLA class I and class II haplotypes in admixed families from several regions of Mexico. Mol Immunol. 2008, 45: 1171-1178. 10.1016/j.molimm.2007.07.042.View ArticlePubMedGoogle Scholar
- Gorodezky C, Alaez C, Vazquez-Garcia MN, de la Rosa G, Infante E, Balladares S, Toribio R, Perez-Luque E, Munoz L: The genetic structure of Mexican Mestizos of different locations: tracking back their origins through MHC genes, blood group systems, and microsatellites. Hum Immunol. 2001, 62: 979-991. 10.1016/S0198-8859(01)00296-8.View ArticlePubMedGoogle Scholar
- Excoffier L, Laval G, Balding D: Gametic phase estimation over large genomic regions using an adaptive window approach. Hum Genomics. 2003, 1: 7-19.PubMed CentralView ArticlePubMedGoogle Scholar
- Storey JD, Tibshirani R: Statistical significance for genomewide studies. Proc Natl Acad Sci USA. 2003, 100: 9440-9445. 10.1073/pnas.1530509100.PubMed CentralView ArticlePubMedGoogle Scholar
- Soto-Ramirez LE: HIV/AIDS in Latin America. Science. 2008, 321: 465-10.1126/science.1162896.View ArticlePubMedGoogle Scholar
- Bastos FI, Caceres C, Galvao J, Veras MA, Castilho EA: AIDS in Latin America: assessing the current status of the epidemic and the ongoing response. Int J Epidemiol. 2008, 37: 729-737.View ArticlePubMedGoogle Scholar
- Magis-Rodríguez C, Hernández-Avila M: Epidemiología del SIDA en México. 25 Años de SIDA en México Logros, desaciertos y retos. Edited by: Córdova-Villalobos J, Ponce de León-Rosales S, Valdespino J. 2008, Cuernavaca, Morelos, México: Instituto Nacional de Salud Pública, 85-103. FirstGoogle Scholar
- Llop E: El sistema mayor de histocompatibilidad humano, HLA. Poblaciones chilenas: cuatro décadas de investigaciones bioantropológicas. Edited by: Rothhammer F, Llop E. 2004, Chile: Editorial Universitaria, 165-188.Google Scholar
- Brumme ZL, Tao I, Szeto S, Brumme CJ, Carlson JM, Chan D, Kadie C, Frahm N, Brander C, Walker B, et al: Human leukocyte antigen-specific polymorphisms in HIV-1 Gag and their association with viral load in chronic untreated infection. Aids. 2008, 22: 1277-1286. 10.1097/QAD.0b013e3283021a8c.View ArticlePubMedGoogle Scholar
- Karlsson AC, Deeks SG, Barbour JD, Heiken BD, Younger SR, Hoh R, Lane M, Sallberg M, Ortiz GM, Demarest JF, et al: Dual pressure from antiretroviral therapy and cell-mediated immune response on the human immunodeficiency virus type 1 protease gene. J Virol. 2003, 77: 6743-6752. 10.1128/JVI.77.12.6743-6752.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- Matthews PC, Prendergast A, Leslie A, Crawford H, Payne R, Rousseau C, Rolland M, Honeyborne I, Carlson J, Kadie C, et al: Central role of reverting mutations in HLA associations with human immunodeficiency virus set point. J Virol. 2008, 82: 8548-8559. 10.1128/JVI.00580-08.PubMed CentralView ArticlePubMedGoogle Scholar
- Kiepiela P, Leslie AJ, Honeyborne I, Ramduth D, Thobakgale C, Chetty S, Rathnavalu P, Moore C, Pfafferott KJ, Hilton L, et al: Dominant influence of HLA-B in mediating the potential co-evolution of HIV and HLA. Nature. 2004, 432: 769-775. 10.1038/nature03113.View ArticlePubMedGoogle Scholar
- Bailey JR, Brennan TP, O'Connell KA, Siliciano RF, Blankson JN: Evidence of CD8+ T-cell-mediated selective pressure on human immunodeficiency virus type 1 nef in HLA-B*57+ elite suppressors. J Virol. 2009, 83: 88-97. 10.1128/JVI.01958-08.PubMed CentralView ArticlePubMedGoogle Scholar
- Goulder PJ, Altfeld MA, Rosenberg ES, Nguyen T, Tang Y, Eldridge RL, Addo MM, He S, Mukherjee JS, Phillips MN, et al: Substantial differences in specificity of HIV-specific cytotoxic T cells in acute and chronic HIV infection. J Exp Med. 2001, 193: 181-194. 10.1084/jem.193.2.181.PubMed CentralView ArticlePubMedGoogle Scholar
- Karlsson AC, Iversen AK, Chapman JM, de Oliviera T, Spotts G, McMichael AJ, Davenport MP, Hecht FM, Nixon DF: Sequential broadening of CTL responses in early HIV-1 infection is associated with viral escape. PLoS ONE. 2007, 2: e225-10.1371/journal.pone.0000225.PubMed CentralView ArticlePubMedGoogle Scholar
- Rozera G, Abbate I, Bruselles A, Vlassi C, D'Offizi G, Narciso P, Chillemi G, Prosperi M, Ippolito G, Capobianchi MR: Massively parallel pyrosequencing highlights minority variants in the HIV-1 env quasispecies deriving from lymphomonocyte sub-populations. Retrovirology. 2009, 6: 15-10.1186/1742-4690-6-15.PubMed CentralView ArticlePubMedGoogle Scholar
- Bailey JR, Zhang H, Wegweiser BW, Yang HC, Herrera L, Ahonkhai A, Williams TM, Siliciano RF, Blankson JN: Evolution of HIV-1 in an HLA-B*57-positive patient during virologic escape. J Infect Dis. 2007, 196: 50-55. 10.1086/518515.View ArticlePubMedGoogle Scholar
- Borderia AV, Codoner FM, Sanjuan R: Selection promotes organ compartmentalization in HIV-1: evidence from gag and pol genes. Evolution. 2007, 61: 272-279. 10.1111/j.1558-5646.2007.00025.x.View ArticlePubMedGoogle Scholar
- Li B, Gladden AD, Altfeld M, Kaldor JM, Cooper DA, Kelleher AD, Allen TM: Rapid reversion of sequence polymorphisms dominates early human immunodeficiency virus type 1 evolution. J Virol. 2007, 81: 193-201. 10.1128/JVI.01231-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Avila-Rios S, Reyes-Teran G, Espinosa E: Cornering HIV: taking advantage of interactions between selective pressures. Med Hypotheses. 2007, 69: 422-431. 10.1016/j.mehy.2006.12.012.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.