HLA-driven convergence of HIV-1 viral subtypes B and F toward the adaptation to immune responses in human populations

Genetic divergence of HIV-1 B subtype in Italy over the years 2003-2016 and impact on CTL escape prevalence

HIV-1 is characterized by high genetic variability, with implications for spread, and immune-escape selection. Here, the genetic modification of HIV-1 B subtype over time was evaluated on 3,328 pol and 1,152 V3 sequences belonging to B subtype and collected from individuals diagnosed in Italy between 2003 and 2016. Sequences were analyzed for genetic-distance from consensus-B (Tajima-Nei), non-synonymous and synonymous rates (dN and dS), CTL escapes, and intra-host evolution over four time-spans (2003–2006, 2007–2009, 2010–2012, 2013–2016). Genetic-distance increased over time for both pol and V3 sequences (P < 0.0001 and 0.0003). Similar results were obtained for dN and dS. Entropy-value significantly increased at 16 pol and two V3 amino acid positions. Seven of them were CTL escape positions (protease: 71; reverse-transcriptase: 35, 162, 177, 202, 207, 211). Sequences with ≥3 CTL escapes increased from 36.1% in 2003–2006 to 54.0% in 2013–2016 (P < 0.0001), and showed better intra-host adaptation than those containing ≤2 CTL escapes (intra-host evolution: 3.0 × 10⁻³ [2.9 × 10⁻³–3.1 × 10⁻³] vs. 4.3 × 10⁻³ [4.0 × 10⁻³–5.0 × 10⁻³], P[LRT] < 0.0001[21.09]). These data provide evidence of still ongoing modifications, involving CTL escape mutations, in circulating HIV-1 B subtype in Italy. These modifications might affect the process of HIV-1 adaptation to the host, as suggested by the slow intra-host evolution characterizing viruses with a high number of CTL escapes.

Next-generation sequencing analyses of the emergence and maintenance of mutations in CTL epitopes in HIV controllers with differential viremia control

Background

Despite the low level of viral replication in HIV controllers (HICs), studies have reported viral mutations related to escape from cytotoxic T-lymphocyte (CTL) response in HIV-1 plasma sequences. Thus, evaluating the dynamics of the emergence of CTL-escape mutants in HICs reservoirs is important for understanding viremia control. To analyze the HIV-1 mutational profile and dynamics of CTL-escape mutants in HICs, we selected 11 long-term non-progressor individuals and divided them into the following groups: (1) viremic controllers (VCs; n = 5) and (2) elite controllers (ECs; n = 6). For each individual, we used HIV-1 proviral DNA from PBMCs related to earliest (V_E) and latest (V_L) visits to obtain gag and nef sequences using the Illumina HiSeq system. The consensus of each mapped gene was used to assess viral divergence, and next-generation sequencing data were employed to identify SNPs and variations within and flanking CTL epitopes.

Results

Divergence analysis showed higher values for nef compared to gag among the HICs. EC and VC groups showed similar divergence rates for both genes. Analysis of the number of SNPs showed that VCs present more variability in both genes. Synonymous/non-synonymous mutation ratios were < 1 for gag among ECs and for nef among ECs and VCs, exhibiting a predominance of non-synonymous mutations. Such mutations were observed in regions encoding CTL-restricted epitopes in all individuals. All ECs presented non-synonymous mutations in CTL epitopes but generally at low frequency (< 1%); all VCs showed a high number of mutations, with significant frequency changes between V_E and V_L visits. A higher frequency of internal mutations was observed for gag epitopes, with significant changes across visits compared to Nef epitopes, indicating a pattern associated with differential genetic pressure.

Conclusions

The high genetic conservation of HIV-1 gag and nef among ECs indicates that the higher level of viremia control restricts the evolution of both genes. Although viral replication levels in HICs are low or undetectable, all individuals exhibited CTL epitope mutations in proviral gag and nef variants, indicating that potential CTL escape mutants are present in HIC reservoirs and that situations leading to a disequilibrium of the host-virus relationship can result in the spread of CTL-escape variants.

Electronic supplementary material

The online version of this article (10.1186/s12977-018-0444-z) contains supplementary material, which is available to authorized users.

Computational comparison of availability in CTL/gag epitopes among patients with acute and chronic HIV-1 infection

BACKGROUND

Recent studies indicate that there is selection bias for transmission of viral polymorphisms associated with higher viral fitness. Furthermore, after transmission and before a specific immune response is mounted in the recipient, the virus undergoes a number of reversions which allow an increase in their replicative capacity. These aspects, and others, affect the viral population characteristic of early acute infection.

METHODS

160 singlegag-gene amplifications were obtained by limiting-dilution RT-PCR from plasma samples of 8 ARV-naïve patients with early acute infection (<30 days, 22 days average) and 8 ARV-naive patients with approximately a year of infection (10 amplicons per patient). Sanger sequencing and NGS SMRT technology (Pacific Biosciences) were implemented to sequence the amplicons. Phylogenetic analysis was performed by using MEGA 6.06. HLA-I (A and B) typing was performed by SSOP-PCR method. The chromatograms were analyzed with Sequencher 4.10. Epitopes and immune-proteosomal cleavages prediction was performed with CBS prediction server for the 30 HLA-A and -B alleles most prevalent in our population with peptide lengths from 8 to 14 mer. Cytotoxic response prediction was performed by using IEDB Analysis Resource.

RESULTS

After implementing epitope prediction analysis, we identified a total number of 325 possible viral epitopes present in two or more acute or chronic patients. 60.3% (n = 196) of them were present only in acute infection (prevalent acute epitopes) while 39.7% (n = 129) were present only in chronic infection (prevalent chronic epitopes). Within p24, the difference was equally dramatic with 59.4% (79/133) being acute epitopes (p < 0.05). This is consistent with progressive viral adaptation to immune response in time and further supported by the fact that cytotoxic responses prediction showed that acute epitopes are more likely to generate immune response than chronic epitopes. Interestingly, only 27.5% of acute epitopes match the population-level consensus sequence of the virus.

CONCLUSIONS

Our results indicate that certain non-consensus viral residues might be transmitted more frequently than consensus-residues when located in immunological relevant positions (epitopes). This observation might be relevant to the rationale behind development of an effective vaccineto reduce viral reservoir and induce functional cure of HIV infection based in prevalent acute epitopes.

Virus-Host Gene Interactions Define HIV-1 Disease Progression

In this chapter, we will review recent research on the virology of HIV-1 transmission and the impact of the transmitted virus genotype on subsequent disease progression. In most instances of HIV-1 sexual transmission, a single genetic variant, or a very limited number of variants from the diverse viral quasi-species present in the transmitting partner establishes systemic infection. Transmission involves both stochastic and selective processes, such that in general a minority variant in the donor is transmitted. While there is clear evidence for selection, the biological properties that mediate transmission remain incompletely defined. Nevertheless, the genotype of the transmitted founder virus, which reflects prior exposure to and escape from host immune responses, clearly influences disease progression. Some escape mutations impact replicative capacity, while others effectively cloak the virus from the newly infected host's immune response by preventing recognition. It is the balance between the impact of escape mutations on viral fitness and susceptibility to the host immunogenetics that defines HIV-1 disease progression.

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

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

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

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.

Genotypic and functional impact of HIV-1 adaptation to its host population during the North American epidemic

HLA-restricted immune escape mutations that persist following HIV transmission could gradually spread through the viral population, thereby compromising host antiviral immunity as the epidemic progresses. To assess the extent and phenotypic impact of this phenomenon in an immunogenetically diverse population, we genotypically and functionally compared linked HLA and HIV (Gag/Nef) sequences from 358 historic (1979-1989) and 382 modern (2000-2011) specimens from four key cities in the North American epidemic (New York, Boston, San Francisco, Vancouver). Inferred HIV phylogenies were star-like, with approximately two-fold greater mean pairwise distances in modern versus historic sequences. The reconstructed epidemic ancestral (founder) HIV sequence was essentially identical to the North American subtype B consensus. Consistent with gradual diversification of a "consensus-like" founder virus, the median "background" frequencies of individual HLA-associated polymorphisms in HIV (in individuals lacking the restricting HLA[s]) were ∼ 2-fold higher in modern versus historic HIV sequences, though these remained notably low overall (e.g. in Gag, medians were 3.7% in the 2000s versus 2.0% in the 1980s). HIV polymorphisms exhibiting the greatest relative spread were those restricted by protective HLAs. Despite these increases, when HIV sequences were analyzed as a whole, their total average burden of polymorphisms that were "pre-adapted" to the average host HLA profile was only ∼ 2% greater in modern versus historic eras. Furthermore, HLA-associated polymorphisms identified in historic HIV sequences were consistent with those detectable today, with none identified that could explain the few HIV codons where the inferred epidemic ancestor differed from the modern consensus. Results are therefore consistent with slow HIV adaptation to HLA, but at a rate unlikely to yield imminent negative implications for cellular immunity, at least in North America. Intriguingly, temporal changes in protein activity of patient-derived Nef (though not Gag) sequences were observed, suggesting functional implications of population-level HIV evolution on certain viral proteins.

HIV, HBV, and HCV molecular epidemiology among trans (transvestites, transsexuals, and transgender) sex workers in Argentina

Commercial sex work is frequent among male-to-female transvestites, transsexuals and transgenders in Argentina, leading to high susceptibility to HIV, HBV, and HCV among other sexually transmitted infections. In a global context of scarce data on the trans sex workers population, this study was aimed to study the genomic characterization of these viruses. Plasma presence of HIV, HBV, and HCV genomic material was evaluated in samples from 273 trans sex workers. Genomic sequences of HIV-gag, pol, and vif-vpu genes, HBV-S gene, and HCV-5'UT and NS5B genes were obtained. Molecular characterization involved phylogenetic analysis and several in silico tools. Resistance-associated mutations in HIV and HBV pol genes were also analyzed. The HIV genomic characterization in 62 trans sex workers samples showed that 54.8% of the isolates corresponded to BF intersubtype recombinants, and 38.7% to subtype B. The remaining were classified as subtypes C (4.8%) and A (1.6%). HBV and HCV co-infection prevalence among HIV positive trans sex workers yielded rates of 3.2% and 6.5% respectively. Drug resistance-associated mutations were found in 12/62 (19%) HIV pol sequences, but none among HBV. Based on phylogenetic relationships, HIV isolates characterized as subtypes BF and B appeared intermingled with those from other high-risk groups. Despite trans sex workers declared not to have received antiviral treatment, complex drug resistance-associated mutation patterns were found in several HIV isolates. Planned prevention, screening, and treatment are needed to reduce further transmission and morbidity.

T cells target APOBEC3 proteins in human immunodeficiency virus type 1-infected humans and simian immunodeficiency virus-infected Indian rhesus macaques

APOBEC3 proteins mediate potent antiretroviral activity by hypermutating the retroviral genome during reverse transcription. To counteract APOBEC3 and gain a replicative advantage, lentiviruses such as human immunodeficiency virus type 1 (HIV-1) and simian immunodeficiency virus (SIV) have evolved the Vif protein, which targets APOBEC3 proteins for proteasomal degradation. However, the proteasome plays a critical role in the generation of T cell peptide epitopes. Whether Vif-mediated destruction of APOBEC3 proteins leads to the generation and presentation of APOBEC3-derived T cell epitopes on the surfaces of lentivirus-infected cells remains unknown. Here, using peptides derived from multiple Vif-sensitive APOBEC3 proteins, we identified APOBEC3-specific T cell responses in both HIV-1-infected patients and SIV-infected rhesus macaques. These results raise the possibility that these T cell responses may be part of the larger antiretroviral immune response.

Estimation of HIV-testing rates to maximize early diagnosis-derived benefits at the individual and population level

Background

In HIV infection, initiation of treatment is associated with improved clinical outcom and reduced rate of sexual transmission. However, difficulty in detecting infection in early stages impairs those benefits. We determined the minimum testing rate that maximizes benefits derived from early diagnosis.

Methods

We developed a mathematical model of HIV infection, diagnosis and treatment that allows studying both diagnosed and undiagnosed populations, as well as determining the impact of modifying time to diagnosis and testing rates. The model’s external consistency was assessed by estimating time to AIDS and death in absence of treatment as well as by estimating age-dependent mortality rates during treatment, and comparing them with data previously reported from CASCADE and DHCS cohorts.

Results

In our model, life expectancy of patients diagnosed before 8 years post infection is the same as HIV-negative population. After this time point, age at death is significantly dependent on diagnosis delay but initiation of treatment increases life expectancy to similar levels as HIV-negative population. Early mortality during HAART is dependent on treatment CD4 threshold until 6 years post infection and becomes dependent on diagnosis delay after 6 years post infection. By modifying testing rates, we estimate that an annual testing rate of 20% leads to diagnosis of 90% of infected individuals within the first 8.2 years of infection and that current testing rate in middle-high income settings stands close to 10%. In addition, many differences between low-income and middle-high incomes can be predicted by solely modifying the diagnosis delay.

Conclusions

To increase testing rate of undiagnosed HIV population by two-fold in middle-high income settings will minimize early mortality during initiation of treatment and global mortality rate as well as maximize life expectancy. Our results highlight the impact of achieving early diagnosis and the importance of strongly work on improving HIV testing rates.

Pyrosequencing reveals restricted patterns of CD8+ T cell escape-associated compensatory mutations in simian immunodeficiency virus

CD8+ T cells play a major role in the containment of human immunodeficiency virus (HIV) and simian immunodeficiency virus (SIV) replication. CD8+ T cell-driven variations in conserved regions under functional constraints result in diminished viral replicative capacity. While compensatory mutations outside an epitope can restore replicative capacity, the kinetics with which they arise remains unknown. Additionally, certain patterns of linked mutations associated with CD8+ T cell epitope escape in these highly conserved regions may lead to variable levels of viral fitness. Here, we used pyrosequencing to investigate the kinetics and patterns of mutations surrounding the Mamu-A1*00101-bound Gag(181-189)CM9 CD8+ T cell epitope. We obtained more than 400 reads for each sequencing time point, allowing us to confidently detect the emergence of viral variants bearing escape mutations with frequencies as low as 1% of the circulating virus. With this level of detail, we demonstrate that compensatory mutations generally arise concomitantly with Gag(181-189)CM9 escape mutations. We observed distinct patterns of linked flanking mutations, most of which were found downstream of Gag(181-189)CM9. Our data indicate that, whereas Gag(181-189)CM9 escape is much more complex that previously appreciated, it occurs in a coordinated fashion, with very specific patterns of flanking mutations required for immune evasion. This is the first detailed report of the ontogeny of compensatory mutations that allow CD8+ T cell epitope escape in infected individuals.

Loss of HIV-1-derived cytotoxic T lymphocyte epitopes restricted by protective HLA-B alleles during the HIV-1 epidemic

OBJECTIVE AND DESIGN

HIV-1 is known to adapt to the human immune system, leading to accumulation of escape mutations during the course of infection within an individual. Cross-sectional studies have shown an inverse correlation between the prevalence of human leukocyte antigen (HLA) alleles in a population and the number of cytotoxic T lymphocyte (CTL) escape mutations in epitopes restricted by those HLA alleles. Recently, it was demonstrated that at a population level HIV-1 is adapting to the humoral immune response, which is reflected in an increase in resistance to neutralizing antibodies over time. Here we investigated whether adaptations to cellular immunity have also accumulated during the epidemic.

METHODS

We compared the number of CTL epitopes in HIV-1 strains isolated from individuals who seroconverted at the beginning of the HIV-1 epidemic and from individuals who seroconverted in recent calendar time.

RESULTS

The number of CTL epitopes in HIV-1 variants restricted by the most common HLA alleles in the population did not change significantly during the epidemic. In contrast, we found a significant loss of CTL epitopes restricted by HLA-B alleles associated with a low relative hazard of HIV-1 disease progression during the epidemic. Such a loss was not observed for CTL epitopes restricted by HLA-A alleles.

CONCLUSION

Despite the large degree of HLA polymorphism, HIV-1 has accumulated adaptations to CTL responses within 20 years of the epidemic. The fact that such adaptations are driven by the HLA-B molecules that provide best protection against HIV-1 disease progression has important implications for our understanding of HIV evolution.

Analysis of HIV type 1 BF recombinant sequences from South America dates the origin of CRF12_BF to a recombination event in the 1970s

HIV-1 epidemics in South America are believed to have originated in part from the subtype B epidemic initiated in the Caribbean/North America region. However, circulation of BF recombinants in similar proportions was extensively reported. Information currently shows that many BF recombinants share a recombination structure similar to that found in the CRF12_BF. In the present study, analyzing a set of 405 HIV sequences, we identified the most likely origin of the BF epidemic in an early event of recombination. We found that the subtype B epidemics in South America analyzed in the present study were initiated by a founder event that occurred in the early 1970s, a few years after the introduction of these strains in the Americas. Regarding the F/BF recombinant epidemics, by analyzing a subtype F genomic segment within the viral gene gag present in the majority of the BF recombinants, we found evidence of a geographic divergence very soon after the introduction of subtype F strains in South America. Moreover, through analysis of a subtype B segment present in all the CRF12_BF-like recombination structure, we estimated the circulation of the subtype B strain that gave rise to that recombinant structure around the same time period estimated for the introduction of subtype F strains. The HIV epidemics in South America were initiated in part through a founder event driven by subtype B strains coming from the previously established epidemic in the north of the continent. A second introduction driven by subtype F strains is likely to have encountered the incipient subtype B epidemic that soon after their arrival recombined with them, originating the BF epidemic in the region. These results may explain why in South America the majority of F sequences are found as BF recombinants.