The association of cardiometabolic, diet and lifestyle parameters with plasma glucagon-like peptide-1: An IMI DIRECT study.

CONTEXT: The role of glucagon-like peptide-1(GLP-1) in Type 2 diabetes (T2D) and obesity is not fully understood. OBJECTIVE: We investigate the association of cardiometabolic, diet and lifestyle parameters on fasting and postprandial GLP-1 in people at risk of, or living with, T2D. METHOD: We analysed cross-sectional data from the two Innovative Medicines Initiative (IMI) Diabetes Research on Patient Stratification (DIRECT) cohorts, cohort 1(n=2127) individuals at risk of diabetes; cohort 2 (n=789) individuals with new-onset of T2D. RESULTS: Our multiple regression analysis reveals that fasting total GLP-1 is associated with an insulin resistant phenotype and observe a strong independent relationship with male sex, increased adiposity and liver fat particularly in the prediabetes population. In contrast, we showed that incremental GLP-1 decreases with worsening glycaemia, higher adiposity, liver fat, male sex and reduced insulin sensitivity in the prediabetes cohort. Higher fasting total GLP-1 was associated with a low intake of wholegrain, fruit and vegetables inpeople with prediabetes, and with a high intake of red meat and alcohol in people with diabetes. CONCLUSION: These studies provide novel insights into the association between fasting and incremental GLP-1, metabolic traits of diabetes and obesity, and dietary intake and raise intriguing questions regarding the relevance of fasting GLP-1 in the pathophysiology T2D.

Treatment effect heterogeneity following type 2 diabetes treatment with GLP1-receptor agonists and SGLT2-inhibitors: a systematic review.

BACKGROUND: A precision medicine approach in type 2 diabetes requires the identification of clinical and biological features that are reproducibly associated with differences in clinical outcomes with specific anti-hyperglycaemic therapies. Robust evidence of such treatment effect heterogeneity could support more individualized clinical decisions on optimal type 2 diabetes therapy. METHODS: We performed a pre-registered systematic review of meta-analysis studies, randomized control trials, and observational studies evaluating clinical and biological features associated with heterogenous treatment effects for SGLT2-inhibitor and GLP1-receptor agonist therapies, considering glycaemic, cardiovascular, and renal outcomes. After screening 5,686 studies, we included 101 studies of SGLT2-inhibitors and 75 studies of GLP1-receptor agonists in the final systematic review. RESULTS: Here we show that the majority of included papers have methodological limitations precluding robust assessment of treatment effect heterogeneity. For SGLT2-inhibitors, multiple observational studies suggest lower renal function as a predictor of lesser glycaemic response, while markers of reduced insulin secretion predict lesser glycaemic response with GLP1-receptor agonists. For both therapies, multiple post-hoc analyses of randomized control trials (including trial meta-analysis) identify minimal clinically relevant treatment effect heterogeneity for cardiovascular and renal outcomes. CONCLUSIONS: Current evidence on treatment effect heterogeneity for SGLT2-inhibitor and GLP1-receptor agonist therapies is limited, likely reflecting the methodological limitations of published studies. Robust and appropriately powered studies are required to understand type 2 diabetes treatment effect heterogeneity and evaluate the potential for precision medicine to inform future clinical care.

This study reviews the available evidence on which patient features (such as age, sex, and blood test results) are associated with different outcomes for two recently introduced type 2 diabetes medications: SGLT2-inhibitors and GLP1-receptor agonists. Understanding what individual characteristics are associated with different response patterns may help clinical providers and people living with diabetes make more informed decisions about which type 2 diabetes treatments will work best for an individual. We focus on three outcomes: blood glucose levels (raised blood glucose is the primary symptom of diabetes and a primary aim of diabetes treatment is to lower this), heart disease, and kidney disease. We identified some potential factors that reduce effects on blood glucose levels, including poorer kidney function for SGLT2-inhibitors and lower production of the glucose-lowering hormone insulin for GLP1-receptor agonists. We did not identify clear factors that alter heart and kidney disease outcomes for either medication. Most of the studies had limitations, meaning more research is needed to fully understand the factors that influence treatment outcomes in type 2 diabetes.

A gene risk score using missense variants in SLCO1B1 is associated with earlier onset statin intolerance.

BACKGROUND AND AIMS: The efficacy of statin therapy is hindered by intolerance to the therapy, leading to discontinuation. Variants in SLCO1B1, which encodes the hepatic transporter OATB1B1, influence statin pharmacokinetics, resulting in altered plasma concentrations of the drug and its metabolites. Current pharmacogenetic guidelines require sequencing of the SLCO1B1 gene, which is more expensive and less accessible than genotyping. In this study, we aimed to develop an easy, clinically implementable functional gene risk score (GRS) of common variants in SLCO1B1 to identify patients at risk of statin intolerance. METHODS AND RESULTS: A GRS was developed from four common variants in SLCO1B1. In statin users from Tayside, Scotland, UK, those with a high-risk GRS had increased odds across three phenotypes of statin intolerance [general statin intolerance (GSI): ORGSI 2.42; 95% confidence interval (CI): 1.29-4.31, P = 0.003; statin-related myopathy: ORSRM 2.51; 95% CI: 1.28-4.53, P = 0.004; statin-related suspected rhabdomyolysis: ORSRSR 2.85; 95% CI: 1.03-6.65, P = 0.02]. In contrast, using the Val174Ala genotype alone or the recommended OATP1B1 functional phenotypes produced weaker and less reliable results. A meta-analysis with results from adjudicated cases of statin-induced myopathy in the PREDICTION-ADR Consortium confirmed these findings (ORVal174Ala 1.99; 95% CI: 1.01-3.95, P = 0.048; ORGRS 1.76; 95% CI: 1.16-2.69, P = 0.008). For those requiring high-dose statin therapy, the high-risk GRS was more consistently associated with the time to onset of statin intolerance amongst the three phenotypes compared with Val174Ala (GSI: HRVal174Ala 2.49; 95% CI: 1.09-5.68, P = 0.03; HRGRS 2.44; 95% CI: 1.46-4.08, P < 0.001). Finally, sequence kernel association testing confirmed that rare variants in SLCO1B1 are associated with the risk of intolerance (P = 0.02). CONCLUSION: We provide evidence that a GRS based on four common SLCO1B1 variants provides an easily implemented genetic tool that is more reliable than the current recommended practice in estimating the risk and predicting early-onset statin intolerance.

Precision medicine in type 2 diabetes: A systematic review of treatment effect heterogeneity for GLP1-receptor agonists and SGLT2-inhibitors.

Background: A precision medicine approach in type 2 diabetes requires identification of clinical and biological features that are reproducibly associated with differences in clinical outcomes with specific anti-hyperglycaemic therapies. Robust evidence of such treatment effect heterogeneity could support more individualized clinical decisions on optimal type 2 diabetes therapy. Methods: We performed a pre-registered systematic review of meta-analysis studies, randomized control trials, and observational studies evaluating clinical and biological features associated with heterogenous treatment effects for SGLT2-inhibitor and GLP1-receptor agonist therapies, considering glycaemic, cardiovascular, and renal outcomes. Results: After screening 5,686 studies, we included 101 studies of SGLT2-inhibitors and 75 studies of GLP1-receptor agonists in the final systematic review. The majority of papers had methodological limitations precluding robust assessment of treatment effect heterogeneity. For glycaemic outcomes, most cohorts were observational, with multiple analyses identifying lower renal function as a predictor of lesser glycaemic response with SGLT2-inhibitors and markers of reduced insulin secretion as predictors of lesser response with GLP1-receptor agonists. For cardiovascular and renal outcomes, the majority of included studies were post-hoc analyses of randomized control trials (including meta-analysis studies) which identified limited clinically relevant treatment effect heterogeneity. Conclusions: Current evidence on treatment effect heterogeneity for SGLT2-inhibitor and GLP1-receptor agonist therapies is limited, likely reflecting the methodological limitations of published studies. Robust and appropriately powered studies are required to understand type 2 diabetes treatment effect heterogeneity and evaluate the potential for precision medicine to inform future clinical care. Plain language summary: This review identifies research that helps understand which clinical and biological factors that are associated with different outcomes for specific type 2 diabetes treatments. This information could help clinical providers and patients make better informed personalized decisions about type 2 diabetes treatments. We focused on two common type 2 diabetes treatments: SGLT2-inhibitors and GLP1-receptor agonists, and three outcomes: blood glucose control, heart disease, and kidney disease. We identified some potential factors that are likely to lessen blood glucose control including lower kidney function for SGLT2-inhibitors and lower insulin secretion for GLP1-receptor agonists. We did not identify clear factors that alter heart and renal disease outcomes for either treatment. Most of the studies had limitations, meaning more research is needed to fully understand the factors that influence treatment outcomes in type 2 diabetes.

Type 2 Diabetes risk alleles in Peptidyl-glycine Alpha-amidating Monooxygenase influence GLP-1 levels and response to GLP-1 Receptor Agonists.

Patients with type 2 diabetes vary in their response to currently available therapeutic agents (including GLP-1 receptor agonists) leading to suboptimal glycemic control and increased risk of complications. We show that human carriers of hypomorphic T2D-risk alleles in the gene encoding peptidyl-glycine alpha-amidating monooxygenase (PAM), as well as Pam-knockout mice, display increased resistance to GLP-1 in vivo. Pam inactivation in mice leads to reduced gastric GLP-1R expression and faster gastric emptying: this persists during GLP-1R agonist treatment and is rescued when GLP-1R activity is antagonized, indicating resistance to GLP-1's gastric slowing properties. Meta-analysis of human data from studies examining GLP-1R agonist response (including RCTs) reveals a relative loss of 44% and 20% of glucose lowering (measured by glycated hemoglobin) in individuals with hypomorphic PAM alleles p.S539W and p.D536G treated with GLP-1R agonist. Genetic variation in PAM has effects on incretin signaling that alters response to medication used commonly for treatment of T2D.

Weight variability and cardiovascular outcomes: a systematic review and meta-analysis.

The association between body weight variability and the risk of cardiovascular disease (CVD) has been investigated previously with mixed findings. However, there has been no extensive study which systematically evaluates the current evidence. Furthermore, the impact of ethnicity and type 2 diabetes on this phenomena has not yet been investigated. Therefore, the aim of this study was to comprehensively evaluate the effect of weight variability on risk of CVD (any cardiovascular (CV) event, composite CV outcome, CV death, Stroke, Myocardial Infarction) and the influence of ethnicity and type 2 diabetes status on the observed association. A systematic review and meta-analysis was performed according to the meta-analyses of observational studies in epidemiology (MOOSE) guidelines. The electronic databases PubMed, Web of Science, and the Cochrane Library were searched for studies that investigated the relationship between body weight or BMI variability and CV diseases using Medical Subject Headings (MeSH) terms and keywords. The relative risks (RRs) for the outcomes were collected from studies, pooled, and analysed using a random-effects model to estimate the overall relative risk. Of 5645 articles screened, 23 studies with a total population of 15,382,537 fulfilled the prespecified criteria and were included. Individuals in the highest strata of body weight variability were found to have significantly increased risk of any CV event (RR = 1.27; 95% Confidence Interval (CI) 1.17-1.38; P < 0.0001; I2 = 97.28%), cardiovascular death (RR = 1.29; 95% CI 1.03-1.60; P < 0.0001; I2 = 55.16%), myocardial infarction (RR = 1.32; 95% CI 1.09-1.59; P = 0.0037; I2 = 97.14%), stroke (RR = 1.21; 95% CI 1.19-1.24; P < 0.0001; I2 = 0.06%), and compound CVD outcomes (RR = 1.36; 95% CI 1.08-1.73; P = 0.01; I2 = 92.41%). Similar RRs were observed regarding BMI variability and per unit standard deviation (SD) increase in body weight variability. Comparable effects were seen in people with and without diabetes, in White Europeans and Asians. In conclusion, body weight variability is associated with increased risk of CV diseases regardless of ethnicity or diabetes status. Future research is needed to prove a causative link between weight variability and CVD risk, as appropriate interventions to maintain stable weight could positively influence CVD.

Precision Medicine in Diabetes.

Tailoring treatment or management to groups of individuals based on specific clinical, molecular, and genomic features is the concept of precision medicine. Diabetes is highly heterogenous with respect to clinical manifestations, disease progression, development of complications, and drug response. The current practice for drug treatment is largely based on evidence from clinical trials that report average effects. However, around half of patients with type 2 diabetes do not achieve glycaemic targets despite having a high level of adherence and there are substantial differences in the incidence of adverse outcomes. Therefore, there is a need to identify predictive markers that can inform differential drug responses at the point of prescribing. Recent advances in molecular genetics and increased availability of real-world and randomised trial data have started to increase our understanding of disease heterogeneity and its impact on potential treatments for specific groups. Leveraging information from simple clinical features (age, sex, BMI, ethnicity, and co-prescribed medications) and genomic markers has a potential to identify sub-groups who are likely to benefit from a given drug with minimal adverse effects. In this chapter, we will discuss the state of current evidence in the discovery of clinical and genetic markers that have the potential to optimise drug treatment in type 2 diabetes.

Identification of Genetic Variation Influencing Metformin Response in a Multiancestry Genome-Wide Association Study in the Diabetes Prevention Program (DPP).

Genome-wide significant loci for metformin response in type 2 diabetes reported elsewhere have not been replicated in the Diabetes Prevention Program (DPP). To assess pharmacogenetic interactions in prediabetes, we conducted a genome-wide association study (GWAS) in the DPP. Cox proportional hazards models tested associations with diabetes incidence in the metformin (MET; n = 876) and placebo (PBO; n = 887) arms. Multiple linear regression assessed association with 1-year change in metformin-related quantitative traits, adjusted for baseline trait, age, sex, and 10 ancestry principal components. We tested for gene-by-treatment interaction. No significant associations emerged for diabetes incidence. We identified four genome-wide significant variants after correcting for correlated traits (P < 9 × 10-9). In the MET arm, rs144322333 near ENOSF1 (minor allele frequency [MAF]AFR = 0.07; MAFEUR = 0.002) was associated with an increase in percentage of glycated hemoglobin (per minor allele, ß = 0.39 [95% CI 0.28, 0.50]; P = 2.8 × 10-12). rs145591055 near OMSR (MAF = 0.10 in American Indians) was associated with weight loss (kilograms) (per G allele, ß = -7.55 [95% CI -9.88, -5.22]; P = 3.2 × 10-10) in the MET arm. Neither variant was significant in PBO; gene-by-treatment interaction was significant for both variants [P(G×T) < 1.0 × 10-4]. Replication in individuals with diabetes did not yield significant findings. A GWAS for metformin response in prediabetes revealed novel ethnic-specific associations that require further investigation but may have implications for tailored therapy.

Pharmacogenomics of GLP-1 receptor agonists: a genome-wide analysis of observational data and large randomised controlled trials.

BACKGROUND: In the treatment of type 2 diabetes, GLP-1 receptor agonists lower blood glucose concentrations, body weight, and have cardiovascular benefits. The efficacy and side effects of GLP-1 receptor agonists vary between people. Human pharmacogenomic studies of this inter-individual variation can provide both biological insight into drug action and provide biomarkers to inform clinical decision making. We therefore aimed to identify genetic variants associated with glycaemic response to GLP-1 receptor agonist treatment. METHODS: In this genome-wide analysis we included adults (aged ≥18 years) with type 2 diabetes treated with GLP-1 receptor agonists with baseline HbA1c of 7% or more (53 mmol/mol) from four prospective observational cohorts (DIRECT, PRIBA, PROMASTER, and GoDARTS) and two randomised clinical trials (HARMONY phase 3 and AWARD). The primary endpoint was HbA1c reduction at 6 months after starting GLP-1 receptor agonists. We evaluated variants in GLP1R, then did a genome-wide association study and gene-based burden tests. FINDINGS: 4571 adults were included in our analysis, of these, 3339 (73%) were White European, 449 (10%) Hispanic, 312 (7%) American Indian or Alaskan Native, and 471 (10%) were other, and around 2140 (47%) of the participants were women. Variation in HbA1c reduction with GLP-1 receptor agonists treatment was associated with rs6923761G→A (Gly168Ser) in the GLP1R (0·08% [95% CI 0·04-0·12] or 0·9 mmol/mol lower reduction in HbA1c per serine, p=6·0 × 10-5) and low frequency variants in ARRB1 (optimal sequence kernel association test p=6·7 × 10-8), largely driven by rs140226575G→A (Thr370Met; 0·25% [SE 0·06] or 2·7 mmol/mol  [SE 0·7] greater HbA1c reduction per methionine, p=5·2 × 10-6). A similar effect size for the ARRB1 Thr370Met was seen in Hispanic and American Indian or Alaska Native populations who have a higher frequency of this variant (6-11%) than in White European populations. Combining these two genes identified 4% of the population who had a 30% greater reduction in HbA1c than the 9% of the population with the worse response. INTERPRETATION: This genome-wide pharmacogenomic study of GLP-1 receptor agonists provides novel biological and clinical insights. Clinically, when genotype is routinely available at the point of prescribing, individuals with ARRB1 variants might benefit from earlier initiation of GLP-1 receptor agonists. FUNDING: Innovative Medicines Initiative and the Wellcome Trust.

Initial Insights into the Genetic Variation Associated with Metformin Treatment Failure in Youth with Type 2 Diabetes.

Metformin is the first-line treatment for type 2 diabetes (T2D) in youth but with limited sustained glycemic response. To identify common variants associated with metformin response, we used a genome-wide approach in 506 youth from the Treatment Options for Type 2 Diabetes in Adolescents and Youth (TODAY) study and examined the relationship between T2D partitioned polygenic scores (pPS), glycemic traits, and metformin response in these youth. Several variants met a suggestive threshold (P < 1 × 10-6), though none including published adult variants reached genome-wide significance. We pursued replication of top nine variants in three cohorts, and rs76195229 in ATRNL1 was associated with worse metformin response in the Metformin Genetics Consortium (n = 7,812), though statistically not being significant after Bonferroni correction (P = 0.06). A higher ß-cell pPS was associated with a lower insulinogenic index (P = 0.02) and C-peptide (P = 0.047) at baseline and higher pPS related to two insulin resistance processes were associated with increased C-peptide at baseline (P = 0.04,0.02). Although pPS were not associated with changes in glycemic traits or metformin response, our results indicate a trend in the association of the ß-cell pPS with reduced ß-cell function over time. Our data show initial evidence for genetic variation associated with metformin response in youth with T2D.

Heterogeneity in phenotype, disease progression and drug response in type 2 diabetes.

Type 2 diabetes (T2D) is a complex chronic disease characterized by considerable phenotypic heterogeneity. In this study, we applied a reverse graph embedding method to routinely collected data from 23,137 Scottish patients with newly diagnosed diabetes to visualize this heterogeneity and used partitioned diabetes polygenic risk scores to gain insight into the underlying biological processes. Overlaying risk of progression to outcomes of insulin requirement, chronic kidney disease, referable diabetic retinopathy and major adverse cardiovascular events, we show how these risks differ by patient phenotype. For example, patients at risk of retinopathy are phenotypically different from those at risk of cardiovascular events. We replicated our findings in the UK Biobank and the ADOPT clinical trial, also showing that the pattern of diabetes drug monotherapy response differs for different drugs. Overall, our analysis highlights how, in a European population, underlying phenotypic variation drives T2D onset and affects subsequent diabetes outcomes and drug response, demonstrating the need to incorporate these factors into personalized treatment approaches for the management of T2D.

Young-onset diabetes in Asian Indians is associated with lower measured and genetically determined beta cell function.

AIMS/HYPOTHESIS: South Asians in general, and Asian Indians in particular, have higher risk of type 2 diabetes compared with white Europeans, and a younger age of onset. The reasons for the younger age of onset in relation to obesity, beta cell function and insulin sensitivity are under-explored. METHODS: Two cohorts of Asian Indians, the ICMR-INDIAB cohort (Indian Council of Medical Research-India Diabetes Study) and the DMDSC cohort (Dr Mohan's Diabetes Specialties Centre), and one of white Europeans, the ESDC (East Scotland Diabetes Cohort), were used. Using a cross-sectional design, we examined the comparative prevalence of healthy, overweight and obese participants with young-onset diabetes, classified according to their BMI. We explored the role of clinically measured beta cell function in diabetes onset in Asian Indians. Finally, the comparative distribution of a partitioned polygenic score (pPS) for risk of diabetes due to poor beta cell function was examined. Replication of the genetic findings was sought using data from the UK Biobank. RESULTS: The prevalence of young-onset diabetes with normal BMI was 9.3% amongst white Europeans and 24-39% amongst Asian Indians. In Asian Indians with young-onset diabetes, after adjustment for family history of type 2 diabetes, sex, insulin sensitivity and HDL-cholesterol, stimulated C-peptide was 492 pmol/ml (IQR 353-616, p<0.0001) lower in lean compared with obese individuals. Asian Indians in our study, and South Asians from the UK Biobank, had a higher number of risk alleles than white Europeans. After weighting the pPS for beta cell function, Asian Indians have lower genetically determined beta cell function than white Europeans (p<0.0001). The pPS was associated with age of diagnosis in Asian Indians but not in white Europeans. The pPS explained 2% of the variation in clinically measured beta cell function, and 1.2%, 0.97%, and 0.36% of variance in age of diabetes amongst Asian Indians with normal BMI, or classified as overweight and obese BMI, respectively. CONCLUSIONS/INTERPRETATION: The prevalence of lean BMI in young-onset diabetes is over two times higher in Asian Indians compared with white Europeans. This phenotype of lean, young-onset diabetes appears driven in part by lower beta cell function. We demonstrate that Asian Indians with diabetes also have lower genetically determined beta cell function.

Diabetes status modifies the long-term effect of lipoprotein-associated phospholipase A2 on major coronary events.

AIMS/HYPOTHESIS: Lipoprotein-associated phospholipase A2 (Lp-PLA2) activity has an independent prognostic association with major coronary events (MCE). However, no study has investigated whether type 2 diabetes status modifies the effect of Lp-PLA2 activity or inhibition on the risk of MCE. We investigate the interaction between diabetes status and Lp-PLA2 activity with risk of MCE. Subsequently, we test the resulting hypothesis that diabetes status will play a role in modifying the efficacy of an Lp-PLA2 inhibitor. METHODS: A retrospective cohort study design was utilised in two study populations. Discovery analyses were performed in the Genetics of Diabetes Audit and Research in Tayside Scotland (GoDARTS) cohort based in Scotland, UK. Participants were categorised by type 2 diabetes control status: poorly controlled (HbA1c ≥ 48 mmol/mol or ≥6.5%) and well-controlled (HbA1c < 48 mmol/mol or <6.5%) diabetes (n = 7420). In a secondary analysis of the Stabilization of Atherosclerotic Plaque by Initiation of Darapladib Therapy (STABILITY) trial of Lp-PLA2 inhibitor (darapladib) efficacy, 15,828 participants were stratified post hoc by type 2 diabetes diagnosis status (diabetes or no diabetes) at time of recruitment. Lp-PLA2 activity was then divided into population-specific quartiles. MCE were determined from linked medical records in GoDARTS and trial records in STABILITY. First, the interaction between diabetes control status and Lp-PLA2 activity on the outcome of MCE was explored in GoDARTS. The effect was replicated in the placebo arm of STABILITY. The effect of Lp-PLA2 on MCE was then examined in models stratified by diabetes status. This helped determine participants at higher risk. Finally, the effect of Lp-PLA2 inhibition was assessed in STABILITY in the higher risk group. Cox proportional hazards models adjusted for confounders were used to assess associations. RESULTS: In GoDARTS, a significant interaction between increased Lp-PLA2 activity (continuous and quartile divided) and diabetes control status was observed in the prediction of MCE (p < 0.0001). These effects were replicated in the placebo arm of STABILITY (p < 0.0001). In GoDARTS, stratified analyses showed that, among individuals with poorly controlled diabetes, the hazards of MCE for those with high (Q4) Lp-PLA2 activity was 1.19 compared with individuals with lower (Q1-3) Lp-PLA2 activity (95% CI 1.11, 1.38; p < 0.0001) and 1.35 (95% CI 1.16, 1.57; p < 0.0001) when compared with those with the lowest activity (Q1). Those in the higher risk group were identified as individuals with the highest Lp-PLA2 activity (Q4) and poorly controlled diabetes or diabetes. Based on these observations in untreated populations, we hypothesised that the Lp-PLA2 inhibitor would have more benefit in this higher risk group. In this risk group, Lp-PLA2 inhibitor use was associated with a 33% reduction in MCE compared with placebo (HR 0.67 [95% CI 0.50, 0.90]; p = 0.008). In contrast, Lp-PLA2 inhibitor showed no efficacy in individuals with low activity, regardless of diabetes status, or among those with no baseline diabetes and high Lp-PLA2 activity. CONCLUSIONS/INTERPRETATION: These results support the hypothesis that diabetes status modifies the association between Lp-PLA2 activity and MCE. These results suggest that cardiovascular morbidity and mortality associated with Lp-PLA2 activity is especially important in patients with type 2 diabetes, particularly those with worse glycaemic control. Further investigation of the effects of Lp-PLA2 inhibition in diabetes appears warranted. DATA AVAILABILITY: STABILITY trial data are available from clinicaltrials.gov repository through the GlaxoSmithKline clinical study register https://clinicaltrials.gov/ct2/show/NCT00799903 . GoDARTS datasets generated during and/or analysed during the current study are available following request to the GoDARTS Access Managements Group https://godarts.org/scientific-community/ .

Genomic editing of metformin efficacy-associated genetic variants in SLC47A1 does not alter SLC47A1 expression.

Several pharmacogenetics studies have identified an association between a greater metformin-dependent reduction in HbA1c levels and the minor A allele at rs2289669 in intron 10 of SLC47A1, encoding multidrug and toxin extrusion 1 (MATE1), a presumed metformin transporter. It is currently unknown if the rs2289669 locus is a cis-eQTL, which would validate its role as predictor of metformin efficacy. We looked at association between common genetic variants in the SLC47A1 gene region and HbA1c reduction after metformin treatment using locus-wise meta-analysis from the MetGen consortium. CRISPR-Cas9 was applied to perform allele editing of, or genomic deletion around, rs2289669 and of the closely linked rs8065082 in HepG2 cells. The genome-edited cells were evaluated for SLC47A1 expression and splicing. None of the common variants including rs2289669 showed significant association with metformin response. Genomic editing of either rs2289669 or rs8065082 did not alter SLC47A1 expression or splicing. Experimental and in silico analyses show that the rs2289669-containing haploblock does not appear to carry genetic variants that could explain its previously reported association with metformin efficacy.

Common Statin Intolerance Variants in ABCB1 and LILRB5 Show Synergistic Effects on Statin Response: An Observational Study Using Electronic Health Records.

Background: Statin intolerance impacts approximately 10% of statin users, with side effects ranging from mild myalgia to extreme intolerance resulting in myopathy and rhabdomyolysis. Statin intolerance results in poor adherence to therapy and can impact statin efficacy. Many genetic variants are associated with statin intolerance. The effect of these variants on statin efficacy has not been systematically explored. Methods: Using longitudinal electronic health records and genetic biobank data from Tayside, Scotland, we examined the effect of seven genetic variants with previously reported associations with simvastatin or atorvastatin intolerance on the outcome of statin response. Statin response was measured by the reduction achieved when comparing pre- and post-statin non-high-density lipoprotein-cholesterol (non-HDL-C). Post-treatment statin response was limited to non-HDL-C measured within 6months of therapy initiation. Univariate and multivariable linear regression models were used to assess the main and adjusted effect of the variants on statin efficacy. Results: Around 9,401 statin users met study inclusion criteria, of whom 8,843 were first prescribed simvastatin or atorvastatin. The average difference in post-treatment compared to pre-treatment non-HDL-cholesterol was 1.45 (±1.04) mmol/L. In adjusted analyses, only two variants, one in the gene ATP-binding cassette transporter B1 (ABCB1; rs1045642), and one in leukocyte immunoglobulin like receptor B5 (LILRB5; rs12975366), were associated with statin efficacy. In ABCB1, homozygous carriers of the C allele at rs1045642 had 0.06mmol/L better absolute reduction in non-HDL-cholesterol than carriers of the T allele (95% CI: 0.01, 0.1). In LILRB5 (rs12975366), carriers of the C allele had 0.04mmol/L better absolute reduction compared to those homozygous for the T allele (95% CI: 0.004, 0.08). When combined into a two-variant risk score, individuals with both the rs1045642-CC genotype and the rs12975366-TC or CC genotype had a 0.11mmol/L greater absolute reduction in non-HDL-cholesterol compared to those with rs1045642-TC or TT genotype and the rs12975366-TT genotype (95% CI: 0.05, 0.16; p<0.001). Conclusion: We report two genetic variants for statin adverse drug reactions (ADRs) that are associated with statin efficacy. While the ABCB1 variant has been shown to have an association with statin pharmacokinetics, no similar evidence for LILRB5 has been reported. These findings highlight the value of genetic testing to deliver precision therapeutics to statin users.

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.

Polymorphism in INSR Locus Modifies Risk of Atrial Fibrillation in Patients on Thyroid Hormone Replacement Therapy.

AIMS: Atrial fibrillation (AF) is a risk for patients receiving thyroid hormone replacement therapy. No published work has focused on pharmacogenetics relevant to thyroid dysfunction and AF risk. We aimed to assess the effect of L-thyroxine on AF risk stratified by a variation in a candidate gene. METHODS AND RESULTS: A retrospective follow-up study was done among European Caucasian patients from the Genetics of Diabetes Audit and Research in Tayside Scotland cohort (Scotland, United Kingdom). Linked data on biochemistry, prescribing, hospital admissions, demographics, and genetic biobank were used to ascertain patients on L-thyroxine and diagnosis of AF. A GWAS-identified insulin receptor-INSR locus (rs4804416) was the candidate gene. Cox survival models and sensitivity analyses by taking competing risk of death into account were used. Replication was performed in additional sample (The Genetics of Scottish Health Research register, GoSHARE), and meta-analyses across the results of the study and replication cohorts were done. We analyzed 962 exposed to L-thyroxine and 5,840 unexposed patients who were rs4804416 genotyped. The rarer G/G genotype was present in 18% of the study population. The total follow-up was up to 20 years, and there was a significant increased AF risk for patients homozygous carriers of the G allele exposed to L-thyroxine (RHR = 2.35, P = 1.6e-02). The adjusted increased risk was highest within the first 3 years of exposure (RHR = 9.10, P = 8.5e-04). Sensitivity analysis yielded similar results. Effects were replicated in GoSHARE (n = 3,190). CONCLUSION: Homozygous G/G genotype at the INSR locus (rs4804416) is associated with an increased risk of AF in patients on L-thyroxine, independent of serum of free thyroxine and thyroid-stimulating hormone serum concentrations.

Utilizing Large Electronic Medical Record Data Sets to Identify Novel Drug-Gene Interactions for Commonly Used Drugs.

Real-world prescribing of drugs differs from the experimental systems, physiological-pharmacokinetic models, and clinical trials used in drug development and licensing, with drugs often used in patients with multiple comorbidities with resultant polypharmacy. The increasing availability of large biobanks linked to electronic healthcare records enables the potential to identify novel drug-gene interactions in large populations of patients. In this study we used three Scottish cohorts and UK Biobank to identify drug-gene interactions for the 50 most commonly used drugs and 162 variants in genes involved in drug pharmacokinetics. We defined two phenotypes based upon prescribing behavior-drug-stop or dose-decrease. Using this approach, we replicate 11 known drug-gene interactions including, for example, CYP2C9/CYP2C8 variants and sulfonylurea/thiazolidinedione prescribing and ABCB1/ABCG2 variants and statin prescribing. We identify eight novel associations after Bonferroni correction, three of which are replicated or validated in the UK Biobank or have other supporting results: The C-allele at rs4918758 in CYP2C9 was associated with a 25% (15-44%) lower odds of dose reduction of quinine, P = 1.6 × 10-5 ; the A-allele at rs9895420 in ABCC3 was associated with a 46% (24-62%) reduction in odds of dose reduction with doxazosin, P = 1.2 × 10-4 , and altered blood pressure response in the UK Biobank; the CYP2D6*2 variant was associated with a 30% (18-40%) reduction in odds of stopping ramipril treatment, P = 1.01 × 10-5 , with similar results seen for enalapril and lisinopril and with other CYP2D6 variants. This study highlights the scope of using large population bioresources linked to medical record data to explore drug-gene interactions at scale.