Impact of HLA Class I Alleles on Timing of HIV Rebound After Antiretroviral Treatment Interruption

Published on November 27, 2017
Read in PDF format

Correction

This is a corrected version of this article, published on June 27, 2018. The original version of this article, published on November 27, 2017, was missing an informed consent statement.

Authors

You Jeong Park1,2, Behzad Etemad1, Hayat Ahmed1, Vivek Naranbhai3,4,5, Evgenia Aga2, Ronald J. Bosch6, John W. Mellors7, Daniel R. Kuritzkes1, Michael Para8, Rajesh T. Gandhi3, Mary Carrington3,4, Jonathan Z. Li1

Affiliated Institutions

1 Brigham and Women’s Hospital, Harvard Medical School, Boston, Massachusetts
2 Harvard College, Cambridge, Massachusetts
3 Massachusetts General Hospital and Ragon Institute, Harvard Medical School, Boston, Massachusetts
4 Cancer and Inflammation Program, Laboratory of Experimental Immunology, Leidos Biomedical Research Institute, Frederick National Laboratory for Cancer Research, Frederick, Maryland
5 Centre for the AIDS Programme of Research in South Africa, Durban, KwaZulu Natal, South Africa
6 Harvard T.H. Chan School of Public Health, Boston, Massachusetts
7 University of Pittsburgh, Pittsburgh, Pennsylvania
8 Ohio State University, Columbus, Ohio

Corresponding Author

Jonathan Li
Division of Infectious Diseases
Brigham and Women’s Hospital
65 Landsdowne St, Rm 421
Cambridge, MA 02139
Tel: (617) 768-8476
Fax: (617) 768-8738
[email protected]

Running Title

Association of HLA alleles with HIV reservoir and viral rebound timing

Suggested Citation

Park Y, Etemad B, Ahmed H, Naranbhai V, Aga E, Bosch RJ, Mellors JW, Kuritzkes DR, Para M, Gandhi RT, Carrington M, Li JZ. Impact of HLA Class I Alleles on Timing of HIV Rebound After Antiretroviral Treatment Interruption. Pathogens and Immunity. 2017;2(3):431-45. doi: 10.20411/pai.v2i3.222

Abstract

Background: Identifying host determinants associated with HIV reservoir size and timing of viral rebound after an analytic treatment interruption (ATI) is an important step in the search for an HIV functional cure. We performed a pooled analysis of 103 participants from 4 AIDS Clinical Trials Group ATI studies to identify the association between HLA class I alleles with HIV reservoir size and viral rebound timing.

Methods: Total HIV DNA and cell-associated HIV RNA (CA-RNA) were quantified in pre-ATI peripheral blood mononuclear cell samples, and residual plasma viremia was measured using the single-copy assay. HLA class I typing was performed, and we generated an odds ratio (OR) of predicted HLA effect on HIV viremia control for each individual and compared this with time to viral rebound, and levels of HIV DNA and CA-RNA.

Results: There was no significant association between the HLA ORs and levels of HIV DNA or CA-RNA, but carriage of protective HLA-B alleles (lower OR scores) was associated with delayed viral rebound (P=0.02). Higher OR scores at the HLA-C locus were associated with longer duration of ART treatment (P = 0.02) and this trend was also seen with the combined OR score (P < 0.01). Individuals with protective HLA-B alleles had delayed viral rebound after treatment interruption that was not explained by differences in baseline reservoir size.

Conclusions: The results indicate the vital role of cellular host immunity in preventing HIV rebound and the importance of taking into account the HLA status of study participants being evaluated in trials for an HIV cure.

Keywords

HLA, HIV reservoir, viral rebound timing, analytic treatment interruption

Introduction

Although antiretroviral therapy (ART) is effective at suppressing HIV replication in infected individuals, the virus persists in a stable pool of resting CD4+ T cells in the form of latent provirus that almost inevitably rebounds when treatment is stopped [1]. Eradication of HIV-1 disease necessitates the elimination of this latent reservoir, and identifying host determinants associated with latency and reservoir size in patients receiving ART is an important step in the search for an intervention that provides sustained, long-term ART-free virological remission [2].

One of the main host genetic determinants that influences the control of HIV and disease progression is the human leukocyte antigen (HLA) class I locus [3]. Specific HLA class I alleles like HLA-B*57 and B*27 have consistently been shown to display protective effects and are associated with viral load control and delayed progression to AIDS [4, 5], while the allelic group HLA-B*35Px is associated with accelerated disease progression [6]. A combined analysis of 9 HIV cohorts has since provided additional details for HLA class I allele-specific effects in maintaining HIV control [7]. Using the results from that study, we examined the association between alleles at 3 HLA class I loci (HLA-A, B, and C) and both reservoir size and viral dynamics after treatment interruption in participants of the AIDS Clinical Trials Group (ACTG) who underwent an analytic treatment interruption (ATI). Identifying host genetic factors influencing the size of the infected cell population will provide insights into the mechanisms behind viral control after treatment interruption with implications for the efforts to induce ART-free HIV remission.

Methods

Study Population

Participants from 4 ACTGATI studies (A5024 [8], A5068 [9], A5170 [10], and A5197 [11]) were included if they were receiving suppressive ART, received no immunological interventions (eg, therapeutic vaccination), had HIV-1 RNA less than 50 copies/mL at the time of ATI, and had samples available for HLA typing. These participants received ART during chronic infection and have been previously described in an analysis showing the virological predictors of HIV rebound timing [12].

Written informed consent was provided by all study participants for use of stored samples in HIV-related research. This study was approved by the Partners Institutional Review Board.

Human leukocyte antigen typing

HLA class I typing was performed following the PCR-SSOP (sequence-specific oligonucleotide probing) and the PCR-SBT (sequence based typing) protocols recommended by the 13th International Histocompatibility Workshop (http://www.ihwg.org). Each HLA allele was assigned an odds ratio (OR) for its effects on HIV disease progression based on a previously published multivariable logistic regression analysis of HIV-1 controllers and non-controllers (Supplementary Table 4) [7]. Controllers refer to HIV-1 infected individuals whose mean viral loads are maintained at less than 2,000 copies/mL of plasma despite not receiving ART. Non-controllers refer to individuals whose mean viral loads surpass 10,000 copies/mL. Based on this model, alleles with OR values less than 1 are associated with the protective phenotype seen in controllers while OR values greater than 1 are associated with the unfavorable phenotype of rapid disease progression. For our analysis, alleles that did not reach significance in the model were given an OR value of 1 and defined as neutral. For each participant, the OR’s of the 2 alleles in each HLA locus were multiplied to calculate ORA, ORB, and ORC for the HLA-A, B, and C locus, respectively. The OR scores for all loci were also multiplied to produce a cumulative OR score (ORABC) to reflect the combined effects of all HLA alleles. All HLA B*35Px group alleles (B*3502, 3503, 3504, and B*5301) [6] were assigned the same OR values as B*3502.

As a sensitivity analysis, participants were also grouped into 3 discrete categories: protective, neutral, or unfavorable groups. Participants with at least 1 protective and no unfavorable alleles were categorized as having “protective” HLA genotypes, those with at least 1 unfavorable allele and no protective alleles were categorized as “unfavorable”, and those with 1 protective and 1 unfavorable allele or those with neither protective nor unfavorable alleles were categorized as “neutral”. Categorization of alleles was based on the previously described OR values [7]. HLA-B*57/27 alleles were considered protective regardless of the last 2 digits, and HLA-B*35Px were considered unfavorable for all races.

HIV-1 reservoir quantification

Cell-associated HIV-1 RNA (CA-RNA) and DNA were isolated from cryopreserved PBMCs using the AllPrep DNA/RNA Mini Kit (Qiagen). Cellular integrity for RNA analysis was assessed by the measurement of total extracted RNA and evaluation of the IPO-8 housekeeping gene [13]. Unspliced CA-RNA and total HIV DNA were quantified using a real-time PCR approach with primers/probes targeting conserved regions of HIV LTR/gag as previously described [14].

CD4+ cell count and pre- and post-ATI viral load measurements

The CD4+ cell count before ATI and plasma HIV-1 RNA were defined as the most recent measurements on or before the date of ART discontinuation. The timing of viral rebound was defined as a confirmed viral load ≥ 200 copies/mL or (2) a single viral load ≥ 1000 copies/mL. Viral load set point was defined as the mean viral load between weeks 12 and 16. The timing of viral rebound was categorized as early (≤ 4 weeks), intermediate (5 to 8 weeks), or delayed (> 8 weeks) after ATI.

Statistical analysis

Associations between OR scores (ORA, ORB, ORC, and ORABC) and reservoir size measurements or viral rebound times were analyzed using Spearman correlation, Kruskal-Wallis, Wilcoxon rank sum, and Fisher’s exact tests. Associations between categorical HLA-groups and reservoir size measurements or viral rebound times were also analyzed using Fisher’s exact tests.

Results

Study participants and baseline characteristics

A total of 103 participants were included from the pooled ACTG studies with ATI. Table 1 lists their baseline characteristics. Participants in this study had a median age of 42 years, 90% were male, and 83% were white. The median (Q1, Q3) CD4+ cell count at the start of ATI was 843 (687, 1042) cells/µL. Participants had received ART for a median of 5.3 years.


Table 1. Baseline characteristics of participants included in the study.

Characteristic

Total (N=103)

Sex, male, n (%)

93 (90%)

Age, median years (Q1, Q3)

42 (38, 50)

Race/Ethnicity, n (%)

White

85 (83%)

Black

18 (17%)

Nadir CD4+ cell count, median cells/µL (Q1, Q3)

422 (351, 544)

Pre-ATI CD4+ cell count, median cells/ µL (Q1, Q3)

850 (686, 1048)

Duration of ART, median years (Q1, Q3)

5.3 (3.2, 6.6)

NNRTI-based ART

65 (63%)

Source study, n (%)

A5024

7 (7%)

A5068

12 (12%)

A5170

62 (60%)

A5197

22 (21%)

NNRTI, non-nucleoside reverse transcriptase inhibitor; ART, antiretroviral therapy


There were no significant associations between background characteristics and the estimated effect of each HLA locus individually (ORA, ORB or ORc) and the overall effect (HLAABC) scores, with the exception of ART duration (Supplementary Tables 1-3). A positive correlation was detected between ORABC and the duration of ART (Spearman r = 0.26, P< 0.01, Figure 1), indicating that participants harboring HLA alleles associated with greater risk of HIV non-control had been treated with ART for a longer period of time. A similar trend could be seen for each of the individual HLA loci, but only the ORc analysis reached statistical significance (Spearman r = 0.24, P = 0.01). These results were consistent with the categorical analysis where protective HLA-C alleles were associated with a shorter duration of ART (protective vs neutral: median 2.5 vs 5.3 years, exact Wilcoxon P < 0.01; protective vs unfavorable: median 2.5 vs 5.8 years, Wilcoxon P < 0.01).


Figure 1. Correlation between ORABC and duration of ART. ORABC represents an overall assessment of protective status and is calculated by multiplying the ORs of HIV control of all 6 alleles at the A, B, and C loci. A lower ORABC is associated with greater control of HIV viremia.


Protective effect of HLA is associated with delayed viral rebound timing

Participants with delayed viral rebound were found to have significantly more protective HLA-B alleles (lower ORB scores) compared to those with early or intermediate rebound times (early vs delayed: median OR B 1 vs 0.6, P = 0.02; intermediate vs delayed: median 1.35 vs 0.6, Figure 2). These results were consistent with the categorical analysis, which showed that there were significant differences in the distribution of viral rebound timing by HLA-B category (P = 0.02, Supplementary Figure 1). The results were largely unchanged with the single ≥ 1,000 HIV-1 RNA copies/mL definition of viral rebound. No significant associations were found in the protective effect of the HLA alleles and either pre-ATI levels of HIV DNA, CA-RNA, or the post-ATI viral load set point.


Figure 2. Association between ORB and timing of viral rebound. ORB represents the protective status at the HLA-B locus and is calculated for each individual by multiplying the OR estimate of HIV non-control for each of the 2 alleles carried by that individual.


Discussion

In this pooled analysis of ACTG studies with ATI, we evaluated the association between HLA class I alleles with HIV reservoir size and post-ATI viral dynamics using a new method of assessing the additive effect of alleles at multiple loci. We found that the presence of protective HLA-B alleles was associated with delayed viral rebound after treatment interruption that was not readily explained by differences in baseline reservoir size. One key immune response associated with control of viral replication is the HIV-1-specific CD8+ T lymphocyte responses which are modulated by the host’s HLA class I molecules [16]. The role that virus-specific T lymphocytes may play in the control of viral replication is suggested by the association between the appearance of HIV-1-specific CD8+ T lymphocytes in peripheral blood with initial decline in viremia during primary infection [17]. This link is further supported by the loss of the ability to control viral replication in macaques infected with simian immunodeficiency virus (SIV) once depleted of CD8+ lymphocytes [18]. Although the complete mechanism underlying HLA modulation of HIV disease progression is not well understood, class I alleles associated with the protective phenotype of slower disease progression strongly contribute to the total HIV-specific CD8+ T lymphocyte response during acute infection [19]. In this study, higher OR scores (ie, more unfavorable HLA alleles) were associated with a longer time receiving ART. One possible explanation is that unfavorable HLA alleles may lead to more rapid disease progression and therefore to patients starting ART earlier.

To our knowledge, this is the first study demonstrating that protective HLA alleles, specifically in the HLA-B loci, are associated with delayed time to viral rebound, indicating the vital role that cellular host immunity has in preventing HIV rebound after treatment interruption. Given the lack of evidence for an association between HLA alleles and pre-ATI reservoir size (HIV DNA) or activity (CA-RNA), this effect may be mediated by the identification and control of reactivated T cells and viral spread after ART discontinuation. The results also suggest that boosting T-cell-mediated immune responses could be an important component of efforts to induce sustained ART-free HIV remission. In the SPARTAC study, participants with unfavorable HLA alleles had higher HIV DNA levels prior to ART initiation, and only pre-ART levels of T-cell exhaustion predicted the time of viral rebound [20]. However, the SPARTAC study only evaluated participants who initiated ART during early HIV infection, and their sample size was less than half that of this study, limiting the authors' ability to assess HLA-mediated effects. Despite the fact that protective class I HLA alleles are associated with the spontaneous control of HIV, these alleles were not found to be enriched in the VISCONTI cohort of post-treatment controllers [21]. One potential explanation is that while the VISCONTI study limited the categorization of protective HLA alleles to HLA-B*27 and B*57, our analysis provided a more comprehensive and nuanced approach to assessing HLA-related effects. It is also possible that non-HLA-mediated immune responses are needed for sustained HIV remission. In summary, individuals with protective HLA-B alleles had delayed viral rebound after treatment interruption that was not explained by differences in baseline reservoir size. The results indicate the vital role of cellular host immunity in preventing HIV rebound and the importance of taking into account the HLA status of study participants being evaluated in cure trials for HIV because an imbalance in HLA types between study arms may have an unexpected effect on the outcome.

Acknowledgments

We thank the participants, staff, and principal investigators of the ACTG studies A5024 (J. Michael Kilby, Ronald Mitsuyasu), A5068 (Jeffrey Jacobson, Ian Frank, Michael Saag, Joseph Eron), A5170 (Daniel Skiest, David Margolis, Diane Havlir), A5187 (Daniel Barouch, Eric Rosenberg, Daniel Kuritzkes), and A5197 (Robert Schooley, Michael Lederman, Diane Havlir).

Financial Support

This work was supported in part by a grant from the Harvard University Center for AIDS Research (to Drs. Li and Gandhi, NIAID 5P30AI060354-08), National Institutes of Health (NIH) grants AI125109 (to Dr. Li), UM1 AI068634 (Statistical and Data Management Center of the AIDS Clinical Trials Group), UM1 AI068636 (AIDS Clinical Trials Group), and a subcontract from UM1 AI068636 to the Harvard Virology Support Laboratory (to Dr. Kuritzkes).

Potential Conflicts of Interest

Dr. Li has received research support and consulted for Gilead Sciences and Merck. Dr. Gandhi has received grant support from Gilead, Merck, and Viiv.

Footnotes

Submitted August 25, 2017 | Accepted September 12, 2017 | Published November 27, 2017

References

1. Wong JK, Hezareh M, Gunthard HF, Havlir DV, Ignacio CC, Spina CA, Richman DD. Recovery of replication-competent HIV despite prolonged suppression of plasma viremia. Science. 1997;278(5341):1291-5. PubMed PMID: 9360926.

2. Abdel-Mohsen M, Wang C, Strain MC, Lada SM, Deng X, Cockerham LR, Pilcher CD, Hecht FM, Liegler T, Richman DD, Deeks SG, Pillai SK. Select host restriction factors are associated with HIV persistence during antiretroviral therapy. AIDS. 2015;29(4):411-20. PubMed PMID: 25602681. Pubmed Central PMCID: PMC4385712. doi: 10.1097/QAD.0000000000000572

3. Bashirova AA, Thomas R, Carrington M. HLA/KIR restraint of HIV: surviving the fittest. Annu Rev Immunol. 2011;29:295-317. PubMed PMID: 21219175. Pubmed Central PMCID: PMC3725604. doi: 10.1146/annurev-immunol-031210-101332

4. Migueles SA, Sabbaghian MS, Shupert WL, Bettinotti MP, Marincola FM, Martino L, Hallahan CW, Selig SM, Schwartz D, Sullivan J, Connors M. HLA B*5701 is highly associated with restriction of virus replication in a subgroup of HIV-infected long term nonprogressors. Proc Natl Acad Sci U S A. 2000;97(6):2709-14. PubMed PMID: 10694578. Pubmed Central PMCID: PMC15994. doi: 10.1073/pnas.050567397

5. Kiepiela P, Leslie AJ, Honeyborne I, Ramduth D, Thobakgale C, Chetty S, Rathnavalu P, Moore C, Pfafferott KJ, Hilton L, Zimbwa P, Moore S, Allen T, Brander C, Addo MM, Altfeld M, James I, Mallal S, Bunce M, Barber LD, Szinger J, Day C, Klenerman P, Mullins J, Korber B, Coovadia HM, Walker BD, Goulder PJ. Dominant influence of HLA-B in mediating the potential co-evolution of HIV and HLA. Nature. 2004;432(7018):769-75. PubMed PMID: 15592417. doi: 10.1038/nature03113

6. Gao X, Nelson GW, Karacki P, Martin MP, Phair J, Kaslow R, Goedert JJ, Buchbinder S, Hoots K, Vlahov D, O'Brien SJ, Carrington M. Effect of a single amino acid change in MHC class I molecules on the rate of progression to AIDS. N Engl J Med. 2001;344(22):1668-75. PubMed PMID: 11386265. doi: 10.1056/NEJM200105313442203

7. Bashirova AA, Martin-Gayo E, Jones DC, Qi Y, Apps R, Gao X, Burke PS, Taylor CJ, Rogich J, Wolinsky S, Bream JH, Duggal P, Hussain S, Martinson J, Weintrob A, Kirk GD, Fellay J, Buchbinder SP, Goedert JJ, Deeks SG, Pereyra F, Trowsdale J, Lichterfeld M, Telenti A, Walker BD, Allen RL, Carrington M, Yu XG. LILRB2 interaction with HLA class I correlates with control of HIV-1 infection. PLoS Genet. 2014;10(3):e1004196. PubMed PMID: 24603468. Pubmed Central PMCID: PMC3945438. doi: 10.1371/journal.pgen.1004196

8. Kilby JM, Bucy RP, Mildvan D, Fischl M, Santana-Bagur J, Lennox J, Pilcher C, Zolopa A, Lawrence J, Pollard RB, Habib RE, Sahner D, Fox L, Aga E, Bosch RJ, Mitsuyasu R, Adult ACTGAPT. A randomized, partially blinded phase 2 trial of antiretroviral therapy, HIV-specific immunizations, and interleukin-2 cycles to promote efficient control of viral replication (ACTG A5024). J Infect Dis. 2006;194(12):1672-6. PubMed PMID: 17109338. doi: 10.1086/509508

9. Jacobson JM, Pat Bucy R, Spritzler J, Saag MS, Eron JJ, Jr., Coombs RW, Wang R, Fox L, Johnson VA, Cu-Uvin S, Cohn SE, Mildvan D, O'Neill D, Janik J, Purdue L, O'Connor DK, Vita CD, Frank I, National Institute of A, Infectious Diseases ACTGPT. Evidence that intermittent structured treatment interruption, but not immunization with ALVAC-HIV vCP1452, promotes host control of HIV replication: the results of AIDS Clinical Trials Group 5068. J Infect Dis. 2006;194(5):623-32. PubMed PMID: 16897661. doi: 10.1086/506364

10. Skiest DJ, Su Z, Havlir DV, Robertson KR, Coombs RW, Cain P, Peterson T, Krambrink A, Jahed N, McMahon D, Margolis DM, Team ACTGS. Interruption of antiretroviral treatment in HIV-infected patients with preserved immune function is associated with a low rate of clinical progression: a prospective study by AIDS Clinical Trials Group 5170. J Infect Dis. 2007;195(10):1426-36. PubMed PMID: 17436222. doi: 10.1086/512681

11. Schooley RT, Spritzler J, Wang H, Lederman MM, Havlir D, Kuritzkes DR, Pollard R, Battaglia C, Robertson M, Mehrotra D, Casimiro D, Cox K, Schock B, Team ACTGS. AIDS clinical trials group 5197: a placebo-controlled trial of immunization of HIV-1-infected persons with a replication-deficient adenovirus type 5 vaccine expressing the HIV-1 core protein. J Infect Dis. 2010;202(5):705-16. PubMed PMID: 20662716. Pubmed Central PMCID: PMC2916952. doi: 10.1086/655468

12. Li JZ, Etemad B, Ahmed H, Aga E, Bosch RJ, Mellors JW, Kuritzkes DR, Lederman MM, Para M, Gandhi RT. The size of the expressed HIV reservoir predicts timing of viral rebound after treatment interruption. AIDS. 2016;30(3):343-53. PubMed PMID: 26588174. doi: 10.1097/QAD.0000000000000953

13. Ledderose C, Heyn J, Limbeck E, Kreth S. Selection of reliable reference genes for quantitative real-time PCR in human T cells and neutrophils. BMC Res Notes. 2011;4:427. PubMed PMID: 22011438. Pubmed Central PMCID: PMC3229292. doi: 10.1186/1756-0500-4-427

14. Malnati MS, Scarlatti G, Gatto F, Salvatori F, Cassina G, Rutigliano T, Volpi R, Lusso P. A universal real-time PCR assay for the quantification of group-M HIV-1 proviral load. Nat Protoc. 2008;3(7):1240-8. PubMed PMID: 18600229. doi: 10.1038/nprot.2008.108

15. Palmer S, Wiegand AP, Maldarelli F, Bazmi H, Mican JM, Polis M, Dewar RL, Planta A, Liu S, Metcalf JA, Mellors JW, Coffin JM. New real-time reverse transcriptase-initiated PCR assay with single-copy sensitivity for human immunodeficiency virus type 1 RNA in plasma. J Clin Microbiol. 2003;41(10):4531-6. PubMed PMID: 14532178. Pubmed Central PMCID: PMC254331.

16. Kosmrlj A, Read EL, Qi Y, Allen TM, Altfeld M, Deeks SG, Pereyra F, Carrington M, Walker BD, Chakraborty AK. Effects of thymic selection of the T-cell repertoire on HLA class I-associated control of HIV infection. Nature. 2010;465(7296):350-4. PubMed PMID: 20445539. Pubmed Central PMCID: PMC3098720. doi: 10.1038/nature08997

17. Borrow P, Lewicki H, Hahn BH, Shaw GM, Oldstone MB. Virus-specific CD8+ cytotoxic T-lymphocyte activity associated with control of viremia in primary human immunodeficiency virus type 1 infection. J Virol. 1994;68(9):6103-10. PubMed PMID: 8057491. Pubmed Central PMCID: PMC237022.

18. Schmitz JE, Kuroda MJ, Santra S, Sasseville VG, Simon MA, Lifton MA, Racz P, Tenner-Racz K, Dalesandro M, Scallon BJ, Ghrayeb J, Forman MA, Montefiori DC, Rieber EP, Letvin NL, Reimann KA. Control of viremia in simian immunodeficiency virus infection by CD8+ lymphocytes. Science. 1999;283(5403):857-60. PubMed PMID: 9933172.

19. Altfeld M, Kalife ET, Qi Y, Streeck H, Lichterfeld M, Johnston MN, Burgett N, Swartz ME, Yang A, Alter G, Yu XG, Meier A, Rockstroh JK, Allen TM, Jessen H, Rosenberg ES, Carrington M, Walker BD. HLA Alleles Associated with Delayed Progression to AIDS Contribute Strongly to the Initial CD8(+) T Cell Response against HIV-1. PLoS Med. 2006;3(10):e403. PubMed PMID: 17076553. Pubmed Central PMCID: PMC1626551. doi: 10.1371/journal.pmed.0030403

20. Hurst J, Hoffmann M, Pace M, Williams JP, Thornhill J, Hamlyn E, Meyerowitz J, Willberg C, Koelsch KK, Robinson N, Brown H, Fisher M, Kinloch S, Cooper DA, Schechter M, Tambussi G, Fidler S, Babiker A, Weber J, Kelleher AD, Phillips RE, Frater J. Immunological biomarkers predict HIV-1 viral rebound after treatment interruption. Nat Commun. 2015;6:8495. PubMed PMID: 26449164. Pubmed Central PMCID: 4633715. doi: 10.1038/ncomms9495

21. Saez-Cirion A, Bacchus C, Hocqueloux L, Avettand-Fenoel V, Girault I, Lecuroux C, Potard V, Versmisse P, Melard A, Prazuck T, Descours B, Guergnon J, Viard JP, Boufassa F, Lambotte O, Goujard C, Meyer L, Costagliola D, Venet A, Pancino G, Autran B, Rouzioux C. Post-Treatment HIV-1 Controllers with a Long-Term Virological Remission after the Interruption of Early Initiated Antiretroviral Therapy ANRS VISCONTI Study. PLoS Pathog. 2013;9(3):e1003211. PubMed PMID: 23516360. Pubmed Central PMCID: 3597518. doi: 10.1371/journal.ppat.1003211

Supplementary Materials

Supplementary Table 1. Baseline characteristics for HLA-A groups

Characteristic

HLA-A

Protective (n=7)

Neutral (n=39)

Unfavorable (n=57)

Sex, male, n (%)

7 (100%)

34 (87%)

51 (89%)

Age, median years (Q1, Q3)

38 (33, 60)

42 (36, 48)

43 (39, 50)

Race/Ethnicity, n(%)

White

7 (100%)

24 (62%)

54 (95%)

Black

0 (0%)

15 (38%)

3 (5%)

Nadir CD4+ cell count, median cells/µL (Q1, Q3)

476 (402, 587)

396 (305, 505)

422 (350, 566)

Pre-ATI CD4+ cell count, median cells/ µL (Q1, Q3)

905 (798, 1161)

863 (617, 1099)

832 (681, 1010)

Pre-ART viral load, median log10 HIV RNA copies/mL (Q1, Q3)

4.0 (3.3, 4.7)

4.4 (3.8, 4.9)

4.4 (4.0, 4.9)

CA-DNA, median log10 copies/106 PBMCs (Q1, Q3)

2.3 (2.2, 2.9)

2.4 (1.9, 3.1)

2.3 (1.7, 2.9)

CA-RNA, median log10 copies/106 PBMCs (Q1, Q3)

2.0 (1.1, 2.6)

1.7 (1.1, 2.3)

1.8 (1.1, 2.3)

Duration of ART, median years (Q1, Q3)

3.8 (2.6, 5.3)

5.4 (3.2, 6.5)

5.4 (3.3, 7.2)

NNRTI-based ART

4 (57%)

19 (49%)

41 (72%)

Source study, n (%)

A5024

0 (0%)

3 (8%)

4 (7%)

A5068

0 (0%)

8 (21%)

4 (7s%)

A5170

5 (71%)

22 (56%)

35 (61%)

A5197

2 (29%)

6 (15%)

14 (25%)


Supplementary Table 2. Baseline characteristics for HLA-B groups

Characteristic

HLA-B

Protective (n=22)

Neutral (n=48)

Unfavorable (n=33)

Sex, male, n (%)

19 (86%)

41 (85%)

32 (97%)

Age, median years (Q1, Q3)

42 (38, 46)

43 (38, 49)

41 (48, 53)

Race/Ethnicity, n (%)

White

18 (82%)

43 (90%)

24 (73%)

Black

4 (18%)

5 (10%)

9 (27%)

Nadir CD4+ cell count, median cells/µL (Q1, Q3)

443 (383, 536)

413 (353, 544)

449 (328, 568)

Pre-ATI CD4+ cell count, median cells/µL (Q1, Q3)

876 (626, 1030)

859 (688, 1061)

800 (687, 1127)

Pre-ART viral load, median log10 HIV RNA copies/mL (Q1, Q3)

4.2 (3.9, 4.6)

4.3 (3.7, 4.9)

4.5 (3.8, 5.0)

CA-DNA, median log10 copies/106 PBMCs (Q1, Q3)

2.1 (1.7, 2.8)

2.4 (1.9, 2.9)

2.4 (1.7, 3.1)

CA-RNA, median log10 copies/106 PBMCs (Q1, Q3)

1.6 (1.1, 2.1)

1.8 (1.1, 2.3)

2.0 (1.1, 2.5)

Duration of ART, median years (Q1, Q3)

4.2 (2.8, 5.8)

5.4 (3.3, 6.7)

5.4 (3.3, 8.3)

NNRTI-based ART

12 (55%)

34 (71%)

18 (55%)

Source study, n (%)

A5024

3 (14%)

3 (6%)

1 (3%)

A5068

0 (0%)

7 (15%)

5 (15%)

A5170

14 (64%)

29 (60%)

19 (58%)

A5197

5 (23%)

9 (19%)

8 (24%)


Supplementary Table 3. Baseline characteristics for HLA-C groups

Characteristic

HLA-C

Protective (n=8)

Neutral (n=79)

Unfavorable (n=16)

Sex, male, n (%)

7 (88%)

71 (90%)

14 (88%)

Age, median years (Q1, Q3)

38 (35, 51)

43 (38, 50)

43 (40, 52)

Race/Ethnicity, n (%)

White

5 (63%)

64 (81%)

16 (100%)

Black

3 (37%)

15 (19%)

0 (0%)

Nadir CD4+ cell count, median cells/µL (Q1, Q3)

451 (402, 510)

413 (325, 582)

422 (363, 524)

Pre-ATI CD4+ cell count, median cells/ µL (Q1, Q3)

799 (601, 1310)

852 (653, 1023)

852 (764, 1186)

Pre-ART viral load, median log10 HIV RNA copies/mL (Q1, Q3)

4.5 (3.8, 5.3)

4.3 (3.8, 4.7)

4.9 (4.0, 5.5)

CA-DNA, median log10 copies/106 PBMCs (Q1, Q3)

2.7 (1.7, 3.3)

2.3 (1.8, 3.1)

2.5 (1.7, 2.7)

CA-RNA, median log10 copies/106 PBMCs (Q1, Q3)

1.8 (1.2, 2.6)

1.8 (1.1, 2.3)

1.8 (1.1, 2.3)

Duration of ART, median years (Q1, Q3)

2.5 (1.6, 4.3)

5.3 (3.4, 6.7)

5.8 (3.5, 7.6)

NNRTI-based ART

5 (63%)

50 (63%)

9 (56%)

Source study, n (%)

A5024

1(13%)

6 (8%)

0 (0%)

A5068

0 (0%)

9 (11%)

3 (19%)

A5170

6 (74%)

46 (58%)

10 (63%)

A5197

1 (13%)

18 (23%)

3 (19%)


Supplementary Table 4. Effect of individual class I alleles on viral control (controllers vs non-controllers) based on a logistic regression model with stepwise selection including all HLA class I alleles with phenotypic frequencies of > 2% (Table S5 in [7]).

White

Black

Allele

P

OR

95% CI

Allele

P

OR

95% CI

A*01:01

2E-11

2.3

1.8-2.9

A*23:01

3E-04

1.9

1.3-2.7

A*02:01

4E-05

1.5

1.2-1.8

A*36:01

4E-04

4.3

1.9-9.8

A*25:01

1E-03

0.5

0.3-0.8

A*03:01

4E-02

0.7

0.5-1.0

A*31:01

5E-03

0.6

0.4-0.9

B*57:03

2E-22

0.1

0.1-0.2

A*68:02

3E-02

0.5

0.3-0.9

B*45:01

8E-05

3.3

1.8-5.9

B*57:01

3E-42

0.1

0.1-0.2

B*35:01

3E-04

2.2

1.4-3.3

B*27:05

3E-13

0.3

0.2-0.4

B*15:10

1E-03

2.6

1.5-4.7

B*07:02

2E-07

2.0

1.5-2.6

B*81:01

2E-03

0.4

0.2-0.7

B*52:01

4E-06

0.3

0.2-0.5

B*58:02

3E-03

2.5

1.4-4.5

B*14:02

2E-05

0.5

0.3-0.7

B*18:01

5E-03

2.5

1.3-4.8

B*13:02

3E-05

0.4

0.3-0.6

B*07:02

1E-02

1.7

1.1-2.5

B*40:01

4E-04

2.0

1.4-2.9

B*14:02

4E-02

0.6

0.3-1.0

B*38:01

5E-04

2.9

1.6-5.2

C*12:03

8E-05

0.3

0.2-0.6

B*18:01

7E-04

2.0

1.3-3.0

C*08:04

4E-03

0.3

0.1-0.7

B*40:02

4E-03

0.5

0.3-0.8

C*05:01

2E-02

0.5

0.3-0.9

B*58:01

2E-02

0.5

0.3-0.9

B*35:02

3E-02

2.9

1.1-7.9

B*55:01

4E-02

2.1

1.0-4.4

C*14:02

1E-04

0.4

0.2-0.6

C*04:01

1E-03

1.5

1.2-2.0


Supplementary Figure 1. Univariate analysis of associations at HLA-B locus. Proportions of participants in HLA-B groups grouped based on their time to viral rebound to ≥ 200 HIV RNA copies/mL.