- Open Access
Lack of a significant impact of Gag-Protease-mediated HIV-1 replication capacity on clinical parameters in treatment-naive Japanese individuals
Retrovirologyvolume 12, Article number: 98 (2015)
The Erratum to this article has been published in Retrovirology 2015 12:100
HLA class I-associated escape mutations in HIV-1 Gag can reduce viral replication, suggesting that associated fitness costs could impact HIV-1 disease progression. Previous studies in North American and African cohorts have reported reduced Gag-Protease mediated viral replication capacity (Gag-Pro RC) in individuals expressing protective HLA class I alleles including HLA-B*57:01, B*27:05, and B*81:01. These studies also reported significant positive associations between Gag-Pro RCs and plasma viral load (pVL). However, these HLA alleles are virtually absent in Japan, and the importance of Gag as an immune target is not clearly defined in this population.
We generated chimeric NL4-3 viruses carrying patient-derived Gag-Protease from 306 treatment-naive Japanese individuals chronically infected with HIV-1 subtype B. We analyzed associations between Gag-Pro RC and clinical markers of HIV-1 infection and host HLA expression. We observed no significant correlation between Gag-Pro RC and pVL in Japan in the overall cohort. However, upon exclusion of individuals expressing Japanese protective alleles HLA-B*52:01 and B*67:01, Gag-Pro RC correlated positively with pVL and negatively with CD4 T-cell count. Our results thus contrast with studies from other global cohorts reporting significantly lower Gag-Pro RC among persons expressing protective HLA alleles, and positive relationships between Gag-Pro RC and pVL in the overall study populations. We also identified five amino acids in Gag-Protease significantly associated with Gag-Pro RC, whose effects on RC were confirmed by site-directed mutagenesis. However, of the four mutations that decreased Gag-Pro RC, none were associated with reductions in pVL in Japan though two were associated with lower pVL in North America.
These data indicate that Gag fitness does not affect clinical outcomes in subjects with protective HLA class I alleles as well as the whole Japanese population. Moreover, the impact of Gag fitness costs on HIV-1 clinical parameters in chronic infection is likely low in Japan compared to other global populations.
Within an infected individual, immune pressure by cytotoxic T lymphocytes (CTLs) induces the selection of HLA-associated escape variants in the HIV-1 genome. The persistence of HLA-associated escape variants upon transmission in turn shapes viral diversity at the population level . CTL escape occurs throughout infection. Moreover, certain escape mutations compromise viral fitness. In some of these cases, their effects on immune evasion can be offset in part by their fitness costs though the effects of such mutations on plasma viral load (pVL) are complex and remain incompletely understood [2–9]. HIV-1 Gag is rich in immunodominant CTL epitopes, and a vital role of anti-Gag CTL responses on HIV-1 immune control has been demonstrated [10, 11; also reviewed in 12, 13]. In particular, associations between fitness-costly Gag escape mutations and decreases in pVL have been described [2, 5–8], notably for the immunodominant HLA-B*57:01-restricted TW10 and B*27:05-restricted KK10 epitopes, both of which are well studied because of their protective roles in HIV-1 disease outcome. Specifically, CTL-driven escape mutations in TW10 impose a mild fitness cost that, along with robust HLA-B*57:01-driven CTL responses, likely contribute to the significantly lower pVL observed in these subjects [5, 7, 14–17]. Notably however, fitness costs of escape in HLA-B*57:01 epitopes are generally rescued by the subsequent appearance of compensatory mutations [18, 19]. In contrast, a single amino acid change at Gag 264 (R264K) in the B*27:05-restricted KK10 epitope substantially compromises in vitro viral replication [20–22], but this mutation generally occurs late in infection due to the requirement for the upstream compensatory mutation S173A [7, 21]. As such, the appearance of R264K in vivo generally heralds HIV-1 progression [23, 24]. These studies support a relationship between Gag escape mutations and pVL and underscore a critical role of Gag in HIV-1 disease outcome.
In vitro studies also highlight the importance of Gag as an immune target. Certain HLA-associated mutations in Gag have been shown to impair Gag-Protease-mediated viral replication capacity in vitro (Gag-Pro RC). Gag-Pro RC is in turn positively associated with pVL in numerous global cohorts [18, 25–30]. In particular, reduced Gag-Pro RC has been observed in predominantly North American cohorts of elite controllers as well as individuals in acute/early infection who express protective alleles B*13, B*27, B*57, and/or B*58:01 [18, 26], though the latter associations are less apparent during the chronic phase . In addition, the impact of HLA-associated changes in Gag-Pro RC on pVL have also been reported in African cohorts chronically infected with HIV-1 subtype C [27–29]. A more recent study compared Gag-Pro RCs in treatment-naive Mexican and Barbadian HIV-1 subtype B cohorts which differ in their frequency of protective HLA alleles (8 % versus 34 % of B*27/57/58:01/81:01 in Mexico versus Barbados). The study revealed a higher burden of HLA-mediated Gag escape mutations and significantly lower Gag-Pro RC in chronically infected Barbadian compared to Mexican individuals . Moreover, Gag-Pro RC correlated positively with pVL, supporting the importance of Gag in HIV-1 disease outcome . The significantly lower Gag-Pro RC and pVL in Barbados versus Mexico was attributed to the higher prevalence of protective HLA alleles in the former cohort, implying that protective alleles lower Gag-Pro RC, viral set points, and therefore, disease progression rates at the population level. These results are consistent with a previous study of nine distinct cohorts worldwide that demonstrates HIV-1 adaptation to HLA class I alleles at a population level, presumably as a result of transmission and persistence of certain escape mutations in the population . Thus, HLA allele frequency shapes viral diversity, which could also affect CTL responses and HIV-1 disease outcome in the population.
HLA class I distributions differ markedly in Asian compared with North American and African cohorts . The “classical” protective alleles B*27:05 and B*57:01 are essentially absent in Japan: a published cohort of 504 treatment-naive HIV+ Japanese featured only one B*27:05+ and no B*57:01 patients [31–37]. Instead, B*52:01 and B*67:01 are associated with HIV-1 control in Japan . Therefore, the role of Gag fitness in HIV-1 pathogenesis, or more broadly the mechanism of protection, is likely to be different in Asian populations. As such, the impact of HLA-mediated immune pressure on viral fitness and its influence on HIV-1 pathogenesis need to be specifically defined for Asian populations. To address these issues, we investigated the effect of HLA-associated changes in Gag-Pro RC in 306 treatment-naive HIV+ Japanese patients chronically infected with HIV-1 subtype B in the presence or absence of protective HLA alleles. We also identified amino acid changes in Gag-Protease that altered Gag-Pro RCs and examined their associations with pVL. In contrast to studies of non-Asian cohorts that supported a role of Gag fitness in HIV-1 disease outcome, our results suggested that Japanese protective alleles did not lower Gag-Pro RC and that Gag fitness did not correlate with clinical markers in the whole Japanese population.
Relationship between Gag-Protease-dependent replication capacity and clinical markers of disease progression
We generated recombinant HIV-1NL4-3 containing Gag-Protease from 306 treatment-naïve Japanese subjects chronically infected with HIV-1 subtype B. Plasma viral sequences were determined for all patients in a previous study , and the patient origin of all recombinant viruses was confirmed by maximum likelihood phylogenetic analysis (Additional file 1). Recombinant viruses were used to infect a reporter CD4+ T-cell line, CEM-GXR, which naturally expresses CXCR4 and has been further engineered to express CCR5 and Tat-inducible GFP . Gag-Protease-dependent replication capacities (Gag-Pro RCs) were determined based on the exponential increase in the % GFP+ cells during days 2–8 following infection and normalized against that of the parental HIV-1NL4-3 as described for North American and African studies [18, 25–30]. Similar to previous reports, Gag-Pro RCs were normally distributed (Fig. 1a), with a mean of 1.12, a standard deviation (SD) of 0.13, and a range of 0.67–1.56. The mean Gag-Pro RC in the Japanese cohort was slightly higher than that observed in the North American cohort  but comparable to those reported more recently for an ART-naïve Mexican cohort .
To examine the impact of Gag-Pro RC on HIV-1 pathogenesis, we investigated the relationship between Gag-Protease RC and clinical markers of HIV-1 infection (pVL and CD4 T-cell count). Overall, Japanese Gag-Pro RCs did not correlate with pVL (R = 0.092, P = 0.11; Fig. 1b) nor CD4 T-cell count (R = −0.078, P = 0.17; Fig. 1c). Not only did these results not achieve statistical significance, but the correlation coefficients for pVL and Gag-Pro RC in Japan were weaker than those reported in African (R = 0.24) [28, 29] and Mexican/Barbadian (R = 0.25)  chronic cohorts. The Japanese correlation coefficient R = 0.092 was also weaker to that reported in a North American chronic cohort (R = 0.27 in the presence of protective HLA alleles; R = 0.12 for the total population) . Together, this result suggested a minimal role for Gag fitness as a major determinant of virologic control in Japan compared to other global cohorts.
Impact of HLA class I alleles on Gag-Pro RCs and relationship with clinical parameters
To investigate the role of Gag as an HLA-mediated immune target, we stratified HLA class I alleles based on their median RCs (Additional file 2A-C). For all alleles with frequencies >3 % (n > 8) in the cohort, we did not find any significant associations between HLA alleles and Gag-Pro RC. In addition, for each HLA allele, we compared the relative rankings of median RCs amongst the HLA alleles with those previously determined for pVL . If Gag-Pro RC affects pVL, HLA alleles with lower Gag-Pro RCs would show a tendency toward lower pVL and vice versa. However, the relative RC ranking (in Additional file 2A-C) and pVL for each HLA allele did not match in most cases. We previously reported that A*33:03, B*44:03, B*40:06, C*14:03, and C*12:02 were associated with lower pVL in Japan , but none of these alleles were associated with lower Gag-Pro RC (Additional file 2A-C). Notably, the median Gag-Pro RCs in individuals expressing Japanese protective alleles HLA-B*52:01 and B*67:01  were higher than the population median (Additional file 2B), indicating that protective-allele expression is not associated with lower Gag fitness in Japan. Indeed, no significant correlation was observed between median Gag-Pro RCs and median pVL for each HLA allele (Additional file 2D).
To further explore the relationship between HLA, Gag-Pro RC, and HIV-1 clinical parameters, we grouped our cohort by HLA-allele expression and analyzed associations between Gag-Pro RCs and HIV-1 clinical parameters in each HLA group (Table 1). In two HLA-A alleles (A*02:01 and A*31:01), one B allele (B*07:02), and two C alleles (C*08:01 and C*08:03), Gag-Pro RC exhibited a positive correlation with pVL (p < 0.05), whereas Gag-Pro RC exhibited negative correlations with CD4 count at this threshold in seven alleles (HLA-A*02:06, HLA-B*46:01, B*48:01, B*59:01, C*01:02, C*08:01, and C*08:03; Table 1). However, only in HLA-C*08:01+ (n = 34) and C*08:03+ (n = 10) subjects, Gag-Pro RC correlated with both clinical parameters: a positive correlation with pVL (R = 0.36 and R = 0.70, respectively) and an inverse correlation with CD4 count (R = −0.46 and R = −0.69, respectively; Fig. 2).
We also examined groups of patients based on non-expression of particular HLA alleles. We observed significant correlations between clinical parameters and Gag-Pro RCs in the absence of eight HLA alleles (Table 1), including HLA-B*52:01. In individuals lacking the Japanese protective B*52:01 allele, Gag-Pro RCs correlated positively with pVL (R = 0.15, P = 0.021, q = 0.24, n = 234) and inversely with CD4 T-cell count (R = −0.19, P = 0.005, q = 0.16; Fig. 2). Similarly, Gag-Pro RC also showed significant associations with both pVL and CD4 count in individuals lacking HLA-C*12:02 (Fig. 2). HLA-B*52:01 and C*12:02 form a haplotype in Japan [31–36]. Indeed, all B*52:01+ subjects had C*12:02, and of 85 C*12:02+ subjects in our cohort, only three did not carry B*52:01. Positive correlations between Gag-Pro RC and pVL were additionally observed in the absence of A*11:01, A*33:03, B*39:01, B*44:03, C*07:02, and C*14:03 (Table 1), but none of these alleles associated with CD4 count. Taken together, our results stand in contrast to observations of lower Gag-Pro RC among individuals expressing classical protective HLA alleles in other global cohorts [18, 25–30]. HIV-1 control by protective HLA alleles is not associated with low Gag fitness in Japan.
Impact of Gag on clinical parameters in the absence of protective alleles
In our cohort, Gag-Pro RCs did not correlate with pVL in subjects carrying either of the Japanese protective alleles, B*52:01 or B*67:01. However, exclusion of B*52:01 (n = 234) resulted in associations of Gag-Pro RCs with clinical markers of infection (R = 0.15, P = 0.021, q = 0.24 for pVL; Table 1). Therefore, we next examined correlations between Gag-Pro RCs and clinical parameters in the absence of these protective alleles (i.e. among B*52:01 and B*67:01 double-negative subjects; n = 224). In persons lacking both of these protective alleles, Gag-Pro RCs exhibited a positive association with pVL (R = 0.16, P = 0.017) and a negative association with CD4 count (R = −0.21, P = 0.0019) (Fig. 3a). These associations are weaker than those observed in Mexican/Barbadian, African, and North American individuals expressing protective HLA alleles (R > 0.2 for these cohorts) but comparable in magnitude to the relationship between Gag-Pro RC and pVL in the North American chronic cohort (R = 0.12, P = 0.0007) [18, 29, 30]. Indeed, after excluding B*52/B*67/C*12 alleles, the associations of Gag-Pro RC and clinical parameters became slightly stronger in each HLA allele- group analyzed in Table 1 (Additional file 3). These observations further support our result that Gag fitness does not associate with clinical parameters in subjects with protective HLA alleles in Japan.
When we further explored the relationship between Gag-Pro RC and HIV clinical parameters in the B*52:01−B*67:01− subpopulation by multivariable analyses, a linear regression model adjusting for pVL and CD4 count indicated that Gag-Pro RC in this subset was driven largely by CD4 count (P = 0.043; β Estimate = −0.012 per 100 cell/mm3 increment) and not pVL (P = 0.47; β Estimate = 0.0091 per log10 increment). In other words, the model indicated that, after controlling for pVL, Gag-Pro RC decreased by 0.012 units for each increment of 100 CD4 cells/mm3. Although the model identified CD4 count as an independent predictor of Gag-Pro RC, CD4 count explains only 3.7 % of population variation in Gag-Pro RC (Table 2). In another linear regression model additionally incorporating three HLA alleles—C*08:01, C*08:03, and C*12:02—that exhibited significant associations of Gag-Pro RC with both pVL and CD4 count in the B*52:01−B*67:01− subpopulation, the addition of HLA information to the model did not improve the model fit (P = 0.12). Taken together, results of the multivariable analyses suggest that apparent HLA associations with Gag-Pro RC, HIV clinical parameters, and HLA alleles are likely to be driven primarily by CD4 count rather than pVL in subjects lacking the protective B*52 and B*67 alleles. Moreover, these observations raise the possibility that Gag-Pro RC may not have a major impact on pVL even in the B*52−B*67− subpopulation.
We next analyzed the impact of HLA associated polymorphisms on pVL in 218 B*52:01− B*67:01− subjects based on the HLA associated polymorphisms determined for the Japanese population . Here, HLA-associated polymorphisms were defined based on a published list derived from plasma virus sequences from the present cohort (see methods and ; these included 94 HLA-associated polymorphisms in Gag and 16 in Protease). For each individual in this analysis, we counted the number of HLA-associated polymorphisms present in their viral sequence (regardless of the host HLAs expressed). In doing so we observed no correlation between the number of polymorphisms in Gag, Protease, or Gag-Protease and Gag-Pro RCs (Fig. 3b–d). On the other hand, statistically significant negative correlations were observed between the number of HLA-associated polymorphisms in Gag (R = −0.15, P = 0.024) as well as Gag-Protease (R = −0.16, P = 0.018) and pVL (Fig. 3e–g). This suggests that HLA-associated polymorphism burden in Gag impacts pVL without significantly affecting Gag-Pro RCs in the B*52−B*67− subpopulation, indicating the possible absence of correlation between Gag-Pro RC and pVL in this subpopulation. A simple explanation for the discrepancy would be that a combinatorial effect by mutations in other parts of the virus may have resulted in the significant change in pVL even though the number of polymorphisms in Gag only was not enough to affect Gag replication per se.
Impact of amino acid changes on Gag RCs in the absence of the protective alleles
We expanded our search for mutations in Gag that would affect Gag-Pro RCs in the B*52:01− B*67:01− subjects. In an exploratory analysis, we statistically identified amino acids in Gag that were associated with Gag-Pro RCs in persons lacking the protective B*52:01 and B*67:01 alleles. We identified amino acid variants at four positions, Gag 79, 228, 286, and 357 that were associated with lower Gag-Pro RC and at one position, Gag 218, which was associated with a higher Gag-Pro RC (Table 3). To evaluate the potential combinatorial effects of these mutations on Gag-Pro RC, we assigned a score to each recombinant Gag-Protease viral sequence by assigning +1 for the presence of the RC-increasing residue and −1 for each of the four RC-decreasing residues and computing this sum. Overall, these scores exhibited a significant positive correlation with Gag-Pro RCs (R = 0.35, P = 0.03 × 10−6; Fig. 4a). This suggests a dose-dependent effect of these substitutions on replication capacity, which we sought to confirm in vitro. To do this, we focused on the amino acid changes that reduced Gag-Pro RC and introduced Gag Y79F and R286K substitutions (the most frequent variations at these positions; Additional file 4), along with Gag M228L into pNL4-3 by site-directed mutagenesis. NL4-3 naturally carries glycine at Gag357 and therefore naturally contains a residue other than S at this site (357G). Note that the combination of these four substitutions 79F+228L+286K+357G was observed in naturally occurring sequences from our cohort. The mutant NL4-3 carrying these four Gag-Pro RC-reducing changes exhibited 20 % decreased replication capacity compared to wild-type (P < 0.0001; Fig. 4b). Thus, the in vitro result confirmed the effect of the amino acid changes on the viral replication capacity.
Even though the presence of these amino acid changes in the HLA-B*52:01−B*67:01− subpopulation reduced viral replication capacity in vitro, they did not significantly correlate with pVL in this subpopulation (Fig. 5a). To investigate the potential in vivo importance of these amino acid changes in a non-Japanese cohort, we analyzed associations between these variants and pVL in a North American cohort that previously reported stronger positive correlations between Gag-Pro RC and pVL . In the North American cohort, B*52:01/C*12:02 prevalence is <2 % and B*67:01 prevalence is 0 %, and thus these alleles are unlikely to influence the present analysis . In contrast to the Japanese cohort, the North American cohort showed a weak association between pVL and mutation scores (R = 0.078, P = 0.049; Fig. 5b), suggesting a minor but significant impact of these polymorphisms in vivo. Individual analysis of mutations at each position also revealed significantly lower pVL among individuals harboring amino acid variants at Gag 79 (Y79 versus non-Y) and 228 (M228L) in the North American cohort (Fig. 5d). In contrast, no such relationships were observed in the HLA-B*52:01−B*67:01− Japanese subpopulation (Fig. 5c). The observation that Gag amino acid variants, at least at position 79 and 228, are associated with pVL in North America but not Japan further supports a lesser role of Gag fitness in modulating clinical outcomes in Japan, even in the HLA-B*52:01−B*67:01− subpopulation.
Mounting evidence suggests the importance of immune responses against Gag in the control of HIV-1 infection [reviewed in 12, 13]. In persons expressing protective HLA-B*57:01 and B*27:05 alleles, targeting critical immunodominant epitopes in Gag slows disease progression. In vitro viral replication studies, notably those using recombinant viruses carrying patient-derived Gag-Protease sequences, have contributed to this understanding by demonstrating significantly lower Gag-dependent viral replication in HIV-1 controllers from North American cohorts [25, 26]. Positive correlations between Gag-mediated replication capacity and pVL have also been reported in HIV-1 subtype C in an African cohort [28, 29] and subtype B in Mexican and Barbadian cohorts  which indicated an association between Gag-Pro RC and the frequency of protective HLA alleles in the population.
In contrast, our results indicate no significant association of Gag-Pro RC with pVL and CD4 in the overall Japanese cohort even though identical experimental methods were employed. This lack of significant correlation between Gag-Pro RC and pVL, possibly even in the B*52-B*67- subpopulation according to the multivariate analysis, suggests that Gag fitness does not substantially affect pVL in Japan and raises the hypothesis that Gag-mediated HIV-1 control is less critical in Japan than other populations. The difference could be due to the unique HLA distribution in Japan , which shapes viral diversity at a population level . Particularly, the protective HLA-B*57:01 and B*27:05 alleles are virtually absent in Japan; instead, HLA-B*52:01 (allele frequency 10.4–11.1 %) and Asian-specific B*67:01 (allele frequency 1.1–1.7 %) are associated with HIV-1 control [Allele Frequencies database (http://www.allelefrequencies.net/default.asp); 32]. Intriguingly, neither B*52:01 nor B*67:01 was associated with low Gag-Pro RCs, and Gag-Pro RC did not correlate with clinical parameters in subjects with these protective alleles. As such, we hypothesize that the replicative consequences of B*52/B*67-driven immune escape mutations in Gag are minimal compared to the replicative effects of escape from protective alleles such as B*27 and B*57 in other populations. In support of this hypothesis, our previous results indicate that the number of HLA-associated amino acid changes in Pol, but not in Gag, negatively correlate with pVL . The correlation was observed in the presence of HLA-B*52:01, suggesting that escape from HLA-B*52:01-restricted responses in Pol may confer fitness costs.
The absence of HLA-C*12:02, an allele in tight linkage disequilibrium with HLA-B*52:01, also showed significant associations between RC and both clinical parameters. Our previous study demonstrated a negative correlation between pVL and B*52:01-associated amino acid changes in Pol whereas C*12:02-associated changes did not show a significant association . However, the effect of B*52:01 cannot be separated in the present study due to the extremely low number of B*52:01 or C*12:02 single positive subjects. The presence of C*08:01 (n = 34) and C*08:03 (n = 10) also had strong associations with both clinical parameters, but these two alleles are not significantly associated with lower or higher pVL in chronic phase in our previous study . Therefore, the clinical importance and a putative role of these associations require further investigation.
We identified amino acid polymorphisms in Gag that were associated with alterations of Gag-Pro RC in the absence of protective HLA-B*52:01 and B*67:01 alleles. Of the four polymorphisms associated with lower RC, one was located in p17 (Y79) and the other three in p24 (M228, R286, G357). Our previous analysis of HLA-associated viral polymorphisms in this cohort indicated that changes at Y79 were associated with HLA-B*51:01, C*14:02, and C*14:03, changes at M228 were associated with C*0304, and changes at G357 were associated with B*07:02 . Of note, R286 had been previously identified as B*52:01-associated even though B*52-expressing subjects were excluded from the present analysis of RC-associated substitutions . It is therefore possible that R286K imposes a fitness cost in B*52:01+ subjects, but its effect on Gag-Pro RC may have been indiscernible due to the presence of other RC-modulating mutations in these individuals. Moreover, the impact of R286K on pVL among B*52:01 individuals may be difficult to detect due to two conflicting forces at play: on one hand, R286K may slightly lower replication capacity, but on the other hand, it confers the ability to escape from CTL responses and thus persist in vivo. For B*52:01− subjects, the B*52-driven R286K immune-escape mutation may have been transmitted to B*52:01− recipients, which may have resulted in the observed lower Gag-Pro RC (similar to the reduced RCs observed when the B*57:01-restricted T242N is transmitted to B*57:01− recipients [8, 9] ). Site-directed mutagenesis confirmed that the four amino acid changes significantly reduced Gag-Pro RC. However, the reduction in Gag-Pro RCs by these amino acid changes did not lead to lower pVL even in the B*52−B*67− subpopulation, suggesting a possible lack of a significant impact of Gag fitness on pVL in Japan.
Even though our data argue against a significant role of Gag fitness as a major determinant of HIV-1 disease outcome in Japan, they do not completely exclude the possibility that Gag fitness subtly influences pVL. The amino acid changes may have significantly decreased pVL, but other mutations in viral genome could have restored viral replication capacity in vivo. These scenarios still suggest that pVL does not reflect change in Gag fitness, and that Gag replication has only a minimal impact on HIV-1 pathogenesis in Japan. Further studies, including longitudinal analyses, will be necessary to identify such compensatory mutations. Alternatively, HLA-associated polymorphisms in Gag could actually be the result of immune pressure. In this case, a high immune pressure in Gag would generate immune-escape mutations, but the reduction of pVL may have occurred through cross-reactive CTL function rather than mutation-induced Gag fitness cost. A recent functional study on CTL responses reported 13 epitopes (eight in Gag and five in Pol) involved in HIV-1 control in Japan . B52* and B*67 targeted eight out of the 13 epitopes, including seven of the eight epitopes located in Gag, suggesting the importance of the Gag as an immune target by the Japanese protective alleles . Surprisingly, however, only one of the eight disease-controlling CTL epitopes targeted by B*52 induced escape mutations . These results suggest that CTL responses directly control pVL by killing infected cells via HLA-mediated recognition of Gag epitopes, rather than via escape-mediated fitness costs. In addition to sustained responses to Gag epitopes that are difficult for the virus to escape , HLA-mediated protective effects in Japan may also be mediated by responses to Pol epitopes that eventually induce fitness-costly escape mutations within that viral protein, that partially offset the advantages gained by escape from CTL responses .
Gag fitness has an effect on CD4 count in the Japanese population who do not express protective alleles, in contrast to global cohorts. At least in this subpopulation, escape mutations in Gag may impose fitness costs that are clinically advantageous. However, viral polymorphisms that weakly affected pVL did not have an impact on Gag fitness even in this B*52−B*67− subpopulation, implying a minor contribution of Gag fitness to HIV-1 clinical outcome. Fitness costs of immune-driven mutations in other viral proteins may have more significant impact on HIV-1 disease progression in Japan . These observations are in contrast to the previous reports which indicated significant correlations between Gag-Pro RCs and clinical parameters in Mexican/Barbadian and African chronic cohorts as well as elite controllers and subjects carrying protective HLA alleles in predominantly Caucasian chronic cohorts [18, 25, 27–30].
Previously, one study analyzed Gag-Pro RC in 156 treatment-naive Japanese individuals sampled from 1994 to 2009 and reported reduced Gag-Pro RCs among more recent isolates . In contrast to our results, they observed a significant association between Gag-Pro RC and pVL. The discrepancy between these two studies may be due to the characteristics of the cohorts and/or study time course. The present study analyzed samples collected from 2008 to 2011 whereas the previous one focused on older samples. Newer viruses may have more chance to accumulate mutations that could affect viral replication dependently or independently of those occurring in Gag-Protease. Mutations located outside of Gag-Protease could affect pVL but not Gag-Pro RCs measured in vitro, which would reduce the correlation between Gag-Pro RC and pVL. Sequence comparisons identified significantly different amino acid frequencies at 35 out of 500 amino acid positions in Gag and eight out of 100 positions in Pol between the two cohorts, which may have affected different outcomes in the two cohorts (data not shown). In addition, cohort-level differences in HLA prevalence (e.g. in the other study, the frequencies of A*24+ and C*03+ individuals were 52.6 and 34.6 %, respectively whereas the frequencies of these alleles in the present study were 62.4 and 42.2 %, respectively) may also have affected the results.
Notably, an association between Pol replication capacity and HIV-1 disease progression has also been reported in North American cohorts [41–46]. One of these studies suggested decreased reverse transcriptase-integrase-dependent RC in HIV-1 elite controllers . In another study B*52:01-associated change at amino acid I220 in integrase was shown to decrease viral replication capacity . Considering that Pol enzymatic functions directly affect viral replication, mutations in Pol may meaningfully impact viral fitness and pVL. Possibly, HLA-associated HIV-1 polymorphisms in Gag-Protease may not affect overall HIV-1 disease outcome in Japan as only a few alleles, including HLA-C*08:01 and C*08:03, had relatively strong correlations between Gag-Pro RC and clinical markers of infection. Additionally, the present study comprised a cross-sectional analysis at the chronic phase of HIV infection as previously done in the other global cohort studies except for the North American cohort. Thus, we cannot exclude a possibility of a significant impact of Gag fitness on clinical outcome at different time points: e.g. at the acute/early phase. Further study, including longitudinal analyses, will be necessary to determine a possible role of Gag fitness at the earlier phase as well as Pol fitness in HIV-1 disease progression in Japan.
Our data support a lack of a significant impact of Gag fitness on pVL in treatment-naive Japanese individuals chronically infected with HIV-1. Moreover, in contrast to non-Asian chronic cohorts where Gag-Pro RC fitness may modulate HIV clinical parameters especially in subjects with protective HLA alleles, the protective effects of Japanese B*52 and B*67 alleles do not seem to be associated with Gag fitness. Thus, the impact of Gag fitness on HIV clinical parameters is population dependent.
A total of 430 treatment-naive Japanese subjects chronically infected with HIV-1 subtype B were recruited for a previous study  and enrolled at the National Center for Global Health and Medicine (NCGHM) from 2008 to 2011. HLA typing and the preparation of cDNA samples from patient plasma were performed as described in the previous study. For the present study, specimens from 306 of the original 430 participants were randomly selected for analysis. The median pVL and CD4 count for this N = 306 cohort were 4.5 log10 RNA copies/ml (IQR: 4.1–5.1) and 299 cells/mm3 (IQR: 151–407), respectively. The study was approved by the Ethical Committees in the Faculty of Life Science, Kumamoto University and the NCGHM and conducted in accordance with the Declaration of Helsinki.
For comparison with the Japanese cohort, a data set derived from a published North American (Canada and USA) cohort comprising 803 untreated chronic subjects infected with HIV-1 subtype B was used for the analysis . Median pVL and CD4 count for the cohort were 5.1 log10 RNA copies/ml (IQR: 4.7–5.5) and 273 cells/mm3 (IQR: 130–420), respectively .
Generation of chimeric viruses and in vitro replication assays
The production of chimeric viruses and Gag-Pro RC determination was performed as described previously [18, 25–30]. Briefly, cDNAs were generated from plasma viral RNA for a separate study . Plasma-derived gag-protease DNA fragments were generated from the cDNAs by nested PCR, using TaKaRa ExTaq Hot Start Version (Takara Bio Inc., Shiga, Japan). Chimeric HIV-1NL4-3 encoding plasma-derived gag-protease was generated by electroporating gag-protease PCR fragments with linearized gag-protease-deficient pNL4-3 (delta gag-pro NL4-3) into a CEM-derived GFP-reporter T-cell line engineered to carry LTR-GFP, CXCR4, and CCR5 (CEM-GXR, a generous gift from T. Miura and M.A. Brockman)  . Stock titers were determined by infecting CEM-GXR and measuring GFP+ cells, using flow cytometry. To determine Gag-Pro RCs, CEM-GXR cells were infected with chimeric viruses at a MOI of 0.03 at day 2, and the percentage of HIV-infected (GFP-expressing) cells was monitored daily for 2–8 days post-infection. Gag-Pro RC was calculated as the natural-log slope of the percentage of GFP-expressing cells during the exponential phase of viral spread and normalized against the mean value obtained with wild-type NL4-3 (i.e. NL4-3 RC = 1.00). Virus titration and RC experiments were done in triplicate.
For the validation of virus stocks, we isolated RNA from all virus stocks, using EZ1 Virus Mini Kit v2.0 (Qiagen, Valencia, CA). Gag-protease DNA fragments were generated by RT-PCR, using SuperScript III First-Strand Synthesis System for RT-PCR (Life Technologies, Carlsbad, CA) for cDNA synthesis and TaKaRa ExTaq Hot Start Version for subsequent nested PCR. Viral sequences were analyzed on the ABI3500 genetic analyzer (Applied Biosystems, Carlsbad, CA) and compared with previously determined gag-protease sequences from patients’ plasma viruses  by constructing a maximum-likelihood phylogenetic tree (PhyML) . Phylogenies were viewed using Figtree (http://tree.bio.ed.ac.uk/software/figtree).
Definition and analysis of HLA-associated polymorphisms in Gag-Protease
A total of 110 HLA-associated polymorphisms, occurring at 61 unique codons in Gag-Protease, were previously defined in a cohort of 430 treatment-naive Japanese subjects using phylogenetically-corrected methods . In the present study we wished to identify HLA-associated polymorphisms that modulate Gag-Pro RC. For this purpose, we compared Gag-Pro RCs in the presence or absence of each of these HLA-associated amino acid substitutions and assessed significance using the Mann–Whitney U test. Figures depicting amino acid variation at Gag-Pro-RC-modulating codons were created using WebLogo 3 (http://weblogo.threeplusone.com) [48, 49] to visualize sequence conservation at each position.
Site-directed mutagenesis and mutant virus stocks
Y79F, M228, and R286K substitutions were engineered into HIV-1 pNL4-3 by PCR overlap extension [50, 51], using PCR primers containing specific nucleotide changes: tat to ttc for Y79F, atg to ctg for M228L and aga to aag for R286K. These specific triplet codons were selected for each substitution based on higher frequencies in human codon usage, to exclude possible negative effects of rare codon usage in vitro. PCR fragments containing specific nucleotide changes were digested with BssHII (at nucleotide position 711 in pNL4-3) and SbfI (at 2838 in pNL4-3) followed by subcloning into pNL4-3 digested with the same restriction enzymes. The subcloned fragments were sequenced to confirm specific nucleotide changes and the absence of other mutations.
pNL4-3 carrying specific mutations were transfected into HEK293T cells by lipofectamine 2000 (Life Technologies, Carlsbad, CA) according to the manufacturer’s instructions. Virus stocks were harvested 48 h post-transfection, and titers were determined by infecting CEM-GXR cells and subsequently measuring GFP+ cells. Replication assays were performed in two sets of triplicate for each sample as described above.
Associations between Gag-Pro RC, clinical parameters and mutation scores were determined by Spearman’s rank correlation. Differences in pVL in context of the presence or absence of amino acid substitutions at four locations in Gag were determined by the Mann–Whitney U test. Multiple comparisons were addressed using q-values.
The relationship between Gag-Pro RC and pVL/CD4 count in the population subset that did not express protective HLA class I alleles (HLA-B*52:01−B*67:01− subset) was also investigated using multivariable linear regression. The primary model included only log10 plasma viral load (per log10 increment) and CD4 count (per 1 cell/mm3 increment) as variables. A second model additionally incorporated specific HLA alleles (C*08:01, C*08:02, C*12:02) that exhibited significant associations between Gag-Pro RC and both pVL and CD4 count as variables; here, model estimates are expressed in terms of individuals harboring (versus the reference group of individuals not harboring) the HLA allele in question. Multivariable analyses were performed using R (version 3.1.2).
The study was approved by the Ethical Committees in the Faculty of Life Science Kumamoto University (540) and the National Center for Global Health and Medicine. Informed consent was obtained from all patients, and the study was conducted in accordance with the Declaration of Helsinki.
human leukocyte antigen
plasma viral load
cytotoxic T lymphocytes
Kawashima Y, Pfafferott K, Frater J, Matthews P, Payne R, Addo M, Gatanaga H, Fujiwara M, Hachiya A, Koizumi H, Kuse N, Oka S, Duda A, Prendergast A, Crawford H, Leslie A, Brumme Z, Brumme C, Allen T, Brander C, Kaslow R, Tang J, Hunter E, Allen S, Mulenga J, Branch S, Roach T, John M, Mallal S, Ogwu A, Shapiro R, Prado JG, Fidler S, Weber J, Pybus OG, Klenerman P, Ndung’u T, Phillips R, Heckerman D, Harrigan PR, Walker BD, Takiguchi M, Goulder P. Adaptation of HIV-1 to human leukocyte antigen class I. Nature. 2009;2009(458):641–5.
Prince JL, Claiborne DT, Carlson JM, Schaefer M, Yu T, Lahki S, Prentice HA, Yue L, Vishwanathan SA, Kilembe W, Goepfert P, Price MA, Gilmour J, Mulenga J, Farmer P, Derdeyn CA, Tang J, Heckerman D, Kaslow RA, Allen SA, Hunter E. Role of transmitted Gag CTL polymorphisms in defining replicative capacity and early HIV-1 pathogenesis. PLoS Pathog. 2012;8:e1003041.
Frater AJ, Brown H, Oxenius A, Günthard HF, Hirschel B, Robinson N, Leslie AJ, Payne R, Crawford H, Prendergast A, Brander C, Kiepiela P, Walker BD, Goulder PJ, McLean A, Phillips RE. Effective T-cell responses select human immunodeficiency virus mutants and slow disease progression. J Virol. 2007;81:6742–51.
Matthews PC, Prendergast A, Leslie A, Crawford H, Payne R, Rousseau C, Rolland M, Honeyborne I, Carlson J, Kadie C, Brander C, Bishop K, Mlotshwa N, Mullins JI, Coovadia H, Ndung’u T, Walker BD, Heckerman D, Goulder PJ. Central role of reverting mutations in HLA associations with human immunodeficiency virus set point. J Virol. 2008;82:8548–59.
Martinez-Picado J, Prado JG, Fry EE, Pfafferott K, Leslie A, Chetty S, Thobakgale C, Honeyborne I, Crawford H, Matthews P, Pillay T, Rousseau C, Mullins JI, Brander C, Walker BD, Stuart DI, Kiepiela P, Goulder P. Fitness cost of escape mutations in p24 Gag in association with control of human immunodeficiency virus type 1. J Virol. 2006;80:3617–23.
Brumme ZL, Tao I, Szeto S, Brumme CJ, Carlson JM, Chan D, Kadie C, Frahm N, Brander C, Walker B, Heckerman D, Harrigan PR. Human leukocyte antigen-specific polymorphisms in HIV-1 Gag and their association with viral load in chronic untreated infection. AIDS. 2008;22:1277–86.
Schneidewind A, Brockman MA, Yang R, Adam RI, Li B, Le Gall S, Rinaldo CR, Craggs SL, Allgaier RL, Power KA, Kuntzen T, Tung CS, LaBute MX, Mueller SM, Harrer T, McMichael AJ, Goulder PJ, Aiken C, Brander C, Kelleher AD, Allen TM. Escape from the dominant HLA-B27-restricted cytotoxic T-lymphocyte response in Gag is associated with a dramatic reduction in human immunodeficiency virus type 1 replication. J Virol. 2007;81:12382–93.
Goepfert PA, Lumm W, Farmer P, Matthews P, Prendergast A, Carlson JM, Derdeyn CA, Tang J, Kaslow RA, Bansal A, Yusim K, Heckerman D, Mulenga J, Allen S, Goulder PJ, Hunter E. Transmission of HIV-1 Gag immune escape mutations is associated with reduced viral load in linked recipients. J Exp Med. 2008;205:1009–17.
Chopera DR, Woodman Z, Mlisana K, Mlotshwa M, Martin DP, Seoighe C, Treurnicht F, de Rosa DA, Hide W, Karim SA, Gray CM, Williamson C, CAPRISA 002 Study Team. Transmission of HIV-1 CTL escape variants provides HLA-mismatched recipients with a survival advantage. PLoS Pathog. 2008;4:e1000033.
Janes H, Friedrich DP, Krambrink A, Smith RJ, Kallas EG, Horton H, Casimiro DR, Carrington M, Geraghty DE, Gilbert PB, McElrath MJ, Frahm N. Vaccine-induced gag-specific T cells are associated with reduced viremia after HIV-1 infection. J Infect Dis. 2013;208:1231–9.
Ferrando-Martínez S, Casazza JP, Leal M, Machmach K, Muñoz-Fernández MÁ, Viciana P, Koup RA, Ruiz-Mateos E. Differential Gag-specific polyfunctional T cell maturation patterns in HIV-1 elite controllers. J Virol. 2012;86:3667–74.
Ahlers JD, Belyakov IM. Lessons learned from natural infection: focusing on the design of protective T cell vaccines for HIV/AIDS. Trends Immunol. 2010;3:120–30.
Goulder PJ, Watkins DI. Impact of MHC class I diversity on immune control of immunodeficiency virus replication. Nat Rev Immunol. 2008;8:619–30.
Huang KH, Goedhals D, Carlson JM, Brockman MA, Mishra S, Brumme ZL, Hickling S, Tang CS, Miura T, Seebregts C, Heckerman D, Ndung’u T, Walker B, Klenerman P, Steyn D, Goulder P, Phillips R, Bloemfontein-Oxford Collaborative G, vanVuuren C, Frater J. Progression to AIDS in South Africa is associated with both reverting and compensatory viral mutations. PloS One. 2011;6:e19018.
Brockman MA, Schneidewind A, Lahaie M, Schmidt A, Miura T, Desouza I, Ryvkin F, Derdeyn CA, Allen S, Hunter E, Mulenga J, Goepfert PA, Walker BD, Allen TM. Escape and compensation from early HLA-B57-mediated cytotoxic T-lymphocyte pressure on human immunodeficiency virus type 1 Gag alter capsid interactions with cyclophilin A. J Virol. 2007;81:12608–18.
Leslie AJ, Pfafferott KJ, Chetty P, Draenert R, Addo MM, Feeney M, Tang Y, Holmes EC, Allen T, Prado JG, Altfeld M, Brander C, Dixon C, Ramduth D, Jeena P, Thomas SA, St John A, Roach TA, Kupfer B, Luzzi G, Edwards A, Taylor G, Lyall H, Tudor-Williams G, Novelli V, Martinez-Picado J, Kiepiela P, Walker BD, Goulder PJ. HIV evolution: CTL escape mutation and reversion after transmission. Nat Med. 2004;10:282–9.
Boutwell CL, Rowley CF, Essex M. Reduced viral replication capacity of human immunodeficiency virus type 1 subtype C caused by cytotoxic-T-lymphocyte escape mutations in HLA-B57 epitopes of capsid protein. J Virol. 2009;83:2460–8.
Brockman MA, Brumme ZL, Brumme CJ, Miura T, Sela J, Rosato PC, Kadie CM, Carlson JM, Markle TJ, Streeck H, Kelleher AD, Markowitz M, Jessen H, Rosenberg E, Altfeld M, Harrigan PR, Heckerman D, Walker BD, Allen TM. Early selection in Gag by protective HLA alleles contributes to reduced HIV-1 replication capacity that may be largely compensated for in chronic infection. J Virol. 2010;84:11937–49.
Gijsbers EF, Feenstra KA, van Nuenen AC, Navis M, Heringa J, Schuitemaker H, Kootstra NA. HIV-1 replication fitness of HLA-B*57/58:01 CTL escape variants is restored by the accumulation of compensatory mutations in Gag. PLoS One. 2013;8:e81235.
Ammaranond P, van Bockel DJ, Petoumenos K, McMurchie M, Finlayson R, Middleton MG, Davenport MP, Venturi V, Suzuki K, Gelgor L, Kaldor JM, Cooper DA, Kelleher AD. HIV immune escape at an immunodominant epitope in HLA-B*27-positive individuals predicts viral load outcome. J. Immunol. 2011;186:479–88.
Schneidewind A, Brockman MA, Sidney J, Wang YE, Chen H, Suscovich TJ, Li B, Adam RI, Allgaier RL, Mothé BR, Kuntzen T, Oniangue-Ndza C, Trocha A, Yu XG, Brander C, Sette A, Walker BD, Allen TM. Structural and functional constraints limit options for cytotoxic T-lymphocyte escape in the immunodominant HLA-B27-restricted epitope in human immunodeficiency virus type 1 capsid. J Virol. 2008;82:5594–605.
Kaslow RA, Carrington M, Apple R, Park L, Muñoz A, Saah AJ, Goedert JJ, Winkler C, O’Brien SJ, Rinaldo C, Detels R, Blattner W, Phair J, Erlich H, Mann DL. Influence of combinations of human major histocompatibility complex genes on the course of HIV-1 infection. Nat Med. 1996;2:405–11.
Goulder PJ, Phillips RE, Colbert RA, McAdam S, Ogg G, Nowak MA, Giangrande P, Luzzi G, Morgan B, Edwards A, McMichael AJ, Rowland-Jones S. Late escape from an immunodominant cytotoxic T-lymphocyte response associated with progression to AIDS. Nat Med. 1997;3:212–7.
Kelleher AD, Long C, Holmes EC, Allen RL, Wilson J, Conlon C, Workman C, Shaunak S, Olson K, Goulder P, Brander C, Ogg G, Sullivan JS, Dyer W, Jones I, McMichael AJ, Rowland-Jones S, Phillips RE. Clustered mutations in HIV-1 gag are consistently required for escape from HLA-B27-restricted cytotoxic T lymphocyte responses. J Exp Med. 2001;193:375–86.
Miura T, Brockman MA, Brumme ZL, Brumme CJ, Pereyra F, Trocha A, Block BL, Schneidewind A, Allen TM, Heckerman D, Walker BD. HLA-associated alterations in replication capacity of chimeric NL4-3 viruses carrying gag-protease from elite controllers of human immunodeficiency virus type 1. J Virol. 2009;83:140–9.
Miura T, Brumme ZL, Brockman MA, Rosato P, Sela J, Brumme CJ, Pereyra F, Kaufmann DE, Trocha A, Block BL, Daar ES, Connick E, Jessen H, Kelleher AD, Rosenberg E, Markowitz M, Schafer K, Vaida F, Iwamoto A, Little S, Walker BD. Impaired replication capacity of acute/early viruses in persons who become HIV controllers. J Virol. 2010;84:7581–91.
Wright JK, Naidoo VL, Brumme ZL, Prince JL, Claiborne DT, Goulder PJ, Brockman MA, Hunter E, Ndung’u T. Impact of HLA-B*81-associated mutations in HIV-1 Gag on viral replication capacity. J Virol. 2012;86:3193–9.
Wright JK, Novitsky V, Brockman MA, Brumme ZL, Brumme CJ, Carlson JM, Heckerman D, Wang B, Losina E, Leshwedi M, van der Stok M, Maphumulo L, Mkhwanazi N, Chonco F, Goulder PJ, Essex M, Walker BD, Ndung’u T. Influence of Gag-protease-mediated replication capacity on disease progression in individuals recently infected with HIV-1 subtype C. J Virol. 2011;85:3996–4006.
Wright JK, Brumme ZL, Carlson JM, Heckerman D, Kadie CM, Brumme CJ, Wang B, Losina E, Miura T, Chonco F, van der Stok M, Mncube Z, Bishop K, Goulder PJ, Walker BD, Brockman MA, Ndung’u T. Gag-protease-mediated replication capacity in HIV-1 subtype C chronic infection: associations with HLA type and clinical parameters. J Virol. 2010;84:10820–31.
Juarez-Molina CI, Payne R, Soto-Nava M, Avila-Rios S, Valenzuela-Ponce H, Adland E, Leitman E, Brener J, Muenchhoff M, Branch S, Landis C, Reyes-Teran G, Goulder P. Impact of HLA selection pressure on HIV fitness at a population level in Mexico and Barbados. J Virol. 2014;88:10392–8.
Chikata T, Carlson JM, Tamura Y, Borghan MA, Naruto T, Hashimoto M, Murakoshi H, Le AQ, Mallal S, John M, Gatanaga H, Oka S, Brumme ZL, Takiguchi M. Host-specific adaptation of HIV-1 subtype B in the Japanese population. J Virol. 2014;88:4764–75.
Naruto T, Gatanaga H, Nelson G, Sakai K, Carrington M, Oka S, Takiguchi M. 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:10870–2.
Hashimoto M, Kinoshita T, Yamasaki M, Tanaka H, Imanishi T, Ihara H, Ichikawa Y, Fukunishi T. Gene frequencies and haplotypic associations within the HLA region in 916 unrelated Japanese individuals. Tissue Antigens. 1994;44:166–73.
Saito S, Ota S, Yamada E, Inoko H, 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:522–9.
Tanaka H, Akaza T, Juji T. Report of the Japanese Central Bone Marrow Data Center. Clin Transpl. 1996;10:139–44.
Tokunaga K, Ishikawa Y, Ogawa A, Wang H, Mitsunaga S, Moriyama S, Lin L, Bannai M, Watanabe Y, Kashiwase K, Tanaka H, Akaza T, Tadokoro K, Juji T. Sequence-based association analysis of HLA class I and II alleles in Japanese supports conservation of common haplotypes. Immunogenetics. 1997;46:199–205.
Itoh Y, Mizuki N, Shimada T, Azuma F, Itakura M, Kashiwase K, Kikkawa E, Kulski JK, Satake M, Inoko H. High-throughput DNA typing of HLA-A, -B, -C, and -DRB1 loci by a PCR-SSOP-Luminex method in the Japanese population. Immunogenetics. 2005;57:717–29.
Brockman MA, Tanzi GO, Walker BD, Allen TM. Use of a novel GFP reporter cell line to examine replication capacity of CXCR4- and CCR5-tropic HIV-1 by flow cytometry. J Virol Methods. 2006;131:134–42.
Murakoshi H, Akahoshi T, Koyanagi M, Chikata T, Naruto T, Maruyama R, Tamura Y, Ishizuka N, Gataganaga H, Oka S, Takiguchi M. Clinical control of HIV-1 by cytotoxic T cells specific for multiple conserved epitopes. J Virol. 2015;89:5330–9.
Nomura S, Hosoya N, Brumme ZL, Brockman MA, Kikuchi T, Koga M, Nakamura H, Koibuchi T, Fujii T, Carlson JM, Heckerman D, Kawana-Tachikawa A, Iwamoto A, Miura T. Significant reductions in Gag-protease-mediated HIV-1 replication capacity during the course of the epidemic in Japan. J Virol. 2013;87:1465–76.
Brumme ZL, Li C, Miura T, Sela J, Rosato PC, Brumme CJ, Markle TJ, Martin E, Block BL, Trocha A, Kadie CM, Allen TM, Pereyra F, Heckerman D, Walker BD, Brockman MA. Reduced replication capacity of NL4-3 recombinant viruses encoding reverse transcriptase-integrase sequences from HIV-1 elite controllers. J Acquir Immune Defic Syndr. 2011;56:100–8.
Brockman MA, Chopera DR, Olvera A, Brumme CJ, Sela J, Markle TJ, Martin E, Carlson JM, Le AQ, McGovern R, Cheung PK, Kelleher AD, Jessen H, Markowitz M, Rosenberg E, Frahm N, Sanchez J, Mallal S, John M, Harrigan PR, Heckerman D, Brander C, Walker BD, Brumme ZL. Uncommon pathways of immune escape attenuate HIV-1 integrase replication capacity. J Virol. 2012;86:6913–23.
Barbour JD, Hecht FM, Wrin T, Segal MR, Ramstead CA, Liegler TJ, Busch MP, Petropoulos CJ, Hellmann NS, Kahn JO, Grant RM. Higher CD4+ T cell counts associated with low viral pol replication capacity among treatment-naive adults in early HIV-1 infection. J Infect Dis. 2004;190:251–6.
Daar ES, Kesler KL, Wrin T, Petropoulo CJ, Bates M, Lail A, Hellmann NS, Gomperts E, Donfield S, Hemophilia Growth and Development Study. HIV-1 pol replication capacity predicts disease progression. AIDS. 2005;19:871–7.
Goetz MB, Leduc R, Wyman N, Kostman JR, Labriola AM, Lie Y, Weidler J, Coakley E, Bates M, Luskin-Hawk R, Long Term Monitoring Study (CPCRA060). Terry Beirn Community Programs for Clinical Research on AIDS. HIV replication capacity is an independent predictor of disease progression in persons with untreated chronic HIV infection. J Acquir Immune Defic Syndr. 2010;53:472–9.
Ng OT, Laeyendecker O, Redd AD, Munshaw S, Grabowski MK, Paquet AC, Evans MC, Haddad M, Huang W, Robb ML, Reynolds SJ, Gray RH, Wawer MJ, Serwadda D, Eshleman SH, Quinn TC. HIV type 1 polymerase gene polymorphisms are associated with phenotypic differences in replication capacity and disease progression. J Infect Dis. 2014;209:66–73.
Guindon S, Gascuel O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003;52:696–704.
Crooks GE, Hon G, Chandonia JM, Brenner SE. WebLogo: a sequence logo generator. Genome Res. 2004;14:1188–90.
Schneider TD, Stephens RM. Sequence logos: a new way to display consensus sequences. Nucleic Acids Res. 1990;18:6097–100.
Ho SN, Hunt HD, Horton RM, Pullen JK, Pease LR. Site-directed mutagenesis by overlap extension using the polymerase chain reaction. Gene. 1989;77:51–9.
Heckman KL, Pease LR. Gene splicing and mutagenesis by PCR-driven overlap extension. Nat Protoc. 2007;2:924–32.
KS designed, acquired, analyzed data and drafted the manuscript. TM made substantial contributions to conception, design, and analysis of data and drafted the manuscript. TC and CJB performed the statistical analysis and contributed to drafting the manuscript. ZL generated and analyzed data and revised the manuscript critically for important intellectual content. HG and SO generated clinical data and samples. All authors read and approved the final manuscript.
We thank Drs. Richard Harrigan and Bruce Walker for their cohort data; Drs. Toshiyuki Miura and Mark Brockman for generously providing CEM-GXR cells and/or RC data; Mari Hasegawa, Yoshiki Eto, Yoshiko Tamura, and Lin Zhansong for technical assistance, and Sachiko Sakai for secretarial assistance.
T.C. is a JSPS Research Fellow. Z.L.B. is the recipient of a New Investigator Award from the Canadian Institutes for Health Research and a Scholar Award from the Michael Smith Foundation for Health Research (MSFHR). This study was supported by the Global COE program Global Education and Research Center Aiming at the Control of AIDS, launched as a project commissioned by the Ministry of Education, Science, Sports, and Culture, Japan (K.S. and M.T.), Grant-in-Aid for Scientific Research (B) Grant Number 26293240 (M.T.), and Grants-in-Aid for AIDS Research from the Ministry of Health, Labor, and Welfare, Japan (M.T.).
The authors declare that they have no competing interests.
An erratum to this article is available at http://dx.doi.org/10.1186/s12977-015-0227-8.