HLA class I-driven evolution of human immunodeficiency virus type 1 subtype c proteome: immune escape and viral load

Insights to HIV-1 coreceptor usage by estimating HLA adaptation with Bayesian generalized linear mixed models

The mechanisms triggering the human immunodeficiency virus type I (HIV-1) to switch the coreceptor usage from CCR5 to CXCR4 during the course of infection are not entirely understood. While low CD4+ T cell counts are associated with CXCR4 usage, a predominance of CXCR4 usage with still high CD4+ T cell counts remains puzzling. Here, we explore the hypothesis that viral adaptation to the human leukocyte antigen (HLA) complex, especially to the HLA class II alleles, contributes to the coreceptor switch. To this end, we sequence the viral gag and env protein with corresponding HLA class I and II alleles of a new cohort of 312 treatment-naive, subtype C, chronically-infected HIV-1 patients from South Africa. To estimate HLA adaptation, we develop a novel computational approach using Bayesian generalized linear mixed models (GLMMs). Our model allows to consider the entire HLA repertoire without restricting the model to pre-learned HLA-polymorphisms. In addition, we correct for phylogenetic relatedness of the viruses within the model itself to account for founder effects. Using our model, we observe that CXCR4-using variants are more adapted than CCR5-using variants (p-value = 1.34e-2). Additionally, adapted CCR5-using variants have a significantly lower predicted false positive rate (FPR) by the geno2pheno[coreceptor] tool compared to the non-adapted CCR5-using variants (p-value = 2.21e-2), where a low FPR is associated with CXCR4 usage. Consequently, estimating HLA adaptation can be an asset in predicting not only coreceptor usage, but also an approaching coreceptor switch in CCR5-using variants. We propose the usage of Bayesian GLMMs for modeling virus-host adaptation in general.

HIV-1 Group M Capsid Amino Acid Variability: Implications for Sequence Quality Control of Genotypic Resistance Testing

BACKGROUND

With the approval of the HIV-1 capsid inhibitor, lenacapavir, capsid sequencing will be required for managing lenacapavir-experienced individuals with detectable viremia. Successful sequence interpretation will require examining new capsid sequences in the context of previously published sequence data.

METHODS

We analyzed published HIV-1 group M capsid sequences from 21,012 capsid-inhibitor naïve individuals to characterize amino acid variability at each position and influence of subtype and cytotoxic T lymphocyte (CTL) selection pressure. We determined the distributions of usual mutations, defined as amino acid differences from the group M consensus, with a prevalence ≥ 0.1%. Co-evolving mutations were identified using a phylogenetically-informed Bayesian graphical model method.

RESULTS

162 (70.1%) positions had no usual mutations (45.9%) or only conservative usual mutations with a positive BLOSUM62 score (24.2%). Variability correlated independently with subtype-specific amino acid occurrence (Spearman rho = 0.83; p < 1 × 10-9) and the number of times positions were reported to contain an HLA-associated polymorphism, an indicator of CTL pressure (rho = 0.43; p = 0.0002).

CONCLUSIONS

Knowing the distribution of usual capsid mutations is essential for sequence quality control. Comparing capsid sequences from lenacapavir-treated and lenacapavir-naïve individuals will enable the identification of additional mutations potentially associated with lenacapavir therapy.

Highly Ambiguous HIV-1 Pol Positions Encoding Multiple Amino Acids Usually Result from Antiviral or Immune Selection Pressure

HIV-1 pol nucleotide ambiguities encoding amino acid mixtures occur commonly during population-based genotypic drug resistance testing. However, few studies have addressed the validity of sequences with fully ambiguous codons (FACs) containing codons translatable to more than four amino acids. We identified 839 published HIV-1 pol sequences with 846 FACs at 131 positions and determined their distribution relative to 215 HLA-associated pol positions (HAPs) and 84 drug-resistance positions. Among HIV-1 reverse transcriptase (RT) and protease sequences from antiretroviral therapy (ART)-naive and -experienced persons, there was a strong correlation between the likelihood a position was a FAC and that it was an HAP (Spearman's correlation coefficient rho >0.40; p < 1e-6). Among HIV-1 RT sequences from ART-experienced persons, there was a correlation between the likelihood that a position was a FAC and that it was a drug-resistance position (rho = 0.2; p = 8e-4). In the context of population-based genotypic resistance testing, FACs usually result from antiviral or immune selection pressure.

HAMdetector: A Bayesian regression model that integrates information to detect HLA-associated mutations

MOTIVATION

A key process in anti-viral adaptive immunity is that the Human Leukocyte Antigen system (HLA) presents epitopes as Major Histocompatibility Complex I (MHC I) protein-peptide complexes on cell surfaces and in this way alerts CD8+ cytotoxic T-Lymphocytes (CTLs). This pathway exerts strong selection pressure on viruses, favoring viral mutants that escape recognition by the HLA/CTL system. Naturally, such immune escape mutations often emerge in highly variable viruses, e.g. HIV or HBV, as HLA-associated mutations (HAMs), specific to the hosts MHC I proteins. The reliable identification of HAMs is not only important for understanding viral genomes and their evolution, but it also impacts the development of broadly effective anti-viral treatments and vaccines against variable viruses. By their very nature, HAMs are amenable to detection by statistical methods in paired sequence/HLA data. However, HLA alleles are very polymorphic in the human host population which makes the available data relatively sparse and noisy. Under these circumstances, one way to optimize HAM detection is to integrate all relevant information in a coherent model. Bayesian inference offers a principled approach to achieve this.

RESULTS

We present a new Bayesian regression model for the detection of HAMs that integrates a sparsity-inducing prior, epitope predictions, and phylogenetic bias assessment, and that yields easily interpretable quantitative information on HAM candidates. The model predicts experimentally confirmed HAMs as having high posterior probabilities, and it performs well in comparison to state-of-the-art models for several data sets from individuals infected with HBV, HDV, and HIV.

AVAILABILITY

The source code of this software is available at https://github.com/HAMdetector/Escape.jl under a permissive MIT license.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Predicted impact of the viral mutational landscape on the cytotoxic response against SARS-CoV-2

The massive assessment of immune evasion due to viral mutations that increase COVID-19 susceptibility can be computationally facilitated. The adaptive cytotoxic T response is critical during primary infection and the generation of long-term protection. Here, potential HLA class I epitopes in the SARS-CoV-2 proteome were predicted for 2,915 human alleles of 71 families using the netMHCIpan EL algorithm. Allele families showed extreme epitopic differences, underscoring genetic variability of protective capacity between humans. Up to 1,222 epitopes were associated with any of the twelve supertypes, that is, allele clusters covering 90% population. Next, from all mutations identified in ~118,000 viral NCBI isolates, those causing significant epitope score reduction were considered epitope escape mutations. These mutations mainly involved non-conservative substitutions at the second and C-terminal position of the ligand core, or total ligand removal by large recurrent deletions. Escape mutations affected 47% of supertype epitopes, which in 21% of cases concerned isolates from two or more sub-continental areas. Some of these changes were coupled, but never surpassed 15% of evaded epitopes for the same supertype in the same isolate, except for B27. In contrast to most supertypes, eight allele families mostly contained alleles with few SARS-CoV-2 ligands. Isolates harboring cytotoxic escape mutations for these families co-existed geographically within sub-Saharan and Asian populations enriched in these alleles according to the Allele Frequency Net Database. Collectively, our findings indicate that escape mutation events have already occurred for half of HLA class I supertype epitopes. However, it is presently unlikely that, overall, it poses a threat to the global population. In contrast, single and double mutations for susceptible alleles may be associated with viral selective pressure and alarming local outbreaks. The integration of genomic, geographical and immunoinformatic information eases the surveillance of variants potentially affecting the global population, as well as minority subpopulations.

HLA-associated preadaptation in HIV Vif is associated with higher set point viral load and faster CD4+ decline in Zambian transmission pairs

OBJECTIVE S

We investigated the relationship between human leukocyte antigen (HLA)-associated preadaptation for the entire subtype C HIV-1 proteome of the transmitted founder virus and subsequent HIV-1 disease progression in a cohort of heterosexual linked transmission pairs in Zambia.

DESIGN

An adaptation model was used to calculate an adaptation score for each virus-HLA combination in order to quantify the degree of preadaptation of the transmitted virus to the linked recipient's HLA alleles. These scores were then assessed for their relationship to viral load and longitudinal CD4+ decline in the recipient.

METHODS

Viral RNA was extracted from the plasma of the donor partner and the linked recipient near the time of transmission, as well as longitudinally from the linked recipient. Viral adaptation scores were calculated for each individual and each protein in the subtype C HIV-1 proteome.

RESULTS

The majority of HLA-associated sites were located in Gag, Pol and Nef; however, proportional to protein length, the accessory and regulatory proteins contained a relatively high proportion of HLA-associated sites. Over the course of infection, HLA-mediated immune adaptation increased for all proteins except Vpu and gp120. Preadaptation was positively associated with higher early set point viral load and faster CD4+ decline. When examined by protein, preadaptation in Pol and Vif were statistically significantly associated with these markers of disease progression.

CONCLUSION

Adaptation in Pol had the greatest impact on viral control. Despite containing a large proportion of HLA-associated sites, Vif was the only regulatory or accessory protein for which preadaptation significantly correlated with disease progression.

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.

Mucosal IL-4R antagonist HIV vaccination with SOSIP-gp140 booster can induce high-quality cytotoxic CD4+/CD8+ T cells and humoral responses in macaques

Inducing humoral, cellular and mucosal immunity is likely to improve the effectiveness of HIV-1 vaccine strategies. Here, we tested a vaccine regimen in pigtail macaques using an intranasal (i.n.) recombinant Fowl Pox Virus (FPV)-gag pol env-IL-4R antagonist prime, intramuscular (i.m.) recombinant Modified Vaccinia Ankara Virus (MVA)-gag pol-IL-4R antagonist boost followed by an i.m SOSIP-gp140 boost. The viral vector-expressed IL-4R antagonist transiently inhibited IL-4/IL-13 signalling at the vaccination site. The SOSIP booster not only induced gp140-specific IgG, ADCC (antibody-dependent cellular cytotoxicity) and some neutralisation activity, but also bolstered the HIV-specific cellular and humoral responses. Specifically, superior sustained systemic and mucosal HIV Gag-specific poly-functional/cytotoxic CD4⁺ and CD8⁺ T cells were detected with the IL-4R antagonist adjuvanted strategy compared to the unadjuvanted control. In the systemic compartment elevated Granzyme K expression was linked to CD4⁺ T cells, whilst Granzyme B/TIA-1 to CD8⁺ T cells. In contrast, the cytotoxic marker expression by mucosal CD4⁺ and CD8⁺ T cells differed according to the mucosal compartment. This vector-based mucosal IL-4R antagonist/SOSIP booster strategy, which promotes cytotoxic mucosal CD4⁺ T cells at the first line of defence, and cytotoxic CD4⁺ and CD8⁺ T cells plus functional antibodies in the blood, may prove valuable in combating mucosal infection with HIV-1 and warrants further investigation.

Common evolutionary features of the envelope glycoprotein of HIV-1 in patients belonging to a transmission chain

The diversity of the HIV-1 envelope glycoproteins (Env) is largely a consequence of the pressure exerted by the adaptive immune response to infection. While it was generally assumed that the neutralizing antibody (NAb) response depended mainly on the infected individual, the concept that virus-related factors could be important in inducing this response has recently emerged. Here, we analyzed the influence of the infecting viral strain in shaping NAb responses in four HIV-1 infected subjects belonging to a transmission chain. We also explored the impact of NAb responses on the functional evolution of the viral quasispecies. The four patients developed a strong autologous neutralizing antibody response that drove viral escape and coincided with a parallel evolution of their infecting quasispecies towards increasing infectious properties, increasing susceptibility to T20 and increasing resistance to both CD4 analogs and V3 loop-directed NAbs. This evolution was associated with identical Env sequence changes at several positions in the V3 loop, the fusion peptide and the HR2 domain of gp41. The common evolutionary pattern of Env in different hosts suggests that the capacity of a given Env to adapt to changing environments may be restricted by functional constraints that limit its evolutionary landscape.

Immune escape mutations in HIV-1 controllers in the Brazilian Amazon region

BACKGROUND

Human immunodeficiency virus (HIV-1) infection is characterized by high viral replication and a decrease in CD4+ T cells (CD4+TC), resulting in AIDS, which can lead to death. In elite controllers and viremia controllers, viral replication is naturally controlled, with maintenance of CD4+TC levels without the use of antiretroviral therapy (ART).

METHODS

The aim of the present study was to describe virological and immunological risk factors among HIV-1-infected individuals according to characteristics of progression to AIDS. The sample included 30 treatment-naive patients classified into three groups based on infection duration (> 6 years), CD4+TC count and viral load: (i) 2 elite controllers (ECs), (ii) 7 viremia controllers (VCs) and (iii) 21 nonviremia controllers (NVCs). Nested PCR was employed to amplify the virus genome, which was later sequenced using the Ion PGM platform for subtyping and analysis of immune escape mutations.

RESULTS

Viral samples were classified as HIV-1 subtypes B and F. Greater selection pressure on mutations was observed in the group of viremia controllers, with a higher frequency of immunological escape mutations in the genes investigated, including two new mutations in gag. The viral sequences of viremia controllers and nonviremia controllers did not differ significantly regarding the presence of immune escape mutations.

CONCLUSION

The results suggest that progression to AIDS is not dependent on a single variable but rather on a set of characteristics and pressures exerted by virus biology and interactions with immunogenetic host factors.

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.

Pol-Driven Replicative Capacity Impacts Disease Progression in HIV-1 Subtype C Infection

CD8⁺ T cell-mediated escape mutations in Gag can reduce HIV-1 replication capacity (RC) and alter disease progression, but less is known about immune-mediated attenuation in other HIV-1 proteins. We generated 487 recombinant viruses encoding RT-integrase from individuals with chronic (n = 406) and recent (n = 81) HIV-1 subtype C infection and measured their in vitro RC using a green fluorescent protein (GFP) reporter T cell assay. In recently infected individuals, reverse transcriptase (RT)-integrase-driven RC correlated significantly with viral load set point (r = 0.25; P = 0.03) and CD4⁺ T cell decline (P = 0.013). Moreover, significant associations between RT integrase-driven RC and viral load (r = 0.28; P < 0.0001) and CD4⁺ T cell count (r = -0.29; P < 0.0001) remained in chronic infection. In early HIV infection, host expression of the protective HLA-B*81 allele was associated with lower RC (P = 0.05), as was expression of HLA-B*07 (P = 0.02), suggesting early immune-driven attenuation of RT-integrase by these alleles. In chronic infection, HLA-A*30:09 (in linkage disequilibrium with HLA-B*81) was significantly associated with lower RC (P = 0.05), and all 6 HLA-B alleles with the lowest RC measurements represented protective alleles, consistent with long-term effects of host immune pressures on lowering RT-integrase RC. The polymorphisms V241I, I257V, P272K, and E297K in reverse transcriptase and I201V in integrase, all relatively uncommon polymorphisms occurring in or adjacent to optimally described HLA-restricted cytotoxic T-lymphocyte epitopes, were associated with reduced RC. Together, our data suggest that RT-integrase-driven RC is clinically relevant and provide evidence that immune-driven selection of mutations in RT-integrase can compromise RC.IMPORTANCE Identification of viral mutations that compromise HIV's ability to replicate may aid rational vaccine design. However, while certain escape mutations in Gag have been shown to reduce HIV replication and influence clinical progression, less is known about the consequences of mutations that naturally arise in other HIV proteins. Pol is a highly conserved protein, but the impact of Pol function on HIV disease progression is not well defined. Here, we generated recombinant viruses using the RT-integrase region of Pol derived from HIV-1C-infected individuals with recent and chronic infection and measured their ability to replicate in vitro We demonstrate that RT-integrase-driven replication ability significantly impacts HIV disease progression. We further show evidence of immune-mediated attenuation in RT-integrase and identify specific polymorphisms in RT-integrase that significantly decrease HIV-1 replication ability, suggesting which Pol epitopes could be explored in vaccine development.

Extensive host immune adaptation in a concentrated North American HIV epidemic

Supplemental Digital Content is available in the text

Objective:

HIV incidence in the Canadian province of Saskatchewan, where Indigenous persons make up 80% of those infected, are among the highest on the continent. Reports of accelerated HIV progression, associated with carriage of certain human leukocyte antigen (HLA) alleles (including the typically protective HLA-B∗51) have also emerged from the region. Given that acquisition of HIV preadapted to host HLA negatively impacts clinical outcome, we hypothesized that HIV-host adaptation may be elevated in Saskatchewan.

Design:

Comparative analysis of population-level HIV sequence datasets from Saskatchewan and elsewhere in Canada/USA.

Methods:

We analyzed 1144 HIV subtype B Pol sequences collected in Saskatchewan between 2000 and 2016, comprising ∼65% of cumulative provincial HIV cases, for the presence of 70 HLA-associated Pol mutations. Sequences from British Columbia (N = 6525) and elsewhere in Canada/USA (N = 6517) were used for comparison. HIV adaptation levels to 34 HLA alleles were also computed. Putative HIV transmission clusters were identified, and the prevalence of HLA-associated adaptations within and outside these clusters was investigated.

Results:

Analyses confirmed significantly elevated and temporally increasing levels of HIV adaptation to commonly expressed HLA alleles, in particular B∗51. Notably, HLA-adapted HIV strains were significantly enriched among phylogenetic clusters in Saskatchewan.

Conclusion:

Extensive circulating HIV adaptation to HLA in Saskatchewan provides a plausible explanation for accelerated progression, while enrichment of adapted variants in phylogenetic clusters suggests they are being widely transmitted. Results highlight the utility of Pol sequences, routinely collected for drug resistance monitoring, for surveillance of HIV-host adaptation, and underscore the urgent need to expand HIV prevention and treatment programmes in Saskatchewan.

Detecting consistent patterns of directional adaptation using differential selection codon models

BACKGROUND

Phylogenetic codon models are often used to characterize the selective regimes acting on protein-coding sequences. Recent methodological developments have led to models explicitly accounting for the interplay between mutation and selection, by modeling the amino acid fitness landscape along the sequence. However, thus far, most of these models have assumed that the fitness landscape is constant over time. Fluctuations of the fitness landscape may often be random or depend on complex and unknown factors. However, some organisms may be subject to systematic changes in selective pressure, resulting in reproducible molecular adaptations across independent lineages subject to similar conditions.

RESULTS

Here, we introduce a codon-based differential selection model, which aims to detect and quantify the fine-grained consistent patterns of adaptation at the protein-coding level, as a function of external conditions experienced by the organism under investigation. The model parameterizes the global mutational pressure, as well as the site- and condition-specific amino acid selective preferences. This phylogenetic model is implemented in a Bayesian MCMC framework. After validation with simulations, we applied our method to a dataset of HIV sequences from patients with known HLA genetic background. Our differential selection model detects and characterizes differentially selected coding positions specifically associated with two different HLA alleles.

CONCLUSION

Our differential selection model is able to identify consistent molecular adaptations as a function of repeated changes in the environment of the organism. These models can be applied to many other problems, ranging from viral adaptation to evolution of life-history strategies in plants or animals.

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.

Early evolution of human leucocyte antigen-associated escape mutations in variable Gag proteins predicts CD4+ decline in HIV-1 subtype C-infected women

OBJECTIVE

HIV-1 escape from cytotoxic T-lymphocytes results in the accumulation of human leucocyte antigen (HLA)-associated mutations in the viral genome. To understand the contribution of early escape to disease progression, this study investigated the evolution and pathogenic implications of cytotoxic T-lymphocyte escape in a cohort followed from infection for 5 years.

METHODS

Viral loads and CD4 cell counts were monitored in 78 subtype C-infected individuals from onset of infection until CD4 cell count decline to less than 350 cells/μl or 5 years postinfection. The gag gene was sequenced and HLA-associated changes between enrolment and 12 months postinfection were mapped.

RESULTS

HLA-associated escape mutations were identified in 48 (62%) of the participants and were associated with CD4 decline to less than 350 cells/μl (P = 0.05). Escape mutations in variable Gag proteins (p17 and p7p6) had a greater impact on disease progression than escape in more conserved regions (p24) (P = 0.03). The association between HLA-associated escape mutations and CD4 decline was independent of protective HLA allele (B57, B58 : 01 and B81) expression.

CONCLUSION

The high frequency of escape contributed to rapid disease progression in this cohort. Although HLA-adaption in both conserved and variable Gag domains in the first year of infection was detrimental to long-term clinical outcome, escape in variable domains had greater impact.

HIV and host immunogenetics: unraveling the role of HLA-C

Host genetic factors play a major role in determining the outcome of many infections including human immunodeficiency virus (HIV). Multiple host factors have been studied till date showing their varied role in susceptibility or resistance to HIV infection. HLA-C, however, has been recently started gaining interest in researchers mind revealing its polymorphisms to have an important effect on viral load set-points, disease progression as well as transmission. In this review report, we have compiled these significant findings of HLA-C in HIV infection, in an attempt to highlight the need for further research in the area in different ethnic population to establish its role in the infection.

Transmitted/Founder HIV-1 Subtype C Viruses Show Distinctive Signature Patterns in Vif, Vpr, and Vpu That Are Under Subsequent Immune Pressure During Early Infection

Viral variants that predominate during early infection may exhibit constrained diversity compared with those found during chronic infection and could contain amino acid signature patterns that may enhance transmission, establish productive infection, and influence early events that modulate the infection course. We compared amino acid distributions in 17 patients recently infected with HIV-1C with patients with chronic infection. We found significantly lower entropy in inferred transmitted/founder (t/f) compared with chronic viruses and identified signature patterns in Vif and Vpr from inferred t/f viruses. We investigated sequence evolution longitudinally up to 500 days postseroconversion and compared the impact of selected substitutions on predicted human leukocyte antigen (HLA) binding affinities of published and predicted cytotoxic T-lymphocyte epitopes. Polymorphisms in Vif and Vpr during early infection occurred more frequently at epitope-HLA anchor residues and significantly decreased predicted epitope-HLA binding. Transmission-associated sequence signatures may have implications for novel strategies to prevent HIV-1 transmission.

Functional characterization of Vif proteins from HIV-1 infected patients with different APOBEC3G haplotypes

OBJECTIVE

The human cytidine deaminase APOBEC3G (A3G) potently restricts HIV-1 but the virus, in turn, expresses a Vif protein which degrades A3G. A natural A3G-H186R variant, common in African populations, has been associated with a more rapid AIDS disease progression, but the underlying mechanism remains unknown. We hypothesized that differences in HIV-1 Vif activity towards A3G wild type and A3G-H186R contribute to the distinct clinical AIDS manifestation.

METHODS

Vif variants were cloned from plasma samples of 26 South African HIV-1 subtype C infected patients, which either express wild type A3G or A3G-H186R. The Vif alleles were assessed for their ability to counteract A3G variants using western blot and single-cycle infectivity assays.

RESULTS

We obtained a total of 392 Vif sequences which displayed an amino acid sequence difference of 6.2-19.2% between patients. The intrapatient Vif diversities from patient groups A3G, A3G and A3G were similar. Vif variants obtained from patients expressing A3G and A3G were capable of counteracting both A3G variants with similar efficiency. However, the antiviral activity of A3G-H186R was significantly reduced in both the presence and absence of Vif, indicating that the A3G-H186R variant intrinsically exerts less antiviral activity.

CONCLUSION

A3G wild type and A3G-H186R are equally susceptible to counteraction by Vif, regardless of whether the Vif variant was obtained from A3G and A3G patients. However, the A3G-H186R variant intrinsically displayed lower antiviral activity, which could explain the higher plasma viral loads and accelerated disease progression reported for patients expressing A3G.

Bayesian Inference of the Evolution of a Phenotype Distribution on a Phylogenetic Tree

The distribution of a phenotype on a phylogenetic tree is often a quantity of interest. Many phenotypes have imperfect heritability, so that a measurement of the phenotype for an individual can be thought of as a single realization from the phenotype distribution of that individual. If all individuals in a phylogeny had the same phenotype distribution, measured phenotypes would be randomly distributed on the tree leaves. This is, however, often not the case, implying that the phenotype distribution evolves over time. Here we propose a new model based on this principle of evolving phenotype distribution on the branches of a phylogeny, which is different from ancestral state reconstruction where the phenotype itself is assumed to evolve. We develop an efficient Bayesian inference method to estimate the parameters of our model and to test the evidence for changes in the phenotype distribution. We use multiple simulated data sets to show that our algorithm has good sensitivity and specificity properties. Since our method identifies branches on the tree on which the phenotype distribution has changed, it is able to break down a tree into components for which this distribution is unique and constant. We present two applications of our method, one investigating the association between HIV genetic variation and human leukocyte antigen and the other studying host range distribution in a lineage of Salmonella enterica, and we discuss many other potential applications.

Effects of a Single Escape Mutation on T Cell and HIV-1 Co-adaptation

The mechanistic basis for the progressive accumulation of Y(135)F Nef mutant viruses in the HIV-1-infected population remains poorly understood. Y(135)F viruses carry the 2F mutation within RW8 and RF10, which are two HLA-A(∗)24:02-restricted superimposed Nef epitopes recognized by distinct and adaptable CD8(+) T cell responses. We combined comprehensive analysis of the T cell receptor repertoire and cross-reactive potential of wild-type or 2F RW8- and RF10-specific CD8(+) T cells with peptide-MHC complex stability and crystal structure studies. We find that, by affecting direct and water-mediated hydrogen bond networks within the peptide-MHC complex, the 2F mutation reduces both TCR and HLA binding. This suggests an advantage underlying the evolution of the 2F variant with decreased CD8(+) T cell efficacy. Our study provides a refined understanding of HIV-1 and CD8(+) T cell co-adaptation at the population level.

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.

Population-Level Immune-Mediated Adaptation in HIV-1 Polymerase during the North American Epidemic

Human leukocyte antigen (HLA) class I-associated polymorphisms in HIV-1 that persist upon transmission to HLA-mismatched hosts may spread in the population as the epidemic progresses. Transmission of HIV-1 sequences containing such adaptations may undermine cellular immune responses to the incoming virus in future hosts. Building upon previous work, we investigated the extent of HLA-associated polymorphism accumulation in HIV-1 polymerase (Pol) through comparative analysis of linked HIV-1/HLA class I genotypes sampled during historic (1979 to 1989; n = 338) and modern (2001 to 2011; n = 278) eras from across North America (Vancouver, BC, Canada; Boston, MA; New York, NY; and San Francisco, CA). Phylogenies inferred from historic and modern HIV-1 Pol sequences were star-like in shape, with an inferred most recent common ancestor (epidemic founder virus) sequence nearly identical to the modern North American subtype B consensus sequence. Nevertheless, modern HIV-1 Pol sequences exhibited roughly 2-fold-higher patristic (tip-to-tip) genetic distances than historic sequences, with HLA pressures likely driving ongoing diversification. Moreover, the frequencies of published HLA-associated polymorphisms in individuals lacking the selecting HLA class I allele was on average ∼2.5-fold higher in the modern than in the historic era, supporting their spread in circulation, though some remained stable in frequency during this time. Notably, polymorphisms restricted by protective HLA alleles appear to be spreading to a greater relative extent than others, though these increases are generally of modest absolute magnitude. However, despite evidence of polymorphism spread, North American hosts generally remain at relatively low risk of acquiring an HIV-1 polymerase sequence substantially preadapted to their HLA profiles, even in the present era.

IMPORTANCE

HLA class I-restricted cytotoxic T-lymphocyte (CTL) escape mutations in HIV-1 that persist upon transmission may accumulate in circulation over time, potentially undermining host antiviral immunity to the transmitted viral strain. We studied >600 experimentally collected HIV-1 polymerase sequences linked to host HLA information dating back to 1979, along with phylogenetically reconstructed HIV-1 sequences dating back to the virus' introduction into North America. Overall, our results support the gradual spread of many-though not all-HIV-1 polymerase immune escape mutations in circulation over time. This is consistent with recent observations from other global regions, though the extent of polymorphism accumulation in North America appears to be lower than in populations with high seroprevalence, older epidemics, and/or limited HLA diversity. Importantly, the risk of acquiring an HIV-1 polymerase sequence at transmission that is substantially preadapted to one's HLA profile remains relatively low in North America, even in the present era.

High frequency of transmitted HIV-1 Gag HLA class I-driven immune escape variants but minimal immune selection over the first year of clade C infection

In chronic HIV infection, CD8+ T cell responses to Gag are associated with lower viral loads, but longitudinal studies of HLA-restricted CD8+ T cell-driven selection pressure in Gag from the time of acute infection are limited. In this study we examined Gag sequence evolution over the first year of infection in 22 patients identified prior to seroconversion. A total of 310 and 337 full-length Gag sequences from the earliest available samples (median = 14 days after infection [Fiebig stage I/II]) and at one-year post infection respectively were generated. Six of 22 (27%) individuals were infected with multiple variants. There was a trend towards early intra-patient viral sequence diversity correlating with viral load set point (p = 0.07, r = 0.39). At 14 days post infection, 59.7% of Gag CTL epitopes contained non-consensus polymorphisms and over half of these (35.3%) comprised of previously described CTL escape variants. Consensus and variant CTL epitope proportions were equally distributed irrespective of the selecting host HLA allele and most epitopes remained unchanged over 12 months post infection. These data suggest that intrapatient diversity during acute infection is an indicator of disease outcome. In this setting, there is a high rate of transmitted CTL escape variants and limited immune selection in Gag during the first year of infection. These data have relevance for vaccine strategies designed to elicit effective CD8+ T cell immune responses.

Possible footprints of APOBEC3F and/or other APOBEC3 deaminases, but not APOBEC3G, on HIV-1 from patients with acute/early and chronic infections

Members of the apolipoprotein B mRNA-editing enzyme-catalytic polypeptide-like-3 (APOBEC3) innate cellular cytidine deaminase family, particularly APOBEC3F and APOBEC3G, can cause extensive and lethal G-to-A mutations in HIV-1 plus-strand DNA (termed hypermutation). It is unclear if APOBEC3-induced mutations in vivo are always lethal or can occur at sublethal levels that increase HIV-1 diversification and viral adaptation to the host. The viral accessory protein Vif counteracts APOBEC3 activity by binding to APOBEC3 and promoting proteasome degradation; however, the efficiency of this interaction varies, since a range of hypermutation frequencies are observed in HIV-1 patient DNA. Therefore, we examined "footprints" of APOBEC3G and APOBEC3F activity in longitudinal HIV-1 RNA pol sequences from approximately 3,000 chronically infected patients by determining whether G-to-A mutations occurred in motifs that were favored or disfavored by these deaminases. G-to-A mutations were more frequent in APOBEC3G-disfavored than in APOBEC3G-favored contexts. In contrast, mutations in APOBEC3F-disfavored contexts were relatively rare, whereas mutations in contexts favoring APOBEC3F (and possibly other deaminases) occurred 16% more often than average G-to-A mutations. These results were supported by analyses of >500 HIV-1 env sequences from acute/early infection.

IMPORTANCE

Collectively, our results suggest that APOBEC3G-induced mutagenesis is lethal to HIV-1, whereas mutagenesis caused by APOBEC3F and/or other deaminases may result in sublethal mutations that might facilitate viral diversification. Therefore, Vif-specific cytotoxic T lymphocyte (CTL) responses and drugs that manipulate the interplay between Vif and APOBEC3 may have beneficial or detrimental clinical effects depending on how they affect the binding of Vif to various members of the APOBEC3 family.

The impact of the interferon-lambda family on the innate and adaptive immune response to viral infections

Type-III interferons (IFN-λ, IFNL) are the most recently described family of IFNs. This family of innate cytokines are increasingly being ascribed pivotal roles in host-pathogen interactions. Herein, we will review the accumulating evidence detailing the immune biology of IFNL during viral infection, and the implications of this novel information on means to advance the development of therapies and vaccines against existing and emerging pathogens. IFNLs exert antiviral effects via induction of IFN-stimulated genes. Common single nucleotide polymorphisms (SNPs) in the IFNL3, IFNL4 and the IFNL receptor α-subunit genes have been strongly associated with IFN-α-based treatment of chronic hepatitis C virus infection. The clinical impact of these SNPs may be dependent on the status of viral infection (acute or chronic) and the potential to develop viral resistance. Another important function of IFNLs is macrophage and dendritic cell polarization, which prime helper T-cell activation and proliferation. It has been demonstrated that IFNL increase Th1- and reduce Th2-cytokines. Therefore, can such SNPs affect the IFNL signaling and thereby modulate the Th1/Th2 balance during infection? In turn, this may influence the subsequent priming of cytotoxic T cells versus antibody-secreting B cells, with implications for the breadth and durability of the host response.

HIV subtype influences HLA-B*07:02-associated HIV disease outcome

Genetic polymorphisms within the MHC encoding region have the strongest impact on HIV disease progression of any in the human genome and provide important clues to the mechanisms of HIV immune control. Few analyses have been undertaken of HLA alleles associated with rapid disease progression. HLA-B*07:02 is an HLA class I molecule that is prevalent in most populations worldwide and that has previously been consistently linked to accelerated disease progression in B-clade infection. This study investigates the observation that HLA-B*07:02 is not associated with a high viral setpoint in C-clade infection. We examine the hypothesis that this clade-specific difference in association with disease outcome may be related to distinct targeting of CD8(+) T cell epitopes. We observed that C-clade-infected individuals with HLA-B*07:02 target a broader range of Gag epitopes, and to higher magnitudes, than do individuals infected with B-clade infection. In particular, a novel p17-Gag (Gag22-30, RPGGKKHYM) epitope is targeted in >50% of HLA-B*07:02-positive C-clade-infected individuals but clade-specific differences in this epitope result in nonimmunogenicity in B-clade infection. Only the C-clade p24-Gag "GL9" (Gag355-363, GPSHKARVL) epitope-specific CD8(+) T cell response out of 16 studied was associated with a low viral setpoint. Although this epitope was also targeted in B-clade infection, the escape mutant S357S is present at higher frequency in B-clade infection than in C-clade infection (70% versus 43% in HLA-B*07:02-negative subjects). These data support earlier studies suggesting that increased breadth of the Gag-specific CD8(+) T cell response may contribute to improved HIV immune control irrespective of the particular HLA molecules expressed.

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.

tat Exon 1 exhibits functional diversity during HIV-1 subtype C primary infection

Human immunodeficiency virus type 1 (HIV-1) Tat is a mediator of viral transcription and is involved in the control of virus replication. However, associations between HIV-1 Tat diversity and functional effects during primary HIV-1 infection are still unclear. We estimated selection pressures in tat exon 1 using the mixed-effects model of evolution with 672 viral sequences generated from 20 patients infected with HIV-1 subtype C (HIV-1C) over 500 days postseroconversion. tat exon 1 residues 3, 4, 21, 24, 29, 39, and 68 were under positive selection, and we established that specific amino acid signature patterns were apparent in primary HIV-1C infection compared with chronic infection. We assessed the impact of these mutations on long terminal repeat (LTR) activity and found that Tat activity was negatively affected by the Ala(21) substitution identified in 13/20 (65%) of patients, which reduced LTR activity by 88% (± 1%) (P < 0.001). The greatest increase in Tat activity was seen with the Gln(35)/Lys(39) double mutant that resulted in an additional 49% (± 14%) production of LTR-driven luciferase (P = 0.012). There was a moderate positive correlation between Tat-mediated LTR activity and HIV-1 RNA in plasma (P = 0.026; r = 0.400) after 180 days postseroconversion that was reduced by 500 days postseroconversion (P = 0.043; r = 0.266). Although Tat activation of the LTR is not a strong predictor of these clinical variables, there are significant linear relationships between Tat transactivation and patients' plasma viral loads and CD4 counts, highlighting the complex interplay between Tat mutations in early HIV-1C infection.

Developments in virus-like particle-based vaccines for HIV

Virus-like particles (VLPs) hold great promise for the development of effective and affordable vaccines. VLPs, indeed, are suitable for presentation and efficient delivery to antigen-presenting cells of linear as well as conformational antigens. This will ultimately result in a crosspresentation with both MHC class I and II molecules to prime CD4(+) T-helper and CD8(+) cytotoxic T cells. This review describes an update on the development and use of VLPs as vaccine approaches for HIV.

Cumulative impact of host and viral factors on HIV-1 viral-load control during early infection

In HIV-1 infection, the early set-point viral load strongly predicts both viral transmission and disease progression. The factors responsible for the wide spectrum of set-point viral loads are complex and likely reflect an interplay between the transmitted virus and genetically defined factors in both the transmitting source partner and the seroconverter. Indeed, analysis of 195 transmission pairs from Lusaka, Zambia, revealed that the viral loads in transmitting source partners contributed only ∼2% of the variance in early set-point viral loads of seroconverters (P = 0.046 by univariable analysis). In multivariable models, early set-point viral loads in seroconverting partners were a complex function of (i) the viral load in the source partner, (ii) the gender of the seroconverter, (iii) specific HLA class I alleles in the newly infected partner, and (iv) sharing of HLA-I alleles between partners in a transmission pair. Each of these factors significantly and independently contributed to the set-point viral load in the newly infected partner, accounting for up to 37% of the variance observed and suggesting that many factors operate in concert to define the early virological phenotype in HIV-1 infection.

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.

Role of transmitted Gag CTL polymorphisms in defining replicative capacity and early HIV-1 pathogenesis

Initial studies of 88 transmission pairs in the Zambia Emory HIV Research Project cohort demonstrated that the number of transmitted HLA-B associated polymorphisms in Gag, but not Nef, was negatively correlated to set point viral load (VL) in the newly infected partners. These results suggested that accumulation of CTL escape mutations in Gag might attenuate viral replication and provide a clinical benefit during early stages of infection. Using a novel approach, we have cloned gag sequences isolated from the earliest seroconversion plasma sample from the acutely infected recipient of 149 epidemiologically linked Zambian transmission pairs into a primary isolate, subtype C proviral vector, MJ4. We determined the replicative capacity (RC) of these Gag-MJ4 chimeras by infecting the GXR25 cell line and quantifying virion production in supernatants via a radiolabeled reverse transcriptase assay. We observed a statistically significant positive correlation between RC conferred by the transmitted Gag sequence and set point VL in newly infected individuals (p = 0.02). Furthermore, the RC of Gag-MJ4 chimeras also correlated with the VL of chronically infected donors near the estimated date of infection (p = 0.01), demonstrating that virus replication contributes to VL in both acute and chronic infection. These studies also allowed for the elucidation of novel sites in Gag associated with changes in RC, where rare mutations had the greatest effect on fitness. Although we observed both advantageous and deleterious rare mutations, the latter could point to vulnerable targets in the HIV-1 genome. Importantly, RC correlated significantly (p = 0.029) with the rate of CD4+ T cell decline over the first 3 years of infection in a manner that is partially independent of VL, suggesting that the replication capacity of HIV-1 during the earliest stages of infection is a determinant of pathogenesis beyond what might be expected based on set point VL alone.

In the majority of HIV-1 cases, a single virus establishes infection. However, mutations in the viral genome accumulate over time in order to avoid recognition by the host immune response. Certain mutations in the main structural protein, Gag, driven by cytotoxic T lymphocytes are detrimental to viral replication, and we showed previously that, upon transmission, viruses with higher numbers of escape mutations in Gag were associated with lower early set point viral loads. We hypothesized that this could be attributed to attenuation of the transmitted virus. Here, we have cloned the gag gene from 149 newly infected individuals from linked transmission pairs into a clade C proviral vector and determined the replicative capacity in vitro. We found that the replicative capacity conferred by the transmitted Gag correlated with set point viral loads in newly infected individuals, as well as with the viral load of the transmitting partner, and we identified previously unrecognized residues associated with increasing and decreasing replicative capacity. Importantly, we demonstrate that transmitted viruses with high replicative capacity cause more rapid CD4+ decline over the first three years, independent of viral load. This suggests that the trajectory of pathogenesis may be affected very early in infection, before adaptive immunity can respond.

HIV-1 drug resistance-associated mutations among antiretroviral-naive Thai patients with chronic HIV-1 infection

Antiretroviral therapy (ART) has increased in resource-limited settings. This study determined the prevalence of HIV-1 drug resistance-associated mutations (DRAMs) among patients with chronic HIV-1 infections and compare DRAMs between CRF01_AE and B subtypes. ART-naive Thai patients who had ART initiation between 2010 and 2011 were enrolled prospectively. Genotypic assays were performed on viral reverse transcriptase and protease genes within 4 weeks before starting ART. DRAMs were assessed using the International AIDS Society-USA 2011 list. A total of 330 patients were included. HIV-1 subtypes included CRF01_AE (73%), B (23.9%), and others (3.1%). Median (IQR) CD4+ was 66 (23-172) cells/mm(3) and median (IQR) HIV-1 RNA was 5.2 (4.6-5.8) log copies/ml. The prevalence of patients with ≥1 DRAMs for any antiretroviral agents was 17.6%. DRAM prevalence was 17% for non-nucleoside reverse transcriptase inhibitors (NNRTIs), 0.6% for NRTIs, and 0.6% for protease inhibitors (PIs). DRAMs to NNRTIs were V106I (7%), V179D (4.2%), V179T (1.8%), E138A (1.5%), V90I (1.2%), K103N (0.9%), Y181C (0.9%), and P225H (0.3%). DRAMs to NRTIs were M184V (0.3%) and T215S (0.3%). The only major DRAM for PIs was M46L (0.6%). Minor DRAMs to PIs including I13V, M36I, H69K, and L89M were observed more frequently in CRF_01 AE. By multivariate analysis, the factors "HIV-1 subtype B" and "low pretreated CD4+ cell count" were associated with a higher rate of DRAMs. HIV-1 DRAMs, especially to NNRTIs, are emerging in a middle-income country after widespread use of NNRTI-based ART. HIV genotypic assays before ART initiation in patients with chronic HIV-1 infection should be considered.

Significant reductions in Gag-protease-mediated HIV-1 replication capacity during the course of the epidemic in Japan

Human immunodeficiency virus type 1 (HIV-1) evolves rapidly in response to host immune selection pressures. As a result, the functional properties of HIV-1 isolates from earlier in the epidemic may differ from those of isolates from later stages. However, few studies have investigated alterations in viral replication capacity (RC) over the epidemic. In the present study, we compare Gag-Protease-associated RC between early and late isolates in Japan (1994 to 2009). HIV-1 subtype B sequences from 156 antiretroviral-naïve Japanese with chronic asymptomatic infection were used to construct a chimeric NL4-3 strain encoding plasma-derived gag-protease. Viral replication capacity was examined by infecting a long terminal repeat-driven green fluorescent protein-reporter T cell line. We observed a reduction in the RC of chimeric NL4-3 over the epidemic, which remained significant after adjusting for the CD4(+) T cell count and plasma virus load. The same outcome was seen when limiting the analysis to a single large cluster of related sequences, indicating that our results are not due to shifts in the molecular epidemiology of the epidemic in Japan. Moreover, the change in RC was independent of genetic distance between patient-derived sequences and wild-type NL4-3, thus ruling out potential temporal bias due to genetic similarity between patient and historic viral backbone sequences. Collectively, these data indicate that Gag-Protease-associated HIV-1 replication capacity has decreased over the epidemic in Japan. Larger studies from multiple geographical regions will be required to confirm this phenomenon.

SIV genome-wide pyrosequencing provides a comprehensive and unbiased view of variation within and outside CD8 T lymphocyte epitopes

Deep sequencing technology is revolutionizing our understanding of HIV/SIV evolution. It is known that acute SIV sequence variation within CD8 T lymphocyte (CD8-TL) epitopes is similar among MHC-identical animals, but we do not know whether this persists into the chronic phase. We now determine whether chronic viral variation in MHC-identical animals infected with clonal SIV is similar throughout the entire coding sequence when using a sensitive deep sequencing approach. We pyrosequenced the entire coding sequence of the SIV genome isolated from a unique cohort of four SIVmac239-infected, MHC-identical Mauritian cynomolgus macaques (MCM) 48 weeks after infection; one MCM in the cohort became an elite controller. Among the three non-controllers, we found that genome-wide sequences were similar between animals and we detected increased sequence complexity within 64% of CD8-TL epitopes when compared to Sanger sequencing methods. When we compared sequences between the MHC-matched controller and the three non-controllers, we found the viral population in the controller was less diverse and accumulated different variants than the viral populations in the non-controllers. Importantly, we found that initial PCR amplification of viral cDNA did not significantly affect the sequences detected, suggesting that data obtained by pyrosequencing PCR-amplified viral cDNA accurately represents the diversity of sequences replicating within an animal. This demonstrates that chronic sequence diversity across the entire SIV coding sequence is similar among MHC-identical animals with comparable viral loads when infected with the same clonal virus stock. Additionally, our approach to genome-wide SIV sequencing accurately reflects the diversity of sequences present in the replicating viral population. In sum, our study suggests that genome-wide pyrosequencing of immunodeficiency viruses captures a thorough and unbiased picture of sequence diversity, and may be a useful approach to employ when evaluating which sequences to include as part of a vaccine immunogen.

Conservation of HIV-1 T cell epitopes across time and clades: validation of immunogenic HLA-A2 epitopes selected for the GAIA HIV vaccine

HIV genomic sequence variability has complicated efforts to generate an effective globally relevant vaccine. Regions of the viral genome conserved in sequence and across time may represent the "Achilles' heel" of HIV. In this study, highly conserved T-cell epitopes were selected using immunoinformatics tools combining HLA-A2 supertype binding predictions with relative global conservation. Analysis performed in 2002 on 10,803 HIV-1 sequences, and again in 2009, on 43,822 sequences, yielded 38 HLA-A2 epitopes. These epitopes were experimentally validated for HLA binding and immunogenicity with PBMCs from HIV-infected patients in Providence, Rhode Island, and/or Bamako, Mali. Thirty-five (92%) stimulated an IFNγ response in PBMCs from at least one subject. Eleven of fourteen peptides (79%) were confirmed as HLA-A2 epitopes in both locations. Validation of these HLA-A2 epitopes conserved across time, clades, and geography supports the hypothesis that such epitopes could provide effective coverage of virus diversity and would be appropriate for inclusion in a globally relevant HIV vaccine.

HIV-1 dynamics: a reappraisal of host and viral factors, as well as methodological issues

The dynamics of HIV-1 viremia is a complex and evolving landscape with clinical and epidemiological (public health) implications. Most studies have relied on the use of set-point viral load (VL) as a readily available proxy of viral dynamics to assess host and viral correlates. This review highlights recent findings from population-based studies of set-point VL, focusing primarily on robust data related to host genetics. A comprehensive understanding of viral dynamics will clearly need to consider both host and viral characteristics, with close attention to (i) the timing of VL measurements, (ii) the biology of viral evolution, (iii) compartments of active viral replication, (iv) the transmission source partner as the immediate past microenvironment, and (v) proper application of statistical models.

Distinct evolutionary pressures underlie diversity in simian immunodeficiency virus and human immunodeficiency virus lineages

Simian immunodeficiency virus (SIV) infection of rhesus macaques causes immune depletion and disease closely resembling human AIDS and is well recognized as the most relevant animal model for the human disease. Experimental investigations of viral pathogenesis and vaccine protection primarily involve a limited set of related viruses originating in sooty mangabeys (SIVsmm). The diversity of human immunodeficiency virus type 1 (HIV-1) has evolved in humans in about a century; in contrast, SIV isolates used in the macaque model evolved in sooty mangabeys over millennia. To investigate the possible consequences of such different evolutionary histories for selection pressures and observed diversity in SIVsmm and HIV-1, we isolated, sequenced, and analyzed 20 independent isolates of SIVsmm, including representatives of 7 distinct clades of viruses isolated from natural infection. We found SIVsmm diversity to be lower overall than HIV-1 M group diversity. Reduced positive selection (i.e., less diversifying evolution) was evident in extended regions of SIVsmm proteins, most notably in Gag p27 and Env gp120. In addition, the relative diversities of proteins in the two lineages were distinct: SIVsmm Env and Gag were much less diverse than their HIV-1 counterparts. This may be explained by lower SIV-directed immune activity in mangabeys relative to HIV-1-directed immunity in humans. These findings add an additional layer of complexity to the interpretation and, potentially, to the predictive utility of the SIV/macaque model, and they highlight the unique features of human and simian lentiviral evolution that inform studies of pathogenesis and strategies for AIDS vaccine design.

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.

Differential clade-specific HLA-B*3501 association with HIV-1 disease outcome is linked to immunogenicity of a single Gag epitope

The strongest genetic influence on immune control in HIV-1 infection is the HLA class I genotype. Rapid disease progression in B-clade infection has been linked to HLA-B*35 expression, in particular to the less common HLA-B*3502 and HLA-B*3503 subtypes but also to the most prevalent subtype, HLA-B*3501. In these studies we first demonstrated that whereas HLA-B*3501 is associated with a high viral set point in two further B-clade-infected cohorts, in Japan and Mexico, this association does not hold in two large C-clade-infected African cohorts. We tested the hypothesis that clade-specific differences in HLA associations with disease outcomes may be related to distinct targeting of critical CD8(+) T-cell epitopes. We observed that only one epitope was significantly targeted differentially, namely, the Gag-specific epitope NPPIPVGDIY (NY10, Gag positions 253 to 262) (P = 2 × 10(-5)). In common with two other HLA-B*3501-restricted epitopes, in Gag and Nef, that were not targeted differentially, a response toward NY10 was associated with a significantly lower viral set point. Nonimmunogenicity of NY10 in B-clade-infected subjects derives from the Gag-D260E polymorphism present in ∼90% of B-clade sequences, which critically reduces recognition of the Gag NY10 epitope. These data suggest that in spite of any inherent HLA-linked T-cell receptor repertoire differences that may exist, maximizing the breadth of the Gag-specific CD8(+) T-cell response, by the addition of even a single epitope, may be of overriding importance in achieving immune control of HIV infection. This distinction is of direct relevance to development of vaccines designed to optimize the anti-HIV CD8(+) T-cell response in all individuals, irrespective of HLA type.

Temporal association of HLA-B*81:01- and HLA-B*39:10-mediated HIV-1 p24 sequence evolution with disease progression

HLA-B*81:01 and HLA-B*39:10 alleles have been associated with viremic control in HIV-1 subtype C infection. Both alleles restrict the TL9 epitope in p24 Gag, and cytotoxic-T-lymphocyte (CTL)-mediated escape mutations in this epitope have been associated with an in vitro fitness cost to the virus. We investigated the timing and impact of mutations in the TL9 epitope on disease progression in five B*81:01- and two B*39:10-positive subtype C-infected individuals. Whereas both B*39:10 participants sampled at 2 months postinfection had viruses with mutations in the TL9 epitope, in three of the five (3/5) B*81:01 participants, TL9 escape mutations were only detected 10 months after infection, taking an additional 10 to 15 months to reach fixation. In the two remaining B*81:01 individuals, one carried a TL9 escape variant at 2 weeks postinfection, whereas no escape mutations were detected in the virus from the other participant for up to 33 months postinfection, despite CTL targeting of the epitope. In all participants, escape mutations in TL9 were linked to coevolving residues in the region of Gag known to be associated with host tropism. Late escape in TL9, together with coevolution of putative compensatory mutations, coincided with a spontaneous increase in viral loads in two individuals who were otherwise controlling the infection. These results provide in vivo evidence of the detrimental impact of B*81:01-mediated viral evolution, in a single Gag p24 epitope, on the control of viremia.

Rapid evolution of HIV-1 to functional CD8⁺ T cell responses in humanized BLT mice

The development of mouse/human chimeras through the engraftment of human immune cells and tissues into immunodeficient mice, including the recently described humanized BLT (bone marrow, liver, thymus) mouse model, holds great promise to facilitate the in vivo study of human immune responses. However, little data exist regarding the extent to which cellular immune responses in humanized mice accurately reflect those seen in humans. We infected humanized BLT mice with HIV-1 as a model pathogen and characterized HIV-1-specific immune responses and viral evolution during the acute phase of infection. HIV-1-specific CD8(+) T cell responses in these mice were found to closely resemble those in humans in terms of their specificity, kinetics, and immunodominance. Viral sequence evolution also revealed rapid and highly reproducible escape from these responses, mirroring the adaptations to host immune pressures observed during natural HIV-1 infection. Moreover, mice expressing the protective HLA-B*57 allele exhibited enhanced control of viral replication and restricted the same CD8(+) T cell responses to conserved regions of HIV-1 Gag that are critical to its control of HIV-1 in humans. These data reveal that the humanized BLT mouse model appears to accurately recapitulate human pathogen-specific cellular immunity and the fundamental immunological mechanisms required to control a model human pathogen, aspects critical to the use of a small-animal model for human pathogens.

Further progress on defining highly conserved immunogenic epitopes for a global HIV vaccine: HLA-A3-restricted GAIA vaccine epitopes

Two major obstacles confronting HIV vaccine design have been the extensive viral diversity of HIV-1 globally and viral evolution driven by escape from CD8(+) cytotoxic T-cell lymphocyte (CTL)-mediated immune pressure. Regions of the viral genome that are not able to escape immune response and that are conserved in sequence and across time may represent the "Achilles' heel" of HIV and would be excellent candidates for vaccine development. In this study, T-cell epitopes were selected using immunoinformatics tools, combining HLA-A3 binding predictions with relative sequence conservation in the context of global HIV evolution. Twenty-seven HLA-A3 epitopes were chosen from an analysis performed in 2003 on 10,803 HIV-1 sequences, and additional sequences were selected in 2009 based on an expanded set of 43,822 sequences. These epitopes were tested in vitro for HLA binding and for immunogenicity with PBMCs of HIV-infected donors from Providence, Rhode Island. Validation of these HLA-A3 epitopes conserved across time, clades, and geography supports the hypothesis that epitopes such as these would be candidates for inclusion in our globally relevant GAIA HIV vaccine constructs.

Cytotoxic T-lymphocyte escape mutations identified by HLA association favor those which escape and revert rapidly

Identifying human immunodeficiency virus (HIV) immune escape mutations has implications for understanding the impact of host immunity on pathogen evolution and guiding the choice of vaccine antigens. One means of identifying cytotoxic-T-lymphocyte (CTL) escape mutations is to search for statistical associations between mutations and host human leukocyte antigen (HLA) class I alleles at the population level. The impact of evolutionary rates on the strength of such associations is not well defined. Here, we address this topic using a mathematical model of within-host evolution and between-host transmission of CTL escape mutants that predicts the prevalence of escape mutants at the population level. We ask how the rates at which an escape mutation emerges in a host who bears the restricting HLA and reverts when transmitted to a host who does not bear the HLA affect the strength of an association. We consider the impact of these factors when using a standard statistical method to test for an association and when using an adaptation of that method that corrects for phylogenetic relationships. We show that with both methods, the average sample size required to identify an escape mutation is smaller if the mutation escapes and reverts quickly. Thus, escape mutations identified as HLA associated systematically favor those that escape and revert rapidly. We also present expressions that can be used to infer escape and reversion rates from cross-sectional escape prevalence data.

Interplay between HIV-1 and Host Genetic Variation: A Snapshot into Its Impact on AIDS and Therapy Response

As of February 2012, 50 circulating recombinant forms (CRFs) have been reported for HIV-1 while one CRF for HIV-2. Also according to HIV sequence compendium 2011, the HIV sequence database is replete with 414,398 sequences. The fact that there are CRFs, which are an amalgamation of sequences derived from six or more subtypes (CRF27_cpx (cpx refers to complex) is a mosaic with sequences from 6 different subtypes besides an unclassified fragment), serves as a testimony to the continual divergent evolution of the virus with its approximate 1% per year rate of evolution, and this phenomena per se poses tremendous challenge for vaccine development against HIV/AIDS, a devastating disease that has killed 1.8 million patients in 2010. Here, we explore the interaction between HIV-1 and host genetic variation in the context of HIV/AIDS and antiretroviral therapy response.

Aminopeptidase substrate preference affects HIV epitope presentation and predicts immune escape patterns in HIV-infected individuals

Viruses evade immune detection partly through immune-associated mutations. Analyses of HIV sequences derived from infected individuals have identified numerous examples of HLA-associated mutations within or adjacent to T cell epitopes, but the potential impact of most mutations on epitope production and presentation remains unclear. The multistep breakdown of proteins into epitopes includes trimming of N-extended peptides into epitopes by aminopeptidases before loading onto MHC class I molecules. Definition of sequence signatures that modulate epitope production would lead to a better understanding of factors driving viral evolution and immune escape at the population level. In this study, we identified cytosolic aminopeptidases cleavage preferences in primary cells and its impact on HIV Ag degradation into epitopes in primary human cell extracts by mass spectrometry and on epitope presentation to CTL. We observed a hierarchy of preferred amino acid cleavage by cytosolic aminopeptidases. We demonstrated that flanking mutations producing more or less cleavable motifs can increase or decrease epitope production and presentation by up to 14-fold. We found that the efficiency of epitope production correlates with cleavability of flanking residues. These in vitro findings were supported by in vivo population-level analyses of clinically derived viral sequences from 1134 antiretroviral-naive HIV-infected individuals: HLA-associated mutations immune pressures drove the selection of residues that are less cleavable by aminopeptidases predominantly at N-flanking sites, leading to reduced epitope production and immune recognition. These results underscore an important and widespread role of Ag processing mutations in HIV immune escape and identify molecular mechanisms underlying impaired epitope presentation.