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

Tomohiro Akahoshi; Hiroyuki Gatanaga; Nozomi Kuse; Takayuki Chikata; Madoka Koyanagi; Naoki Ishizuka; Chanson J. Brumme; Hayato Murakoshi; Zabrina L. Brumme; Shinichi Oka; Masafumi Takiguchi

doi:10.1371/journal.ppat.1009177

Abstract

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.

Author summary

HIV-1 strains harboring immune escape mutations can accumulate in circulation, but it remains unclear to what extent these 'escaped' HIV-1 strains continue to evolve under ongoing population-level immune pressures. We investigated population-level changes at HIV-1 Reverse Transcriptase codon 135 (RT135), which is under pressure by T cells restricted by HLA-B*51:01 and B*52:01, highly frequent alleles in Japan. While threonine initially accumulated at RT135, RT135L has subsequently increased markedly. Our findings revealed that RT135V selection by HLA-B*52:01-restricted T-cells led to the creation of a new epitope restricted by HLA-C*12:02, an allele in strong linkage disequilibrium with HLA-B*52:01. HLA-C*12:02-restricted T cells in turn suppressed replication of RT135T virus and selected RT135L. Notably, population-level shifts at this codon are particular to Japan, where HLA-B*52:01-C*12:02 represents the most prevalent HLA haplotype. Our findings highlight multiple virus-specific T cells as dynamic drivers of population-level − and host population-specific − HIV-1 evolution over the long term.

Citation: Akahoshi T, Gatanaga H, Kuse N, Chikata T, Koyanagi M, Ishizuka N, et al. (2020) T-cell responses to sequentially emerging viral escape mutants shape long-term HIV-1 population dynamics. PLoS Pathog 16(12): e1009177. https://doi.org/10.1371/journal.ppat.1009177

Editor: Charles R. M. Bangham, Imperial College London, UNITED KINGDOM

Received: September 25, 2020; Accepted: November 18, 2020; Published: December 28, 2020

Copyright: © 2020 Akahoshi et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: All relevant data are within the manuscript and its Supporting Information files.

Funding: This research was supported in part by a grant-in-aid for AIDS Research (15fk0410019h0001, 16fk0410202h0002, and 17fk0410302h0003) from AMED to MT and by a grant-in-aid (26293240 and 26860335) for scientific research from the Ministry of Education, Science, Sports and Culture, Japan, to MT and TA, and by a project grant from the Canadian Institutes of Health Research (PJT-148621) to ZLB. ZLB is supported by a Scholar Award from the Michael Smith Foundation for Health Research. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

HIV-1 mutations, many of which are selected by HIV-1-specific cytotoxic T lymphocytes (CTLs), can accumulate at the population level [1–5]. Some of these mutations may revert to wild-type upon transmission if the new hosts cannot elicit epitope-specific CTLs and/or if these mutant viruses have reduced viral fitness [6–9]. Many mutations however do not carry a fitness cost, and can therefore persist and accumulate [3,10–12]. Studies of HLA-associated polymorphisms in HIV-1 suggest that most of these are selected by CTLs [13–17], and a recent report that compared historic and recent HIV-1 strains in North America demonstrated that, on average, circulating HLA-associated mutation frequencies have approximately doubled since the 1980s [18,19]. Thus, at the population level, HIV-1 strains are continually evolving to adapt to antiviral HLA-mediated immune responses in their host population(s).

As immune escape mutations compromise the recognition of HIV-1 by virus-specific T cells, their accumulation at the population level compromises T cell-mediated control of HIV-1 [10,20–23]. These negative implications also extend to protective HLA alleles such as HLA-B57, HLA-B27, and HLA-B51 [6,24–28]. For example, the RT I135T mutation selected by HLA-B*51:01-restricted T cells specific for the TI8 epitope (RT codons 128–135) completely abrogates HIV-1 recognition by these cells [2,29]. Consequently, although HLA-B*51:01 was a protective allele in Japan more than 30 years ago [29,30], HLA-B*51:01 is no longer protective in individuals infected more recently [31,32]. This is likely because the RT135T mutation accumulated to particularly high levels in Japan [2,33].

In Japan however, RT135 mutations also frequently occur in individuals harboring HLA-B*52:01[33], which is common in the region [34]. Specifically, RT135T is predominantly observed in HLA-B*51:01⁺ individuals whereas RT135T, RT135L, and RT135V are observed at relatively equal frequencies in HLA-B*52:01⁺ individuals [33]. Curiously however, while HLA-B*51:01-restricted and HLA-B*52:01-restricted TI8-specific CTLs fail to recognize target cells infected with RT135T or RT135L mutant virus, they can still recognize cells infected with RT135V virus [2,33]. The mechanisms underlying differential RT135 mutation patterns in HLA-B*51:01⁺ and HLA-B*52:01⁺ individuals however remain unclear.

Two hypotheses are plausible. Undiscovered HLA-B*52:01-restricted epitopes in the region could explain RT135 mutations, but our previous intensive searches did not identify any additional HLA-B*52:01-restricted epitopes in this region that contained RT135 [35]. Alternatively, T-cells restricted by allele(s) in strong linkage disequilibrium with HLA-B*52:01, that target undiscovered epitopes in the region, could select RT135 mutations. In Japan, HLA-C*12:02 is in very strong linkage disequilibrium (LD) with HLA-B*52:01 (LD = 0.99) [34,36], and the HLA-B*52:01-C*12:02 haplotype, at 20% frequency, represents the most prevalent haplotype in Japan. As such, we hypothesized that HLA-C*12:02-restricted T cells may contribute to the selection of different mutations in individuals carrying this haplotype.

In the present study, we searched for novel HLA-C*12:02-restricted T cell epitopes covering RT135 and investigated the role of these T cells in selecting mutations at RT135 in individuals harboring the HLA-B*52:01-C*12:02 haplotype. We further explored the possibility that such T cells could underpin the dynamic changes in population-level RT135 mutation frequencies that have occurred over the course of the HIV-1 epidemic in Japan.

Results

Accumulation of RT135V mutation in HLA-B52:01⁺HLA-C12:02⁺ hemophiliacs with HIV

To exclude the effect of ongoing domestic HIV-1 transmission on the accumulation of mutant viruses in the population, we analyzed RT135 amino acid frequencies in 95 Japanese hemophiliacs who acquired HIV-1 via transfusion with contaminated blood products imported from the USA around 1983 [37,38]. Overall, 49% acquired HIV-1 carrying RT135I (wild-type virus), but this percentage rose to 87% when restricting to B*51:01^-/B*52:01^- hemophiliacs (Fig 1A). This supports our previous observation that the majority of hemophiliacs acquired wild-type (RT135I) virus [2]. The frequency of RT135V was significantly higher among HLA-B*52:01⁺ individuals than those lacking this allele (p = 8.9x10^-5, q = 5.3x10^-4), whereas RT135T was significantly higher among HLA-B*51:01⁺ individuals than those lacking this allele (p = 3.8X10^-5, q = 3.4x10^-4) (Fig 1A). However, no other RT135 mutations were enriched among HLA-B*52:01⁺ or HLA-B*51:01⁺ individuals. These results clearly demonstrate that RT135V and RT135T mutations were overrepresented among HLA-B*52:01⁺ and HLA-B*51:01⁺ hemophiliacs, respectively.

Download:

Fig 1. Amino acid variation at RT135 in Japanese individuals with chronic HIV-1.

A. HIV-1 RT135 amino acid variation in 95 Japanese hemophiliac individuals with HIV-1. RT135 amino acid frequencies were compared between groups using fisher's exact test (two-tailed). Multiple tests were addressed using q values (indicated below the P values). B. RT135 amino acid variation in 83 Japanese male individuals diagnosed with HIV-1 before 1997. RT135 amino acid frequencies were compared as in Fig 1A.

https://doi.org/10.1371/journal.ppat.1009177.g001

We next analyzed RT135 mutation frequencies in non-hemophiliac Japanese individuals who had been diagnosed with HIV-1 before 1997, and compared these to the frequencies observed in the hemophiliacs. The results revealed that RT135T had accumulated to comparable levels in HLA-B*51:01⁺ non-hemophiliac and hemophiliac individuals (Fig 1B). In contrast, RT135V mutation frequencies were markedly lower in HLA-B*52:01⁺ non-hemophiliac individuals compared to HLA-B*52:01⁺ hemophiliacs (Fig 1B). Taken together our results suggest that, not only does RT135 mutation selection differ between HLA-B*52:01⁺ and HLA-B*51:01⁺ individuals, but that it may also differ between HLA-B*52:01⁺ hemophiliac and non-hemophiliac individuals, where RT135 mutant frequencies in the latter group may be further influenced by ongoing domestic HIV-1 transmission over time.

Generation of novel HLA-C12:02-restricted epitope after selection of RT135V mutation by HLA-B52:01-restricted TI8-specific CTLs

Only one nucleotide substitution between I (ata) and T (aca) or V (gta) at RT135 was observed in hemophiliacs and non-hemophiliacs (S1 Fig and S1 Table), suggesting that RT135T and RT135V each evolved from RT135I. It is not clear however why RT135V predominated in HLA-B*52:01⁺ hemophiliacs while RT135T predominated in HLA-B*51:01⁺ ones. We hypothesized that the differential RT135 mutation patterns observed in HLA-B*51:01⁺ and HLA-B*52:01⁺ individuals may be mediated by HLA-C*12:02-restricted T cell pressures in the latter group, since HLA-B*52:01 and HLA-C*12:02 are in strong linkage disequilibrium. We therefore sought to identify novel HLA-C*12:02-restricted CTL epitopes spanning RT135 using a cocktail of overlapping 11-mer peptides including RT135 (Pol cocktail 15). We identified one individual KI-793 who exhibited weak T cell responses to Pol cocktail 15 (S2A Fig). We further found that cultured T cells stimulated with KN11 peptide exhibited a weak HLA-C*12:02-restricted recognition of KN11 peptide, but no recognition by HLA-B*52:01-restricted cells (S2B Fig). Since this individual harbored RT135V mutant virus, we speculated that his T cells more effectively recognized a peptide containing the RT135V mutation than the wild-type peptide. Indeed, when PBMCs from this individual were cultured after stimulation with KN11 or KN11-10V peptide, the latter induced far stronger responses (S2C Fig). We next analyzed the responses of KN11-10V-induced T cells to truncated peptides of KN11-10V. The T cells recognized TN9-8V most effectively (Fig 2A), identifying this as the optimal HLA-C*12:02-restricted epitope. Additional analysis using an HLA-C*12:02-TN9-8V tetramer revealed a high frequency of HLA-C*12:02-restricted TN9-8V-specific CD8⁺ T cells in PBMCs from KI-793 (1.82% of CD8⁺ T cells) (Figs 2B and S3A). These results together indicate that HIV-1 strains harboring the RT135V mutant contain a novel HLA-C*12:02-restricted TN9-8V T cell epitope.

Download:

Fig 2. Identification and characterization of a novel HLA-C*12:02-restricted epitope TN9-8V.

A. Identification of an optimal HLA-C*12:02-restricted epitope, TN9-8V. Recognition of KN11-10V-specific bulk CD8⁺ T cells from an individual harboring RT135V virus (participant KI-793) by C1R-C1202 cells prepulsed with KN11-10V peptide or its truncated derivatives was analyzed by ICS assay. B. Detection of HLA-C*12:02-restricted TN9-8V-specific CD8⁺ T cells in PBMCs from KI-793 by tetramer staining. C. Frequencies of HLA-B*52:01-restricted TI8-specific and HLA-C*12:02-restricted TN9-8V-specific CD8⁺ T cells in HLA-B*52:01⁺C*12:02⁺ individuals with chronic HIV-1, by staining of TN9-8V/C1202 and TI8/B5201 tetramers. RT135 amino acid variation in these individuals was determined by bulk sequencing. D. Representative participants harboring TN9-8V-specific and TI8-specific CD8⁺ T cells.

https://doi.org/10.1371/journal.ppat.1009177.g002

The results further suggested that TN9-8V-specific CD8⁺ T cells would be elicited in HLA-B*52:01⁺ HLA-C*12:02⁺ individuals who acquired RT135V virus, but not in those who acquired wild-type virus or other mutants at this position. Indeed, TN9-8V-specific CD8⁺ T cells were detected in 5 of 8 HLA-B*52:01⁺HLA-C*12:02⁺ individuals harboring RT135V mutant virus, whereas these T cells were not detected in any of the 7 HLA-B*52:01⁺HLA-C*12:02⁺ individuals harboring wild-type or other mutant virus (Fig 2C). This directly supports the idea that TN9-8V-specific CD8⁺ T cells are specifically elicited in HLA-B*52:01⁺HLA-C*12:02⁺ individuals infected with RT135V virus.

Of note, a previous study showed that TI8-specific HLA-B*52:01-restricted CD8⁺ T cells suppressed replication of the 8V mutant virus less than that of wild-type virus [33], suggesting that HLA-B*52:01-restricted CD8⁺ T cells select RT135V. Indeed, TI8-specific and TN9-8V-specific T cells were detected only in individuals harboring RT135V mutant virus, but not in those harboring other mutant viruses (Fig 2C and 2D). This supports the notion that the C*12:02-TN9-8V epitope is "created" following B*52:01-mediated selection of RT135V in vivo.

Recognition of target cells infected with RT135 mutant viruses by TN9-8V-specific T cells

To clarify the ability of TN9-8V-specific T cells to recognize RT135 variants, we established TN9-8V-specific CD8⁺ T cell lines from 3 individuals (KI-069, KI-661, and KI-793) and assessed the ability of these T cell lines to recognize four peptides (TN9-8T, -8R, -8L, and -8I) versus the cognate TN9-8V. These T cell lines cross-recognized TN9-8T, though to a lesser extent than TN9-8V (Fig 3A). We measured the binding affinity of TN9-8V and the mutant peptides to HLA-C*12:02 using the HLA class I stabilization assay. Among all peptides tested, TN9-8V and TN9-8T exhibited the highest and second highest affinities to HLA-C*12:02, respectively. Of note, these affinities were stronger than those of the known HLA-C*12:02-restricted HIV-1 epitope KY9. In contrast, TN9-8L and TN9-8R did not bind to HLA-C*12:02 (Fig 3B), indicating that RT135L/R mutations critically affect peptide recognition by HLA-C*12:02-specific T cells.

Download:

Fig 3. Recognition of RT135 mutant virus-infected cells by HLA-C*12:02-restricted TN9-8V-specific CTL lines.

A. Recognition of TN9-8V and its mutant peptides by TN9-8V-specific T cell lines derived from HLA-B*52:01⁺C*12:02⁺ individuals with HIV-1. The T cell responses to C1R-C1202 cells prepulsed with TN9-8V or its mutant peptides were analyzed by ICS. B. Binding affinity of TN9-8V and its mutant peptides to HLA-C*12:02 molecules by an HLA class I stabilization assay using RMA-S-C1202 cells. KY9 is a reported HLA-C*12:02-restricted epitope. Statistical analysis was conducted using Student's t-test (two-tailed). Two independent experiments showed the same result. C. Recognition of cells infected with RT135I virus (NL43; NL43-8I) or the mutant virus (NL43-8V, -8T or -8L) by TN9-8V-specific CTL lines. The T-cell responses to 721.221-CD4-C1202 cells infected with NL43-8I or the mutant virus were measured by ICS. The frequencies of p24⁺ cells among 721.221-CD4-C1202 cells infected with NL43-8I, -8V, -8T, and -8L were 72.5%, 85.6%, 73.6%, and 59.8%, respectively, for KI-069 and KI-661, while they were 63.4%, 67.9%, 64.6%, and 40.4%, respectively, for KI-793. See also S4 Fig. D. Suppression of viral replication by TN9-8V-specfic CTL lines, at an E:T ratio of 1:1. The concentration of p24 Ag in the culture supernatants collected on Day 5 was measured by ELISA (left) and the percentage of suppression was calculated (right). Data are shown as the means and ±SD of triplicate assays (A-D).

https://doi.org/10.1371/journal.ppat.1009177.g003

We next investigated the ability of these T cell lines to recognize cells infected with the mutant viruses and to suppress HIV-1 replication. All three TN9-8V-specific CD8⁺ T cell lines recognized 721.221-CD4-C1202 cells infected with RT135V and RT135T mutant virus, though recognition of RT135V was consistently stronger than RT135T. In contrast, recognition of RT135I virus was low or nonexistent (Figs 3C and S4). Moreover, while the T cell lines from KI-069 and KI-793 effectively suppressed the replication of both RT135V and RT135T viruses, they only partially suppressed RT135I virus (KI-069) or failed to suppress it (KI-793) (Fig 3D). These results together suggest that TN9-8V-specific T cells effectively suppress replication of both RT135V and RT135T viruses in vivo.

Population-level shifts in HIV-1 RT135 mutation frequencies Japan over two decades

Among Japanese hemophiliacs, RT135T was predominantly selected in HLA-B*51:01⁺ individuals while RT135V was predominantly selected in HLA-B*52:01⁺ individuals. RT135T was also predominantly selected in HLA-B*51:01⁺ non-hemophiliacs. Somewhat in contrast however, RT135T, RT135L and RT135V occurred at relatively equal frequencies among HLA-B*52:01⁺ non-hemophiliacs [33]. These findings suggest that B*52:01-associated selection pressures on RT135 differed between the hemophiliac and non-hemophiliac cohorts. As the HIV-1 epidemic in Japanese non-hemophiliacs is primarily fueled by ongoing domestic transmission, we investigated trends in RT135 frequencies over time.

To do so, we analyzed RT135 mutation frequencies in 1,478 non-hemophilic Japanese individuals diagnosed with HIV-1 between 1997 and 2015. Analysis was performed in the cohort as a whole, and stratified by allele carriage (HLA-B*51:01⁺B*52:01/C*12:02^-, HLA-B*51:01^-B*52:01/C*12:02⁺, and HLA-B*51:01^-B*52:01/C*12:02^-). RT135 amino acid frequencies were compared across four periods according to year of HIV-1 diagnosis: before 1997 (<1997), 1998–2003, 2004–2009, and 2010–2015 (Fig 4A). The frequencies of V, T, R, M, L, K and I changed significantly across these periods, overall and in stratified analyses (all P<0.001). Overall, the frequency of wild-type RT135I among non-hemophiliacs declined from 25% before 1997 to 12% in 2010–2015, confirming gradual accumulation of RT135 mutants in the population (P<0.0001). Notably, the frequency of RT135L increased during these periods (22% to 44%, P<0.0001), whereas that of RT135T decreased (41% to 25%, P = 0.0005) (S2 Table). Stratification by HLA allele carriage further revealed that RT135L increased in all groups (P<0.0001). In contrast, RT135T decreased in only 2 groups: HLA-B*51:01⁺B*52:01^-/C*12:02^- individuals and HLA-B*51:01^-B*52:01^-/C*12:02^- individuals (P = 0.001) (S2 Table).

Download:

Fig 4. Amino acid variation at RT135 in Japanese and Canadian individuals with subtype B HIV-1 over 2 decades.

A. HIV-1 RT135 amino acid frequencies in 1,478 Japanese non-hemophiliac individuals with chronic HIV-1 subtype B before 1997, and between 1998–2003, 2004–2009, and 2010–2015. The cohort is analyzed overall, and stratified by HLA allele carriage as follows: HLA-B*51:01⁺B*52:01^-C*12:02^- (n = 212); HLA-B*51:01^-B*52:01⁺C*12:02⁺ (n = 302); and HLA-B*51:01^-B*52:01^-C*12:02^- (n = 964). Statistical analysis was conducted using the Cochran-Mantel-Haenszel test (See also S2 Table). B. RT135 amino acid frequencies in 6,258 Canadian individuals with chronic HIV-1 subtype B recruited from 2000 to 2015. "X" indicates a mixture of two or more amino acids at that residue (as bulk HIV-1 sequencing was performed).

https://doi.org/10.1371/journal.ppat.1009177.g004

To demonstrate that these significant increases in RT135L mutation frequencies are particular to Japan, and not mirrored in global regions where B*52:01/C*12:02 carriage is uncommon, we analyzed RT135 amino acid frequencies in individuals with HIV subtype B recruited from 2000 to 2015 in British Columbia (BC), Canada. Of 6,258 individuals analyzed, 56% and 27% harbored RT135I or RT135T, respectively, but only 1% of them harbored RT135L. Moreover, RT135 amino acid frequencies remained unchanged over three periods: 2000–2003, 2004–2009, and 2010–2015 (Fig 4B). In British Columbia, population frequency of HLA-B*51:01 is 11.7% whereas that of HLA-B*52:01 or -C*12:02 is only 1.4% [24], suggesting that RT135V selection by HLA-B*52:01-restricted TI8-specific T cells and elicitation of HLA-C*12:02-restricted TN9-8V-specific T cells would occur only rarely in this population.

Selection of RT135L by TN9-8V-specific HLA-C*12:02-restricted T cells

Our observation that TN9-8V-specific CD8⁺ T cells could not recognize target cells pulsed with TN9-8L peptide nor cells infected with RT135L mutant virus (Fig 3A and 3C) suggests that RT135L confers escape from these T cells. To investigate this, we analyzed HIV-1 evolution in a hemophiliac individual homozygous for the HLA-B*52:01⁺HLA-C*12:02⁺ haplotype (KI-528). Longitudinal bulk and deep sequencing of plasma HIV-1 RNA revealed that this individual acquired both wild-type (RT135I) and RT135V viruses in May 2003 (Fig 5A). By March 2006, RT135V had disappeared and RT135L had emerged as the dominant variant (68%). Moreover, within-host viral variants harboring RT135L comprised 2 groups having cta and tta at RT135, suggesting that both variants could have evolved from either RT135I or RT135V sequences via only one nucleotide change. By June 2007, RT135L had increased to 92% frequency, with both cta and tta variants in equal proportions.

Download:

Fig 5. Longitudinal analysis of RT135 mutant selection in an HLA-B*52:01⁺C*12:02⁺ individual.

A. Longitudinal sequence analysis at RT135 in an HLA-B*52:01⁺C*12:02⁺ hemophiliac individual with HIV-1, KI-528. Frequencies of amino acids at RT135 were determined by both bulk and deep sequencing. Both analyses showed similar results. The amino acids at RT135, their codon usage, and their frequency are shown at each sampling date. Estimated pathways of viral evolution during May 2003-June 2007 are indicated using arrows. B. Detection of TI8- and TN9-8V-specific CD8⁺ T cells in PBMCs from KI-528 at two different time points. These T cells were detected using TI8/B5201 and TN9-8V/C1202 tetramers. C. Recognition of TN9-8V and TN9-8L peptides by TN9-8V-specific bulk T cells from KI-528. T-cell responses to 721.221-CD4-C1202 cells prepulsed with TN9-8V or TN9-8L peptide were analyzed by ICS, and representative flow cytometric results are shown (left). The relative percent of IFN-γ⁺ cells among CD8⁺ cells was calculated by dividing their response frequency to each peptide by their response frequency to TN9-8V (right). Data are shown as the means and SD of triplicate assays. Statistical analysis was conducted using Student's t-test (two-tailed). ***P<0.001.

https://doi.org/10.1371/journal.ppat.1009177.g005

We next investigated the emergence of TI8-specific HLA-B*52:01-restricted and TN9-8V-specific HLA-C*12:02-restricted CD8⁺ T cells in this individual. Analysis using B*52:01-TI8 and C*12:02-TN9-8V tetramers confirmed the presence of TI8-specific and TN9-8V-specific CD8⁺ T cells in May 2004 (Figs 5B and S3B). The frequency of these T cells however decreased in April 2006, after RT135L had become the dominant variant. We stimulated PBMCs from this individual with TN9-8V peptide and cultured them to induce TN9-8V-specific T cells. The bulk T cells responded to TN9-8V peptide but not TN9-8L peptide (Fig 5C), suggesting that the TN9-8V-specific T cells can select RT135L mutant virus.

Suppression of RT135T mutant virus by TN9-8T-specific HLA-C*12:02-restricted T cells

The observation that TN9-8T has a strong affinity for HLA-C*12:02 raises the possibility that TN9-8T-specific T cells can target RT135T virus-infected cells. To investigate this, we analyzed an HLA-B*52:01⁺C*12:02⁺ individual, KI-638, who had acquired both RT135I and RT135T viruses in April 2009 (Fig 6A). By May 2013, RT135T and RT135I had markedly decreased whereas RT135V predominated (85%). Analysis using two tetramers confirmed TI8-specific T cells in PBMCs from this individual in July 2009, and TN9-8V-specific T cells had emerged by May 2013 (Figs 6B left and S3C). These results support the idea that C*12:02-restricted TN9-8V-specific T cells are elicited after the emergence of 8V mutant virus, which is selected by TI8-specific HLA-B*52:01-restricted CD8⁺ T cells.

Download:

Fig 6. Viral suppression of RT135T mutant virus by TN9-8T-specific CTL line.

A. Longitudinal sequence analysis at RT135 in an HLA-B*52:01⁺C*12:02⁺ individual with HIV-1, KI-638. Frequencies of amino acids at RT135 were determined by both bulk and deep sequencing. Both analyses showed similar results. The amino acids at RT135, their codon usage, and their frequencies are summarized by sampling date. B. Detection of TI8-, TN9-8V- and TN9-8T-specific CD8⁺ T cells in PBMCs from KI-638. These T cells were detected by TN9-8V/C1202 and TI8/B5201 tetramers or both TN9-8T/C1202 and TN9-8V/C1202 tetramers. C. Recognition of TN9-8T and its mutant peptides by TN9-8T-specific T cell line derived from KI-638. T-cell responses to 721.221-CD4-C1202 cells prepulsed with TN9-8T or its mutant peptides were analyzed by ICS. D. Binding of TN9-8V or TN9-8T tetramer to TN9-8T-specific T cell line. Data are shown as the mean and SD of triplicate assays E. Recognition of RT135 mutant virus-infected cells by a TN9-8T-specific T cell line. T-cell responses to 721.221-CD4-C1202 cells infected with NL43 (NL43-8I) or the mutant virus (NL43-8V, -8T or -8L) were analyzed by ICS. The frequencies of p24⁺ cells among 721.221-CD4-C1202 cells infected with NL43-8I, -8V, -8T, and -8L were 31.8%, 27.0%, 31.8%, and 28.5%, respectively. Representative flow cytometric results (left) and summarized results for the frequency of IFN-γ⁺ cells among CD8⁺ cells (right) is shown. F. Suppression of viral replication by TN9-8T-specfic CTL line, where E:T ratio was 1:1. The concentration of p24 Ag in the culture supernatants collected at Day 5 was measured by ELISA (left) and the percentage of suppression was calculated (right). G. Detection of TN9-8V- and TN9-8T-specific CD8⁺ T cells in PBMCs from 7 individuals infected with RT135T mutant virus. HLA-C*12:02-restricted TN9-8V and TN9-8T-specific CD8⁺ T cells were analyzed by using TN9-8T/C1202- or TN9-8V/C1202-tetramer. Statistical analysis was performed using a paired t-test (two-tailed). Data are shown as the means and SD of triplicate assays (C, D, E, and F).

https://doi.org/10.1371/journal.ppat.1009177.g006

Analysis using HLA-C*12:02-tetramers specific for the TN9-8V or TN9-8T epitopes revealed that, in July 2009, this individual had HLA-B*52:01-restricted TI8-specific T cells and HLA-C*12:02-restricted TN9-8T-specific T cells, but no HLA-C*12:02-restricted TN9-8V-specific ones (Figs 6B right and S3C). This suggests that TN9-8T-specific T cells were elicited by RT135T virus-infected cells, not RT135V virus-infected ones. In contrast, in May 2013, three types of HLA-C*12:02-restricted T cells − TN9-8T-specific, TN9-8V-specific, and TN9-8V/8T cross-reactive ones − were observed in this participant (Fig 6B). Together this suggests that TN9-8V-specific and TN9-8V/8T cross-reactive T cells were elicited after the emergence of RT135V virus.

To further characterize the TN9-8T-specific T cells found in July 2009, we established a TN9-8T-specific T cell line. This T cell line cross-recognized the TN9-8V peptide to the same level as TN9-8T peptide (Fig 6C). The binding affinity of the TN9-8T tetramer to the T cell line was much higher than that of the TN9-8V tetramer (Fig 6D). These results indicated that the T cell expressed a T cell receptor (TCR) that could recognize both epitopes, though its affinity for HLA-C*12:02-TN9-8T was higher than its affinity for HLA-C*12:02-TN9-8V. The TN9-8T-specific T cell line also strongly recognized target cells infected with RT135T or RT135V virus (Fig 6E), and strongly suppressed the replication of HIV harboring these mutations (Fig 6F). Together, these results support the idea that TN9-8T-specific T cells suppressed the replication of RT135T virus in this individual between July 2009 and May 2013.

Finally, to clarify whether TN9-8T-specific T cells are frequently elicited in HLA-B*52:01⁺C*12:02⁺ individuals harboring RT135T virus, we stained PBMCs from seven additional HLA-B*52:01⁺C*12:02⁺ individuals harboring this mutant virus using HLA-C*12:02/TN9-8T or HLA-C*12:02/TN9-8V tetramers. We detected TN9-8T-specific T cells in four of them (Figs 6G and S5), suggesting that TN9-8T-specific T cells are frequently elicited in these individuals. However, the frequency of HLA-C*12:02/TN9-8V tetramer⁺ T cells was significantly lower than that of HLA-C*12:02/TN9-8T tetramer⁺ ones (Figs 6G and S5). This indicates that TN9-8T-specific T cells are also preferentially elicited in these individuals.

Discussion

HIV-1 escape mutations reduce CTL recognition of HIV-infected cells by reducing or abrogating epitope processing, HLA binding and/or TCR recognition of the HLA-epitope complex [10,39,40]. Although many escape mutant epitopes cannot elicit specific T cells, some of them can to some extent [41–46]. For example, CTLs specific for two superimposed HLA-A*24:02-restricted epitopes, NefRF10 and NefRW8, select the Nef135F mutation [47]. Whereas the 10-mer epitope having this mutation (RF10-2F) further elicited HLA-A*24:02-restricted CTLs with high affinity TCR specific for RF10-2F, the 8-mer mutant epitope failed to elicit HLA-A*24:02-restricted T cells specific for RW8-2F [48]. We also speculated that immune escape could generate a new epitope recognized by another HLA class I-restricted CTL, which in turn could drive further escape, but such a phenomenon had not been reported to date. Here, we demonstrated that a new HLA-C*12:02-restricted epitope TN9-8V was generated by the RT135V mutation which was selected by HLA-B*52:01-restricted TI8-specific T cells in individuals carrying the HLA-B*52:01-C*12:02 haplotype. Moreover, due to strong linkage disequilibrium between these alleles and the high prevalence (20%) of this haplotype in Japan [34,36], this novel TN9-8V epitope was generated in the Japanese population far more so than in other global populations.

Japanese hemophiliacs acquired HIV-1 strains circulating in North America via contaminated blood products imported from the USA around 1983. To exclude the effect of ongoing HIV-1 domestic transmission on population-level accumulation of escape mutations, we began by studying these individuals. We demonstrated that RT135T and RT135V were significantly enriched among hemophiliacs carrying HLA-B*51:01⁺ and HLA-B*52:01⁺, respectively, suggesting that these mutations are selected by HLA-B*51:01- and HLA-B*52:01-restricted TI8-specific T cells, respectively, in this group. Previous studies demonstrated that neither HLA-B*51:01- nor HLA-B*52:01-restricted TI8-specific CTLs could recognize cells infected with HIV-1 carrying RT135T, RT135L, or RT135R mutations, whereas they could recognize those infected with RT135V mutant virus [2,30,33]. HLA-B*51:01-restricted TI8-specific T cells from individuals harboring wild-type (RT135I) virus suppressed replication of RT135V and wild-type viruses to a similar extent, whereas those from individuals harboring RT135V virus suppressed RT135V mutant virus more readily than wild-type virus in a competitive suppression assay [30]. In contrast, HLA-B*52:01-restricted TI8-specific T cells suppressed replication of wild-type virus more effectively than RT135V virus [33]. These findings together partially explain why HLA-B*52:01-restricted TI8-specific T cells select RT135V more often than HLA-B*51:01-restricted ones. It does not explain why in the Japanese hemophiliac group however, HLA-B*52:01-restricted T cells did not select other RT135 mutants such as RT135T, L or R as was observed in the non-hemophiliac group. Furthermore, although RT135V was enriched in HLA-B*52:01⁺ hemophiliacs, it was not enriched in HLA-B*52:01⁺ non-hemophiliacs. This may be explained by the source of transmitted HIV-1: whereas hemophiliacs acquired HIV-1 from North America around 1983, non-hemophiliacs predominantly acquired HIV-1 strains circulating in Japan. Since non-hemophiliac Japanese individuals diagnosed before 1997 predominantly harbored RT135T/L mutations, it is likely that these mutant viruses were initially transmitted and subsequently accumulated in domestically-acquired HIV cases, who were mostly men who have sex with men.

RT135 mutation frequencies have shifted dramatically over the past 20 years in non-hemophiliac Japanese individuals with HIV-1. Among those diagnosed before 1997, RT135T dominated (41%), but its frequency decreased to 25% among those diagnosed between 2010–2015. In contrast, RT135L dominated at this time (44%), whereas its frequency was only 22% among those diagnosed before 1997. HLA-B*51:01⁺ and HLA-B*51:01^-HLA-B*52:01^- non-hemophiliac individuals showed similar patterns of change in mutation frequencies over the past 20 years: a decrease in RT135T frequency and an increase in RT135L. This may be explained by viral transmission from HLA-B*52:01⁺C*12:02⁺ individuals who selected RT135L but suppressed replication of RT135T. The frequency of RT135T in HLA-B*52:01⁺C*12:02⁺ individuals was lower than among HLA-B*51:01⁺ ones, and comparable to that in HLA-B*51:01^-HLA-B*52:01^- ones who do not have HLA-B*52:01-restricted TI8-specific/HLA-C*12:02-restricted TN9-8V-specific T cells. This raises the possibility that RT135T is not directly selected by HLA-B*52:01-restricted TI8-specific T cells, or is eliminated by TN9-8T/8V-specific T cells in HLA-B*52:01⁺C*12:02⁺ individuals, as outlined in our proposed mechanistic model shown in Fig 7. On the other hand, in British Columbia, Canada, where the frequencies of HLA-B*51:01 and HLA-B*52:01 or HLA-C*12:02 were 11.7% and 1.4%, respectively, amino acid variation at RT135 remained extremely stable over a >15 year period with RT135T consistently representing the dominant mutation (26–28%) [49]. These findings further support the idea that HLA-B*52:01- and/or HLA-C*12:02-mediated immune responses drove the shifts in RT135 mutation frequencies over the last 20 years in Japan.

Download:

Fig 7. Proposed mechanism for coevolution of HIV-1 with HIV-1-specific T cells.

RT135V and RT135T viruses were selected by TI8-specific HLA-B*52:01-restricted and TI8-specific HLA-B*51:01-restricted T cells, respectively. In the case of selection of RT135V virus in HLA-B*52:01⁺HLA-C*12:02⁺ persons or transmission of this mutant virus to them, HLA-C*12:02-restricted TN9-8V-specific T cells were subsequently elicited, which led to the selection of RT135L virus rather than RT135T, since these T cells can suppress replication of RT135T virus but not RT135L. In the case of RT135T transmission to HLA-B*52:01⁺HLA-C*12:02⁺ individuals, HLA-C*12:02-restricted TN9-8T-specific T cells were subsequently elicited, which effectively suppressed replication of this mutant virus. Since RT135L virus cannot be recognized by HLA-C*12:02-restricted TN9-8T-specific T cells, it can be selected in HLA-B*52:01⁺HLA-C*12:02⁺ individuals infected with RT135T virus.

https://doi.org/10.1371/journal.ppat.1009177.g007

The frequency of HLA-B*51:01-restricted TI8-specific T cells correlated inversely with pVL whereas that of HLA-B*51:01-restricted T cells specific for three other epitopes did not [29], suggesting that TI8-specific T cells contribute to HLA-B*51:01-restricted HIV-1 control. Consistent with this, HLA-B*51:01 was previously a protective allele in Japanese hemophiliacs [29] and other populations [26,50]. However, this allele is no longer protective in non-hemophiliac Japanese individuals [31,32]. Since 80% of HLA-B*51:01⁺ Japanese individuals harbor RT135T/L, it is unlikely that HLA-B*51:01⁺ TI8-specific T cells can control replication of RT135T/L mutant virus in these persons. Thus, accumulation of RT135T/L at the population level critically compromised HLA-B*51:01-mediated control of HIV-1. Instead, HLA-B*52:01-C*12:02 represents a protective haplotype in non-hemophiliac Japanese individuals [32]. Consistent with this, a genome-wide association study of HIV-1 control in White (Caucasian) and Black individuals identified HLA-B*52:01, but not HLA-B*51:01, as a protective allele together with HLA-B*57 and HLA-B*27 [51]. Previous studies described four HLA-B*52:01-restricted and two HLA-C*12:02-restricted protective epitopes [35,52] but TI8 was not among these. In contrast, our recent study of HLA-B*52:01⁺ TI8 responders diagnosed from 2011 to 2017 in Japan reported good clinical outcome in these individuals [53], suggesting that HLA-B*52:01-restricted TI8-specific T cells may have the ability to control HIV-1 in individuals who acquired the virus in more recent years. The role of TN9-8V-specific T cells in HIV-1 control in vivo remains unclear, and merits further study.

In summary, we demonstrated dynamic shifts in HIV-1 RT135 mutation frequency over the past two decades in Japan, where RT135T frequency significantly decreased among non-hemophiliacs with HIV-1, and RT135L significantly increased. We further describe a plausible mechanism for these shifts, summarized in Fig 7. Specifically, we describe novel T cell responses restricted by the most common HLA-B/HLA-C allele haplotype expressed in Japan, and demonstrate that these likely contributed to dynamic population-level changes in RT135 mutation frequencies in Japan that are not mirrored in other parts of the globe. Our observations underscore HIV-1's dynamic ability to co-evolve with virus-specific T cells in a host population-specific manner.

Methods

Ethics statement

This study was approved by the ethics committee of Kumamoto University (RINRI-1340 and GENOME-342) and the National Center for Global Health and Medicine (NCGM-A-000172-01). Written informed consent was obtained from all individuals according to the Declaration of Helsinki.

Participants

Plasma and peripheral blood mononuclear cells (PBMCs) were separated from whole blood. HLA genotypes were determined by the HLA Foundation Laboratory, Kyoto, Japan. Plasma and PBMCs from 23 HLA-B*52:01⁺HLA-C*12:02⁺ individuals with chronic HIV-1 subtype B were characterized with respect to HIV-1 sequence and CTL response. HIV-1 sequences from 95 Japanese hemophiliac patients, who acquired HIV-1 subtype B via contaminated blood products imported from the USA around 1983, were analyzed for HIV-1 amino acid variation at RT135. For comparison, amino acid frequencies at RT135 derived from 1,478 Japanese individuals with subtype B HIV-1 enrolled in NCGM from before 1997 to 2015, and 6,258 British Columbians with subtype B HIV-1 collected between 2000 and 2015 [49] were also analyzed.

Cells

C1R cells expressing HLA-B*52:01 (C1R-B5201) and HLA-C*12:02 (C1R-C1202), and 721.221 cells expressing CD4 and HLA-C*12:02 (721.221-CD4-C1202) were previously generated [33,54]. These cells were cultured in RPMI-1640 medium (Invitrogen) containing with 5–10% fetal bovine serum (R5-R10) and 0.15 mg/ml hygromycin B. RMA-S-C1202 cells, a TAP deficient mouse-derived cell line that expresses HLA-C*12:02, were previously generated and cultured in R10 containing 0.15 mg/ml hygromycin B [55]. MAGIC-5 cells (CCR5-transfected HeLa-CD4/LTR-β-gal cells) were cultured in R10 as described previously [56].

Induction of peptide-specific bulk T cells

PBMCs from HLA-B*52:01⁺ and HLA-C*12:02⁺ expressing individual with chronic HIV-1, KI-793, were stimulated with KN11 or KN11-10V peptide (1 μM) in R10 supplemented with 200 U/ml human recombinant interleukin-2 (rIL-2) to induce peptide-specific bulk T cells for 14 days. After induction, bulk T cells were examined for IFN-γ production using an intracellular cytokine staining (ICS) assay.

ICS assay

C1R-C1202 or 721.221-CD4-C1202 cells were pre-incubated with or without their cognate peptide(s) at concentrations from 0.1 to 10 nM at 37°C for 1 hr, and washed twice with R10. Cultured bulk T cells or established CTL lines were incubated with peptide-pulsed or virus-infected 721.221-CD4-C1202 cells in R10 containing Brefeldin A (10 μg/ml) in a 96-U plate (Nunc) at 37°C for 4 hrs. Subsequently, the cells were stained with allophycocyanin (APC)-conjugated anti-CD8 monoclonal antibody (mAb) (Dako) and 7-amino-actinomycin D (7-AAD) (BD Bioscience) at 4°C for 30 min, after which the cells were fixed with 4% paraformaldehyde solution and permeabilized with permeabilization buffer (0.1% saponin and 5% FBS in phosphate-buffered saline) at 4°C for 10 min. Thereafter the cells were stained with fluorescein isothiocyanate (FITC)-conjugated anti-IFN-γ mAb (BD Biosciences) at room temperature for 30 min and then washed twice with the permeabilization buffer. The percentage of CD8⁺ cells producing IFN-γ was analyzed by flow cytometry (FACS Canto II).

Generation of TN9-8V- and TN9-8T-specific CTL lines

TN9-8V-specific CTL lines were generated from PBMCs of three chronically HIV-1-infected individuals, KI-069, -661 and -793, while TN9-8T-specific CTL lines were generated from PBMCs from chronically HIV-1-infected individual KI-638. PBMCs were stained with TN9-8V/C1202- or TN9-8T/C1202-tetramer at 37°C for 30 min and washed twice with R10, after which they were stained with FITC-conjugated anti-CD3, Pacific blue-conjugated CD8 and 7-Aminoactinomycin D (7AAD) viability dye at 4°C for 30 min and washed twice with R10. After staining, TN9-8V- or TN9-8T-specific CD8⁺ T cells were sorted by FACSAria into a 96-well U-bottom plate and cultured with irradiated feeder and stimulator cells in cloning medium containing 100 nM of TN9-8V or TN9-8T peptide. Half of the medium was replaced every 2–3 days with fresh medium without peptide. After 14–21 days in culture, the growing cells were tested for their ability to recognize TN9-8V or TN9-8T peptide by ICS assay.

HIV-1 clones

NL43-RT135X mutant viruses were previously generated by site-directed mutagenesis [2].

Generation of HLA-peptide tetrameric complexes

HLA class I-peptide tetrameric complexes (tetramers) were synthesized as previously described [57]. TI8/B5201 monomer complex was previously generated [33]. TN9-8V or TN9-8T peptide was added to the refolding solution containing the biotinylation-sequence-tagged extracellular domain of HLA-C*12:02 molecule and β2 microglobulin. The purified monomer complexes were mixed with Phycoerythrin (PE)- or APC labeled streptavidin (Molecular Probes) at a molar ratio of 4:1.

Tetramer binding assay

PBMCs and the TN9-8T-specific CTL line were stained with PE or APC-conjugated tetramer at 100 nM and 1–100 nM, respectively at 37°C for 30 min. After washing twice with R10, the cells were stained with FITC-conjugated anti-CD3 mAb, Pacific blue-conjugated anti-CD8 mAb and 7-AAD at 4°C for 30 min. Thereafter, the cells were washed twice with R10 and then analyzed by flow cytometry (FACS Canto II).

Viral replication suppression assay

The ability of TN9-8V-specific and TN9-8T-specific CTL lines to suppress HIV-1 replication was examined as previously described [58]. Activated CD4⁺ T cells (1x10⁶ or 5x10⁵ cells) were isolated from PBMCs of healthy HLA-B*52:01⁺ C*12:02⁺ donors and incubated with a given HIV-1 clone at 37°C for 6 hrs. After washing twice with R10, the cells (1x10⁴/well) were co-cultured with TN9-8V- or TN9-8T-specific CTL lines at E:T ratios of 1:1 in culture media containing IL2. From day 3 to day 7 post-infection, 30 μl of culture supernatant was collected and replaced with fresh medium. The concentration of p24 Ag was measured using an enzyme-linked immunosorbent assay (ELISA) (ZeptoMetrix or Rimco). Percent suppression was calculated as follows: (concentration of p24 Ag without the CTLs—concentration of p24 Ag with the CTLs)/concentration of p24 Ag without the CTLs x 100.

Sequence of autologous virus

For longitudinal HIV-1 sequence analysis, data were obtained as follows. Viral RNA was extracted from plasma using a QIAamp MinElute virus spin kit (QIAGEN). cDNA was synthesized from the RNA with SuperScript III first-strand kit (Invitrogen) and then the Pol region was amplified by nested PCR with Taq DNA polymerase (Promega). DNA sequencing was performed using a BigDye Terminator cycle sequencing kit (Applied Biosystems) on an ABI PRISM genetic analyzer.

To determine the frequencies of minority RT135 variants in studied samples, deep sequencing was performed using Miseq (Illumina). PCR products were purified by PCR-cleanup kit (SIGMA). Tagmentation and normalization were performed using the Nextera XT sample preparation kit (Illumina). Normalized samples were pooled and loaded into Miseq reagent kit v3 600 cycles and sequenced following the manufacturer’s instructions. Miseq reads were quality trimmed using Sickle (https://github.com/najoshi/sickle). The reads with high-quality scores (Q>30) were mapped to the HIV-1 NL4-3 reference genome (GenBank: AF324493) using Bowtie2. After sorting and converting using samtools, mapped reads were visualized using the Integrative Genomics Viewer (IGV) after which amino acid frequencies were calculated.

HLA class I stabilization assay

Peptide-binding activity to HLA-C*12:02 was examined using RMA-S-C1202 cells as previously described [59]. Briefly, RMA-S-C1202 cells were cultured at 26°C for 16 h and then pulsed with peptides at concentrations of 1 to 100 μM at 26°C for 1 h. Subsequently, the cells were incubated at 37°C for 3 h, followed by staining with anti-HLA-Bw6 antibody (SFR8-B6) which cross-binds to HLA-C [60] and FITC-conjugated sheep anti-mouse IgG (Jackson Immuno Research). The MFI was measured by flow cytometry (FACS Canto II). Relative MFI was calculated by subtracting the MFI of peptide-unpulsed cells from that of the peptide-pulsed ones.

Statistical analysis

The frequency of the amino acids between HLA⁺ and HLA^- individuals was analyzed by using Fisher’s exact test. Multiple tests were addressed using q values, the p-value analogue of the false discovery rate. A significance threshold of q <0.2 was employed. The differences in the peptide binding affinities (Fig 3B) or in the T-cell responses (Fig 5C) between TN9-8V and other peptides were statistically analyzed by performing the two-tailed unpaired t test. The difference between tetramer⁺ T-cell frequencies was compared using the paired t test (Fig 6G). All tests were two-tailed and p values of <0.05 were considered significant.

The association between V, T, R, M, L, K, I frequencies across the four periods was assessed using the Cochran–Mantel–Haenszel test, where the mutations were treated as nominal variables and the four periods were treated as ordered categorical data (Fig 4A). Analyses were performed overall, and stratified by allele carriage. Statistical analyses were performed using SAS (SAS Institute), version 9.4.

Supporting information

S1 Fig. Nucleotide substitution at RT135 in Japanese individuals with chronic HIV-1, Related to Fig 1.

Codon usages for I, V, and T are shown, along with their frequencies observed in Japanese individuals with chronic HIV-1 (Hemo; 95 hemophiliacs shown in Fig 1A, Non-hemo; 83 non-hemophiliac individuals shown in Fig 1B and S1 Table).

https://doi.org/10.1371/journal.ppat.1009177.s001

(TIF)

S2 Fig. Identification of a novel HLA-C*12:02-restricted CTL epitope including RT135, Related to Fig 2A.

A. T-cell responses to 11-mer overlapping Pol peptides containing RT135. T-cell responses of PBMCs from an HLA-B*52:01⁺C*12:02⁺ individual (KI-793) to five 11-mer overlapping Pol peptides containing RT135 position and Pol peptide cocktail 15 including the 5 overlapping peptides were analyzed at a concentration of 100 nM by ELISPOT assay. B. Identification of HLA-restriction of the response to the Pol 11–141 (KN11) peptide. IFN-γ production from KN11-induced bulk T cells stimulated with C1R cells expressing either HLA-B*52:01 or -C*12:02 molecule pre-pulsed with the KN11 peptide at concentration of 100 nM was analyzed by ICS assay. C. Comparison of induction efficiency between bulk T cells induced with the KN11 peptide and those with KN11-10V mutant one. IFN-γ production from KN11-induced or KN11-10V-induced bulk T cells stimulated with C1R-C1202 cells pre-pulsed with the KN11 or KN11-10V peptide was analyzed by ICS assay. Relative % of IFN-γ⁺ cells among CD8⁺ T cells was calculated as follows: % of IFN-γ⁺ cells with peptide–% of IFN-γ⁺ cells without peptide.

https://doi.org/10.1371/journal.ppat.1009177.s002

(TIF)

S3 Fig. Negative controls of ex vivo tetramer staining assay, Related to Figs 2B, 5B and 6B.

A. Staining of KI-793 PBMCs without tetramer, related to Fig 2B. B. Staining of KI-528 PBMCs without tetramers, related to Fig 5B. C. Staining of KI-638 PBMCs without tetramers (upper) and staining of PBMCs derived from a healthy donor with tetramers (bottom), related to Fig 6B.

https://doi.org/10.1371/journal.ppat.1009177.s003

(TIF)

S4 Fig. Recognition of virus-infected cells by TN9-8V-specific T cell lines, Related to Fig 3C.

Responses by TN9-8V-specific CTL lines were established from 3 HLA-B*52:01⁺C*12:02⁺ individuals with chronic HIV-1. The ability of these T cells to recognize 721.221-CD4-C1202 cells infected with NL43, or NL43-RT135X mutant viruses was analyzed by ICS assay. Frequency of IFN-γ⁺ cells among CD8⁺ cells was indicated.

https://doi.org/10.1371/journal.ppat.1009177.s004

(TIF)

S5 Fig. Ex vivo detection of TN9-8T-specific CD8⁺ T cells, Related to Fig 6G.

PBMCs from seven HLA-B*52:01⁺ C*12:02⁺ individuals harboring HIV-1 RT135T virus were stained with TN9-8V/C1202 tetramer or TN9-8T one at concentration of 100 nM. Representative cases corresponding to Fig 6G were shown. Frequency of tetramer⁺ cells among CD3⁺CD8⁺ T cells is indicated in red.

https://doi.org/10.1371/journal.ppat.1009177.s005

(TIF)

S1 Table. Codon usages of amino acids at RT135 observed in Japanese hemophiliacs and non-hemophiliacs with chronic HIV-1

https://doi.org/10.1371/journal.ppat.1009177.s006

(DOCX)

S2 Table. Statistical analysis using Cochran-Mantel-Haenszel test on the association between the frequency of 7 amino acids and the four periods.

https://doi.org/10.1371/journal.ppat.1009177.s007

(DOCX)

Acknowledgments

The authors thank Jeffrey Joy, Julio Montaner and the BC Centre for Excellence in HIV/AIDS laboratory for data access and support.

References

1. Allen TM, Altfeld M, Geer SC, Kalife ET, Moore C, O'Sullivan K M, et al. Selective escape from CD8+ T-cell responses represents a major driving force of human immunodeficiency virus type 1 (HIV-1) sequence diversity and reveals constraints on HIV-1 evolution. J Virol. 2005;79(21):13239–49. pmid:16227247
- View Article
- PubMed/NCBI
- Google Scholar
2. Kawashima Y, Pfafferott K, Frater J, Matthews P, Payne R, Addo M, et al. Adaptation of HIV-1 to human leukocyte antigen class I. Nature. 2009;458(7238):641–5. pmid:19242411
- View Article
- PubMed/NCBI
- Google Scholar
3. Leslie A, Kavanagh D, Honeyborne I, Pfafferott K, Edwards C, Pillay T, et al. Transmission and accumulation of CTL escape variants drive negative associations between HIV polymorphisms and HLA. J Exp Med. 2005;201(6):891–902. pmid:15781581
- View Article
- PubMed/NCBI
- Google Scholar
4. Liu Y, McNevin J, Cao J, Zhao H, Genowati I, Wong K, et al. Selection on the human immunodeficiency virus type 1 proteome following primary infection. J Virol. 2006;80(19):9519–29. pmid:16973556
- View Article
- PubMed/NCBI
- Google Scholar
5. Moore CB, John M, James IR, Christiansen FT, Witt CS, and Mallal SA. Evidence of HIV-1 adaptation to HLA-restricted immune responses at a population level. Science. 2002;296(5572):1439–43. pmid:12029127
- View Article
- PubMed/NCBI
- Google Scholar
6. Crawford H, Lumm W, Leslie A, Schaefer M, Boeras D, Prado JG, et al. Evolution of HLA-B*5703 HIV-1 escape mutations in HLA-B*5703-positive individuals and their transmission recipients. J Exp Med. 2009;206(4):909–21. pmid:19307327
- View Article
- PubMed/NCBI
- Google Scholar
7. Friedrich TC, Dodds EJ, Yant LJ, Vojnov L, Rudersdorf R, Cullen C, et al. Reversion of CTL escape-variant immunodeficiency viruses in vivo. Nat Med. 2004;10(3):275–81. pmid:14966520
- View Article
- PubMed/NCBI
- Google Scholar
8. Leslie AJ, Pfafferott KJ, Chetty P, Draenert R, Addo MM, Feeney M, et al. HIV evolution: CTL escape mutation and reversion after transmission. Nat Med. 2004;10(3):282–9. pmid:14770175
- View Article
- PubMed/NCBI
- Google Scholar
9. Li B, Gladden AD, Altfeld M, Kaldor JM, Cooper DA, Kelleher AD, et al. Rapid reversion of sequence polymorphisms dominates early human immunodeficiency virus type 1 evolution. J Virol. 2007;81(1):193–201. pmid:17065207
- View Article
- PubMed/NCBI
- Google Scholar
10. Goulder PJ, Brander C, Tang Y, Tremblay C, Colbert RA, Addo MM, et al. Evolution and transmission of stable CTL escape mutations in HIV infection. Nature. 2001;412(6844):334–8. pmid:11460164
- View Article
- PubMed/NCBI
- Google Scholar
11. Fryer HR, Frater J, Duda A, Roberts MG, Investigators ST, Phillips RE, et al. Modelling the evolution and spread of HIV immune escape mutants. PLoS Pathog. 2010;6(11):e1001196. pmid:21124991
- View Article
- PubMed/NCBI
- Google Scholar
12. Navis M, Matas DE, Rachinger A, Koning FA, van Swieten P, Kootstra NA, et al. Molecular evolution of human immunodeficiency virus type 1 upon transmission between human leukocyte antigen disparate donor-recipient pairs. PLoS One. 2008;3(6):e2422. pmid:18560583
- View Article
- PubMed/NCBI
- Google Scholar
13. Avila-Rios S, Ormsby CE, Carlson JM, Valenzuela-Ponce H, Blanco-Heredia J, Garrido-Rodriguez D, et al. Unique features of HLA-mediated HIV evolution in a Mexican cohort: a comparative study. Retrovirology. 2009;6:72. pmid:19664284
- View Article
- PubMed/NCBI
- Google Scholar
14. Brumme ZL, Brumme CJ, Heckerman D, Korber BT, Daniels M, Carlson J, et al. Evidence of differential HLA class I-mediated viral evolution in functional and accessory/regulatory genes of HIV-1. PLoS Pathog. 2007;3(7):e94. pmid:17616974
- View Article
- PubMed/NCBI
- Google Scholar
15. Brumme ZL, John M, Carlson JM, Brumme CJ, Chan D, Brockman MA, et al. HLA-associated immune escape pathways in HIV-1 subtype B Gag, Pol and Nef proteins. PLoS One. 2009;4(8):e6687. pmid:19690614
- View Article
- PubMed/NCBI
- Google Scholar
16. Rousseau CM, Daniels MG, Carlson JM, Kadie C, Crawford H, Prendergast A, et al. HLA class I-driven evolution of human immunodeficiency virus type 1 subtype c proteome: immune escape and viral load. J Virol. 2008;82(13):6434–46. pmid:18434400
- View Article
- PubMed/NCBI
- Google Scholar
17. Carlson JM, Le AQ, Shahid A, and Brumme ZL. HIV-1 adaptation to HLA: a window into virus-host immune interactions. Trends Microbiol. 2015;23(4):212–24. pmid:25613992
- View Article
- PubMed/NCBI
- Google Scholar
18. Cotton LA, Kuang XT, Le AQ, Carlson JM, Chan B, Chopera DR, et al. Genotypic and functional impact of HIV-1 adaptation to its host population during the North American epidemic. PLoS Genet. 2014;10(4):e1004295. pmid:24762668
- View Article
- PubMed/NCBI
- Google Scholar
19. Kinloch NN, MacMillan DR, Le AQ, Cotton LA, Bangsberg DR, Buchbinder S, et al. Population-Level Immune-Mediated Adaptation in HIV-1 Polymerase during the North American Epidemic. J Virol. 2016;90(3):1244–58. pmid:26559841
- View Article
- PubMed/NCBI
- Google Scholar
20. Martin E, Carlson JM, Le AQ, Chopera DR, McGovern R, Rahman MA, et al. Early immune adaptation in HIV-1 revealed by population-level approaches. Retrovirology. 2014;11:64. pmid:25212686
- View Article
- PubMed/NCBI
- Google Scholar
21. Murakoshi H, Koyanagi M, Akahoshi T, Chikata T, Kuse N, Gatanaga H, et al. Impact of a single HLA-A*24:02-associated escape mutation on the detrimental effect of HLA-B*35:01 in HIV-1 control. EBioMedicine. 2018;36:103–12. pmid:30249546
- View Article
- PubMed/NCBI
- Google Scholar
22. Monaco DC, Dilernia DA, Fiore-Gartland A, Yu T, Prince JL, Dennis KK, et al. Balance between transmitted HLA preadapted and nonassociated polymorphisms is a major determinant of HIV-1 disease progression. J Exp Med. 2016;213(10):2049–63. pmid:27551154
- View Article
- PubMed/NCBI
- Google Scholar
23. Carlson JM, Du VY, Pfeifer N, Bansal A, Tan VY, Power K, et al. Impact of pre-adapted HIV transmission. Nat Med. 2016;22(6):606–13. pmid:27183217
- View Article
- PubMed/NCBI
- Google Scholar
24. Carlson JM, Brumme CJ, Martin E, Listgarten J, Brockman MA, Le AQ, et al. Correlates of protective cellular immunity revealed by analysis of population-level immune escape pathways in HIV-1. J Virol. 2012;86(24):13202–16. pmid:23055555
- View Article
- PubMed/NCBI
- Google Scholar
25. Murakoshi H, Koyanagi M, Chikata T, Rahman MA, Kuse N, Sakai K, et al. 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. J Virol. 2017;91(4). pmid:27903797
- View Article
- PubMed/NCBI
- Google Scholar
26. Zhang Y, Peng Y, Yan H, Xu K, Saito M, Wu H, et al. Multilayered defense in HLA-B51-associated HIV viral control. J Immunol. 2011;187(2):684–91. pmid:21670313
- View Article
- PubMed/NCBI
- Google Scholar
27. Goulder PJ, Phillips RE, Colbert RA, McAdam S, Ogg G, Nowak MA, et al. Late escape from an immunodominant cytotoxic T-lymphocyte response associated with progression to AIDS. Nat Med. 1997;3(2):212–7. pmid:9018241
- View Article
- PubMed/NCBI
- Google Scholar
28. Kelleher AD, Long C, Holmes EC, Allen RL, Wilson J, Conlon C, et al. Clustered mutations in HIV-1 gag are consistently required for escape from HLA-B27-restricted cytotoxic T lymphocyte responses. J Exp Med. 2001;193(3):375–86. pmid:11157057
- View Article
- PubMed/NCBI
- Google Scholar
29. Kawashima Y, Kuse N, Gatanaga H, Naruto T, Fujiwara M, Dohki S, et al. Long-term control of HIV-1 in hemophiliacs carrying slow-progressing allele HLA-B*5101. J Virol. 2010;84(14):7151–60. pmid:20410273
- View Article
- PubMed/NCBI
- Google Scholar
30. Kuse N, Akahoshi T, Gatanaga H, Ueno T, Oka S, and Takiguchi M. Selection of TI8-8V mutant associated with long-term control of HIV-1 by cross-reactive HLA-B*51:01-restricted cytotoxic T cells. J Immunol. 2014;193(10):4814–22. pmid:25305317
- View Article
- PubMed/NCBI
- Google Scholar
31. Koga M, Kawana-Tachikawa A, Heckerman D, Odawara T, Nakamura H, Koibuchi T, et al. Changes in impact of HLA class I allele expression on HIV-1 plasma virus loads at a population level over time. Microbiol Immunol. 2010;54(4):196–205. pmid:20377748
- View Article
- PubMed/NCBI
- Google Scholar
32. Naruto T, Gatanaga H, Nelson G, Sakai K, Carrington M, Oka S, et al. HLA class I-mediated control of HIV-1 in the Japanese population, in which the protective HLA-B*57 and HLA-B*27 alleles are absent. J Virol. 2012;86(19):10870–2. pmid:22811530
- View Article
- PubMed/NCBI
- Google Scholar
33. Yagita Y, Kuse N, Kuroki K, Gatanaga H, Carlson JM, Chikata T, et al. Distinct HIV-1 escape patterns selected by cytotoxic T cells with identical epitope specificity. J Virol. 2013;87(4):2253–63. pmid:23236061
- View Article
- PubMed/NCBI
- Google Scholar
34. Ikeda N, Kojima H, Nishikawa M, Hayashi K, Futagami T, Tsujino T, et al. Determination of HLA-A, -C, -B, -DRB1 allele and haplotype frequency in Japanese population based on family study. Tissue Antigens. 2015;85(4):252–9. pmid:25789826
- View Article
- PubMed/NCBI
- Google Scholar
35. Murakoshi H, Akahoshi T, Koyanagi M, Chikata T, Naruto T, Maruyama R, et al. Clinical Control of HIV-1 by Cytotoxic T Cells Specific for Multiple Conserved Epitopes. J Virol. 2015;89(10):5330–9. pmid:25741000
- View Article
- PubMed/NCBI
- Google Scholar
36. Saito S, Ota S, Yamada E, Inoko H, and Ota M. Allele frequencies and haplotypic associations defined by allelic DNA typing at HLA class I and class II loci in the Japanese population. Tissue antigens. 2000;56(6):522–9. pmid:11169242
- View Article
- PubMed/NCBI
- Google Scholar
37. Iizuka A, Nagao T, Shohji H, Nitta A, Yoshihara N, and Matsuhashi T. A longitudinal study on seroconversion of HTLV-III/LAV in Japanese haemophiliacs. Br J Haematol. 1987;65(2):249–50. pmid:3644665
- View Article
- PubMed/NCBI
- Google Scholar
38. Kashiwagi S, Hayashi J, Ikematsu H, Noguchi A, Ikeda K, Kishida K, et al. Prevalence of human immunodeficiency virus (HIV) infection among hemophiliacs in Fukuoka, Japan. Microbiol Immunol. 1988;32(8):851–5. pmid:3143898
- View Article
- PubMed/NCBI
- Google Scholar
39. Altman JD, and Feinberg MB. HIV escape: there and back again. Nat Med. 2004;10(3):229–30. pmid:14991039
- View Article
- PubMed/NCBI
- Google Scholar
40. Bailey J, Blankson JN, Wind-Rotolo M, and Siliciano RF. Mechanisms of HIV-1 escape from immune responses and antiretroviral drugs. Curr Opin Immunol. 2004;16(4):470–6. pmid:15245741
- View Article
- PubMed/NCBI
- Google Scholar
41. Allen TM, Yu XG, Kalife ET, Reyor LL, Lichterfeld M, John M, et al. De novo generation of escape variant-specific CD8+ T-cell responses following cytotoxic T-lymphocyte escape in chronic human immunodeficiency virus type 1 infection. J Virol. 2005;79(20):12952–60. pmid:16188997
- View Article
- PubMed/NCBI
- Google Scholar
42. Feeney ME, Tang Y, Pfafferott K, Roosevelt KA, Draenert R, Trocha A, et al. HIV-1 viral escape in infancy followed by emergence of a variant-specific CTL response. J Immunol. 2005;174(12):7524–30. pmid:15944251
- View Article
- PubMed/NCBI
- Google Scholar
43. Melhem NM, Smith KN, Huang XL, Colleton BA, Jiang W, Mailliard RB, et al. The impact of viral evolution and frequency of variant epitopes on primary and memory human immunodeficiency virus type 1-specific CD8(+) T cell responses. Virology. 2014;450–451:34–48. pmid:24503065
- View Article
- PubMed/NCBI
- Google Scholar
44. Roider J, Kalteis AL, Vollbrecht T, Gloning L, Stirner R, Henrich N, et al. Adaptation of CD8 T cell responses to changing HIV-1 sequences in a cohort of HIV-1 infected individuals not selected for a certain HLA allele. PLoS One. 2013;8(12):e80045. pmid:24312453
- View Article
- PubMed/NCBI
- Google Scholar
45. Bailey JR, Williams TM, Siliciano RF, and Blankson JN. Maintenance of viral suppression in HIV-1-infected HLA-B*57+ elite suppressors despite CTL escape mutations. J Exp Med. 2006;203(5):1357–69. pmid:16682496
- View Article
- PubMed/NCBI
- Google Scholar
46. Qin K, Boppana S, Du VY, Carlson JM, Yue L, Dilernia DA, et al. CD8 T cells targeting adapted epitopes in chronic HIV infection promote dendritic cell maturation and CD4 T cell trans-infection. PLoS Pathog. 2019;15(8):e1007970. pmid:31398241
- View Article
- PubMed/NCBI
- Google Scholar
47. Sun X, Fujiwara M, Shi Y, Kuse N, Gatanaga H, Appay V, et al. Superimposed epitopes restricted by the same HLA molecule drive distinct HIV-specific CD8+ T cell repertoires. J Immunol. 2014;193(1):77–84. pmid:24899498
- View Article
- PubMed/NCBI
- Google Scholar
48. Sun X, Shi Y, Akahoshi T, Fujiwara M, Gatanaga H, Schonbach C, et al. Effects of a Single Escape Mutation on T Cell and HIV-1 Co-adaptation. Cell Rep. 2016;15(10):2279–91. pmid:27239036
- View Article
- PubMed/NCBI
- Google Scholar
49. Brumme ZL, Kinloch NN, Sanche S, Wong A, Martin E, Cobarrubias KD, et al. Extensive host immune adaptation in a concentrated North American HIV epidemic. AIDS. 2018;32(14):1927–38. pmid:30048246
- View Article
- PubMed/NCBI
- Google Scholar
50. O'Brien SJ, Gao X, and Carrington M. HLA and AIDS: a cautionary tale. Trends Mol Med. 2001;7(9):379–81. pmid:11530315
- View Article
- PubMed/NCBI
- Google Scholar
51. International HIVCS, Pereyra F, Jia X, McLaren PJ, Telenti A, de Bakker PI, et al. The major genetic determinants of HIV-1 control affect HLA class I peptide presentation. Science. 2010;330(6010):1551–7. pmid:21051598
- View Article
- PubMed/NCBI
- Google Scholar
52. Chikata T, Murakoshi H, Koyanagi M, Honda K, Gatanaga H, Oka S, et al. Control of HIV-1 by an HLA-B*52:01-C*12:02 Protective Haplotype. J Infect Dis. 2017;216(11):1415–24. pmid:28968792
- View Article
- PubMed/NCBI
- Google Scholar
53. Zou C, Murakoshi H, Kuse N, Akahoshi T, Chikata T, Gatanaga H, et al. Effective Suppression of HIV-1 Replication by Cytotoxic T Lymphocytes Specific for Pol Epitopes in Conserved Mosaic Vaccine Immunogens. J Virol. 2019;93(7). pmid:30674626
- View Article
- PubMed/NCBI
- Google Scholar
54. Honda K, Zheng N, Murakoshi H, Hashimoto M, Sakai K, Borghan MA, et al. Selection of escape mutant by HLA-C-restricted HIV-1 Pol-specific cytotoxic T lymphocytes carrying strong ability to suppress HIV-1 replication. Eur J Immunol. 2011;41(1):97–106. pmid:21182081
- View Article
- PubMed/NCBI
- Google Scholar
55. Kuse N, Rahman MA, Murakoshi H, Tran GV, Chikata T, Koyanagi M, et al. Different Effects of Nonnucleoside Reverse Transcriptase Inhibitor Resistance Mutations on Cytotoxic T Lymphocyte Recognition between HIV-1 Subtype B and Subtype A/E Infections. J Virol. 2015;89(14):7363–72. pmid:25972553
- View Article
- PubMed/NCBI
- Google Scholar
56. Gatanaga H, Hachiya A, Kimura S, and Oka S. Mutations other than 103N in human immunodeficiency virus type 1 reverse transcriptase (RT) emerge from K103R polymorphism under non-nucleoside RT inhibitor pressure. Virology. 2006;344(2):354–62. pmid:16219331
- View Article
- PubMed/NCBI
- Google Scholar
57. Altman JD, Moss PA, Goulder PJ, Barouch DH, McHeyzer-Williams MG, Bell JI, et al. Phenotypic analysis of antigen-specific T lymphocytes. Science. 1996;274(5284):94–6. pmid:8810254
- View Article
- PubMed/NCBI
- Google Scholar
58. Tomiyama H, Akari H, Adachi A, and Takiguchi M. Different effects of Nef-mediated HLA class I down-regulation on human immunodeficiency virus type 1-specific CD8(+) T-cell cytolytic activity and cytokine production. J Virol. 2002;76(15):7535–43. pmid:12097566
- View Article
- PubMed/NCBI
- Google Scholar
59. Lin Z, Kuroki K, Kuse N, Sun X, Akahoshi T, Qi Y, et al. HIV-1 Control by NK Cells via Reduced Interaction between KIR2DL2 and HLA-C (*)12:02/C (*)14:03. Cell Rep. 2016;17(9):2210–20. pmid:27880898
- View Article
- PubMed/NCBI
- Google Scholar
60. Radka SF, Kostyu DD, and Amos DB. A monoclonal antibody directed against the HLA-Bw6 epitope. J Immunol. 1982;128(6):2804–6. pmid:6176653
- View Article
- PubMed/NCBI
- Google Scholar

[ref1] 1. Allen TM, Altfeld M, Geer SC, Kalife ET, Moore C, O'Sullivan K M, et al. Selective escape from CD8+ T-cell responses represents a major driving force of human immunodeficiency virus type 1 (HIV-1) sequence diversity and reveals constraints on HIV-1 evolution. J Virol. 2005;79(21):13239–49. pmid:16227247
View Article
PubMed/NCBI
Google Scholar

[2] View Article

[3] PubMed/NCBI

[4] Google Scholar

[ref2] 2. Kawashima Y, Pfafferott K, Frater J, Matthews P, Payne R, Addo M, et al. Adaptation of HIV-1 to human leukocyte antigen class I. Nature. 2009;458(7238):641–5. pmid:19242411
View Article
PubMed/NCBI
Google Scholar

[6] View Article

[7] PubMed/NCBI

[8] Google Scholar

[ref3] 3. Leslie A, Kavanagh D, Honeyborne I, Pfafferott K, Edwards C, Pillay T, et al. Transmission and accumulation of CTL escape variants drive negative associations between HIV polymorphisms and HLA. J Exp Med. 2005;201(6):891–902. pmid:15781581
View Article
PubMed/NCBI
Google Scholar

[10] View Article

[11] PubMed/NCBI

[12] Google Scholar

[ref4] 4. Liu Y, McNevin J, Cao J, Zhao H, Genowati I, Wong K, et al. Selection on the human immunodeficiency virus type 1 proteome following primary infection. J Virol. 2006;80(19):9519–29. pmid:16973556
View Article
PubMed/NCBI
Google Scholar

[14] View Article

[15] PubMed/NCBI

[16] Google Scholar

[ref5] 5. Moore CB, John M, James IR, Christiansen FT, Witt CS, and Mallal SA. Evidence of HIV-1 adaptation to HLA-restricted immune responses at a population level. Science. 2002;296(5572):1439–43. pmid:12029127
View Article
PubMed/NCBI
Google Scholar

[18] View Article

[19] PubMed/NCBI

[20] Google Scholar

[ref6] 6. Crawford H, Lumm W, Leslie A, Schaefer M, Boeras D, Prado JG, et al. Evolution of HLA-B*5703 HIV-1 escape mutations in HLA-B*5703-positive individuals and their transmission recipients. J Exp Med. 2009;206(4):909–21. pmid:19307327
View Article
PubMed/NCBI
Google Scholar

[22] View Article

[23] PubMed/NCBI

[24] Google Scholar

[ref7] 7. Friedrich TC, Dodds EJ, Yant LJ, Vojnov L, Rudersdorf R, Cullen C, et al. Reversion of CTL escape-variant immunodeficiency viruses in vivo. Nat Med. 2004;10(3):275–81. pmid:14966520
View Article
PubMed/NCBI
Google Scholar

[26] View Article

[27] PubMed/NCBI

[28] Google Scholar

[ref8] 8. Leslie AJ, Pfafferott KJ, Chetty P, Draenert R, Addo MM, Feeney M, et al. HIV evolution: CTL escape mutation and reversion after transmission. Nat Med. 2004;10(3):282–9. pmid:14770175
View Article
PubMed/NCBI
Google Scholar

[30] View Article

[31] PubMed/NCBI

[32] Google Scholar

[ref9] 9. Li B, Gladden AD, Altfeld M, Kaldor JM, Cooper DA, Kelleher AD, et al. Rapid reversion of sequence polymorphisms dominates early human immunodeficiency virus type 1 evolution. J Virol. 2007;81(1):193–201. pmid:17065207
View Article
PubMed/NCBI
Google Scholar

[34] View Article

[35] PubMed/NCBI

[36] Google Scholar

[ref10] 10. Goulder PJ, Brander C, Tang Y, Tremblay C, Colbert RA, Addo MM, et al. Evolution and transmission of stable CTL escape mutations in HIV infection. Nature. 2001;412(6844):334–8. pmid:11460164
View Article
PubMed/NCBI
Google Scholar

[38] View Article

[39] PubMed/NCBI

[40] Google Scholar

[ref11] 11. Fryer HR, Frater J, Duda A, Roberts MG, Investigators ST, Phillips RE, et al. Modelling the evolution and spread of HIV immune escape mutants. PLoS Pathog. 2010;6(11):e1001196. pmid:21124991
View Article
PubMed/NCBI
Google Scholar

[42] View Article

[43] PubMed/NCBI

[44] Google Scholar

[ref12] 12. Navis M, Matas DE, Rachinger A, Koning FA, van Swieten P, Kootstra NA, et al. Molecular evolution of human immunodeficiency virus type 1 upon transmission between human leukocyte antigen disparate donor-recipient pairs. PLoS One. 2008;3(6):e2422. pmid:18560583
View Article
PubMed/NCBI
Google Scholar

[46] View Article

[47] PubMed/NCBI

[48] Google Scholar

[ref13] 13. Avila-Rios S, Ormsby CE, Carlson JM, Valenzuela-Ponce H, Blanco-Heredia J, Garrido-Rodriguez D, et al. Unique features of HLA-mediated HIV evolution in a Mexican cohort: a comparative study. Retrovirology. 2009;6:72. pmid:19664284
View Article
PubMed/NCBI
Google Scholar

[50] View Article

[51] PubMed/NCBI

[52] Google Scholar

[ref14] 14. Brumme ZL, Brumme CJ, Heckerman D, Korber BT, Daniels M, Carlson J, et al. Evidence of differential HLA class I-mediated viral evolution in functional and accessory/regulatory genes of HIV-1. PLoS Pathog. 2007;3(7):e94. pmid:17616974
View Article
PubMed/NCBI
Google Scholar

[54] View Article

[55] PubMed/NCBI

[56] Google Scholar

[ref15] 15. Brumme ZL, John M, Carlson JM, Brumme CJ, Chan D, Brockman MA, et al. HLA-associated immune escape pathways in HIV-1 subtype B Gag, Pol and Nef proteins. PLoS One. 2009;4(8):e6687. pmid:19690614
View Article
PubMed/NCBI
Google Scholar

[58] View Article

[59] PubMed/NCBI

[60] Google Scholar

[ref16] 16. Rousseau CM, Daniels MG, Carlson JM, Kadie C, Crawford H, Prendergast A, et al. HLA class I-driven evolution of human immunodeficiency virus type 1 subtype c proteome: immune escape and viral load. J Virol. 2008;82(13):6434–46. pmid:18434400
View Article
PubMed/NCBI
Google Scholar

[62] View Article

[63] PubMed/NCBI

[64] Google Scholar

[ref17] 17. Carlson JM, Le AQ, Shahid A, and Brumme ZL. HIV-1 adaptation to HLA: a window into virus-host immune interactions. Trends Microbiol. 2015;23(4):212–24. pmid:25613992
View Article
PubMed/NCBI
Google Scholar

[66] View Article

[67] PubMed/NCBI

[68] Google Scholar

[ref18] 18. Cotton LA, Kuang XT, Le AQ, Carlson JM, Chan B, Chopera DR, et al. Genotypic and functional impact of HIV-1 adaptation to its host population during the North American epidemic. PLoS Genet. 2014;10(4):e1004295. pmid:24762668
View Article
PubMed/NCBI
Google Scholar

[70] View Article

[71] PubMed/NCBI

[72] Google Scholar

[ref19] 19. Kinloch NN, MacMillan DR, Le AQ, Cotton LA, Bangsberg DR, Buchbinder S, et al. Population-Level Immune-Mediated Adaptation in HIV-1 Polymerase during the North American Epidemic. J Virol. 2016;90(3):1244–58. pmid:26559841
View Article
PubMed/NCBI
Google Scholar

[74] View Article

[75] PubMed/NCBI

[76] Google Scholar

[ref20] 20. Martin E, Carlson JM, Le AQ, Chopera DR, McGovern R, Rahman MA, et al. Early immune adaptation in HIV-1 revealed by population-level approaches. Retrovirology. 2014;11:64. pmid:25212686
View Article
PubMed/NCBI
Google Scholar

[78] View Article

[79] PubMed/NCBI

[80] Google Scholar

[ref21] 21. Murakoshi H, Koyanagi M, Akahoshi T, Chikata T, Kuse N, Gatanaga H, et al. Impact of a single HLA-A*24:02-associated escape mutation on the detrimental effect of HLA-B*35:01 in HIV-1 control. EBioMedicine. 2018;36:103–12. pmid:30249546
View Article
PubMed/NCBI
Google Scholar

[82] View Article

[83] PubMed/NCBI

[84] Google Scholar

[ref22] 22. Monaco DC, Dilernia DA, Fiore-Gartland A, Yu T, Prince JL, Dennis KK, et al. Balance between transmitted HLA preadapted and nonassociated polymorphisms is a major determinant of HIV-1 disease progression. J Exp Med. 2016;213(10):2049–63. pmid:27551154
View Article
PubMed/NCBI
Google Scholar

[86] View Article

[87] PubMed/NCBI

[88] Google Scholar

[ref23] 23. Carlson JM, Du VY, Pfeifer N, Bansal A, Tan VY, Power K, et al. Impact of pre-adapted HIV transmission. Nat Med. 2016;22(6):606–13. pmid:27183217
View Article
PubMed/NCBI
Google Scholar

[90] View Article

[91] PubMed/NCBI

[92] Google Scholar

[ref24] 24. Carlson JM, Brumme CJ, Martin E, Listgarten J, Brockman MA, Le AQ, et al. Correlates of protective cellular immunity revealed by analysis of population-level immune escape pathways in HIV-1. J Virol. 2012;86(24):13202–16. pmid:23055555
View Article
PubMed/NCBI
Google Scholar

[94] View Article

[95] PubMed/NCBI

[96] Google Scholar

[ref25] 25. Murakoshi H, Koyanagi M, Chikata T, Rahman MA, Kuse N, Sakai K, et al. 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. J Virol. 2017;91(4). pmid:27903797
View Article
PubMed/NCBI
Google Scholar

[98] View Article

[99] PubMed/NCBI

[100] Google Scholar

[ref26] 26. Zhang Y, Peng Y, Yan H, Xu K, Saito M, Wu H, et al. Multilayered defense in HLA-B51-associated HIV viral control. J Immunol. 2011;187(2):684–91. pmid:21670313
View Article
PubMed/NCBI
Google Scholar

[102] View Article

[103] PubMed/NCBI

[104] Google Scholar

[ref27] 27. Goulder PJ, Phillips RE, Colbert RA, McAdam S, Ogg G, Nowak MA, et al. Late escape from an immunodominant cytotoxic T-lymphocyte response associated with progression to AIDS. Nat Med. 1997;3(2):212–7. pmid:9018241
View Article
PubMed/NCBI
Google Scholar

[106] View Article

[107] PubMed/NCBI

[108] Google Scholar

[ref28] 28. Kelleher AD, Long C, Holmes EC, Allen RL, Wilson J, Conlon C, et al. Clustered mutations in HIV-1 gag are consistently required for escape from HLA-B27-restricted cytotoxic T lymphocyte responses. J Exp Med. 2001;193(3):375–86. pmid:11157057
View Article
PubMed/NCBI
Google Scholar

[110] View Article

[111] PubMed/NCBI

[112] Google Scholar

[ref29] 29. Kawashima Y, Kuse N, Gatanaga H, Naruto T, Fujiwara M, Dohki S, et al. Long-term control of HIV-1 in hemophiliacs carrying slow-progressing allele HLA-B*5101. J Virol. 2010;84(14):7151–60. pmid:20410273
View Article
PubMed/NCBI
Google Scholar

[114] View Article

[115] PubMed/NCBI

[116] Google Scholar

[ref30] 30. Kuse N, Akahoshi T, Gatanaga H, Ueno T, Oka S, and Takiguchi M. Selection of TI8-8V mutant associated with long-term control of HIV-1 by cross-reactive HLA-B*51:01-restricted cytotoxic T cells. J Immunol. 2014;193(10):4814–22. pmid:25305317
View Article
PubMed/NCBI
Google Scholar

[118] View Article

[119] PubMed/NCBI

[120] Google Scholar

[ref31] 31. Koga M, Kawana-Tachikawa A, Heckerman D, Odawara T, Nakamura H, Koibuchi T, et al. Changes in impact of HLA class I allele expression on HIV-1 plasma virus loads at a population level over time. Microbiol Immunol. 2010;54(4):196–205. pmid:20377748
View Article
PubMed/NCBI
Google Scholar

[122] View Article

[123] PubMed/NCBI

[124] Google Scholar

[ref32] 32. Naruto T, Gatanaga H, Nelson G, Sakai K, Carrington M, Oka S, et al. HLA class I-mediated control of HIV-1 in the Japanese population, in which the protective HLA-B*57 and HLA-B*27 alleles are absent. J Virol. 2012;86(19):10870–2. pmid:22811530
View Article
PubMed/NCBI
Google Scholar

[126] View Article

[127] PubMed/NCBI

[128] Google Scholar

[ref33] 33. Yagita Y, Kuse N, Kuroki K, Gatanaga H, Carlson JM, Chikata T, et al. Distinct HIV-1 escape patterns selected by cytotoxic T cells with identical epitope specificity. J Virol. 2013;87(4):2253–63. pmid:23236061
View Article
PubMed/NCBI
Google Scholar

[130] View Article

[131] PubMed/NCBI

[132] Google Scholar

[ref34] 34. Ikeda N, Kojima H, Nishikawa M, Hayashi K, Futagami T, Tsujino T, et al. Determination of HLA-A, -C, -B, -DRB1 allele and haplotype frequency in Japanese population based on family study. Tissue Antigens. 2015;85(4):252–9. pmid:25789826
View Article
PubMed/NCBI
Google Scholar

[134] View Article

[135] PubMed/NCBI

[136] Google Scholar

[ref35] 35. Murakoshi H, Akahoshi T, Koyanagi M, Chikata T, Naruto T, Maruyama R, et al. Clinical Control of HIV-1 by Cytotoxic T Cells Specific for Multiple Conserved Epitopes. J Virol. 2015;89(10):5330–9. pmid:25741000
View Article
PubMed/NCBI
Google Scholar

[138] View Article

[139] PubMed/NCBI

[140] Google Scholar

[ref36] 36. Saito S, Ota S, Yamada E, Inoko H, and Ota M. Allele frequencies and haplotypic associations defined by allelic DNA typing at HLA class I and class II loci in the Japanese population. Tissue antigens. 2000;56(6):522–9. pmid:11169242
View Article
PubMed/NCBI
Google Scholar

[142] View Article

[143] PubMed/NCBI

[144] Google Scholar

[ref37] 37. Iizuka A, Nagao T, Shohji H, Nitta A, Yoshihara N, and Matsuhashi T. A longitudinal study on seroconversion of HTLV-III/LAV in Japanese haemophiliacs. Br J Haematol. 1987;65(2):249–50. pmid:3644665
View Article
PubMed/NCBI
Google Scholar

[146] View Article

[147] PubMed/NCBI

[148] Google Scholar

[ref38] 38. Kashiwagi S, Hayashi J, Ikematsu H, Noguchi A, Ikeda K, Kishida K, et al. Prevalence of human immunodeficiency virus (HIV) infection among hemophiliacs in Fukuoka, Japan. Microbiol Immunol. 1988;32(8):851–5. pmid:3143898
View Article
PubMed/NCBI
Google Scholar

[150] View Article

[151] PubMed/NCBI

[152] Google Scholar

[ref39] 39. Altman JD, and Feinberg MB. HIV escape: there and back again. Nat Med. 2004;10(3):229–30. pmid:14991039
View Article
PubMed/NCBI
Google Scholar

[154] View Article

[155] PubMed/NCBI

[156] Google Scholar

[ref40] 40. Bailey J, Blankson JN, Wind-Rotolo M, and Siliciano RF. Mechanisms of HIV-1 escape from immune responses and antiretroviral drugs. Curr Opin Immunol. 2004;16(4):470–6. pmid:15245741
View Article
PubMed/NCBI
Google Scholar

[158] View Article

[159] PubMed/NCBI

[160] Google Scholar

[ref41] 41. Allen TM, Yu XG, Kalife ET, Reyor LL, Lichterfeld M, John M, et al. De novo generation of escape variant-specific CD8+ T-cell responses following cytotoxic T-lymphocyte escape in chronic human immunodeficiency virus type 1 infection. J Virol. 2005;79(20):12952–60. pmid:16188997
View Article
PubMed/NCBI
Google Scholar

[162] View Article

[163] PubMed/NCBI

[164] Google Scholar

[ref42] 42. Feeney ME, Tang Y, Pfafferott K, Roosevelt KA, Draenert R, Trocha A, et al. HIV-1 viral escape in infancy followed by emergence of a variant-specific CTL response. J Immunol. 2005;174(12):7524–30. pmid:15944251
View Article
PubMed/NCBI
Google Scholar

[166] View Article

[167] PubMed/NCBI

[168] Google Scholar

[ref43] 43. Melhem NM, Smith KN, Huang XL, Colleton BA, Jiang W, Mailliard RB, et al. The impact of viral evolution and frequency of variant epitopes on primary and memory human immunodeficiency virus type 1-specific CD8(+) T cell responses. Virology. 2014;450–451:34–48. pmid:24503065
View Article
PubMed/NCBI
Google Scholar

[170] View Article

[171] PubMed/NCBI

[172] Google Scholar

[ref44] 44. Roider J, Kalteis AL, Vollbrecht T, Gloning L, Stirner R, Henrich N, et al. Adaptation of CD8 T cell responses to changing HIV-1 sequences in a cohort of HIV-1 infected individuals not selected for a certain HLA allele. PLoS One. 2013;8(12):e80045. pmid:24312453
View Article
PubMed/NCBI
Google Scholar

[174] View Article

[175] PubMed/NCBI

[176] Google Scholar

[ref45] 45. Bailey JR, Williams TM, Siliciano RF, and Blankson JN. Maintenance of viral suppression in HIV-1-infected HLA-B*57+ elite suppressors despite CTL escape mutations. J Exp Med. 2006;203(5):1357–69. pmid:16682496
View Article
PubMed/NCBI
Google Scholar

[178] View Article

[179] PubMed/NCBI

[180] Google Scholar

[ref46] 46. Qin K, Boppana S, Du VY, Carlson JM, Yue L, Dilernia DA, et al. CD8 T cells targeting adapted epitopes in chronic HIV infection promote dendritic cell maturation and CD4 T cell trans-infection. PLoS Pathog. 2019;15(8):e1007970. pmid:31398241
View Article
PubMed/NCBI
Google Scholar

[182] View Article

[183] PubMed/NCBI

[184] Google Scholar

[ref47] 47. Sun X, Fujiwara M, Shi Y, Kuse N, Gatanaga H, Appay V, et al. Superimposed epitopes restricted by the same HLA molecule drive distinct HIV-specific CD8+ T cell repertoires. J Immunol. 2014;193(1):77–84. pmid:24899498
View Article
PubMed/NCBI
Google Scholar

[186] View Article

[187] PubMed/NCBI

[188] Google Scholar

[ref48] 48. Sun X, Shi Y, Akahoshi T, Fujiwara M, Gatanaga H, Schonbach C, et al. Effects of a Single Escape Mutation on T Cell and HIV-1 Co-adaptation. Cell Rep. 2016;15(10):2279–91. pmid:27239036
View Article
PubMed/NCBI
Google Scholar

[190] View Article

[191] PubMed/NCBI

[192] Google Scholar

[ref49] 49. Brumme ZL, Kinloch NN, Sanche S, Wong A, Martin E, Cobarrubias KD, et al. Extensive host immune adaptation in a concentrated North American HIV epidemic. AIDS. 2018;32(14):1927–38. pmid:30048246
View Article
PubMed/NCBI
Google Scholar

[194] View Article

[195] PubMed/NCBI

[196] Google Scholar

[ref50] 50. O'Brien SJ, Gao X, and Carrington M. HLA and AIDS: a cautionary tale. Trends Mol Med. 2001;7(9):379–81. pmid:11530315
View Article
PubMed/NCBI
Google Scholar

[198] View Article

[199] PubMed/NCBI

[200] Google Scholar

[ref51] 51. International HIVCS, Pereyra F, Jia X, McLaren PJ, Telenti A, de Bakker PI, et al. The major genetic determinants of HIV-1 control affect HLA class I peptide presentation. Science. 2010;330(6010):1551–7. pmid:21051598
View Article
PubMed/NCBI
Google Scholar

[202] View Article

[203] PubMed/NCBI

[204] Google Scholar

[ref52] 52. Chikata T, Murakoshi H, Koyanagi M, Honda K, Gatanaga H, Oka S, et al. Control of HIV-1 by an HLA-B*52:01-C*12:02 Protective Haplotype. J Infect Dis. 2017;216(11):1415–24. pmid:28968792
View Article
PubMed/NCBI
Google Scholar

[206] View Article

[207] PubMed/NCBI

[208] Google Scholar

[ref53] 53. Zou C, Murakoshi H, Kuse N, Akahoshi T, Chikata T, Gatanaga H, et al. Effective Suppression of HIV-1 Replication by Cytotoxic T Lymphocytes Specific for Pol Epitopes in Conserved Mosaic Vaccine Immunogens. J Virol. 2019;93(7). pmid:30674626
View Article
PubMed/NCBI
Google Scholar

[210] View Article

[211] PubMed/NCBI

[212] Google Scholar

[ref54] 54. Honda K, Zheng N, Murakoshi H, Hashimoto M, Sakai K, Borghan MA, et al. Selection of escape mutant by HLA-C-restricted HIV-1 Pol-specific cytotoxic T lymphocytes carrying strong ability to suppress HIV-1 replication. Eur J Immunol. 2011;41(1):97–106. pmid:21182081
View Article
PubMed/NCBI
Google Scholar

[214] View Article

[215] PubMed/NCBI

[216] Google Scholar

[ref55] 55. Kuse N, Rahman MA, Murakoshi H, Tran GV, Chikata T, Koyanagi M, et al. Different Effects of Nonnucleoside Reverse Transcriptase Inhibitor Resistance Mutations on Cytotoxic T Lymphocyte Recognition between HIV-1 Subtype B and Subtype A/E Infections. J Virol. 2015;89(14):7363–72. pmid:25972553
View Article
PubMed/NCBI
Google Scholar

[218] View Article

[219] PubMed/NCBI

[220] Google Scholar

[ref56] 56. Gatanaga H, Hachiya A, Kimura S, and Oka S. Mutations other than 103N in human immunodeficiency virus type 1 reverse transcriptase (RT) emerge from K103R polymorphism under non-nucleoside RT inhibitor pressure. Virology. 2006;344(2):354–62. pmid:16219331
View Article
PubMed/NCBI
Google Scholar

[222] View Article

[223] PubMed/NCBI

[224] Google Scholar

[ref57] 57. Altman JD, Moss PA, Goulder PJ, Barouch DH, McHeyzer-Williams MG, Bell JI, et al. Phenotypic analysis of antigen-specific T lymphocytes. Science. 1996;274(5284):94–6. pmid:8810254
View Article
PubMed/NCBI
Google Scholar

[226] View Article

[227] PubMed/NCBI

[228] Google Scholar

[ref58] 58. Tomiyama H, Akari H, Adachi A, and Takiguchi M. Different effects of Nef-mediated HLA class I down-regulation on human immunodeficiency virus type 1-specific CD8(+) T-cell cytolytic activity and cytokine production. J Virol. 2002;76(15):7535–43. pmid:12097566
View Article
PubMed/NCBI
Google Scholar

[230] View Article

[231] PubMed/NCBI

[232] Google Scholar

[ref59] 59. Lin Z, Kuroki K, Kuse N, Sun X, Akahoshi T, Qi Y, et al. HIV-1 Control by NK Cells via Reduced Interaction between KIR2DL2 and HLA-C (*)12:02/C (*)14:03. Cell Rep. 2016;17(9):2210–20. pmid:27880898
View Article
PubMed/NCBI
Google Scholar

[234] View Article

[235] PubMed/NCBI

[236] Google Scholar

[ref60] 60. Radka SF, Kostyu DD, and Amos DB. A monoclonal antibody directed against the HLA-Bw6 epitope. J Immunol. 1982;128(6):2804–6. pmid:6176653
View Article
PubMed/NCBI
Google Scholar

[238] View Article

[239] PubMed/NCBI

[240] Google Scholar

Figures

Abstract

Author summary

Introduction

Results

Accumulation of RT135V mutation in HLA-B*52:01+HLA-C*12:02+ hemophiliacs with HIV

Generation of novel HLA-C*12:02-restricted epitope after selection of RT135V mutation by HLA-B*52:01-restricted TI8-specific CTLs

Recognition of target cells infected with RT135 mutant viruses by TN9-8V-specific T cells

Population-level shifts in HIV-1 RT135 mutation frequencies Japan over two decades

Selection of RT135L by TN9-8V-specific HLA-C*12:02-restricted T cells

Suppression of RT135T mutant virus by TN9-8T-specific HLA-C*12:02-restricted T cells

Discussion

Methods

Ethics statement

Participants

Cells

Induction of peptide-specific bulk T cells

ICS assay

Generation of TN9-8V- and TN9-8T-specific CTL lines

HIV-1 clones

Generation of HLA-peptide tetrameric complexes

Tetramer binding assay

Viral replication suppression assay

Sequence of autologous virus

HLA class I stabilization assay

Statistical analysis

Supporting information

S1 Fig. Nucleotide substitution at RT135 in Japanese individuals with chronic HIV-1, Related to Fig 1.

S2 Fig. Identification of a novel HLA-C*12:02-restricted CTL epitope including RT135, Related to Fig 2A.

S3 Fig. Negative controls of ex vivo tetramer staining assay, Related to Figs 2B, 5B and 6B.

S4 Fig. Recognition of virus-infected cells by TN9-8V-specific T cell lines, Related to Fig 3C.

S5 Fig. Ex vivo detection of TN9-8T-specific CD8+ T cells, Related to Fig 6G.

S1 Table. Codon usages of amino acids at RT135 observed in Japanese hemophiliacs and non-hemophiliacs with chronic HIV-1

S2 Table. Statistical analysis using Cochran-Mantel-Haenszel test on the association between the frequency of 7 amino acids and the four periods.

Acknowledgments

References

Accumulation of RT135V mutation in HLA-B52:01⁺HLA-C12:02⁺ hemophiliacs with HIV

Generation of novel HLA-C12:02-restricted epitope after selection of RT135V mutation by HLA-B52:01-restricted TI8-specific CTLs

S5 Fig. Ex vivo detection of TN9-8T-specific CD8⁺ T cells, Related to Fig 6G.