Previous case–control studies have reported aberrations of the gut microbiota in individuals with prediabetes. The primary objective of the present study was to explore the dynamics of the gut microbiota of individuals with prediabetes over 4 years with a secondary aim of relating microbiota dynamics to temporal changes of metabolic phenotypes. The study included 486 European patients with prediabetes. Gut microbiota profiling was conducted using shotgun metagenomic sequencing and the same bioinformatics pipelines at study baseline and after 4 years. The same phenotyping protocols and core laboratory analyses were applied at the two timepoints. Phenotyping included anthropometrics and measurement of fasting plasma glucose and insulin levels, mean plasma glucose and insulin under an oral glucose tolerance test (OGTT), 2-h plasma glucose after an OGTT, oral glucose insulin sensitivity index, Matsuda insulin sensitivity index, body mass index, waist circumference, and systolic and diastolic blood pressure. Measures of the dynamics of bacterial microbiota were related to concomitant changes in markers of host metabolism. Over 4 years, significant declines in richness were observed in gut bacterial and viral species and microbial pathways accompanied by significant changes in the relative abundance and the genetic composition of multiple bacterial species. Additionally, bacterial-viral interactions diminished over time. Despite the overall reduction in bacterial richness and microbial pathway richness, 80 dominant core bacterial species and 78 core microbial pathways were identified at both timepoints in 99

AIMS:Heterogeneity in the rate of β-cell loss in newly diagnosed type 1 diabetes patients is poorly understood and creates a barrier to designing and interpreting disease-modifying clinical trials. Integrative analyses of baseline multi-omics data obtained after the diagnosis of type 1 diabetes may provide mechanistic insight into the diverse rates of disease progression after type 1 diabetes diagnosis. METHODS:We collected samples in a pan-European consortium that enabled the concerted analysis of five different omics modalities in data from 97 newly diagnosed patients. In this study, we used Multi-Omics Factor Analysis to identify molecular signatures correlating with post-diagnosis decline in β-cell mass measured as fasting C-peptide. RESULTS:Two molecular signatures were significantly correlated with fasting C-peptide levels. One signature showed a correlation to neutrophil degranulation, cytokine signalling, lymphoid and non-lymphoid cell interactions and G-protein coupled receptor signalling events that were inversely associated with a rapid decline in β-cell function. The second signature was related to translation and viral infection was inversely associated with change in β-cell function. In addition, the immunomics data revealed a Natural Killer cell signature associated with rapid β-cell decline. CONCLUSIONS:Features that differ between individuals with slow and rapid decline in β-cell mass could be valuable in staging and prediction of the rate of disease progression and thus enable smarter (shorter and smaller) trial designs for disease modifying therapies as well as offering biomarkers of therapeutic effect.

AIMS/HYPOTHESIS:There is a growing need for markers that could help indicate the decline in beta cell function and recognise the need and efficacy of intervention in type 1 diabetes. Measurements of suitably selected serum markers could potentially provide a non-invasive and easily applicable solution to this challenge. Accordingly, we evaluated a broad panel of proteins previously associated with type 1 diabetes in serum from newly diagnosed individuals during the first year from diagnosis. To uncover associations with beta cell function, comparisons were made between these targeted proteomics measurements and changes in fasting C-peptide levels. To further distinguish proteins linked with the disease status, comparisons were made with measurements of the protein targets in age- and sex-matched autoantibody-negative unaffected family members (UFMs). METHODS:Selected reaction monitoring (SRM) mass spectrometry analyses of serum, targeting 85 type 1 diabetes-associated proteins, were made. Sera from individuals diagnosed under 18 years (n=86) were drawn within 6 weeks of diagnosis and at 3, 6 and 12 months afterwards (288 samples in total). The SRM data were compared with fasting C-peptide/glucose data, which was interpreted as a measure of beta cell function. The protein data were further compared with cross-sectional SRM measurements from UFMs (n=194). RESULTS:Eleven proteins had statistically significant associations with fasting C-peptide/glucose. Of these, apolipoprotein L1 and glutathione peroxidase 3 (GPX3) displayed the strongest positive and inverse associations, respectively. Changes in GPX3 levels during the first year after diagnosis indicated future fasting C-peptide/glucose levels. In addition, differences in the levels of 13 proteins were observed between the individuals with type 1 diabetes and the matched UFMs. These included GPX3, transthyretin, prothrombin, apolipoprotein C1 and members of the IGF family. CONCLUSIONS/INTERPRETATION:The association of several targeted proteins with fasting C-peptide/glucose levels in the first year after diagnosis suggests their connection with the underlying changes accompanying alterations in beta cell function in type 1 diabetes. Moreover, the direction of change in GPX3 during the first year was indicative of subsequent fasting C-peptide/glucose levels, and supports further investigation of this and other serum protein measurements in future studies of beta cell function in type 1 diabetes.

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.

BACKGROUND:Type 1 diabetes is a complex heterogenous autoimmune disease without therapeutic interventions available to prevent or reverse the disease. This study aimed to identify transcriptional changes associated with the disease progression in patients with recent-onset type 1 diabetes. METHODS:Whole-blood samples were collected as part of the INNODIA study at baseline and 12 months after diagnosis of type 1 diabetes. We used linear mixed-effects modelling on RNA-seq data to identify genes associated with age, sex, or disease progression. Cell-type proportions were estimated from the RNA-seq data using computational deconvolution. Associations to clinical variables were estimated using Pearson's or point-biserial correlation for continuous and dichotomous variables, respectively, using only complete pairs of observations. FINDINGS:We found that genes and pathways related to innate immunity were downregulated during the first year after diagnosis. Significant associations of the gene expression changes were found with ZnT8A autoantibody positivity. Rate of change in the expression of 16 genes between baseline and 12 months was found to predict the decline in C-peptide at 24 months. Interestingly and consistent with earlier reports, increased B cell levels and decreased neutrophil levels were associated with the rapid progression. INTERPRETATION:There is considerable individual variation in the rate of progression from appearance of type 1 diabetes-specific autoantibodies to clinical disease. Patient stratification and prediction of disease progression can help in developing more personalised therapeutic strategies for different disease endotypes. FUNDING:A full list of funding bodies can be found under Acknowledgments.

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.

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.

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.

Background: Development of disease-modifying therapy for type 1 diabetes (T1D) is hampered by the limited understanding of pathogenesis, heterogeneity, lack of disease biomarkers and stratifiers. Using the INNODIA consortium pan-European infrastructure to collect prospective clinical data from newly diagnosed subjects combined with multi-omics data, we mined integrated datasets using deep learning to identify novel relationships that could transform the disease monitoring landscape. Method: Samples collected <6 weeks from diagnosis from 100 islet autoantibody positive participants were prioritised for omics analysis and the data integrated with clinical parameters and C-peptide (CP) from mixed meal tolerance tests (at 3, 6 and 12 months). SNP array data, plasma metabolomics, lipidomics and miRNAs, serum proteomics, whole blood RNA-sequencing and multiFACS immunomics from PBMCs were collected. We applied a variational autoencoder (VAE) to simultaneously integrate the data and learn CP progression over 12 months. The latent (lower dimensional) representation of the VAE integrated data was clustered to identify clinically meaningful subgroups. Results: Clustering identified five age-independent subgroups with significantly different progression patterns (p<0.05). Cluster 1 and 3 associated with lower fasted and stimulated CP. Cluster 2 and 5 associated with higher fasted CP, while cluster 4 most strongly associated with higher stimulated CP. Cluster 1 had the lowest stimulated CP at baseline, cluster 5 the lowest fasted CP and cluster 4 the highest fasted and stimulated CP at baseline. Investigation of the most predictive signals from the omics data is ongoing. Discussion: The preliminary results suggest that the VAE can learn meaningful progression patterns from integrated clinical and multi-omics data, which holds great promise to uncover important data structures not readily discovered when analyzing single data types. Disclosure C. Brorsson: Employee; Spouse/Partner; Novo Nordisk A/S. J. Almagro armenteros: None. G. Mazzoni: Employee; Self; Novo Nordisk A/S. S. Kaur: None. A. M. Schulte: None. C. Mathieu: Advisory Panel; Self; Novo Nordisk, Sanofi, Merck Sharp and Dohme Ltd., Eli Lilly and Company, Novartis, AstraZeneca, Boehringer Ingelheim, Roche, Medtronic, ActoBio Therapeutics, Pfizer, Insulet and Zealand Pharma, Research Support; Self; Medtronic, Novo Nordisk, Sanofi and ActoBio Therapeutics, Speaker’s Bureau; Self; Novo Nordisk, Sanofi, Eli Lilly and Company, Boehringer Ingelheim, AstraZeneca and Novartis. S. Brunak: Board Member; Self; Intomics A/S, Proscion A/S, Stock/Shareholder; Self; Hoba Therapeutics Aps, Novo Nordisk A/S. On behalf of the innodia consortium: n/a. Funding INNODIA (115797)

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.

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

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

The study aimed to explore the role of ERBB3 in type 1 diabetes (T1D). We examined whether genetic variation of ERBB3 (rs2292239) affects residual β-cell function in T1D cases. Furthermore, we examined the expression of ERBB3 in human islets, the effect of ERBB3 knockdown on apoptosis in insulin-producing INS-1E cells and the genetic and regulatory architecture of the ERBB3 locus to provide insights to how rs2292239 may confer disease susceptibility. rs2292239 strongly correlated with residual β-cell function and metabolic control in children with T1D. ERBB3 locus associated lncRNA (NONHSAG011351) was found to be expressed in human islets. ERBB3 was expressed and down-regulated by pro-inflammatory cytokines in human islets and INS-1E cells; knockdown of ERBB3 in INS-1E cells decreased basal and cytokine-induced apoptosis. Our data suggests an important functional role of ERBB3 and its potential regulators in the β-cells and may constitute novel targets to prevent β-cell destruction in T1D.

Genome-wide association studies (GWAS) have identified over 40 type 1 diabetes risk loci. The clinical impact of these loci on β-cell function during disease progression is unknown. We aimed at testing whether a genetic risk score could predict glycemic control and residual β-cell function in type 1 diabetes (T1D). As gene expression may represent an intermediate phenotype between genetic variation and disease, we hypothesized that genes within T1D loci which are expressed in islets and transcriptionally regulated by proinflammatory cytokines would be the best predictors of disease progression. Two-thirds of 46 GWAS candidate genes examined were expressed in human islets, and 11 of these significantly changed expression levels following exposure to proinflammatory cytokines (IL-1β + IFNγ + TNFα) for 48 h. Using the GWAS single nucleotide polymorphisms (SNPs) from each locus, we constructed a genetic risk score based on the cumulative number of risk alleles carried in children with newly diagnosed T1D. With each additional risk allele carried, HbA1c levels increased significantly within first year after diagnosis. Network and gene ontology (GO) analyses revealed that several of the 11 candidate genes have overlapping biological functions and interact in a common network. Our results may help predict disease progression in newly diagnosed children with T1D which can be exploited for optimizing treatment.