Unique features of HLA-mediated HIV evolution in a Mexican cohort: a comparative study

Effect of Difference in Consensus Sequence between HIV-1 Subtype A/E and Subtype B Viruses on Elicitation of Gag-Specific CD8+ T Cells and Accumulation of HLA-Associated Escape Mutations

The Gag280 mutation is associated with HLA-C*01:02 but not with HLA-B*52:01 in subtype A/E-infected individuals, whereas this mutation is associated with HLA-B*52:01 but not with HLA-C*01:02 in subtype B infections. Although it is known that the Gag280 mutant is selected by HLA-B*52:01-restricted GagRI8 (Gag275-282)-specific T cells in subtype B infections, it remains unknown why this Gag280 mutation is associated with HLA-C*01:02 rather than HLA-B*52:01 in subtype A/E infections. The subtype B and A/E viruses have different consensus sequence, with Thr and Val at Gag280, respectively. To clarify the effect of this difference in Gag280 consensus sequence, we investigated the role of HLA-C*01:02-restricted GagYI9 (Gag277-285)-specific T cells in selection of Gag280 mutations in subtype A/E-infected Vietnamese and subtype B-infected Japanese individuals. GagYI9-4V-specific T cells, which were frequently elicited in Vietnamese individuals infected with the consensus-type A/E virus, failed to recognize GagV280T mutant A/E virus-infected cells. GagYI9-4T mutant epitope-specific T cells, which were weakly elicited in individuals infected with the mutant A/E virus, had weak or no ability to recognize the mutant virus. These results account for the mechanism for selection and accumulation of GagV280T mutants in the case of subtype A/E infections. In contrast, HLA-C*01:02-restricted GagYI9-4T-specific T cells were weakly elicited in Japanese individuals infected with the subtype B virus, explaining why HLA-C*01:02-restricted Gag280 mutations are not accumulated in the case of a subtype B infection. The present study demonstrated that a difference in the Gag280 consensus sequence influenced the elicitation of the GagYI9-specific T cells involved in the accumulation of HLA-C*01:02-associated Gag280 mutations.IMPORTANCE HIV-1 mutations escaped from HIV-specific CD8⁺ T cells are mostly detected as HLA-associated mutations. A diversity of HLA-associated mutations is somewhat distinct to each race and region, since HLA allele distribution differs among them. A difference in the consensus sequence among HIV-1 subtypes may also influence the diversity of HLA-associated mutations. HLA-C*01:02-associated GagV280T and HLA-B*52:01-associated GagT280A/S mutations were previously identified in HIV-1 subtype A/E-infected and subtype B-infected individuals, respectively, though these subtype viruses have a different consensus sequence at Gag280. We demonstrated that the GagV280T mutant virus was selected by HLA-C*01:02-restricted GagYI9-4V-specific T cells in subtype A/E-infected Vietnamese but that HLA-C*01:02-restricted GagYI9-4T-specific T cells were weakly elicited in subtype B-infected Japanese. Together with our recent study which demonstrated the mechanism for the accumulation of HLA-B*52:01-associated mutations, we clarified the mechanism for the accumulation of different Gag280 mutations and the effect of the difference in the consensus sequence on the accumulation of escape mutations.

T-cell responses to sequentially emerging viral escape mutants shape long-term HIV-1 population dynamics

HIV-1 strains harboring immune escape mutations can persist in circulation, but the impact of selection by multiple HLA alleles on population HIV-1 dynamics remains unclear. In Japan, HIV-1 Reverse Transcriptase codon 135 (RT135) is under strong immune pressure by HLA-B*51:01-restricted and HLA-B*52:01-restricted T cells that target a key epitope in this region (TI8; spanning RT codons 128-135). Major population-level shifts have occurred at HIV-1 RT135 during the Japanese epidemic, which first affected hemophiliacs (via imported contaminated blood products) and subsequently non-hemophiliacs (via domestic transmission). Specifically, threonine accumulated at RT135 (RT135T) in hemophiliac and non-hemophiliac HLA-B*51:01+ individuals diagnosed before 1997, but since then RT135T has markedly declined while RT135L has increased among non-hemophiliac individuals. We demonstrated that RT135V selection by HLA-B*52:01-restricted TI8-specific T-cells led to the creation of a new HLA-C*12:02-restricted epitope TN9-8V. We further showed that TN9-8V-specific HLA-C*12:02-restricted T cells selected RT135L while TN9-8T-specific HLA-C*12:02-restricted T cells suppressed replication of the RT135T variant. Thus, population-level accumulation of the RT135L mutation over time in Japan can be explained by initial targeting of the TI8 epitope by HLA-B*52:01-restricted T-cells, followed by targeting of the resulting escape mutant by HLA-C*12:02-restricted T-cells. We further demonstrate that this phenomenon is particular to Japan, where the HLA-B*52:01-C*12:02 haplotype is common: RT135L did not accumulate over a 15-year longitudinal analysis of HIV sequences in British Columbia, Canada, where this haplotype is rare. Together, our observations reveal that T-cell responses to sequentially emerging viral escape mutants can shape long-term HIV-1 population dynamics in a host population-specific manner.

High-Resolution HLA Typing of HLA-A, -B, -C, -DRB1, and -DQB1 in Kinh Vietnamese by Using Next-Generation Sequencing

Human leukocyte antigen (HLA) genotyping displays the particular characteristics of HLA alleles and haplotype frequencies in each population. Although it is considered the current gold standard for HLA typing, high-resolution sequence-based HLA typing is currently unavailable in Kinh Vietnamese populations. In this study, high-resolution sequence-based HLA typing (3-field) was performed using an amplicon-based next-generation sequencing platform to identify the HLA-A, -B, -C, -DRB1, and -DQB1 alleles of 101 unrelated healthy Kinh Vietnamese individuals from southern Vietnam. A total of 28 HLA-A, 41 HLA-B, 21 HLA-C, 26 HLA-DRB1, and 25 HLA-DQB1 alleles were identified. The most frequently occurring HLA alleles were A^∗11:01:01, B^∗15:02:01, C^∗07:02:01, DRB1^∗12:02:01, and DQB1^∗03:01:01. Haplotype calculation showed that A^∗29:01:01∼B^∗07:05:01, DRB1^∗12:02:01∼DQB1^∗3:01:01, A^∗29:01:01∼C^∗15:05:02∼B^∗07:05:01, A^∗33:03:01∼B^∗58:01:01∼DRB1^∗03:01:01, and A^∗29:01:01∼C^∗15:05:02∼B^∗07:05:01∼DRB1^∗10:01:01∼DQB1^∗05:01:01 were the most common haplotypes in the southern Kinh Vietnamese population. Allele distribution and haplotype analyses demonstrated that the Vietnamese population shares HLA features with South-East Asians but retains unique characteristics. Data from this study will be potentially applicable in medicine and anthropology.

Clinical and evolutionary consequences of HIV adaptation to HLA: implications for vaccine and cure

PURPOSE OF REVIEW

The purpose of this review is to summarize recent advances in our understanding of HIV adaptation to human leukocyte antigen (HLA)-associated immune pressures and its relevance to HIV prevention and cure research.

RECENT FINDINGS

Recent research has confirmed that HLA is a major driver of individual and population-level HIV evolution, that HIV strains are adapting to the immunogenetic profiles of the different human ethnic groups in which they circulate, and that HIV adaptation has substantial clinical and immunologic consequences. As such, adaptation represents a major challenge to HIV prevention and cure. At the same time, there are opportunities: Studies of HIV adaptation are revealing why certain HLA alleles are protective in some populations and not others; they are identifying immunogenic viral epitopes that harbor high mutational barriers to escape, and they may help illuminate novel, vaccine-relevant HIV epitopes in regions where circulating adaptation is extensive. Elucidation of HLA-driven adapted and nonadapted viral forms in different human populations and HIV subtypes also renders 'personalized' immunogen selection, as a component of HIV cure strategies, conceptually feasible.

SUMMARY

Though adaptation represents a major challenge to HIV prevention and cure, achieving an in-depth understanding of this phenomenon can help move the design of such strategies forward.

Genotypic and Mechanistic Characterization of Subtype-Specific HIV Adaptation to Host Cellular Immunity

The extent to which viral genetic context influences HIV adaptation to human leukocyte antigen (HLA) class I-restricted immune pressures remains incompletely understood. The Ugandan HIV epidemic, where major pandemic group M subtypes A1 and D cocirculate in a single host population, provides an opportunity to investigate this question. We characterized plasma HIV RNA gag, pol, and nef sequences, along with host HLA genotypes, in 464 antiretroviral-naive individuals chronically infected with HIV subtype A1 or D. Using phylogenetically informed statistical approaches, we identified HLA-associated polymorphisms and formally compared their strengths of selection between viral subtypes. A substantial number (32%) of HLA-associated polymorphisms identified in subtype A1 and/or D had previously been reported in subtype B, C, and/or circulating recombinant form 01_AE (CRF01_AE), confirming the shared nature of many HLA-driven escape pathways regardless of viral genetic context. Nevertheless, 34% of the identified HLA-associated polymorphisms were significantly differentially selected between subtypes A1 and D. Experimental investigation of select examples of subtype-specific escape revealed distinct underlying mechanisms with important implications for vaccine design: whereas some were attributable to subtype-specific sequence variation that influenced epitope-HLA binding, others were attributable to differential mutational barriers to immune escape. Overall, our results confirm that HIV genetic context is a key modulator of viral adaptation to host cellular immunity and highlight the power of combined bioinformatic and mechanistic studies, paired with knowledge of epitope immunogenicity, to identify appropriate viral regions for inclusion in subtype-specific and universal HIV vaccine strategies.IMPORTANCE The identification of HIV polymorphisms reproducibly selected under pressure by specific HLA alleles and the elucidation of their impact on viral function can help identify immunogenic viral regions where immune escape incurs a fitness cost. However, our knowledge of HLA-driven escape pathways and their functional costs is largely limited to HIV subtype B and, to a lesser extent, subtype C. Our study represents the first characterization of HLA-driven adaptation pathways in HIV subtypes A1 and D, which dominate in East Africa, and the first statistically rigorous characterization of differential HLA-driven escape across viral subtypes. The results support a considerable impact of viral genetic context on HIV adaptation to host HLA, where HIV subtype-specific sequence variation influences both epitope-HLA binding and the fitness costs of escape. Integrated bioinformatic and mechanistic characterization of these and other instances of differential escape could aid rational cytotoxic T-lymphocyte-based vaccine immunogen selection for both subtype-specific and universal HIV vaccines.

Novel HLA class I associations with HIV-1 control in a unique genetically admixed population

Associations between HLA class I alleles and HIV progression in populations exhibiting Amerindian and Caucasian genetic admixture remain understudied. Using univariable and multivariable analyses we evaluated HLA associations with five HIV clinical parameters in 3,213 HIV clade B-infected, ART-naïve individuals from Mexico and Central America (MEX/CAM cohort). A Canadian cohort (HOMER, n = 1622) was used for comparison. As expected, HLA allele frequencies in MEX/CAM and HOMER differed markedly. In MEX/CAM, 13 HLA-A, 24 HLA-B, and 14 HLA-C alleles were significantly associated with at least one clinical parameter. These included previously described protective (e.g. B*27:05, B*57:01/02/03 and B*58:01) and risk (e.g. B*35:02) alleles, as well as novel ones (e.g. A*03:01, B*15:39 and B*39:02 identified as protective, and A*68:03/05, B*15:30, B*35:12/14, B*39:01/06, B*39:05~C*07:02, and B*40:01~C*03:04 identified as risk). Interestingly, both protective (e.g. B*39:02) and risk (e.g. B*39:01/05/06) subtypes were identified within the common and genetically diverse HLA-B*39 allele group, characteristic to Amerindian populations. While HLA-HIV associations identified in MEX and CAM separately were similar overall (Spearman's rho = 0.33, p = 0.03), region-specific associations were also noted. The identification of both canonical and novel HLA/HIV associations provides a first step towards improved understanding of HIV immune control among unique and understudied Mestizo populations.

Weaker HLA Footprints on HIV in the Unique and Highly Genetically Admixed Host Population of Mexico

HIV circumvents HLA class I-restricted CD8⁺ T-cell responses through selection of escape mutations that leave characteristic mutational “footprints,” also known as HLA-associated polymorphisms (HAPs), on HIV sequences at the population level. While many HLA footprints are universal across HIV subtypes and human populations, others can be region specific as a result of the unique immunogenetic background of each host population. Using a published probabilistic phylogenetically informed model, we compared HAPs in HIV Gag and Pol (PR-RT) in 1,612 subtype B-infected, antiretroviral treatment-naive individuals from Mexico and 1,641 individuals from Canada/United States. A total of 252 HLA class I allele subtypes were represented, including 140 observed in both cohorts, 67 unique to Mexico, and 45 unique to Canada/United States. At the predefined statistical threshold of a q value of <0.2, 358 HAPs (201 in Gag, 157 in PR-RT) were identified in Mexico, while 905 (534 in Gag and 371 in PR-RT) were identified in Canada/United States. HAPs identified in Mexico included both canonical HLA-associated escape pathways and novel associations, in particular with HLA alleles enriched in Amerindian and mestizo populations. Remarkably, HLA footprints on HIV in Mexico were not only fewer but also, on average, significantly weaker than those in Canada/United States, although some exceptions were noted. Moreover, exploratory analyses suggested that the weaker HLA footprint on HIV in Mexico may be due, at least in part, to weaker and/or less reproducible HLA-mediated immune pressures on HIV in this population. The implications of these differences for natural and vaccine-induced anti-HIV immunity merit further investigation.

IMPORTANCE HLA footprints on HIV identify viral regions under intense and consistent pressure by HLA-restricted immune responses and the common mutational pathways that HIV uses to evade them. In particular, HLA footprints can identify novel immunogenic regions and/or epitopes targeted by understudied HLA alleles; moreover, comparative analyses across immunogenetically distinct populations can illuminate the extent to which HIV immunogenic regions and escape pathways are shared versus population-specific pathways, information which can in turn inform the design of universal or geographically tailored HIV vaccines. We compared HLA-associated footprints on HIV in two immunogenetically distinct North American populations, those of Mexico and Canada/United States. We identify both shared and population-specific pathways of HIV adaptation but also make the surprising observation that HLA footprints on HIV in Mexico overall are fewer and weaker than those in Canada/United States, raising the possibility that HLA-restricted antiviral immune responses in Mexico are weaker, and/or escape pathways somewhat less consistent, than those in other populations.

Accumulation of Pol Mutations Selected by HLA-B*52:01-C*12:02 Protective Haplotype-Restricted Cytotoxic T Lymphocytes Causes Low Plasma Viral Load Due to Low Viral Fitness of Mutant Viruses

HLA-B*52:01-C*12:02, which is the most abundant haplotype in Japan, has a protective effect on disease progression in HIV-1-infected Japanese individuals, whereas HLA-B*57 and -B*27 protective alleles are very rare in Japan. A previous study on HLA-associated polymorphisms demonstrated that the number of HLA-B*52:01-associated mutations at four Pol positions was inversely correlated with plasma viral load (pVL) in HLA-B*52:01-negative individuals, suggesting that the transmission of HIV-1 with these mutations could modulate the pVL in the population. However, it remains unknown whether these mutations were selected by HLA-B*52:01-restricted CTLs and also reduced viral fitness. In this study, we identified two HLA-B*52:01-restricted and one HLA-C*12:02-restricted novel cytotoxic T-lymphocyte (CTL) epitopes in Pol. Analysis using CTLs specific for these three epitopes demonstrated that these CTLs failed to recognize mutant epitopes or more weakly recognized cells infected with mutant viruses than wild-type virus, supporting the idea that these mutations were selected by the HLA-B*52:01- or HLA-C*12:02-restricted T cells. We further showed that these mutations reduced viral fitness, although the effect of each mutation was weak. The present study demonstrated that the accumulation of these Pol mutations selected by HLA-B*52:01- or HLA-C*12:02-restricted CTLs impaired viral replication capacity and thus reduced the pVL. The fitness cost imposed by the mutations partially accounted for the effect of the HLA-B*52:01-C*12:02 haplotype on clinical outcome, together with the effect of HLA-B*52:01-restricted CTLs on viral replication, which had been previously demonstrated.

IMPORTANCE

Numerous population-based studies identified HLA-associated HIV-1 mutations to predict HIV-1 escape mutations from cytotoxic T lymphocytes (CTLs). However, the majority of these HLA-associated mutations have not been identified as CTL escape mutations. Our previous population-based study showed that five HLA-B*52:01-associated mutations at four Pol positions were inversely correlated with the plasma viral load in HLA-B*52:01-negative Japanese individuals. In the present study, we demonstrated that these mutations were indeed selected by CTLs specific for novel B*52:01- and C*12:02-restricted epitopes and that the accumulation of these mutations reduced the viral fitness in vitro This study elucidated the mechanism by which the accumulation of these CTL escape mutations contributed to the protective effect of the HLA-B*52:01-HLA-C*12:02 haplotype on disease progression in HIV-1-infected Japanese individuals.

Balance between transmitted HLA preadapted and nonassociated polymorphisms is a major determinant of HIV-1 disease progression

In a cohort of Zambian heterosexual transmission pairs, the authors show that HIV-1–transmitted variants already exhibit a significant degree of preadaptation to the new host's HLA alleles, which, modulated by polymorphisms that decrease viral fitness, determines early set-point VL and the rate of disease progression in the newly infected individual.

HIV-1 adapts to a new host through mutations that facilitate immune escape. Here, we evaluate the impact on viral control and disease progression of transmitted polymorphisms that were either preadapted to or nonassociated with the new host’s HLA. In a cohort of 169 Zambian heterosexual transmission pairs, we found that almost one-third of possible HLA-linked target sites in the transmitted virus Gag protein are already adapted, and that this transmitted preadaptation significantly reduced early immune recognition of epitopes. Transmitted preadapted and nonassociated polymorphisms showed opposing effects on set-point VL and the balance between the two was significantly associated with higher set-point VLs in a multivariable model including other risk factors. Transmitted preadaptation was also significantly associated with faster CD4 decline (<350 cells/µl) and this association was stronger after accounting for nonassociated polymorphisms, which were linked with slower CD4 decline. Overall, the relative ratio of the two classes of polymorphisms was found to be the major determinant of CD4 decline in a multivariable model including other risk factors. This study reveals that, even before an immune response is mounted in the new host, the balance of these opposing factors can significantly influence the outcome of HIV-1 infection.

A strong association of human leukocyte antigen-associated Pol and Gag mutations with clinical parameters in HIV-1 subtype A/E infection

OBJECTIVES

Identification of human leukocyte antigen-associated HIV-1 polymorphisms (HLA-APs) in different global populations furthers our understanding of HIV-1 pathogenesis and may help identify candidate immunogens for HIV vaccines targeted to these populations. Although numerous population-based studies identifying HLA-APs have been conducted in HIV-1 subtype B- and subtype C-infected cohorts, few have focused on subtype A/E.

DESIGN

We investigated HLA-APs in a cohort of chronically HIV-1 subtype A/E-infected Vietnamese individuals.

METHODS

HLA-APs in HIV-1 Gag, Pol, and Nef regions from 388 treatment-naive individuals chronically infected with HIV-1 subtype A/E were analyzed using phylogenetically informed approaches.

RESULTS

A total of 303 HLA-APs were identified. HLA-APs occurring at six positions in Gag and six positions in Pol were significantly associated with higher plasma viral load (pVL), whereas HLA-APs occurring at two positions in Gag and 13 positions in Pol were significantly associated with lower CD4 T-cell counts. Furthermore, the proportion of Pol codons harboring an HLA-AP specific to the host's HLA correlated positively with HIV-1 pVL (R = 0.22; P < 0.0001) and inversely with CD4 T-cell counts (R = -0.32; P < 0.0001). Similarly, the proportion of HLA-associated Gag codons harboring host-specific HLA-AP correlated inversely with CD4 T-cell counts (R = -0.13; P = 0.01).

CONCLUSION

These significant associations between HIV-1 amino acids adapted to Vietnamese HLA alleles and higher pVL and lower CD4 T-cell counts suggests that accumulation of cytotoxic T cells escape mutations may influence clinical outcomes in HIV-1 subtype A/E-infected Vietnamese individuals.

Identification of Immunogenic Cytotoxic T Lymphocyte Epitopes Containing Drug Resistance Mutations in Antiretroviral Treatment-Naïve HIV-Infected Individuals

Background

Therapeutic HIV vaccines may prove helpful to intensify antiretroviral treatment (ART) efficacy and may be an integral part of future cure strategies.

Methods

We examined IFN-gamma ELISpot responses to a panel of 218 HIV clade B consensus-based HIV protease-reverse transcriptase peptides, designed to mimic previously described and predicted cytotoxic T lymphocyte epitopes overlapping drug resistance (DR) positions, that either included the consensus sequence or the DR variant sequence, in 49 ART-naïve HIV-infected individuals. Next generation sequencing was used to assess the presence of minority DR variants in circulating viral populations.

Results

Although a wide spectrum of differential magnitudes of response to DR vs. WT peptide pairs was observed, responses to DR peptides were frequent and strong in the study cohort. No difference between the median magnitudes of response to DR vs. WT peptides was observed. Interestingly, of the 22 peptides that were recognized by >15% of the participants, two-thirds (64%) corresponded to DR peptides. When analysing responses per peptide pair per individual, responses to only WT (median 4 pairs/individual) or DR (median 6 pairs/individual) were more common than responses to both WT and DR (median 2 pairs/individual; p<0.001). While the presence of ELISpot responses to WT peptides was frequently associated with the presence of the corresponding peptide sequence in the patient’s virus (mean 68% of cases), responses to DR peptides were generally not associated with the presence of DR mutations in the viral population, even at low frequencies (mean 1.4% of cases; p = 0.0002).

Conclusions

Our data suggests that DR peptides are frequently immunogenic and raises the potential benefit of broadening the antigens included in a therapeutic vaccine approach to immunogenic epitopes containing common DR sequences. Further studies are needed to assess the quality of responses elicited by DR peptides.

HLA-B allele and haplotype diversity among Thai patients identified by PCR-SSOP: evidence for high risk of drug-induced hypersensitivity

Background: There are 3 classes of HLA molecules; HLA class I, II and III, of which different classes have different functions. HLA-B gene which belongs to HLA class I play an important role predicting drug hypersensitivity.

Materials and Methods: Nine hundred and eighty-six Thai subjects who registered at a pharmacogenomics laboratory were determined for HLA-B genotype using a two-stage sequence-specific oligonucleotide probe system (PCR-SSOP).

Results: In this study, HLA-B alleles did not deviate from Hardy-Weinberg equilibrium (P > 0.05). The most common HLA-B alleles observed in this population were HLA-B^*46:01 (11.51%), HLA-B^*58:01 (8.62%), HLA-B^*40:01 (8.22%), HLA-B^*15:02 (8.16%) and HLA-B^*13:01 (6.95%). This finding revealed that HLA-B allele frequency in the Thai population was consistent with the Chinese population (p > 0.05), however, differed from the Malaysian population (p < 0.05). The top five HLA-B genotypes were HLA-B^*40:01/46:01 (2.13%), HLA-B^*46:01/46:01 (2.03%), HLA-B^*40:01/58:01 (2.03%), HLA-B^*46:01/58:01 (1.93%) and HLA-B^*15:02/46:01 (1.83%). This study found that 15.92% of Thai subjects carry HLA-B^*15:02, which has been associated with carbamazepine-induced severe cutaneous adverse drug reactions (SCARs). Moreover, 16.33% of Thai subjects carry the HLA-B^*58:01 allele, which has been associated with allopurinol-induced SCARs.

Conclusion: This study demonstrates a high diversity of HLA-B polymorphisms in this Thai population. The high frequency of HLA-B pharmacogenomic markers in the population emphasizes the importance of such screening to predict/avoid drug hypersensitivity.

HIV-1 adaptation to HLA: a window into virus-host immune interactions

HIV-1 develops specific mutations within its genome that allow it to escape detection by human leukocyte antigen (HLA) class I-restricted immune responses, notably those of CD8(+) cytotoxic T lymphocytes (CTL). HLA thus represents a major force driving the evolution and diversification of HIV-1 within individuals and at the population level. Importantly, the study of HIV-1 adaptation to HLA also represents an opportunity to identify what qualities constitute an effective immune response, how the virus in turn adapts to these pressures, and how we may harness this information to design HIV-1 vaccines that stimulate effective cellular immunity.

Impact of HLA selection pressure on HIV fitness at a population level in Mexico and Barbados

Previous studies have demonstrated that effective cytotoxic T lymphocyte (CTL) responses drive the selection of escape mutations that reduce viral replication capacity (VRC). Escape mutations, including those with reduced VRC, can be transmitted and accumulate in a population. Here we compared two antiretroviral therapy (ART)-naive HIV clade B-infected cohorts, in Mexico and Barbados, in which the most protective HLA alleles (HLA-B*27/57/58:01/81:01) are differentially expressed, at 8% and 34%, respectively. Viral loads were significantly higher in Mexico than in Barbados (median, 40,774 versus 14,200; P < 0.0001), and absolute CD4⁺ T-cell counts were somewhat lower (median, 380/mm³ versus 403/mm³; P = 0.007). We tested the hypothesis that the disparate frequencies of these protective HLA alleles would be associated with a higher VRC at the population level in Mexico. Analysis of VRC in subjects in each cohort, matched for CD4⁺ T-cell count, revealed that the VRC was indeed higher in the Mexican cohort (mean, 1.13 versus 1.03; P = 0.0025). Although CD4 counts were matched, viral loads remained significantly higher in the Mexican subjects (P = 0.04). This VRC difference was reflected by a significantly higher frequency in the Barbados cohort of HLA-B*27/57/58:01/81:01-associated Gag escape mutations previously shown to incur a fitness cost on the virus (P = 0.004), a difference between the two cohorts that remained statistically significant even in subjects not expressing these protective alleles (P = 0.01). These data suggest that viral set points and disease progression rates at the population level may be significantly influenced by the prevalence of protective HLA alleles such as HLA-B*27/57/58:01/81:01 and that CD4 count-based guidelines to initiate antiretroviral therapy may need to be modified accordingly, to optimize the effectiveness of treatment-for-prevention strategies and reduce HIV transmission rates to the absolute minimum.

IMPORTANCE Immune control of HIV at an individual level is strongly influenced by the HLA class I genotype. HLA class I molecules mediating effective immune control, such as HLA-B*27 and HLA-B*57, are associated with the selection of escape mutants that reduce viral replicative capacity. The escape mutants selected in infected patients can be transmitted and affect the viral load and CD4 count in the recipient. These findings prompt the hypothesis that the frequency of protective alleles in a population may affect viral set points and rates of disease progression in that population. These studies in Mexico and Barbados, where the prevalence rates of protective HLA alleles are 8% and 34%, respectively, support this hypothesis. These data suggest that antiretroviral therapy (ART) treatment-for-prevention strategies will be less successful in populations such as those in Mexico, where viral loads are higher for a given CD4 count. Consideration may therefore usefully be given to ART initiation at higher absolute CD4 counts in such populations to optimize the impact of ART for prevention.

Host-specific adaptation of HIV-1 subtype B in the Japanese population

The extent to which HIV-1 clade B strains exhibit population-specific adaptations to host HLA alleles remains incompletely known, in part due to incomplete characterization of HLA-associated HIV-1 polymorphisms (HLA-APs) in different global populations. Moreover, it remains unknown to what extent the same HLA alleles may drive significantly different escape pathways across populations. As the Japanese population exhibits distinctive HLA class I allele distributions, comparative analysis of HLA-APs between HIV-1 clade B-infected Japanese and non-Asian cohorts could shed light on these questions. However, HLA-APs remain incompletely mapped in Japan. In a cohort of 430 treatment-naive Japanese with chronic HIV-1 clade B infection, we identified 284 HLA-APs in Gag, Pol, and Nef using phylogenetically corrected methods. The number of HLA-associated substitutions in Pol, notably those restricted by HLA-B*52:01, was weakly inversely correlated with the plasma viral load (pVL), suggesting that the transmission and persistence of B*52:01-driven Pol mutations could modulate the pVL. Differential selection of HLA-APs between HLA subtype members, including those differing only with respect to substitutions outside the peptide-binding groove, was observed, meriting further investigation as to their mechanisms of selection. Notably, two-thirds of HLA-APs identified in Japan had not been reported in previous studies of predominantly Caucasian cohorts and were attributable to HLA alleles unique to, or enriched in, Japan. We also identified 71 cases where the same HLA allele drove significantly different escape pathways in Japan versus predominantly Caucasian cohorts. Our results underscore the distinct global evolution of HIV-1 clade B as a result of host population-specific cellular immune pressures.

IMPORTANCE

Cytotoxic T lymphocyte (CTL) escape mutations in HIV-1 are broadly predictable based on the HLA class I alleles expressed by the host. Because HLA allele distributions differ among worldwide populations, the pattern and diversity of HLA-associated escape mutations are likely to be somewhat distinct to each race and region. HLA-associated polymorphisms (HLA-APs) in HIV-1 have previously been identified at the population level in European, North American, Australian, and African cohorts; however, large-scale analyses of HIV-1 clade B-specific HLA-APs in Asians are lacking. Differential intraclade HIV-1 adaptation to global populations can be investigated via comparative analyses of HLA-associated polymorphisms across ethnic groups, but such studies are rare. Here, we identify HLA-APs in a large Japanese HIV-1 clade B cohort using phylogenetically informed methods and observe that the majority of them had not been previously characterized in predominantly Caucasian populations. The results highlight HIV's unique adaptation to cellular immune pressures imposed by different global populations.

Distinct HIV-1 escape patterns selected by cytotoxic T cells with identical epitope specificity

Pol283-8-specific, HLA-B*51:01-restricted, cytotoxic T cells (CTLs) play a critical role in the long-term control of HIV-1 infection. However, these CTLs select for the reverse transcriptase (RT) I135X escape mutation, which may be accumulating in circulating HIV-1 sequences. We investigated the selection of the I135X mutation by CTLs specific for the same epitope but restricted by HLA-B*52:01. We found that Pol283-8-specific, HLA-B*52:01-restricted CTLs were elicited predominantly in chronically HIV-1-infected individuals. These CTLs had a strong ability to suppress the replication of wild-type HIV-1, though this ability was weaker than that of HLA-B*51:01-restricted CTLs. The crystal structure of the HLA-B*52:01-Pol283-8 peptide complex provided clear evidence that HLA-B*52:01 presents the peptide similarly to HLA-B*51:01, ensuring the cross-presentation of this epitope by both alleles. Population level analyses revealed a strong association of HLA-B*51:01 with the I135T mutant and a relatively weaker association of HLA-B*52:01 with several I135X mutants in both Japanese and predominantly Caucasian cohorts. An in vitro viral suppression assay revealed that the HLA-B*52:01-restricted CTLs failed to suppress the replication of the I135X mutant viruses, indicating the selection of these mutants by the CTLs. These results suggest that the different pattern of I135X mutant selection may have resulted from the difference between these two CTLs in the ability to suppress HIV-1 replication.

Correlates of protective cellular immunity revealed by analysis of population-level immune escape pathways in HIV-1

HLA class I-associated polymorphisms identified at the population level mark viral sites under immune pressure by individual HLA alleles. As such, analysis of their distribution, frequency, location, statistical strength, sequence conservation, and other properties offers a unique perspective from which to identify correlates of protective cellular immunity. We analyzed HLA-associated HIV-1 subtype B polymorphisms in 1,888 treatment-naïve, chronically infected individuals using phylogenetically informed methods and identified characteristics of HLA-associated immune pressures that differentiate protective and nonprotective alleles. Over 2,100 HLA-associated HIV-1 polymorphisms were identified, approximately one-third of which occurred inside or within 3 residues of an optimally defined cytotoxic T-lymphocyte (CTL) epitope. Differential CTL escape patterns between closely related HLA alleles were common and increased with greater evolutionary distance between allele group members. Among 9-mer epitopes, mutations at HLA-specific anchor residues represented the most frequently detected escape type: these occurred nearly 2-fold more frequently than expected by chance and were computationally predicted to reduce peptide-HLA binding nearly 10-fold on average. Characteristics associated with protective HLA alleles (defined using hazard ratios for progression to AIDS from natural history cohorts) included the potential to mount broad immune selection pressures across all HIV-1 proteins except Nef, the tendency to drive multisite and/or anchor residue escape mutations within known CTL epitopes, and the ability to strongly select mutations in conserved regions within HIV's structural and functional proteins. Thus, the factors defining protective cellular immune responses may be more complex than simply targeting conserved viral regions. The results provide new information to guide vaccine design and immunogenicity studies.

Intersubtype differences in the effect of a rare p24 gag mutation on HIV-1 replicative fitness

Certain immune-driven mutations in HIV-1, such as those arising in p24(Gag), decrease viral replicative capacity. However, the intersubtype differences in the replicative consequences of such mutations have not been explored. In HIV-1 subtype B, the p24(Gag) M250I mutation is a rare variant (0.6%) that is enriched among elite controllers (7.2%) (P = 0.0005) and appears to be a rare escape variant selected by HLA-B58 supertype alleles (P < 0.01). In contrast, in subtype C, it is a relatively common minor polymorphic variant (10 to 15%) whose appearance is not associated with a particular HLA allele. Using site-directed mutant viruses, we demonstrate that M250I reduces in vitro viral replicative capacity in both subtype B and subtype C sequences. However, whereas in subtype C downstream compensatory mutations at p24(Gag) codons 252 and 260 reduce the adverse effects of M250I, fitness costs in subtype B appear difficult to restore. Indeed, patient-derived subtype B sequences harboring M250I exhibited in vitro replicative defects, while those from subtype C did not. The structural implications of M250I were predicted by protein modeling to be greater in subtype B versus C, providing a potential explanation for its lower frequency and enhanced replicative defects in subtype B. In addition to accounting for genetic differences between HIV-1 subtypes, the design of cytotoxic-T-lymphocyte-based vaccines may need to account for differential effects of host-driven viral evolution on viral fitness.

Adaptive evolution of HIV at HLA epitopes is associated with ethnicity in Canada

Host immune selection pressure influences the development of mutations that allow for HIV escape. Mutation patterns induced in HIV by the human leukocyte antigen (HLA) are HLA-allele specific. As ethnic groups have distinct and characteristic HLA allele frequencies, we can expect divergent viral evolution within ethnicities. Here, we have sequenced and analyzed the HIV pol gene from 1248 subtype B infected, treatment-naïve individuals in Canada. Phylogenetic analysis showed no separation between pol sequences from five self-identified ethnic groups, yet fixation index (F_ST) values showed significant divergence between ethnicities. A total of 17 amino acid sites showed an ethnic-specific fixation pattern (0.015<F_ST <0.060, p<0.01), and 27 codons were inferred to be under positive selection (p<0.01), with each set of sites strongly associated with HLA sites (p = 1.78×10⁻⁶ and p = 1.91×10⁻⁷, respectively). Within the pol gene, eight sites under HLA selective pressure were correlated with ethnicity, indicating ‘adaptive divergence’ between the groups studied. Our findings highlight challenges in HIV vaccine design in ethnically diverse countries with subtype B epidemics.

Widespread impact of HLA restriction on immune control and escape pathways of HIV-1

The promiscuous presentation of epitopes by similar HLA class I alleles holds promise for a universal T-cell-based HIV-1 vaccine. However, in some instances, cytotoxic T lymphocytes (CTL) restricted by HLA alleles with similar or identical binding motifs are known to target epitopes at different frequencies, with different functional avidities and with different apparent clinical outcomes. Such differences may be illuminated by the association of similar HLA alleles with distinctive escape pathways. Using a novel computational method featuring phylogenetically corrected odds ratios, we systematically analyzed differential patterns of immune escape across all optimally defined epitopes in Gag, Pol, and Nef in 2,126 HIV-1 clade C-infected adults. Overall, we identified 301 polymorphisms in 90 epitopes associated with HLA alleles belonging to shared supertypes. We detected differential escape in 37 of 38 epitopes restricted by more than one allele, which included 278 instances of differential escape at the polymorphism level. The majority (66 to 97%) of these resulted from the selection of unique HLA-specific polymorphisms rather than differential epitope targeting rates, as confirmed by gamma interferon (IFN-γ) enzyme-linked immunosorbent spot assay (ELISPOT) data. Discordant associations between HLA alleles and viral load were frequently observed between allele pairs that selected for differential escape. Furthermore, the total number of associated polymorphisms strongly correlated with average viral load. These studies confirm that differential escape is a widespread phenomenon and may be the norm when two alleles present the same epitope. Given the clinical correlates of immune escape, such heterogeneity suggests that certain epitopes will lead to discordant outcomes if applied universally in a vaccine.

High-avidity, high-IFNγ-producing CD8 T-cell responses following immune selection during HIV-1 infection

HIV-1 mutations, which reduce or abolish CTL responses against virus-infected cells, are frequently selected in acute and chronic HIV infection. Among population HIV-1 sequences, immune selection is evident as human leukocyte antigen (HLA) allele-associated substitutions of amino acids within or near CD8 T-cell epitopes. In these cases, the non-adapted epitope is susceptible to immune recognition until an escape mutation renders the epitope less immunogenic. However, several population-based studies have independently identified HLA-associated viral changes, which lead to the formation of a new T-cell epitope, suggesting that the immune responses that these variants or 'neo-epitopes' elicit provide an evolutionary advantage to the virus rather than the host. Here, we examined the functional characteristics of eight CD8 T-cell responses that result from viral adaptation in 125 HLA-genotyped individuals with chronic HIV-1 infection. Neo-epitopes included well-characterized immunodominant epitopes restricted by common HLA alleles, and in most cases the T-cell responses against the neo-epitope showed significantly greater functional avidity and higher IFNγ production than T cells for non-adapted epitopes, but were not more cytotoxic. Neo-epitope formation and emergence of cognate T-cell response coincident with a rise in viral load was then observed in vivo in an acutely infected individual. These findings show that HIV-1 adaptation not only abrogates the immune recognition of early targeted epitopes, but may also increase immune recognition to other epitopes, which elicit immunodominant but non-protective T-cell responses. These data have implications for immunodominance associated with polyvalent vaccines based on the diversity of chronic HIV-1 sequences.

Translation of HLA-HIV associations to the cellular level: HIV adapts to inflate CD8 T cell responses against Nef and HLA-adapted variant epitopes

Strong statistical associations between polymorphisms in HIV-1 population sequences and carriage of HLA class I alleles have been widely used to identify possible sites of CD8 T cell immune selection in vivo. However, there have been few attempts to prospectively and systematically test these genetic hypotheses arising from population-based studies at a cellular, functional level. We assayed CD8 T cell epitope-specific IFN-γ responses in 290 individuals from the same cohort, which gave rise to 874 HLA-HIV associations in genetic analyses, taking into account autologous viral sequences and individual HLA genotypes. We found immunological evidence for 58% of 374 associations tested as sites of primary immune selection and identified up to 50 novel HIV-1 epitopes using this reverse-genomics approach. Many HLA-adapted epitopes elicited equivalent or higher-magnitude IFN-γ responses than did the nonadapted epitopes, particularly in Nef. At a population level, inclusion of all of the immunoreactive variant CD8 T cell epitopes in Gag, Pol, Nef, and Env suggested that HIV adaptation leads to an inflation of Nef-directed immune responses relative to other proteins. We concluded that HLA-HIV associations mark viral epitopes subject to CD8 T cell selection. These results can be used to guide functional studies of specific epitopes and escape mutations, as well as to test, train, and evaluate analytical models of viral escape and fitness. The inflation of Nef and HLA-adapted variant responses may have negative effects on natural and vaccine immunity against HIV and, therefore, has implications for diversity coverage approaches in HIV vaccine design.

Early selection in Gag by protective HLA alleles contributes to reduced HIV-1 replication capacity that may be largely compensated for in chronic infection

Mutations that allow escape from CD8 T-cell responses are common in HIV-1 and may attenuate pathogenesis by reducing viral fitness. While this has been demonstrated for individual cases, a systematic investigation of the consequence of HLA class I-mediated selection on HIV-1 in vitro replication capacity (RC) has not been undertaken. We examined this question by generating recombinant viruses expressing plasma HIV-1 RNA-derived Gag-Protease sequences from 66 acute/early and 803 chronic untreated subtype B-infected individuals in an NL4-3 background and measuring their RCs using a green fluorescent protein (GFP) reporter CD4 T-cell assay. In acute/early infection, viruses derived from individuals expressing the protective alleles HLA-B*57, -B*5801, and/or -B*13 displayed significantly lower RCs than did viruses from individuals lacking these alleles (P < 0.05). Furthermore, acute/early RC inversely correlated with the presence of HLA-B-associated Gag polymorphisms (R = -0.27; P = 0.03), suggesting a cumulative effect of primary escape mutations on fitness during the first months of infection. At the chronic stage of infection, no strong correlations were observed between RC and protective HLA-B alleles or with the presence of HLA-B-associated polymorphisms restricted by protective alleles despite increased statistical power to detect these associations. However, RC correlated positively with the presence of known compensatory mutations in chronic viruses from B*57-expressing individuals harboring the Gag T242N mutation (n = 50; R = 0.36; P = 0.01), suggesting that the rescue of fitness defects occurred through mutations at secondary sites. Additional mutations in Gag that may modulate the impact of the T242N mutation on RC were identified. A modest inverse correlation was observed between RC and CD4 cell count in chronic infection (R = -0.17; P < 0.0001), suggesting that Gag-Protease RC could increase over the disease course. Notably, this association was stronger for individuals who expressed B*57, B*58, or B*13 (R = -0.27; P = 0.004). Taken together, these data indicate that certain protective HLA alleles contribute to early defects in HIV-1 fitness through the selection of detrimental mutations in Gag; however, these effects wane as compensatory mutations accumulate in chronic infection. The long-term control of HIV-1 in some persons who express protective alleles suggests that early fitness hits may provide lasting benefits.

Adaptive interactions between HLA and HIV-1: highly divergent selection imposed by HLA class I molecules with common supertype motifs

Currently, 1.1 million individuals in the United States are living with HIV-1 infection. Although this is a relatively small proportion of the global pandemic, the remarkable mix of ancestries in the United States, drawn together over the past two centuries of continuous population migrations, provides an important and unique perspective on adaptive interactions between HIV-1 and human genetic diversity. HIV-1 is a rapidly adaptable organism and mutates within or near immune epitopes that are determined by the HLA class I genotype of the infected host. We characterized HLA-associated polymorphisms across the full HIV-1 proteome in a large, ethnically diverse national United States cohort of HIV-1-infected individuals. We found a striking divergence in the immunoselection patterns associated with HLA variants that have very similar or identical peptide-binding specificities but are differentially distributed among racial/ethnic groups. Although their similarity in peptide binding functionally clusters these HLA variants into supertypes, their differences at sites within the peptide-binding groove contribute to race-specific selection effects on circulating HIV-1 viruses. This suggests that the interactions between the HLA/HIV peptide complex and the TCR vary significantly within HLA supertype groups, which, in turn, influences HIV-1 evolution.