Tracking the culprit: HIV-1 evolution and immune selection revealed by single-genome amplification

Analysis of temporal changes in HIV-1 CRF01_AE gag genetic variability and CD8 T-cell epitope evolution

Currently, little is known about the time-dependent evolution of human immunodeficiency virus-1 (HIV-1) circulating recombinant forms (CRF) 01_AE, a dominant recombinant form associated with HIV-1 epidemics worldwide. Since gag is a highly immunodominant HIV-1 protein, we performed a comparative analysis of the CRF01_AE gag protein’s time-dependent changes and evolution. A total of 3105 HIV-1 CRF01_AE gag sequences representing 17 countries from the timeline 1990–2017 were obtained. The sequences’ phylogenetic relationship and epidemic dynamics were analyzed through a Maximum Likelihood tree and Bayesian Skyline plot, respectively. Genomic variability was measured through Shannon entropy and time-dependent immunoevolution was analyzed using changes in proteasomal degradation pattern, cytotoxic T lymphocytes (CTL) epitopes, and Human leukocyte antigens (HLA) restriction profile. The most recent common ancestor of the HIV CRF01_AE epidemic was estimated to be 1974±1. A period of exponential growth in effective population size began in 1982, fluctuated, and then stabilized in 1999. Genetic variability (entropy) consistently increased, however, epitope variability remained comparable; the highest number of novel CTL epitopes were present in 1995–1999, which were lost over time. The spread of the HIV-1 CRF01_AE epidemic is predominant in countries within Asia. Population immunogenetic pressures in the region may have contributed to the initial changes and following adaptation/stabilization of epitope diversity within gag sequences.

Increased Frequency of Inter-Subtype HIV-1 Recombinants Identified by Near Full-Length Virus Sequencing in Rwandan Acute Transmission Cohorts

Most studies of HIV-1 transmission have focused on subtypes B and C. In this study, we determined the genomic sequences of the transmitted founder (TF) viruses from acutely infected individuals enrolled between 2005 and 2011 into IAVI protocol C in Rwanda and have compared these isolates to viruses from more recent (2016–2019) acute/early infections in three at risk populations – MSM, high risk women (HRW), and discordant couples (DC). For the Protocol C samples, we utilized near full-length single genome (NFLG) amplification to generate 288 HIV-1 amplicons from 26 acutely infected seroconverters (SC), while for the 21 recent seroconverter samples (13 from HRW, two from DC, and six from MSM), we PCR amplified overlapping half-genomes. Using PacBio SMRT technology combined with the MDPseq workflow, we performed multiplex sequencing to obtain high accuracy sequences for each amplicon. Phylogenetic analyses indicated that the majority of recent transmitted viruses from DC and HRW clustered within those of the earlier Protocol C cohort. However, five of six sequences from the MSM cohort branched together and were greater than 97% identical. Recombination analyses revealed a high frequency (6/26; 23%) of unique inter-subtype recombination in Protocol C with 19% AC and 4% CD recombinant viruses, which contrasted with only 6.5% of recombinants defined by sequencing of the pol gene previously. The frequency of recombinants was significantly higher (12/21; 57%) in the more recent isolates, although, the five related viruses from the MSM cohort had identical recombination break points. While major drug resistance mutations were absent from Protocol C viruses, 4/21 of recent isolates exhibited transmitted nevirapine resistance. These results demonstrate the ongoing evolution and increased prevalence of recombinant and drug resistant transmitted viruses in Rwanda and highlight the importance of defining NFLG sequences to fully understand the nature of TF viruses and in particular the prevalence of unique recombinant forms (URFs) in transmission cohorts.

Guidelines for cloning, expression, purification and functional characterization of primary HIV-1 envelope glycoproteins

The trimeric HIV-1 envelope (Env) glycoproteins gp120 and gp41 mediate virus entry into target cells by engaging CD4 and the coreceptors CCR5 or CXCR4 at the cell surface and driving membrane fusion. Receptor/gp120 interactions regulate the virus life cycle, HIV infection transmission and pathogenesis. Env is also the target of neutralizing antibodies. Efforts have thus been made to produce soluble HIV-1 glycoproteins to develop vaccines and study the role and mechanisms of HIV/receptor interactions. However, production and purification of Env glycoproteins and their functional assessment has to cope with multiple obstacles. These include difficulties in amplifying and cloning env sequences and setting up receptor binding assays that are suitable for studies on large collections of glycoproteins, flexible enough to adapt to Env and receptor structural heterogeneities, and allow recapitulating the receptor binding properties of virion-associated Env trimers. Here we identify these difficulties and present protocols to produce primary gp120 and determination of their binding properties to receptors. The receptor binding assays confirmed that the produced glycoproteins are competent for binding CD4 and undergo proper CD4-induced conformational changes required for interaction with CCR5. These assays may help elucidate the role of gp120/receptor interactions in the pathophysiology of HIV infection and develop HIV-1 entry inhibitors.

HIV-1 subtype A gag variability and epitope evolution

Objective

The aim of this study was to examine the course of time-dependent evolution of HIV-1 subtype A on a global level, especially with respect to the dynamics of immunogenic HIV gag epitopes.

Methods

We used a total of 1,893 HIV-1 subtype A gag sequences representing a timeline from 1985 through 2010, and 19 different countries in Africa, Europe and Asia. The phylogenetic relationship of subtype A gag and its epidemic dynamics was analysed through a Maximum Likelihood tree and Bayesian Skyline plot, genomic variability was measured in terms of G→A substitutions and Shannon entropy, and the time-dependent evolution of HIV subtype A gag epitopes was examined. Finally, to confirm observations on globally reported HIV subtype A sequences, we analysed the gag epitope data from our Kenyan, Pakistani, and Afghan cohorts, where both cohort-specific gene epitope variability and HLA restriction profiles of gag epitopes were examined.

Results

The most recent common ancestor of the HIV subtype A epidemic was estimated to be 1956±1. A period of exponential growth began about 1980 and lasted for approximately 7 years, stabilized for 15 years, declined for 2–3 years, then stabilized again from about 2004. During the course of evolution, a gradual increase in genomic variability was observed that peaked in 2005–2010. We observed that the number of point mutations and novel epitopes in gag also peaked concurrently during 2005–2010.

Conclusion

It appears that as the HIV subtype A epidemic spread globally, changing population immunogenetic pressures may have played a role in steering immune-evolution of this subtype in new directions. This trend is apparent in the genomic variability and epitope diversity of HIV-1 subtype A gag sequences.

Population-specific evolution of HIV Gag epitopes in genetically diverged patients

BACKGROUND

Under the host selection pressure HIV evolves rapidly to override crucial steps in the antigen presentation pathway. This allows the virus to escape binding and recognition by cytotoxic T lymphocytes. Selection pressures on HIV can be unique depending on the immunogenetics of host populations. It is therefore logical to hypothesize that the virus evolving in a given population will carry signature mutations that will allow it to survive in that particular host milieu.

OBJECTIVES

The aim of this study was to perform a comparative analysis of HIV-1 Gag subtype A sequences from two genetically diverged populations, namely, Kenyan and Pakistani. We analyzed unique mutations that could intercept the antigen processing pathway and potentially change the repertoire of Gag epitopes in each study group.

METHODS

Twenty-nine Kenyan and 56 Pakistani samples from HIV-1 subtype A-infected patients were used in this study. The HIV-1 gag region p24 and p2p7p1p6 was sequenced and mutations affecting proteasomal degradation, TAP binding, HLA binding and CTL epitope generation, were analyzed using the in silico softwares NetChop and MAPPP, TAPPred, nHLAPred and CTLPred, respectively.

RESULTS

Certain mutations unique to either Pakistani or Kenyan patients were observed to affect sites for proteasomal degradation, TAP binding, and HLA binding. As a consequence of these mutations, epitope pattern in these populations was altered.

CONCLUSION

Unique selection pressures can steer the direction of viral epitope evolution in the host populations. Population-specific HIV epitopes have to be taken into account while designing treatment as well as vaccine for HIV.

HIV-1 Tat B-cell epitope vaccination was ineffectual in preventing viral rebound after ART cessation: HIV rebound with current ART appears to be due to infection with new endogenous founder virus and not to resurgence of pre-existing Tat-dependent viremia

CD4 T cell activation, essential for productive HIV infection, is provided initially in acute HIV infection by innate immune system secretion of activating cytokines. This cytokine response wanes with time and long-term activation of CD4 cells is maintained by HIV Tat protein secreted by HIV infected cells. Structured treatment interruption (STI) in well-controlled antiretroviral-treated (ART) subjects was explored a decade ago with a consensus finding that, in most subjects, HIV levels rebounded within four weeks to pre-ART levels. Based on these observations we initiated a randomized placebo-controlled study of a universal anti-Tat epitope vaccine, TUTI-16, to determine if immunological blockade of Tat would prevent HIV rebound after ART cessation. TUTI-16 immunization was safe, with predominantly mild local and systemic injection-related adverse reactions. TUTI-16 was also immunogenic, with high levels of anti-Tat antibodies compared with levels previously shown to reduce HIV replication in vivo. Of 21 subjects analyzed, 13 (62%) had HIV rebounds vs. 8 (38%) that remained aviremia, but this distribution was not vaccine-related (p = 0.61 log-rank (Mantel-Cox) test), nullifying our hypothesis that anti-Tat antibodies would block rebound of Tat-dependent set-point HIV viremia after ART cessation. Our present findings are consistent with recent molecular findings that rebounding virus following STI is homogeneous and unrelated to previous circulating HIV, suggesting that rebounding HIV represents new founder virus, akin to the original acute HIV infection. We propose, therefore, that STI may have potential as a practical and economical approach to testing the safety and efficacy of candidate prophylactic HIV vaccines.

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.

Exploratory clinical studies of a synthetic HIV-1 Tat epitope vaccine in asymptomatic treatment-naïve and antiretroviral-controlled HIV-1 infected subjects plus healthy uninfected subjects

TUTI-16 is a synthetic universal HIV-1 Tat epitope vaccine, designed to induce anti-Tat antibodies that block the function of circulating Tat, an HIV encoded protein secreted by HIV-1 infected cells. Circulating Tat activates CD4 T cells, permitting HIV replication and sustained viremia. Safety, immunogenicity and antiretroviral potential of TUTI-16 were explored in a randomized double-blind dose-escalating study in asymptomatic treatment-naïve HIV-1 infected subjects. TUTI-16 was safe, with mild local and systemic injection-related adverse reactions, but the antibody response was barely detectable. Surprisingly, a highly statistically significant reduction of HIV-1 viral load was found in the lowest 30 μg vaccine dose group (p < 0.01) but not at the higher doses. We posited that an anti-Tat antibody response below the limit of detection inhibited HIV viral load at this dose, an effect nullified at higher vaccine doses by activating cytokines induced by adjuvant components in TUTI-16. To clarify this immunogenicity/activation conundrum open label immunogenicity studies were performed in healthy HIV uninfected and aviremic ART-controlled HIV-infected subjects. These established that (1) healthy HIV negative subjects had robust antibody responses, maximal with 1 mg TUTI-16, (2) ART-controlled aviremic HIV infected subjects had similarly robust antibody responses, and (3) adjuvant-induced increases of HIV viral load did not occur in the presence of ART. These studies provided us a basis for the design of a protocol to explore the therapeutic potential of TUTI-16 vaccination to provide drug free control of HIV-1 viremia.

Molecular characterization of HIV-1 CRF01_AE in Mekong Delta, Vietnam, and impact of T-cell epitope mutations on HLA recognition (ANRS 12159)

Background

To date, 11 HIV-1 subtypes and 48 circulating recombinant forms have been described worldwide. The underlying reason why their distribution is so heterogeneous is not clear. Host genetic factors could partly explain this distribution. The aim of this study was to describe HIV-1 strains circulating in an unexplored area of Mekong Delta, Vietnam, and to assess the impact of optimal epitope mutations on HLA binding.

Methods

We recruited 125 chronically antiretroviral-naive HIV-1-infected subjects from five cities in the Mekong Delta. We performed high-resolution DNA typing of HLA class I alleles, sequencing of Gag and RT-Prot genes and phylogenetic analysis of the strains. Epitope mutations were analyzed in patients bearing the HLA allele restricting the studied epitope. Optimal wild-type epitopes from the Los Alamos database were used as reference. T-cell epitope recognition was predicted using the immune epitope database tool according to three different scores involved in antigen processing (TAP and proteasome scores) and HLA binding (MHC score).

Results

All sequences clustered with CRF01_AE. HLA class I genotyping showed the predominance of Asian alleles as A*11:01 and B*46:01 with a Vietnamese specificity held by two different haplotypes. The percentage of homology between Mekong and B consensus HIV-1 sequences was above 85%. Divergent epitopes had TAP and proteasome scores comparable with wild-type epitopes. MHC scores were significantly lower in divergent epitopes with a mean of 2.4 (±0.9) versus 2 (±0.7) in non-divergent ones (p<0.0001).

Conclusions

Our study confirms the wide predominance of CRF01_AE in the Mekong Delta where patients harbor a specific HLA pattern. Moreover, it demonstrates the lower MHC binding affinity among divergent epitopes. This weak immune pressure combined with a narrow genetic diversity favors immune escape and could explain why CRF01_AE is still predominant in Vietnam, particularly in the Mekong area.

MHC heterozygote advantage in simian immunodeficiency virus-infected Mauritian cynomolgus macaques

The importance of a broad CD8 T lymphocyte (CD8-TL) immune response to HIV is unknown. Ex vivo measurements of immunological activity directed at a limited number of defined epitopes provide an incomplete portrait of the actual immune response. We examined viral loads in simian immunodeficiency virus (SIV)-infected major histocompatibility complex (MHC)-homozygous and MHC-heterozygous Mauritian cynomolgus macaques. Chronic viremia in MHC-homozygous macaques was 80 times that in MHC-heterozygous macaques. Virus from MHC-homozygous macaques accumulated 11 to 14 variants, consistent with escape from CD8-TL responses after 1 year of SIV infection. The pattern of mutations detected in MHC-heterozygous macaques suggests that their epitope-specific CD8-TL responses are a composite of those present in their MHC-homozygous counterparts. These results provide the clearest example of MHC heterozygote advantage among individuals infected with the same immunodeficiency virus strain, suggesting that broad recognition of multiple CD8-TL epitopes should be a key feature of HIV vaccines.

A universal anti-HIV-1 Tat epitope vaccine that is fully synthetic and self-adjuvanting

Circulating HIV-1 Tat protein is essential for maintenance of the chronic HIV replication that predicates both HIV transmission and clinical progression to AIDS/death. A synthetic anti-Tat epitope vaccine (TUTI-16) was designed to induce neutralizing antibodies to Tat and, hopefully, provide immunological control of HIV replication. TUTI-16 is composed of (1) a conserved Tat B cell epitope (Tat 4-12), rendered universal by introducing known variant amino acids at variable positions 7, 9 and 12 during solid phase synthesis, (2) a promiscuous T helper sequence and (3) a lipopeptide toll-like receptor 2 (TLR2) agonist. TUTI-16 induced high titer antibodies against all 8 known variants of the Tat epitope.

Timing constraints of in vivo gag mutations during primary HIV-1 subtype C infection

Background

Aiming to answer the broad question “When does mutation occur?” this study examined the time of appearance, dominance, and completeness of in vivo Gag mutations in primary HIV-1 subtype C infection.

Methods

A primary HIV-1C infection cohort comprised of 8 acutely and 34 recently infected subjects were followed frequently up to 500 days post-seroconversion (p/s). Gag mutations were analyzed by employing single-genome amplification and direct sequencing. Gag mutations were determined in relation to the estimated time of seroconversion. Time of appearance, dominance, and completeness was compared for different types of in vivo Gag mutations.

Results

Reverse mutations to the wild type appeared at a median (IQR) of 62 (44;139) days p/s, while escape mutations from the wild type appeared at 234 (169;326) days p/s (p<0.001). Within the subset of mutations that became dominant, reverse and escape mutations appeared at 54 (30;78) days p/s and 104 (47;198) days p/s, respectively (p<0.001). Among the mutations that reached completeness, reverse and escape mutations appeared at 54 (30;78) days p/s and 90 (44;196) days p/s, respectively (p = 0.006). Time of dominance for reverse mutations to and escape mutations from the wild type was 58 (44;105) days p/s and 219 (90;326) days p/s, respectively (p<0.001). Time of completeness for reverse and escape mutations was 152 (100;176) days p/s and 243 (101;370) days p/s, respectively (p = 0.001). Fitting a Cox proportional hazards model with frailties confirmed a significantly earlier time of appearance (hazard ratio (HR): 2.6; 95% CI: 2.3–3.0), dominance (4.8 (3.4–6.8)), and completeness (3.6 (2.3–5.5)) of reverse mutations to the wild type Gag than escape mutations from the wild type. Some complex mutational pathways in Gag included sequential series of reversions and escapes.

Conclusions

The study identified the timing of different types of in vivo Gag mutations in primary HIV-1 subtype C infection in relation to the estimated time of seroconversion. Overall, the in vivo reverse mutations to the wild type occurred significantly earlier than escape mutations from the wild type. This shorter time to incidence of reverse mutations remained in the subsets of in vivo Gag mutations that reached dominance or completeness.