Variants in the CYP7B1 gene region do not affect natural resistance to HIV-1 infection
- Manuela Sironi†1,
- Mara Biasin†2Email author,
- Chiara Pontremoli1,
- Rachele Cagliani1,
- Irma Saulle2,
- Daria Trabattoni2,
- Francesca Vichi3,
- Sergio Lo Caputo3,
- Francesco Mazzotta3,
- Wbeimar Aguilar-Jimenez4,
- Maria Teresa Rugeles4,
- Samandhy Cedeno5, 6,
- Jorge Sanchez7,
- Christian Brander5, 6 and
- Mario Clerici8, 9
© Sironi et al. 2015
Received: 21 July 2015
Accepted: 14 September 2015
Published: 24 September 2015
The genetic bases of natural resistance to HIV-1 infection remain largely unknown. Recently, two genome-wide association studies suggested a role for variants within or in the vicinity of the CYP7B1 gene in modulating HIV susceptibility. CYP7B1 is an appealing candidate for this due to its contribution to antiviral immune responses. We analyzed the frequency of two previously described CYP7B1 variants (rs6996198 and rs10808739) in three independent cohorts of HIV-1 infected subjects and HIV-1 exposed seronegative individuals (HESN).
rs6996198 and rs10808739 were genotyped in three case/control cohorts of sexually-exposed HESN and HIV-1-infected individuals from Italy, Peru and Colombia. Comparison of the allele and genotype frequencies of the two SNPs under different models showed that the only significant difference was seen for rs6996198 in the Peruvian sample (nominal p = 0.048, dominant model). For this variant, a random-effect meta-analysis yielded non-significant results (dominant model, p = 0.78) and revealed substantial heterogeneity among cohorts. No significant effect of the rs10808739 allelic status on HIV-1 infection susceptibility (additive model, p = 0.30) emerged from the meta-analysis.
Although our study had limited power to detect association due to the small sample size, comparisons among the three cohorts revealed very similar allelic and genotypic frequencies in HESN and HIV-1 positive subjects. Overall, these data indicate that the two GWAS-defined variants in the CYP7B1 region do not strongly influence HIV-1 infection susceptibility.
Individual variation in susceptibility to HIV-1 infection is well documented [1–3], and a minority of individuals (15 % of the HIV exposed people), usually referred to as HIV-1 exposed seronegative individuals (HESN), does not become productively infected despite multiple exposures to the virus [2, 4, 5]. Because only a minority of these subjects is homozygous for the CCR5Δ32 deletion , it is logical to postulate that other, additional genetic factors can modulate susceptibility to HIV-1 infection. Some of such factors were identified through candidate gene approaches; only few of them however were replicated across different samples . In recent years, genome-wide association studies (GWAS) have been used for the identification of common variants that underlie complex phenotypes. Among these, two GWASs for HIV-1 infection susceptibility have suggested a role of variants within or in the vicinity of the CYP7B1 gene in this phenomenon [7, 8]. In one study, HIV-1 infected individuals were compared with HIV-1 negative subjects . Using meta-analysis of two European cohorts and further validation in European Americans, a single signal with genome-wide significance was detected at rs6996198, a SNP that is located downstream the CYP7B1 gene and possibly modulates its expression . Another variant in CYP7B1 (rs10808739, in intron 1) was described in a GWAS of HIV-1 serodiscordant couples from Eastern and Southern Africa. Although the variant did not reach genome-wide significance, it was one of the top signals in the study .
The CYP7B1 gene encodes an enzyme belonging to the cytochrome P450 superfamily that is involved in cholesterol catabolism by inactivation of oxysterols and their subsequent conversion into bile salts, mainly in the liver, brain, and reproductive tract . Controlling intracellular cholesterol metabolism is essential for the correct entry, assembly and budding of virions, and many viruses are dependent on cellular lipid metabolisms to ensure the correct budding of functional virions from infected cells . HIV alters cholesterol trafficking switching it from a physiological efflux to a virus-controlled transport, thus reducing the ability of a cell to export cholesterol. CYP7B1 also modulates several immune functions, as well including proinflammatory cytokine release  and, via its catabolite 25-hydroxycholesterol , programmed cell death  and the synthesis of IgA , a class of antibodies that has been associated with resistance to HIV infection . Hence, CYP7B1 represents an appealing candidate as a modifier of infection susceptibility due to its contribution to the synthesis of virions and the initiation of antiviral defense.
To verify the role played by CYP7B1 in resistance to HIV infection we recruited 125 Italian HESN exposed through unprotected sexual intercourse (SexExp-HESN). Inclusion criteria were a history of multiple unprotected sexual episodes for more than 4 years at the time of enrolment, with at least three 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 . Infection in HESN was ruled-out by plasma HIV RNA and proviral DNA analyses. HESN and 114 seropositive (SP) partners were recruited at the S. M. Annunziata Hospital, Florence; all of them were Italian of Caucasian origin.
Sixty-two Colombian SexExp-HESN and 51 SP partners were also included. The inclusion criteria for these HESN subjects were previously reported  and included a negative HIV-1/2 ELISA test within 1 month of sample taking. The similar ancestry component and pair-wise fixation index (FST) values in the Colombian cohort  indicated no intra-cohort stratification by ethnicity.
The third HESN cohort was recruited in Peru and has been described [18, 19] both in terms of host genetics (HLA and KIR) and immune reactivity to HIV and viral co-pathogens. The similar frequency of HLA and KIR alleles  suggested no major intra-cohort stratification. For this study, samples were available for a total of 130 HESN and 95 SP individuals, all of whom were recruited at IMPACTA clinics across Lima. HESN were tested on a 3-monthly basis for newly acquired HIV infection. Risk criteria for the HESN cohort were more than five different sexual partners over the last 3 months, reported STI over the last 6 months, sexual intercourse with a known HIV SP partner in the last 6 months and having accepted money for sex as described previously [J. Coll et al. Bangkok AIDS Vaccine meeting 2011, Factors influencing recruitment, retention and seroconversion rates in MSM at high risk for HIV infection in Lima and Barcelona].
The study was designed and performed according to the Helsinki declaration (1975 revised in 2000) and was approved by the Ethics Committee of the participating units. All subjects provided written informed consent to participate in this study.
Genotyping of CYP7B1 rs6996198 and rs10808739 was performed by TaqMan probe assays (TaqMan SNP genotyping assay, Applied Biosystems, Foster City, CA, USA) using the allelic discrimination real-time PCR method. Genotyping rate was >0.90 for both SNPs in all samples. The two SNPs complied to Hardy–Weinberg equilibrium in all samples with the only exception of rs6966198 in the SP Peruvian sample (uncorrected p value = 0.02). Association analysis for single variants was performed using logistic regressions under different models using the PLINK software . Meta-analysis was performed using a random effect model using the R package “meta” (R package version 4.3-0. http://CRAN.R-project.org/package=meta). Linkage analysis was performed using HaploReg v2 .
Demographic status of the populations examined
SP (n = 51)
HESN (n = 62)
SP (n = 114)
HESN (n = 125)
SP (n = 95)
HESN (n = 130)
Age, mean years ± SD
33.9 ± 7.5
35.1 ± 10.6
42.4 ± 8.8
41.6 ± 9.1
30.8 ± 6.7
31.2 ± 10.7
Males, n (%)
Viral load, median copies/mL (interquartile range)
CD4+ T cell/µL count, median (interquartile range)
Monthly unprotected sexual episodes, mean (range)c
Previous history of sexually transmitted diseases and/or AIDS-defining illnesses (%)
Heterosexual orientation (%)
Homosexual orientation (%)
Bisexual orientation (%)
Ethnicity, ancestryb (%)
European (Tuscan): 100
European (Tuscan): 100
Association of rs6966198 and rs10808739 with HIV-1 infection susceptibility
p additive b
p dominant c
p recessive d
CT + TTa
p additive b
p dominant c
p recessive d
CYP7B1 plays a major role in the synthesis of HIV and in immunity and is thus a likely candidate in the modulation of susceptibility to HIV. Results, however, did not confirm the previously reported associations, suggesting that these particular CYP7B1 genetic variants do not play a role in HIV-1 infection susceptibility, and indicating that other variants will need to be analyzed. A limitation of our study is the relatively small sample sizes of the HESN and HIV-1 infected cohorts, resulting in low power to detect potentially existing associations. Nonetheless, the combined sample size was comparable to that analyzed in the serodiscordant couple study that detected an association at rs10808739 . As for rs6966198, it reached genome-wide significance in a relatively large two-stage study that employed a different design from the one applied herein . Thus, in that case, HIV-1-infected individuals were compared to HIV-1 negative subjects, irrespective of their exposure status. This may have determined a proportion of HIV-1 negative controls to be misclassified—i.e. most of with a consequent reduction of power. Our analyses were performed in well characterized HESN cohorts hence overcoming this problem; a degree of misclassification may nevertheless still remain for HIV-1 positive subjects because their exposure history may not always be properly and completely recorded and their group is likely to include individuals in a wide susceptibility range.
Taking into account these limitations, comparisons in the three cohorts showed very similar allelic and genotypic frequencies in HESN and HIV-1-infected individuals, with the only exception being the Peruvian sample that reached nominal significance for rs6996198. For the same variant, the Italian cohort also showed some difference in frequency, albeit not significant, that was in the opposite direction to that reported by the original association study  and to that observed in Peruvians. The original study that identified rs6996198 as potentially associated with HIV-1 infection susceptibility estimated an OR of 0.64 with narrow confidence intervals (95 % IC: 0.54–075). Herein, also due to heterogeneity among samples, we obtained an OR very close to 1 with wide uncertainty (95 % IC: 0.5727–2.1016). Thus, our results do not imply that rs6996198 has no role in modulating HIV-1 infection susceptibility, but indicate that its effect (if any) is small. Overall, these data lead us to conclude that the two genotyped variants in the CYP7B1 region do not strongly influence HIV-1 infection susceptibility.
MS, MB and MC, conceived the study and wrote the paper; IS, RC and CP performed the experiments and analyzed the data; SLC, FM, FV and DT enrolled and selected the cohort of HESN subjects in Florence; WA-J, and MTR enrolled and selected the cohort of HESN subjects in Colombia; CB enrolled and selected the cohort of HESN subjects in Peru. All authors read and approved the final manuscript.
This work was supported by: Estrategia de Sostenibilidad 2014–2015 de la Universidad de Antioquia; NIH-NIDCR R01 DE018925-04; the HIVACAT program and CUTHIVAC 241904. Thanks to “Grupo de Vida” foundation in Medellin—Colombia, to the staff of HERES Health, from Santa Marta—Colombia for their collaboration in patient recruitment.
Compliance with ethical guidelines
Competing interests The authors declare that they have no competing interests.
Open AccessThis 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.
- Plummer FA, Ball TB, Kimani J, Fowke KR. Resistance to HIV-1 infection among highly exposed sex workers in Nairobi: what mediates protection and why does it develop? Immunol Lett. 1999;66(1–3):27–34.View ArticlePubMedGoogle Scholar
- Fowke KR, Nagelkerke NJ, Kimani J, Simonsen JN, Anzala AO, Bwayo JJ, et al. Resistance to HIV-1 infection among persistently seronegative prostitutes in Nairobi, Kenya. Lancet. 1996;348(9038):1347–51.View ArticlePubMedGoogle Scholar
- Piacentini L, Biasin M, Fenizia C, Clerici M. Genetic correlates of protection against HIV infection: the ally within. J Intern Med. 2009;265(1):110–24.View ArticlePubMedGoogle Scholar
- Clerici M, Giorgi JV, Chou CC, Gudeman VK, Zack JA, Gupta P, et al. Cell mediated immune response to human immunodeficiency virus type 1 (HIV-1) in seronegative homosexuals with recent sexual exposure to HIV-1. JID. 1992;165:1012–9.View ArticlePubMedGoogle Scholar
- Miyazawa M, Lopalco L, Mazzotta F, Lo Caputo S, Veas F, Clerici M. The, “immunologic advantage” of HIV-exposed seronegative individuals. AIDS. 2009;23:161–76.View ArticlePubMedGoogle Scholar
- McLaren PJ, Fellay J. Human genetic variation in HIV disease: beyond genome-wide association studies. Curr Opin HIV AIDS. 2015;10(2):110–5.View ArticlePubMedGoogle Scholar
- Limou S, Delaneau O, van Manen D, An P, Sezgin E, Le Clerc S, et al. Multicohort genomewide association study reveals a new signal of protection against HIV-1 acquisition. J Infect Dis. 2012;205(7):1155–62.PubMed CentralView ArticlePubMedGoogle Scholar
- Lingappa JR, Petrovski S, Kahle E, Fellay J, Shianna K, McElrath MJ, et al. Genomewide association study for determinants of HIV-1 acquisition and viral set point in HIV-1 serodiscordant couples with quantified virus exposure. PLoS One. 2011;6(12):e28632.PubMed CentralView ArticlePubMedGoogle Scholar
- Yantsevich AV, Dichenko YV, Mackenzie F, Mukha DV, Baranovsky AV, Gilep AA, et al. Human steroid and oxysterol 7α-hydroxylase CYP7B1: substrate specificity, azole binding and misfolding of clinically relevant mutants. FEBS J. 2014;281(6):1700–13.View ArticlePubMedGoogle Scholar
- Chukkapalli V, Heaton NS, Randall G. Lipids at the interface of virus–host interactions. Curr Opin Microbiol. 2012;15:512–8.PubMed CentralView ArticlePubMedGoogle Scholar
- Dulos J, van der Vleuten MA, Kavelaars A, Heijnen CJ, Boots AM. CYP7B expression and activity in fibroblast-like synoviocytes from patients with rheumatoid arthritis: regulation by proinflammatory cytokines. Arthritis Rheum. 2005;52(3):770–8.View ArticlePubMedGoogle Scholar
- Martin C, Bean R, Rose K, Habib F, Seckl J. Cyp7b1 catalyses the 7alpha-hydroxylation of dehydroepiandrosterone and 25-hydroxycholesterol in rat prostate. Biochem J. 2001;355(Pt 2):509–15.PubMed CentralView ArticlePubMedGoogle Scholar
- Ayala-Torres S, Moller PC, Johnson BH, Thompson EB. Characteristics of 25-hydroxycholesterol-induced apoptosis in the human leukemic cell line CEM. Exp Cell Res. 1997;235(1):35–47.View ArticlePubMedGoogle Scholar
- Bauman DR, Bitmansour AD, McDonald JG, Thompson BM, Liang G, Russell DW. 25-Hydroxycholesterol secreted by macrophages in response to Toll-like receptor activation suppresses immunoglobulin A production. Proc Natl Acad Sci USA. 2009;106(39):16764–9.PubMed CentralView ArticlePubMedGoogle Scholar
- Broliden K, Hinkula J, Devito C, Kiama P, Kimani J, Trabbatoni D, et al. Functional HIV-1 specific IgA antibodies in HIV-1 exposed, persistently IgG seronegative female sex workers. Immunol Lett. 2001;79(1–2):29–36.View ArticlePubMedGoogle Scholar
- Zapata W, Rodriguez B, Weber J, Estrada H, Quiñones-Mateu ME, Zimermman PA, et al. Increased levels of human beta-defensins mRNA in sexually HIV-1 exposed but uninfected individuals. Curr HIV Res. 2008;6(6):531–8.PubMed CentralView ArticlePubMedGoogle Scholar
- Zapata W, Aguilar-Jiménez W, Pineda-Trujillo N, Rojas W, Estrada H, Rugeles MT. Influence of CCR5 and CCR2 genetic variants in the resistance/susceptibility to HIV in serodiscordant couples from Colombia. AIDS Res Hum Retroviruses. 2013;29(12):1594–603.PubMed CentralView ArticlePubMedGoogle Scholar
- Olvera A, Pérez-Álvarez S, Ibarrondo J, Ganoza C, Lama JR, Lucchetti A, et al. The HLA-C*04: 01/KIR2DS4 gene combination and human leukocyte antigen alleles with high population frequency drive rate of HIV disease progression. AIDS. 2015;29(5):507–17.PubMedGoogle Scholar
- Xu GJ, Kula T, Xu Q, Li MZ, Vernon SD, Ndung’u T, et al. Comprehensive serological profiling of human populations using a synthetic human virome. Science. 2015;348(6239):aaa0698.View ArticlePubMedGoogle Scholar
- Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81:559–75.PubMed CentralView ArticlePubMedGoogle Scholar
- Ward LD, Kellis M. HaploReg: a resource for exploring chromatin states, conservation, and regulatory motif alterations within sets of genetically linked variants. Nucleic Acids Res. 2012;40(Database issue):D930–4.PubMed CentralView ArticlePubMedGoogle Scholar