PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of jacLink to Publisher's site
 
J Antimicrob Chemother. 2010 September; 65(9): 1889–1893.
Published online 2010 July 17. doi:  10.1093/jac/dkq260
PMCID: PMC2920179

Haplotype structure of CYP2B6 and association with plasma efavirenz concentrations in a Chilean HIV cohort

Abstract

Objectives

Efavirenz is extensively metabolized by CYP2B6, and associations between CYP2B6 polymorphisms and plasma efavirenz exposure have been reported. The objective of this study was to investigate CYP2B6 haplotype structure and functional consequences in a Latin American population.

Patients and methods

Two hundred and nineteen patients were recruited at Fundación Arriarán, Chile, between September and December 2008. Plasma efavirenz concentrations were determined using liquid chromatography with mass spectrometry. Genotyping for 30 single nucleotide polymorphisms (SNPs) with a minor allele frequency of >0.05 in the HapMap CEU population at intervals of ~1 kb across the CYP2B6 locus was conducted using Sequenom iPLEX MALDI-TOF.

Results

Thirteen SNPs passed quality control and, of these, statistically significant associations (P < 0.001) with plasma efavirenz concentrations were observed for 11. Pairwise tagging SNP analysis (R2 > 0.8) identified 3 SNPs (rs10403955, rs2279345 and rs8192719) representative of the 11 associated SNPs. A composite genetic model of these three alleles was constructed, and an association between carriers of four to six of these alleles and the risk of efavirenz plasma concentrations >4 µg/mL was identified with an odds ratio of 48.1 (95% confidence interval: 13.5–207.7). This represents a positive predictive value of 80.9% and a negative predictive value of 91.8%, with sensitivity of 57.9% and specificity of 97.2%.

Conclusions

A composite genetic model of CYP2B6 SNPs in a Chilean HIV-positive cohort may have value in predicting concentrations of efavirenz associated with a higher likelihood of CNS toxicity. Further investigation of the functional basis of these associations is now required.

Keywords: antiretrovirals, pharmacokinetics, genetics, single nucleotide polymorphisms

Introduction

Efavirenz is a first-line therapy for the treatment of HIV globally. Despite the potency and tolerability of efavirenz-containing regimens, some patients can fail to achieve durable viral suppression or develop treatment-limiting toxicity.1,2 Efavirenz plasma concentrations >4 µg/mL have been associated with an increased risk of CNS adverse effects, while concentrations <1 µg/mL have been associated with virological failure.3 Although other enzymes are involved, efavirenz is predominantly metabolized by CYP2B6.4

The most commonly studied CYP2B6 single nucleotide polymorphism (SNP) is c.516G→T (rs3745274), which encodes a Gln172His amino acid substitution. This polymorphism is associated with lower hepatic expression and enzymatic activity for CYP2B6.5,6 Many studies have shown an association between the c.516G→T allele and higher efavirenz exposure in different ethnic groups.710 In addition to c.516G→T, CYP2B6 polymorphisms, such as the c.983T→C SNP, have also been shown to impact on plasma efavirenz concentrations in black populations.1113

To date, there is little information on the frequency and functional consequences of CYP2B6 polymorphisms in Latin American populations. Therefore, the aim of this study was to investigate the CYP2B6 haplotype structure in a Chilean HIV cohort and to identify associations with efavirenz plasma concentrations. The ability of CYP2B6 genotypes to predict individuals at risk of high efavirenz plasma concentrations was also assessed.

Methods

Patients

Two hundred and nineteen HIV-positive patients were recruited at Fundación Arriarán, Santiago, Chile, between September and December 2008. Of these individuals, 214 were self-described white or mestizos (mixed European and Amerindian ancestry), with 5 of native ancestry. All individuals were >18 years of age and taking 600 mg of efavirenz once daily for ≥2 weeks. The median [interquartile range (IQR)] duration of therapy was 3.6 years (1.6–5.8 years). The median (IQR) time since last dose was 11.7 h (10.8–12.7 h) and, as such, all efavirenz plasma concentration measurements can be considered as Ctrough. Patients were excluded if informed consent could not be obtained or if rifampicin was used as a co-medication. Ethical approval for the study was obtained from the Comité de Ética de Investigación en Seres Humanos de la Facultad de Medicina Universidad de Chile.

Plasma efavirenz analysis

Sample preparation consisted of protein precipitation with 0.1% acetic acid in methanol followed by centrifugation at 2000 g for 5 min. An aliquot of 250 µL of supernatant was used and 100 µL of an internal standard consisting of 2 mg/L 6,7-dimethyl-2,3-di(2-pyridyl)-quinoxaline (Aldrich, Milwaukee, WI, USA) was added. All solvents and chemicals were obtained from Sigma (Oakville, ON, Canada) and were of HPLC grade. The liquid chromatography with mass spectrometry system consisted of a HP1100 LC system (Agilent Technologies, Wilmington, DE, USA) with a Supelcosil™ ABZ+ (15 cm × 4.6 mm, 3 µm) C18 column (Supelco, Bellefonte, PA, USA) coupled to an API-2000 mass spectrometer (AB/MDS/Sciex, Concord, ON, Canada) with a turbo ion-spray source. The efavirenz calibration curve ranged from 0.05 to 5 mg/L. The intraday accuracy and precision were −2.6%–2.0% and 5.4%–11.2%, respectively. The interday accuracy and precision were −5.3% to −0.9% and 8.4%–13.0%, respectively. More details on this assay can be found elsewhere.14

SNP selection

SNPs with a minor allele frequency (MAF) of >0.05 were identified within dbSNP (http://www.ncbi.nlm.nih.gov/SNP) at intervals of ~1 kb across the CYP2B6 locus at chr19q13.2 and an additional 5 kb up- and downstream of the gene. Where no SNPs with validated MAFs were available, the nearest positional SNP was included.

DNA extraction and genotyping

DNA was extracted using a QIAamp whole blood mini kit (Qiagen, UK), according to the manufacturer's protocol. Thirteen DNA samples were of insufficient quality or quantity to acquire reliable genotype data and so were omitted. Therefore, a total of 206 individuals were genotyped for 32 SNPs using the MALDI-TOF-based Sequenom iPLEX system (Sequenom Inc., San Diego, CA, USA), according to the manufacturer's protocol. Sequence-specific PCR and extend reaction oligonucleotides were obtained from Metabion GmbH (Martinsried, Germany). Multiplex assays were designed using the software available at https://mysequenom.com/default.aspx. SNPs were divided into a 21- and an 11-plex assay.

Statistical analysis

SNPs with MAF <0.05 or call rate <90% were omitted. Statistical analysis was undertaken using SPSS 17.0. Continuous variable analysis by genotype was by one-way analysis of variance (ANOVA). Binary outcomes by genotype were analysed by Fisher's exact test. P < 0.05 was regarded as significant. Tagging SNP analysis was carried out using the Tagger function within Haploview 4.1 software (www.broad.mit.edu/mpg/haploview/).

Results

Association of individual CYP2B6 SNPs with plasma efavirenz concentrations

Thirteen SNPs were carried forward for further analysis (Table 1) and all were in Hardy–Weinburg equilibrium. Ten SNPs were rejected for MAF <0.05 and seven for call rate <90%. One-way ANOVA of efavirenz plasma concentration with genotype for individual SNPs indicated statistically significant associations (P < 0.001) for 11/13 SNPs within the CYP2B6 gene. The exceptions were rs4802100 (5′ upstream) (P = 0.08) and rs34083050 (intron 4) (P = 0.91). The statistical significance of the genotype associations with efavirenz plasma concentration ranged from P = 2.6 × 10−8 (rs2279345) to P = 3.6 × 10−22 (rs8192719). Only one associated SNP was exonic [c.516G→T (rs3745274)] (P = 5.6 × 10−20).

Table 1
Genotype and allelic frequencies of CYP2B6 single nucleotide polymorphisms and median plasma efavirenz concentrations

CYP2B6 haplotype structure in Chilean HIV-positive population

Linkage disequilibrium (LD) analysis of the 13 CYP2B6 polymorphisms showed a significant degree of LD (Figure 1a). Pairwise tagging SNP analysis (R2 > 0.8) identified 3 SNPs (rs10403955, rs2279345 and rs8192719) that were representative of the 11 SNPs associated with plasma efavirenz concentration. These three SNPs were therefore used to construct a composite genetic model. Inferred haplotype analysis of the three SNPs showed that the GCT haplotype was more frequent in Chilean patients (0.34) than in the HapMap (http://hapmap.ncbi.nlm.nih.gov/) Caucasian population (0.27), but less frequent than in the HapMap Yoruba population (0.40). No data were available for rs10403955 in the HapMap Mexican ancestry in Los Angeles cohort and, as such, a direct comparison of inferred haplotype frequency was not possible. However, no significant differences were observed between Chilean patients and the HapMap Mexican population for the frequency of other genotyped SNPs.

Figure 1
LD and association of composite CYP2B6 genotype with efavirenz plasma concentrations. (a) Data are shown as represented in Haploview software. White squares, D < 1 and logarithm of odds (LOD) < 2; ...

Association of CYP2B6 composite genetic model with plasma efavirenz concentrations

The number of alleles associated with high efavirenz concentration (rs10403955 G, rs2297345 C and rs8192719 T) possessed by each individual was tallied and analysed in a composite genetic model (Figure 1b). Individuals possessing zero of six associated alleles (n = 69) had a median (range) efavirenz plasma concentration of 2.26 µg/mL (0.13–4.43 µg/mL). Median (range) plasma efavirenz concentrations were significantly higher in those with one to three alleles (n = 87) at 2.76 µg/mL (0.72–9.18 µg/mL) (P < 0.01). At 5.32 µg/mL (1.37–10.1 µg/mL), individuals with four to six associated alleles had a significantly higher median efavirenz plasma concentration than either the zero or one to three allele groups (P < 0.001 in both instances).

An association between carriers of four to six risk alleles with the risk of an efavirenz plasma concentration greater than the proposed minimal toxic concentration (MTC) (4 µg/mL) was also identified. Of the 38 individuals with an efavirenz plasma concentration >MTC, 22 possessed four to six associated alleles compared with 4/140 of those with an efavirenz plasma concentration <MTC. This represents an odds ratio of 48.1 (95% confidence interval: 13.5–207.7) (P < 0.0001). If the prevalence of >MTC is taken to be 17% (38/219) then the positive predictive value for this association was 80.9% and the negative predictive value was 91.8%, with a sensitivity of 57.9% and a specificity of 97.2%. The positive and negative predictive values of homozygous (TT) carriers of the top associated SNP (rs8192719) were 80.8% and 91.1%, though the sensitivity was marginally lower than the composite genetic model at 53.6%.

Discussion

Many studies have identified associations between polymorphisms in CYP2B6 and efavirenz pharmacokinetics.710 The data reported show that strong LD exists across the CYP2B6 locus. Our report is the first to show the haplotype structure of CYP2B6 in Chile. The frequency of indigenous Chilean ethnicity in our cohort (2.3%) is comparable to the Chilean 2002 census (www.ine.cl/cd2002/sintesiscensal.pdf; 4.6%), suggesting that our cohort is ethnically representative. Given that >97% of the study population are of white European or mestizo ancestry (mixed European and Amerindian), comparisons with Los Angeles residents with Mexican ancestry are valid. Allelic frequencies for functional CYP2B6 polymorphisms appear comparable to those of Mexican ancestry (MEX) in the HapMap. For example, the 516T allele frequency in our cohort was 0.35 versus 0.29 in the MEX cohort (P > 0.05). However, the 516T frequency was statistically significantly higher (P < 0.05) than that of the Caucasian (CEU) HapMap cohort (0.27). It should, however, be noted that the HapMap MEX cohort is smaller (n = 100) than the Chilean cohort in the present study (and the HapMap CEU cohort).

This study identified 11 CYP2B6 SNPs associated with increased efavirenz exposure in Chilean patients. Based on LD across these 11 SNPs, 3 tagging SNPs were identified, which, when combined in a composite model, strengthened the observed association and had higher positive and negative predictive values for efavirenz concentrations above the proposed MTC (80.9% and 91.9%, respectively) than the c.516G→T TT genotype alone (71.1% and 90.2%, respectively) in this cohort. Given that individuals with a plasma concentration >4 µg/mL are three times more likely to develop CNS toxicity,3 a test with a high positive predictive value may be of clinical value. One of the three SNPs within the composite genetic model (rs10403955) tagged the c.516G→T SNP, which goes some way to explaining its contribution. However, it is less clear what the contribution of the other two SNPs within the composite model is. Both rs8192719 and rs2279345 are located in intronic enhancer sites, but have a low likelihood of directly exerting functional effects. Therefore, it is likely that their contribution is indirect, by tagging other functional variants within the CYP2B6 gene.

The frequency of the GCT inferred haplotype was higher in the Chilean patients than in the HapMap Caucasian population, but of lower frequency than in the HapMap Yoruba (Ibadan, Nigeria) population. It is therefore interesting to note that the plasma efavirenz concentrations in the present study (median = 2.6; range = 0.7–9.8 µg/mL) were ~23% higher than those reported previously in a Caucasian population (P = 0.008; median = 2.1; range = 0.7–17.4 µg/mL; n = 113) recruited in Germany.12 Curiously, however, concentrations were also higher than in a Black population (P = 0.04; median = 2.1; range = 0.5–26.0 µg/mL; n = 54) recruited as part of the same study,12 indicating that genetic factors other than CYP2B6 also contribute to the observed differences. Aside from CYP2B6, recent studies have identified the influence on efavirenz pharmacokinetics of polymorphisms within the CYP2A6 and UGT2B7 genes, which also warrant further investigation.7,15,16

One limitation of this study worthy of mention is the potential for the introduction of selection bias towards those patients more likely to need antiretroviral treatment and, crucially, those more likely to tolerate efavirenz. This may also have implications for the predictive values for MTC >4 µg/mL, since patients maintained on standard-dose efavirenz may not be representative of an antiretroviral-naive population. Unpublished data from C. P. Cortes from the national Chilean AIDS Cohort (ChiAC) from October 2001 to March 2008 identified 3107 patients starting an antiretroviral regimen with efavirenz. Of these, 639 (20.6%) had an adverse event requiring a change of drug, 3% of which were CNS related. The mean number of days to efavirenz discontinuation due to CNS issues for Fundación Arriarán, the major contributing centre of ChiAC, was 99 days (13–200 days) (C. P. Cortes, unpublished data).

In conclusion, these data indicate that a composite genetic model that includes multiple CYP2B6 SNPs is more strongly associated with efavirenz plasma concentrations than the c.516G→T polymorphism alone. These observations now warrant confirmation in other studies.

Funding

This work was supported by the NHS R&D Health Technology Assessment (HTA) Programme (project number 04/35/08) for the position of D. F. C. The National Institute of Health Research (NIHR–Department of Health) and the Northwest Development Agency (NWDA) provided infrastructural and project support.

Transparency declarations

None to declare.

References

1. Arribas JR, Pozniak AL, Gallant JE, et al. Tenofovir disoproxil fumarate, emtricitabine, and efavirenz compared with zidovudine/lamivudine and efavirenz in treatment-naive patients: 144-week analysis. J Acquir Immune Defic Syndr. 2008;47:74–8. [PubMed]
2. Gallant JE, DeJesus E, Arribas JR, et al. Tenofovir DF, emtricitabine, and efavirenz vs. zidovudine, lamivudine, and efavirenz for HIV. N Engl J Med. 2006;354:251–60. doi:10.1056/NEJMoa051871. [PubMed]
3. Marzolini C, Telenti A, Decosterd LA, et al. Efavirenz plasma levels can predict treatment failure and central nervous system side effects in HIV-1-infected patients. AIDS. 2001;15:71–5. doi:10.1097/00002030-200101050-00011. [PubMed]
4. Ward BA, Gorski JC, Jones DR, et al. The cytochrome P450 2B6 (CYP2B6) is the main catalyst of efavirenz primary and secondary metabolism: implication for HIV/AIDS therapy and utility of efavirenz as a substrate marker of CYP2B6 catalytic activity. J Pharmacol Exp Ther. 2003;306:287–300. doi:10.1124/jpet.103.049601. [PubMed]
5. Desta Z, Saussele T, Ward B, et al. Impact of CYP2B6 polymorphism on hepatic efavirenz metabolism in vitro. Pharmacogenomics. 2007;8:547–58. doi:10.2217/14622416.8.6.547. [PubMed]
6. Hofmann MH, Blievernicht JK, Klein K, et al. Aberrant splicing caused by single nucleotide polymorphism c.516G>T [Q172H], a marker of CYP2B6*6, is responsible for decreased expression and activity of CYP2B6 in liver. J Pharmacol Exp Ther. 2008;325:284–92. doi:10.1124/jpet.107.133306. [PubMed]
7. Kwara A, Lartey M, Sagoe KW, et al. CYP2B6 (c.516G→T) and CYP2A6 (*9B and/or *17) polymorphisms are independent predictors of efavirenz plasma concentrations in HIV-infected patients. Br J Clin Pharmacol. 2009;67:427–36. doi:10.1111/j.1365-2125.2009.03368.x. [PMC free article] [PubMed]
8. Ramachandran G, Hemanth Kumar AK, Rajasekaran S, et al. CYP2B6 G516T polymorphism but not rifampin coadministration influences steady-state pharmacokinetics of efavirenz in human immunodeficiency virus-infected patients in South India. Antimicrob Agents Chemother. 2009;53:863–8. doi:10.1128/AAC.00899-08. [PMC free article] [PubMed]
9. Rotger M, Tegude H, Colombo S, et al. Predictive value of known and novel alleles of CYP2B6 for efavirenz plasma concentrations in HIV-infected individuals. Clin Pharmacol Ther. 2007;81:557–66. doi:10.1038/sj.clpt.6100072. [PubMed]
10. Tsuchiya K, Gatanaga H, Tachikawa N, et al. Homozygous CYP2B6 *6 (Q172H and K262R) correlates with high plasma efavirenz concentrations in HIV-1 patients treated with standard efavirenz-containing regimens. Biochem Biophys Res Commun. 2004;319:1322–6. doi:10.1016/j.bbrc.2004.05.116. [PubMed]
11. Gatanaga H, Hayashida T, Tsuchiya K, et al. Successful efavirenz dose reduction in HIV type 1-infected individuals with cytochrome P450 2B6 *6 and *26. Clin Infect Dis. 2007;45:1230–7. doi:10.1086/522175. [PubMed]
12. Wyen C, Hendra H, Vogel M, et al. Impact of CYP2B6 983T>C polymorphism on non-nucleoside reverse transcriptase inhibitor plasma concentrations in HIV-infected patients. J Antimicrob Chemother. 2008;61:914–8. doi:10.1093/jac/dkn029. [PubMed]
13. Haas DW, Gebretsadik T, Mayo G, et al. Associations between CYP2B6 polymorphisms and pharmacokinetics after a single dose of nevirapine or efavirenz in African Americans. J Infect Dis. 2009;199:872–80. doi:10.1086/597125. [PMC free article] [PubMed]
14. Lorello G, la Porte C, Pilon R, et al. Discordance in HIV-1 viral loads and antiretroviral drug concentrations comparing semen and blood plasma. HIV Med. 2009;10:548–54. doi:10.1111/j.1468-1293.2009.00725.x. [PubMed]
15. Arab-Alameddine M, Di Iulio J, Buclin T, et al. Pharmacogenetics-based population pharmacokinetic analysis of efavirenz in HIV-1-infected individuals. Clin Pharmacol Ther. 2009;85:485–94. doi:10.1038/clpt.2008.271. [PubMed]
16. di Iulio J, Fayet A, Arab-Alameddine M, et al. In vivo analysis of efavirenz metabolism in individuals with impaired CYP2A6 function. Pharmacogenet Genomics. 2009;19:300–9. [PubMed]

Articles from Journal of Antimicrobial Chemotherapy are provided here courtesy of Oxford University Press