Early immune adaptation in HIV-1 revealed by population-level approaches
Retrovirology volume 11, Article number: 64 (2014)
The reproducible nature of HIV-1 escape from HLA-restricted CD8+ T-cell responses allows the identification of HLA-associated viral polymorphisms “at the population level” – that is, via analysis of cross-sectional, linked HLA/HIV-1 genotypes by statistical association. However, elucidating their timing of selection traditionally requires detailed longitudinal studies, which are challenging to undertake on a large scale. We investigate whether the extent and relative timecourse of immune-driven HIV adaptation can be inferred via comparative cross-sectional analysis of independent early and chronic infection cohorts.
Similarly-powered datasets of linked HLA/HIV-1 genotypes from individuals with early (median < 3 months) and chronic untreated HIV-1 subtype B infection, matched for size (N > 200/dataset), HLA class I and HIV-1 Gag/Pol/Nef diversity, were established. These datasets were first used to define a list of 162 known HLA-associated polymorphisms detectable at the population level in cohorts of the present size and host/viral genetic composition. Of these 162 known HLA-associated polymorphisms, 15% (occurring at 14 Gag, Pol and Nef codons) were already detectable via statistical association in the early infection dataset at p ≤ 0.01 (q < 0.2) – identifying them as the most consistently rapidly escaping sites in HIV-1. Among these were known rapidly-escaping sites (e.g. B*57-Gag-T242N) and others not previously appreciated to be reproducibly rapidly selected (e.g. A*31:01-associated adaptations at Gag codons 397, 401 and 403). Escape prevalence in early infection correlated strongly with first-year escape rates (Pearson’s R = 0.68, p = 0.0001), supporting cross-sectional parameters as reliable indicators of longitudinally-derived measures. Comparative analysis of early and chronic datasets revealed that, on average, the prevalence of HLA-associated polymorphisms more than doubles between these two infection stages in persons harboring the relevant HLA (p < 0.0001, consistent with frequent and reproducible escape), but remains relatively stable in persons lacking the HLA (p = 0.15, consistent with slow reversion). Published HLA-specific Hazard Ratios for progression to AIDS correlated positively with average escape prevalence in early infection (Pearson’s R = 0.53, p = 0.028), consistent with high early within-host HIV-1 adaptation (via rapid escape and/or frequent polymorphism transmission) as a correlate of progression.
Cross-sectional host/viral genotype datasets represent an underutilized resource to identify reproducible early pathways of HIV-1 adaptation and identify correlates of protective immunity.
HIV-1 escape from Human Leukocyte-Antigen (HLA) class I-restricted CD8+ T-lymphocytes (CTL) occurs in a broadly predictable manner based on the HLA alleles expressed by the host . Reversion of escape mutations, usually to consensus, upon HIV-1 transmission to an individual lacking the restricting HLA also occurs reproducibly in many -, though not all -, cases. The reproducible nature of viral adaptation allows us to identify HLA-associated polymorphisms in HIV-1 (that is, viral polymorphisms that are significantly over- or under- represented among persons expressing a given HLA allele) “at the population level” (that is, via the analysis of cross-sectional, linked HLA/HIV-1 genotypes via statistical association approaches that additionally correct for various potential confounders -). Such studies are normally undertaken in chronic infection, as the virus has undergone a majority of its within-host adaptation by this stage. For example, a recent population-level study of >1800 chronically HIV-1 subtype B-infected persons identified over >2000 HLA-associated polymorphisms across HIV-1, with a majority occurring in Gag, Pol and Nef .
Though HLA-associated polymorphisms in HIV-1 can be identified using cross-sectional approaches, their timing of selection cannot be directly determined by these methods. Rather, temporal information is ideally established via detailed longitudinal study of untreated individuals recently infected with HIV-1 (e.g.: ,,-). However, identifying large numbers of recently-infected persons is challenging. Another consideration is that, given the current evidence and clinical recommendations supporting HIV-1 treatment initiation in early infection , prospective longitudinal observational study of untreated HIV-1 infection may no longer be feasible nor ethical moving forward.
As such, cross-sectional pretreatment host/viral genotype datasets from individuals at different HIV-1 infection stages enrolled in established (or future) cohorts could potentially provide alternate data sources to infer the extent and time course of immune-driven HIV-1 adaptation, including the earliest events post-infection, using population-level approaches. Though such approaches have been investigated ,, they remain underutilized in this context. Notably, population-level approaches offer one key advantage in that - by definition - they specifically identify HIV-1 adaptations that occur reproducibly in persons expressing the restricting HLA  (as opposed to longitudinal studies that characterize immune escape dynamics in individual persons, but cannot elucidate the extent to which such pathways are shared between persons, e.g.,-). As such, population-level studies may be particularly useful in identifying the HLA-restricted CTL escape mutations that are most rapidly and reproducibly selected following HIV-1 infection.
In an attempt to achieve these goals, we undertook a proof-of-concept study that compared the prevalence of known HLA-associated polymorphisms in HIV-1 Gag, Pol and Nef  in identically-sized cross-sectional early and chronic infection cohorts that were matched as closely as possible for their HLA allele distributions and their total HIV-1 diversity. Our main goals were: 1) to assess the utility of population-level approaches to identify the most reproducibly rapid escape mutations in HIV-1; 2) to estimate the extent of escape and reversion between early and chronic infection; and 3) to investigate whether features related to population-level early immune escape signal can discriminate protective from non-protective HLA alleles.
Results and discussion
Assembling early and chronic infection cohorts matched for size, HLA and HIV-1 diversity
Our study sought to demonstrate that the extent, reproducibility and relative timing (early versus later) of HLA-driven escape in HIV-1 can be inferred via comparative analysis of independent cross-sectional host/virus genotype datasets from different infection stages. This strategy ideally requires cross-sectional datasets that are identically powered with respect to host and viral genetic diversity (i.e. datasets that mimic longitudinal data as closely as possible, in that they differ only with respect to infection stage of the participants). As such, our first step was to assemble early and chronic HIV-1 subtype B cohorts of identical size that were matched as closely as possible for HLA class I allele distribution and HIV-1 diversity. We did so by drawing upon host and viral genotype data from early and chronic infection cohorts in North America, Europe and Australia (methods and ,-). Our final early and chronic datasets comprised 221 Gag, 203 Pol and 219 Nef HIV-1 subtype B sequences per cohort, for which linked HLA class I types were available. Early cohort patients were recruited a median of 88 [IQR 63–120] days following infection. All early and >75% of chronic patients were antiretroviral naïve; the remainder were untreated at time of sampling.
A total of 59 HLA class I alleles, classified at subtype-level (4-digit) resolution, were observed at a frequency >1% in the early and/or chronic cohorts; these comprised 17 HLA-A, 23 HLA-B and 19 HLA-C alleles (Figure 1). Of these 59 alleles, the frequencies of 56 (94.9%) were comparable between cohorts; only three alleles (HLA-A*02:06, A*30:02 and B*39:01) exhibited significantly different frequencies between cohorts (all p < 0.05; 0.08 < q < 0.33) (Figure 1A-C). As such, our early and chronic cohorts were generally well-matched with respect to host HLA diversity.
HIV-1 Gag, Pol and Nef diversity was also generally comparable between the two cohorts. Mean patristic (pairwise) genetic distances between HIV-1 sequences in early versus chronic datasets, measured in units of substitutions per nucleotide site, were 0.076 (Standard Deviation [SD] ± 0.011) versus 0.071 (SD ± 0.010) respectively for Gag (Figure 1D left and middle panels), 0.057 (SD ± 0.008) versus 0.053 (SD ± 0.008) for Pol, and 0.119 (SD ± 0.018) versus 0.120 (SD ± 0.021) for Nef (not shown). Moreover, no gross inter-cohort segregation was observed in a combined HIV-1 Gag phylogeny (Figure 1D, right), indicating that neither cohort was dominated by large epidemiologically linked clusters nor exhibited evidence of recent descent from distinct ancestors. Together, these data suggest that our early and chronic datasets are similarly powered with respect to host and viral genetic diversity, and thus differ only with respect to infection stage.
Defining the list of HLA-associated polymorphisms for investigation in cohorts of the present size and composition
A total of 453 HLA-associated polymorphisms in Gag/Pol/Nef had previously been identified at q < 0.05 in an independent cohort of N > 1800 individuals with chronic HIV-1 subtype B infection , which contained no overlap with the cohorts studied here. These HLA-associated polymorphisms comprise “adapted” associations (HIV-1 amino acids significantly over-represented in the presence of the HLA allele in question) as well as “nonadapted” associations (HIV-1 amino acids significantly under-represented in the presence of the HLA allele). For example, at Gag codon 242 the nonadapted amino acid associated with HLA-B*57:01 is the subtype B consensus “T” whereas the B*57:01 adapted form is “N”, denoted as “B*57:01-Gag-T242N”. The cohort wherein these HLA-associations were originally defined however  was more than seven times larger than the cohorts presently studied. Therefore, we do not have sufficient statistical power to interrogate all of them in the present study. As such, our next step was to define, from the published list , the subset of known HLA-associated polymorphisms that is appropriate for study in cohorts of the present size and host/viral genetic composition.
Theoretically, if all immune escape mutations, once selected, persisted for the remainder of the host’s lifetime, and if we had achieved perfect genetic matching between our early and chronic cohorts, then we could define an appropriate subset of HLA-associated polymorphisms by interrogating our chronic cohort for the presence of these N = 453 known HLA-associated polymorphisms. Those detectable at the population level in chronic infection, a stage when a majority of within-host adaptation has already occurred, would represent an appropriate subset for study in cohorts of the present size. Thus we first interrogated our chronic cohort for the presence of these 453 published HLA-associated polymorphisms using statistical association with phylogenetic correction (see Methods and ,), and in doing so identified 157 (35%) “adapted” and “nonadapted” HLA associations at p < 0.01 (corresponding to q < 0.01 in this analysis) (Figure 2 and Additional file 1). These comprised 54, 52 and 51 HLA-associated polymorphisms in Gag, Pol and Nef respectively.
However, if genetic matching of early and chronic cohorts was imperfect, and/or if immune escape mutations occasionally arose transiently during infection, then defining HLA-associated polymorphisms based on the chronic cohort only could potentially lead us to miss some associations. Thus, we interrogated our early cohort for this same list of published HLA-associated polymorphisms using identical methods. In doing so we identified an additional 5 HLA-associated polymorphisms at p < 0.01, two occurring at Gag codon 401 and one each occurring at Nef codons 38, 102 and 133, that had not been identified in the chronic analysis (Figure 2 and Additional file 1).
We therefore took the union of these results, totaling N = 162 HLA-associated polymorphisms, as our definition of known HLA-associated polymorphisms appropriate for cohorts of the present size and composition (Figure 2 and Additional file 1). To further validate this subset, we applied a published phylogenetically-corrected interaction test (see Methods and ,) to compare the strengths of selection of these individual HLA-associated polymorphisms in early versus chronic cohorts. Given that within-host HIV-1 adaptation increases over the infection course via the selection of immune escape mutations, we would expect higher strengths of association between the restricting HLA and the HIV-1 polymorphism at later versus earlier stages: indeed, of the 162 HLA-associated polymorphisms studied, the strengths of 101 (62%) of them differed significantly (p < 0.01, q ≤ 0.1) between early and chronic infection (Additional file 1). Note also that none of the N = 162 polymorphisms were restricted by the three HLA alleles (A*02:06, A*30:02 and B*39:01) whose frequencies differed significantly between cohorts. All downstream analyses therefore focused on this list of N = 162 HLA-associated polymorphisms.
Fifteen percent of known HLA-associated polymorphisms are already detectable at the population level in early infection
Our first objective was to assess the extent of population-level signal for HLA-driven escape in early HIV-1 infection. Of the N = 162 HLA-associated polymorphisms identified for study, 24 (15%), occurring at 14 unique codons in Gag, Pol and Nef, were detectable at the population level in the early cohort at a threshold of p ≤ 0.01 (Table 1). In total, these 24 associations comprised 16% of those investigated in Gag (9 of 56), 6% of those investigated in Pol (3 of 52) and 22% of those investigated in Nef (12 of 54).
As expected, among these were escape mutations known or previously observed to occur in the first year of infection, including B*57:01 Gag-T242N, B*51:01 RT-I135X and Int-L28I, C*03:04 Nef-V85L, A*11:01 Nef-K92R and A*24:02 Nef-Y135F ,,,,,,. Our findings therefore provide proof-of-concept that the most consistently rapid host adaptations in HIV-1 can be identified using cross-sectional methods. It is notable that, by < 3 months post-infection, the magnitude of statistical association between certain HLA alleles and their associated viral polymorphisms is already very high, and in some cases not significantly different from their magnitudes of association in chronic infection. For example, the Odds Ratio [OR] of association between B*57 and Gag-T242N is 33 in early infection (p = 3 × 10-9) compared to 151 (p = 5 × 10-16) in chronic infection, which, though stronger during the latter stage, does not represent a statistically significant difference (inter-cohort comparison p = 0.3) (Table 1). This observation underscores the rapid and highly reproducible nature of certain HLA-driven within-host adaptations in HIV-1, where, for certain mutations such as Gag-T242N, escape (and by extension our ability to detect this association via population-level methods) is already near maximal in B*57:01-expressing persons < 3 months post-infection.
Of note, HLA-associated HIV-1 polymorphisms with strong early population-level escape signal also included understudied viral sites. Notable among these were A*31:01-associated polymorphisms at Gag codons 397, 401 and 403, the first of which represented the second strongest p-value detected in the early cohort (Odds Ratio = 6.2, p = 1 × 10-7, Table 1). These associations are located within the novel A*31:01-restricted CR9 CD8+ epitope originally characterized via detailed longitudinal analysis of a single HIV-1 subtype B-infected individual ,. By definition, population-level studies identify viral adaptations that occur reproducibly in persons expressing the restricting HLA; as such, the present results extend those of the original individual-level study  by indicating that escape within CR9 is both rapid and highly consistent in HLA-A*31:01-expressing persons. By extension, a lack of population-level early escape signal does not necessarily mean that a given site never escapes early: rather, it indicates that a given site does not reproducibly escape early (or at least does not do so to an extent that achieves statistical significance in a dataset of the present size). For example, very rapid (<30 days) escape was previously documented within the A*01-restricted GY9 epitope (Gag codons 71–79) in two HIV-1 subtype C-infected persons , but no evidence of reproducible early escape Gag codon 79 in A*01-expressing persons was observed in our early dataset, suggesting that rapid A*01-driven escape at this position is atypical in HIV-1 subtype B.
Taken together, population-level analyses extend those of individual-level studies by identifying escape mutations that are rapidly and reproducibly selected across patients. Our observation that 15% of known HLA-associated polymorphisms, notably those in Nef and Gag, are already detectable < 3 months post-infection, underscores the predictable and rapid nature of HIV-1 adaptation despite each individual’s unique combination of host HLA and transmitted virus genetics. Further, the detection of substantial population-level escape signal within unknown or understudied CD8+ epitopes in HIV-1 (Table 1) argues for continued efforts to map novel epitopes commonly targeted during this critical infection stage.
Can population-level approaches identify transient early escape pathways?
Recent longitudinal studies have revealed that immune escape is often characterized by the initial appearance of transient mutant forms that often retain some ability to be targeted by existing (or de novo) CTL ,, which then drive the selection of more effective escape variants that ultimately become fixed within the host ,,. If such “transient” escape pathways are reproducible across hosts, we wondered whether population-level approaches could theoretically be used as exploratory tools to identify them. If so, we reasoned that such transient escapes would display stronger population-level escape signal in early compared to chronic infection (since, in some persons, the early variant would be subsequently replaced with another, thereby reducing population-level signal in later stages). Although our inter-cohort comparative analysis revealed no HLA-associated polymorphisms that displayed significantly stronger signal in early versus chronic infection (Additional file 1 and data not shown), we were nevertheless intrigued by the five HLA-associated polymorphisms in Gag and Nef that exhibited population-level escape signal of p < 0.01 in our early cohort but p ≥ 0.01 in the chronic cohort (Table 1), suggesting these as possible transient escape pathways.
Indeed, analysis of available longitudinal bulk plasma HIV-1 RNA Nef sequences from seven B*57:01 expressing individuals identified one case where an individual harbored the Nef-133I adapted mutation at the earliest sampled timepoint 30 days post-infection, which was replaced by V at 86 days post-infection and then by a mixture of I/V at 228 days-post infection (Table 2). Similarly, analysis of available bulk plasma Gag sequences from seven A*30:01 expressing persons identified one case where an individual harbored the adapted Gag-401 L variant at the earliest timepoint post-infection, that was subsequently replaced by a non-adapted form (and/or a mixture of the two) within a year of infection. In both of these cases the association is located at position 7 within the epitope, which is consistent with transient early escape mutations representing incomplete TCR repertoire escape variants ,. The idea that escape mutations, once selected, may not always persist for the lifetime of the host is also supported by within-host reversion of certain escape mutations in very advanced disease . We thus cautiously interpret the data to suggest that cross-sectional approaches could theoretically be used to identify reproducible HLA-driven adaptations that represent “transient” early escape variants in some individuals, though such findings would require validation in independent cohorts, as well as experimentally.
Escape prevalence in early infection correlates with longitudinal first-year escape rates
Another objective was to investigate to what extent cross-sectional data could be used to infer the extent and time course of immune-driven HIV-1 adaptation. As such, we first wished to demonstrate that early escape frequencies calculated cross-sectionally predict rates of immune escape calculated longitudinally. Published first-year rates of escape were available for 27 optimally-described CD8+ T-cell epitopes  which contained one or more HLA-associated polymorphisms investigated in the present study. For example, the estimated first year escape rate for the Gag-TW10 epitope (Gag240-249) is 38.36% per person-month , while the prevalence of the Gag-T242N mutation among B*57-expressing persons in our early dataset is 67% (Figure 3). As expected, longitudinal first-year CD8+ epitope escape rates correlated significantly with HIV-1 polymorphism prevalence among persons expressing the relevant HLA in our early infection dataset (Pearson’s R = 0.68, p = 0.0001; Figure 3). Because ~40% of the patients in the present early infection cohort were included in the published longitudinal study , we re-analyzed our data with these overlapping patients removed, and observed that the correlation remained strong (Pearson’s R = 0.55, p = 0.0035, not shown). This supports HLA-associated escape mutation prevalence calculated cross-sectionally at < 3 months post-infection as a reliable surrogate marker of first year escape rates calculated longitudinally.
Inferring the extent of host adaptation via comparative analysis of cross-sectional data from early and chronic infection
We next wished to use our cross-sectional early and chronic cohorts to quantify the extent of HLA-driven escape occurring between these two infection stages. For this analysis, we specifically defined “escape” as the specific adapted viral form associated with a given HLA allele at a given HIV-1 codon - for example, Gag 242 N is the B*57:01-associated adapted form at this position. This adapted list comprised N = 74 HLA-associated polymorphisms (25, 24, and 25 in Gag, Pol and Nef respectively) (Figure 2 and Additional file 1). We calculated the prevalence of each of these polymorphisms in persons expressing the relevant HLA allele in our early versus chronic cohorts, thus allowing us to estimate the extent of within-host HIV-1 adaptation between these two stages. Overall, the median “percentage escaped” (defined as the % of individuals expressing the relevant HLA and harboring the HIV-1 polymorphism of interest) was 23.8% [Interquartile range (IQR) 5.3-44.4%] in early infection versus 55.1% [IQR 28.4-73.0%] in chronic infection (p < 0.0001; Figure 4A). This indicates that, on average, escape prevalence in persons expressing the restricting HLA allele more than doubles between these infection stages. Breaking the analysis down by HIV-1 protein, the median early versus chronic escape prevalence was 23.5% [IQR 15.1-49.7%] vs. 55.6% [IQR 29.7-85.7%] in Gag, 11.3% [IQR 1.2-33.3%] vs. 50.5% [IQR 25.2-69.2%] in Pol, and 31.3% [7.5-63.0%] vs. 54.6% [21.4-73.3%] in Nef (all p ≤ 0.001, not shown). This is consistent with early escape occurring predominantly in Gag and Nef ,,, while escape in Pol is generally slower but nevertheless approaches comparable levels by chronic infection.
Though summary statistics are informative, individual polymorphisms differ widely in their timing and extent of selection over the infection course. For this reason, details on polymorphism prevalence in persons expressing vs. not expressing the relevant HLA in early and chronic infection, along with their statistical measures of association, are provided in the Additional file 1. We highlight some examples here. First, for a substantial minority of polymorphisms (notably those in Table 1), escape is rapid, reproducible and largely complete within < 3 months post-infection. For example, 67% of B*57:01-expressing individuals already harbored Gag-242 N in early infection, a proportion that increased to 83% in the chronic phase (Figure 4A). Noting that Gag-242 N frequency was only 5.8% among persons lacking B*57:01 in early infection (Additional file 1), these results are consistent with escape having already occurred in over two-thirds of B*57-expressing persons by < 3 months post-infection ,, with an additional minority escaping somewhat later.
While the prevalence of Gag-242 N is low in the general population (5.8% among B*57:01-negative persons and ~1% among persons lacking an allele belonging to the B58 supertype), other polymorphisms are quite prevalent in circulation, but are nevertheless significantly enriched among HLA-expressing persons in early infection (Figure 4A). In this case, their high early prevalence is attributable to both frequent transmission and reproducible early escape. For example, both the A*31:01-associated Gag-403 K and C*03:04-associated Nef-85 L polymorphisms are observed at >40% prevalence in HIV-1 subtype B sequences, but their prevalence is ~55% and ~72% respectively among persons expressing the relevant HLA < 3 months post-infection (Odds Ratios 1.5 and 2.4 respectively, p < 0.01, Table 1 and Additional file 1). The observation that population-level approaches are capable of detecting strong escape signals despite high polymorphism background frequencies has previously been demonstrated in high-powered chronic infection cohorts ; the present study extends this to demonstrate such signals can also be detected very early in infection, in more modestly-powered datasets. Overall, if one uses the original criterion of early population-level statistical signal of p ≤ 0.01 to define HIV-1 sites that predominantly escape early, 15% (11 of 74) of adapted polymorphisms fall into this category; using a more liberal threshold of p < 0.05, this increases to 21.6% (16 of 74) (Figure 4B).
The remaining 78.4% (58 of 74) polymorphisms generally reproducibly escape later than 3 months following infection (Figure 4B). It is important to note that later escape can occur because CTL responses against these regions generally arise later during infection (i.e. there is no immune pressure on these epitopes in early infection), or because CTL responses arise relatively early, but escape does not reproducibly occur rapidly in a significant proportion of individuals expressing the relevant HLA allele. Among these later-escaping polymorphisms are those whose population background (transmission) frequencies are generally low, and those whose background frequencies are generally high. The B*51:01-associated Integrase-32I polymorphism at position 5 of the B*51-restricted LI9 epitope (Integrase28-36) provides an example of the former. In early infection, its frequency in B*51:01-expressing persons is 5%, not significantly different from background, but this rises to 64% by chronic infection (Figure 4A and Additional file 1). The LI9 epitope is known to be consistently targeted in B*51-expressing persons early after infection ,. The observation that this epitope ultimately escapes via Int-32I in >60% of B*51:01-expressing persons suggests this epitope is under strong, sustained and reproducible CD8+ T-cell pressure by B*51 in vivo, where delayed escape is likely explained by a combination of mutational/fitness constraints and both intra-individual (“vertical”) and inter-individual (“horizontal”) CD8+ T-cell immunodominance hierarchies ,,.
An example of a later-escaping polymorphism with high population background frequency is B*44:02-Gag-312E. Located at position 7 of the B*44:02-restricted AW11 epitope (Gag306-317), it represents the HIV-1 subtype B consensus residue at this codon. Its >60% frequency in both B*44:02 and non-B*44:02-expressing persons in early infection reflects its high transmission frequency, rather than early selection by B*44:02. Nevertheless, by chronic infection, 83% of B*44:02-expressing persons harbored Gag-312E, consistent with later escape (Figure 4A). A full categorization of HLA-associated “adapted” polymorphisms in terms of “earlier” vs. “later” escaping (defined as early p < 0.05 vs. p ≥ 0.05 respectively) and “lower” vs. “higher” background (estimated transmission) frequency (defined as <30% vs. p ≥ 30% respectively), is provided in the Additional file 1. A graphic depicting the proportion of HLA-associated polymorphisms in each of these categories, broken down by HIV protein, is provided in Figure 4B.
The extent of reversion of HLA-associated polymorphisms over time can be similarly estimated by calculating their prevalence in persons lacking the relevant HLA in the early versus chronic cohorts. The overall median percentage of individuals harboring a given polymorphism in the absence of the restricting HLA allele was comparable in early (13.7% [IQR 4.9-34.2%]) and chronic (14.9% [4.2-34.9%]) infection (p = 0.15; Figure 4C), consistent with slow reversion reported for many polymorphisms ,,,. Note that inferred reversion frequencies merit cautious interpretation in cases where polymorphisms are selected by multiple alleles (e.g. the seemingly stable prevalence of Gag-147 L in A*25:01-negative individuals is likely due in part to its selection by B*13:02 and B*57:01 , associations that were not investigated in the present study). Nevertheless, results confirmed that HLA-B*57:01-Gag-T242N reverts between early and chronic infection (though RT-245E, Int-122I, Int-124N, or Nef-116N revert slowly or not at all, as reported previously ,,). The reversion analysis additionally revealed novel sites of potential interest. For example, the early escaping A*31:01-Gag 397R polymorphism (Table 1 and ,) displayed evidence of reversion, suggesting that this mutation may have a high fitness cost.
Host adaptation-related features distinguish protective and non-protective HLA class I alleles
Lastly, we wished to identify adaptation-related features that discriminate protective from non-protective HLA alleles, defined here as their published hazard ratios for progression to AIDS [HR-AIDS] in natural history studies . Although the timecourse of viral escape is influenced by complex factors including epitope immunodominance hierarchies, strength of selection, mutational/fitness constraints and transmitted virus characteristics ,,,-, we reasoned that HLA alleles that restrict polymorphisms that are already highly prevalent in early infection (due to rapid escape and/or frequent transmission) would be generally unfavorable for HIV-1 control. Thus, for all HLA alleles for which ≥2 adapted polymorphisms were investigated in the present study (N = 17 alleles total), we computed their mean “percentage escaped” in early infection. That is, we took the prevalence of each of these adapted polymorphisms in persons expressing the relevant HLA allele in our early cohort (displayed in Figure 4A), and, for each HLA allele, computed the mean of these values.
Consistent with our hypothesis, we observed a positive correlation between an HLA allele’s average extent of adaptation in early infection, and its HR-AIDS (Pearson’s R = 0.53, p = 0.028; Figure 5A). Of note, HLA-B*57:01 appears as somewhat of an outlier, exhibiting higher than expected escape prevalence in early infection given its protective nature. We hypothesize that the reason B*57 can maintain sustained HIV-1 control despite rapid escape in some epitopes (Table 1 and ,,,,) is because the early B*57-restricted CD8+ response often simultaneously targets more than one epitope, notably in p24Gag,-, where escape is accompanied by fitness costs -.
HLA-associated polymorphisms identified at the population level mark viral sites under strong, reproducible in vivo immune pressure by individual HLA alleles . We thus further hypothesized that HLA alleles for which escape was substantial (i.e. selected in a high proportion of persons expressing the relevant HLA) but generally delayed (i.e. selected post early-phase) would tend to be more protective. As such, for the same set of HLA-associated polymorphisms we computed their fold-change in escape between early and chronic infection, and calculated the mean of these values per HLA allele. By this measure, alleles for which the majority of escape had already occurred in early infection would exhibit low subsequent fold-changes, whereas alleles selecting escape mutations that generally occurred later in infection would exhibit fold-changes reflecting the extent of selection pressure on these sites in persons expressing the relevant HLA. Consistent with our hypothesis, we observed an inverse correlation between an HLA allele’s HR-AIDS and its average fold-increase in escape in chronic versus early infection (Pearson’s R = 0.54, p = 0.025; Figure 5B). Of note, B*27 appears as an additional outlier in this analysis, possibly due to escape in the critical Gag-KK10 epitope requiring nearly a decade in some individuals ,, due to its high fitness/mutational barrier ,.
Overall, our findings are consistent with a high early burden of adaptation to host HLA (either via rapid escape or frequent polymorphism transmission) as a correlate of HLA-associated progression risk. Conversely, HLA alleles from which HIV-1 escape is substantial and reproducible yet occurs on a delayed timescale appears to be a correlate of protection. Taken together with observations that protective HLA alleles contribute substantially to the total HIV-specific CD8+ response in early infection ,, that they impose broad selection pressures on HIV-1 ; and that some CD8+ epitopes escape slowly despite sustained CD8+ targeting ,, results suggest that the capacity to exert consistent, substantial and sustained pressure, ideally on multiple epitopes, from which the virus can only escape on a relatively delayed timescale, is a correlate of protection.
Some limitations of our study merit mention. Firstly, our early and chronic datasets are relatively modestly powered by association testing standards, so it was not possible to examine all published HLA-associated polymorphisms in HIV-1 subtype B. Secondly, due to the lack of information on duration of infection for chronic patients, it is likely that the chronic cohort comprised patients at a range of infection stages. Inclusion of some chronic patients with less advanced infection could underestimate the extent of escape at this stage. Finally, although care was taken to match our early and chronic datasets as closely as possible for HIV and HLA genetic diversity and distribution, it is essentially impossible to achieve perfectly matched datasets. As such, we cannot rule out small differences in statistical power between cohorts, and therefore should interpret candidate “transient” early escape results with some caution.
In conclusion, our results provide proof-of-concept that statistical association approaches can be applied to cross-sectional host/viral genetic datasets to identify the most rapidly selected HLA-associated polymorphisms in HIV-1 that are also reproducibly selected across patients. As such, results from these types of population-level studies complement those of individual-level longitudinal analyses that cannot assess inter-patient reproducibility. Furthermore, the extent and relative timing of individual escape events (in terms of early versus later in the infection course) can also be inferred from cross-sectional data. In particular, we demonstrate that high escape prevalence in early infection (either due to rapid selection and/or high transmission frequency) is a correlate of HLA-associated progression risk while reproducible later escape (a surrogate of consistent immune selection on a given site in persons expressing the relevant HLA) is a correlate of protection.
Given that longitudinal observational studies of untreated persons are incompatible with current recommendations for early HIV-1 treatment initiation  and treatment as prevention , cross-sectional analyses of pretreatment host/viral genotypes could provide relevant alternative tools to advance our knowledge of HIV-1 adaptation, including the earliest events post-infection. We suggest that studies such as the present one be undertaken with expanded cross-sectional cohorts, comprised of individuals at various clinical stages of infection, including from different HIV-1 subtypes.
Early and chronic infection cohorts
The early cohort was comprised of HIV-1 subtype B infected patients recruited through various observational seroconverter studies including the Acute Infection and Early Disease Research Program (AIEDRP) sites in Boston and New York (USA), Sydney (Australia), a private medical clinic in Berlin (Germany), and observational cohort studies in Montreal and Vancouver (Canada) ,,. Infection dates for the patients in the early cohort were estimated as described in ,. Briefly, for patients with positive HIV RNA (>5,000 copies/ml) or detectable serum p24 antigen but a negative HIV-1 enzyme immunoassay (EIA), 4 weeks were subtracted from the negative EIA date. For patients with positive EIA but negative/indeterminate Western blot, 6 weeks were subtracted from the positive EIA date. For patients with negative detuned EIA, 4 months were subtracted from this date. For the remainder, infection dates were estimated as the midpoint between the last negative and the first positive HIV test. Clinical histories were incorporated into infection date estimates where available.
To maximize our power to detect HLA-associated polymorphisms in the early infection stage, all available early infection patients were included in the present study. For each of these patients, the sample closest to ~3-months following the estimated date of infection was selected, yielding a median sampling distribution of 88 days [Interquartile Range 63–120 days] post-infection for early samples. In contrast, the chronic cohort was assembled from the baseline (pre-therapy) timepoint from a total of more than 300 HIV-1 subtype B infected individuals initiating antiretroviral therapy in British Columbia, Canada, and untreated HIV-1 subtype B infected individuals in Boston, USA ,. Time since infection is unknown for individuals in the chronic cohort, however the median CD4 count at sampling was 250 [IQR 147–360] cells/mm3 for this group. To create HIV-1 gene-specific chronic infection datasets of equal size to the early cohort, that were also matched as closely as possible for HLA class I and HIV-1 diversity of the early cohort, chronic patients were selected from the total group using an iterative process to achieve the closest matching of HIV-1 and HLA distributions (Figure 1). All early and >75% of chronic patients were antiretroviral naïve; the remainder were untreated at time of sampling.
All patients provided written informed consent. Ethical approval was obtained through the institutional review boards at the Massachusetts General Hospital, the BC Centre for Excellence in HIV/AIDS and Simon Fraser University.
HIV-1 and host (HLA class I) genotyping
HIV-1 RNA was extracted from plasma using standard methods. Gag, Pol (including protease, codons 1–400 of Reverse Transcriptase, and Integrase), and Nef were amplified in separate nested RT-PCR reactions using HIV-1 subtype B-specific primers. Amplicons were bulk-sequenced bidirectionally on a 3130xl and/or 3730xl automated DNA sequencer (Applied Biosystems). Chromatograms were analyzed using Sequencher v5.0 (Genecodes) or RECall  with nucleotide mixtures called if the height of the secondary peak exceeded 25% of the height of the dominant peak (Sequencher) or 20% of the dominant peak area (RECall). HIV-1 sequences were confirmed as subtype B using the recombinant identification program (RIP; http://www.hiv.lanl.gov/content/sequence/RIP/RIP.html) and aligned to the HIV-1 subtype B reference strain HXB2. Phylogenetic trees were constructed using PhyML  and visualized using FigTree (http://tree.bio.ed.ac.uk/software/figtree/). Pairwise genetic distances were computed from newick treefiles using PATRISTIC . HLA class I typing was performed using sequence-based methods  and imputed where necessary to high resolution using a machine learning algorithm (; http://research.microsoft.com/en-us/projects/bio/mbt.aspx#HLA-Completion). HIV-1 sequences from cohorts where REBs allow public sequence deposition have been deposited in GenBank: accession numbers are Gag (KJ869442 - KJ869609), Protease-RT (KJ869900 - KJ870015), Integrase (KJ869610 - KJ869735), Nef (KJ869736 - KJ869899). A full summary of polymorphism frequencies, broken down by HLA allele carriage and infection stage is provided as Additional file 1.
Definition and identification of HLA-associated polymorphisms
The published reference list of N = 453 HLA-associated polymorphisms in HIV-1 subtype B Gag, Pol and Nef sequences was defined in an independent international cohort of >1800 individuals chronically infected with HIV-1 subtype-B using phylogenetically-informed methods at q < 0.05 . The cohort used to define these associations had no overlap with the early and chronic cohorts studied here . Briefly, to identify HLA-associated polymorphisms in linked HIV/HLA datasets, maximum likelihood phylogenetic trees (one per HIV-1 gene) are constructed, and a model of conditional adaptation is inferred for each observed HIV-1 amino acid at each codon. The amino acid is assumed to evolve independently along the tree until it reaches the tips, representing the present host. Selection via host HLA-mediated pressures and HIV-1 amino acid covariation is directly modeled using a weighted logistic regression, in which the individual’s HLA repertoire and covarying HIV-1 amino acids are used as predictors, and the bias is determined by the inferred possible transmitted sequences (as inferred via reconstruction of the amino acid frequencies at the penultimate internal nodes in the phylogeny) . Here, the null hypothesis is that the observed amino acids at the tree tips are explained by the phylogeny only, whereas the alternative hypothesis is that they are better explained by the presence of a specific HLA (or covarying HIV-1 amino acid) in the present host. To identify which factors (HLA and/or HIV-1 covariation) contribute to the selection pressure, a forward selection procedure is employed where the most significant association is added to the model in an iterative fashion, with p-values computed using the likelihood ratio test. Statistical significance is reported using q-values , the p-value analogue of the false discovery rate (FDR). Q-values denote the expected proportion of false positives among results deemed significant at a given p-value threshold; for example, at q ≤ 0.05, we expect 5% of identified associations to be false positives.
HLA-associated polymorphisms are classified into two categories: (1) “Adapted forms”, amino acids significantly overrepresented in the presence of the HLA allele in question, which represent the putative escape forms associated with that HLA at that codon, and (2) “Nonadapted forms”, amino acids significantly underrepresented in the presence of the HLA allele in question, which represent the immunologically susceptible form associated with that HLA at that codon. In most cases, HLA-associated nonadapted forms represent the subtype consensus amino acid while adapted forms represent polymorphic variants – but exceptions exist.
To identify an appropriate subset of known HLA-associated HIV-1 polymorphisms that are appropriate for study in datasets of the present size (N ~ 200) and host/viral genetic distribution, we interrogated our early and chronic infection cohorts for these N = 453 published polymorphisms  using the phylogenetically-corrected methods described above. As described in the results, this yielded a subset of N = 162 HLA-associated polymorphisms detectable in our early and/or chronic datasets (Figure 2 and Additional file 1).
Our analyses also featured comparisons of the strength of selection of HLA-associated polymorphisms between early and chronic cohorts, undertaken using a previously-described phylogenetically-corrected interaction test ,,. Briefly, we took the union of all HLA-associated polymorphisms detectable at the population level in either the early or chronic cohorts (N = 162). For each association on the list, we constructed a phylogenetically-corrected logistic regression model using the restricting HLA as a predictor. Using a likelihood ratio test, we compare this model to a more expressive one that includes an additional interaction term that assigns “1” if the individual expresses the restricting HLA allele and is in the chronic cohort, or “0” otherwise. This allows us to obtain a p-value testing the null hypothesis that HLA-associated selection at that site is not significantly different in early versus chronic cohorts.
Fisher’s exact test was used to compare HLA class I allele frequencies between cohorts. The Mann–Whitney paired test was used to compare the prevalence of HLA-associated polymorphisms in the presence/absence of their restricting HLA, in early versus chronic cohorts, as these data were non-normally distributed. Pearson’s correlation was used to investigate the relationship between early escape prevalence and published first-year rates of escape , as well as with published HLA allele-specific Hazard Ratios for progression to AIDS , as these data did not significantly violate the assumption that values were drawn from a normal distribution. In single analyses, significance is denoted by p < 0.05. In the case of multiple tests, q-values are used ; thresholds are defined throughout the paper. All tests of significance were two-tailed.
Goulder PJ, Watkins DI: HIV and SIV CTL escape: implications for vaccine design. Nat Rev Immunol. 2004, 4: 630-640. 10.1038/nri1417.
Leslie AJ, Pfafferott KJ, Chetty P, Draenert R, Addo MM, Feeney M, Tang Y, Holmes EC, Allen T, Prado JG, Altfeld M, Brander C, Dixon C, Ramduth D, Jeena P, Thomas SA, St John A, Roach TA, Kupfer B, Luzzi G, Edwards A, Taylor G, Lyall H, Tudor-Williams G, Novelli V, Martinez-Picado J, Kiepiela P, Walker BD, Goulder PJ: HIV evolution: CTL escape mutation and reversion after transmission. Nat Med. 2004, 10: 282-289. 10.1038/nm992.
Henn MR, Boutwell CL, Charlebois P, Lennon NJ, Power KA, Macalalad AR, Berlin AM, Malboeuf CM, Ryan EM, Gnerre S, Zody MC, Erlich RL, Green LM, Berical A, Wang Y, Casali M, Streeck H, Bloom AK, Dudek T, Tully D, Newman R, Axten KL, Gladden AD, Battis L, Kemper M, Zeng Q, Shea TP, Gujja S, Zedlack C, Gasser O, et al: Whole genome deep sequencing of HIV-1 reveals the impact of early minor variants upon immune recognition during acute infection. PLoS Pathog. 2012, 8: e1002529-10.1371/journal.ppat.1002529.
Crawford H, Prado JG, Leslie A, Hue S, Honeyborne I, Reddy S, van der Stok M, Mncube Z, Brander C, Rousseau C, Mullins JI, Kaslow R, Goepfert P, Allen S, Hunter E, Mulenga J, Kiepiela P, Walker BD, Goulder PJ: Compensatory mutation partially restores fitness and delays reversion of escape mutation within the immunodominant HLA-B*5703-restricted Gag epitope in chronic human immunodeficiency virus type 1 infection. J Virol. 2007, 81: 8346-8351. 10.1128/JVI.00465-07.
Fryer HR, Frater J, Duda A, Roberts MG, Phillips RE, McLean AR: Modelling the evolution and spread of HIV immune escape mutants. PLoS Pathog. 2010, 6: e1001196-10.1371/journal.ppat.1001196.
Leslie A, Kavanagh D, Honeyborne I, Pfafferott K, Edwards C, Pillay T, Hilton L, Thobakgale C, Ramduth D, Draenert R, Le Gall S, Luzzi G, Edwards A, Brander C, Sewell AK, Moore S, Mullins J, Moore C, Mallal S, Bhardwaj N, Yusim K, Phillips R, Klenerman P, Korber B, Kiepiela P, Walker B, Goulder P: 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.
Schneidewind A, Brumme ZL, Brumme CJ, Power KA, Reyor LL, O'Sullivan K, Gladden A, Hempel U, Kuntzen T, Wang YE, Oniangue-Ndza C, Jessen H, Markowitz M, Rosenberg ES, Sekaly RP, Kelleher AD, Walker BD, Allen TM: Transmission and long-term stability of compensated CD8 escape mutations. J Virol. 2009, 83: 3993-3997. 10.1128/JVI.01108-08.
Cornelissen M, Hoogland FM, Back NK, Jurriaans S, Zorgdrager F, Bakker M, Brinkman K, Prins M, van der Kuyl AC: Multiple transmissions of a stable human leucocyte antigen-B27 cytotoxic T-cell-escape strain of HIV-1 in The Netherlands. AIDS. 2009, 23: 1495-1500. 10.1097/QAD.0b013e32832d9267.
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.
Brumme ZL, John M, Carlson JM, Brumme CJ, Chan D, Brockman MA, Swenson LC, Tao I, Szeto S, Rosato P, Sela J, Kadie CM, Frahm N, Brander C, Haas DW, Riddler SA, Haubrich R, Walker BD, Harrigan PR, Heckerman D, Mallal S: HLA-associated immune escape pathways in HIV-1 subtype B Gag, Pol and Nef proteins. PLoS ONE. 2009, 4: e6687-10.1371/journal.pone.0006687.
Carlson JM, Brumme CJ, Martin E, Listgarten J, Brockman MA, Le AQ, Chui CK, Cotton LA, Knapp DJ, Riddler SA, Haubrich R, Nelson G, Pfeifer N, Deziel CE, Heckerman D, Apps R, Carrington M, Mallal S, Harrigan PR, John M, Brumme ZL: Correlates of protective cellular immunity revealed by analysis of population-level immune escape pathways in HIV-1. J Virol. 2012, 86: 13202-13216. 10.1128/JVI.01998-12.
Carlson JM, Listgarten J, Pfeifer N, Tan V, Kadie C, Walker BD, Ndung'u T, Shapiro R, Frater J, Brumme ZL, Goulder PJ, Heckerman D: Widespread Impact of HLA Restriction on Immune Control and Escape Pathways of HIV-1. J Virol. 2012, 86: 5230-5243. 10.1128/JVI.06728-11.
Brumme ZL, Brumme CJ, Carlson J, Streeck H, John M, Eichbaum Q, Block BL, Baker B, Kadie C, Markowitz M, Jessen H, Kelleher AD, Rosenberg E, Kaldor J, Yuki Y, Carrington M, Allen TM, Mallal S, Altfeld M, Heckerman D, Walker BD: Marked epitope- and allele-specific differences in rates of mutation in human immunodeficiency type 1 (HIV-1) Gag, Pol, and Nef cytotoxic T-lymphocyte epitopes in acute/early HIV-1 infection. J Virol. 2008, 82: 9216-9227. 10.1128/JVI.01041-08.
Duda A, Lee-Turner L, Fox J, Robinson N, Dustan S, Kaye S, Fryer H, Carrington M, McClure M, McLean AR, Fidler S, Weber J, Phillips RE, Frater AJ: HLA-associated clinical progression correlates with epitope reversion rates in early human immunodeficiency virus infection. J Virol. 2009, 83: 1228-1239. 10.1128/JVI.01545-08.
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.
Yager N, Robinson N, Brown H, Flanagan P, Frater J, Fidler S, Weber J, Phillips R: Longitudinal analysis of an HLA-B*51-restricted epitope in integrase reveals immune escape in early HIV-1 infection. AIDS. 2013, 27: 313-323. 10.1097/QAD.0b013e32835b8cf5.
Fischer W, Ganusov VV, Giorgi EE, Hraber PT, Keele BF, Leitner T, Han CS, Gleasner CD, Green L, Lo CC, Nag A, Wallstrom TC, Wang S, McMichael AJ, Haynes BF, Hahn BH, Perelson AS, Borrow P, Shaw GM, Bhattacharya T, Korber BT: Transmission of single HIV-1 genomes and dynamics of early immune escape revealed by ultra-deep sequencing. PLoS ONE. 2010, 5: e12303-10.1371/journal.pone.0012303.
Goonetilleke N, Liu MK, Salazar-Gonzalez JF, Ferrari G, Giorgi E, Ganusov VV, Keele BF, Learn GH, Turnbull EL, Salazar MG, Weinhold KJ, Moore S, Letvin N, Haynes BF, Cohen MS, Hraber P, Bhattacharya T, Borrow P, Perelson AS, Hahn BH, Shaw GM, Korber BT, McMichael AJ: The first T cell response to transmitted/founder virus contributes to the control of acute viremia in HIV-1 infection. J Exp Med. 2009, 206: 1253-1272. 10.1084/jem.20090365.
Liu MK, Hawkins N, Ritchie AJ, Ganusov VV, Whale V, Brackenridge S, Li H, Pavlicek JW, Cai F, Rose-Abrahams M, Treurnicht F, Hraber P, Riou C, Gray C, Ferrari G, Tanner R, Ping LH, Anderson JA, Swanstrom R, Cohen M, Karim SS, Haynes B, Borrow P, Perelson AS, Shaw GM, Hahn BH, Williamson C, Korber BT, Gao F, Self S, et al: Vertical T cell immunodominance and epitope entropy determine HIV-1 escape. J Clin Invest. 2013, 123: 380-393.
Panel on Antiretroviral Guidelines for Adults and Adolescents. Guidelines for the use of antiretroviral agents in HIV-1-infected adults and adolescents. Bethesda (MD): Department of Health and Human Services (DHHS); 2013.
Fryer HR, Frater J, Duda A, Palmer D, Phillips RE, McLean AR: Cytotoxic T-lymphocyte escape mutations identified by HLA association favor those which escape and revert rapidly. J Virol. 2012, 86: 8568-8580. 10.1128/JVI.07020-11.
Huang KH, Goedhals D, Carlson JM, Brockman MA, Mishra S, Brumme ZL, Hickling S, Tang CS, Miura T, Seebregts C, Heckerman D, Ndung'u T, Walker B, Klenerman P, Steyn D, Goulder P, Phillips R, van Vuuren C, Frater J: Progression to AIDS in South Africa Is Associated with both Reverting and Compensatory Viral Mutations. PLoS One. 2011, 6: e19018-10.1371/journal.pone.0019018.
Brockman MA, Chopera DR, Olvera A, Brumme CJ, Sela J, Markle TJ, Martin E, Carlson JM, Le AQ, McGovern R, Cheung PK, Kelleher AD, Jessen H, Markowitz M, Rosenberg E, Frahm N, Sanchez J, Mallal S, John M, Harrigan PR, Heckerman D, Brander C, Walker BD, Brumme ZL: Uncommon pathways of immune escape attenuate HIV-1 integrase replication capacity. J Virol. 2012, 86: 6913-6923. 10.1128/JVI.07133-11.
Poon AF, McGovern RA, Mo T, Knapp DJ, Brenner B, Routy JP, Wainberg MA, Harrigan PR: Dates of HIV infection can be estimated for seroprevalent patients by coalescent analysis of serial next-generation sequencing data. AIDS. 2011, 25: 2019-2026. 10.1097/QAD.0b013e32834b643c.
Cotton LA, Kuang XT, Le AQ, Carlson JM, Chan B, Chopera DR, Brumme CJ, Markle TJ, Martin E, Shahid A, Anmole G, Mwimanzi P, Nassab P, Penney KA, Rahman MA, Milloy MJ, Schechter MT, Markowitz M, Carrington M, Walker BD, Wagner T, Buchbinder S, Fuchs J, Koblin B, Mayer KH, Harrigan PR, Brockman MA, Poon AF, Brumme ZL: Genotypic and Functional Impact of HIV-1 Adaptation to Its Host Population during the North American Epidemic. PLoS Genet. 2014, 10: e1004295-10.1371/journal.pgen.1004295.
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.
Chikata T, Carlson JM, Tamura Y, Borghan MA, Naruto T, Hashimoto M, Murakoshi H, Le AQ, Mallal S, John M, Gatanaga H, Oka S, Brumme ZL, Takiguchi M: Host-specific adaptation of HIV-1 subtype B in the Japanese population. J Virol. 2014, 88: 4764-4775. 10.1128/JVI.00147-14.
Turnbull EL, Baalwa J, Conrod KE, Wang S, Wei X, Wong M, Turner J, Pellegrino P, Williams I, Shaw GM, Borrow P: Escape is a more common mechanism than avidity reduction for evasion of CD8+ T cell responses in primary human immunodeficiency virus type 1 infection. Retrovirology. 2011, 8: 41-10.1186/1742-4690-8-41.
Brackenridge S, Evans EJ, Toebes M, Goonetilleke N, Liu MK, di Gleria K, Schumacher TN, Davis SJ, McMichael AJ, Gillespie GM: An early HIV mutation within an HLA-B*57-restricted T cell epitope abrogates binding to the killer inhibitory receptor 3DL1. J Virol. 2011, 85: 5415-5422. 10.1128/JVI.00238-11.
Tomiyama H, Sakaguchi T, Miwa K, Oka S, Iwamoto A, Kaneko Y, Takiguchi M: Identification of multiple HIV-1 CTL epitopes presented by HLA-B*5101 molecules. Hum Immunol. 1999, 60: 177-186. 10.1016/S0198-8859(98)00113-X.
Altfeld M, Kalife ET, Qi Y, Streeck H, Lichterfeld M, Johnston MN, Burgett N, Swartz ME, Yang A, Alter G, Yu XG, Meier A, Rockstroh JK, Allen TM, Jessen H, Rosenberg ES, Carrington M, Walker BD: 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.
Goulder PJ, Phillips RE, Colbert RA, McAdam S, Ogg G, Nowak MA, Giangrande P, Luzzi G, Morgan B, Edwards A, McMichael AJ, Rowland-Jones S: Late escape from an immunodominant cytotoxic T-lymphocyte response associated with progression to AIDS. Nat Med. 1997, 3: 212-217. 10.1038/nm0297-212.
Schneidewind A, Brockman MA, Yang R, Adam RI, Li B, Le Gall S, Rinaldo CR, Craggs SL, Allgaier RL, Power KA, Kuntzen T, Tung CS, LaBute MX, Mueller SM, Harrer T, McMichael AJ, Goulder PJ, Aiken C, Brander C, Kelleher AD, Allen TM: 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.
Herbeck JT, Rolland M, Liu Y, McLaughlin S, McNevin J, Zhao H, Wong K, Stoddard JN, Raugi D, Sorensen S, Genowati I, Birditt B, McKay A, Diem K, Maust BS, Deng W, Collier AC, Stekler JD, McElrath MJ, Mullins JI: Demographic processes affect HIV-1 evolution in primary infection before the onset of selective processes. J Virol. 2011, 85: 7523-7534. 10.1128/JVI.02697-10.
Brumme ZL, Brumme CJ, Heckerman D, Korber BT, Daniels M, Carlson J, Kadie C, Bhattacharya T, Chui C, Szinger J, Mo T, Hogg RS, Montaner JS, Frahm N, Brander C, Walker BD, Harrigan PR: 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.
O'Brien SJ, Gao X, Carrington M: HLA and AIDS: a cautionary tale. Trends Mol Med. 2001, 7: 379-381. 10.1016/S1471-4914(01)02131-1.
Liu Y, McNevin J, Cao J, Zhao H, Genowati I, Wong K, McLaughlin S, McSweyn MD, Diem K, Stevens CE, Maenza J, He H, Nickle DC, Shriner D, Holte SE, Collier AC, Corey L, McElrath MJ, Mullins JI: Selection on the human immunodeficiency virus type 1 proteome following primary infection. J Virol. 2006, 80: 9519-9529. 10.1128/JVI.00575-06.
Streeck H, Jolin JS, Qi Y, Yassine-Diab B, Johnson RC, Kwon DS, Addo MM, Brumme C, Routy JP, Little S, Jessen HK, Kelleher AD, Hecht FM, Sekaly RP, Rosenberg ES, Walker BD, Carrington M, Altfeld M: Human immunodeficiency virus type 1-specific CD8+ T-cell responses during primary infection are major determinants of the viral set point and loss of CD4+ T cells. J Virol. 2009, 83: 7641-7648. 10.1128/JVI.00182-09.
Prince JL, Claiborne DT, Carlson JM, Schaefer M, Yu T, Lahki S, Prentice HA, Yue L, Vishwanathan SA, Kilembe W, Goepfert P, Price MA, Gilmour J, Mulenga J, Farmer P, Derdeyn CA, Tang J, Heckerman D, Kaslow RA, Allen SA, Hunter E: Role of transmitted Gag CTL polymorphisms in defining replicative capacity and early HIV-1 pathogenesis. PLoS Pathog. 2012, 8: e1003041-10.1371/journal.ppat.1003041.
Song H, Pavlicek JW, Cai F, Bhattacharya T, Li H, Iyer SS, Bar KJ, Decker JM, Goonetilleke N, Liu MK, Berg A, Hora B, Drinker MS, Eudailey J, Pickeral J, Moody MA, Ferrari G, McMichael A, Perelson AS, Shaw GM, Hahn BH, Haynes BF, Gao F: Impact of immune escape mutations on HIV-1 fitness in the context of the cognate transmitted/founder genome. Retrovirology. 2012, 9: 89-10.1186/1742-4690-9-89.
Borghans JA, Molgaard A, de Boer RJ, Kesmir C: HLA alleles associated with slow progression to AIDS truly prefer to present HIV-1 p24. PLoS ONE. 2007, 2: e920-10.1371/journal.pone.0000920.
Streeck H, Lichterfeld M, Alter G, Meier A, Teigen N, Yassine-Diab B, Sidhu HK, Little S, Kelleher A, Routy JP, Rosenberg ES, Sekaly RP, Walker BD, Altfeld M: 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.
Brennan CA, Ibarrondo FJ, Sugar CA, Hausner MA, Shih R, Ng HL, Detels R, Margolick JB, Rinaldo CR, Phair J, Jacobson LP, Yang OO, Jamieson BD: Early HLA-B*57-restricted CD8+ T lymphocyte responses predict HIV-1 disease progression. J Virol. 2012, 86: 10505-10516. 10.1128/JVI.00102-12.
Brockman MA, Schneidewind A, Lahaie M, Schmidt A, Miura T, Desouza I, Ryvkin F, Derdeyn CA, Allen S, Hunter E, Mulenga J, Goepfert PA, Walker BD, Allen TM: Escape and compensation from early HLA-B57-mediated cytotoxic T-lymphocyte pressure on human immunodeficiency virus type 1 Gag alter capsid interactions with cyclophilin A. J Virol. 2007, 81: 12608-12618. 10.1128/JVI.01369-07.
Martinez-Picado J, Prado JG, Fry EE, Pfafferott K, Leslie A, Chetty S, Thobakgale C, Honeyborne I, Crawford H, Matthews P, Pillay T, Rousseau C, Mullins JI, Brander C, Walker BD, Stuart DI, Kiepiela P, Goulder P: Fitness cost of escape mutations in p24 Gag in association with control of human immunodeficiency virus type 1. J Virol. 2006, 80: 3617-3623. 10.1128/JVI.80.7.3617-3623.2006.
Crawford H, Lumm W, Leslie A, Schaefer M, Boeras D, Prado JG, Tang J, Farmer P, Ndung'u T, Lakhi S, Gilmour J, Goepfert P, Walker BD, Kaslow R, Mulenga J, Allen S, Goulder PJ, Hunter E: Evolution of HLA-B*5703 HIV-1 escape mutations in HLA-B*5703-positive individuals and their transmission recipients. J Exp Med. 2009, 206: 909-921.
Feeney ME, Tang Y, Roosevelt KA, Leslie AJ, McIntosh K, Karthas N, Walker BD, Goulder PJ: Immune escape precedes breakthrough human immunodeficiency virus type 1 viremia and broadening of the cytotoxic T-lymphocyte response in an HLA-B27-positive long-term-nonprogressing child. J Virol. 2004, 78: 8927-8930. 10.1128/JVI.78.16.8927-8930.2004.
Gao X, Bashirova A, Iversen AK, Phair J, Goedert JJ, Buchbinder S, Hoots K, Vlahov D, Altfeld M, O'Brien SJ, Carrington M: AIDS restriction HLA allotypes target distinct intervals of HIV-1 pathogenesis. Nat Med. 2005, 11: 1290-1292. 10.1038/nm1333.
Kelleher AD, Long C, Holmes EC, Allen RL, Wilson J, Conlon C, Workman C, Shaunak S, Olson K, Goulder P, Brander C, Ogg G, Sullivan JS, Dyer W, Jones I, McMichael AJ, Rowland-Jones S, Phillips RE: Clustered mutations in HIV-1 gag are consistently required for escape from HLA-B27-restricted cytotoxic T lymphocyte responses. J Exp Med. 2001, 193: 375-386. 10.1084/jem.193.3.375.
Montaner JS: Treatment as prevention: toward an AIDS-free generation. Top Antivir Med. 2013, 21: 110-114.
Woods CK, Brumme CJ, Liu TF, Chui CK, Chu AL, Wynhoven B, Hall TA, Trevino C, Shafer RW, Harrigan PR: Automating HIV drug resistance genotyping with RECall, a freely accessible sequence analysis tool. J Clin Microbiol. 2012, 50: 1936-1942. 10.1128/JCM.06689-11.
Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O: New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010, 59: 307-321. 10.1093/sysbio/syq010.
Fourment M, Gibbs MJ: PATRISTIC: a program for calculating patristic distances and graphically comparing the components of genetic change. BMC Evol Biol. 2006, 6: 1-10.1186/1471-2148-6-1.
Cotton LA, Rahman MA, Ng C, Le AQ, Milloy MJ, Mo T, Brumme ZL: HLA class I sequence-based typing using DNA recovered from frozen plasma. J Immunol Methods. 2012, 382: 40-47. 10.1016/j.jim.2012.05.003.
Listgarten J, Brumme Z, Kadie C, Xiaojiang G, Walker B, Carrington M, Goulder P, Heckerman D: Statistical resolution of ambiguous HLA typing data. PLoS Comput Biol. 2008, 4: e1000016-10.1371/journal.pcbi.1000016.
Storey JD, Tibshirani R: Statistical significance for genomewide studies. Proc Natl Acad Sci U S A. 2003, 100: 9440-9445. 10.1073/pnas.1530509100.
We thank Colin Shen and Zhixing (Samuel) Tan for laboratory assistance. We thank Chanson Brumme, Conan Woods and Daniel MacMillan for database assistance. We thank Richard Harrigan and Bruce D. Walker for data access and mentorship. We thank Mark Brockman for helpful discussions.
This work was supported by operating grants from the Canadian Institutes for Health Research (CIHR) MOP-93536 and HOP-115700 to ZLB. The VIDUS and ACCESS projects are funded by the National Institute on Drug Abuse, NIH (RO1DA011591 and RO1DA021525). This project has been funded in whole or in part with federal funds from the Frederick National Laboratory for Cancer Research, under Contract No. HHSN261200800001E. The content of this publication does not necessarily reflect the views or policies of the Department of Health and Human Services, nor does mention of trade names, commercial products, or organizations imply endorsement by the U.S. Government. This Research was supported in part by the Intramural Research Program of the NIH, Frederick National Lab, Center for Cancer Research. EM was supported by a Master’s Scholarship from the Canadian Association of HIV Research and Abbott Virology. AQL is the recipient of a CIHR Frederick Banting and Charles Best Masters award. DRC was the recipient of a CIHR CANADA-HOPE fellowship. M-JM is supported by post-doctoral fellowships from CIHR and the Michael Smith Foundation for Health Research (MSFHR). ZLB is a recipient of a CIHR New Investigator Award and a MSFHR Scholar Award.
The funding bodies played no role in the design, collection, analysis, or interpretation of data, nor in the writing of the manuscript or the decision to submit it for publication.
The authors declare that they have no competing interests.
EM and JMC designed and executed data analyses. AQL, DRC, RM, MAR, and CN performed HLA class I and/or plasma HIV-1 genotyping. HJ, ADK, MM, TMA, MJM, MC, and MAW contributed patient specimens and/or data. ZLB conceived and designed the study, analyzed data, and wrote the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: A full summary of HLA-associated polymorphisms, their HIV genomic locations and directions of association (adapted vs nonadapted), and their observed frequencies in early versus chronic infection in persons harboring vs. not harboring the restricting HLA allele are provided in columns A-L and O-T. P-values and q-values of association are also provided for early infection (Columns M-N) and chronic infection (Columns U-V). P- and q-values comparing the strength of selection of individual HLA-associated polymorphisms between early and chronic stages are provided in Columns W-X. Finally, all HLA-associated “adapted” polymorphisms are categorized with respect to their relative timescale of escape (early vs. later; Column Y) as well as their relative background frequency in the population (lower vs. higher; Column Z). (XLSX 87 KB)
About this article
Cite this article
Martin, E., Carlson, J.M., Le, A.Q. et al. Early immune adaptation in HIV-1 revealed by population-level approaches. Retrovirology 11, 64 (2014). https://doi.org/10.1186/s12977-014-0064-1
- Human immunodeficiency virus type-1 (HIV-1)
- Human leukocyte antigen (HLA) class I
- CD8+ cytotoxic T-lymphocytes (CTL)
- Immune escape
- HLA-associated polymorphism
- Acute/early infection
- Population-level analysis
- Statistical association with phylogenetic correction