Abstract
Abstract
The use of pharmacogenomics (PGx) knowledge in treatment of individual patients is becoming a common phenomenon in the developed world. However, poorly resourced countries have thus far been constrained for three main reasons. First, the cost of whole genome sequencing is still considerably high in comparison to other (non-genomics) diagnostics in the developing world where both science and social dynamics create a dynamic and fragile healthcare ecosystem. Second, studies correlating genomic differences with drug pharmacokinetics and pharmacodynamics have not been consistent, and more importantly, often not indexed to impact on societal end-points, beyond clinical practice. Third, ethics regulatory frames over PGx testing require improvements based on nested accountability systems and in ways that address the user community needs. Thus, CYP2B6 is a crucial enzyme in the metabolism of antiretroviral drugs, efavirenz and nevirapine. More than 40 genetic variants have been reported, but only a few contribute to differences in plasma EFV and NVP concentrations. The most widely reported CYP2B6 variants affecting plasma drug levels include c.516G>T, c.983T>C, and to a lesser extent, g.15582C>T, which should be considered in future PGx tests. While the first two variants are easily characterized, the g.15582C>T detection has been performed primarily by sequencing, which is costly, labor intensive, and requires access to barely available expertise in the developing world. We report here on a simple, practical PCR-RFLP method with vast potentials for use in resource-constrained world regions to detect the g.15582C>T variation among South African and Cameroonian persons. The effects of CYP2B6 g.15582C>T on plasma EFV concentration were further evaluated among HIV/AIDS patients. We report no differences in the frequency of the g.15582T variant between the South African (0.08) and Cameroonian (0.06) groups, which are significantly lower than reported in Asians (0.39) and Caucasians (0.31). The g.15582C/T and T/T genotypes were associated with significantly reduced EFV levels (p=0.006). This article additionally presents the policy relevance of the PGX global health diagnostics and therefore, collectively makes an original interdisciplinary contribution to the field of integrative biology and personalized medicine in developing world. Such studies are, in fact, broadly important because resource-constrained regions exist not only in developing world but also in major geographical parts of the G20 nations and the developed countries.
Introduction
I
A second factor of importance for PGx testing in resource-constrained settings is that, although there are recognized examples of PGx applications, studies that correlate genomic differences with drug pharmacokinetics and pharmacodynamics have not been uniformly consistent, and more importantly, often not indexed to changes in public health or societal end-points, beyond clinical practice.
A third notable factor is the ethics frames over PGx testing. They, too, need to be improved but based on nested accountability systems whereby each innovation actor and stakeholder can crosscheck transparently each other's behavior, claims of novel PGx tests and their value for society (Dove and Ozdemir, 2014; Ostrom et al. 1999), and in ways that address the local community needs, bearing in mind the inter-dependency of “local” and “global” health. Antiretroviral drugs, we suggest, are one of the prime examples where the real-life value of genomics/genetics testing for personalized medicine can be examined in resource constrained world regions. Such studies are important since resource-constrained regions exist not only in developing world but also in major geographical parts of the G20 and the developed countries (Hotez, 2013).
Genetic variation in CYP2B6 has become important because of the increasing number of different substrates it metabolizes, as well as the elusive cure for HIV/AIDS, which means that more and more people will be put on antiretroviral drugs as part of chronic medication. The two main drugs that form the backbone of antiretroviral therapy (ART), efavirenz (EFV) and nevirapine (NVP), are substrates of CYP2B6 (di Iulio et al., 2009). CYP2B6 gene polymorphisms are particularly relevant in sub-Sahara Africa, where HIV/AIDS is continuing to affect the lives of millions of patients. More than 40 SNPs have been reported in the CYP2B6 gene and its regulatory regions. The most important are the variants that affect CYP2B6 expression and activity, which have a bearing on the systemic exposure to drugs that are principally metabolized by CYP2B6 (Arab-Alameddine et al., 2009; di Iulio et al., 2009).
Several studies have reported on the contribution of the different CYP2B6 variants, culminating in the observation of qualitative and quantitative differences in the distribution of CYP2B6 variants in different geographical populations (Dandara et al., 2011; Lamba et al., 2003; Swart et al., 2013). For example, the variant c.1459C>T, which is associated with decreased CYP2B6 activity (Lang et al., 2001) and occurs in 2%–3%, 9%–14% and 1% of Africans (observed in West Africans), Caucasians, and Chinese, respectively (Crettol et al., 2005; Guan et al., 2006; Jacob et al., 2004; Klein et al., 2005; Mehlotra et al., 2006). This example clearly shows that some variants are important in some population groups and not in others.
Recent studies have narrowed down the variants that contribute most of the observed variation in CYP2B6 activity to CYP2B6 c.516G>T, c.983T>C, and the intron 3 variant, g.15582C>T (Holzinger et al., 2012; Lamba et al., 2003). The first two are exonic mutations and affect protein structure, whereas the latter is an intronic mutation and has been associated with CYP2B6 functionality. The g.15582T variant disrupts a cryptic splice acceptor and also lies within an exonic splicing enhancer (ESE), thus, enhancing the skipping of exons 4–6 in individuals who also carry the exon 4 mutation, c.516G>T (g.15631G>T) (Lamba et al., 2003).
Since CYP2B6 is involved in the metabolism of two of the most extensively used ARV drugs, EFV and NVP, it is important that inexpensive and easy-to-use genetic tests are developed. Although the genotyping for CYP2B6 c.516G>T and c.983T>C has been presented using less expensive methods, evaluation of the intron 3 variation, g.15582C>T, has been carried out using sequencing, which is more costly and more challenging to analyze and interpret. The aim of this study was to design a simple, practical PCR-RFLP method with vast potentials for use in resource-constrained world regions to detect the g.15582C>T SNP, to report on the frequency of the CYP2B6 g.15582T-allele in South Africans and Cameroonians, and finally to evaluate the effects of the CYP2B6 g.15582C>T single nucleotide polymorphism on plasma EFV levels.
Materials and Methods
The effects of CYP2B6 g.15582C>T SNP were investigated in a cohort of 321 South Africans (104 healthy individuals and 217 HIV/AIDS patients) previously described by Swart et al. (2012; 2013) and the distribution of its frequency was investigated in these plus an additional 99 healthy Cameroonian individuals (Swart et al., 2012). The study received ethics clearance and approval from the Human Research Ethics Committee, Faculty of Health Sciences, University of Cape Town (HREC REF 103/2009).
DNA isolation was performed from whole blood (400 μL) or buffy coat (200 μL) samples, depending on availability, by using either the QIAmp Blood Midi Kit (QIAGEN, Hilden, Germany) or GENelute Kit (Sigma-Aldrich, St Louis, MO, USA). Primers flanking the CYP2B6 (NCBI accession number NC_000019.10) intron 3 variation were adapted from Lang et al. (2001) and synthesized by Integrated DNA Technologies (IDT). Primer sequences are CYP2B6-F: 5′-GGTCTGCCCATCTATAAAC-3′ and CYP2B6-4AS: 5′-TCCCTCTCCGTCTCCCTG-3′. The primer set was used to amplify a 411 bp PCR fragment encompassing two restriction sites for the Mbo I restriction enzyme. The PCR reaction mixture contained 50–100 ng of DNA, 1X Green GoTaq Reaction Buffer (Promega, Madison, USA), 0.4 mM dNTPs (BioLine, London, UK), 1.5 mM MgCl2, 0.4 μM of each primer, and 0.5 U GoTaq DNA polymerase. The “T100™ thermal cycler from Bio-Rad was used for PCR amplification. The PCR amplification conditions included an initial denaturation at 94°C for 3 min, followed by 30 cycles of denaturation at 94°C for 30 sec, annealing at 54.9°C for 30 sec, and extension at 72°C for 40s. This was followed by a final extension at 72°C for 10 min.
After PCR, the PCR product was digested overnight at 37°C using the restriction enzyme, Mbo I. The enzyme restriction digestion mix contained the following; 10 μL of PCR product, 1X Buffer R, and 3 U Mbo I (New England BioLabs, Inc., Ipswich, MA). PCR digests were electrophoresed on a 2.5% agarose gel for 120 min at 100V in the presence of GRGreen dye (InnoVita, Inc., Gaithersburg, MD). Digestion of the g.15582C-allele resulted in two fragments (92 bp and 319 bp), while in the presence of the g.15582T-allele, three fragments of 92 bp, 109 bp, and 210 bp were observed (Fig. 1). This PCR-RFLP method makes it easy for screening of the CYP2B6 g.15582C>T SNP and analyzing the possible association with particular phenotypes, for example, measures of response to EFV or NVP therapy.

Agarose gel electrophoresis (2.5% agarose gel, 80V, 2 h) for the detection of CYP2B6 g.15582C>T SNP. M=Molecular weight marker; Lanes 1, 4, 5, 7–9=homozygous C/C genotype (319 and 92 bp); Lane 10=T/T genotype (210 bp, 109 bp, 92 bp); Lanes 2 and 6=heterozygous C/T genotype (319 bp, 210 bp, 109 bp, and 92 bp).
Allele and genotype frequencies for different populations were obtained from the dbSNP database (http://www.ncbi.nlm.nih.gov/SNP/) and compared to the frequencies observed in this study among the South Africans. Pearson's χ2-test and Fisher's exact test were used to assess if differences in frequency distribution between populations were significant (p<0.05). Plasma EFV levels were log-transformed and ANOVA was used to determine the effect of CYP2B6 g.155282C>T SNP alone and also with or without a c.516T- or c.983C-allele. The genotyping results for CYP2B6 c.516G>T and c.983T>C SNPs in this cohort and the effects on EFV levels have previously been reported on by Swart et al. (2013).
Results and Discussion
The CYP2B6 g.15582T-allele was observed among 104 healthy South Africans and 99 Cameroonians, at a frequency of 8% and 6%, respectively. The allele and genotype frequencies were further compared to those of Caucasians, Asians, and other Africans using data from HapMap and 1000 genomes project (Table 1). No significant differences were observed between the allele and genotype frequencies when comparing South Africans to Cameroonians or Nigerians (i.e., Yoruba). Despite the similarity in allele and genotype frequencies for the CYP2B6 g.15582C>T SNP between South Africans, Cameroonians, Yoruba, and African-Americans, previous studies have shown that differences in minor allele frequencies for SNPs exist between these populations (Swart et al., 2012). African populations are, however, consistently shown to be significantly different from European and Asian populations (Dandara et al., 2014b) and this has major implications for personalised medicine in Low- and Middle-Income Countries (LMIC).
Indicates overall significant differences (p=0.0001) for genotype and allele frequencies compared to frequencies among South Africans.
Linkage disequilibrium analysis was performed between CYP2B6 g.15582C>T, c.516G>T, and c.983T>C. Linkage disequilibrium parameters (D' and r2, respectively) were as follows, g.15582C>T vs. c.516G>T (0.74 and 0.04), g.15582C>T vs c.983T>C (1.00 and 0.00), and c.516G>T vs. c.983T>C (0.95 and 0.01). Holzinger et al. (2012) reported that the CYP2B6 g.15582C>T and c.983T>C SNPs were in strong LD in a Caucasian population, thus, indicating that whereas in one population, testing for one SNP may give strong clues about the other SNP, in another population, there may be a need to test for each SNP independently.
For example, in a GWAS study that aimed to identify contributors to variable EFV levels, it was observed that after taking into account the c.516G>T and the c.983T>C variations, the g.15582C>T loci was significant (Holzinger et al., 2012), whereas in this study, among the 217 South African individuals with data on EFV levels, the CYP2B6 g.15582T-allele was associated with decreased plasma EFV levels (p=0.0059) (Fig. 2), but significance was lost upon adjusting for the c.516G>T and c.983T>C SNPs (Figs. 3 and 4) and this is in agreement with a recent report by Sinxadi et al. (2015) in a separate South African HIV/AIDS cohort.

Effect of the CYP2B6 g.15582C>T variation on plasma EFV levels. This analysis was done without stratifying for either CYP2B6 c.516G>T or c.983T>C single nucleotide polymorphisms.

Effect of CYP2B6 g.15582C>T on plasma EFV levels stratified for CYP2B6 c.516G>T and 983T>C SNPs. Only individuals with homozygous 516G/G and 983T/T genotypes were included in the analysis.

Effect of CYP2B6 g.15582C>T on plasma EFV levels among individuals with CYP2B6 c.516G/T or c.983T/C genotypes or their combination.
This is in contrast to the results reported by Holzinger et al. (2012) and Lamba et al. (2003), showing an association between the CYP2B6 g.15582T-allele alone with lower CYP2B6 expression and increased plasma EFV levels even after correcting for c.516G>T and the c.983T>C variations. Lamba et al. (2003) reported LD between the CYP2B6 g.15582C>T SNP and the CYP2B6 g.-2320T>C SNP in a putative HNF4 binding site within the promoter region, and the combination of g.-2320T/15582C was associated with high levels of CYP2B6 protein in Caucasian females, while the g.-2320C/15582T combination was associated with low levels of CYP2B6 protein (Lamba et al., 2003).
The effect of the CYP2B6 g.-2320T>C SNP should be further investigated, as this SNP could potentially be involved in the high levels of CYP2B6 that gives rise to the decreased EFV levels observed in this study. The c.516T- and c.983C-alleles have consistently been associated with CYP2B6 loss-of-function and significantly high plasma EFV concentrations (Thorn et al., 2010). However, for most of the other SNPs, associations with altered plasma EFV levels have been conflicting in different populations. For example, Kwara et al. (2009) reported an association between CYP2A6 -48T>G SNP and altered plasma EFV levels, but this association was not confirmed in two other studies in different populations (Elens et al., 2010; Leger et al., 2009). These conflicting reports make translation of PGx approach to clinical healthcare setting complicated, but certainly more ethnic specific.
The availability of this quick PCR-RFLP genotyping method will enable further genotyping of the CYP2B6 g.15582C>T SNP in other African populations and also in HIV/AIDS patients receiving either EFV or NVP to establish the importance of the CYP2B6 g.15582C>T SNP conclusively. This genotyping method can contribute to addressing the need for cost-effective PGx tests in resource limited settings, while new next generation sequencing (NGS) technologies still remains pricy and not easily available in developing countries.
We report on the effect of CYP2B6 g.15582C>T, alone and in combination with c.516G>T and c.983T>C on plasma EFV levels by using an easy method for quick screening of the CYP2B6 g.15582C>T SNP. We also provide data that the frequency of CYP2B6 g.15582C>T varies between different global population groups.
Implications for policy
The use of PGx for the treatment of individual patients in low and middle income countries has thus far been limited for two main reasons. First, genotyping has proven useful, but impractical and costly within resource constrained settings. This has been the case with other notable examples such as glucose-6-phosphate dehydrogenase (G6PD) deficiency (Relling et al., 2014). Second, studies that correlate genetic differences in absorption, distribution, metabolism, and excretion (ADME), with treatment efficacy in ART therapies, have generally been inconsistent.
Drug response phenotype is a complex phenotype influenced by environmental and genetic factors. In addition to gene–environment interactions, gene–gene interactions also play an important role in drug response. The variation in microbiome composition within and between individuals and how that variability affects drug response, the new area of pharmacomicrobiomics (ElRakaiby et al., 2014), also needs to be taken into account to understand the full spectrum of PGx. The interaction between these factors and their contribution to explaining differences in ADME is not fully understood, with only a portion of the variability in plasma EFV levels explained. Failure to take into account all variables likely to affect PGx could be the reason for inconsistencies in replication of PGx association studies in various studies or populations.
Although difficult to measure, a systems biology approach supported by bioinformatics and appropriate statistical analysis is the most likely to be able to incorporate all the PGx variables (ElRakaiby et al., 2014; Pinelli et al., 2012). The varied and different approaches used in PGx have led to generation of inadequate and inconsistent evidence to warrant adoption of a PGx approach (Pavlos et al., 2012; Pavlos and Phillips, 2012). This study provides not only evidence to support a PGx approach, but a methodology that can be adapted for a more cost effective and practical application in the clinical setting.
Policy that allows translation of this work is somehow limited. The greater proportion of PGx research is carried out by industry, and is targeted towards new drug candidates. Studies concerning existing drugs such as EFV and NVP have received less funding. Studies such as this one need to be supported by further research on translation and how findings may be implemented in clinical practice. Day to day PGx practice entails the collection, storage, and use of genetic data collected from patients. Managing this data requires specialist software and database systems for digital curation, processing, and collation of data, and qualified personnel with appropriate expertise.
The addition of genetic counsellors may also be a consideration as they have proved useful in supporting families affected by inherited genetic conditions. As the use of genomic technologies within routine care increases, many more patients and their families will need support to understand the implications of particular genetic findings, and perhaps to help them make important decisions regarding treatment.
Not only do primary health care providers and public health professionals need to be “genetically literate”, but also the public, if PGx is to be used in treatment. A consistent and coordinated approach to engaging with the public is needed to promote understanding of these new technologies and what it means for healthcare. Public engagement, education, or consultation exercises can be costly, particularly if the target population is large and widely distributed, or if there are several languages in use. High levels of literacy can contribute positively to public engagement efforts. Public engagement to introduce new health interventions, may however face an uphill struggle if there are high levels of population displacement and instability, as can happen in conflict zones.
Clear and stable regulatory and governance frameworks will assist in translation. A study by Hopkins et al. (2006) looking at Europe, the US, and Japan found that there is room for further harmonization and rationalization of regulatory frameworks governing clinical research in PGx. Furthermore, controls on the application of diagnostic tests remain a neglected area. For example, Hopkins et al. (2006) recommend assessment of the social and ethical impacts of PGx testing on a case-by-case basis. Earlier, De Vries et al. (De Vries, 2004; De Vries et al., 2004) have suggested to take into account both ethics and ethics of ethics in evaluation of new knowledge emergence in society. There has been insufficient regulatory analysis done to determine whether regulatory frameworks in the global south will help or hinder the development and commercialization of the type of test demonstrated in this article.
Also underexplored is how the assessment of social and ethical implications of testing is provided for by governance frameworks. More widely, the governance of PGx testing will need to cover aspects as far reaching as marketing and direct to consumer advertising, patents, access and intellectual property, liability, and exploitation (of patients, health care professionals, and governments from exaggerating the efficacy, scope, or necessity of tests), and addressing the controversy over the use of race and ethnicity in PGx research (World Health Organization, 2007).
This article therefore collectively makes an original contribution to the emerging field of “global OMICS” using the example of PGx testing in resource-constrained settings. Such studies with the aim of addressing the needs of global health diagnostics, could conceivably be extended to the testing of noncommunicable diseases (NCDs) that are now impacting the developed and developing world alike (Agirbasli et al., 2013).
Footnotes
Acknowledgments
We thank the South African Medical Research Council (SA MRC), the National Research Foundation (NRF) of South Africa, and the University of Cape Town for funding.
Author Disclosure Statement
No conflicting financial interests exist.
