A Regulatory Polymorphism in HAVCR2 Modulates Susceptibility to HIV-1 Infection

Abstract

The HAVCR2 gene encodes TIM-3, an immunoglobulin superfamily member expressed by exhausted CD8+ T cells during chronic viral infection. We investigated whether genetic variation at HAVCR2 modulates the susceptibility to HIV-1 acquisition; specifically we focused on a 3′ UTR variant (rs4704846, A/G) that represents a natural selection target. We genotyped rs4704846 in three independent cohorts of HIV-1 exposed seronegative (HESN) individuals with different geographic origin (Italy and Spain) and distinct route of exposure to HIV-1 (sexual and injection drug use). Matched HIV-1 positive subjects and healthy controls were also analyzed. In all case-control cohorts the minor G allele at rs4704846 was more common in HIV-1 infected individuals than in HESN, with healthy controls showing intermediate frequency. Results from the three association analyses were combined through a random effect meta-analysis, which revealed no heterogeneity among samples (Cochrane's Q, p value =  0.89, I² =  0) and yielded a p value of 6.8 ×10⁻⁴. The minor G allele at rs4704846 was found to increase HAVCR2 expression after in vitro HIV-1 infection. Thus, a positively selected polymorphism in the 3′ UTR, which modulates HAVCR2 expression, is associated with the susceptibility to HIV-1 infection. These data warrant further investigation into the role of TIM-3 in the prevention and treatment of HIV-1/AIDS.

About the Authors:

Manuela Sironi

Contributed equally to this work with: Manuela Sironi, Mara Biasin

* E-mail: [email protected]

Affiliation: Bioinformatics, Scientific Institute IRCCS E. MEDEA, Bosisio Parini, Italy

Mara Biasin

Contributed equally to this work with: Manuela Sironi, Mara Biasin

Affiliation: Department of Biomedical and Clinical Sciences, University of Milan, Milan, Italy

Federica Gnudi

Affiliation: Department of Biomedical and Clinical Sciences, University of Milan, Milan, Italy

Rachele Cagliani

Affiliation: Bioinformatics, Scientific Institute IRCCS E. MEDEA, Bosisio Parini, Italy

Irma Saulle

Affiliation: Department of Biomedical and Clinical Sciences, University of Milan, Milan, Italy

Diego Forni

Affiliation: Bioinformatics, Scientific Institute IRCCS E. MEDEA, Bosisio Parini, Italy

Veronica Rainone

Affiliation: Department of Biomedical and Clinical Sciences, University of Milan, Milan, Italy

Daria Trabattoni

Affiliation: Department of Biomedical and Clinical Sciences, University of Milan, Milan, Italy

Micaela Garziano

Affiliation: Department of Biomedical and Clinical Sciences, University of Milan, Milan, Italy

Francesco Mazzotta

Affiliation: Infectious Disease Unit, S. Maria Annunziata Hospital Florence, Florence, Italy

Luis Miguel Real

Affiliation: Infectious Diseases and Microbiology Clinical Unit, Valme Hospital, Seville, Spain

Antonio Rivero-Juarez

Affiliation: Maimonides Institut for Biomedical Research (IMIBIC)-Reina Sofia Universitary Hospital-University of Cordoba, Cordoba, Spain

Antonio Caruz

Affiliation: Immunogenetics Unit, Department of Experimental Biology, University of Jaen, Jaen, Spain

Sergio Lo Caputo

Affiliation: Infectious Disease Unit, S. Maria Annunziata Hospital Florence, Florence, Italy

Mario Clerici

Affiliations Chair of Immunology, Department of Physiopathology and Transplantation, University of Milan, Milan, Italy, Don C. Gnocchi Foundation ONLUS, IRCCS, Milan, Italy

Introduction

TIM-3 (T cell immunoglobulin and mucin domain-containing molecule 3) is an immunoglobulin superfamily member encoded in humans by the HAVCR2 (hepatitis A virus cellular receptor 2) gene. Initially identified as a marker of IFN-γ-producing CD4+ Th1 and CD8+ Tc1 cells [1], TIM-3 was more recently shown to be expressed by several other immune cell types including NK/NTK, macrophages/monocytes, and dendritic cells [2]. The protein product of HAVCR2 can bind both galectin 9 (Gal-9) and phosphatidylserine [2]. Engagement of TIM-3 by Gal-9 on T cells induces cell death and promotes peripheral tolerance [2]. Thus, TIM-3 plays an important role in the negative regulation of T-cell mediated responses, and abrogation of its signaling increases the secretion of IFN-γ by activated human T cells [3]. Recent evidences have indicated that expression of TIM-3 marks a population of exhausted CD8+ T cells during chronic viral infection [2]. Specifically, in progressive HIV-1 infection TIM-3 defines an abundant population of CD8+ T cells and its expression correlates positively with viral load and inversely with CD4+ T cell counts [4]. The loss of proliferative activity of HIV-specific TIM-3-expressing CD8+ cells is partially mediated by the interaction with Gal-9 on T_reg cells and is modulated by HLA-B allelic status [5]. Despite these observations, the role of TIM-3 in HIV-1 acquisition has never been analyzed.

We have reported that a variant located in the 3′UTR of HAVCR2 (rs4704846) has been a target of natural selection in human populations and suggested that the selective pressure is accounted for by infectious agents [6]. In line with this view, a SNP (rs3087616) located 62 bp apart and in full linkage disequilibrium with rs4704846 (r² =  1 in Europeans) has recently been shown to act as an expression QTL (eQTL) in CD14+ monocytes [7]. Given the central role of TIM-3 in viral infection [2], and because evolutionary and eQTL analyses point to rs4704846 (or a closely linked variant) as a functional polymorphism, we investigated whether this SNP modulates the susceptibility to HIV-1 infection.

Materials and Methods

Ethics statement

The study was designed and performed according to the Helsinki declaration and was approved by the Ethics Committees of following Institutions: University of Jaen, Valme Hospital (Seville), Reina Sofia Universitary Hospital (Cordoba), and S. Maria Annunziata Hospital (Florence). All patients and healthy blood donors provided written informed consent to participate in this study.

Subject cohorts

Ninety-three Italian HESN that had been exposed to the virus through unprotected sexual intercourse (SexExp-HESN) and 87 HIV-1-infected subjects were recruited at the S. M. Annunziata Hospital in Florence, Italy; all of them were Italian of European origin. Inclusion criteria for HESN were a history of multiple unprotected sexual episodes for more than 4 years at the time of the enrolment, with at least 3 episodes of at-risk intercourse within 4 months prior to study entry, and an average of 30 (range, 18 to >100) reported unprotected sexual contacts per year [8]. All individuals (SexExp-HESN and HIV-1 infected) had been longitudinally followed for >4 years before the study by the Department of Obstetrics and Gynecology of the S. M. Annunziata Hospital. This allowed us to exclude from the study HESN and HIV-1 infected subjects in whom sexually transmitted diseases or any other pathology were reported during that time period. The range of CD4 counts in HIV-1 infected patients were 36–850 cells/ml, and viral loads were >50–750000 copies/ml. All of the patients were receiving highly active antiretroviral therapy (HAART) at the time of the study.

Thirty-eight Spanish HESN that had been exposed to the virus through unprotected sexual intercourse (SexExp-HESN) were recruited as well. These subjects are female partners of HIV-1 infected patients that were treatment-naive and viremic. In this case, mean number of unprotected sexual intercourse per year was 110 and the mean number of years of unprotected sex was 5 (range 3–17 years). Healthy controls (HC, n = 77) were anonymous blood donors from the City of Jaen Hospital in Jaen, Spain.

Finally, we recruited 190 males exposed to HIV-1 infection by injection drug use (IDU) and enrolled in prospective cohort studies in Spain (Valme Hospital, Sevilla) who had shared needles for >3 months. Concurrent markers of hepatitis C virus (HCV) infection, the most chronic viral infection transmitted by sharing needles, were present in 100% of IDU subjects. These values are significantly higher than the reported HCV prevalence of 1%–2% for the general population in Spain. Ninety-three of these subjects were HIV-1 negative (IDU-HESN), 97 were HIV-1 positive (IDU-CTR). The mean of CD4 cells in HIV-1 infected patients was 648±408 mm³, and viral loads were undetectable in 81% of the patients. 87% of these patients were receiving HAART at the time of the study.

All subjects were Spanish of European origin. The main epidemiological characteristics of the populations studied are detailed in Table 1.

[Figure omitted. See PDF.]

Table 1. Characteristic of the study cohorts.

https://doi.org/10.1371/journal.pone.0106442.t001

Genotyping and statistical analysis

Genomic DNA was used as template for PCR amplification using TaqMan probes specifically designed to perform a SNP genotyping assay for rs4704846 (G/A) and using the allelic discrimination real-time PCR method. Assays were performed in 10 µl reactions, using TaqMan Genotyping Master Mix on 96-well plates using a ABI 7000 instrument (Applied Biosystems Foster City, CA, USA). The variant complied to Hardy-Weinberg equilibrium in all samples. Genetic association analyses were performed by logistic regression using an additive model, and results from the three cohorts were combined using a random-effect meta-analysis; all analyses were performed using PLINK [9]. Linkage disequilibrium analyses were performed using Haploview (v. 4.1) [10] and blocks were identified through the confidence interval algorithm implemented in the software [11].

PBMC isolation, HIV-1 infection, and transcript quantification

Whole blood was collected from 40 healthy volunteers by venipuncture in Vacutainer tubes containing EDTA (Becton Dickinson), and PBMCs were separated on lymphocyte separation medium (Organon Teknica, Malvern, PA). PBMCs (10×10⁶ cells/ml) were cultured for 2 days at 37°C and 5% CO2 in RPMI 1640 containing FBS (20%), PHA (7.5 ìg/ml), and IL-2 (15 ng/ml). After viability assessment, 2.5×10⁶ cells were resuspended in medium containing or not containing 1 ng HIV-1Ba-L/1×10⁶ PBMC and incubated for 3 h at 37°C. Cells were then washed and resuspended in 3 ml complete medium with IL-2 (15 ng/ml). Cells were plated in 24-well tissue culture plates and incubated at 37°C and 5% CO2. After 3 days, 2.5×10⁵ PBMCs were collected for gene expression analysis. RNA was extracted from cultured PBMCs and from HIV-1-infected PBMCs by using the acid guanidium thiocyanate–phenol–chloroform method. The RNA was treated with RNase-free DNase (New England Biolabs, Ipswich, MA). One microgram of RNA was reverse transcribed into first-strand cDNA in a 20-ìl final volume containing 1 ìM random hexanucleotide primers, 1 ìM oligonucleotide, and 200 U Moloney murine leukemia virus reverse transcriptase (Promega, Madison, WI). cDNA quantification for HAVCR2 and GAPDH was performed by real-time PCR (DNA Engine Opticon 2; MJ Research, Ramsey, NJ). Reactions were performed using a SYBR Green PCR mix (5 prime, Gaithersburg, MD). Results were expressed as ÄÄCt and presented as ratios between the target gene and the GAPDH housekeeping mRNA.

Results

As mentioned above, rs4704846 is located in the 3′UTR of HAVCR2. The ancestral minor G allele displays a frequency of 0.39, 0.20, and 0.01 in Africans (Yoruba), Europeans and Japanese plus Chinese, respectively, as assessed by the 1000 Genomes Project data [12]. Analysis of linkage disequilibrium (LD) in Europeans along HAVCR2 indicated that rs4704846 lies within an LD block that also includes rs3087616, an eQTL in naive CD14+ monocytes. The variant is also in full LD with several SNPs that were described as eQTLs in lymphoblastoid cell lines [13], [14] (Fig. 1). The 3′UTR of HAVCR2 is moderately conserved in mammals, as assessed through the GERP (Genomic Evolutionary Rate Profiling) score [15] (Fig. 1); prediction of regulatory motifs using rSNPBase [16] indicated that rs4704846 maps to regions showing H3K27ac histone marks (usually associated with active regulatory elements) in CD14+ monocytes and several microRNA binding sites are observed across the entire 3′ UTR (microRNa.org, http://www.microrna.org/microrna/home.do). The presence of regulatory motifs was further analyzed using HaploReg [17]: rs4704846 was found to affect a Smad4 binding site (Fig. 1). Smad4 is a mediator of TGF-beta signaling and its expression after SIV infection differs in rhesus macaque and African green monkeys, which are differentially susceptible to SIV-mediated immunopathology [18].

[Figure omitted. See PDF.]

Figure 1. LD analysis and functional annotation.

The HAVCR2 gene region is shown within the UCSC Genome Browser view. Variants that have been described as eQTLs in monocytes and lymphoblastoid cell lines are in green and blue, respectively. The variant we analyzed is in red. The location of predicted functional elements and microRNA binding sites is also shown. LD analysis is reported in the bottom panel; data refer to Europeans and derive from the 1000 Genomes Projects Phase I data for Europeans (CEU) [12]. r² was calculated with the Haploview software and blocks were identified through the implemented confidence interval algorithm (see methods).

https://doi.org/10.1371/journal.pone.0106442.g001

To explore the role of rs4704846 in HIV-1 acquisition, we genotyped this variant in a well characterized cohort of 93 heterosexual Italian subjects who have a history of unprotected sex with their seropositive partners (sex-exposed HESN, SexExp-HESN). The allele frequency of rs4704846 in these subjects was compared to that observed in a sample of 87 Italian HIV-1 positive individuals. A significant difference was observed, with the minor G allele being much more common in HIV-1 infected (0.25) compared to SexExp-HESN (0.14) (Table 2). The frequency of the G allele in Italians is 0.21, as determined by the 1000 Genomes Phase I Project. We genotyped 300 Italian healthy subjects and obtained a similar frequency of 0.20, which is therefore intermediate between HIV-1 infected and SexExp-HESN. A very similar result was obtained in a second and smaller cohort of SexExp-HESN (n = 38) from Spain: in these subjects the frequency of the G allele was 0.12, much lower than in a sample of 77 Spanish healthy controls (HC, frequency = 0.20); due to the small sample size the association p value did not reach statistical significance (Table 2). Finally, to replicate these results, a third HESN population with a different route of exposure to HIV-1 was analyzed. In particular, we recruited 190 Spanish injection drug users (IDU): all of them were HCV-positive, but, whereas 97 subjects were HIV-1 infected (IDU-HIV-1 infected) as well, the remaining 93 individuals tested HIV-1 negative despite multiple exposures through needle sharing (IDU-HESN). Again, the G allele was more common in IDU-HIV-1 infected compared to IDU-HESN with a borderline significance of 0.064.

[Figure omitted. See PDF.]

Table 2. Association of rs4704846 with HIV-1 infection susceptibility.

https://doi.org/10.1371/journal.pone.0106442.t002

Results from the three association analyses were combined through a random effect meta-analysis, which revealed no heterogeneity among samples (Cochrane's Q p value =  0.89, I² =  0) and yielded a p value of 6.8 ×10⁻⁴ (Table 2). Overall, these results strongly suggest that the minor G allele of rs4704846 is associated with increased risk of HIV-1 acquisition.

To assess whether the modulatory effect of rs4704846 on HIV-1 susceptibility is exerted through changes in HAVCR2 expression, we performed an in vitro infection assay. Specifically, PBMCs from 40 healthy subjects (27 AA, 13 AG) were infected with HIV-1 and HAVCR2 transcript abundance was quantified by real-time PCR after 3 days. As shown in figure 2, significantly higher HAVCR2 expression was observed following in vitro HIV-1 infection in AG heterozygotes compared to AA homozygotes (Student's t-test, p = 0.028) (Fig. 2); the same trend was observed in uninfected PBMCs, although the difference was less marked and did not reach statistical significance (Fig. 2).

[Figure omitted. See PDF.]

Figure 2. Box-and-whisker plot of HAVCR2 expression depending on rs4704846 genotype.

Data derive from PBMCs from 40 healthy volunteers uninfected or infected with HIV-1. HAVCR2 transcript levels are log-transformed and shown in standard box-and-whisker plot representation (thick line: median; box: quartiles; whiskers: 1.5 × interquartile range); p values are calculated using the Student's t-test.

https://doi.org/10.1371/journal.pone.0106442.g002

Discussion

We analyzed three independent HESN cohorts and in all cases we observed the same trend: HIV-1 infected subjects are more likely to carry the G allele than HESN, with healthy controls showing intermediate frequency. Therefore, although statistical analysis was fully significant for the Italian cohort only, the meta-analysis we performed revealed no heterogeneity in the effect of rs4704846 among the three samples and yielded a highly significant association result. Indeed, the small sample size of HESN cohorts, which is conceivable given the characteristics of these subjects, make replication across populations and meta-analysis a powerful approach to discover susceptibility variants for HIV-1 infection. We should add that the single case-control samples used either HIV-1 infected or HC subjects as a comparison to HESN; nonetheless, we believe meta-analysis to be appropriate as HIV-1 infection susceptibility is a common condition to most humans [19], [20].

The interaction between TIM-3 and Gal-9 acts to limit the extent of CD8⁺ T cell immunity to HSV infection [21] and over-expression of TIM-3 on CD4+ and CD8+ T cells correlates with diseases progression in chronic hepatitis B infection [22]. Likewise, the chronic persistence of HIV-1 is associated with the increased expression of TIM-3 on CD4+ and virus-specific CD8+ T cells [4]. Results herein fit within this scenario by showing that the expression level of HAVCR2 following infection is at least partially determined by a polymorphism that also affects susceptibility to HIV-1. Although protein and RNA levels do not correlate perfectly in humans[23], it is conceivable that allelic status at rs4704846 also affects TIM-3 abundance in PBMCs.

Natural selection targets variants with a phenotypic effect and acts in response to specific selective pressures. The major allele of rs4704846, which associates with lower susceptibility to HIV-1 acquisition, has been driven to high frequency in human populations by natural selection [6]. Because its appearance as a human pathogen is recent, HIV-1 cannot be regarded as the underlying selective pressure. Nonetheless, as mentioned above, TIM-3 dysregulation has been associated with chronic and acute infections with other viral species [4], [22], [24], suggesting that extant or extinct pathogens drove the frequency increase of the protective allele. The location of rs4704846 in the 3′ UTR and its full linkage-disequilibrium with eQTLs in monocytes and lymphoblastoid cell lines are in agreement with the effect we observed on HAVCR2 expression; nonetheless, the underlying molecular mechanism(s) (miRNA-mediated regulation, transcription factor binding site alteration or other) remains to be determined. Although replication in additional independent samples will be necessary, data herein warrant further investigation into the role of HAVCR2 in the prevention and treatment of HIV-1/AIDS.

Author Contributions

Conceived and designed the experiments: MS MB MC. Performed the experiments: FG RC IS DF VR MG. Analyzed the data: MS MB DT. Contributed reagents/materials/analysis tools: FM LMR ARJ AC SLC. Contributed to the writing of the manuscript: MS MB MC.

Citation: Sironi M, Biasin M, Gnudi F, Cagliani R, Saulle I, Forni D, et al. (2014) A Regulatory Polymorphism in HAVCR2 Modulates Susceptibility to HIV-1 Infection. PLoS ONE 9(9): e106442. https://doi.org/10.1371/journal.pone.0106442

References

1. Monney L, Sabatos CA, Gaglia JL, Ryu A, Waldner H, et al. (2002) Th1-specific cell surface protein tim-3 regulates macrophage activation and severity of an autoimmune disease. Nature 415: 536–541.

2. Freeman GJ, Casasnovas JM, Umetsu DT, DeKruyff RH (2010) TIM genes: A family of cell surface phosphatidylserine receptors that regulate innate and adaptive immunity. Immunol Rev 235: 172–189.

3. Koguchi K, Anderson DE, Yang L, O'Connor KC, Kuchroo VK, et al. (2006) Dysregulated T cell expression of TIM3 in multiple sclerosis. J Exp Med 203: 1413–1418.

4. Jones RB, Ndhlovu LC, Barbour JD, Sheth PM, Jha AR, et al. (2008) Tim-3 expression defines a novel population of dysfunctional T cells with highly elevated frequencies in progressive HIV-1 infection. J Exp Med 205: 2763–2779.

5. Elahi S, Dinges WL, Lejarcegui N, Laing KJ, Collier AC, et al. (2011) Protective HIV-specific CD8+ T cells evade treg cell suppression. Nat Med 17: 989–995.

6. Forni D, Cagliani R, Pozzoli U, Colleoni M, Riva S, et al. (2013) A 175 million year history of T cell regulatory molecules reveals widespread selection, with adaptive evolution of disease alleles. Immunity 38: 1129–1141.

7. Fairfax BP, Humburg P, Makino S, Naranbhai V, Wong D, et al. (2014) Innate immune activity conditions the effect of regulatory variants upon monocyte gene expression. Science 343: 1246949.

8. Miyazawa M, Lopalco L, Mazzotta F, Lo Caputo S, Veas F, et al. (2009) The ‘immunologic advantage’ of HIV-exposed seronegative individuals. AIDS 23: 161–175.

9. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, et al. (2007) PLINK: A tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 81: 559–575.

10. Barrett JC, Fry B, Maller J, Daly MJ (2005) Haploview: Analysis and visualization of LD and haplotype maps. Bioinformatics 21: 263–265 10.1093/bioinformatics/bth457.

11. Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, et al. (2002) The structure of haplotype blocks in the human genome. Science 296: 2225–2229 10.1126/science.1069424.

12. Genomes Project Consortium, Abecasis GR, Auton A, Brooks LD, DePristo MA, et al (2012) An integrated map of genetic variation from 1,092 human genomes. Nature 491: 56–65.

13. Lappalainen T, Sammeth M, Friedlander MR, 't Hoen PA, Monlong J, et al. (2013) Transcriptome and genome sequencing uncovers functional variation in humans. Nature 501: 506–511.

14. Grundberg E, Small KS, Hedman AK, Nica AC, Buil A, et al. (2012) Mapping cis- and trans-regulatory effects across multiple tissues in twins. Nat Genet 44: 1084–1089.

15. Cooper GM, Stone EA, Asimenos G (2005) NISC Comparative Sequencing Program, Green ED, et al (2005) Distribution and intensity of constraint in mammalian genomic sequence. Genome Res 15: 901–913.

16. Guo L, Du Y, Chang S, Zhang K, Wang J (2014) rSNPBase: A database for curated regulatory SNPs. Nucleic Acids Res 42: D1033–9.

17. Ward LD, Kellis M (2012) HaploReg: A resource for exploring chromatin states, conservation, and regulatory motif alterations within sets of genetically linked variants. Nucleic Acids Res 40: D930–4.

18. Ploquin MJ, Desoutter JF, Santos PR, Pandrea I, Diop OM, et al. (2006) Distinct expression profiles of TGF-beta1 signaling mediators in pathogenic SIVmac and non-pathogenic SIVagm infections. Retrovirology 3: 37.

19. Plummer FA, Ball TB, Kimani J, Fowke KR (1999) Resistance to HIV-1 infection among highly exposed sex workers in nairobi: What mediates protection and why does it develop? Immunol Lett 66: 27–34.

20. Fowke KR, Nagelkerke NJ, Kimani J, Simonsen JN, Anzala AO, et al. (1996) Resistance to HIV-1 infection among persistently seronegative prostitutes in nairobi, kenya. Lancet 348: 1347–1351.

21. Reddy PB, Sehrawat S, Suryawanshi A, Rajasagi NK, Mulik S, et al. (2011) Influence of galectin-9/Tim-3 interaction on herpes simplex virus-1 latency. J Immunol 187: 5745–5755.

22. Wu W, Shi Y, Li J, Chen F, Chen Z, et al. (2011) Tim-3 expression on peripheral T cell subsets correlates with disease progression in hepatitis B infection. Virol J 8: 113–422X-8-113.

23. Wu L, Candille SI, Choi Y, Xie D, Jiang L, et al. (2013) Variation and genetic control of protein abundance in humans. Nature 499: 79–82.

24. Golden-Mason L, Palmer BE, Kassam N, Townshend-Bulson L, Livingston S, et al. (2009) Negative immune regulator tim-3 is overexpressed on T cells in hepatitis C virus infection and its blockade rescues dysfunctional CD4+ and CD8+ T cells. J Virol 83: 9122–9130.