An omics-based machine learning approach to predict diabetes progression: a RHAPSODY study

AIMS/HYPOTHESIS

People with type 2 diabetes are heterogeneous in their disease trajectory, with some progressing more quickly to insulin initiation than others. Although classical biomarkers such as age, HbA1c and diabetes duration are associated with glycaemic progression, it is unclear how well such variables predict insulin initiation or requirement and whether newly identified markers have added predictive value.

METHODS

In two prospective cohort studies as part of IMI-RHAPSODY, we investigated whether clinical variables and three types of molecular markers (metabolites, lipids, proteins) can predict time to insulin requirement using different machine learning approaches (lasso, ridge, GRridge, random forest). Clinical variables included age, sex, HbA1c, HDL-cholesterol and C-peptide. Models were run with unpenalised clinical variables (i.e. always included in the model without weights) or penalised clinical variables, or without clinical variables. Model development was performed in one cohort and the model was applied in a second cohort. Model performance was evaluated using Harrel's C statistic.

RESULTS

Of the 585 individuals from the Hoorn Diabetes Care System (DCS) cohort, 69 required insulin during follow-up (1.0-11.4 years); of the 571 individuals in the Genetics of Diabetes Audit and Research in Tayside Scotland (GoDARTS) cohort, 175 required insulin during follow-up (0.3-11.8 years). Overall, the clinical variables and proteins were selected in the different models most often, followed by the metabolites. The most frequently selected clinical variables were HbA1c (18 of the 36 models, 50%), age (15 models, 41.2%) and C-peptide (15 models, 41.2%). Base models (age, sex, BMI, HbA1c) including only clinical variables performed moderately in both the DCS discovery cohort (C statistic 0.71 [95% CI 0.64, 0.79]) and the GoDARTS replication cohort (C 0.71 [95% CI 0.69, 0.75]). A more extensive model including HDL-cholesterol and C-peptide performed better in both cohorts (DCS, C 0.74 [95% CI 0.67, 0.81]; GoDARTS, C 0.73 [95% CI 0.69, 0.77]). Two proteins, lactadherin and proto-oncogene tyrosine-protein kinase receptor, were most consistently selected and slightly improved model performance.

CONCLUSIONS/INTERPRETATION

Using machine learning approaches, we show that insulin requirement risk can be modestly well predicted by predominantly clinical variables. Inclusion of molecular markers improves the prognostic performance beyond that of clinical variables by up to 5%. Such prognostic models could be useful for identifying people with diabetes at high risk of progressing quickly to treatment intensification.

DATA AVAILABILITY

Summary statistics of lipidomic, proteomic and metabolomic data are available from a Shiny dashboard at https://rhapdata-app.vital-it.ch .

Phenotype-based targeted treatment of SGLT2 inhibitors and GLP-1 receptor agonists in type 2 diabetes

AIMS/HYPOTHESIS

A precision medicine approach in type 2 diabetes could enhance targeting specific glucose-lowering therapies to individual patients most likely to benefit. We aimed to use the recently developed Bayesian causal forest (BCF) method to develop and validate an individualised treatment selection algorithm for two major type 2 diabetes drug classes, sodium-glucose cotransporter 2 inhibitors (SGLT2i) and glucagon-like peptide-1 receptor agonists (GLP1-RA).

METHODS

We designed a predictive algorithm using BCF to estimate individual-level conditional average treatment effects for 12-month glycaemic outcome (HbA1c) between SGLT2i and GLP1-RA, based on routine clinical features of 46,394 people with type 2 diabetes in primary care in England (Clinical Practice Research Datalink; 27,319 for model development, 19,075 for hold-out validation), with additional external validation in 2252 people with type 2 diabetes from Scotland (SCI-Diabetes [Tayside & Fife]). Differences in glycaemic outcome with GLP1-RA by sex seen in clinical data were replicated in clinical trial data (HARMONY programme: liraglutide [n=389] and albiglutide [n=1682]). As secondary outcomes, we evaluated the impacts of targeting therapy based on glycaemic response on weight change, tolerability and longer-term risk of new-onset microvascular complications, macrovascular complications and adverse kidney events.

RESULTS

Model development identified marked heterogeneity in glycaemic response, with 4787 (17.5%) of the development cohort having a predicted HbA1c benefit >3 mmol/mol (>0.3%) with SGLT2i over GLP1-RA and 5551 (20.3%) having a predicted HbA1c benefit >3 mmol/mol with GLP1-RA over SGLT2i. Calibration was good in hold-back validation, and external validation in an independent Scottish dataset identified clear differences in glycaemic outcomes between those predicted to benefit from each therapy. Sex, with women markedly more responsive to GLP1-RA, was identified as a major treatment effect modifier in both the UK observational datasets and in clinical trial data: HARMONY-7 liraglutide (GLP1-RA): 4.4 mmol/mol (95% credible interval [95% CrI] 2.2, 6.3) (0.4% [95% CrI 0.2, 0.6]) greater response in women than men. Targeting the two therapies based on predicted glycaemic response was also associated with improvements in short-term tolerability and long-term risk of new-onset microvascular complications.

CONCLUSIONS/INTERPRETATION

Precision medicine approaches can facilitate effective individualised treatment choice between SGLT2i and GLP1-RA therapies, and the use of routinely collected clinical features for treatment selection could support low-cost deployment in many countries.

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.

Body Mass Index's influence on arterial hypertension in Type 1 diabetes - A brief report from IMI-SOPHIA study

Information on BMI and risk of developing hypertension in type 1 diabetes (T1D) is scarce, and it comes mostly from cross-sectional analyses. This study underscores a risk of developing hypertension in T1D individuals with high BMI, and this risk appears to be higher than in those with type 2 diabetes.

Trajectories of clinical characteristics, complications and treatment choices in data-driven subgroups of type 2 diabetes

AIMS/HYPOTHESIS

This study aimed to explore the added value of subgroups that categorise individuals with type 2 diabetes by k-means clustering for two primary care registries (the Netherlands and Scotland), inspired by Ahlqvist's novel diabetes subgroups and previously analysed by Slieker et al. METHODS: We used two Dutch and Scottish diabetes cohorts (N=3054 and 6145; median follow-up=11.2 and 12.3 years, respectively) and defined five subgroups by k-means clustering with age at baseline, BMI, HbA1c, HDL-cholesterol and C-peptide. We investigated differences between subgroups by trajectories of risk factor values (random intercept models), time to diabetes-related complications (logrank tests and Cox models) and medication patterns (multinomial logistic models). We also compared directly using the clustering indicators as predictors of progression vs the k-means discrete subgroups. Cluster consistency over follow-up was assessed.

RESULTS

Subgroups' risk factors were significantly different, and these differences remained generally consistent over follow-up. Among all subgroups, individuals with severe insulin resistance faced a significantly higher risk of myocardial infarction both before (HR 1.65; 95% CI 1.40, 1.94) and after adjusting for age effect (HR 1.72; 95% CI 1.46, 2.02) compared with mild diabetes with high HDL-cholesterol. Individuals with severe insulin-deficient diabetes were most intensively treated, with more than 25% prescribed insulin at 10 years of diagnosis. For severe insulin-deficient diabetes relative to mild diabetes, the relative risks for using insulin relative to no common treatment would be expected to increase by a factor of 3.07 (95% CI 2.73, 3.44), holding other factors constant. Clustering indicators were better predictors of progression variation relative to subgroups, but prediction accuracy may improve after combining both. Clusters were consistent over 8 years with an accuracy ranging from 59% to 72%.

CONCLUSIONS/INTERPRETATION

Data-driven subgroup allocations were generally consistent over follow-up and captured significant differences in risk factor trajectories, medication patterns and complication risks. Subgroups serve better as a complement rather than as a basis for compressing clustering indicators.

Pharmacogenetics at scale in real-world bioresources: CYP2C19 and clopidogrel outcomes in UK Biobank

OBJECTIVE

The impact of CYP2C19 genotype on clopidogrel outcomes is one of the most well established pharmacogenetic interactions, supported by robust evidence and recommended by the Food and Drug Administration and clinical pharmacogenetics implementation consortium. However, there is a scarcity of large-scale real-world data on the extent of this pharmacogenetic effect, and clinical testing for the CYP2C19 genotype remains infrequent. This study utilizes the UK Biobank dataset, including 10 365 patients treated with clopidogrel, to offer the largest observational analysis of these pharmacogenetic effects to date.

METHODS

Incorporating time-varying drug exposure and repeated clinical outcome, we adopted semiparametric frailty models to detect and quantify exposure-based effects of CYP2C19 (*2,*17) variants and nongenetic factors on the incidence risks of composite outcomes of death or recurrent hospitalizations due to major adverse cardiovascular events (MACE) or hemorrhage in the entire cohort of clopidogrel-treated patients.

RESULTS

Out of the 10 365 clopidogrel-treated patients, 40% (4115) experienced 10 625 MACE events during an average follow-up of 9.23 years. Individuals who received clopidogrel (coverage >25%) with a CYP2C19*2 loss-of-function allele had a 9.4% higher incidence of MACE [incidence rate ratios (IRR), 1.094; 1.044-1.146], but a 15% lower incidence of hemorrhage (IRR, 0.849; 0.712-0.996). These effects were stronger with high clopidogrel exposure. Conversely, the gain-of-function CYP2C19*17 variant was associated with a 5.3% lower incidence of MACE (IRR, 0.947; 0.903-0.983). Notably, there was no evidence of *2 or *17 effects when clopidogrel exposure was low, confirming the presence of a drug-gene interaction.

CONCLUSION

The impact of CYP2C19 on clinical outcomes in clopidogrel-treated patients is substantial, highlighting the importance of incorporating genotype-based prescribing into clinical practice, regardless of the reason for clopidogrel use or the duration of treatment. Moreover, the methodology introduced in this study can be applied to further real-world investigations of known drug-gene and drug-drug interactions and the discovery of novel interactions.

Multi-omics subgroups associated with glycaemic deterioration in type 2 diabetes: an IMI-RHAPSODY Study

Introduction

Type 2 diabetes (T2D) onset, progression and outcomes differ substantially between individuals. Multi-omics analyses may allow a deeper understanding of these differences and ultimately facilitate personalised treatments. Here, in an unsupervised "bottom-up" approach, we attempt to group T2D patients based solely on -omics data generated from plasma.

Methods

Circulating plasma lipidomic and proteomic data from two independent clinical cohorts, Hoorn Diabetes Care System (DCS) and Genetics of Diabetes Audit and Research in Tayside Scotland (GoDARTS), were analysed using Similarity Network Fusion. The resulting patient network was analysed with Logistic and Cox regression modelling to explore relationships between plasma -omic profiles and clinical characteristics.

Results

From a total of 1,134 subjects in the two cohorts, levels of 180 circulating plasma lipids and 1195 proteins were used to separate patients into two subgroups. These differed in terms of glycaemic deterioration (Hazard Ratio=0.56;0.73), insulin sensitivity and secretion (C-peptide, p=3.7e-11;2.5e-06, DCS and GoDARTS, respectively; Homeostatic model assessment 2 (HOMA2)-B; -IR; -S, p=0.0008;4.2e-11;1.1e-09, only in DCS). The main molecular signatures separating the two groups included triacylglycerols, sphingomyelin, testican-1 and interleukin 18 receptor.

Conclusions

Using an unsupervised network-based fusion method on plasma lipidomics and proteomics data from two independent cohorts, we were able to identify two subgroups of T2D patients differing in terms of disease severity. The molecular signatures identified within these subgroups provide insights into disease mechanisms and possibly new prognostic markers for T2D.

Analysis of type 2 diabetes heterogeneity with a tree-like representation: insights from the prospective German Diabetes Study and the LURIC cohort

BACKGROUND

Heterogeneity in type 2 diabetes can be represented by a tree-like graph structure by use of reversed graph-embedded dimensionality reduction. We aimed to examine whether this approach can be used to stratify key pathophysiological components and diabetes-related complications during longitudinal follow-up of individuals with recent-onset type 2 diabetes.

METHODS

For this cohort analysis, 927 participants aged 18-69 years from the German Diabetes Study (GDS) with recent-onset type 2 diabetes were mapped onto a previously developed two-dimensional tree based on nine simple clinical and laboratory variables, residualised for age and sex. Insulin sensitivity was assessed by a hyperinsulinaemic-euglycaemic clamp, insulin secretion was assessed by intravenous glucose tolerance test, hepatic lipid content was assessed by 1 H magnetic resonance spectroscopy, serum interleukin (IL)-6 and IL-18 were assessed by ELISA, and peripheral and autonomic neuropathy were assessed by functional and clinical measures. Participants were followed up for up to 16 years. We also investigated heart failure and all-cause mortality in 794 individuals with type 2 diabetes undergoing invasive coronary diagnostics from the Ludwigshafen Risk and Cardiovascular Health (LURIC) cohort.

FINDINGS

There were gradients of clamp-measured insulin sensitivity (both dimensions: p<0·0001) and insulin secretion (pdim1<0·0001, pdim2=0·00097) across the tree. Individuals in the region with the lowest insulin sensitivity had the highest hepatic lipid content (n=205, pdim1<0·0001, pdim2=0·037), pro-inflammatory biomarkers (IL-6: n=348, pdim1<0·0001, pdim2=0·013; IL-18: n=350, pdim1<0·0001, pdim2=0·38), and elevated cardiovascular risk (nevents=143, pdim1=0·14, pdim2<0·00081), whereas individuals positioned in the branch with the lowest insulin secretion were more prone to require insulin therapy (nevents=85, pdim1=0·032, pdim2=0·12) and had the highest risk of diabetic sensorimotor polyneuropathy (nevents=184, pdim1=0·012, pdim2=0·044) and cardiac autonomic neuropathy (nevents=118, pdim1=0·0094, pdim2=0·06). In the LURIC cohort, all-cause mortality was highest in the tree branch showing insulin resistance (nevents=488, pdim1=0·12, pdim2=0·0032). Significant gradients differentiated individuals having heart failure with preserved ejection fraction from those who had heart failure with reduced ejection fraction.

INTERPRETATION

These data define the pathophysiological underpinnings of the tree structure, which has the potential to stratify diabetes-related complications on the basis of routinely available variables and thereby expand the toolbox of precision diabetes diagnosis.

FUNDING

German Diabetes Center, German Federal Ministry of Health, Ministry of Culture and Science of the state of North Rhine-Westphalia, German Federal Ministry of Education and Research, German Diabetes Association, German Center for Diabetes Research, European Community, German Research Foundation, and Schmutzler Stiftung.

Associations of offspring birthweight and placental weight with subsequent parental coronary heart disease: survival regression using the walker cohort

Low birth weight (BW) is consistently correlated with increased parental risk of subsequent cardiovascular disease, but the links with offspring placental weight (PW) are mostly unexplored. We have investigated the associations between parental coronary heart disease (CHD) and offspring BW and PW using the Walker cohort, a collection of 48,000 birth records from Dundee, Scotland, from the 1950s and 1960s. We linked the medical history of 13,866 mothers and 8,092 fathers to their offspring's records and performed Cox survival analyses modelling maternal and paternal CHD risk by their offspring's BW, PW, and the ratio between both measurements. We identified negative associations between offspring BW and both maternal (hazard ratio [HR]: 0.91, 95% confidence interval [CI]: 0.88-0.95) and paternal (HR: 0.96, 95% CI: 0.93-1.00) CHD risk, the stronger maternal correlation being consistent with previous reports. Offspring PW to BW ratio was positively associated with maternal CHD risk (HR: 1.14, 95% CI: 1.08-1.21), but the associations with paternal CHD were not significant. These analyses provide additional evidence for intergenerational associations between early growth and parental disease, identifying directionally opposed correlations of maternal CHD with offspring BW and PW, and highlight the importance of the placenta as a determinant of early development and adult disease.

Genome-wide association study of placental weight identifies distinct and shared genetic influences between placental and fetal growth

A well-functioning placenta is essential for fetal and maternal health throughout pregnancy. Using placental weight as a proxy for placental growth, we report genome-wide association analyses in the fetal (n = 65,405), maternal (n = 61,228) and paternal (n = 52,392) genomes, yielding 40 independent association signals. Twenty-six signals are classified as fetal, four maternal and three fetal and maternal. A maternal parent-of-origin effect is seen near KCNQ1. Genetic correlation and colocalization analyses reveal overlap with birth weight genetics, but 12 loci are classified as predominantly or only affecting placental weight, with connections to placental development and morphology, and transport of antibodies and amino acids. Mendelian randomization analyses indicate that fetal genetically mediated higher placental weight is causally associated with preeclampsia risk and shorter gestational duration. Moreover, these analyses support the role of fetal insulin in regulating placental weight, providing a key link between fetal and placental growth.

Metformin increases the uptake of glucose into the gut from the circulation in high-fat diet-fed male mice, which is enhanced by a reduction in whole-body Slc2a2 expression

OBJECTIVES

Metformin is the first line therapy recommended for type 2 diabetes. However, the precise mechanism of action remains unclear and up to a quarter of patients show some degree of intolerance to the drug, with a similar number showing poor response to treatment, limiting its effectiveness. A better understanding of the mechanism of action of metformin may improve its clinical use. SLC2A2 (GLUT2) is a transmembrane facilitated glucose transporter, with important roles in the liver, gut and pancreas. Our group previously identified single nucleotide polymorphisms in the human SLC2A2 gene, which were associated with reduced transporter expression and an improved response to metformin treatment. The aims of this study were to model Slc2a2 deficiency and measure the impact on glucose homoeostasis and metformin response in mice.

METHODS

We performed extensive metabolic phenotyping and 2-deoxy-2-[18F]fluoro-d-glucose ([18F]FDG)-positron emission tomography (PET) analysis of gut glucose uptake in high-fat diet-fed (HFD) mice with whole-body reduced Slc2a2 (Slc2a2+/-) and intestinal Slc2a2 KO, to assess the impact of metformin treatment.

RESULTS

Slc2a2 partial deficiency had no major impact on body weight and insulin sensitivity, however mice with whole-body reduced Slc2a2 expression (Slc2a2+/-) developed an age-related decline in glucose homoeostasis (as measured by glucose tolerance test) compared to wild-type (Slc2a2+/+) littermates. Glucose uptake into the gut from the circulation was enhanced by metformin exposure in Slc2a2+/+ animals fed HFD and this action of the drug was significantly higher in Slc2a2+/- animals. However, there was no effect of specifically knocking-out Slc2a2 in the mouse intestinal epithelial cells.

CONCLUSIONS

Overall, this work identifies a differential metformin response, dependent on expression of the SLC2A2 glucose transporter, and also adds to the growing evidence that metformin efficacy includes modifying glucose transport in the gut. We also describe a novel and important role for this transporter in maintaining efficient glucose homoeostasis during ageing.

Understanding diabetes heterogeneity: key steps towards precision medicine in diabetes

Diabetes is a highly heterogeneous condition; yet, it is diagnosed by measuring a single blood-borne metabolite, glucose, irrespective of aetiology. Although pragmatically helpful, disease classification can become complex and limit advances in research and medical care. Here, we describe diabetes heterogeneity, highlighting recent approaches that could facilitate management by integrating three disease models across all forms of diabetes, namely, the palette model, the threshold model and the gradient model. Once diabetes has developed, further worsening of established diabetes and the subsequent emergence of diabetes complications are kept in check by multiple processes designed to prevent or circumvent metabolic dysfunction. The impact of any given disease risk factor will vary from person-to-person depending on their background, diabetes-related propensity, and environmental exposures. Defining the consequent heterogeneity within diabetes through precision medicine, both in terms of diabetes risk and risk of complications, could improve health outcomes today and shine a light on avenues for novel therapy in the future.

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.

Second international consensus report on gaps and opportunities for the clinical translation of precision diabetes medicine

Precision medicine is part of the logical evolution of contemporary evidence-based medicine that seeks to reduce errors and optimize outcomes when making medical decisions and health recommendations. Diabetes affects hundreds of millions of people worldwide, many of whom will develop life-threatening complications and die prematurely. Precision medicine can potentially address this enormous problem by accounting for heterogeneity in the etiology, clinical presentation and pathogenesis of common forms of diabetes and risks of complications. This second international consensus report on precision diabetes medicine summarizes the findings from a systematic evidence review across the key pillars of precision medicine (prevention, diagnosis, treatment, prognosis) in four recognized forms of diabetes (monogenic, gestational, type 1, type 2). These reviews address key questions about the translation of precision medicine research into practice. Although not complete, owing to the vast literature on this topic, they revealed opportunities for the immediate or near-term clinical implementation of precision diabetes medicine; furthermore, we expose important gaps in knowledge, focusing on the need to obtain new clinically relevant evidence. Gaps include the need for common standards for clinical readiness, including consideration of cost-effectiveness, health equity, predictive accuracy, liability and accessibility. Key milestones are outlined for the broad clinical implementation of precision diabetes medicine.

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.

Genetic analysis of blood molecular phenotypes reveals common properties in the regulatory networks affecting complex traits

We evaluate the shared genetic regulation of mRNA molecules, proteins and metabolites derived from whole blood from 3029 human donors. We find abundant allelic heterogeneity, where multiple variants regulate a particular molecular phenotype, and pleiotropy, where a single variant associates with multiple molecular phenotypes over multiple genomic regions. The highest proportion of share genetic regulation is detected between gene expression and proteins (66.6%), with a further median shared genetic associations across 49 different tissues of 78.3% and 62.4% between plasma proteins and gene expression. We represent the genetic and molecular associations in networks including 2828 known GWAS variants, showing that GWAS variants are more often connected to gene expression in trans than other molecular phenotypes in the network. Our work provides a roadmap to understanding molecular networks and deriving the underlying mechanism of action of GWAS variants using different molecular phenotypes in an accessible tissue.

Competing risks analysis for neutrophil to lymphocyte ratio as a predictor of diabetic retinopathy incidence in the Scottish population

BACKGROUND

Diabetic retinopathy (DR) is a major sight-threatening microvascular complication in individuals with diabetes. Systemic inflammation combined with oxidative stress is thought to capture most of the complexities involved in the pathology of diabetic retinopathy. A high level of neutrophil-lymphocyte ratio (NLR) is an indicator of abnormal immune system activity. Current estimates of the association of NLR with diabetes and its complications are almost entirely derived from cross-sectional studies, suggesting that the nature of the reported association may be more diagnostic than prognostic. Therefore, in the present study, we examined the utility of NLR as a biomarker to predict the incidence of DR in the Scottish population.

METHODS

The incidence of DR was defined as the time to the first diagnosis of R1 or above grade in the Scottish retinopathy grading scheme from type 2 diabetes diagnosis. The effect of NLR and its interactions were explored using a competing risks survival model adjusting for other risk factors and accounting for deaths. The Fine and Gray subdistribution hazard model (FGR) was used to predict the effect of NLR on the incidence of DR.

RESULTS

We analysed data from 23,531 individuals with complete covariate information. At 10 years, 8416 (35.8%) had developed DR and 2989 (12.7%) were lost to competing events (death) without developing DR and 12,126 individuals did not have DR. The median (interquartile range) level of NLR was 2.04 (1.5 to 2.7). The optimal NLR cut-off value to predict retinopathy incidence was 3.04. After accounting for competing risks at 10 years, the cumulative incidence of DR and deaths without DR were 50.7% and 21.9%, respectively. NLR was associated with incident DR in both Cause-specific hazard (CSH = 1.63; 95% CI: 1.28-2.07) and FGR models the subdistribution hazard (sHR = 2.24; 95% CI: 1.70-2.94). Both age and HbA1c were found to modulate the association between NLR and the risk of DR.

CONCLUSIONS

The current study suggests that NLR has a promising potential to predict DR incidence in the Scottish population, especially in individuals less than 65 years and in those with well-controlled glycaemic status.

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.

Common and Distinct Genetic Architecture of Age at Diagnosis of Diabetes in South Indian and European Populations

OBJECTIVE

South Asians are diagnosed with type 2 diabetes (T2D) more than a decade earlier in life than seen in European populations. We hypothesized that studying the genomics of age of diagnosis in these populations may give insight into the earlier age diagnosis of T2D among individuals of South Asian descent.

RESEARCH DESIGN AND METHODS

We conducted a meta-analysis of genome-wide association studies (GWAS) of age at diagnosis of T2D in 34,001 individuals from four independent cohorts of European and South Asian Indians.

RESULTS

We identified two signals near the TCF7L2 and CDKAL1 genes associated with age at the onset of T2D. The strongest genome-wide significant variants at chromosome 10q25.3 in TCF7L2 (rs7903146; P = 2.4 × 10-12, β = -0.436; SE 0.02) and chromosome 6p22.3 in CDKAL1 (rs9368219; P = 2.29 × 10-8; β = -0.053; SE 0.01) were directionally consistent across ethnic groups and present at similar frequencies; however, both loci harbored additional independent signals that were only present in the South Indian cohorts. A genome-wide signal was also obtained at chromosome 10q26.12 in WDR11 (rs3011366; P = 3.255 × 10-8; β = 1.44; SE 0.25), specifically in the South Indian cohorts. Heritability estimates for the age at diagnosis were much stronger in South Indians than Europeans, and a polygenic risk score constructed based on South Indian GWAS explained ∼2% trait variance.

CONCLUSIONS

Our findings provide a better understanding of ethnic differences in the age at diagnosis and indicate the potential importance of ethnic differences in the genetic architecture underpinning T2D.

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.

A publication-wide association study (PWAS), historical language models to prioritise novel therapeutic drug targets

Most biomedical knowledge is published as text, making it challenging to analyse using traditional statistical methods. In contrast, machine-interpretable data primarily comes from structured property databases, which represent only a fraction of the knowledge present in the biomedical literature. Crucial insights and inferences can be drawn from these publications by the scientific community. We trained language models on literature from different time periods to evaluate their ranking of prospective gene-disease associations and protein-protein interactions. Using 28 distinct historical text corpora of abstracts published between 1995 and 2022, we trained independent Word2Vec models to prioritise associations that were likely to be reported in future years. This study demonstrates that biomedical knowledge can be encoded as word embeddings without the need for human labelling or supervision. Language models effectively capture drug discovery concepts such as clinical tractability, disease associations, and biochemical pathways. Additionally, these models can prioritise hypotheses years before their initial reporting. Our findings underscore the potential for extracting yet-to-be-discovered relationships through data-driven approaches, leading to generalised biomedical literature mining for potential therapeutic drug targets. The Publication-Wide Association Study (PWAS) enables the prioritisation of under-explored targets and provides a scalable system for accelerating early-stage target ranking, irrespective of the specific disease of interest.

Proteomic Pathways across Ejection Fraction Spectrum in Heart Failure: an EXSCEL Substudy

Background

Ejection fraction (EF) is a key component of heart failure (HF) classification, including the increasingly codified HF with mildly reduced EF (HFmrEF) category. However, the biologic basis of HFmrEF as an entity distinct from HF with preserved EF (HFpEF) and reduced EF (HFrEF) has not been well characterized.

Methods

The EXSCEL trial randomized participants with type 2 diabetes (T2DM) to once-weekly exenatide (EQW) vs. placebo. For this study, profiling of ∼5000 proteins using the SomaLogic SomaScan platform was performed in baseline and 12-month serum samples from N=1199 participants with prevalent HF at baseline. Principal component analysis (PCA) and ANOVA (FDR p<0.1) were used to determine differences in proteins between three EF groups, as previously curated in EXSCEL (EF>55% [HFpEF], EF 40-55% [HFmrEF], EF<40% [HFrEF]). Cox proportional hazards was used to assess association between baseline levels of significant proteins, and changes in protein level between baseline and 12-month, with time-to-HF hospitalization. Mixed models were used to assess whether significant proteins changed differentially with exenatide vs. placebo therapy.

Results

Of N=1199 EXSCEL participants with prevalent HF, 284 (24%), 704 (59%) and 211 (18%) had HFpEF, HFmrEF and HFrEF, respectively. Eight PCA protein factors and 221 individual proteins within these factors differed significantly across the three EF groups. Levels of the majority of proteins (83%) demonstrated concordance between HFmrEF and HFpEF, but higher levels in HFrEF, predominated by the domain of extracellular matrix regulation, e.g. COL28A1 and tenascin C [TNC]; p<0.0001. Concordance between HFmrEF and HFrEF was observed in a minority of proteins (1%) including MMP-9 (p<0.0001). Biologic pathways of epithelial mesenchymal transition, ECM receptor interaction, complement and coagulation cascades, and cytokine receptor interaction demonstrated enrichment among proteins with the dominant pattern, i.e. HFmrEF-HFpEF concordance. Baseline levels of 208 (94%) of the 221 proteins were associated with time-to-incident HF hospitalization including domains of extracellular matrix (COL28A1, TNC), angiogenesis (ANG2, VEGFa, VEGFd), myocyte stretch (NT-proBNP), and renal function (cystatin-C). Change in levels of 10 of the 221 proteins from baseline to 12 months (including increase in TNC) predicted incident HF hospitalization (p<0.05). Levels of 30 of the 221 significant proteins (including TNC, NT-proBNP, ANG2) were reduced differentially by EQW compared with placebo (interaction p<0.0001).

Conclusions

In this HF substudy of a large clinical trial of people with T2DM, we found that serum levels of most proteins across multiple biologic domains were similar between HFmrEF and HFpEF. HFmrEF may be more biologically similar to HFpEF than HFrEF, and specific related biomarkers may offer unique data on prognosis and pharmacotherapy modification with variability by EF.

Potential Value of Identifying Type 2 Diabetes Subgroups for Guiding Intensive Treatment: A Comparison of Novel Data-Driven Clustering With Risk-Driven Subgroups

OBJECTIVE

To estimate the impact on lifetime health and economic outcomes of different methods of stratifying individuals with type 2 diabetes, followed by guideline-based treatment intensification targeting BMI and LDL in addition to HbA1c.

RESEARCH DESIGN AND METHODS

We divided 2,935 newly diagnosed individuals from the Hoorn Diabetes Care System (DCS) cohort into five Risk Assessment and Progression of Diabetes (RHAPSODY) data-driven clustering subgroups (based on age, BMI, HbA1c, C-peptide, and HDL) and four risk-driven subgroups by using fixed cutoffs for HbA1c and risk of cardiovascular disease based on guidelines. The UK Prospective Diabetes Study Outcomes Model 2 estimated discounted expected lifetime complication costs and quality-adjusted life-years (QALYs) for each subgroup and across all individuals. Gains from treatment intensification were compared with care as usual as observed in DCS. A sensitivity analysis was conducted based on Ahlqvist subgroups.

RESULTS

Under care as usual, prognosis in the RHAPSODY data-driven subgroups ranged from 7.9 to 12.6 QALYs. Prognosis in the risk-driven subgroups ranged from 6.8 to 12.0 QALYs. Compared with homogenous type 2 diabetes, treatment for individuals in the high-risk subgroups could cost 22.0% and 25.3% more and still be cost effective for data-driven and risk-driven subgroups, respectively. Targeting BMI and LDL in addition to HbA1c might deliver up to 10-fold increases in QALYs gained.

CONCLUSIONS

Risk-driven subgroups better discriminated prognosis. Both stratification methods supported stratified treatment intensification, with the risk-driven subgroups being somewhat better in identifying individuals with the most potential to benefit from intensive treatment. Irrespective of stratification approach, better cholesterol and weight control showed substantial potential for health gains.

Identification of biomarkers for glycaemic deterioration in type 2 diabetes

We identify biomarkers for disease progression in three type 2 diabetes cohorts encompassing 2,973 individuals across three molecular classes, metabolites, lipids and proteins. Homocitrulline, isoleucine and 2-aminoadipic acid, eight triacylglycerol species, and lowered sphingomyelin 42:2;2 levels are predictive of faster progression towards insulin requirement. Of ~1,300 proteins examined in two cohorts, levels of GDF15/MIC-1, IL-18Ra, CRELD1, NogoR, FAS, and ENPP7 are associated with faster progression, whilst SMAC/DIABLO, SPOCK1 and HEMK2 predict lower progression rates. In an external replication, proteins and lipids are associated with diabetes incidence and prevalence. NogoR/RTN4R injection improved glucose tolerance in high fat-fed male mice but impaired it in male db/db mice. High NogoR levels led to islet cell apoptosis, and IL-18R antagonised inflammatory IL-18 signalling towards nuclear factor kappa-B in vitro. This comprehensive, multi-disciplinary approach thus identifies biomarkers with potential prognostic utility, provides evidence for possible disease mechanisms, and identifies potential therapeutic avenues to slow diabetes progression.

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.

Metformin: evidence from preclinical and clinical studies for potential novel applications in cardiovascular disease

INTRODUCTION

For a long time, metformin has been the first-line treatment for glycaemic control in type 2 diabetes, however, the results of recent cardiovascular outcome trials of sodium-glucose co-transporter 2 inhibitors and glucagon-like peptide 1 receptor agonists have caused many to question metformin's position in the guidelines. Although there are several plausible mechanisms by which metformin might have beneficial cardiovascular effects, for example its anti-inflammatory effects and metabolic properties, and numerous observational data suggesting improved cardiovascular outcomes with metformin use, the main randomised clinical trial data for metformin was published over 20 years ago. Nevertheless, the overwhelming majority of participants in contemporary type 2 diabetes trials were prescribed metformin.

AREAS COVERED

In this review we will summarise the potential mechanisms of cardiovascular benefit with metformin, before discussing clinical data in individuals with or without diabetes.

EXPERT OPINION

Metformin may have some cardiovascular benefit in patients with and without diabetes, however the majority of clinical trials were small and are before the use SGLT2 inhibitors and GLP1-RAs. Larger contemporary randomised trials with metformin evaluating its cardiovascular benefit are warranted.

Cardiovascular Safety in Type 2 Diabetes With Sulfonylureas as Second-Line Drugs: A Nation-Wide Population-Based Comparative Safety Study

OBJECTIVE

To assess the real-world cardiovascular (CV) safety for sulfonylureas (SU), in comparison with dipeptidylpeptidase-4 inhibitors (DPP4i) and thiazolidinediones (TZD), through development of robust methodology for causal inference in a whole nation study.

RESEARCH DESIGN AND METHODS

A cohort study was performed including people with type 2 diabetes diagnosed in Scotland before 31 December 2017, who failed to reach HbA1c 48 mmol/mol despite metformin monotherapy and initiated second-line pharmacotherapy (SU/DPP4i/TZD) on or after 1 January 2010. The primary outcome was composite major adverse cardiovascular events (MACE), including hospitalization for myocardial infarction, ischemic stroke, heart failure, and CV death. Secondary outcomes were each individual end point and all-cause death. Multivariable Cox proportional hazards regression and an instrumental variable (IV) approach were used to control confounding in a similar way to the randomization process in a randomized control trial.

RESULTS

Comparing SU to non-SU (DPP4i/TZD), the hazard ratio (HR) for MACE was 1.00 (95% CI: 0.91-1.09) from the multivariable Cox regression and 1.02 (0.91-1.13) and 1.03 (0.91-1.16) using two different IVs. For all-cause death, the HR from Cox regression and the two IV analyses was 1.03 (0.94-1.13), 1.04 (0.93-1.17), and 1.03 (0.90-1.17).

CONCLUSIONS

Our findings contribute to the understanding that second-line SU for glucose lowering are unlikely to increase CV risk or all-cause mortality. Given their potent efficacy, microvascular benefits, cost effectiveness, and widespread use, this study supports that SU should remain a part of the global diabetes treatment portfolio.

Higher mitochondrial DNA copy number is associated with metformin-induced weight loss

BACKGROUND

Considerable variability exists in response to metformin with few effective biomarkers to guide the treatment. Here we evaluated whether whole blood derived mitochondrial DNA copy number (mtDNA-CN) is a biomarker of metformin response as measured by glucose reduction or weight loss.

METHODS

Using data from the trial of Metformin (n = 304) and AcaRbose (n = 300) in Chinese as the initial Hypoglycaemic treatment (MARCH), we examined the association between mtDNA-CN and two metformin response outcomes of HbA1c reduction and weight loss. The acarbose arm was used as a comparator group. Whole blood mtDNA-CN was estimated by deep whole genome sequencing with adjustments for confounders. Multiple linear regression and repeated measurement analyses were used to evaluate the association between mtDNA-CN and drug response outcomes.

RESULTS

Here we show that glucose reduction is not significantly associated with mtDNA-CN and in either treatment arm. In the metformin arm, each increase of 1 SD in mtDNA-CN is significantly (P = 0.006) associated with a 0.43 kg more weight loss. Repeated measurement analysis shows that after 16 weeks of metformin monotherapy, patients in the top tertile of mtDNA-CN consistently lost 1.21 kg more weight than those in the bottom tertile (P < 0.001). In comparison, mtDNA-CN is not significantly associated with acarbose-induced weight loss.

CONCLUSIONS

Patients with higher mtDNA-CN are likely to lose more weight upon metformin treatment, suggesting mtDNA-CN as a potential novel biomarker for more effective weight management in type 2 diabetes.

Patient stratification for determining optimal second-line and third-line therapy for type 2 diabetes: the TriMaster study

Precision medicine aims to treat an individual based on their clinical characteristics. A differential drug response, critical to using these features for therapy selection, has never been examined directly in type 2 diabetes. In this study, we tested two hypotheses: (1) individuals with body mass index (BMI) > 30 kg/m², compared to BMI ≤ 30 kg/m², have greater glucose lowering with thiazolidinediones than with DPP4 inhibitors, and (2) individuals with estimated glomerular filtration rate (eGFR) 60-90 ml/min/1.73 m², compared to eGFR >90 ml/min/1.73 m², have greater glucose lowering with DPP4 inhibitors than with SGLT2 inhibitors. The primary endpoint for both hypotheses was the achieved HbA1c difference between strata for the two drugs. In total, 525 people with type 2 diabetes participated in this UK-based randomized, double-blind, three-way crossover trial of 16 weeks of treatment with each of sitagliptin 100 mg once daily, canagliflozin 100 mg once daily and pioglitazone 30 mg once daily added to metformin alone or metformin plus sulfonylurea. Overall, the achieved HbA1c was similar for the three drugs: pioglitazone 59.6 mmol/mol, sitagliptin 60.0 mmol/mol and canagliflozin 60.6 mmol/mol (P = 0.2). Participants with BMI > 30 kg/m², compared to BMI ≤ 30 kg/m², had a 2.88 mmol/mol (95% confidence interval (CI): 0.98, 4.79) lower HbA1c on pioglitazone than on sitagliptin (n = 356, P = 0.003). Participants with eGFR 60-90 ml/min/1.73 m², compared to eGFR >90 ml/min/1.73 m², had a 2.90 mmol/mol (95% CI: 1.19, 4.61) lower HbA1c on sitagliptin than on canagliflozin (n = 342, P = 0.001). There were 2,201 adverse events reported, and 447/525 (85%) randomized participants experienced an adverse event on at least one of the study drugs. In this precision medicine trial in type 2 diabetes, our findings support the use of simple, routinely available clinical measures to identify the drug class most likely to deliver the greatest glycemic reduction for a given patient. (ClinicalTrials.gov registration: NCT02653209 ; ISRCTN registration: 12039221 .).

Development and external validation of multivariable risk models to predict incident and resolved neuropathic pain: a DOLORisk Dundee study

Neuropathic pain is difficult to treat, and an understanding of the risk factors for its onset and resolution is warranted. This study aimed to develop and externally validate two clinical risk models to predict onset and resolution of chronic neuropathic pain. Participants of Generation Scotland: Scottish Family Health Study (GS; general Scottish population; n = 20,221) and Genetic of Diabetes Audit and Research in Tayside Scotland (GoDARTS; n = 5236) were sent a questionnaire on neuropathic pain and followed- -up 18 months later. Chronic neuropathic pain was defined using DN4 scores (≥ 3/7) and pain for 3 months or more. The models were developed in GS using logistic regression with backward elimination based on the Akaike information criterion. External validation was conducted in GoDARTS and assessed model discrimination (ROC and Precision-Recall curves), calibration and clinical utility (decision curve analysis [DCA]). Analysis revealed incidences of neuropathic pain onset (6.0% in GS [236/3903] and 10.7% in GoDARTS [61/571]) and resolution (42.6% in GS [230/540] and 23.7% in GoDARTS [56/236]). Psychosocial and lifestyle factors were included in both onset and resolved prediction models. In GoDARTS, these models showed adequate discrimination (ROC = 0.636 and 0.699), but there was evidence of miscalibration (Intercept = - 0.511 and - 0.424; slope = 0.623 and 0.999). The DCA indicated that the models would provide clinical benefit over a range of possible risk thresholds. To our knowledge, these are the first externally validated risk models for neuropathic pain. The findings are of interest to patients and clinicians in the community, who may take preventative or remedial measures.

Patient preference for second- and third-line therapies in type 2 diabetes: a prespecified secondary endpoint of the TriMaster study

Patient preference is very important for medication selection in chronic medical conditions, like type 2 diabetes, where there are many different drugs available. Patient preference balances potential efficacy with potential side effects. As both aspects of drug response can vary markedly between individuals, this decision could be informed by the patient personally experiencing the alternative medications, as occurs in a crossover trial. In the TriMaster (NCT02653209, ISRCTN12039221), randomized double-blind, three-way crossover trial patients received three different second- or third-line once-daily type 2 diabetes glucose-lowering drugs (pioglitazone 30 mg, sitagliptin 100 mg and canagliflozin 100 mg). As part of a prespecified secondary endpoint, we examined patients' drug preference after they had tried all three drugs. In total, 448 participants were treated with all three drugs which overall showed similar glycemic control (HbA1c on pioglitazone 59.5 sitagliptin 59.9, canagliflozin 60.5 mmol mol^-1, P = 0.19). In total, 115 patients (25%) preferred pioglitazone, 158 patients (35%) sitagliptin and 175 patients (38%) canagliflozin. The drug preferred by individual patients was associated with a lower HbA1c (mean: 4.6; 95% CI: 3.9, 5.3) mmol mol^-1 lower versus nonpreferred) and fewer side effects (mean: 0.50; 95% CI: 0.35, 0.64) fewer side effects versus nonpreferred). Allocating therapy based on the individually preferred drugs, rather than allocating all patients the overall most preferred drug (canagliflozin), would result in more patients achieving the lowest HbA1c for them (70% versus 30%) and the fewest side effects (67% versus 50%). When precision approaches do not predict a clear optimal therapy for an individual, allowing patients to try potential suitable medications before they choose long-term therapy could be a practical alternative to optimizing treatment for type 2 diabetes.

A phenome-wide comparative analysis of genetic discordance between obesity and type 2 diabetes

Obesity and type 2 diabetes are causally related, yet there is considerable heterogeneity in the consequences of both conditions and the mechanisms of action are poorly defined. Here we show a genetic-driven approach defining two obesity profiles that convey highly concordant and discordant diabetogenic effects. We annotate and then compare association signals for these profiles across clinical and molecular phenotypic layers. Key differences are identified in a wide range of traits, including cardiovascular mortality, fat distribution, liver metabolism, blood pressure, specific lipid fractions and blood levels of proteins involved in extracellular matrix remodelling. We find marginal differences in abundance of Bacteroidetes and Firmicutes bacteria in the gut. Instrumental analyses reveal prominent causal roles for waist-to-hip ratio, blood pressure and cholesterol content of high-density lipoprotein particles in the development of diabetes in obesity. We prioritize 17 genes from the discordant signature that convey protection against type 2 diabetes in obesity, which may represent logical targets for precision medicine approaches.

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; I² = 97.28%), cardiovascular death (RR = 1.29; 95% CI 1.03-1.60; P < 0.0001; I² = 55.16%), myocardial infarction (RR = 1.32; 95% CI 1.09-1.59; P = 0.0037; I² = 97.14%), stroke (RR = 1.21; 95% CI 1.19-1.24; P < 0.0001; I² = 0.06%), and compound CVD outcomes (RR = 1.36; 95% CI 1.08-1.73; P = 0.01; I² = 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.

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 HbA_1c 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 HbA_1c 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 HbA_1c 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 HbA_1c 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 HbA_1c 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 HbA_1c 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.

Using Data to Improve the Management of Diabetes: The Tayside Experience

Tayside is a region in the East of Scotland and forms one of nine local government regions in the country. It is home to approximately 416,000 individuals who fall under the National Health Service (NHS) Tayside health board, which provides health care services to the population. In Tayside, Scotland, a comprehensive informatics network for diabetes care and research has been established for over 25 years. This has expanded more recently to a comprehensive Scotland-wide clinical care system, Scottish Care Information - Diabetes (SCI-Diabetes). This has enabled improved diabetes screening and integrated management of diabetic retinopathy, neuropathy, nephropathy, cardiovascular health, and other comorbidities. The regional health informatics network links all of these specialized services with comprehensive laboratory testing, prescribing records, general practitioner records, and hospitalization records. Not only do patients benefit from the seamless interconnectedness of these data, but also the Tayside bioresource has enabled considerable research opportunities and the creation of biobanks. In this article we describe how health informatics has been used to improve care of people with diabetes in Tayside and Scotland and, through anonymized data linkage, our understanding of the phenotypic and genotypic etiology of diabetes and associated complications and comorbidities.

Development of a treatment selection algorithm for SGLT2 and DPP-4 inhibitor therapies in people with type 2 diabetes: a retrospective cohort study

BACKGROUND

Current treatment guidelines do not provide recommendations to support the selection of treatment for most people with type 2 diabetes. We aimed to develop and validate an algorithm to allow selection of optimal treatment based on glycaemic response, weight change, and tolerability outcomes when choosing between SGLT2 inhibitor or DPP-4 inhibitor therapies.

METHODS

In this retrospective cohort study, we identified patients initiating SGLT2 and DPP-4 inhibitor therapies after Jan 1, 2013, from the UK Clinical Practice Research Datalink (CPRD). We excluded those who received SGLT2 or DPP-4 inhibitors as first-line treatment or insulin at the same time, had estimated glomerular filtration rate (eGFR) of less than 45 mL/min per 1·73 m², or did not have a valid baseline glycated haemoglobin (HbA_1c) measure (<53 or ≥120 mmol/mol). The primary efficacy outcome was the HbA_1c value reached 6 months after drug initiation, adjusted for baseline HbA_1c. Clinical features associated with differential HbA_1c outcome on the two therapies were identified in CPRD (n=26 877), and replicated in reanalysis of 14 clinical trials (n=10 414). An algorithm to predict individual-level differential HbA_1c outcome on the two therapies was developed in CPRD (derivation; n=14 069) and validated in head-to-head trials (n=2499) and CPRD (independent validation; n=9376). In CPRD, we further explored heterogeneity in 6-month weight change and treatment discontinuation.

FINDINGS

Among 10 253 patients initiating SGLT2 inhibitors and 16 624 patients initiating DPP-4 inhibitors in CPRD, baseline HbA_1c, age, BMI, eGFR, and alanine aminotransferase were associated with differential HbA_1c outcome with SGLT2 inhibitor and DPP-4 inhibitor therapies. The median age of participants was 62·0 years (IQR 55·0-70·0). 10 016 (37·3%) were women and 16 861 (62·7%) were men. An algorithm based on these five features identified a subgroup, representing around four in ten CPRD patients, with a 5 mmol/mol or greater observed benefit with SGLT2 inhibitors in all validation cohorts (CPRD 8·8 mmol/mol [95% CI 7·8-9·8]; CANTATA-D and CANTATA-D2 trials 5·8 mmol/mol [3·9-7·7]; BI1245.20 trial 6·6 mmol/mol [2·2-11·0]). In CPRD, predicted differential HbA_1c response with SGLT2 inhibitor and DPP-4 inhibitor therapies was not associated with weight change. Overall treatment discontinuation within 6 months was similar in patients predicted to have an HbA_1c benefit with SGLT2 inhibitors over DPP-4 inhibitors (median 15·2% [13·2-20·3] vs 14·4% [12·9-16·7]). A smaller subgroup predicted to have greater HbA_1c reduction with DPP-4 inhibitors were twice as likely to discontinue SGLT2 inhibitors than DPP-4 inhibitors (median 26·8% [23·4-31·0] vs 14·8% [12·9-16·8]).

INTERPRETATION

A validated treatment selection algorithm for SGLT2 inhibitor and DPP-4 inhibitor therapies can support decisions on optimal treatment for people with type 2 diabetes.

FUNDING

BHF-Turing Cardiovascular Data Science Award and the UK Medical Research Council.

Associations between parental type 2 diabetes risk and offspring birthweight and placental weight: a survival analysis using the Walker cohort

AIMS/HYPOTHESIS

Low birthweight (BW) is associated with the development of type 2 diabetes. Genome-wide analyses have identified a strong genetic component to this association, with many BW-associated loci also involved in glucose metabolism. We hypothesised that offspring BW and placental weight (PW) are correlated with parental type 2 diabetes risk, reflecting the inheritance of diabetes risk alleles that also influence fetal growth.

METHODS

The Walker cohort, a collection of birth records from Dundee, Scotland, from the 1950s and the 1960s was used to test this hypothesis by linking BW and PW measurements to parental health outcomes. Using data from SCI-Diabetes and the national death registry, we obtained health records for over 20,000 Walker parents. We performed Fine-Gray survival analyses of parental type 2 diabetes risk with competing risk of death, and Cox regression analyses of risk of death, independently in the maternal and paternal datasets, modelled by offspring BW and PW.

RESULTS

We found significant associations between increased paternal type 2 diabetes risk and reduced offspring BW (subdistribution hazard ratio [SHR] 0.92 [95% CI 0.87, 0.98]) and PW (SHR 0.87 [95% CI 0.81, 0.94]). The association of maternal type 2 diabetes risk with offspring BW or PW was not significant. Lower offspring BW was also associated with increased risk of death in both mothers (HR 0.91 [95% CI 0.89, 0.94]) and fathers (HR 0.95 [95% CI 0.92, 0.98]), and higher offspring PW was associated with increased maternal mortality risk (HR 1.08 [95% CI 1.04, 1.13]) when adjusted for BW.

CONCLUSIONS/INTERPRETATION

We identified associations between offspring BW and reduced paternal type 2 diabetes risk, most likely resulting from the independent effects of common type 2 diabetes susceptibility alleles on fetal growth, as described by the fetal insulin hypothesis. Moreover, we identified novel associations between offspring PW and reduced paternal type 2 diabetes risk, a relationship that might also be caused by the inheritance of diabetes predisposition variants. We found differing associations between offspring BW and PW and parental risk of death. These results provide novel epidemiological support for the use of offspring BW and PW as predictors for future risk of type 2 diabetes and death in mothers and fathers.

A roadmap to achieve pharmacological precision medicine in diabetes

Current pharmacological treatment of diabetes is largely algorithmic. Other than for cardiovascular disease or renal disease, where sodium-glucose cotransporter 2 inhibitors and/or glucagon-like peptide-1 receptor agonists are indicated, the choice of treatment is based upon overall risks of harm or side effect and cost, and not on probable benefit. Here we argue that a more precise approach to treatment choice is necessary to maximise benefit and minimise harm from existing diabetes therapies. We propose a roadmap to achieve precision medicine as standard of care, to discuss current progress in relation to monogenic diabetes and type 2 diabetes, and to determine what additional work is required. The first step is to identify robust and reliable genetic predictors of response, recognising that genotype is static over time and provides the skeleton upon which modifiers such as clinical phenotype and metabolic biomarkers can be overlaid. The second step is to identify these metabolic biomarkers (e.g. beta cell function, insulin sensitivity, BMI, liver fat, metabolite profile), which capture the metabolic state at the point of prescribing and may have a large impact on drug response. Third, we need to show that predictions that utilise these genetic and metabolic biomarkers improve therapeutic outcomes for patients, and fourth, that this is cost-effective. Finally, these biomarkers and prediction models need to be embedded in clinical care systems to enable effective and equitable clinical implementation. Whilst this roadmap is largely complete for monogenic diabetes, we still have considerable work to do to implement this for type 2 diabetes. Increasing collaborations, including with industry, and access to clinical trial data should enable progress to implementation of precision treatment in type 2 diabetes in the near future.

Cohort profile: the Scottish Diabetes Research Network national diabetes cohort - a population-based cohort of people with diabetes in Scotland

PURPOSE

The Scottish Diabetes Research Network (SDRN)-diabetes research platform was established to combine disparate electronic health record data into research-ready linked datasets for diabetes research in Scotland. The resultant cohort, 'The SDRN-National Diabetes Dataset (SDRN-NDS)', has many uses, for example, understanding healthcare burden and socioeconomic trends in disease incidence and prevalence, observational pharmacoepidemiology studies and building prediction tools to support clinical decision making.

PARTICIPANTS

We estimate that >99% of those diagnosed with diabetes nationwide are captured into the research platform. Between 2006 and mid-2020, the cohort comprised 472 648 people alive with diabetes at any point in whom there were 4 million person-years of follow-up. Of the cohort, 88.1% had type 2 diabetes, 8.8% type 1 diabetes and 3.1% had other types (eg, secondary diabetes). Data are captured from all key clinical encounters for diabetes-related care, including diabetes clinic, primary care and podiatry and comprise clinical history and measurements with linkage to blood results, microbiology, prescribed and dispensed drug and devices, retinopathy screening, outpatient, day case and inpatient episodes, birth outcomes, cancer registry, renal registry and causes of death.

FINDINGS TO DATE

There have been >50 publications using the SDRN-NDS. Examples of recent key findings include analysis of the incidence and relative risks for COVID-19 infection, drug safety of insulin glargine and SGLT2 inhibitors, life expectancy estimates, evaluation of the impact of flash monitors on glycaemic control and diabetic ketoacidosis and time trend analysis showing that diabetic ketoacidosis (DKA) remains a major cause of death under age 50 years. The findings have been used to guide national diabetes strategy and influence national and international guidelines.

FUTURE PLANS

The comprehensive SDRN-NDS will continue to be used in future studies of diabetes epidemiology in the Scottish population. It will continue to be updated at least annually, with new data sources linked as they become available.

Visit-to-visit glycated hemoglobin A1c variability in adults with type 2 diabetes: a systematic review and meta-analysis

BACKGROUND

Current practice uses the latest measure of glycated hemoglobin (HbAlc) to facilitate clinical decision-making. Studies have demonstrated that HbAlc variability links the risk of death and complications of diabetes. However, the role of HbAlc variability is unclear in clinical practice. This systematic review summarized the evidence of visit-to-visit HbAlc variability regarding different metrics in micro- and macro-vascular complications and death in people with type 2 diabetes.

METHODS

We searched PubMed, EMBASE (via OVID), and Cochrane Central Register (CENTRAL, via OVID) for studies investigating the association between HbAlc variability and adverse outcomes in patients with type 2 diabetes and performed random-effects meta-analysis stratified by HbAlc variability metrics in terms of standard deviation (SD), coefficient of variation (CV), and HbAlc variability score (HVS).

RESULTS

In people with type 2 diabetes, the highest quantile of all three HbAlc variability metrics (HbAlc-standard deviation [HbAlc-SD], HbAlc-coefficient of variance [HbAlc-CV], and HVS) is associated with increased risks of all-cause mortality, cardiovascular events, progression to chronic kidney disease, amputation, and peripheral neuropathy. For example, the hazard ratio of HbAlc-SD on all-cause mortality was l.89 with 95% confidence interval (95% CI) l.46-2.45 (HbAlc-CV l.47, 95% CI l.26-l.72; HVS l.67, 95% CI l.34-2.09).

CONCLUSIONS

High HbAlc variability leads to micro- and macro-vascular complications of type 2 diabetes and related death. People with type 2 diabetes and high HbAlc variability need additional attention and care for the potential adverse outcomes.

The impact of birthweight on subsequent phenotype of type 2 diabetes in later life

AIMS

It is well established that low birthweight is associated with subsequent risk of type 2 diabetes (T2DM). The aim of our study was to use a large birth cohort linked to a national diabetes registry to investigate how birthweight impacts the phenotype at diagnosis of T2DM and the subsequent rate of glycaemic deterioration.

METHODS

We linked the Walker Birth Cohort (48,000 births, 1952-1966, Tayside, Scotland) to the national diabetes registry in Scotland (SCI-Diabetes). Birthweight was adjusted for gestational age. Simple linear regression was performed to assess the impact of the adjusted birthweight on the diabetes phenotype at diagnosis. This was then built up into a multiple regression model to allow for the adjustment of confounding variables. A cox proportional hazards model was then used to evaluate the impact of birthweight on diabetes progression.

RESULTS

Lower birthweights were associated with a 293 day younger age of diagnosis of T2DM per 1 kg reduction in birthweight, p = 0.005; and a 1.29 kg/m² lower BMI at diagnosis per 1 kg reduction in birthweight, p < 0.001. There was no significant association of birthweight on diabetes progression.

CONCLUSION

For the first time, we have shown that a lower birthweight is associated with younger onset of T2DM, with those with lower birthweight also being slimmer at diagnosis. These results suggest that lower birthweight impacts on T2DM phenotype via reduced beta-cell function rather than insulin resistance.

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.

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.

Metabolic dysfunction-related liver disease as a risk factor for cancer

Objective

The aim of this study was to investigate the association between obesity, diabetes and metabolic related liver dysfunction and the incidence of cancer.

Design

This study was conducted with health record data available from the National Health Service in Tayside and Fife. Genetics of Diabetes Audit and Research Tayside, Scotland (GoDARTS), Scottish Health Research Register (SHARE) and Tayside and Fife diabetics, three Scottish cohorts of 13 695, 62 438 and 16 312 patients, respectively, were analysed in this study. Participants in GoDARTS were a volunteer sample, with half having type 2 diabetes mellitus(T2DM). SHARE was a volunteer sample. Tayside and Fife diabetics was a population-level cohort. Metabolic dysfunction-related liver disease (MDLD) was defined using alanine transaminase measurements, and individuals with alternative causes of liver disease (alcohol abuse, viruses, etc) were excluded from the analysis.

Results

MDLD associated with increased cancer incidence with a HR of 1.31 in a Cox proportional hazards model adjusted for sex, type 2 diabetes, body mass index(BMI), and smoking status (95% CI 1.27 to 1.35, p<0.0001). This was replicated in two further cohorts, and similar associations with cancer incidence were found for Fatty Liver Index (FLI), Fibrosis-4 Index (FIB-4) and non-alcoholic steatohepatitis (NASH). Homozygous carriers of the common non-alcoholic fatty liver disease (NAFLD) risk-variant PNPLA3 rs738409 had increased risk of cancer. (HR=1.27 (1.02 to 1.58), p=3.1×10⁻²). BMI was not independently associated with cancer incidence when MDLD was included as a covariate.

Conclusion

MDLD, FLI, FIB-4 and NASH associated with increased risk of cancer incidence and death. NAFLD may be a major component of the relationship between obesity and cancer incidence.

Prediction of Major Adverse Cardiovascular Events From Retinal, Clinical, and Genomic Data in Individuals With Type 2 Diabetes: A Population Cohort Study

OBJECTIVE

Improved identification of individuals with type 2 diabetes at high cardiovascular (CV) risk could help in selection of newer CV risk-reducing therapies. The aim of this study was to determine whether retinal vascular parameters, derived from retinal screening photographs, alone and in combination with a genome-wide polygenic risk score for coronary heart disease (CHD PRS) would have independent prognostic value over traditional CV risk assessment in patients without prior CV disease.

RESEARCH DESIGN AND METHODS

Patients in the Genetics of Diabetes Audit and Research Tayside Scotland (GoDARTS) study were linked to retinal photographs, prescriptions, and outcomes. Retinal photographs were analyzed using VAMPIRE (Vascular Assessment and Measurement Platform for Images of the Retina) software, a semiautomated artificial intelligence platform, to compute arterial and venous fractal dimension, tortuosity, and diameter. CHD PRS was derived from previously published data. Multivariable Cox regression was used to evaluate the association between retinal vascular parameters and major adverse CV events (MACE) at 10 years compared with the pooled cohort equations (PCE) risk score.

RESULTS

Among 5,152 individuals included in the study, a MACE occurred in 1,017 individuals. Reduced arterial fractal dimension and diameter and increased venous tortuosity each independently predicted MACE. A risk score combining these parameters significantly predicted MACE after adjustment for age, sex, PCE, and the CHD PRS (hazard ratio 1.11 per SD increase, 95% CI 1.04-1.18, P = 0.002) with similar accuracy to PCE (area under the curve [AUC] 0.663 vs. 0.658, P = 0.33). A model incorporating retinal parameters and PRS improved MACE prediction compared with PCE (AUC 0.686 vs. 0.658, P < 0.001).

CONCLUSIONS

Retinal parameters alone and in combination with genome-wide CHD PRS have independent and incremental prognostic value compared with traditional CV risk assessment in type 2 diabetes.