Genome-Wide Meta-analysis Identifies Genetic Variants Associated With Glycemic Response to Sulfonylureas.

OBJECTIVE: Sulfonylureas, the first available drugs for the management of type 2 diabetes, remain widely prescribed today. However, there exists significant variability in glycemic response to treatment. We aimed to establish heritability of sulfonylurea response and identify genetic variants and interacting treatments associated with HbA1c reduction. RESEARCH DESIGN AND METHODS: As an initiative of the Metformin Genetics Plus Consortium (MetGen Plus) and the DIabetes REsearCh on patient straTification (DIRECT) consortium, 5,485 White Europeans with type 2 diabetes treated with sulfonylureas were recruited from six referral centers in Europe and North America. We first estimated heritability using the generalized restricted maximum likelihood approach and then undertook genome-wide association studies of glycemic response to sulfonylureas measured as HbA1c reduction after 12 months of therapy followed by meta-analysis. These results were supported by acute glipizide challenge in humans who were naïve to type 2 diabetes medications, cis expression quantitative trait loci (eQTL), and functional validation in cellular models. Finally, we examined for possible drug-drug-gene interactions. RESULTS: After establishing that sulfonylurea response is heritable (mean ± SEM 37 ± 11%), we identified two independent loci near the GXYLT1 and SLCO1B1 genes associated with HbA1c reduction at a genome-wide scale (P < 5 × 10-8). The C allele at rs1234032, near GXYLT1, was associated with 0.14% (1.5 mmol/mol), P = 2.39 × 10-8), lower reduction in HbA1c. Similarly, the C allele was associated with higher glucose trough levels (ß = 1.61, P = 0.005) in healthy volunteers in the SUGAR-MGH given glipizide (N = 857). In 3,029 human whole blood samples, the C allele is a cis eQTL for increased expression of GXYLT1 (ß = 0.21, P = 2.04 × 10-58). The C allele of rs10770791, in an intronic region of SLCO1B1, was associated with 0.11% (1.2 mmol/mol) greater reduction in HbA1c (P = 4.80 × 10-8). In 1,183 human liver samples, the C allele at rs10770791 is a cis eQTL for reduced SLCO1B1 expression (P = 1.61 × 10-7), which, together with functional studies in cells expressing SLCO1B1, supports a key role for hepatic SLCO1B1 (encoding OATP1B1) in regulation of sulfonylurea transport. Further, a significant interaction between statin use and SLCO1B1 genotype was observed (P = 0.001). In statin nonusers, C allele homozygotes at rs10770791 had a large absolute reduction in HbA1c (0.48 ± 0.12% [5.2 ± 1.26 mmol/mol]), equivalent to that associated with initiation of a dipeptidyl peptidase 4 inhibitor. CONCLUSIONS: We have identified clinically important genetic effects at genome-wide levels of significance, and important drug-drug-gene interactions, which include commonly prescribed statins. With increasing availability of genetic data embedded in clinical records these findings will be important in prescribing glucose-lowering drugs.

Variation in the Plasma Membrane Monoamine Transporter (PMAT) (Encoded by SLC29A4) and Organic Cation Transporter 1 (OCT1) (Encoded by SLC22A1) and Gastrointestinal Intolerance to Metformin in Type 2 Diabetes: An IMI DIRECT Study.

OBJECTIVE: Gastrointestinal adverse effects occur in 20-30% of patients with metformin-treated type 2 diabetes, leading to premature discontinuation in 5-10% of the cases. Gastrointestinal intolerance may reflect localized high concentrations of metformin in the gut. We hypothesized that reduced transport of metformin via the plasma membrane monoamine transporter (PMAT) and organic cation transporter 1 (OCT1) could increase the risk of severe gastrointestinal adverse effects. RESEARCH DESIGN AND METHODS: The study included 286 severe metformin-intolerant and 1,128 metformin-tolerant individuals from the IMI DIRECT (Innovative Medicines Initiative: DIabetes REsearCh on patient straTification) consortium. We assessed the association of patient characteristics, concomitant medication, and the burden of mutations in the SLC29A4 and SLC22A1 genes on odds of intolerance. RESULTS: Women (P < 0.001) and older people (P < 0.001) were more likely to develop metformin intolerance. Concomitant use of transporter-inhibiting drugs increased the odds of intolerance (odds ratio [OR] 1.72, P < 0.001). In an adjusted logistic regression model, the G allele at rs3889348 (SLC29A4) was associated with gastrointestinal intolerance (OR 1.34, P = 0.005). rs3889348 is the top cis-expression quantitative trait locus for SLC29A4 in gut tissue where carriers of the G allele had reduced expression. Homozygous carriers of the G allele treated with transporter-inhibiting drugs had more than three times higher odds of intolerance compared with carriers of no G allele and not treated with inhibiting drugs (OR 3.23, P < 0.001). Use of a genetic risk score derived from rs3889348 and SLC22A1 variants found that the odds of intolerance were more than twice as high in individuals who carry three or more risk alleles compared with those carrying none (OR 2.15, P = 0.01). CONCLUSIONS: These results suggest that intestinal metformin transporters and concomitant medications play an important role in the gastrointestinal adverse effects of metformin.

Integrative network analysis highlights biological processes underlying GLP-1 stimulated insulin secretion: A DIRECT study.

Glucagon-like peptide 1 (GLP-1) stimulated insulin secretion has a considerable heritable component as estimated from twin studies, yet few genetic variants influencing this phenotype have been identified. We performed the first genome-wide association study (GWAS) of GLP-1 stimulated insulin secretion in non-diabetic individuals from the Netherlands Twin register (n = 126). This GWAS was enhanced using a tissue-specific protein-protein interaction network approach. We identified a beta-cell protein-protein interaction module that was significantly enriched for low gene scores based on the GWAS P-values and found support at the network level in an independent cohort from Tübingen, Germany (n = 100). Additionally, a polygenic risk score based on SNPs prioritized from the network was associated (P < 0.05) with glucose-stimulated insulin secretion phenotypes in up to 5,318 individuals in MAGIC cohorts. The network contains both known and novel genes in the context of insulin secretion and is enriched for members of the focal adhesion, extracellular-matrix receptor interaction, actin cytoskeleton regulation, Rap1 and PI3K-Akt signaling pathways. Adipose tissue is, like the beta-cell, one of the target tissues of GLP-1 and we thus hypothesized that similar networks might be functional in both tissues. In order to verify peripheral effects of GLP-1 stimulation, we compared the transcriptome profiling of ob/ob mice treated with liraglutide, a clinically used GLP-1 receptor agonist, versus baseline controls. Some of the upstream regulators of differentially expressed genes in the white adipose tissue of ob/ob mice were also detected in the human beta-cell network of genes associated with GLP-1 stimulated insulin secretion. The findings provide biological insight into the mechanisms through which the effects of GLP-1 may be modulated and highlight a potential role of the beta-cell expressed genes RYR2, GDI2, KIAA0232, COL4A1 and COL4A2 in GLP-1 stimulated insulin secretion.

Metabolite ratios as potential biomarkers for type 2 diabetes: a DIRECT study.

AIMS/HYPOTHESIS: Circulating metabolites have been shown to reflect metabolic changes during the development of type 2 diabetes. In this study we examined the association of metabolite levels and pairwise metabolite ratios with insulin responses after glucose, glucagon-like peptide-1 (GLP-1) and arginine stimulation. We then investigated if the identified metabolite ratios were associated with measures of OGTT-derived beta cell function and with prevalent and incident type 2 diabetes. METHODS: We measured the levels of 188 metabolites in plasma samples from 130 healthy members of twin families (from the Netherlands Twin Register) at five time points during a modified 3 h hyperglycaemic clamp with glucose, GLP-1 and arginine stimulation. We validated our results in cohorts with OGTT data (n = 340) and epidemiological case-control studies of prevalent (n = 4925) and incident (n = 4277) diabetes. The data were analysed using regression models with adjustment for potential confounders. RESULTS: There were dynamic changes in metabolite levels in response to the different secretagogues. Furthermore, several fasting pairwise metabolite ratios were associated with one or multiple clamp-derived measures of insulin secretion (all p < 9.2 × 10-7). These associations were significantly stronger compared with the individual metabolite components. One of the ratios, valine to phosphatidylcholine acyl-alkyl C32:2 (PC ae C32:2), in addition showed a directionally consistent positive association with OGTT-derived measures of insulin secretion and resistance (p ≤ 5.4 × 10-3) and prevalent type 2 diabetes (ORVal_PC ae C32:2 2.64 [ß 0.97 ± 0.09], p = 1.0 × 10-27). Furthermore, Val_PC ae C32:2 predicted incident diabetes independent of established risk factors in two epidemiological cohort studies (HRVal_PC ae C32:2 1.57 [ß 0.45 ± 0.06]; p = 1.3 × 10-15), leading to modest improvements in the receiver operating characteristics when added to a model containing a set of established risk factors in both cohorts (increases from 0.780 to 0.801 and from 0.862 to 0.865 respectively, when added to the model containing traditional risk factors + glucose). CONCLUSIONS/INTERPRETATION: In this study we have shown that the Val_PC ae C32:2 metabolite ratio is associated with an increased risk of type 2 diabetes and measures of insulin secretion and resistance. The observed effects were stronger than that of the individual metabolites and independent of known risk factors.

HbA1c is associated with altered expression in blood of cell cycle- and immune response-related genes.

AIMS/HYPOTHESIS: Individuals with type 2 diabetes are heterogeneous in their glycaemic control as tracked by blood HbA1c levels. Here, we investigated the extent to which gene expression levels in blood reflect current and future HbA1c levels. METHODS: HbA1c levels at baseline and 1 and 2 year follow-up were compared with gene expression levels in 391 individuals with type 2 diabetes from the Hoorn Diabetes Care System Cohort (15,564 genes, RNA sequencing). The functions of associated baseline genes were investigated further using pathway enrichment analysis. Using publicly available data, we investigated whether the genes identified are also associated with HbA1c in the target tissues, muscle and pancreas. RESULTS: At baseline, 220 genes (1.4%) were associated with baseline HbA1c. Identified genes were enriched for cell cycle and complement system activation pathways. The association of 15 genes extended to the target tissues, muscle (n = 113) and pancreatic islets (n = 115). At follow-up, expression of 25 genes (0.16%) associated with 1 year HbA1c and nine genes (0.06%) with 2 year HbA1c. Five genes overlapped across all time points, and 18 additional genes between baseline and 1 year follow-up. After adjustment for baseline HbA1c, the number of significant genes at 1 and 2 years markedly decreased, suggesting that gene expression levels in whole blood reflect the current glycaemic state and but not necessarily the future glycaemic state. CONCLUSIONS/INTERPRETATION: HbA1c levels in individuals with type 2 diabetes are associated with expression levels of genes that link to the cell cycle and complement system activation.

The Hoorn Diabetes Care System (DCS) cohort. A prospective cohort of persons with type 2 diabetes treated in primary care in the Netherlands.

PURPOSE: People with type 2 diabetes (T2D) have a doubled morbidity and mortality risk compared with persons with normal glucose tolerance. Despite treatment, clinical targets for cardiovascular risk factors are not achieved. The Hoorn Diabetes Care System cohort (DCS) is a prospective cohort representing a comprehensive dataset on the natural course of T2D, with repeated clinical measures and outcomes. In this paper, we describe the design of the DCS cohort. PARTICIPANTS: The DCS consists of persons with T2D in primary care from the West-Friesland region of the Netherlands. Enrolment in the cohort started in 1998 and this prospective dynamic cohort currently holds 12 673 persons with T2D. FINDINGS TO DATE: Clinical measures are collected annually, with a high internal validity due to the centrally organised standardised examinations. Microvascular complications are assessed by measuring kidney function, and screening feet and eyes. Information on cardiovascular disease is obtained by 1) self-report, 2) electrocardiography and 3) electronic patient records. In subgroups of the cohort, biobanking and additional measurements were performed to obtain information on, for example, lifestyle, depression and genomics. Finally, the DCS cohort is linked to national cancer and all-cause mortality registers. A selection of published findings from the DCS includes identification of subgroups with distinct development of haemoglobin A1c, blood pressure and retinopathy, and their predictors; validation of a prediction model for personalised retinopathy screening; the assessment of the role of genetics in development and treatment of T2D, providing options for personalised medicine. FUTURE PLANS: We will continue with the inclusion of persons with newly diagnosed T2D, follow-up of persons in the cohort and linkage to morbidity and mortality registries. Currently, we are involved in (inter)national projects on, among others, biomarkers and prediction models for T2D and complications and we are interested in collaborations with external researchers. TRIAL REGISTRATION: ISRCTN26257579.

Variation in the glucose transporter gene SLC2A2 is associated with glycemic response to metformin.

Metformin is the first-line antidiabetic drug with over 100 million users worldwide, yet its mechanism of action remains unclear. Here the Metformin Genetics (MetGen) Consortium reports a three-stage genome-wide association study (GWAS), consisting of 13,123 participants of different ancestries. The C allele of rs8192675 in the intron of SLC2A2, which encodes the facilitated glucose transporter GLUT2, was associated with a 0.17% (P = 6.6 × 10(-14)) greater metformin-induced reduction in hemoglobin A1c (HbA1c) in 10,577 participants of European ancestry. rs8192675 was the top cis expression quantitative trait locus (cis-eQTL) for SLC2A2 in 1,226 human liver samples, suggesting a key role for hepatic GLUT2 in regulation of metformin action. Among obese individuals, C-allele homozygotes at rs8192675 had a 0.33% (3.6 mmol/mol) greater absolute HbA1c reduction than T-allele homozygotes. This was about half the effect seen with the addition of a DPP-4 inhibitor, and equated to a dose difference of 550 mg of metformin, suggesting rs8192675 as a potential biomarker for stratified medicine.

The CTRB1/2 locus affects diabetes susceptibility and treatment via the incretin pathway.

The incretin hormone glucagon-like peptide 1 (GLP-1) promotes glucose homeostasis and enhances ß-cell function. GLP-1 receptor agonists (GLP-1 RAs) and dipeptidyl peptidase-4 (DPP-4) inhibitors, which inhibit the physiological inactivation of endogenous GLP-1, are used for the treatment of type 2 diabetes. Using the Metabochip, we identified three novel genetic loci with large effects (30-40%) on GLP-1-stimulated insulin secretion during hyperglycemic clamps in nondiabetic Caucasian individuals (TMEM114; CHST3 and CTRB1/2; n = 232; all P ≤ 8.8 × 10(-7)). rs7202877 near CTRB1/2, a known diabetes risk locus, also associated with an absolute 0.51 ± 0.16% (5.6 ± 1.7 mmol/mol) lower A1C response to DPP-4 inhibitor treatment in G-allele carriers, but there was no effect on GLP-1 RA treatment in type 2 diabetic patients (n = 527). Furthermore, in pancreatic tissue, we show that rs7202877 acts as expression quantitative trait locus for CTRB1 and CTRB2, encoding chymotrypsinogen, and increases fecal chymotrypsin activity in healthy carriers. Chymotrypsin is one of the most abundant digestive enzymes in the gut where it cleaves food proteins into smaller peptide fragments. Our data identify chymotrypsin in the regulation of the incretin pathway, development of diabetes, and response to DPP-4 inhibitor treatment.

Human mineralocorticoid receptor (MR) gene haplotypes modulate MR expression and transactivation: implication for the stress response.

Stress causes activation of the hypothalamic-pituitary-adrenal (HPA) axis, resulting in secretion of corticosteroids which facilitate behavioural adaptation. These effects exerted by corticosteroids are mediated by two brain corticosteroid receptor types, the mineralocorticoid receptor (MR), with a high affinity already occupied under basal conditions and the glucocorticoid receptor (GR), with a low affinity only activated during stress. Here, we studied MR gene haplotypes constituted by the two single nucleotide polymorphisms MR-2G/C (rs2070951) and MRI180V (rs5522). The haplotypes showed differences in cortisol-induced gene transcription and protein expression while the structural variant MRI180V did not affect ligand binding. Moreover, in a well characterized cohort of 166 school teachers these haplotypes have been associated with perceived chronic stress (Trier Inventory for the Assessment of Chronic Stress, TICS) and, in a subgroup of 47 subjects, with ACTH, cortisol and heart rate responses to acute psychosocial stress (Trier Social Stress Test, TSST). MR haplotypes were significantly associated with the TICS scales "excessive demands at work" and "social overload". Subjects homozygous for haplotype MR-2C/MRI180, which in vitro showed highest expression and transactivational activity, displayed the highest salivary cortisol (p<0.001), plasma cortisol (p=0.010), plasma ACTH (p=0.003) and heart rate (p=0.018) responses. It is concluded that the investigated MR haplotypes modulate cortisol-induced gene transcription in vitro. Moreover, these haplotypes may contribute to individual differences in perceived chronic stress as well as neuroendocrine and cardiovascular stress responses.

Mineralocorticoid receptor gene variants as determinants of HPA axis regulation and behavior.

We tested if common mineralocorticoid receptor (MR) gene variants contribute to the variability in neuroendocrine control and behavioral reactivity as observed in humans. For that purpose we screened for genetic variability and tested functionality of the identified human MR gene variants in vitro. Four haplotypes were tested for transactivational capacity in vitro and showed profound significant differences when stimulated with cortisol. The MR gene variants were associated with basal levels of cortisol, cortisol levels after dexamethasone administration and with stress-induced hypothalamic-pituitary-adrenal axis and autonomic reactivity. In an elderly cohort, one of the functional MR gene variants, MR-I180V, associated with higher feelings of depression. Moreover, we found an association with neuroticism in a second cohort consisting of depressed patients. In conclusion, we report here new findings on common functional human MR gene variants which reveal a hitherto unknown role of these variants in neuroticism conferring vulnerability to stress-related mental disorders, such as depression and posttraumatic stress syndrome.

The functional c.-2G>C variant of the mineralocorticoid receptor modulates blood pressure, renin, and aldosterone levels.

The mineralocorticoid receptor (MR) is essential in the regulation of volemia and blood pressure. Rare mutations in the MR gene cause type 1 pseudohypoaldosteronism and hypertension. In this study we characterized the common MR polymorphism c.-2G>C (rs2070951) in vitro and tested its influence on parameters related to blood pressure regulation and the renin-angiotensin system. In vitro studies showed that the G allele was associated with decreased MR protein levels and reduced transcriptional activation compared with the C allele. Association studies were performed with several outcome variables in 3 independent cohorts: a mild hypertensive group subjected to a salt-sensitivity test, a healthy normotensive group included in a crossover study to receive both a high and low Na/K diet, and a large cohort (The Netherlands Study of Depression and Anxiety), in which blood pressure was measured. Subjects with the GG genotype had significantly higher plasma renin levels both in the mild hypertensive group and in normal volunteers compared with homozygous C carriers. The GG genotype was also correlated with higher plasma aldosterone levels in healthy subjects. In both the mild hypertensive group and The Netherlands Study of Depression and Anxiety cohort the genotype GG was associated with higher systolic blood pressure in males. In conclusion, the G allele of the common functional genetic polymorphism c.-2G>C in the MR gene associates with increased activation of the renin-angiotensin-aldosterone axis and with increased blood pressure, probably related to decreased MR expression.

Functional mineralocorticoid receptor (MR) gene variation influences the cortisol awakening response after dexamethasone.

Stress causes activation of the hypothalamic-pituitary-adrenal (HPA) axis and results in the secretion of corticosteroids, which facilitate behavioral adaptation and promote the termination of the stress response. These actions exerted by cortisol are mediated by two brain corticosteroid receptor types: the high affinity mineralocorticoid (MR) and the low affinity glucocorticoid receptor (GR). Dexamethasone is a potent GR agonist with affinity to MR. Administration of dexamethasone in the evening results in a significant suppression of the morning cortisol awakening response (CAR). Here we tested the involvement of MR variants in this effect of dexamethasone in 218 young healthy subjects (125 females, all using oral contraceptives). For this purpose we determined two single nucleotide polymorphisms (SNPs) in the MR gene, the previously described MRI180V (rs5522) and the MR-2G/C (rs2070951), which both affect in vitro the transactivational capacity of the MR in response to either cortisol or dexamethasone. Administration of a low dose dexamethasone (0.25mg) at 2300h resulted in a significant suppression of the cortisol awakening response (CAR). Both SNPs modulated the suppression of the CAR after dexamethasone significantly and in a sex specific manner. Suppression of the CAR was highest in the female MR-2G/C GG subjects while in male GG subjects the dexamethasone suppression of the CAR was attenuated compared to the MR-2G/C GC and CC groups. For the MRI180V, male AA subjects showed after dexamethasone a higher CAR than AG subjects while this effect was not observed in females. The SNPs had no significant influence on the CAR without prior dexamethasone treatment. The association of the CAR with functional MR gene variants only in dexamethasone treated subjects suggests the involvement of MR in dexamethasone induced suppression of morning cortisol.

Corticosteroid receptor-gene variants: modulators of the stress-response and implications for mental health.

The stress-response, including autonomic and hypothalamic-pituitary-adrenal (HPA) axis reactivity, is essential for maintaining homeostasis during a challenge. Brain mineralocorticoid receptors and glucocorticoid receptors operate in balance to coordinate the stress-response. Genetic variants in both the human mineralocorticoid and glucocorticoid receptor-genes have been functionally characterized. In vitro effects of these genetic variants on transactivation and mRNA stability have been described. In vivo, two mineralocorticoid receptor-gene SNPs (-2 G/C (allele frequency: 50%), MR I180V (11%)) and four glucocorticoid receptor-gene SNPs (ER22/23EK (3%), N363S (4%), BclI (37%), A3669G (15%)) are associated with changes in hypothalamic-pituitary-adrenal (HPA) axis reactivity. Importantly, the two mineralocorticoid receptor-gene variants (but none of the glucocorticoid receptor-gene variants) also associate with changes in autonomic output as measured as increased heart beat following a psychosocial stress (TSST). Moreover, several of these mineralocorticorticoid receptor- and glucocorticoid receptor variants have been found associated with stress-related disorders, including depression. These data indicate that dysregulation of mineralocorticoid- and glucocorticoid receptor are causative in the pathogenesis of depression. Moreover, these mineralocorticoid- and glucocorticoid receptor-gene variants constitute part of the genetic make up that determines individual stress-responsiveness inducing vulnerability to disease. Furthermore, mineralocorticoid- and glucocorticoid receptors are drug targets, thereby aiming at the underlying mechanisms of stress-related disorders.