Pharmacogenomics in type II diabetes mellitus management: Steps toward personalized medicine

Abstract

Advances in genotype technology in the last decade have put the pharmacogenomics revolution at the forefront of future medicine in clinical practice. Discovery of novel gene variations in drug transporters, drug targets, effector proteins and metabolizing enzymes in the form of single-nucleotide polymorphisms (SNPs) continue to provide insight into the biological phenomena that govern drug efficacy and toxicity. To date, novel gene discoveries extracted from genome-wide association scans and candidate gene studies in at least four antidiabetic drug classes have helped illuminate possible causes of interindividual variability in response. Inadequate protocol guidelines for pharmacogenomics studies often leads to poorly designed studies, making it hard to formulate a definitive conclusion regarding the clinical relevance of the information at hand. These issues, along with the ethical, social, political, legislative, technological, and economic challenges associated with pharmacogenomics have only delayed its entry to mainstream clinical practice. On the other hand, these issues are being actively pursued and rapid progress is being made in each area which assures the possibility of gaining widespread acceptance in clinical practice.

Keywords:

• pharmacogenomics
• genetics
• pharmacokinetics
• pharmacodynamics
• personalized medicine
• type 2 diabetes
• pharmacotherapy
• antidiabetic drugs
• efficacy, and safety

Introduction

In the previous decade a major component of discussions involving disease prevention, therapy and prognosis has shifted from genetics to genomics.^1–3 The completion of the Human Genome Project in 2003 and advances in genotyping technology have enabled researchers to explore the entire genome, updating their approach from looking at just genes to the genome as a whole, in what are referred to as genome-wide association studies (GWAS).⁴ Novel information regarding individual disease susceptibility, response to therapy and progression path has emerged from this approach at a relatively rapid pace, especially in the past two years.^4,5 This is because a determinant of an individual’s disease susceptibility, response to therapy, and progression are attributed, at least partly, to common variations in certain genes; mainly in the form of single nucleotide polymorphisms (SNPs).⁵ This approach has infiltrated many analytical fields to produce many hybridized genomic fields, which include toxicogenomics, proteomics, and pharmacogenomics.^1,3 However, the focus here will be on how polymorphisms in an individual’s genome can alter a response to drug therapy, an area referred to as pharmacogenomics.

Pharmacogenomics is an approach that has evolved from pharmacogenetics and become a new scope for the pharmaceutical and biomedical fields.^1,3 Pharmacogenomics promises a new drug selection process, which takes into account variations in an individual’s genetic makeup to optimize pharmacokinetics and pharmacodynamics to ultimately increase drug efficacy and safety profile.⁶ In other words, creating genetically tailored drug regimens to optimize an individual’s response.^1,6 It is estimated that about 95% of the variability in drug response is due to genetic differences; accounting for these differences then would be highly beneficial, not only for the health care industry, but for patients themselves, decreasing the burden of treatment failures and adverse events on society.^7,8

The pharmacogenomics approach has already supplied researchers with candidate genes and their translational ramifications on drug response in many complex disease states.² A popular focus of several GWAS and candidate gene analysis has been dedicated to understanding the implications of genetic variation in individuals diagnosed with type 2 diabetes mellitus (T2D).⁹ In recent years, a great deal of progress has been made in the form of establishing a genetic explanation for etiological mechanisms by which T2D develops and interindividual variability in response to standard therapy.¹⁰

According to 2007 surveys, the prevalence of diabetes is estimated at 246 million, which is roughly a 63% increase since the turn of the 21st century. Perhaps more alarming is of the prediction of an increased prevalence in years to come, with an estimate of 380 million cases in the year 2025.¹¹ In the United States alone, the prevalence of diagnosed and undiagnosed diabetes cases was 7.8% (23.6 million) of the population, 90% to 95% of which are T2D.¹¹ Additionally, DM is the leading cause of blindness, end-stage renal failure and lower extremity amputations, and increases the risk for many other life-threatening conditions such as heart disease and stroke by two to four times that of a nondiabetic.¹² Financially speaking, an estimated one out of every five health care dollars are spent on diabetic care in the United States and about 174 billion total health care dollars spent annually.¹³

These numbers can be partly attributed to an aging population, unhealthy diet, increasing rates of obesity and sedentary lifestyles.¹¹ Although obesity rates are increasing among all age groups, a significant contributor to these overwhelming prevalence trends is the increase in childhood obesity rates in developing countries.^11,14 Therefore, it comes as no surprise that behavioral modifications (diet and exercise) are the initial recommendation for managing T2D.^15,16 But when euglycemia cannot be achieved through dietary and exercise modification alone, supplementation with oral hypoglycemics and insulin, in some cases, is recommended.¹⁶ Type 2 DM is treatable with ten drug classes: biguanides, sulfonylureas (SUs), thiazolidinediones (TZD), meglitinides, α-glucosidase inhibitors, amylin mimetics, glucagon-like peptide-1 (GLP-1) mimetics, dipeptidyl peptidase-4 (DPP-4) inhibitors, and insulin.^16,17 However, there is a significant amount of variability in responses to these antidiabetic agents. Recent pharmacogenomic studies have linked anti-diabetic drug response variability to genetic polymorphisms in transporters, targets, effector proteins and metabolizing enzymes. The purpose of this article is to explain how pharmacogenomics aims to achieve the best therapeutic outcome based on genomic insights, to ultimately lead to the development of “personalized medicine”.⁷

Biguanide class (Metformin)

Metformin’s place in therapy and mechanism of action

Introduced to the United States market in 1995, metformin has become a first-line therapy for the management of T2D and is widely prescribed in the United States.¹⁷ By 2006, metformin accounted for 37% of noninsulin prescriptions in the United States alone.¹⁷ Although the full pharmacological action profile remains unclear at this point, evidence suggests that the major action of metformin is exerted in the liver, primarily from the activation of adenosine monophosphate-activated protein kinase (AMPK). Activation (phosphorylation) of AMPK is done by serine-threonine kinase (LKB1), which suppresses hepatic gluconeogenesis, thereby, reducing glucagon-mediated glucose output by the liver.^17,18 Metformin is able to exert glucose-lowering actions with a low risk of hypoglycemia, as well as reduce the likelihood of developing macrovascular (hypertension, atherosclerosis) and microvascular (neuropathy, retinopathy, nephropathy) complications, these qualities are what make it an attractive first-line therapy.^17–19 However, metformin is not without adverse events such as diarrhea and nausea that occur in about 30% of patients; or a more serious but very rare side effect, lactic acidosis.^17,18 Despite an exceptional efficacy and safety profile, several T2Ds (about 38%) still fail to reach glycemic goals in metformin therapy.¹⁷

Recent studies suggest that interpatient variability in response to metformin therapy could be related to polymorphisms in the organic cation transporter (OCT) genes and/or the multidrug and toxin extrusion (MATE) genes (Table 1).^{17,18,20–22} The SLC22A1 and SLC22A2 genes encode for OCT1 and OCT2, respectively.^22,23 OCT1 is primarily expressed in hepatocytes and enterocytes and mediate the uptake of metformin in these cells.^17,24 Expressed predominately in the renal distal tubule, OCT2 functions to facilitate urinary elimination many xenobiotics, including metformin.¹⁸ The SLC47A1 and SLC47A2 genes encode for the MATE1 and MATE2-K proteins, respectively.²⁵ Originally, the MATE proteins were found in bacteria and function to sustain cell equilibrium by excreting metabolic wastes and xenobiotics, but recently they have been shown to play a similar role in humans.²⁶ The MATE1 protein is expressed in the bile canalicular membrane of hepatocytes and in the renal epithelium, whereas MATE2-K is expressed only in the renal epithelium.²⁵ It has been demonstrated that metformin is an exceptional substrate for MATE1 and MATE2-K and that these proteins play a role in the elimination of metformin into the bile (MATE1) and urine (MATE1 and MATE2-K).^25,27

Table 1 Possible role and effect of variations on metformin therapy

	SLC22A1	SLC22A2	SLC47A1	SLC47A2
Codes	OCT1	OCT2	MATE1	MATE2-K
Site	Hepatocytes, enterocytes	Renal distal tubule	Bile canalicular membrane of hepatocytes and in renal epithelium	Renal epithelium
Action	Mediates metformin uptake, accumulation and pharmacological action in the liver (AMPK activation)	Facilitates urinary elimination of metformin	Metformin secretion in bile and urine	Metformin excretion in urine
Variants	Decreased hepatic and intestinal metformin uptake and accumulation. Decreased AMPK activation. Increased plasma glucose levels and lower insulin levels	Decreased metformin clearance and increased plasma concentration	May affect glucose lowering effects of metformin. With each A allele there is an associated 0.3% HbA1C reduction. G allele has no association with HbA1C reduction

Role of SLC22A1 polymorphisms in metformin response

Wang and colleagues observed a 30-fold decrease in hepatic uptake and a 3- to 7–fold decrease in intestinal (duodenum, jejunum, ileum) uptake of metformin in Oct1^−/− mice compared to Oct1^+/+ mice.²⁸ However, the renal distribution and elimination of metformin was found to be nearly equal when the two mice groups were compared.²⁸ These results demonstrate that OCT1 is involved in hepatocyte and enterocyte uptake and that kidney uptake relied minimally, if at all, on OCT1.²⁸ Since kidney uptake was not affected by the presence or absence of OCT1, the authors concluded that distribution to the kidney was achieved primarily by another transport mechanism.²⁸ The major findings of this study were that hepatic and intestinal uptake was accomplished through OCT1, but whether OCT1 had any influence on metformin’s response or toxicity (ie, nausea, diarrhea, lactic acidosis) remained in question.²⁸

In 2007, Shu and colleagues sought to discover whether OCT1 plays a role in the actions of metformin and the associated implications on drug response.²¹ Therefore, investigators initially looked at whether AMPK activation was associated with OCT1 expression.²¹ In order to evaluate this, metformin-stimulated AMPK activation of cells transfected with human OCT1 (HEK-OCT1 cells) were compared to control cells (empty vector-transfected cells).²¹ Metformin accumulated considerably more in the HEK-OCT1 cells, which also corresponded with a higher AMPK activation than that of the empty vector cells.²¹ This demonstrates that OCT1 is needed for metformin uptake and pharmacological actions in the liver in vitro. Next, investigators compared the in vivo uptake of metformin in Oct1^−/− mice with cells containing a functional Oct1 allele (Oct1^+/+ and Oct1^+/−) after a single oral dose of metformin.²¹ The results from this experiment were consistent with Wang and colleagues and other research that showed the hepatic uptake in wild-type mice (Oct1^+/+ and Oct1^+/−) was significantly (4.2-fold) higher in comparison to the Oct1^−/− mice.^21,28 Additionally, metformin suppressed glucagon-stimulated glucose production in the wild-type mice, but not in the Oct1^−/− mice.²¹ Subsequently, the Oct1^−/− and wild-type mice were fed the same high-fat diet for eight weeks to increase fasting plasma glucose (FPG) levels and were subjected to a five-day treatment with either saline or metformin.²¹ Wild-type mice receiving metformin therapy illustrated a reduction in FPG by over 30%, but Oct1^−/− mice that received metformin did not have a notable reduction in FPG levels.²¹ In summary, Shu and colleagues revealed that metformin accumulation and AMPK activation are dependent on functional OCT1 expression and that the level of AMPK activation corresponded to suppression in glucagon-stimulated glucose production and reduced FPG levels.²¹ The final assessment was to determine whether OCT1 polymorphisms altered metformin response in humans by measuring metformin uptake in 12 known OCT1 variants (S14F, R61C, F160L, S189L, G220V, P341L, R342H, G401S, V408M, 420del, G465R, R488M) and comparing them to an OCT1-reference using HEK293 cells. This resulted in the identification of seven reduced function OCT1 variants (S14F, R61C, S189L, G220V, G401S, 420del, G465R; p < 0.001) in terms of metformin uptake and AMPK activation in transfected HEK293 cells. Next, humans expressing these genotypes were studied by measuring plasma glucose levels using a three-hour oral glucose tolerance test (OGTT), after receiving oral glucose and metformin treatment. The results showed that metformin in reference alleles produced a significant reduction in plasma glucose levels relative to the reduced function alleles. In addition, reference alleles were associated with significantly (p < 0.05) higher insulin levels two hours after the OGTT administration when compared to the reduced function alleles. Overall, the information derived from this study demonstrates that polymorphisms of OCT1 genes may modulate an individual’s clinical metformin response in humans. However, further studies are needed to replicate this effect in larger populations and in various ethnic groups to determine if the effects can be generalized to the overall population of the reduced function allele carriers. It may also be beneficial if studies to correlate SLC22A1 polymorphisms with metformin toxicity.

On the other hand, a study examining whether OCT1 or OCT2 polymorphisms had any impact on the efficacy of metformin had different outcomes than other work in the area.²⁰ In the study, Japanese subjects (n = 33) were categorized based on response to greater than three months of metformin therapy (24 responders and 9 nonresponders) and no remarkable differences in allelic frequency and hepatic OCT1 mRNA levels among the two groups were revealed, with the exception of the 408 Met allele.²⁰ The 408Met allele showed a tendency to be higher in the nonresponders and also corresponded to a lower level of hepatic OCT1 mRNA levels, but this difference was deemed not significant.²⁰ When OCT2 was analyzed, investigators found that OCT2 gene variation occurs less often than with other known drug transporters (ie, MDR1, MRP1, MRP2, OATP-C), which is an observation that is consistent with other research results.²⁰ There were no significant differences in the prevalence of OCT2 gene variation among the two groups either.²⁰ Possible reasons for the disagreement between Shikata and colleagues and Shu and colleagues may be contributed ethnic differences, small sample size (n = 33) and other concurrent medications (antihyperglycemics as well as others not listed) may have altered the participants response to metformin. The limitations of this study make it hard to arrive at a definitive conclusion, which necessitates further work in this area to confirm whether OCT1 is a mediator of metformin response.

Role of OCT2 polmorphisms in determining metformin response

Polymorphisms in the OCT2 gene have also been identified and shown to cause inter-patient variability in the pharmacokinetic and pharmacodynamic profile of metformin.^18,23 Leabman and colleagues showed that four gene variants (M165I, A270S, R400C, K432Q) showed altered uptake of 1-methyl-4-phenylpyridinium (MPP+), which is a high-affinity substrate for OCT2 [¹⁸, ²³]. Three (M165I, A270S, R400C) of the four showed revealed a reduced MPP + uptake capacity and the other (K432Q) showed an increased uptake activity.¹⁸ Two studies have been done to show the relationship between metformin clearance in individuals carrying three variants (c.596C > T, c.602C > T, c.808G > T) of the SLC22A2 gene. One study was done in Korean subjects (n = 26) and the other was done in Chinese subjects (n = 15), but both studies showed that metformin clearance was reduced and plasma metformin concentrations were increased in carriers of these variants.^23,29 However, the study using Chinese subjects also demonstrated that drug–drug interactions between metformin and cimetidine also depends on these polymorphisms.²⁹ Results with concomitant cimetidine use showed that renal secretion and clearance of metformin were reduced in for all variants, but particularly reduced for the homozygous TT variant carriers (18.7% of the normal metformin clearance).²⁹ Together, the results from these studies exemplify that the clearance and plasma concentrations of metformin is influenced by OCT2 polymorphisms and dose modifications based on the particular variants may be beneficial to maximize the efficacy and minimize the toxicity of metformin.^23,29 However, there are no studies that relate SLC22A2 polymorphisms with metformin toxicity (ie, lactic acidosis) nor is there enough evidence to develop a dose-adjustment scheme to manage patients with these OCT2 variants. Although the results are consistent between the two studies, further investigation with larger and more diverse sample populations need to be done to confirm these results and determine the need for developing a dose-adjustment algorithm for managing these patients. Additional studies to determine the specific structural-activity relationships associated with these polymorphisms may provide useful information on how to manage these polymorphisms to optimize metformin therapy.

Role of MATE1 polymorphisms in metformin response

Becker and colleagues. looked at one polymorphism of the SLC47A1 gene (rs2289669G > A) and the associated reduction in HbA_1C in subjects (n = 116) treated with metformin between 1997 and 2008.²² The results showed that with each A-allele came a 0.3% larger reduction in HbA_1C (0.6% larger reduction than expected for homozygous rs2289669AA) after a 30-day treatment with metformin compared to the reference allele.²² In contrast, the homozygous G allele showed no reduction in the same context.²² HbA_1C changes were based upon a value within 90 days (with an average time prior metformin initiation of 12 days and a standard deviation of ±16 days) of starting metformin therapy and 30 days after continuous metformin use.²² This demonstrates that SNPs in the SLC47A1 may affect the glucose-lowering affects of metformin.²² However, this was the first study to assess this association and further studies are needed to replicate this information in larger and diverse populations to determine if there are any associated clinical implications.

Sulfonylureas

Role of K_ATP channels in stimulating insulin secretion in pancreatic β-cells

Pancreatic β-cells are equipped with ATP-sensitive potassium channels (K_ATP channels) that close in response to elevated cytosolic ATP/ADP concentrations, produced from glucose metabolism in the setting of elevated blood glucose. These K_ATP channels are composed of a Kir6.2 pore and SUR subunits, which regulate the open and closure status of the Kir6.2 pore based on ATP levels and is the target for the SU drug class.^7,30 Closure leads to membrane depolarization and opening of intracellular voltage-gated Ca²⁺ channels that results in an increase in intracellular [Ca²⁺], which stimulates insulin vesicle exocytosis and insulin release to reduce blood glucose.^7,10,30,31

Sulfonylurea mechanism of action

Sulfonylurea exert their pharmacological actions through specifically binding the regulatory SU receptor (SUR1) subunit of the K_ATP channel to induce transformation to the closed state; thereby promoting insulin secretion and reduction in blood glucose.³¹ In addition, research has also proposed that SUs may also bind to the SUR2A subunit on peripheral K_ATP channels to increase glucose uptake in skeletal muscle and adipose tissues, but this action remains to be defined.³¹

Sulfonylureas’ place in therapy and possible reasons for interindividual variability in response

The SUs are an alternative to metformin as a first-line option for the management of T2D.⁷ However, between 10% and 20% of patients experience a condition referred to as primary failure of SU treatment, which is characterized by a reduction in insulin secretion within the first three months of initiation.¹⁶ Explanations for this variability include polymorphisms in the genes that encode for the Kir6.2 and SUR1 subunits (KCNJ11 and ABCC8, respectively) of the K_ATP channel and in metabolizing enzyme, CYP2C9.³⁰

Several studies demonstrate genetic variations among the genes that encode for the Kir6.2 pore (KCNJ11 gene) and the SUR1 subunits (ABCC8 gene) of pancreatic K_ATP channels can alter response to SU therapy, but they were mainly observed in the neonatal diabetic population because that is where these polymorphisms are predominately found.^30,31 Neonatal diabetes is a monogenic form of diabetes typically diagnosed during the first three months of life and subdivided into a form that that usually reconciles after 18 months, transient neonatal diabetes (TND), and a lifelong form, persistent neonatal diabetes (PND).³¹ These forms of diabetes can be attributed to genetic variations in the K_ATP channel, which results in a decrease in ATP-sensitivity and/or a pore fixed in the open conformation.³¹ In general, KCNJ11and ABCC8 gene polymorphisms alter the K_ATP channels sensitivity to ATP or favor an open state leading to a reduction in insulin secretion and subsequent hyperglycemia in these individuals (Table 2).^31–33

Table 2 Possible role and effects of variations on sulfonylurea treatment goals

	KCNJ11	ABCC8	CYP2C9
Codes	Kir6.2	SUR1	Enzyme in cytochrome P450
Site	KATP channels of pancreas		Liver
Action	SU blocks conductance of KATP channels by interacting with SUR1 and Kir6.2 subunits		Metabolize SU
Variants	Decreased ATP sensitivity and/or pore fixed in the open conformation leading to decreased insulin secretion and subsequent hyperglycemia. Different variations result in different sensitivities to SU	Ser1369Ala variant can influence SU efficacy	Reduced clearance which may correlate to SU toxicity

Abbreviations: ATP, adenosine triphosphate; SU, sulfonylurea.

KCNJ11 and ABCC8 polymorphisms in neonatal diabetes and the associated effects on sulfonylurea response

One study categorizes the variations to those that directly affect ATP sensitivity and those that indirectly affect ATP sensitivity and examine the associated response to SU therapy.³¹ Genetic variations that directly affect ATP sensitivity (I182V; ABCC8) do so by decreasing K_ATP channel affinity for ATP, which causes the K_ATP channel to remain open that translates to a reduction in insulin secretion.³¹ Indirect genetic variations (Q52R, V59M, I296L; KCNJ11) increases the stability of the K_ATP channel in the open conformation, which has been previously shown to induce functional detachment of the SUR1 subunit from the pore-forming Kir6.2 subunits.³¹ SUs block the conductance of K_ATP channels by interacting with the SUR1, referred to as the high-affinity site, as well as Kir6.2 subunits referred to as the low-affinity site.³¹ The functional detachment of the SUR1 from the Kir6.2 subunit leads to a decreased sensitivity to SUs (in this study, glibenclamide and tolbutamide) by reducing the high-affinity block.³¹ The results show that glibenclamide blocked the conductance K_ATP channel almost completely in the wild-type, but the variant cells were relatively insensitive to glibenclamide and I296L carriers essentially unresponsive to glibenclamide treatment. With regard to tolbutamide treatment, mild reductions in current were seen in Q52R, V59M, and I183V, but I296L was completely insensitive to the effects of tolbutamide. These results indicate that individuals carrying the I296L variant are less sensitive to SU therapy and may require higher doses of SU therapy to achieve the same goals compared to the other three variants or may even require insulin therapy. However, the other three variants (Q52R, V59M, I183V) may simply need a higher SU dose to achieved goals and the dose would be dependent on the variant expressed. Although, a definitive dosing guide is not available presently, the rank of sensitivities to either SU (I182V > V59M > Q52R > I296L) may provide a general dosing algorithm to help adapt the dose to the particular variant.

There are several other KCNJ11 variants (which include F333I, F35V, R201H, R201C, Q52R, I296L, L164P, G53S, G53R) and ABCC8 variants (which include I182V, H1023Y, I1424Y, F132L) that increase susceptibility to neonatal diabetes, but in the majority of the variant genotypes seem to have a preserved response to SU therapy with the exception of Q52R, I296L, and L164P variants.^7,31,33–38 KCNJ11 variant carriers may only need higher doses of SU therapy to achieve therapeutic goals. Since, SUs are not effective for all variant carriers genetic screening for these variants would be a valuable tool for selecting SU therapy candidates.⁷ Alternatively, further investigations are needed to develop a dosing algorithm for those with KCNJ11 and ABCC8 polymorphisms that respond to SU therapy, but require dose modifications.

Patients with permanent neonatal DM and mutations affecting Kir6.2 response to SU.^33,39 suggest that many of these patients may successfully be switched to oral drugs (such as SU) instead of insulin. Pearson and colleagues reported a 90% success rate (44 of 49 patients) in switching from subcutaneous insulin injections to oral glyburide.³⁹ Babenko and colleagues had five of nine patients successfully switch from insulin injection to oral SU therapy with permanent neonatal DM, in these cases caused by mutations in the SUR1 subunit.³² This instance of patients with permanent neonatal DM with mutant Kir6.2 or SUR1 to SU’s illustrates just how powerful genetics may be in identifying patients who may or may not benefit from a certain treatment. It may be of value to test all those with neonatal diabetes for mutations in Kir6.2 and SUR1.

Role of KCNJ11 and ABCC8 polymorphisms in sulfonylurea response in type 2 diabetics

A study by Feng and colleagues with a total of 1,268 Chinese T2D patients was conducted. These patients were all recently diagnosed diabetics (within the past five years) and had had no recent hypoglycemic treatment. They were treated with gliclazide for eight weeks, and fasting and oral glucose tolerance test two-hour plasma glucose, fasting insulin and A1C were measured at the start and after eight weeks of treatment. After eight weeks of treatment it was found that patients with Ser1369Ala of the ABCC8 gene and rs5210 of the KCNJ11 gene had a significant decrease in FPG (p = 0.002). This was also confirmed in cohort 2 of the study and in the pooled analysis subjects with the Ala/Ala genotype had a 7.7% greater decrease in FPG (p < 0.001), a 3.5 greater decrease in A1C (p = 0.06) and a 11.9% greater decrease in two-hour plasma glucose (p = 0.003) than those with the Ser/Ser genotype after eight weeks of gliclazide treatment.⁴⁰

Role of CYP2C9 polymorphisms in sulfonylurea response

While genetic alterations in the transporting and target proteins can produce variability in response, drug metabolizing enzyme activity may also contribute the response variability. With respect to four United States available (Tolbutamide [Orinase^®], Glipizide [Glucotrol^®], Glibenclamide [Glyburide^®] and Glimepiride (Amaryl^®]) SUs, the major metabolizing enzyme is cytochrome P450 2C9 (CYP2C9).¹⁰

Kirchheiner and colleagues investigated three allelic forms (CYP2C9*3/*2/*1) of CYP2C9 and compared there hepatic metabolism to the wild-type (CYP2C9*1/*1).¹⁰ The CYP2C9*3/*3 revealed lower clearances of tolbutamide, glibemclamide, and glipizide (16%, 50%, and 20%, respectively) compared to that of the wild-type.¹⁰ Homozygous CYP2C9*2/*2 showed a mildly reduced clearance of tolbutamide and glibemclamide (75% and 90%, respectively) compared to the wild-type, but this difference was not deemed statistically significant.¹⁰ Heterozygous condition CYP2C9*1/*2 showed a similar clearance with respect to the wild-type for tolbutamide and allele *2 carriers showed a mildly reduced clearance for glibenclamide compared to the wild-type (homozygous *2/*2 carriers showed 90% clearance compared to wild-type). Both heterozygous CYP2C9*1/*3 and CYP2C9*2/*3 showed a 50% reduction in the clearance of tolbutamide, but and the reduction was less pronounced for glibenclamide. Overall, the wild-type, CYP2C9*2/*2, and *1/*2 genotypes showed roughly normal (CYP2C9*2/*2) to extensive (CYP2C9*1/*2) clearance rates, CYP2C9*1/*3 and *2/*3 genotypes showed intermediate clearance rate and the CYP2C9*3/*3 genotype clearance rate was significantly reduced in participants. This data suggests that pre-screening individuals for these polymorphisms may help guide dosing, which may limit toxicity in these individuals and maximize efficacy; although, correlations between enzyme activity and toxicity was not demonstrated making it hard to rationalize dose adjustments. Therefore, more studies replicating this data in larger and ethnically diverse populations are needed.

In a second study, Caucasian volunteers (n = 475) expressing CYP2C9 variable function alleles and receiving SU therapy (tolbutamide 62.3%, glibenclamide 16.2%, glimepiride 16%, glicazide 5.5%) were analyzed based on changes in dose from their first prescription to the tenth prescription and changes in blood glucose in a period 90 days before and 180 days after beginning SU therapy (n = 79).⁴¹ In contrast to the previous study, clinically relevant variation in SU pharmacokinetics and response in participants with different CYP2C9 polymorphisms were only found in the tolbutamide group and not in any other treatment groups.⁴¹ The tolbutamide group carrying the *3 allele showed only a minute change in dose between the first and the tenth prescriptions compared to that of the wild-type (CYP2C9*1/*1; 12 mg versus 279 mg, respectively).⁴¹ Additionally, 15 of the 20 individuals expressing the *3 allele receiving tolbutamide showed no change in dose from the first prescription to the tenth prescription.⁴¹ No differences between tolbutamide group expressing the *1/*2 and *2/*2 genotypes were found compared to the wild-type group, indicating that that metabolism was not necessary in these individuals, which is consistent with Kirchheiner and colleagues. In a supplementary investigation, the change in dose from the sixth to twentieth prescription was compared to the first prescription using repeated measurements.⁴¹ These results showed that *3 allele carriers receiving tolbutamide showed a daily prescribed dose of 316 mg lower compared to the wild-type tolbutamide group (p = 0.0008).⁴¹ Lastly, in the secondary objective analysis (that measured changes in FPG between a period 90 days prior and 180 after staring SU therapy) no reductions in the tolbutamide groups, regardless of allele, compared to the wild-type reached statistical significance. The findings of this data clearly show that there is a relationship between the CYP2C9*3 allele carriers and their corresponding prescribed dose of tolbutamide, but this relationship may not be extended to the glucose lower potential of SU therapy.⁴¹ However, the relationship between CYP2C9 polymorphism and SU toxicity may have clinical relevance and should be considered in future studies.⁴¹

Role of TCF7L2 variation in sulfonylurea response

Transcription factor 7-like 2 (TCF7L2) variants are associated with T2D risk, it is thought due to decreased β-cell function. A study by Pearson and colleagues in which 4,469 T2D patients were determined to have the TCF7L2 genotype, of these 901 were incidentally being treated with SU and 945 with metformin. A1C pretreatment was available for a subset of these patients (579 of the patients treated with SU and 755 of those on metformin). A logistic regression was used and treatment failure was defined as A1C of >7% within 3–12 months of treatment initiation. The study found that carriers of the risk allele were less likely to respond to SU’s (odds ratio [OR] for failure 1.95%; p = 0.005, comparing rs12255372 T/T vs G/G). When the A1C baseline was included the association was strengthened (OR 2.16; p = 0.009). In conclusion the researchers found that TCF7L2 variants did in fact influence patient response to SUs, but not metformin.⁴²

Thiazolidinedione

Thiazolidinedione place in therapy

The TZDs, otherwise referred to as glitazones, are included in the battery of type 2 antidiabetic agents that are prescribed as an adjunct to diet and exercise (monotherapy) or in combination with SUs, metformin, or insulin.⁴³ The TZDs currently available to the United States market are pioglitazone (Actos^®) and rosiglitazone (Avandia^®).⁴³ TZDs have shown to increase insulin sensitivity of myocytes, adipocytes, and hepatocytes as well as reducing FPG and HbA_1C levels in several, but not all individuals.⁴³ An estimated 20% and 30% of individuals fail to reach goal insulin sensitivity and glycemic measures following TZD therapy.⁴³ Studies have shown that the interindividual variability in TZD response may be attributed, at least partially, to genetic inheritance.⁴³ Genetic polymorphisms in TZD target proteins, effector proteins, and metabolizing enzymes have been shown to contribute to variability in response to TZDs in a number of studies (Table 3).⁴³

Table 3 Possible role and effects of variations on TZD therapy

TZD therapy	ACDC	S447S
	Adiponectin	LPL
Action	Protective, plays role in glucose uptake and fatty acid oxidation. Decrease in adiponectin associated with insulin resistance and atherosclerosis	Hydrolyze triglycerides and mediate clearance of atherogenic lipoproteins
Variants	+45GG allele smaller decreases in FPG level, A1C and increase in adiponectin than those with +45T allele Similarly the +276GG genotype when compared to the +276T allele	S447X shows reduced response

Thiazolidinedione mechanism of action

TZDs pharmacologically act by agonizing the nuclear perixisome proliferator-activated receptor gamma (PPAR-γ) to form heterodimers with the retinoid X receptors (RXR), which binds the PPAR-response elements and regulates target gene transcription.⁴³ The heterodimer (PPAR-γ + RXR) complex plays a role in regulating the transcription of genes involved in fatty acid uptake and storage, glucose uptake, adipocyte differentiation, and cytokine production and was recently shown to regulate the expression OCT1 via the fatty acid binding protein 2 (FABP2) gene.^43,44

Role of PPAR-γ in determining type 2 diabetes susceptibility

Discovery of a polymorphism in the PPAR-γ2, Pro12Ala, found mainly in the adipocytes has been studied extensively in attempt to provide some plausible insight into TZD interindividual response variability.^43,45 It is estimated that Pro12Ala has a frequency of roughly 12% to 15% in the general population, making it the most common PPAR-γ variant.⁴⁶ Two functional studies show that the TZD-activated PPAR-γ2 Ala12 allele exhibited a reduced binding affinity for the PPAR promoter region.^47,48 A meta-analysis showed that Ala12 allele displayed a protective effect with regard to developing T2D, reducing the risk in carriers by 19%.⁴⁵ However, the risk reductions were more pronounced in Asians (35%) than in North Americans (18%) or Europeans (15%), which progressively diminished from northern (10%) to southern Europe (0%). Furthermore, the protective effect of the Ala12 allele appeared to be inversely related to body mass index.⁴⁵

Role of PPAR-γ polymorphisms in thiazolidinedione response

The question of how PPAR-γ polymorphisms affects response to TZD therapy was addressed by Bluher and colleagues and Kang and colleagues, but yielded inconsistent results.^46,49 Bluher and colleagues conducted a study in German T2D subjects (n = 131) using pioglitazone and determined that there was no significant difference in response to piogltiazone between the heterozygous genotype (Pro12Ala) and the homozygous genotype (Pro12Pro).⁴⁶ With regard to the homozygous carriers, the study was underpowered to detect a difference between the response rates compared to Pro12Ala and Pro12Pro genotypes.⁴⁶ Additionally, individuals with heterozygous and homozygous Pro12 genotypes showed a significant (p < 0.05) reduction in triglycerides and carriers of the Ala12Ala genotype showed a significant (p < 0.05) elevation in high-density lipoprotein levels after six months of pioglitazone therapy.⁴⁶ On the other hand, Kang and colleagues conducted a similar study in Korean T2Ds receiving rosiglitazone for 12 weeks and found the response rate of the Pro12Al genotype was significantly (p = 0.002) higher in comparison to the homozygous Pro12 genotype (86.7% and 43.7%, respectively).⁴⁹ As in the previous study, the homozygous Ala12 was not observed.⁴⁹ Based on the results of each study, it would be hard to arrive at a convincing conclusion to explain if response rates depend on Pro12Al or not. Potential factors leading to inconsistent results may be attributed to differences in ethnicity, environmental factors, selected TZDs, and differences in baseline characteristics.⁴⁹ Additionally, insulin sensitivity was restored to a greater degree in the heterozygous (Pro12Ala) versus the homozygous (Pro12Pro) participants, which was inconsistent with other studies using troglitazone for 12 weeks and one year.^49–51 Possible factors leading to inconsistencies between these three trials include differences in the ethnicity of study populations, treatment drug and duration and method for measuring insulin sensitivity.

Role of adiponectin in thiazolidinedione response

Alternatively, effector protein, adiponectin, was shown to be increased in individuals receiving TZD treatment and found to play a role in the insulin-sensitizing effect of TZDs.⁵² Encoded by the ACDC gene, adiponectin is secreted into circulation from the adipose tissue and has been seen to be associated with insulin resistance and atherosclerosis in low systemic concentrations and is protective at high concentrations.⁵² This may be contributed to the fact that adiponectin modulates metabolic processes, including glucose uptake and fatty acid oxidation.⁴³ Studies have revealed several SNPs in the encoding gene, ACDC.⁵² Of those that were discovered, three (SNP45, SNP276, SNP11377) are associated with T2D, circulating adiponectin levels, and insulin resistance in Japanese individuals and two (SNP11391, SNP11377) in French individuals.⁵² One study examined the effects of two SNPs (SNP45 [TT/GG/TG], SNP276 [TT/GG/TG]) on the response rate in Korean individuals receiving rosiglitazone 4 mg for 12 weeks.⁵² The result indicated rosiglitazone significantly reduced FPG levels, HbA_1C levels and significantly raised adiponectin levels after 12 weeks, without considering genotype.⁵² Taking into account the individual SNPs and genotypes of the individuals, investigators found that carriers of the +45GG genotype had significantly smaller decreases in FPG level, HbA_1C and increases in adiponectin levels compared to those who carried the +45T allele.⁵² Similarly, the +276GG genotype showed the same profile when compared to the +276T allele participants.⁵² Furthermore, the 45/276GG haplotype followed suit in showing smaller reductions in FPG and HbA_1C and smaller increases in adiponectin levels relative to the other haplotypes following rosiglitazone therapy in the 166 subjects tested.⁵² In conclusion, the 45/276GG haplotype can reduce an individual’s response to rosiglitazone therapy and may be a clinically relevant consideration for rosiglitazone initiation, but the question of whether these results can be extrapolated to pioglitazone remains to be elucidated along with the functional mechanism by which this occurs.

Role of lipoprotein lipase in thiazolidinedione response

Lipoprotein lipase (LPL) is another effector protein of TZD response and the expression has been shown to increase after pioglitazone treatment.⁴³ LPL is an enzyme that functions to hydrolyze triglyceride and mediate the clearance of atherogenic lipoproteins. A C→G polymorphism at the position 1595 in the LPL gene results in a premature stop codon at number 447, referred to as S447X.⁴³ In one study, the association between the S447X polymorphism and pioglitazone response was analyzed in Chinese individuals with T2D.⁴³ The results of this study indicated that the homozygous genotype, S447S, was associated with a higher percentage of patients achieving the response criteria compared to the other genotypes.⁴³ In contrast, the S447X show a 0.538-fold reduction in response to pioglitazone treatment compared to the S447S genotype.⁴³ This suggests that LPL polymorphisms are a useful predictor in determining how well an individual will respond to pioglitazone therapy.

Common polymorphisms associated with thiazolidinedione toxicity

While efficacy is an important part of drug response, toxicity may negate the efficacy benefits in some cases making it imperative that it is considered along with the efficacy component. TZDs have displeasing adverse effects, which often outweigh the utility of their benefits and limit their use. Adverse events include weight gain, edema, heart failure, bone fractures in women, and idiosyncratic hepatotoxicity.⁴³

Perilipin and thiazolidinedione-induced weight gain

Weight gain was addressed in a study evaluating the SNPs in the candidate gene, perilipin.⁴³ Four common polymorphisms (6209G/A, 11482G/A, 13041A/G, 14995A/T) were examined in this study.⁴³ Results showed that the 11482G/A variant is significantly associated with weight gain following 12 weeks of rosiglitazone therapy.⁴³ The most pronounced weight gain occurred in the homozygous 11482G/G carriers relative to G/A and A/A (1.33, 0.85, and 0.03 kg, respectively).⁴³ This suggests that A/A genotype carriers are less likely to exhibit weight gain.⁴³

Role of Pro12Ala genotype in thiazolidinedione-induced edema

In terms of edema, there is one study that association between the Pro12Ala polymorphism and edema risk in T2Ds.⁴³ Ragalitazar is a dual PPAR-γ and PPAR-α activator that improves FPG and HbA_1C comparable to pioglitazone, but displays more positive effects on the lipid profile on a study performed on 177 subjects.⁵³ Data from this study showed a 4.4-fold decrease in the risk of edema in participants carrying the Pro12Pro genotype compared to those that had the Ala12 allele, suggesting that the Ala12 may be less prone to TZD-induced edema. However, whether these results may be extrapolated to the TZDs remains to be determined, even when it is known that ragaglitazar acts via a similar mechanism.⁴³

Role of CYP2C9 polymorphisms in thiazolidinedione-induced hepatotoxicity

A Japanese retrospective study looked at the role of CYP2C19 or CYP2D6 enzyme polymorphisms in troglitazone-induced hepatotoxicity in a study that aimed to define the cause of TZD-induced hepatotoxicity.⁴³ Retrospective analyses revealed that the group with troglitazone-induced hepatotoxicity was comprised of a larger percentage of CYP2C19-poor metabolizers compared to the group who did not experience hepatotoxicity.⁴³ However, rosiglitazone and pioglitazone have a much lower likelihood of causing TZD-induced hepatotoxicity making it hard to generalize the results to them.⁴³ Lastly, the role of CYP2C19 in troglitazone metabolism was also undetermined.⁴³

A chief contributor in the development of any to concentration-dependent toxicity is the integrity of an individual’s metabolizing enzymes. TZDs are mainly metabolized by CYP2C8, but CYP2C9 and CYP3A4 also contribute to lesser extents.⁴³ Rosiglitazone and pioglitazone are primarily metabolized by CYP2C8, but CYP2C9 and CYP3A4 contribute to a lesser extent in the metabolism of rosiglitazone and pioglitazone, respectively.^10,43 The studies that looked at the association between CYP2C8 polymorphisms and TZD toxicity are inconsistent and generally did not produce statistically significant results.¹⁰ Therefore, it can only be speculated that polymorphisms in TZD-metabolizing enzymes are associated with toxicity.

Meglitindes

Use of meglitinides and their mechanism of action

The antidiabetic drug class of meglitinides is comprised of nateglinide (Starlix^®) and rapeglinide (Prandin^®).⁵⁴ Similar to the SUs, the meglitinides act on the K_ATP channel, but at a distinct binding site, to induce insulin secretion as well as on calcium voltage-gated channels.¹⁰ They stimulate early-phase insulin production to reduce postprandial glucose and their effects are short-lived and leave late-phase and fasting insulin untouched.¹⁰ Possible reasons for interindividual variability in response to this class of antidiabetic medications can be contributed to polymorphisms in organic anion transporting polypeptide 1B1 (OATP1B1) gene (SLCO1B1) or the metabolizing enzyme, CYP2C9 (Table 4).^55,56

Table 4 Possible role and effects of variations on meglitinide therapy

	SLCO1B1	CYP2C9
Action	Possible role in hepatic uptake (transporter)	Metabolization
Variants	521TC, 521CC significantly higher plasma concentrations compared to homozygous TT carriers. Half-life longer in 521CC SNP compared to 521TT	Reduced metabolizing activity in some variations

Role of SLCO1B1 polymorphisms in meglitinide response

Zhang and colleagues studied three SNPs (521TT, 521TC, 521CC) of the SLCO1B1 gene and their associated effect on the pharmacokinetics of nateglinide.⁵⁵ The results showed carriers of heterozygous 521TC and homozygous 521CC had significantly higher plasma concentrations (83% and 76%, respectively) of nateglinide compared to the homozygous TT carriers.⁵⁵ The half-life of nateglinide was 78% longer (p = 0.036) in individuals carrying the 521CC SNP compared to the 521TT carriers.⁵⁵ Lastly, the AUC in the 521TC and 521CC carriers were significantly higher (82% and 108%, respectively) than that of the 521TT carriers.⁵⁵ These results suggest that SLCO1B1 polymorphisms have a significant impact on the plasma concentrations of nateglinide, which may be contributed to a reduced hepatic uptake.⁵⁵ However, the hepatic uptake, efficacy, and toxicity associated with the elevated plasma levels of nateglinide were not assessed, which makes it hard to apply this data clinically, but may act as a starting point for genetics-based nateglinide-prescribing practice.

Another study that assessed the effect of SLCO1B1 polymorphisms and the associated effects on nateglinide and rapeglinide pharmacokinetics produced conflicting results with respect to nateglinide.⁵⁶ Results showed that the 521T > C polymorphisms of SLCO1B1 does not alter the pharmacokinetics and efficacy (ability to reduce blood glucose) of nateglinide or its metabolite (M7), regardless of the individual’s genotype.⁵⁶ In contrast, the pharmacokinetics of repaglinide and its M2 and M4 metabolites were significantly affected by the 521T > C polymorphism.⁵⁶ Individuals expressing the 521CC genotype tended to have higher plasma concentrations of repaglinide and its metabolites, which also correlated with a greater ability to reduce blood glucose.⁵⁶ This indicates that the 521CC genotype may encode for a reduced function OATP1B1 transporter. Although the incidence of adverse events was not assessed in either study, this may be an important avenue for future studies to determine the clinical relevance this polymorphism.

Role of CYP2C9 in meglitinde metabolism and associated effect on response

Nateglinide is metabolized by CYP2C9 and repaglinide is metabolized by CYP2C8. Polymorphism of CYP2C9 were addressed in the same study that addressed the SU and was conducted using the same allelic forms (CYP2C9*1/*2/*3).¹⁰ Heterozygous forms, CYP2C9*1/*3 and CYP2C9*2/*3, and homozygous forms, CYP2C9*2/*2, exhibited intermediately reduced metabolizing activity compared to the wild-type (CYP2C9*1/*1).¹⁰ However, similar to the SUs, CYP2C9*3/*3 showed a significantly reduced (50%) metabolizing activity compared to the wild–type.¹⁰ A population pharmacokinetic/pharmacodynamics model was used to extrapolate associations between nateglinide and plasma insulin and glucose concentrations. Results showed that the risk for hypoglycemia increases for doses of 180 mg and higher and individuals with the CYP2C9*3/*3 genotype demonstrated an enhanced risk.¹⁰ Although the recommended upper limit is 120 mg, caution should be used in the homozygous CYP2C9*3 carriers.¹⁰ This information is inadequate to make dose-adjustment recommendations and carriers of these genotypes should be evaluated further to assess associated response and toxicity. On the other hand, a study investigating the CYP2C8 polymorphisms revealed highly contradictory results, which makes it hard to see a definitive effect with repaglinide.¹⁰ Other studies are warranted to examine the association between repaglinide and CYP2C8 polymorphisms.

Discussion

Pharmacogenomics has become a revolutionary approach for defining an individual’s clinical response to drug therapy by permitting researchers to investigate most, if not all, patient genomes for insight. However, there are still many ethical, social, economic, legislative, and research protocol issues that need to be sorted out before the promise of “personalized medicine” become a customary approach clinical practice. Furthermore, a great deal of education for the public and health care professionals in the area is necessary before gaining overall acceptance. Nevertheless, advances are rapidly moving this approach closer to clinical practice to be implicated in the management of many complex diseases (ie, cardiovascular disease, cancer, and DM).

Pharmacogenomic investigations of T2D have yielded variations in genomic loci containing common variants with translational consequences of altering disease susceptibility risk, clinical course, and response to standard therapy. Discoveries of gene polymorphisms in drug transporter, targets, effector proteins, and metabolizing enzymes have only served to augment discoveries of possible causal genes. Therefore, it is important to keep pushing the pharmacogenomic approach to advance our knowledge in disease state prevention and management to attain optimal health outcomes.

Executive summary

Inquiry regarding disease prevention, therapy and prognosis adopted an approach that interrogates the entire genome (genomics) for answers versus single genes (genetics), which has facilitated the discovery of gene variations in the form of SNPs in several genes associated with complex disease states.

One of the many areas in complex disease state management that the genomic approach has offered great insight is in T2D pharmacotherapy. Specifically, pharmacogenomics has provided a greater understanding of biological mechanisms that cause or contribute to interindividual variability in response of oral antidiabetic therapy. OCT1, OCT2, and MATE1 SNPs have shown to alter metformin’s pharmacokinetics and ability to inhibit glucose output from the liver and may help explain variability in response. SNPs in genes encoding the subunits of the K_ATP channel and hepatic metabolizing enzyme CYP2C9 have been shown to alter individual sensitivity to SU therapy, whereas polymorphisms in CYP2C9 has been shown to alter the pharmacokinetics.

Efficacy of thiazolidinedione therapy appears to be enhanced in individuals expressing genes resulting in higher levels of adiponectin and lipoprotein lipase. Additionally, SNPs in PPAR-γ have the potential to either enhance or reduce patient response to TZD therapy. Altered response to meglitinide therapy is largely due to SNPs that result in a lower functioning OATP1B1 transporter, resulting in reduced hepatic uptake or polymorphisms in metabolizing enzymes.

Ethical, political, legislative, economic, and technological issues have delayed widespread acceptance of pharmacogenomics into clinical practice. However, these issues are being resolved at a fairly rapid pace, which puts pharmacogenomics at the forefront of future medicine in clinical practice.

Disclosures

The authors report no conflicts of interest in this work.