Skip to main content

Rapid HIV disease progression following superinfection in an HLA-B*27:05/B*57:01-positive transmission recipient

Abstract

Background

The factors determining differential HIV disease outcome among individuals expressing protective HLA alleles such as HLA-B*27:05 and HLA-B*57:01 remain unknown. We here analyse two HIV-infected subjects expressing both HLA-B*27:05 and HLA-B*57:01. One subject maintained low-to-undetectable viral loads for more than a decade of follow up. The other progressed to AIDS in < 3 years.

Results

The rapid progressor was the recipient within a known transmission pair, enabling virus sequences to be tracked from transmission. Progression was associated with a 12% Gag sequence change and 26% Nef sequence change at the amino acid level within 2 years. Although next generation sequencing from early timepoints indicated that multiple CD8+ cytotoxic T lymphocyte (CTL) escape mutants were being selected prior to superinfection, < 4% of the amino acid changes arising from superinfection could be ascribed to CTL escape. Analysis of an HLA-B*27:05/B*57:01 non-progressor, in contrast, demonstrated minimal virus sequence diversification (1.1% Gag amino acid sequence change over 10 years), and dominant HIV-specific CTL responses previously shown to be effective in control of viraemia were maintained. Clonal sequencing demonstrated that escape variants were generated within the non-progressor, but in many cases were not selected. In the rapid progressor, progression occurred despite substantial reductions in viral replicative capacity (VRC), and non-progression in the elite controller despite relatively high VRC.

Conclusions

These data are consistent with previous studies demonstrating rapid progression in association with superinfection and that rapid disease progression can occur despite the relatively the low VRC that is typically observed in the setting of multiple CTL escape mutants.

Background

Outcome in HIV infection is strongly influenced by the particular HLA alleles that are expressed, the most protective of which are HLA-B*27:05 and HLA-B*57:01 [1, 2]. The mechanism by which these HLA alleles confer protection is related to the presentation of a broad array of epitopes that enable CD8+ cytotoxic T lymphocytes (CTL) to recognise and kill HIV-infected target cells [2, 3]. As HIV-specific CTL are likely to play a critical role in strategies designed to eradicate viral reservoirs [4], it remains important to define what features affect the ability of CTL to be effective in killing virus-infected cells. We here contrast HIV disease outcomes in two HIV-infected subjects expressing both HLA-B*27:05 and HLA-B*57:01. One subject is a typical HIV non-progressor, maintaining viral loads below the level of detection for more than a decade of follow up. The other subject initially presented the characteristic phenotype of an HLA-B*27:05/57:01-positive HIV-infected subject, with a low viral load (65 copies/ml) and high absolute CD4+ T cell count (950 cells/mm3), but within 2.7 years had progressed rapidly to AIDS (CD4+ T cell count < 200 cells/mm3).

Methods

Study subjects

The adult Caucasian transmission pair R096 (HLA-A*02:01/A*11:01 B*13:02/B*51:01 C*03:03 C*06:02) and R097 (HLA-A*11:01/A*24:01 B*27:05/B*57:01 C*02:02 C*06:02) and non-progressor RI088 (HLA-A*01:01/A*03:01 B*27:05/B*57:01 C*02:02 C*06:02) were recruited from the Thames Valley Cohort, UK, previously described [5].

DNA and RNA extraction

Proviral and genomic DNA were extracted for proviral sequencing and HLA-typing respectively, using PureGene reagents according to the manufacturers instructions. Viral plasma RNA was extracted using the Qiamp Viral RNA Mini Kit (Qiagen, UK) according to the manufacturers instructions. For RNA extraction for ultra-deep sequencing, minor modifications to the manufacturers protocol were made, as previously described [6].

HLA-typing

Four digit HLA Class I typing was performed at the CLIA/ASHI accredited laboratory of William Hildebrand, PhD, D (ABHI), University of Oklahoma Health Sciences Center. Locus specific PCR amplification of exons 2 and 3 and heterozygous sequencing were performed from genomic DNA. Ambiguities were resolved by homozygous sequencing [7].

IFN-γ ELISpot assays

Peripheral blood mononuclear cells were tested for responses to 410 18mer overlapping peptides spanning the B clade HIV proteome as previously described [8] using pools of 11–12 peptides followed by confirmatory assays using individual 18mer peptides.

Sanger sequencing of proviral DNA

The gag gene was amplified from proviral DNA using a nested touchdown PCR with BioTaq DNA polymerase (Bioline, UK) using the following primers: 5′-CTCTAGCAGTGGCGCCCGAA-3′ and 5′-TCCTTTCCACATTTCCAACAGCC-3′ for the first round PCR and 5′-ACTCGGCTTGCTGAAGTGC-3′ and 5′ CAATTTCTGGCTATGTGCCC-3′ for the second round PCR. Twenty cycles of denaturation at 94 °C for 15 s, annealing at 60 °C for 30 s and elongation at 72 °C for 1 min were performed followed by another 20 cycles with an annealing temperature of 57 °C. Purified PCR product was used to prepare sequencing templates using BigDye Terminator v3.1 reaction mix (Applied Biosystems, UK) and sequenced on an ABI 3730xl DNA Analyzer (Applied Biosystems, UK) by the Department of Zoology Sequencing Facility, University of Oxford. For clonal sequencing, purified PCR product was cloned into TOPO vectors using Zero Blunt TOPO PCR Cloning Kit (Invitrogen, UK) and used to transform chemically competent One Shot TOP10 E. coli cells (Invitrogen, UK) according to the manufacturers instructions. Where possible, clones were selected from two independent PCRs to avoid PCR amplification bias. Selected colonies were cultured overnight in LB broth before proceeding with mini-prep plasmid DNA extraction using Montage 96-well plasmid preparation kits (Millipore, US) according to the manufacturers instructions. Plasmid DNA was used for sequencing template preparation as described above.

Ultra-deep sequencing de novo assembly of consensus sequences and minor variant haplotype analysis

The full-length HIV genome was amplified in four fragments from plasma RNA using Superscript III One-Step RT PCR Kit with Platinum Taq High Fidelity enzyme (Invitrogen, UK) as previously described [9]. Sequencing of pooled amplicons was performed using Illumina MiSeq 250 bp paired-end technology. Quality control of reads was performed using QUASR (http://sourceforge.net/projects/quasr/) as previously described [6, 10, 11]. A de novo assembly was constructed using SPAdes version 2.4.0 [12] and a consensus sequence was generated using Abacas version 1.3.1 and MUMmer version 3.2 [13]. Haplotypes in the epitope regions were determined using Vprofiler [12] by selecting reads that span the epitope region and which contain only accepted variants.

Phylogenetic analysis and recombination detection

Maximum likelihood phylogenetic trees were constructed using Mega 6.06 software under the General Time Reversible model of nucleotide substitution as determined by jModelTest version 0.1.1 [14] with 1000 bootstrap replicates and viewed using FigTree v1.4.0 software. Recombination analysis was performed using RDP 4.46 (Recombination Detection Programme) [15] and SimPlot 3.5.1 [16].

Generation of recombinant gag-protease HIV virions and viral replicative capacity assays (VRC)

Recombinant viral vectors expressing autologous-gag-protease sequences from HLA-B*27:05/HLA-B*57:01+ individuals were constructed by transfecting gag-protease amplicons into the Gag-Pro deleted HIV backbone ΔGag-Pro NL4-3 as previously described [17, 18]. Gag-protease was amplified by nested PCR from plasma RNA using previously published primers [18]. The first round amplification was performed using Superscript III One-Step RT PCR Kit with Platinum Taq High Fidelity enzyme (Invitrogen, UK) followed by second round PCR using High Fidelity Platinum (Invitrogen, UK) according to the manufacturers’ instructions. PCR product was purified using QIAquick PCR Purification Kit (Qiagen, UK) and sequenced by Sanger sequencing to confirm accurate amplification of the patient-derived gag-protease sequence.

The ΔGag-Pro NL4-3 was produced by introducing BstEII restriction enzyme sites on either side of the gag-protease-gene in NL4-3, followed by BstEII digest to delete gag-protease and allow self ligation of the plasmid [18]. For transfection, the plasmid was linearised by BstEII digest at 60 °C for 2 h. For each transfection, 2.5 μg of purified gag-protease amplicon and 10 μg of linearised ΔGag-Pro NL4-3 were co-transfected into 2 million CEM-GXR reporter cells [17, 19] by electroporation (300 V, 500 μF capacitance and infinite resistance). Recombinant viral cultures were expanded over 1 month and viral spread measure by GFP expression detected by flow cytometry. Viral supernatant was collected when > 30% of live cells were infected. Viral RNA extraction and gag-protease sequencing confirmed the correct sequence amplification of viral stocks. Viral titres were determined by infection of 1 million CEM-GXR cells with a standardized volume of viral supernatant and examining GFP expression after 48 h.

Viral replicative capacity assays (VRC): CEM-GXR cells were infected in triplicate at a low multiplicity of infection (0.025% GFP+ cells at day 2). GFP expression was monitored daily by flow cytometry over 12 days.

Logistic curve modelling of viral replicative capacity data

Viral growth curves were modeled on the logistic curve function:

$$\begin{aligned} y & = m\left( {x,\theta } \right) + \varepsilon \\ & = \theta _{1} /1 + \exp \left( { - \theta _{2} + \theta _{3} x} \right) \\ \end{aligned}$$

where θ1 represents the asymptotic y point, θ2 the midpoint, and θ3 the scaling factor (rate of change) of the curve. Curves were fitted using the nonlinear regression nls function in R version 3.2.1, plotting the best fit logistic growth curve. To assess whether each parameter was statistically different when comparing curves, models were built where one of the parameters was fixed while the others were allowed to vary. A model where all parameters were allowed to vary was then fitted and a statistical comparisons of parameters was made using ANOVA. The parameters used for the curves in Fig. 4 are shown in Additional file 1: Table 1, together with the statistics for the overall fit and for the separate parameters.

Sequence accession numbers

Sequence data has been deposited in Genbank under the accession numbers MF039091-MF039203.

Results

Rapid progression to AIDS in a transmission recipient expressing HLA-B*27:05/B*57:01

The starting point for this study was the MSM transmission pair R096/R097 (Fig. 1a, b). The transmission recipient, R097, expressed both HLA-B*27:05 and HLA-B*57:01, both of which are strongly associated with slow progression to AIDS. The viral load in R097 at diagnosis (time 0) was 65 HIV RNA copies/ml plasma and the absolute CD4+ T cell count was 950 cells/mm3, reflecting successful immune control of HIV. Over the next 2–3 years, however, despite the co-expression of HLA-B*27:05 and HLA-B*57:01, progression to AIDS (absolute CD4+ T cell < 200 cells/mm3) occurred rapidly, in association with viral loads increasing to > 105 copies/ml. The donor R096 at time 0 had a viral load of 221,839 copies/ml and an absolute CD4+ T cell count of 260 cells/mm3 and initiated antiretroviral therapy (ART) at this time.

Fig. 1
figure 1

Plasma HIV RNA viral load and CD4+ T cell counts for HLA-B*27:05/57:01+ subject R097 and his transmission partner R096 and phylogenetic analysis of HIV sequences showing epidemiological linkage and sequence diversity in transmission pair R096/R097. a Transmission pair donor, R096. b Transmission pair recipient, R097. ‘Early’ and ‘late’ time points of sampling for sequencing are shown (see c). ‘Time 0’ represents the time of diagnosis of recipient R097. The horizontal dotted line represents the limit of detection (LOD) of the viral load assay (40 HIV RNA copies/ml). Grey shading indicates the period during which the subject received antiretroviral (ART) therapy. c Maximum likelihood phylogenetic tree of 1091 bp alignment of RNA (ultra-deep consensus) sequences across the Gag p17 and p24 genes. The donor R096 sequence is shown in blue. Early and late sequencing time points are shown in a. The early R097 sequences are shown in red (n = 2). The late (post-2 year) sequences for R097 are shown in green (n = 11). The early and late sequences form two distinct clusters indicating considerable intra-host sequence diversification over time. 143 B clade reference sequences from the US and UK collected between 2003 and 2011 from the Los Alamos database (http://www.hiv.lanl.gov/) are shown in black. Bootstrap values for the R097 sequence clusters based on 1000 bootstrap replicates are shown in italics

To investigate why disease had progressed so rapidly in R097, and to seek evidence to support the clinical history suggesting transmission from R096 to R097, we ultra-deep sequenced the virus in both individuals (Fig. 1c). Tight phylogenetic clustering of the viral sequences at ‘time 0’ in R096 and R097 strongly supported the notion that R096/R097 were a transmission pair. Viral sequences from R097 post-progression (2.0–3.4 years after ‘time 0’) clustered separately from the time 0 sequences. Gag sequences in R097 differed by 12% at the amino acid level in the 2 years from ‘time 0’ and Nef sequences by 26% at the amino acid level over this short time. These changes occurred in the context of a 12.5% nucleotide change across the full-length genome.

The direction of transmission from R096 to R097 was supported by the presence of virus encoding escape mutants in five HLA-B*51:01-restricted CTL epitopes in both subjects (Table 1). Since HLA-B*51:01 is expressed by R096 and not R097, this is consistent with the variant being selected initially in R096 and being transmitted to R097. In four of the HLA-B*51:01-restricted epitopes, the variant we observed was the one most commonly seen in association with HLA-B*51:01, a so-called HLA footprint [for example, I135T within the epitope TAFTIPSI (RT 128-135)] [20,21,22]. In contrast, no HLA-B*57:01 or HLA-B*27:05-associated mutants were observed in the ‘time 0’ sequences in R096 and R097 (Table 2). One variant L268I was shared by R096 and R097 at ‘time 0’, which lies within the HLA-B*27:05-restricted epitope KRWIILGLNK (Gag 263-272). However, this mutant is an HLA-B*08:01-associated variant flanking the overlapping HLA-B*08:01 epitope EIYKRWII (Gag 260-267) [22] and thus is likely to have been originally selected in an HLA-B*08:01-positive subject prior to transmission to R096.

Table 1 HLA-B*51:01 and HLA-B*13:02 footprints in the recipient R097 sequence
Table 2 Donor R096 and recipient R097 viral sequences encoding HLA-B*27:05/B*57:01-restricted epitopes

Mechanism of distinct phylogenetic clustering of R097 sequences following progression

Comparison of the viral sequences in R097 at ‘time 0’ and the cluster following progression 2.0–3.4 years later showed dramatic differences, prompting the hypothesis that R097 might have progressed as a consequence of HIV superinfection. To determine whether the dramatic sequence changes observed were likely to be the result of superinfection, or of multiple escape mutations being selected over a short period of time to give the appearance of superinfection, we first sought evidence of escape within HLA-B*27:05 and HLA-B*57:01-restricted epitopes selected between ‘time 0’ and 2.0 years later (Table 2). In patient R097 between ‘time 0’ and 2.0 years later, we examined 17 HLA-B*27:05/57:01 epitopes previously shown to contain well-characterised footprints [22], and found escape mutants in 12 of these (one mutation being contained in two overlapping HLA-B*27:05 and HLA-B*57:01 epitopes in Vpr). Additional mutants were observed that were selected between ‘time 0’ and 2.0 years later in HLA-B*27:05/57:01 restricted epitopes, but where HLA footprints have not been described (for example in the HLA-B*27:05 epitope GRRGWEALKY). Thus a proportion of the amino acid differences between the ‘time 0’ sequence and the viral sequences post progression can be explained by the selection of escape mutants in R097. However this did not distinguish rapid selection of escape mutants from the possibility of superinfection with an HIV variant carrying multiple HLA-B*27:05 and HLA-B*57:01-associated escape variants.

To examine this further, we analysed the timing of selection of escape mutants in R097 using ultra-deep sequencing (Table 3). At ‘time 0’ in R097, many of the escape variants that had reached fixation 2.0 years later were already present at low frequencies. For example, Gag-A146P was present at 19%, Gag-T242N at 4%, Gag-R264K at 3%, and Pol-T839N at 50% in R097 at ‘time 0’, and all were present at 99 or 100% 2.0 years later. Each of these variants were present in 0% of viruses in R096 at ‘time 0’. These data are consistent with the rapid increase in frequency of escape mutations in multiple epitopes restricted by HLA-B*27:05 and HLA-B*57:01 in the 2 years after ‘time 0’, but clearly selection of escape mutants had already been initiated by ‘time 0’.

Table 3 Ultra-deep RNA sequencing of HLA-B*57:01 and HLA-B*27:05-restricted epitopes and associated compensatory positions from R097 and his transmission partner R096

However, two features of the sequence data argue strongly in favour of superinfection having occurred in R097. First, some of the escape mutants selected at ‘time 0’ and not shared at the ‘time 2.0’ timepoint, and vice versa, are relatively rare. Examples here are A248E (87% at ‘time 0’ and 0% at ‘time 2.0’), which is a highly effective escape mutant [18], and S148T (0% at ‘time 0’ and 100% at ‘time 2.0’). Both of these mutants are within HLA-B*57:01-restricted CTL epitopes, but they are not selected sufficiently often to be included as HLA-B*57:01-associated ‘footprints’ [22]. In contrast, many of the HLA-B*27:05/B*57:01-associated escape mutants present both at ‘time 0’ and 0% at ‘time 2.0’, such as T242N and R264K, are commonly observed footprints for these HLA alleles [22], and could easily have arisen independently in the two viruses in different subjects expressing HLA-B*27:05/57:01. Second, the number of amino acid differences between the viruses at ‘time 0’ and at ‘time 2.0’ across the full-length proteomes numbered 810. Less than 3% of these amino acid changes (19/810) can be ascribed to escape arising within epitopes that are restricted by the HLA alleles expressed by R097. Thus the dramatic sequence changes observed in R097 between ‘time 0’ and 0% at ‘time 2.0’ cannot have been driven by CTL escape alone, but are consistent with HIV superinfection.

HCV co-infection as an additional explanation for rapid immunological escape

Having established that HIV superinfection was likely to account for the rapid evolutionary changes in HIV sequences and disease progression in subject R097, we considered alternative explanations for the uncharacteristic speed of progression in an HLA B*27/B*57-positive individual. At the time of peak HIV viraemia (time 2.0 years) the recipient was diagnosed with hepatitis C virus (HCV) genotype 1a co-infection by PCR (viral load 4.03 log copies/ml), although his HCV antibody test remained negative for a further 2.5 years before he finally seroconverted. An HCV RNA test 10 months earlier (that is, 14 months after ‘time 0’) had been negative. He went on to receive treatment with IFNα/Ribavirin for 24 weeks, starting 5 months after the HCV diagnosis was made, although this did not result in clearance. Although HCV coinfection typically accelerates disease progression in HIV infection [23], in subject R097, as described above, escape mutants in multiple HIV-specific HLA-B*27:05- and HLA-B*57:01-restricted epitopes had already been selected, and the viral load was on an upward trajectory, more than a year prior to the HCV infection. Nonetheless, it is evident that HIV superinfection occurred either simultaneous with, or weeks-to-months before, HCV co-infection. In either case, HCV co-infected is likely to have accelerated the speed of progression to AIDS in this subject.

Lack of sequence diversification in an HLA-B*27:05/57:01-positive non-progressor

To contrast the speed and extent of viral escape in the HLA-B*27:05/57:01-positive rapid progressor, R097, even before superinfection had occurred, with the sequence changes that might be expected in a non-progressor, we analysed gag viral sequences over a 10-year period in an ART-naive HLA-B*27:05/57:01-positive individual RI088, whose viral loads remained at < 1000 copies/ml, and at < 40 copies/ml for much of this period (Fig. 2). In this instance, viral sequence diversification was negligible with 1.1% difference in Gag at the amino acid level over 10 years. In contrast, we saw a 12% difference at the amino acid level in Gag over 2 years in the rapid progressor R097. Analysis of viral sequence at the clonal level in RI088 (Table 4, Additional file 2: Figure 1) shows that escape mutants were generated periodically but not always selected. For example, the most commonly selected variant, T242N, within the HLA-B*57:01 epitope TSTLQEQIGW (Gag 240–249), reached a frequency of 55% 1.4 years after diagnosis, but subsequently declined to 0% (at 3.4 years) before being selected later. Variants within KF11 and KK10 were observed transiently and at low frequencies. The relative lack of sequence change in the virus in RI088 is also reflected in persistence of CTL responses, such as to the HLA-B*27:05-restricted Gag epitope KK10 and the HLA-B*57:01-restricted Nef epitope HW9 (Fig. 3). In contrast, although these same responses are present at ‘time 0’ in R097, they declined in association with the selection of escape (Fig. 3).

Fig. 2
figure 2

Clinical course of infection and phylogenetic analysis of viral sequences from HLA-B*27:05/57:01-positive controller RI088. a Plasma HIV viral load and CD4+ T cell count for HLA-B*27:05/B*57:01-positive recipient RI088. ‘Time 0’ represents the time of diagnosis. The horizontal dotted line represents the limit of detection (LOD) of the viral load assay (40 copies/ml). b Maximum likelihood phylogenetic tree of a 1089 bp alignment of the Gag gene sequenced from proviral DNA (consensus of clonal sequences shown in Table 4). Subject RI088 (n = 6) sequences spanning 10 years are shown in red. 180 B clade reference sequences from the US and UK collected between 2003 and 2011 from the Los Alamos database (http://www.hiv.lanl.gov/) are shown in black. Bootstrap values for the RI088 sequence cluster based on 1000 bootstrap replicates are shown in italics

Table 4 Clonal Sequencing of HIV proviral DNA at HLA‐B*57:01 and HLA-B*27:05-restricted Gag epitopes from controller RI088
Fig. 3
figure 3

HIV-specific CTL responses in progressor R097 and in non-progressor RI088. EliSpot assays tested recognition of overlapping 18mer peptides, which together spanned the B clade proteome. The figure includes the full sequence of each 18mer peptide recognised and the HLA restriction where known, with known optimal epitopes shown in bold. a EliSpot responses detected in progressor R097 at time 0 (diagnosis) and 2.0 years later. Data shown at 2.0 years are the median of 3 assays undertaken at 1.96, 1.98 and 2.15 years after time 0. No responses were detected to the 18mer peptides containing the epitopes ISPRTLNAW or STTVKAACWW at either timepoint. b Responses detected in non-progressor RI088 at the times shown post diagnosis. Selected peptides that were strongly recognized in R097 but not in RI088 are included in panel B

Progression to AIDS in subject R097 associated with increasing viral replicative capacity

One of the mechanisms that has been proposed for HLA-B*27:05 and HLA-B*57:01-mediated protection against progression to HIV disease is the observation that the Gag escape mutations associated with these HLA molecules significantly reduce viral replicative capacity [5, 24,25,26,27]. Thus, even if the virus succeeds in evading the effective Gag-specific CTL responses restricted by HLA-B*27:05/57:01, the cost of reduced VRC facilitates maintained control of HIV through other immune responses from which the virus has not escaped. However, it is also well-described that reductions in viral replicative capacity resulting from selection of certain escape mutants are typically mitigated, to a varying degree, by selection of alternative escape mutants [28] and the co-selection of compensatory mutants [27, 29,30,31]. To investigate whether VRC was substantially reduced in the face of the selection of, or superinfection by a virus carrying, multiple Gag escape mutants in R097, we generated chimeric viruses comprising NL4-3 and the autologous gag-pro sequence from R097 (Fig. 4a). At ‘time 0’, when the viral load was low and the CD4+ T cell count was high, the VRC in R097 was substantially lower than both NL4-3 wildtype and ‘time 0’ R096 virus. However, during the course of progression, the VRC ultimately had increased to an intermediate level (at 3.4 years after ‘time 0’), reflecting some successful correction of fitness costs of escape variants by compensatory mutants or by the selection of alternative, higher fitness escape mutants. In contrast, although in the case of RI088 the transmitted virus was not known, the lack of multiple Gag escape mutants in RI088 appeared to have little impact on VRC compared to wildtype NL4-3 (Fig. 4b). These data are consistent with previous findings that the HLA-B*27:05- and HLA-B*57:01-associated escape mutants in Gag reduce viral replicative capacity, but even a relatively disabled virus can drive rapid progression to AIDS in the absence of effective HIV-specific CD8+ T cell activity.

Fig. 4
figure 4

Viral replicative capacity (VRC) of recombinant viruses produced from autologous gag-pro from HLA-B*27:05/B*57:01-positive subjects. VRC of gag-pro chimeric virus derived from a longitudinal sampling of progressor R097 and his transmission partner R096; and b controller RI088. VRC is given by the percentage of infected (GFP-expressing) target cells

Discussion

We here describe two HIV-infected adults who express both HLA-B*27:05 and HLA-B*57:01, the most protective HLA class molecules against HIV disease progression, one of whom is a typical long-term non-progressor, maintaining low-to-undetectable viral loads over more than a decade of follow up, whilst the other is a rapid progressor, whose CD4+ T cell decline from immune control of HIV to AIDS occurs within a period of 3 years. In the rapid progressor studied here, we have shown using ultra-deep sequencing that escape mutations were selected in early infection across several of the known HLA-B*27:05 and HLA-B*57:01 epitopes, from the earliest timepoint when the CD4+ T cell count was 950 cells/mm3 and viral load was 65 copies/ml, preceding loss of immune control. However, the most dramatic viral sequence changes observed in this subject, coincident with rapid progression to AIDS, were associated with superinfection with a virus whose ancestry is substantially different from the original transmitted strain of HIV. Although many of the classical HLA-B*27:05 and HLA-B*57:01 ‘footprints’ [22] were present both pre- and post-superinfection, these commonly arise predictably and independently in individuals expressing these HLA molecules, and > 96% of amino acid changes were unrelated to CTL-driven escape.

In contrast, very little viral sequence change occurred over more than a decade in the non-progressor. Clonal sequencing reveals that escape mutants were generated all of this time, but in the majority of instances these were not selected. Finally, analysis of the viral replicative capacity (VRC) of the viruses in the two subjects revealed that differential rates of progression were not explained by lower VRC in the non-progressor. On the contrary, rapid progression occurred despite reductions in VRC in the subject with progressive disease.

Approximately 16% of Caucasians express HLA-B*27:05 and/or HLA-B*57:01 and ~ 60% of Caucasian elite controllers—approximately 0.3% of infected individuals overall [18]—express one or the other of these alleles [18]. Therefore, only a minority of individuals expressing HLA-B*27:05 and/or HLA-B*57:01—approximately 2%—become elite controllers, for reasons which remain unclear. Factors proposed to explain differential outcome in subjects with matched HLA-B allele expression may be categorised as immunological, genetic and virological. Immunological factors include breadth of the Gag-specific CTL response [3], avidity of the T cell receptor for the peptide-MHC complex [32], and ability of the TCR repertoire to recognise potential escape variants [33]. Genetic factors include the HLA-A and HLA-C [34, 35] and KIR alleles expressed [36], and the presence of CCR5 variants such as the delta-32 mutant [37]. Virological factors include replicative capacity of the transmitted virus and pre-adaptation of the transmitted virus to the HLA class I molecules expressed in the transmission recipient [38,39,40,41].

A further influence on HIV disease outcome is the introduction of co-infection with another blood-borne virus. HCV may skew the immunological response to HIV, for example by increasing immune activation [42]. The immunological influence of HCV therapy on HIV is complex, as exemplified by the APRICOT study in which HCV/HIV co-infected patients receiving IFNα had a 1.0 log reduction in HIV viral load, but also a reduction in absolute CD4+ T cell counts [43]. In the patient we describe here, it is likely that intercurrent acute HCV infection (and/or IFN-based therapy) may have contributed to HIV progression. However, this co-infection cannot explain the whole picture as changes in CD4+ T cell count and HIV viral load that herald progression to AIDS were already evident before the HCV infection occurred, and also the extent of the sequence changes observed can only be explained by HIV superinfection.

In the case of the rapid progressor, the observation that several escape mutations were selected very early in the course of infection (prior to ‘time 0’) illustrates that the presence of a broad HIV-specific immune response does not provide complete protection against escape. Indeed, at the initial timepoint, when the CD4+ T cell count was 950 cells/mm3 and the HIV viral load was 65 copies/ml, five CTL responses were detectable (Fig. 3). Perhaps more significant is the observation of variants being generated in viral sequences derived from the non-progressor RI088 that did not get selected. Taken together with the data showing the relatively high viral replicative capacity of virus derived from RI088 compared to the R097 viruses, these findings are consistent with the quality of individual key HIV-specific responses being critical to outcome. These two case reports cannot provide definitive information on which specificities are most important to immune control, but are consistent with numerous previous studies indicating that the HLA-B*27:05-KK10 response plays a vital role in immune control of HIV in subjects expressing HLA-B*27:05.

The influence of genetic factors other than the HLA-B type has not been evaluated in these two study subjects. Although HLA-A and HLA-C alleles can influence immune control of HIV, HLA-B alleles have a greater impact than HLA-A or HLA-C [1, 35, 44]. The HLA-A alleles expressed by the progressor (A*11:01/24:02) and non-progressor (A*01:01/03:01) have not been associated with differential disease outcome and the HLA-C alleles expressed (C*02:02/06:02) were the same in both subjects.

The changes in viral replicative capacity over time in R097 are of interest, showing an initial significant decrease in VRC as the B27/57 escape mutants are selected, and then a partial recovery towards the VRC of the original infecting strain from R096. These findings are consistent with many studies that have shown the ability of the virus to limit the cost of escape mutants either via the selection of alternative escape mutants [28] and/or the co-selection of compensatory mutants [27, 30, 31, 45,46,47,48,49,50]. However, although the VRC returns towards normal, it remains well below that of the original transmitted virus. In contrast, the VRC of the virus in the elite controller remains very close to the reference NL4-3 despite more than a decade of strong selection pressure on the virus.

Conclusions

The description here of two HIV-infected subjects expressing both HLA-B*27:05 and HLA-B*57:01 with differential disease outcomes demonstrates the rapidity with which progression can occur, even when HIV-specific CTL responses restricted by these protective alleles are present. In this report rapid progression occurred following HIV superinfection. However, only a small minority of subjects expressing protective HLA alleles remain non-progressors. Consistent with previous studies, qualitative aspects of the HIV-specific CTL response are likely to play a critical part in maintaining long-term immune control of HIV infection.

Abbreviations

AIDS:

acquired immunodeficiency syndrome

ART:

antiretroviral therapy

CTL:

cytotoxic T lymphocyte

DNA:

deoxyribonucleic acid

ELISpot:

enzyme-linked immunosorbent spot assay

Gag:

group-specific antigen (capsid protein)

HCV:

hepatitis C virus

HIV:

human immunodeficiency virus

HLA:

human leukocyte antigen

IFN:

interferon

LOD:

limit of detection

MHC:

major histocompatibility complex

MSM:

men who have sex with men

Nef:

negative regulatory factor protein

PCR:

polymerase chain reaction

RNA:

ribonucleic acid

RT:

reverse transcription

VRC:

viral replicative capacity

References

  1. Pereyra F, Jia X, McLaren PJ, Telenti A, de Bakker PIW, Walker BD, et al. The major genetic determinants of HIV-1 control affect HLA class I peptide presentation. Science (80-). 2010;330:1551–7.

    Article  PubMed Central  Google Scholar 

  2. Goulder PJR, Walker BD. HIV and HLA class I: an evolving relationship. Immunity. 2012;37:426–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Kiepiela P, Ngumbela K, Thobakgale C, Ramduth D, Honeyborne I, Moodley E, et al. CD8+ T-cell responses to different HIV proteins have discordant associations with viral load. Nat Med. 2007;13:46–53.

    Article  CAS  PubMed  Google Scholar 

  4. Deng K, Pertea M, Rongvaux A, Wang L, Durand CM, Ghiaur G, et al. Broad CTL response is required to clear latent HIV-1 due to dominance of escape mutations. Nature. 2015;517:381–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Payne RP, Kløverpris H, Sacha JB, Brumme Z, Brumme C, Buus S, et al. Efficacious early antiviral activity of HIV Gag- and Pol-specific HLA-B 2705-restricted CD8+ T cells. J Virol. 2010;84:10543–57.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Brener J, Gall A, Batorsky R, Riddell L, Buus S, Leitman E, et al. Disease progression despite protective HLA expression in an HIV-infected transmission pair. Retrovirology. 2015;12:55.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Cano P, Klitz W, Mack SJ, Maiers M, Marsh SGE, Noreen H, et al. Common and well-documented HLA alleles: report of the ad-hoc committee of the American Society for Histocompatiblity and Immunogenetics. Hum Immunol. 2007;68:392–417.

    Article  CAS  PubMed  Google Scholar 

  8. Addo MM, Yu XG, Rathod A, Cohen D, Eldridge RL, Strick D, et al. Comprehensive epitope analysis of human immunodeficiency virus type 1 (HIV-1)-specific T-cell responses directed against the entire expressed HIV-1 genome demonstrate broadly directed responses, but no correlation to viral load. J Virol. 2003;77:2081–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Gall A, Ferns B, Morris C, Watson S, Cotten M, Robinson M, et al. Universal amplification, next-generation sequencing and assembly of HIV-1 genomes. J Clin Microbiol. 2012;50:3838–44.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Bentley DR, Balasubramanian S, Swerdlow HP, Smith GP, Milton J, Brown CG, et al. Accurate whole human genome sequencing using reversible terminator chemistry. Nature. 2008;456:53–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Gall A, Morris C, Kellam P. Complete genome sequence of the WHO international standard for HIV-1 RNA determined by deep sequencing. Genome Announc. 2014;2:10–1.

    Article  Google Scholar 

  12. Henn MR, Boutwell CL, Charlebois P, Lennon NJ, Power K, Macalalad AR, et al. Whole genome deep sequencing of HIV-1 reveals the impact of early minor variants upon immune recognition during acute infection. PLoS Pathog. 2012;8:e1002529.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Kurtz S, Phillippy A, Delcher AL, Smoot M, Shumway M, Antonescu C, et al. Versatile and open software for comparing large genomes. Genome Biol. 2004;5:R12.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Posada D. jModelTest: phylogenetic model averaging. Mol Biol Evol. 2008;25:1253–6.

    Article  CAS  PubMed  Google Scholar 

  15. Martin DP, Murrell B, Golden M, Khoosal A, Muhire B. RDP4: detection and analysis of recombination patterns in virus genomes. Virus Evol. 2015;1:1.

    Article  Google Scholar 

  16. Lole KS, Bollinger RC, Paranjape RS, Gadkari D, Kulkarni SS, Novak NG, et al. Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J Virol. 1999;73:152–60.

    CAS  PubMed  PubMed Central  Google Scholar 

  17. Brockman M, 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.

    Article  CAS  PubMed  Google Scholar 

  18. Miura T, Brockman MA, Brumme ZL, Brumme CJ, Pereyra F, Trocha A, et al. 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.

    Article  PubMed  Google Scholar 

  19. Gervaix A, West D, Leoni LM, Richman DD, Wong-Staal F, Corbeil J. A new reporter cell line to monitor HIV infection and drug susceptibility in vitro. Proc Natl Acad Sci. 1997;94:4653–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Tomiyama H, Sakaguchi T, Miwa K, Oka S, Iwamoto A, Kaneko Y, et al. Identification of multiple HIV-1 CTL epitopes presented by HLA-B*5101 molecules. Hum Immunol. 1999;60:177–86.

    Article  CAS  PubMed  Google Scholar 

  21. 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:641–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Carlson JM, Brumme CJ, Martin E, Listgarten J, Brockman M, 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:13202–16.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Matthews PC, Geretti AM, Goulder PJR, Klenerman P. Epidemiology and impact of HIV coinfection with hepatitis B and hepatitis C viruses in Sub-Saharan Africa. J Clin Virol. 2014;61:1–14.

    Article  Google Scholar 

  24. Martinez-picado J, Prado JG, Fry EE, Pfafferott K, Leslie A, Chetty S, et al. Fitness cost of escape mutations in p24 gag in association with control of human immunodeficiency virus type 1. J Virol. 2006;80:3617–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. 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:282–9.

    Article  CAS  PubMed  Google Scholar 

  26. 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:909–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Schneidewind A, Brockman MA, Yang R, Adam RI, Li B, Le Gall S, et al. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Goonetilleke N, Liu MKP, Salazar-Gonzalez JF, Ferrari G, Giorgi E, Ganusov VV, et al. The first T cell response to transmitted/founder virus contributes to the control of acute viremia in HIV-1 infection. J Exp Med. 2009;206:1253–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Schneidewind A, Brumme ZL, Brumme CJ, Power KA, Reyor LL, O’Sullivan K, et al. Transmission and long-term stability of compensated CD8 escape mutations. J Virol. 2009;83:3993–7.

    Article  CAS  PubMed  Google Scholar 

  30. Crawford H, Prado JG, Leslie A, Hué S, Honeyborne I, Reddy S, et al. Compensatory mutation partially restores fitness and delays reversion of escape mutation within the immunodominant HLA-B*5703-restricted Gag epitope in chronic human immunodeficiency virus type 1 infection. J Virol. 2007;81:8346–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Wright JK, Naidoo VL, Brumme ZL, Prince JL, Claiborne DT, Goulder PJR, et al. Impact of HLA-B*81-associated mutations in HIV-1 gag on viral replication capacity. J Virol. 2012;86:3193–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Almeida JR, Price D, Papagno L, Arkoub ZA, Sauce D, Bornstein E, et al. Superior control of HIV-1 replication by CD8+ T cells is reflected by their avidity, polyfunctionality, and clonal turnover. J Exp Med. 2007;204:2473–85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Chen H, Ndhlovu ZM, Liu D, Porter LC, Fang JW, Darko S, et al. TCR clonotypes modulate the protective effect of HLA class I molecules in HIV-1 infection. Nat Immunol. 2012;13:691–700.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Pereyra F, Addo MM, Kaufmann DE, Liu Y, Miura T, Rathod A, et al. Genetic and immunologic heterogeneity among persons who control HIV infection in the absence of therapy. J Infect Dis. 2008;197:563–71.

    Article  PubMed  Google Scholar 

  35. Apps R, Qi Y, Carlson JM, Chen H, Gao X, Thomas R, et al. Influence of HLA-C expression level on HIV control. Science. 2013;340:87–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Martin MP, Qi Y, Gao X, Yamada E, Martin JN, Pereyra F, et al. Innate partnership of HLA-B and KIR3DL1 subtypes against HIV-1. Nat Genet. 2007;39:733–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Dean M, Carrington M, Winkler C, Huttley GA, Smith MW, Allikmets R, et al. Genetic restriction of HIV-1 infection and progression to AIDS by a deletion allele of the CKR5 structural gene. Science (80-). 1996;273:1856–62.

    Article  CAS  Google Scholar 

  38. Carlson JM, Du VY, Pfeifer N, Bansal A, Tan VYF, Power K, et al. Impact of pre-adapted HIV transmission. Nat Med. 2016;22:606–13.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Claiborne DT, Prince JL, Scully E, Macharia G, Micci L, Lawson B, et al. Replicative fitness of transmitted HIV-1 drives acute immune activation, proviral load in memory CD4+ T cells, and disease progression. Proc Natl Acad Sci. 2015;112:E1480–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  40. Goepfert PA, Lumm W, Farmer P, Matthews P, Prendergast A, Carlson JM, et al. Transmission of HIV-1 Gag immune escape mutations is associated with reduced viral load in linked recipients. J Exp Med. 2008;205:1009–17.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Prince JL, Claiborne DT, Carlson JM, Schaefer M, Yu T, Lahki S, et al. Role of transmitted Gag CTL polymorphisms in defining replicative capacity and early HIV-1 pathogenesis. PLoS Pathog. 2012;8:e1003041.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Gonzalez VD, Falconer K, Blom KG, Reichard O, Morn B, Laursen AL, et al. High levels of chronic immune activation in the T-cell compartments of patients coinfected with hepatitis C virus and human immunodeficiency virus type 1 and on highly active antiretroviral therapy are reverted by alpha interferon and ribavirin treatment. J Virol. 2009;83:11407–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Torriani FJ, Rodriguez-Torres M, Rockstroh JK, Lissen E, Gonzalez-Garcia J, Lazzarin A, et al. Peginterferon Alfa-2a plus ribavirin for chronic hepatitis C virus infection in HIV-infected patients. N Engl J Med. 2004;351:438–50.

    Article  CAS  PubMed  Google Scholar 

  44. Kiepiela P, Leslie AJ, Honeyborne I, Ramduth D, Thobakgale C, Chetty S, et al. Dominant influence of HLA-B in mediating the potential co-evolution of HIV and HLA. Nature. 2004;432:769–75.

    Article  CAS  PubMed  Google Scholar 

  45. Schneidewind A, Brockman MA, Sidney J, Wang YE, Chen H, Suscovich TJ, et al. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Rolland M, Carlson JM, Manocheewa S, Swain JV, Lanxon-Cookson E, Deng W, et al. Amino-acid co-variation in HIV-1 Gag subtype C: HLA-mediated selection pressure and compensatory dynamics. PLoS ONE. 2010;5:e12463.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Brockman MA, Brumme ZL, Brumme CJ, Miura T, Sela J, Rosato PC, et al. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Crawford H, Matthews PC, Schaefer M, Carlson JM, Leslie A, Kilembe W, et al. The hypervariable HIV-1 capsid protein residues comprise HLA-driven CD8+ T-cell escape mutations and covarying HLA-independent polymorphisms. J Virol. 2011;85:1384–90.

    Article  CAS  PubMed  Google Scholar 

  49. Huang KHG, Goedhals D, Carlson JM, Brockman MA, Mishra S, Brumme ZL, et al. Progression to AIDS in South Africa is associated with both reverting and compensatory viral mutations. PLoS ONE. 2011;6:e19018.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Tsai M-H, Muenchhoff M, Adland E, Carlqvist A, Roider J, Cole DK, et al. Paediatric non-progression following grandmother-to-child HIV transmission. Retrovirology. 2016;13:65.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Authors’ contributions

JB carried out experimental work, data analysis and contributed to manuscript preparation. Ultra-deep sequencing was performed in collaboration with AG and PK. AG performed the de novo assembly of consensus sequences. RB performed minor variant sequence haplotype analysis. JH designed the logistic curve models for analysing VRC data. NL contributed to clonal sequencing. FC and AE recruited the study participants. CB and RD contributed to sequence data analysis. OP interpreted phylogenetic data. PM supervised the study and contributed to manuscript preparation. PG conceived and supervised the study and prepared manuscript text. All authors read and approved the final manuscript.

Acknowledgements

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

Sequence data generated during this study is available from Genbank (accession numbers MF039091-MF039203). All other data generated is included in this publication or available from the corresponding author on reasonable request.

Consent for publication

Consent to publish has been obtained as part of the written consent mentioned above. Study participants have been anonymous and cannot be identified from information in this manuscript.

Ethics approval and consent to participants

All subjects gave written informed consent for their participation. Ethics approval was given by the Oxfordshire Rec A Committee (Study Number 06/Q1604/12).

Funding

This work was funded by grants from the National Institutes of Health (RO1AI46995 to PJRG), the Wellcome Trust (WT104748MA to PJRG), PCM received salary support from NIHR and the Wellcome Trust (Grant Number 110110/Z/15/Z), JB received support from the Commonwealth Scholarship Commission and Oppenheimer Memorial Trust. OP was funded by the ERC under European Research Council Grant Agreement Number 614725-PATHPHYLODYN.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Jacqui Brener.

Additional files

Additional file 1: Table 1.

The parameters used for the curves shown in Figure 4, together with the statistics for the overall fit and for the separate parameters.

Additional file 2: Fig. 1.

Phylogenetic analysis of clonal HIV sequences from RI088. Maximum likelihood phylogenetic tree of 1091bp alignment of clonal DNA sequences across the Gag p17 and p24 genes. RI088 sequences are shown in colour (0.6 years in purple, 1.2 years in pink, 1.4 years in green, 3.4 years in orange, 10.8 years in blue, 11.3 years in red). 143 B B clade reference sequences from the US and UK collected between 2003 and 2011 from the Los Alamos database (https://www.hiv.lanl.gov/) are shown in black. Bootstrap values based on 1000 bootstrap replicates are shown in italics.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Brener, J., Gall, A., Hurst, J. et al. Rapid HIV disease progression following superinfection in an HLA-B*27:05/B*57:01-positive transmission recipient. Retrovirology 15, 7 (2018). https://0-doi-org.brum.beds.ac.uk/10.1186/s12977-018-0390-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12977-018-0390-9

Keywords