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.

Discovery of drug-omics associations in type 2 diabetes with generative deep-learning models.

The application of multiple omics technologies in biomedical cohorts has the potential to reveal patient-level disease characteristics and individualized response to treatment. However, the scale and heterogeneous nature of multi-modal data makes integration and inference a non-trivial task. We developed a deep-learning-based framework, multi-omics variational autoencoders (MOVE), to integrate such data and applied it to a cohort of 789 people with newly diagnosed type 2 diabetes with deep multi-omics phenotyping from the DIRECT consortium. Using in silico perturbations, we identified drug-omics associations across the multi-modal datasets for the 20 most prevalent drugs given to people with type 2 diabetes with substantially higher sensitivity than univariate statistical tests. From these, we among others, identified novel associations between metformin and the gut microbiota as well as opposite molecular responses for the two statins, simvastatin and atorvastatin. We used the associations to quantify drug-drug similarities, assess the degree of polypharmacy and conclude that drug effects are distributed across the multi-omics modalities.

Four groups of type 2 diabetes contribute to the etiological and clinical heterogeneity in newly diagnosed individuals: An IMI DIRECT study.

The presentation and underlying pathophysiology of type 2 diabetes (T2D) is complex and heterogeneous. Recent studies attempted to stratify T2D into distinct subgroups using data-driven approaches, but their clinical utility may be limited if categorical representations of complex phenotypes are suboptimal. We apply a soft-clustering (archetype) method to characterize newly diagnosed T2D based on 32 clinical variables. We assign quantitative clustering scores for individuals and investigate the associations with glycemic deterioration, genetic risk scores, circulating omics biomarkers, and phenotypic stability over 36 months. Four archetype profiles represent dysfunction patterns across combinations of T2D etiological processes and correlate with multiple circulating biomarkers. One archetype associated with obesity, insulin resistance, dyslipidemia, and impaired ß cell glucose sensitivity corresponds with the fastest disease progression and highest demand for anti-diabetic treatment. We demonstrate that clinical heterogeneity in T2D can be mapped to heterogeneity in individual etiological processes, providing a potential route to personalized treatments.

SteC and the intracellular Salmonella-induced F-actin meshwork.

Salmonella enterica serovars infect a broad range of mammalian hosts including humans, causing both gastrointestinal and systemic diseases. Following uptake into host cells, bacteria replicate within vacuoles (Salmonella-containing vacuoles; SCVs). Clusters of SCVs are frequently associated with a meshwork of F-actin. This meshwork is dependent on the Salmonella pathogenicity island 2 encoded type III secretion system and its effector SteC. SteC contains a region with weak similarity to conserved subdomains of eukaryotic kinases and has kinase activity that is required for the formation of the F-actin meshwork. Several substrates of SteC have been identified. In this mini-review, we attempt to integrate these findings and propose a more unified model to explain SCV-associated F-actin: SteC (i) phosphorylates the actin sequestering protein Hsp27, which increases the local G-actin concentration (ii) binds to and phosphorylates formin family FMNL proteins, which enables actin polymerisation and (iii) phosphorylates MEK, resulting in activation of the MEK/ERK/MLCK/Myosin II pathway, leading to F-actin bundling. We also consider the possible physiological functions of SCV-associated F-actin and similar structures produced by other intracellular bacterial pathogens.

Genome-Wide Association Analysis of Pancreatic Beta-Cell Glucose Sensitivity.

CONTEXT: Pancreatic beta-cell glucose sensitivity is the slope of the plasma glucose-insulin secretion relationship and is a key predictor of deteriorating glucose tolerance and development of type 2 diabetes. However, there are no large-scale studies looking at the genetic determinants of beta-cell glucose sensitivity. OBJECTIVE: To understand the genetic determinants of pancreatic beta-cell glucose sensitivity using genome-wide meta-analysis and candidate gene studies. DESIGN: We performed a genome-wide meta-analysis for beta-cell glucose sensitivity in subjects with type 2 diabetes and nondiabetic subjects from 6 independent cohorts (n = 5706). Beta-cell glucose sensitivity was calculated from mixed meal and oral glucose tolerance tests, and its associations between known glycemia-related single nucleotide polymorphisms (SNPs) and genome-wide association study (GWAS) SNPs were estimated using linear regression models. RESULTS: Beta-cell glucose sensitivity was moderately heritable (h2 ranged from 34% to 55%) using SNP and family-based analyses. GWAS meta-analysis identified multiple correlated SNPs in the CDKAL1 gene and GIPR-QPCTL gene loci that reached genome-wide significance, with SNP rs2238691 in GIPR-QPCTL (P value = 2.64 × 10-9) and rs9368219 in the CDKAL1 (P value = 3.15 × 10-9) showing the strongest association with beta-cell glucose sensitivity. These loci surpassed genome-wide significance when the GWAS meta-analysis was repeated after exclusion of the diabetic subjects. After correction for multiple testing, glycemia-associated SNPs in or near the HHEX and IGF2B2 loci were also associated with beta-cell glucose sensitivity. CONCLUSION: We show that, variation at the GIPR-QPCTL and CDKAL1 loci are key determinants of pancreatic beta-cell glucose sensitivity.

Processes Underlying Glycemic Deterioration in Type 2 Diabetes: An IMI DIRECT Study.

OBJECTIVE: We investigated the processes underlying glycemic deterioration in type 2 diabetes (T2D). RESEARCH DESIGN AND METHODS: A total of 732 recently diagnosed patients with T2D from the Innovative Medicines Initiative Diabetes Research on Patient Stratification (IMI DIRECT) study were extensively phenotyped over 3 years, including measures of insulin sensitivity (OGIS), ß-cell glucose sensitivity (GS), and insulin clearance (CLIm) from mixed meal tests, liver enzymes, lipid profiles, and baseline regional fat from MRI. The associations between the longitudinal metabolic patterns and HbA1c deterioration, adjusted for changes in BMI and in diabetes medications, were assessed via stepwise multivariable linear and logistic regression. RESULTS: Faster HbA1c progression was independently associated with faster deterioration of OGIS and GS and increasing CLIm; visceral or liver fat, HDL-cholesterol, and triglycerides had further independent, though weaker, roles (R 2 = 0.38). A subgroup of patients with a markedly higher progression rate (fast progressors) was clearly distinguishable considering these variables only (discrimination capacity from area under the receiver operating characteristic = 0.94). The proportion of fast progressors was reduced from 56% to 8-10% in subgroups in which only one trait among OGIS, GS, and CLIm was relatively stable (odds ratios 0.07-0.09). T2D polygenic risk score and baseline pancreatic fat, glucagon-like peptide 1, glucagon, diet, and physical activity did not show an independent role. CONCLUSIONS: Deteriorating insulin sensitivity and ß-cell function, increasing insulin clearance, high visceral or liver fat, and worsening of the lipid profile are the crucial factors mediating glycemic deterioration of patients with T2D in the initial phase of the disease. Stabilization of a single trait among insulin sensitivity, ß-cell function, and insulin clearance may be relevant to prevent progression.

Whole blood co-expression modules associate with metabolic traits and type 2 diabetes: an IMI-DIRECT study.

BACKGROUND: The rising prevalence of type 2 diabetes (T2D) poses a major global challenge. It remains unresolved to what extent transcriptomic signatures of metabolic dysregulation and T2D can be observed in easily accessible tissues such as blood. Additionally, large-scale human studies are required to further our understanding of the putative inflammatory component of insulin resistance and T2D. Here we used transcriptomics data from individuals with (n = 789) and without (n = 2127) T2D from the IMI-DIRECT cohorts to describe the co-expression structure of whole blood that mainly reflects processes and cell types of the immune system, and how it relates to metabolically relevant clinical traits and T2D. METHODS: Clusters of co-expressed genes were identified in the non-diabetic IMI-DIRECT cohort and evaluated with regard to stability, as well as preservation and rewiring in the cohort of individuals with T2D. We performed functional and immune cell signature enrichment analyses, and a genome-wide association study to describe the genetic regulation of the modules. Phenotypic and trans-omics associations of the transcriptomic modules were investigated across both IMI-DIRECT cohorts. RESULTS: We identified 55 whole blood co-expression modules, some of which clustered in larger super-modules. We identified a large number of associations between these transcriptomic modules and measures of insulin action and glucose tolerance. Some of the metabolically linked modules reflect neutrophil-lymphocyte ratio in blood while others are independent of white blood cell estimates, including a module of genes encoding neutrophil granule proteins with antibacterial properties for which the strongest associations with clinical traits and T2D status were observed. Through the integration of genetic and multi-omics data, we provide a holistic view of the regulation and molecular context of whole blood transcriptomic modules. We furthermore identified an overlap between genetic signals for T2D and co-expression modules involved in type II interferon signaling. CONCLUSIONS: Our results offer a large-scale map of whole blood transcriptomic modules in the context of metabolic disease and point to novel biological candidates for future studies related to T2D.

Post-load glucose subgroups and associated metabolic traits in individuals with type 2 diabetes: An IMI-DIRECT study.

AIM: Subclasses of different glycaemic disturbances could explain the variation in characteristics of individuals with type 2 diabetes (T2D). We aimed to examine the association between subgroups based on their glucose curves during a five-point mixed-meal tolerance test (MMT) and metabolic traits at baseline and glycaemic deterioration in individuals with T2D. METHODS: The study included 787 individuals with newly diagnosed T2D from the Diabetes Research on Patient Stratification (IMI-DIRECT) Study. Latent class trajectory analysis (LCTA) was used to identify distinct glucose curve subgroups during a five-point MMT. Using general linear models, these subgroups were associated with metabolic traits at baseline and after 18 months of follow up, adjusted for potential confounders. RESULTS: At baseline, we identified three glucose curve subgroups, labelled in order of increasing glucose peak levels as subgroup 1-3. Individuals in subgroup 2 and 3 were more likely to have higher levels of HbA1c, triglycerides and BMI at baseline, compared to those in subgroup 1. At 18 months (n = 651), the beta coefficients (95% CI) for change in HbA1c (mmol/mol) increased across subgroups with 0.37 (-0.18-1.92) for subgroup 2 and 1.88 (-0.08-3.85) for subgroup 3, relative to subgroup 1. The same trend was observed for change in levels of triglycerides and fasting glucose. CONCLUSIONS: Different glycaemic profiles with different metabolic traits and different degrees of subsequent glycaemic deterioration can be identified using data from a frequently sampled mixed-meal tolerance test in individuals with T2D. Subgroups with the highest peaks had greater metabolic risk.

Predicting and elucidating the etiology of fatty liver disease: A machine learning modeling and validation study in the IMI DIRECT cohorts.

BACKGROUND: Non-alcoholic fatty liver disease (NAFLD) is highly prevalent and causes serious health complications in individuals with and without type 2 diabetes (T2D). Early diagnosis of NAFLD is important, as this can help prevent irreversible damage to the liver and, ultimately, hepatocellular carcinomas. We sought to expand etiological understanding and develop a diagnostic tool for NAFLD using machine learning. METHODS AND FINDINGS: We utilized the baseline data from IMI DIRECT, a multicenter prospective cohort study of 3,029 European-ancestry adults recently diagnosed with T2D (n = 795) or at high risk of developing the disease (n = 2,234). Multi-omics (genetic, transcriptomic, proteomic, and metabolomic) and clinical (liver enzymes and other serological biomarkers, anthropometry, measures of beta-cell function, insulin sensitivity, and lifestyle) data comprised the key input variables. The models were trained on MRI-image-derived liver fat content (<5% or ≥5%) available for 1,514 participants. We applied LASSO (least absolute shrinkage and selection operator) to select features from the different layers of omics data and random forest analysis to develop the models. The prediction models included clinical and omics variables separately or in combination. A model including all omics and clinical variables yielded a cross-validated receiver operating characteristic area under the curve (ROCAUC) of 0.84 (95% CI 0.82, 0.86; p < 0.001), which compared with a ROCAUC of 0.82 (95% CI 0.81, 0.83; p < 0.001) for a model including 9 clinically accessible variables. The IMI DIRECT prediction models outperformed existing noninvasive NAFLD prediction tools. One limitation is that these analyses were performed in adults of European ancestry residing in northern Europe, and it is unknown how well these findings will translate to people of other ancestries and exposed to environmental risk factors that differ from those of the present cohort. Another key limitation of this study is that the prediction was done on a binary outcome of liver fat quantity (<5% or ≥5%) rather than a continuous one. CONCLUSIONS: In this study, we developed several models with different combinations of clinical and omics data and identified biological features that appear to be associated with liver fat accumulation. In general, the clinical variables showed better prediction ability than the complex omics variables. However, the combination of omics and clinical variables yielded the highest accuracy. We have incorporated the developed clinical models into a web interface (see: https://www.predictliverfat.org/) and made it available to the community. TRIAL REGISTRATION: ClinicalTrials.gov NCT03814915.

Discovery of biomarkers for glycaemic deterioration before and after the onset of type 2 diabetes: descriptive characteristics of the epidemiological studies within the IMI DIRECT Consortium.

AIMS/HYPOTHESIS: Here, we describe the characteristics of the Innovative Medicines Initiative (IMI) Diabetes Research on Patient Stratification (DIRECT) epidemiological cohorts at baseline and follow-up examinations (18, 36 and 48 months of follow-up). METHODS: From a sampling frame of 24,682 adults of European ancestry enrolled in population-based cohorts across Europe, participants at varying risk of glycaemic deterioration were identified using a risk prediction algorithm (based on age, BMI, waist circumference, use of antihypertensive medication, smoking status and parental history of type 2 diabetes) and enrolled into a prospective cohort study (n = 2127) (cohort 1, prediabetes risk). We also recruited people from clinical registries with type 2 diabetes diagnosed 6-24 months previously (n = 789) into a second cohort study (cohort 2, diabetes). Follow-up examinations took place at ~18 months (both cohorts) and at ~48 months (cohort 1) or ~36 months (cohort 2) after baseline examinations. The cohorts were studied in parallel using matched protocols across seven clinical centres in northern Europe. RESULTS: Using ADA 2011 glycaemic categories, 33% (n = 693) of cohort 1 (prediabetes risk) had normal glucose regulation and 67% (n = 1419) had impaired glucose regulation. Seventy-six per cent of participants in cohort 1 was male. Cohort 1 participants had the following characteristics (mean ± SD) at baseline: age 62 (6.2) years; BMI 27.9 (4.0) kg/m2; fasting glucose 5.7 (0.6) mmol/l; 2 h glucose 5.9 (1.6) mmol/l. At the final follow-up examination the participants' clinical characteristics were as follows: fasting glucose 6.0 (0.6) mmol/l; 2 h OGTT glucose 6.5 (2.0) mmol/l. In cohort 2 (diabetes), 66% (n = 517) were treated by lifestyle modification and 34% (n = 272) were treated with metformin plus lifestyle modification at enrolment. Fifty-eight per cent of participants in cohort 2 was male. Cohort 2 participants had the following characteristics at baseline: age 62 (8.1) years; BMI 30.5 (5.0) kg/m2; fasting glucose 7.2 (1.4) mmol/l; 2 h glucose 8.6 (2.8) mmol/l. At the final follow-up examination, the participants' clinical characteristics were as follows: fasting glucose 7.9 (2.0) mmol/l; 2 h mixed-meal tolerance test glucose 9.9 (3.4) mmol/l. CONCLUSIONS/INTERPRETATION: The IMI DIRECT cohorts are intensely characterised, with a wide-variety of metabolically relevant measures assessed prospectively. We anticipate that the cohorts, made available through managed access, will provide a powerful resource for biomarker discovery, multivariate aetiological analyses and reclassification of patients for the prevention and treatment of type 2 diabetes.

Sustained influence of metformin therapy on circulating glucagon-like peptide-1 levels in individuals with and without type 2 diabetes.

AIMS: To investigate, in the Carotid Atherosclerosis: Metformin for Insulin Resistance (CAMERA) trial (NCT00723307), whether the influence of metformin on the glucagon-like peptide (GLP)-1 axis in individuals with and without type 2 diabetes (T2DM) is sustained and related to changes in glycaemia or weight, and to investigate basal and post-meal GLP-1 levels in patients with T2DM in the cross-sectional Diabetes Research on Patient Stratification (DIRECT) study. MATERIALS AND METHODS: CAMERA was a double-blind randomized placebo-controlled trial of metformin in 173 participants without diabetes. Using 6-monthly fasted total GLP-1 levels over 18 months, we evaluated metformin's effect on total GLP-1 with repeated-measures analysis and analysis of covariance. In the DIRECT study, we examined active and total fasting and 60-minute post-meal GLP-1 levels in 775 people recently diagnosed with T2DM treated with metformin or diet, using Student's t-tests and linear regression. RESULTS: In CAMERA, metformin increased total GLP-1 at 6 (+20.7%, 95% confidence interval [CI] 4.7-39.0), 12 (+26.7%, 95% CI 10.3-45.6) and 18 months (+18.7%, 95% CI 3.8-35.7), an overall increase of 23.4% (95% CI 11.2-36.9; P < .0001) vs placebo. Adjustment for changes in glycaemia and adiposity, individually or combined, did not attenuate this effect. In the DIRECT study, metformin was associated with higher fasting active (39.1%, 95% CI 21.3-56.4) and total GLP-1 (14.1%, 95% CI 1.2-25.9) but not post-meal incremental GLP-1. These changes were independent of potential confounders including age, sex, adiposity and glycated haemoglobin. CONCLUSIONS: In people without diabetes, metformin increases total GLP-1 in a sustained manner and independently of changes in weight or glycaemia. Metformin-treated patients with T2DM also have higher fasted GLP-1 levels, independently of weight and glycaemia.

Age-related mitochondrial DNA depletion and the impact on pancreatic Beta cell function.

Type 2 diabetes is characterised by an age-related decline in insulin secretion. We previously identified a 50% age-related decline in mitochondrial DNA (mtDNA) copy number in isolated human islets. The purpose of this study was to mimic this degree of mtDNA depletion in MIN6 cells to determine whether there is a direct impact on insulin secretion. Transcriptional silencing of mitochondrial transcription factor A, TFAM, decreased mtDNA levels by 40% in MIN6 cells. This level of mtDNA depletion significantly decreased mtDNA gene transcription and translation, resulting in reduced mitochondrial respiratory capacity and ATP production. Glucose-stimulated insulin secretion was impaired following partial mtDNA depletion, but was normalised following treatment with glibenclamide. This confirms that the deficit in the insulin secretory pathway precedes K+ channel closure, indicating that the impact of mtDNA depletion is at the level of mitochondrial respiration. In conclusion, partial mtDNA depletion to a degree comparable to that seen in aged human islets impaired mitochondrial function and directly decreased insulin secretion. Using our model of partial mtDNA depletion following targeted gene silencing of TFAM, we have managed to mimic the degree of mtDNA depletion observed in aged human islets, and have shown how this correlates with impaired insulin secretion. We therefore predict that the age-related mtDNA depletion in human islets is not simply a biomarker of the aging process, but will contribute to the age-related risk of type 2 diabetes.

Discovery of biomarkers for glycaemic deterioration before and after the onset of type 2 diabetes: rationale and design of the epidemiological studies within the IMI DIRECT Consortium.

AIMS/HYPOTHESIS: The DIRECT (Diabetes Research on Patient Stratification) Study is part of a European Union Framework 7 Innovative Medicines Initiative project, a joint undertaking between four industry and 21 academic partners throughout Europe. The Consortium aims to discover and validate biomarkers that: (1) predict the rate of glycaemic deterioration before and after type 2 diabetes onset; (2) predict the response to diabetes therapies; and (3) help stratify type 2 diabetes into clearly definable disease subclasses that can be treated more effectively than without stratification. This paper describes two new prospective cohort studies conducted as part of DIRECT. METHODS: Prediabetic participants (target sample size 2,200-2,700) and patients with newly diagnosed type 2 diabetes (target sample size ~1,000) are undergoing detailed metabolic phenotyping at baseline and 18 months and 36 months later. Abdominal, pancreatic and liver fat is assessed using MRI. Insulin secretion and action are assessed using frequently sampled OGTTs in non-diabetic participants, and frequently sampled mixed-meal tolerance tests in patients with type 2 diabetes. Biosamples include venous blood, faeces, urine and nail clippings, which, among other biochemical analyses, will be characterised at genetic, transcriptomic, metabolomic, proteomic and metagenomic levels. Lifestyle is assessed using high-resolution triaxial accelerometry, 24 h diet record, and food habit questionnaires. CONCLUSIONS/INTERPRETATION: DIRECT will yield an unprecedented array of biomaterials and data. This resource, available through managed access to scientists within and outside the Consortium, will facilitate the development of new treatments and therapeutic strategies for the prevention and management of type 2 diabetes.