Disease Markers
Volume 2015 (2015), Article ID 542543, 8 pages
http://dx.doi.org/10.1155/2015/542543

Research Article

The Affymetrix DMET Plus Platform Reveals Unique Distribution of ADME-Related Variants in Ethnic Arabs

Salma M. Wakil,¹ Cao Nguyen,² Nzioka P. Muiya,¹ Editha Andres,¹ Agnieszka Lykowska-Tarnowska,¹ Batoul Baz,¹ Asma I. Tahir,¹ Brian F. Meyer,¹ Grant Morahan,² and Nduna Dzimiri¹

¹King Faisal Specialist Hospital and Research Centre, Riyadh 11211, Saudi Arabia
²The University of Western Australia, Perth, WA, Australia

Received 7 October 2014; Revised 9 February 2015; Accepted 9 February 2015

Academic Editor: Helge Frieling

Copyright © 2015 Salma M. Wakil et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Background. The Affymetrix Drug Metabolizing Enzymes and Transporters (DMET) Plus Premier Pack has been designed to genotype 1936 gene variants thought to be essential for screening patients in personalized drug therapy. These variants include the cytochrome P450s (CYP450s), the key metabolizing enzymes, many other enzymes involved in phase I and phase II pharmacokinetic reactions, and signaling mediators associated with variability in clinical response to numerous drugs not only among individuals, but also between ethnic populations. Materials and Methods. We genotyped 600 Saudi individuals for 1936 variants on the DMET platform to evaluate their clinical potential in personalized medicine in ethnic Arabs. Results. Approximately 49% each of the 437 CYP450 variants, 56% of the 581 transporters, 56% of 419 transferases, 48% of the 104 dehydrogenases, and 58% of the remaining 390 variants were detected. Several variants, such as rs3740071, rs6193, rs258751, rs6199, rs11568421, and rs8187797, exhibited significantly either higher or lower minor allele frequencies (MAFs) than those in other ethnic groups. Discussion. The present study revealed some unique distribution trends for several variants in Arabs, which displayed partly inverse allelic prevalence compared to other ethnic populations. The results point therefore to the need to verify and ascertain the prevalence of a variant as a prerequisite for engaging it in clinical routine screening in personalized medicine in any given population.

1. Introduction

The response of an individual to drug therapy is determined by many variables, including food intake, age, and most importantly genetic factors. It is now well acknowledged that alterations in genes encoding drug metabolizing enzymes, ion transporters, and receptors, among others, have a great influence on the pharmacokinetics and pharmacodynamics of therapeutic agents, triggering variations in patient response to drug therapy. By far the largest of such gene families is that encoding the cytochrome P450 (CYP450) superfamily of enzymes, a diverse group of enzymes that catalyze the oxidation of organic substrates, including metabolic intermediates, such as lipids and steroidal hormones, and xenobiotics, such as drugs and other toxic substances [1–4]. The CYP450s are the major enzymes involved in about 75% of drug metabolism and bioactivation. Thereby, the CYP450 3A (CYP3A) isoforms are thought to be responsible for metabolizing approximately 55%, CYP2D6 about 30%, and CYP2C about 15–20% of the drugs [5].

Apart from the CYP450 metabolizers, other protein families involved in the absorption, distribution, metabolism, and elimination (ADME) of therapeutic agents include transporters, transferases, dehydrogenases, monooxygenases, reductases, and receptors. Transporter gene families include the adenosine triphosphate-binding cassette (ABC) and solute transporters (SLC) [6–8], while transferases comprise the glucuronosyltransferases (UGTs), glutathione S transferases (GSTs), sulfotransferases (SULTs), and catechol-O-methyl transferase (COMT) that may also determine the fate of a drug [9–11]. Other ADME-related gene families include those encoding dehydrogenases (DHO), monooxygenases, and reductases [12, 13]. Besides, changes in some receptor families, such as the nuclear receptor subfamilies, β-adrenergic receptors, and peroxisome proliferator-activated receptor subtypes, just to name a couple, also contribute to variations in the ways by which patients respond to drug therapy [14–17].

It is now well established that patient response to drug therapy will depend on the genetic structure of the enzyme or protein involved in its metabolism, and the allelic frequencies and phenotypic consequences may vary considerably between ethnic groups [18–22]. Hence, great effort has been directed at exploiting our knowledge of changes in these genes for clinical purposes in personalized medicine. One of these endeavors is the development of the Affymetrix DMET Plus array which facilitates the highly multiplexed genotyping of known polymorphisms of a panel of markers from 225 ADME-related genes on a single array. These variants have been documented to be of importance in phase I and phase II drug metabolism and disease manifestation. While this platform has found its place in personalized medicine, it has also become increasingly apparent that the influence of these genetic changes varies not only among individuals, but also between ethnical populations. This renders it necessary to acquire adequate prior knowledge of the likely impact of such entities, if they were to be considered for routine clinical purposes in any given society. Currently only isolated data is available on the prevalence of the gene variants in ethnic Arabs. Hence, this study was designed to characterize the prevalence of these variants with the focus on establishing their relevance in personalized management of various indications with different types of drugs in ethnic Arabs, employing the Saudi Arab population as a homogeneous study model.

2. Methods

2.1. Patient Sample Collection

The study candidates comprised 600 individuals randomly drawn from our coronary artery disease (CAD) registry who were subjected to genotyping by the DMET Plus chip. This registry contains affected as well as nonaffected individuals. All individuals signed an informed consent, and the study protocol was approved by the Institutional Review Board (IRB) at the King Faisal Specialist Hospital and Research Centre. Genomic DNA was extracted by a standard phenol extraction procedure.

2.2. Genotyping by Affymetrix DMET Plus Array

The genotyping was accomplished by the DMET (Drug Metabolizing Enzymes and Transporters) Plus Premier Pack, a microarray assay developed by Affymetrix (Affymetrix, Santa Clara, CA, USA) designed specifically to test drug metabolism associations. The DMET array contains 1936 (1931 SNPs and 5 CNVs) drug metabolism markers in 225 genes including 47 phase I enzymes, 80 phase II enzymes, 52 transporters, and 46 other genes. These genetic variants were multiplex genotyped using the molecular inversion probe (MIP) technology [23, 24]. Briefly, some markers from regions containing pseudogenes and close homologs were first preamplified using a multiplex polymerase chain reaction (mPCR) (Qiagen, Valencia, CA, USA). Genomic sequences that contain the polymorphic markers of interest were then preferentially amplified through the use of highly selective MIPs. A first quality control (QC) gel was run to determine the quality of amplified MIPs, which should be a single band represented on a gel in the range of 100–150 base pairs. Smaller DNA fragments were generated by adding fragmentation reagents to improve sample hybridization with the DMET Plus array, and DNA fragment size was checked on the second QC gel, in which the fragment length should be less than 120 base pairs with a smear centered at approximately 50 base pairs. The resulting target DNA was then labelled and hybridized to the DMET Plus array to obtain genotypes using a single color detection format. The profiles for the genotyping call rates and concordance comparisons were generated by the DMET Console software which is based on the BRLMM (Bayesian Robust Linear Model with Mahalanobis distance classifier) algorithm. Fixed genotype boundaries were used as the algorithms for all genotyping configurations. Genotypes were reported as homozygous wild type, heterozygous, homozygous variant or “no call.” CNV markers and SNPS with call rate less than 100% were excluded from the subsequent analysis.

3. Results

The present study genotyped 600 Saudi individuals for the 1936 gene variants involved in drug absorption, distribution, metabolism, and elimination (ADME), with documented functional significance in phase I and II drug metabolism, as well as pharmacodynamics of several therapeutic agents. The gene variants, which include transporters, ion transferases, and receptors, displayed various distribution profiles. Figure 1 summarizes the profiles of the different gene variants by functional groups. As indicated in this figure, overall some 877 (45%) displayed no allelic change at all in our study population.

Figure 1: Comparison of the detected versus non-detected variants by the functional groups of the gene families. ATPT, ATP transporters; OTT, other transporters; TFR, transferases; DHO, dehydrogenases; MOX, monooxygenases; RDTS, reductases; RCPT, receptors; MSC, miscellaneous gene families.

By far the largest number of SNPs on the DMET platform belongs to the drug metabolizing superfamilies, whereby the CYPs constitute the majority. This superfamily comprised about 437 SNPs, of which the CYP2C (53), CYP3A (53), CYP1A (30), and CYP2D (30) constitute the major subfamilies (Table 1). When classified by subfamilies, the data revealed that about 51% of the CYP2C, 32% of the CYP3A, 33% of CYP1A, and 53% of CYP2D variants (Table 1) displayed detectable minor alleles. Furthermore, while, in the majority of the cases, the minor allele frequencies (MAFs) fell in the range of 0.001–0.5, noticeable variations were observed among the family members. To begin with, as depicted in Figure 2, several SNPs exhibited an inverse distributional profile compared to available databases on other populations, such as the Caucasians or Chinese (see DMET Supplementary Data in Supplementary Material available online at http://dx.doi.org/10.1155/2015/542543). Thus, our population displays a whole profile range from variants, such as rs2072200_C>G or rs1573496_C>G in which the minor allele in other ethnic populations turned to be the major allele or vice versa in our population, to those lacking any genetic change such as rs3740071_C>G, rs8187797_C>G, or rs11568421_G>A in our population as opposed to others (Figure 2). Besides, several of the variants displayed either marginally or significantly greater MAFs than those in other ethnical populations, with sizeable number exhibiting MAFs >0.49 (DMET Supplementary Data). Notably, while several variants found in other populations in some smaller CYP gene subfamilies, such as CYP11A, CYP21A, and CYP46A, were not at all detected, those in the CYP2F, CYP4B, CYP7B, CYP8B, and CYP26A subsets displayed discernible minor alleles (Table 1; DMET Supplementary Data).

Table 1: Distribution of some of the common CYP450 variants.

Group
Total
Undetected
<0.01
>0.01

All
437
222 (0.51)
55 (0.125)
160 (0.366)

CYP1A
30
20 (0.666)
1 (0.033)
9 (0.3)

CYP1B
20
16 (0.80)
1 (0.5)
3 (0.15)

CYP2A
34
20 (0.58)
3 (0.09)
11 (0.32)

CYP2B
24
10 (0.42)
2 (0.08)
12 (0.5)

CYP2C
53
26 (0.49)
14 (0.26)
13 (0.25)

CYP2D
30
14 (0.47)
4 (0.13)
12 (0.4)

CYP2E
10
2 (0.20)
3 (0.3)
5 (0.5)

CYP2F
5
0 (0.0)
2 (0.4)
3 (0.6)

CYP2J
9
8 (0.89)
0 (0.0)
1 (0.11)

CYP2S
4
1 (0.25)
1 (0.25)
2 (0.5)

CYP3A
53
36 (0.68)
3 (0.06)
14 (0.26)

CYP4A
4
1 (0.25)
1 (0.25)
2 (0.5)

CYP4B
7
0 (0.0)
2 (0.28)
5 (0.72)

CYP4F
35
12 (0.34)
2 (0.06)
21 (0.6)

CYP4Z
4
1 (0.25)
0 (0.0)
3 (0.75)

CYP7A
7
1 (0.142)
1 (0.142)
5 (0.714)

CYP7B
4
0 (0.0)
1 (0.25)
3 (0.75)

CYP8B
7
0 (0.0)
4 (0.571)
3 (0.428)

CYP11A
7
7 (1.00)
0 (0.0)
0 (0.0)

CYP11B
24
15 (0.63)
0 (0.0)
9 (0.36)

CYP17A
6
3 (0.50)
0 (0.0)
3 (0.50)

CYP19A
11
6 (0.54)
0 (0.0)
5 (0.46)

CYP20A
5
2 (0.40)
0 (0.0)
3 (0.60)

CYP21A
3
3 (1.00)
0 (0.0)
0 (0.0)

CYP24A
8
4 (0.50)
0 (0.0)
4 (0.50)

CYP26A
5
0 (0.0)
2 (0.4)
3 (0.6)

CYP26C
1
0 (0.0)
1 (1.0)
0 (0.0)

CYP27A
6
4 (0.67)
2 (0.33)
0 (0.0)

CYP27B
9
5 (0.56)
2 (0.22)
2 (0.22)

CYP39A
4
1 (0.25)
1 (1.25)
2 (0.50)

CYP46A
1
1 (1.00)
0 (0.0)
0 (0.0)

CYP51A
7
3 (0.43)
2 (0.285)
2 (0.285)

The table displays the relative minor allele frequency distribution of the CYP450 variants determined by DMET Plus chip in 600 individuals.

Figure 2: Examples of SNPs displaying an inverse relationship of the allele frequencies in the minor alleles in Caucasians (CEU) and Chinese (CHN) and ethnic Saudi Arab populations (KSA).

In addition to the metabolizing CYPs, the other large groups of ADME-related variants included the transporters such as the SLCs comprising 322 variants, of which 58% were detectable and the ABCs comprising 242 variants of which 56% were detected in our population. The partial lack of change was also evident among other superfamilies of transporters, transferases, dehydrogenases, monooxygenases, reductases, receptors, and other signalling entities (Figure 1; Table 2). In summary, we were also unable to detect approximately 41% of the variants in other major ADME gene families, including the ABCs, SLCs, SULTs, GSTs, and PPARs.

Table 2: Distribution levels for non-CYP450 variants on the DMET Plus platform.

Group
Type
All
MAF distribution

Undetected
<0.01
>0.01

Transporters
ABC transporters
242
108 (0.44)
33 (0.14)
101 (0.42)

Solute transporters
322
137 (0.42)
47 (0.15)
138 (0.43)

ATPase, Cu²⁺ transporters, alpha peptide
17
7 (0.41)
1 (0.06)
9 (0.53)

Transferases
Glucuronosyltransferases (UGTs)
116
58 (0.50)
11 (0.10)
47 (0.40)

Sulfotransferases
138
45 (0.33)
18 (0.13)
75 (0.54)

Glutathione S-transferase
88
42 (0.48)
4 (0.04)
42 (0.48)

N-Acetyl glucosamine transferase
37
20 (0.54)
5 (0.13)
12 (0.32)

Others (COMT, HNMT, MAT, NNMT, TPMT, QPRT, PNMT)
40
19 (0.48)
5 (0.12)
16 (0.40)

Dehydrogenases
ADHs; ALDHs, XDHs, G6PD
104
54 (0.52)
7 (0.07)
43 (0.41)

Monooxygenases
FMOs, MAOs, AOX
92
38 (0.41)
10 (0.11)
44 (0.48)

Reductases
VKORC1, HMGCRs, POR, CBR
43
16 (0.37)
7 (0.16)
20 (0.47)

Receptors
NR1, members 2 and 3
21
11 (0.52)
3 (0.14)
7 (0.33)

PPAR, gamma, delta
53
13 (0.25)
2 (0.03)
38 (0.72)

Nuclear receptor subfamily 3
12
5 (0.42)
5 (0.42)
2 (0.16)

Others (RALBP1, SPN, APOA, RXRA, AHR, ARNT)
25
11 (0.44)
4 (0.16)
10 (0.40)

Miscellaneous
ALB, ARSA, CA5P AKAP, CDA, CES, PTGIS, CCDC, CROT, SPG, FAAH, EPHX, DCK, PON, TBXAS, TPSG, ORM, PGAP3, SERPINA, PPP1R9, PRSS53, RPL, TYMS
144
71 (0.49)
21 (0.15)
52 (0.36)

ABC, ATP-binding cassette transporters; ALDH, aldehyde dehydrogenase; ADH, alcohol dehydrogenase class 3; ALB, albumin; ARSA, arylsulfatase; AHR, aryl hydrocarbon receptor; AKAP, A-kinase anchor proteins; APOA, apolipoprotein receptor, subtype A; AOX, alternative oxidase; CA5P, carbonic anhydrase; CBR, carbonyl reductase; CCDC, coiled-coil domain-containing protein; CDA, cytidine deaminase; CES, carboxylesterase; COMT, catechol-O-methyl transferase; CROT, peroxisomal carnitine O-octanoyltransferase; DCK, deoxycytidine kinase; EPHX, microsomal epoxide hydrolase; FAAH, fatty acid amide hydrolase; FMO, flavin-containing monooxygenase; G6PD, glucose 6-phosphate dehydrogenase; HNMT, histamine N-methyltransferase; MAO, monoamine oxidase; MAT, methionine adenosyl transferases; NNMT, nicotinamide N-methyltransferase; NR1, nucleic acid receptor, group 1; ORM, orosomucoid; PGAP3, post-GPI attachment to proteins 3; PON, serum paraoxonase/arylesterase; POR, cytochrome P450 reductase; PPAR, peroxisome proliferator-activated receptors; PNMT, phenylethanolamine N-methyltransferase; PTGIS, prostaglandin-I synthase; QPRT, quinolinate phosphoribosyltransferase; PPP1R9, protein phosphatase 1, regulatory subunit 9B; PRSS53; serine protease, 53; RALBP1, ralA binding protein 1; RPL, ribosomal protein-binding receptor; RXRA, retinoic acid receptor; SERPINA, serpin peptidase inhibitor, clade A (alpha-1 antiproteinase, antitrypsin); SPG, spastic paraplegia; SPN, sialophorin receptor; TBXAS, thromboxane A synthase 1 (platelet); TYMS, thymidylate synthetase; TPMT, thiopurine S-methyltransferase; TPSG, tryptase; XDH, xanthine dehydrogenases.

Apart from the transporters, the platform also carries several families of transferases, including the uridine diphosphate glucuronosyltransferases (UGTs), SULTs, glutathione S-transferase alpha (GSTA), glutathione S-transferase Mu (GSTM), glutathione S-transferase omega (GSTO), glutathione S-transferase pi (GSTP), histamine N-methyltransferase, methionine adenosyl transferase (MAT), N-acetyl glucosamine transferase, nicotinamide N-methyltransferase, thiopurine S-methyltransferase (NNMT), COMT, phenylethanolamine N-methyltransferase (PNMT), and quinolinate phosphoribosyltransferase (QPRT). The largest transferase gene superfamily studied was that of the microsomal UGTs, composed of at least 16 genes responsible for the elimination of a myriad of xenobiotics and endogenous compounds. Of the 116 UGT variants, only 50% showed detectable changes (Table 2). Furthermore, we were also able to detect 67% of the 138 SULTs, 52% of the 88 GSTs, and majority of the others.

Aldehyde dehydrogenases (ALDH; 50) and alcohol dehydrogenases (ADHs; 40) similarly constituted the greater part of the 111 DHO family. Other DHOs included six glucose 6-phosphate dehydrogenase (G6PD), seven nicotinamide phosphate (quinone) dehydrogenases (NQOs), and eight xanthine dehydrogenases (XDHs). Altogether 54 (52%) of the DHOs did not show any changes. Furthermore, monooxygenases, such as FMOs or MAOs, constitute a group of 92 variants, of which about 54% were detectable (Table 2). Besides, about 64% of the studied receptors and other signaling gene variants were detectable, including 21 for ralA binding protein 1 (RALBP1), 12 for nuclear receptor subfamily 1, group I, members 2 and 3 (NR1l2 and NR1l3), 53 for PPARs, and 12 for nuclear receptor subfamily 3, group C, member 1 (NR3C), retinoid X receptor, alpha (RXRA) nuclear receptor subfamily 2, group B member, to name a few. Also, approximately 61% of these SNPs showed detectable minor alleles in the present study (Table 2).

Put together, we were able to detect minor alleles for 49% of the CYPs, 57% of the transporters, and 56% of the transferases on the DMET platform in our study population. Of the rest of the gene variants, we were unable to detect minor alleles for 43%. Among these SNPs, 72 showed MAF values of >0.45, whereby 18 displayed values of >0.49.

3.1. Statistical Analysis

Comparison of genotypes and alleles between different groups for continuous dependent variables was accomplished by analysis of variance (ANOVA) or Student’s t-test as appropriate. Categorical variables were analyzed by Chi-Square test, and logistic regression analysis was used to compute odds ratios and their 95% confidence intervals. All other statistical analyses were performed using the SPSS software version 14 (SPSS Inc., Chicago, USA). Associations with a two-tailed value < 0.05 were considered statistically significant.

4. Discussion

The present study established the prevalence of the 1936 DMET Plus platform variants in several gene families involved in the pharmacokinetics and pharmacodynamics of several important therapeutic agents for different ailments. We detected approximately 55% of the SNPs on this platform, pointing to the fact that only a portion of them is likely to be economically worthwhile pursuing in personalized medicine in ethnic Arabs. Currently, there is great lack of data on the distribution of these ADME variants in this population. In fact, to the best of our knowledge, this is the first and largest study reporting their prevalence in an Arab population. This data should therefore serve as a basis for evaluating the usefulness of routinely assaying these SNPs for clinical purposes in this ethnic group.

Perhaps the most widely studied family of metabolizing enzymes is the CYP superfamily. These enzymes display a wide range of phenotypes from poor, rapid, to ultrarapid metabolizers for several important agents, due to the variations in the combinations of their encoding alleles. An example is that of the CYP2C19 with more than 19 variants encoding the nonfunctional CYP2C12 and inactive enzyme CYP2C13, on one hand, and an ultrarapid metabolizing CYP2C117 and extensive metabolizer phenotype CYP2C19 1 on the other hand [4, 25, 26]. Accordingly, poor metabolizers of drugs that are processed through the CYP2C19 pathways frequently experience dramatic changes in drug responses and side effects when they receive standard doses. Thus, for example, the CYP2C12 loss-of-function allele has been associated with a decreased activation of clopidogrel [27, 28], attenuation of its antiplatelet effect [29–35], and contributing to 3- to 6-fold incidence of stent thrombosis in patients treated with percutaneous coronary intervention (PCI) [29–35], while the presence of any gain-of-function CYP2C19 17 has also been linked to increased risk of bleeding [36]. Given these potential clinical consequences of harbouring the CYP2C19 gene variants that may affect therapeutic modalities, it is not surprising that a surge of attempts has grown exponentially in recent years to employ this knowledge clinically in personalized medicine. Besides, some researchers have suggested a link between CYP2C19 polymorphisms and diseases, such as digestive tract cancer [37] and essential hypertension [38]. However, their phenotypic expression has been studied primarily in Caucasians [20] and some other ethnic populations, but only poorly so in Arabs. In fact, very limited information is currently available on the prevalence of these variants in the Saudi population, with only a couple of studies appearing recently on two variants, the CYP2C12 and CYP2C13, albeit involving very small study populations [39, 40]. Hence, the establishment of their prevalence in the present study can be viewed as an important step in identifying the clinically relevant SNPs in this population. Specifically, our results indicate that it is worthwhile screening for the different CYP2C19 variants, for example, for such purposes in our population.

Like the CYP450s, several gene variants encoding other ADME-related proteins also exhibited diversity in their distribution, ranging from those that showed no changes, such as the rs6193_A>G and rs258751_G>A, to those that exhibited inverse profiles, such as rs3740071_C>G and rs17216887_C>G in the ABCs. To date, mutations in the ABCP have been associated with cancer chemotherapy drug resistance [41, 42], atherosclerosis, inflammation [43], and several other diseases [43–45], while disorders linked to ATP7A include Menkes disease and occipital horn syndrome [46, 47]. Hence, in our population, these variants may be relevant not only with respect to drug response, but also disease manifestation, and further studies are necessary to elucidate the extent of their impact on disease in this ethnic population.

Although much remains to be learnt about UGTs, a number of polymorphisms are thought to be of toxicological significance [48] or have been associated with diseases, such as Crigler–Najjar’s and Gilbert’s syndrome [49]. Thus, several of its SNPs, including the rs7586110 (UGT1A12_c.-57T>G) (MAF = 0.417) and rs8175347 (UGT1A28_c.TATA-box) (MAF = 0.271), have been previously linked to different types of diseases, including cancer, cardiovascular diseases, and irinotecan toxicity in patients with Gilbert’s syndrome, to name a few [50–53]. Perhaps one of the best-studied transferase gene families is the GST, which also constitutes one of the largest groups of variants on the platform. The prevalence of the wide majority of the GSTs was similar to that described in other ethnic groups, suggesting that their impact on disease is likely to be global. Furthermore, the study also revealed significant changes in the studied SULTs, which constitute the third largest family of transferases on the platform. Since SULTs can activate procarcinogens to reactive electrophiles [11], enzymes such as the steroid sulfatase and estrogen SULTs have been implicated in human carcinomas [54].

In addition to transporters and transferases, other ADME families including DHOs, monooxygenases, reductases, various receptors, and signal transducers on the platform also equally displayed diversity in the variant profiles. Put together, the data indicates that approximately 45% variants could not be detected in our study population. Furthermore, the majority of those that were detectable presented with MAF >0.01, with a sizeable portion being at variance with the data in the literature.

Since almost 50% of the loci were unchanged in the present population, it was of interest for us to compare the discovered profiles with those of other ethnic groups, as a test for the robustness of the DMET platform as a potential global clinical tool. As might be expected, our analysis revealed some similarities with other ethnic populations in the alterations of a number of variants. Thus, for example, the results point to relatively similar frequencies for the CYP2C12 (rs4244285) to those in several other ethnic groups, including the Romanian (0.12), Lebanese (0.13) [40], Turkish (0.12) [55–57], Jewish (0.15%) [39], Russian (0.11) [58], and Italian (0.12) [59], but slightly lower than those in the Chinese (0.25) [60, 61], North Indian (0.26) [62], and Thai (0.29) [22] populations. We also found low MAFs for CYP2C12 (0.093) and CYP2C13 (0.001) which were comparatively lower than in Africans, while that of the CYP2C19 17 (0.256) matched those of the European populations but was higher than those in African and other Asian populations (see also Supplementary Data). More importantly, in depth analysis pointed to several SNPs, featuring conversions in which a minor allele in the European/Caucasian populations not only turned to be the major one but also exhibited no change in the Saudi population. On the other end of the spectrum were also several major alleles in European or Asian populations that could not be detected in our population. Thus, put together, the study demonstrates that although the distribution of most of these variants was within similar ranges with those in other populations, some distinct interethnical differences in the prevalence of many others were also evident between ethnic Arabs and other ethnic populations.

The important question arises as to the clinical relevance of these findings with respect to targeting the variants for personalized medicine. First, our observations stress the fact that not all ADME variants constitute therapeutically meaningful targets in the ethnic Arab population, as reflected by their absence in our study population. The wide interethnical variations in the prevalence of several of the variants supports the notion that the depth of involvement of these variants will also vary among different ethnic groups. Thus, for example, several genotypes that are otherwise lowly distributed in other ethnic population might be of great significance and vice versa, in this regard. In particular, the observations of inverse genotype relationships in which the alleles displayed literally the opposite level of expression, such as rs3740071 or rs4699735, might also imply, for example, that these variants will exert opposing effects on drug response in different ethnic groups. This, in turn, renders it practically impossible to generalize the mode by which such variants may influence therapeutic modalities globally and therefore necessitates acquisition of adequate knowledge of their prevalence in any given community prior to engaging them in targeted genotyping for clinical purposes in personalized medicine. Besides, our current findings further open the door to also critically evaluate the role of the studied gene variants in disease. Hence, their actual clinical impact on disease management needs to be revisited more closely.

In summary, the present study utilized the availability of the DMET Plus platform to estimate the prevalence of ADME-related variants of potential therapeutic relevance, using the Saudi population, as a basis for informed targeting of these variants in personalized medicine in ethnic Arabs. We were able to detect approximately half of the variants on this platform, not only reaffirming the prevalence of some important variants in our population, but also furnishing some support for the usefulness of the procedure in routinely detecting the presence of these genotypes for clinical purposes. More importantly, we observed some significant differences in the expression of several variants in comparison to other ethnic populations, laying the foundation for adopting evidence-based approaches to personalized medicine in ethnic Arabs.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was supported through a Cardiovascular Research Grant 2010020 of the King Faisal Specialist Hospital and Research Centre. The authors wish to express their gratitude for the support.

References

H.-G. Xie, R. B. Kim, A. J. J. Wood, and C. M. Stein, “Molecular basis of ethnic differences in drug disposition and response,” Annual Review of Pharmacology and Toxicology, vol. 41, pp. 815–850, 2001.

O. Tribut, Y. Lessard, J.-M. Reymann, H. Allain, and D. Bentué-Ferrer, “Pharmacogenomics,” Medical Science Monitor, vol. 8, no. 7, pp. RA152–RA163, 2002.

M. Kazui, Y. Nishiya, T. Ishizuka et al., “Identification of the human cytochrome P450 enzymes involved in the two oxidative steps in the bioactivation of clopidogrel to its pharmacologically active metabolite,” Drug Metabolism and Disposition, vol. 38, no. 1, pp. 92–99, 2010.

J. L. Fayer Rehmel, J. A. Eckstein, N. A. Farid et al., “Interactions of two major metabolites of prasugrel, a thienopyridine antiplatelet agent, with the cytochromes P450,” Drug Metabolism and Disposition, vol. 34, no. 4, pp. 600–607, 2006.

J. G. Hardman and L. E. Limbird, Goodman and Gilman's Pharmacological Basis of Therapeuits, The McGraw-Hill, New York, NY, USA, 9th edition, 1996.

A. L. Davidson, E. Dassa, C. Orelle, and J. Chen, “Structure, function, and evolution of bacterial ATP-binding cassette systems,” Microbiology and Molecular Biology Reviews, vol. 72, no. 2, pp. 317–364, 2008.

P. M. Jones and A. M. George, “The ABC transporter structure and mechanism: perspectives on recent research,” Cellular and Molecular Life Sciences, vol. 61, no. 6, pp. 682–699, 2004.

C. Cobbold, J. Coventry, S. Ponnambalam, and A. P. Monaco, “The Menkes disease ATPase (ATP7A) is internalized via a Rac1-regulated, clathrin- and caveolae-independent pathway,” Human Molecular Genetics, vol. 12, no. 13, pp. 1523–1533, 2003.

F. W. Smith, M. J. Hawkesford, I. M. Prosser, and D. T. Clarkson, “Isolation of a cDNA from Saccharomyces cerevisiae that encodes a high affinity sulphate transporter at the plasma membrane,” Molecular and General Genetics, vol. 247, no. 6, pp. 709–715, 1995.

M. H. Grossman, B. S. Emanuel, and M. L. Budarf, “Chromosomal mapping of the human catechol-O-methyltransferase gene to 22q11.1→q11.2,” Genomics, vol. 12, no. 4, pp. 822–825, 1992.

N. Gamage, A. Barnett, N. Hempel et al., “Human sulfotransferases and their role in chemical metabolism,” Toxicological Sciences, vol. 90, no. 1, pp. 5–22, 2006.

S. Harayama, M. Kok, and E. L. Neidle, “Functional and evolutionary relationships among diverse oxygenases,” Annual Review of Microbiology, vol. 46, pp. 565–601, 1992.

M. H. M. Eppink, C. Bunthof, H. A. Schreuder, and W. J. H. van Berkel, “Phe¹⁶¹ and Arg¹⁶⁶ variants of p-hydroxybenzoate hydroxylase: implications for NADPH recognition and structural stability,” FEBS Letters, vol. 443, no. 3, pp. 251–255, 1999.

D. J. Mangelsdorf, C. Thummel, M. Beato et al., “The nuclear receptor super-family: the second decade,” Cell, vol. 83, no. 6, pp. 835–839, 1995.

V. Chandra, P. Huang, Y. Hamuro et al., “Structure of the intact PPAR-γ-RXR-α nuclear receptor complex on DNA,” Nature, vol. 456, no. 7220, pp. 350–356, 2008.

A. Elbrecht, Y. Chen, C. A. Cullinan et al., “Molecular cloning, expression and characterization of human peroxisome proliferator activated receptors γ1 and γ2,” Biochemical and Biophysical Research Communications, vol. 224, no. 2, pp. 431–437, 1996.

J. R. Jones, C. Barrick, K.-A. Kim et al., “Deletion of PPARγ in adipose tissues of mice protects against high fat diet-induced obesity and insulin resistance,” Proceedings of the National Academy of Sciences of the United States of America, vol. 102, no. 17, pp. 6207–6212, 2005.

J. McGraw and D. Waller, “Cytochrome P450 variations in different ethnic populations,” Expert Opinion on Drug Metabolism and Toxicology, vol. 8, no. 3, pp. 371–382, 2012.

C. Jaja, W. Burke, K. Thummel, K. Edwards, and D. L. Veenstra, “Cytochrome P450 enzyme polymorphism frequency in indigenous and native American populations: a systematic review,” Community Genetics, vol. 11, no. 3, pp. 141–149, 2008.

K. Kurose, E. Sugiyama, and Y. Saito, “Population differences in major functional polymorphisms of pharmacokinetics/pharmacodynamics-related genes in Eastern Asians and Europeans: implications in the clinical trials for novel drug development,” Drug Metabolism and Pharmacokinetics, vol. 27, no. 1, pp. 9–54, 2012.

K. Jazdzewski, J. Boguslawska, J. Jendrzejewski et al., “Thyroid hormone receptor beta (THRB) is a major target gene for microRNAs deregulated in papillary thyroid carcinoma (PTC),” The Journal of Clinical Endocrinology and Metabolism, vol. 96, no. 3, pp. E546–E553, 2011.

R. Polimanti, S. Piacentini, D. Manfellotto, and M. Fuciarelli, “Human genetic variation of CYP450 superfamily: analysis of functional diversity in worldwide populations,” Pharmacogenomics, vol. 13, no. 16, pp. 1951–1960, 2012.

P. Hardenbol, J. Banér, M. Jain et al., “Multiplexed genotyping with sequence-tagged molecular inversion probes,” Nature Biotechnology, vol. 21, no. 6, pp. 673–678, 2003.

L. Dunbar, W. Miles, A. Wheeler, J. Sheridan, J. Pulford, and R. Butler, “The CYP2D6 metaboliser status of patients prescribed risperidone for the treatment of psychosis,” New Zealand Medical Journal, vol. 122, no. 1296, pp. 29–34, 2009.

R. J. Ferguson, S. M. F. De Morais, S. Benhamou et al., “A new genetic defect in human CYP2C19: mutation of the initiation codon is responsible for poor metabolism of S-mephenytoin,” Journal of Pharmacology and Experimental Therapeutics, vol. 284, no. 1, pp. 356–361, 1998.

M. Ingelman-Sundberg, S. C. Sim, A. Gomez, and C. Rodriguez-Antona, “Influence of cytochrome P450 polymorphisms on drug therapies: pharmacogenetic, pharmacoepigenetic and clinical aspects,” Pharmacology and Therapeutics, vol. 116, no. 3, pp. 496–526, 2007.

G. C. Ibeanu, J. Blaisdell, B. I. Ghanayem et al., “An additional defective allele, CYP2C19*5, contributes to the S-mephenytoin poor metabolizer phenotype in Caucasians,” Pharmacogenetics, vol. 8, no. 2, pp. 129–135, 1998.

J. G. Bramness, S. Skurtveit, L. Fauske et al., “Association between blood carisoprodol:meprobamate concentration ratios and CYP2C19 genotype in carisoprodol-drugged drivers: decreased metabolic capacity in heterozygous CYP2C19∗1/CYP2C19∗2 subjects?” Pharmacogenetics, vol. 13, no. 7, pp. 383–388, 2003.

K. A. Kim, P. W. Park, S. J. Hong, and J.-Y. Park, “The effect of CYP2C19 polymorphism on the pharmacokinetics and pharmacodynamics of clopidogrel: a possible mechanism for clopidogrel resistance,” Clinical Pharmacology and Therapeutics, vol. 84, no. 2, pp. 236–242, 2008.

J. L. Mega, S. L. Close, S. D. Wiviott et al., “Cytochrome P450 genetic polymorphisms and the response to prasugrel: relationship to pharmacokinetic, pharmacodynamic, and clinical outcomes,” Circulation, vol. 119, no. 19, pp. 2553–2560, 2009.

T. Simon, C. Verstuyft, M. Mary-Krause et al., “Genetic determinants of response to clopidogrel and cardiovascular events,” The New England Journal of Medicine, vol. 360, no. 4, pp. 363–375, 2009.

J.-P. Collet, J.-S. Hulot, A. Pena et al., “Cytochrome P450 2C19 polymorphism in young patients treated with clopidogrel after myocardial infarction: a cohort study,” The Lancet, vol. 373, no. 9660, pp. 309–317, 2009.

D. Sibbing, J. Stegherr, W. Latz et al., “Cytochrome P450 2C19 loss-of-function polymorphism and stent thrombosis following percutaneous coronary intervention,” European Heart Journal, vol. 30, no. 8, pp. 916–922, 2009.

B. Giusti, A. M. Gori, R. Marcucci et al., “Cytochrome P450 2C19 loss-of-function polymorphism, but not CYP3A4 IVS10+12G/A and P2Y12 T744C polymorphisms, is associated with response variability to dual antiplatelet treatment in high-risk vascular patients,” Pharmacogenetics and Genomics, vol. 17, no. 12, pp. 1057–1064, 2007.

J. T. Brandt, S. L. Close, S. J. Iturria et al., “Common polymorphisms of CYP2C19 and CYP2C9 affect the pharmacokinetic and pharmacodynamic response to clopidogrel but not prasugrel,” Journal of Thrombosis and Haemostasis, vol. 5, no. 12, pp. 2429–2436, 2007.

P. Fontana, J.-S. Hulot, P. de Moerloose, and P. Gaussem, “Influence of CYP2C19 and CYP3A4 gene polymorphisms on clopidogrel responsiveness in healthy subjects,” Journal of Thrombosis and Haemostasis, vol. 5, no. 10, pp. 2153–2155, 2007.

L. Yao, H. C. Wang, J. Z. Liu, and Z. M. Xiong, “Quantitative assessment of the influence of cytochrome P450 2C19 gene polymorphisms and digestive tract cancer risk,” Tumor Biology, vol. 34, no. 5, pp. 3083–3091, 2013.

D. J. Shin, J. Kwon, A.-R. Park et al., “Association of CYP2C19∗2 and ∗3 genetic variants with essential hypertension in Koreans,” Yonsei Medical Journal, vol. 53, no. 6, pp. 1113–1119, 2012.

S. Sviri, S. Shpizen, E. Leitersdorf, M. Levy, and Y. Caraco, “Phenotypic-genotypic analysis of CYP2C19 in the Jewish Israeli population,” Clinical Pharmacology and Therapeutics, vol. 65, no. 3, pp. 275–282, 1999.

M. R. G. Taylor, D. Slavov, K. Humphrey et al., “Pharmacogenetic effect of an endothelin-1 haplotype on response to bucindolol therapy in chronic heart failure,” Pharmacogenetics and Genomics, vol. 19, no. 1, pp. 35–43, 2009.

O. Harismendy, D. Notani, X. Song et al., “9p21 DNA variants associated with coronary artery disease impair interferon-gamma 3 signalling response,” Nature, vol. 470, no. 7333, pp. 264–270, 2011.

C. Holohan, S. van Schaeybroeck, D. B. Longley, and P. G. Johnston, “Cancer drug resistance: an evolving paradigm,” Nature Reviews Cancer, vol. 13, no. 10, pp. 714–726, 2013.

A. Motterle, X. Pu, H. Wood et al., “Functional analyses of coronary artery disease associated variation on chromosome 9p21 in vascular smooth muscle cells,” Human Molecular Genetics, vol. 21, no. 18, pp. 4021–4029, 2012.

A. Congrains, K. Kamide, R. Oguro et al., “Genetic variants at the 9p21 locus contribute to atherosclerosis through modulation of ANRIL and CDKN2A/B,” Atherosclerosis, vol. 220, no. 2, pp. 449–455, 2012.

Y. Li, L. O. Lévesque, and M. B. Anand-Srivastava, “Epidermal growth factor receptor transactivation by endogenous vasoactive peptides contributes to hyperproliferation of vascular smooth muscle cells of SHR,” The American Journal of Physiology—Heart and Circulatory Physiology, vol. 299, no. 6, pp. H1959–H1967, 2010.

G. León-García, A. Santana, N. Villegas-Sepúlveda et al., “The T1048I mutation in ATP7A gene causes an unusual Menkes disease presentation,” BMC Pediatrics, vol. 12, article 150, 2012.

Z. Tümer, “An overview and update of ATP7A mutations leading to Menkes disease and occipital horn syndrome,” Human Mutation, vol. 34, no. 3, pp. 417–429, 2013.

J. O. Miners, R. A. McKinnon, and P. I. Mackenzie, “Genetic polymorphisms of UDP-glucuronosyltransferases and their functional significance,” Toxicology, vol. 181-182, pp. 453–456, 2002.

R. H. Tukey and C. P. Strassburg, “Human UDP-glucuronosyltransferases: metabolism, expression, and disease,” Annual Review of Pharmacology and Toxicology, vol. 40, pp. 581–616, 2000.

Q. Cai, T. Wu, W. Zhang et al., “Genetic polymorphisms in glutathione S-transferases P1 (GSTP1) Ile105Val and prostate cancer risk: a systematic review and meta-analysis,” Tumor Biology, vol. 34, no. 6, pp. 3913–3922, 2013.

D. Liu, Y. Liu, L. Ran, H. Shang, and D. Li, “GSTT1 and GSTM1 polymorphisms and prostate cancer risk in Asians: a systematic review and meta-analysis,” Tumor Biology, vol. 34, no. 5, pp. 2539–2544, 2013.

N. Masood, A. Yasmin, and M. A. Kayani, “Genetic deletions of GSTM1 and GSTT1 in head and neck cancer: review of the literature from 2000 to 2012,” Asian Pacific Journal of Cancer Prevention, vol. 14, no. 6, pp. 3535–3539, 2013.

T. O. Lankisch, C. Schulz, T. Zwingers et al., “Gilbert's syndrome and irinotecan toxicity: combination with UDP-glucuronosyltransferase 1A7 variants increases risk,” Cancer Epidemiology Biomarkers & Prevention, vol. 17, no. 3, pp. 695–701, 2008.

T. Suzuki, Y. Miki, Y. Nakamura, K. Ito, and H. Sasano, “Steroid sulfatase and estrogen sulfotransferase in human carcinomas,” Molecular and Cellular Endocrinology, vol. 340, no. 2, pp. 148–153, 2011.

K. Riemann, H. Struwe, A. Eisenhardt, B. Obermaier, K. W. Schmid, and W. Siffert, “Characterization of intron-1 haplotypes of the G protein beta 4 subunit gene—association with survival and progression in patients with urothelial bladder carcinoma,” Pharmacogenetics and Genomics, vol. 18, no. 11, pp. 999–1008, 2008.

S.-J. Tsai, C.-J. Hong, and Y.-J. Liou, “Brain-derived neurotrophic factor and antidepressant action: another piece of evidence from pharmacogenetics,” Pharmacogenomics, vol. 9, no. 9, pp. 1353–1358, 2008.

R. Tamm, K. Oselin, K. Kallassalu et al., “Thiopurine S-methyltransferase (TPMT) pharmacogenetics: three new mutations and haplotype analysis in the Estonian population,” Clinical Chemistry and Laboratory Medicine, vol. 46, no. 7, pp. 974–979, 2008.

A. Matakidou, N. Hamel, S. Popat et al., “Risk of non-medullary thyroid cancer influenced by polymorphic variation in the thyroglobulin gene,” Carcinogenesis, vol. 25, no. 3, pp. 369–373, 2004.

A. Mihaljević-Peleš, M. Šagud, N. Božina, M. Živković, and N. Jovanović, “Pharmacogenetics and antipsychotics in the light of personalized pharmacotherapy,” Psychiatria Danubina, vol. 22, no. 2, pp. 335–337, 2010.

S. H. Preskorn, “Clinically relevant pharmacology of selective serotonin reuptake inhibitors: an overview with emphasis on pharmacokinetics and effects on oxidative drug metabolism,” Clinical Pharmacokinetics, vol. 32, supplement 1, pp. 1–21, 1997.

M. Vaswani, F. K. Linda, and S. Ramesh, “Role of selective serotonin reuptake inhibitors in psychiatric disorders: a comprehensive review,” Progress in Neuro-Psychopharmacology and Biological Psychiatry, vol. 27, no. 1, pp. 85–102, 2003.

R. T. Coutts, “Polymorphism in the metabolism of drugs, including antidepressant drugs: comments on phenotyping,” Journal of Psychiatry and Neuroscience, vol. 19, no. 1, pp. 30–44, 1994.