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

HIV-1 subtype diversity, phylogenetic analysis and study of drug resistance in strains circulating in the Ural Federal District

Introduction. Ural Federal District (UFD) has been one of the most HIV-affected areas in the Russian Federation during past 20 years. The total number of people living with HIV/AIDS (PLWH) and receiving antiretroviral therapy (ART) exceeds 100,000 (61.7% of all PLWH in the UFD), which creates opportunities for the wide spread of resistant HIV strains.Research aim was to determine the distribution of HIV-1 subtypes, evaluate the genetic heterogeneity of HIV-1 strains, and analyze the prevalence of HIV-1 drug resistance mutations (DRM) and the incidence of acquired resistance to antiretroviral drugs (ARVDs) in PLWH receiving ART in the UFD.Materials and methods. 223 patients receiving ART at stage III–IV of HIV infection living in the UFD were examined. To determine the subtypes and the DRM in the HIV-1 pol gene, molecular genetic studies were performed using the AmpliSense® HIV-Resist-Seq kit by Sanger sequencing on the Applied Biosystems 3500 Genetic Analyzer. The genetic heterogeneity was evaluated by calculating the identity of the genome region of the isolated strains in comparison with the genomes of foreign HIV strains, as well as using phylogenetic analysis.Results. In the studied group of patients, 5 subtypes of HIV-1 were identified: subtype A6 prevalence was 91.03%, that of subtype B was 2.69%, 3 recombinant subtypes (CRF03_A6B, CRF02_AG, CRF63_02A6) accounted for 6.28%. Among analyzed HIV-1 strains, 43.9% had a significant genetic similarity (identity of at least 97%) with the strains isolated from patients from neighboring countries (Belarus, Kazakhstan, Kyrgyzstan, Uzbekistan, Lithuania), 35.9% were similar to the strains isolated from patients from far-abroad countries (USA, China, South Korea, Australia, Sweden, Germany). A high heterogeneity of the circulating genetic variants of HIV-1 strains in the territory of the UFD region was established, which is an unfavorable factor for the diagnosis and treatment of HIV. The most common DRMs to both nucleoside reverse transcriptase inhibitors (NRTI) and non-nucleoside reverse transcriptase inhibitors (NNRTI) were detected in 81 specimens (36.3%). NRTI resistance-forming M184V DRM was more common than any other DRM with statistical significance (p = 0,0008) and was detected in 88 specimens (39.5%).Conclusion. In the subtype structure of HIV-1, the dominant subtype was subtype A6, the most common in the countries that were formerly part of the USSR. The heterogeneity of the HIV-1 strains circulating in the UFD suggests that HIV-1 infection continues to be introduced into the UFD from populations outside the Russian Federation. The findings confirm the high prevalence of DRMs (62.8%) and secondary drug resistance of HIV-1 (60.1%) among PLWH in the territory of the UFD. At the same time, high-level resistance was detected in 56.5% of patients, which requires increasing the coverage of HIV resistance testing, including the introduction of monitoring for primary resistance, in order to optimize first-line ART regimens. 

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.

HIV-1 p24Gag adaptation to modern and archaic HLA-allele frequency differences in ethnic groups contributes to viral subtype diversification

Pathogen-driven selection and past interbreeding with archaic human lineages have resulted in differences in human leukocyte antigen (HLA)-allele frequencies between modern human populations. Whether or not this variation affects pathogen subtype diversification is unknown. Here we show a strong positive correlation between ethnic diversity in African countries and both human immunodeficiency virus (HIV)-1 p24gag and subtype diversity. We demonstrate that ethnic HLA-allele differences between populations have influenced HIV-1 subtype diversification as the virus adapted to escape common antiviral immune responses. The evolution of HIV Subtype B (HIV-B), which does not appear to be indigenous to Africa, is strongly affected by immune responses associated with Eurasian HLA variants acquired through adaptive introgression from Neanderthals and Denisovans. Furthermore, we show that the increasing and disproportionate number of HIV-infections among African Americans in the USA drive HIV-B evolution towards an Africa-centric HIV-1 state. Similar adaptation of other pathogens to HLA variants common in affected populations is likely.

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.

Evidence for a recombinant origin of HIV-1 Group M from genomic variation

Reconstructing the early dynamics of the HIV-1 pandemic can provide crucial insights into the socioeconomic drivers of emerging infectious diseases in human populations, including the roles of urbanization and transportation networks. Current evidence indicates that the global pandemic comprising almost entirely of HIV-1/M originated around the 1920s in central Africa. However, these estimates are based on molecular clock estimates that are assumed to apply uniformly across the virus genome. There is growing evidence that recombination has played a significant role in the early history of the HIV-1 pandemic, such that different regions of the HIV-1 genome have different evolutionary histories. In this study, we have conducted a dated-tip analysis of all near full-length HIV-1/M genome sequences that were published in the GenBank database. We used a sliding window approach similar to the 'bootscanning' method for detecting breakpoints in inter-subtype recombinant sequences. We found evidence of substantial variation in estimated root dates among windows, with an estimated mean time to the most recent common ancestor of 1922. Estimates were significantly autocorrelated, which was more consistent with an early recombination event than with stochastic error variation in phylogenetic reconstruction and dating analyses. A piecewise regression analysis supported the existence of at least one recombination breakpoint in the HIV-1/M genome with interval-specific means around 1929 and 1913, respectively. This analysis demonstrates that a sliding window approach can accommodate early recombination events outside the established nomenclature of HIV-1/M subtypes, although it is difficult to incorporate the earliest available samples due to their limited genome coverage.

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.

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.

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

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

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

The impact of HIV-1 within-host evolution on transmission dynamics

The adaptive potential of HIV-1 is a vital mechanism to evade host immune responses and antiviral treatment. However, high evolutionary rates during persistent infection can impair transmission efficiency and alter disease progression in the new host, resulting in a delicate trade-off between within-host virulence and between-host infectiousness. This trade-off is visible in the disparity in evolutionary rates at within-host and between-host levels, and preferential transmission of ancestral donor viruses. Understanding the impact of within-host evolution for epidemiological studies is essential for the design of preventive and therapeutic measures. Herein, we review recent theoretical and experimental work that generated new insights into the complex link between within-host evolution and between-host fitness, revealing temporal and selective processes underlying the structure and dynamics of HIV-1 transmission.