An interactive atlas of genomic, proteomic, and metabolomic biomarkers promotes the potential of proteins to predict complex diseases.

Multiomics analyses have identified multiple potential biomarkers of the incidence and prevalence of complex diseases. However, it is not known which type of biomarker is optimal for clinical purposes. Here, we make a systematic comparison of 90 million genetic variants, 1453 proteins, and 325 metabolites from 500,000 individuals with complex diseases from the UK Biobank. A machine learning pipeline consisting of data cleaning, data imputation, feature selection, and model training using cross-validation and comparison of the results on holdout test sets showed that proteins were most predictive, followed by metabolites, and genetic variants. Only five proteins per disease resulted in median (min-max) areas under the receiver operating characteristic curves for incidence of 0.79 (0.65-0.86) and 0.84 (0.70-0.91) for prevalence. In summary, our work suggests the potential of predicting complex diseases based on a limited number of proteins. We provide an interactive atlas (macd.shinyapps.io/ShinyApp/) to find genomic, proteomic, or metabolomic biomarkers for different complex diseases.

Transcript and protein signatures derived from shared molecular interactions across cancers are associated with mortality.

BACKGROUND: Characterization of shared cancer mechanisms have been proposed to improve therapy strategies and prognosis. Here, we aimed to identify shared cell-cell interactions (CCIs) within the tumor microenvironment across multiple solid cancers and assess their association with cancer mortality. METHODS: CCIs of each cancer were identified by NicheNet analysis of single-cell RNA sequencing data from breast, colon, liver, lung, and ovarian cancers. These CCIs were used to construct a shared multi-cellular tumor model (shared-MCTM) representing common CCIs across cancers. A gene signature was identified from the shared-MCTM and tested on the mRNA and protein level in two large independent cohorts: The Cancer Genome Atlas (TCGA, 9185 tumor samples and 727 controls across 22 cancers) and UK biobank (UKBB, 10,384 cancer patients and 5063 controls with proteomics data across 17 cancers). Cox proportional hazards models were used to evaluate the association of the signature with 10-year all-cause mortality, including sex-specific analysis. RESULTS: A shared-MCTM was derived from five individual cancers. A shared gene signature was extracted from this shared-MCTM and the most prominent regulatory cell type, matrix cancer-associated fibroblast (mCAF). The signature exhibited significant expression changes in multiple cancers compared to controls at both mRNA and protein levels in two independent cohorts. Importantly, it was significantly associated with mortality in cancer patients in both cohorts. The highest hazard ratios were observed for brain cancer in TCGA (HR [95%CI] = 6.90[4.64-10.25]) and ovarian cancer in UKBB (5.53[2.08-8.80]). Sex-specific analysis revealed distinct risks, with a higher mortality risk associated with the protein signature score in males (2.41[1.97-2.96]) compared to females (1.84[1.44-2.37]). CONCLUSION: We identified a gene signature from a comprehensive shared-MCTM representing common CCIs across different cancers and revealed the regulatory role of mCAF in the tumor microenvironment. The pathogenic relevance of the gene signature was supported by differential expression and association with mortality on both mRNA and protein levels in two independent cohorts.

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.

Incidence of LADA and Type 2 Diabetes in Relation to Tobacco Use and Genetic Susceptibility to Type 2 Diabetes and Related Traits: Findings From a Swedish Case-Control Study and the Norwegian HUNT Study.

OBJECTIVE: Smoking and Swedish smokeless tobacco (snus) are associated with latent autoimmune diabetes in adults (LADA) and type 2 diabetes (T2D). Our aim was to investigate whether genetic susceptibility to T2D, insulin resistance (IR), and insulin secretion (IS) aggravate these associations. RESEARCH DESIGN AND METHODS: We used data from two population-based Scandinavian studies with case subjects with LADA (n = 839) and T2D (n = 5,771), matched control subjects (n = 3,068), and 1,696,503 person-years at risk. Pooled, multivariate relative risks (RR) with 95% CI were estimated for smoking/genetic risk scores (T2D-GRS, IS-GRS, and IR-GRS), and ORs for snus or tobacco/GRS (case-control data). We estimated additive (proportion attributable to interaction [AP]) and multiplicative interaction between tobacco use and GRS. RESULTS: The RR of LADA was elevated in high IR-GRS heavy smokers (≥15 pack-years; RR 2.01 [CI 1.30, 3.10]) and tobacco users (≥15 box/pack-years; RR 2.59 [CI 1.54, 4.35]) compared with low IR-GRS individuals without heavy use, with evidence of additive (AP 0.67 [CI 0.46, 0.89]; AP 0.52 [CI 0.21, 0.83]) and multiplicative (P = 0.003; P = 0.034) interaction. In heavy users, there was additive interaction between T2D-GRS and smoking, snus, and total tobacco use. The excess risk conferred by tobacco use did not differ across GRS categories in T2D. CONCLUSIONS: Tobacco use may confer a higher risk of LADA in individuals with genetic susceptibility to T2D and insulin resistance, whereas genetic susceptibility does not seem to influence the increased T2D incidence associated with tobacco use.

Islet Gene View-a tool to facilitate islet research.

Characterization of gene expression in pancreatic islets and its alteration in type 2 diabetes (T2D) are vital in understanding islet function and T2D pathogenesis. We leveraged RNA sequencing and genome-wide genotyping in islets from 188 donors to create the Islet Gene View (IGW) platform to make this information easily accessible to the scientific community. Expression data were related to islet phenotypes, diabetes status, other islet-expressed genes, islet hormone-encoding genes and for expression in insulin target tissues. The IGW web application produces output graphs for a particular gene of interest. In IGW, 284 differentially expressed genes (DEGs) were identified in T2D donor islets compared with controls. Forty percent of DEGs showed cell-type enrichment and a large proportion significantly co-expressed with islet hormone-encoding genes; glucagon (<i>GCG</i>, 56%), amylin (<i>IAPP</i>, 52%), insulin (<i>INS</i>, 44%), and somatostatin (<i>SST</i>, 24%). Inhibition of two DEGs, <i>UNC5D</i> and <i>SERPINE2</i>, impaired glucose-stimulated insulin secretion and impacted cell survival in a human ß-cell model. The exploratory use of IGW could help designing more comprehensive functional follow-up studies and serve to identify therapeutic targets in T2D.

Genome-wide association analyses highlight etiological differences underlying newly defined subtypes of diabetes.

Type 2 diabetes has been reproducibly clustered into five subtypes with different disease progression and risk of complications; however, etiological differences are unknown. We used genome-wide association and genetic risk score (GRS) analysis to compare the underlying genetic drivers. Individuals from the Swedish ANDIS (All New Diabetics In Scania) study were compared to individuals without diabetes; the Finnish DIREVA (Diabetes register in Vasa) and Botnia studies were used for replication. We show that subtypes differ with regard to family history of diabetes and association with GRS for diabetes-related traits. The severe insulin-resistant subtype was uniquely associated with GRS for fasting insulin but not with variants in the TCF7L2 locus or GRS reflecting insulin secretion. Further, an SNP (rs10824307) near LRMDA was uniquely associated with mild obesity-related diabetes. Therefore, we conclude that the subtypes have partially distinct genetic backgrounds indicating etiological differences.

Distinct Molecular Signatures of Clinical Clusters in People With Type 2 Diabetes: An IMI-RHAPSODY Study.

Type 2 diabetes is a multifactorial disease with multiple underlying aetiologies. To address this heterogeneity, investigators of a previous study clustered people with diabetes according to five diabetes subtypes. The aim of the current study is to investigate the etiology of these clusters by comparing their molecular signatures. In three independent cohorts, in total 15,940 individuals were clustered based on five clinical characteristics. In a subset, genetic (N = 12,828), metabolomic (N = 2,945), lipidomic (N = 2,593), and proteomic (N = 1,170) data were obtained in plasma. For each data type, each cluster was compared with the other four clusters as the reference. The insulin-resistant cluster showed the most distinct molecular signature, with higher branched-chain amino acid, diacylglycerol, and triacylglycerol levels and aberrant protein levels in plasma were enriched for proteins in the intracellular PI3K/Akt pathway. The obese cluster showed higher levels of cytokines. The mild diabetes cluster with high HDL showed the most beneficial molecular profile with effects opposite of those seen in the insulin-resistant cluster. This study shows that clustering people with type 2 diabetes can identify underlying molecular mechanisms related to pancreatic islets, liver, and adipose tissue metabolism. This provides novel biological insights into the diverse aetiological processes that would not be evident when type 2 diabetes is viewed as a homogeneous disease.

Replication and cross-validation of type 2 diabetes subtypes based on clinical variables: an IMI-RHAPSODY study.

AIMS/HYPOTHESIS: Five clusters based on clinical characteristics have been suggested as diabetes subtypes: one autoimmune and four subtypes of type 2 diabetes. In the current study we replicate and cross-validate these type 2 diabetes clusters in three large cohorts using variables readily measured in the clinic. METHODS: In three independent cohorts, in total 15,940 individuals were clustered based on age, BMI, HbA1c, random or fasting C-peptide, and HDL-cholesterol. Clusters were cross-validated against the original clusters based on HOMA measures. In addition, between cohorts, clusters were cross-validated by re-assigning people based on each cohort's cluster centres. Finally, we compared the time to insulin requirement for each cluster. RESULTS: Five distinct type 2 diabetes clusters were identified and mapped back to the original four All New Diabetics in Scania (ANDIS) clusters. Using C-peptide and HDL-cholesterol instead of HOMA2-B and HOMA2-IR, three of the clusters mapped with high sensitivity (80.6-90.7%) to the previously identified severe insulin-deficient diabetes (SIDD), severe insulin-resistant diabetes (SIRD) and mild obesity-related diabetes (MOD) clusters. The previously described ANDIS mild age-related diabetes (MARD) cluster could be mapped to the two milder groups in our study: one characterised by high HDL-cholesterol (mild diabetes with high HDL-cholesterol [MDH] cluster), and the other not having any extreme characteristic (mild diabetes [MD]). When these two milder groups were combined, they mapped well to the previously labelled MARD cluster (sensitivity 79.1%). In the cross-validation between cohorts, particularly the SIDD and MDH clusters cross-validated well, with sensitivities ranging from 73.3% to 97.1%. SIRD and MD showed a lower sensitivity, ranging from 36.1% to 92.3%, where individuals shifted from SIRD to MD and vice versa. People belonging to the SIDD cluster showed the fastest progression towards insulin requirement, while the MDH cluster showed the slowest progression. CONCLUSIONS/INTERPRETATION: Clusters based on C-peptide instead of HOMA2 measures resemble those based on HOMA2 measures, especially for SIDD, SIRD and MOD. By adding HDL-cholesterol, the MARD cluster based upon HOMA2 measures resulted in the current clustering into two clusters, with one cluster having high HDL levels. Cross-validation between cohorts showed generally a good resemblance between cohorts. Together, our results show that the clustering based on clinical variables readily measured in the clinic (age, HbA1c, HDL-cholesterol, BMI and C-peptide) results in informative clusters that are representative of the original ANDIS clusters and stable across cohorts. Adding HDL-cholesterol to the clustering resulted in the identification of a cluster with very slow glycaemic deterioration.

Genetic factors affect the susceptibility to bacterial infections in diabetes.

Diabetes increases the risk of bacterial infections. We investigated whether common genetic variants associate with infection susceptibility in Finnish diabetic individuals. We performed genome-wide association studies and pathway analysis for bacterial infection frequency in Finnish adult diabetic individuals (FinnDiane Study; N = 5092, Diabetes Registry Vaasa; N = 4247) using national register data on antibiotic prescription purchases. Replication analyses were performed in a Swedish diabetic population (ANDIS; N = 9602) and in a Finnish non-diabetic population (FinnGen; N = 159,166). Genome-wide data indicated moderate but significant narrow-sense heritability for infection susceptibility (h2 = 16%, P = 0.02). Variants on chromosome 2 were associated with reduced infection susceptibility (rs62192851, P = 2.23 × 10-7). Homozygotic carriers of the rs62192851 effect allele (N = 44) had a 37% lower median annual antibiotic purchase rate, compared to homozygotic carriers of the reference allele (N = 4231): 0.38 [IQR 0.22-0.90] and 0.60 [0.30-1.20] respectively, P = 0.01). Variants rs6727834 and rs10188087, in linkage disequilibrium with rs62192851, replicated in the FinnGen-cohort (P < 0.05), but no variants replicated in the ANDIS-cohort. Pathway analysis suggested the IRAK1 mediated NF-κB activation through IKK complex recruitment-pathway to be a mediator of the phenotype. Common genetic variants on chromosome 2 may associate with reduced risk of bacterial infections in Finnish individuals with diabetes.

Genetic analysis of obstructive sleep apnoea discovers a strong association with cardiometabolic health.

There is currently limited understanding of the genetic aetiology of obstructive sleep apnoea (OSA). We aimed to identify genetic loci associated with OSA risk, and to test if OSA and its comorbidities share a common genetic background.We conducted the first large-scale genome-wide association study of OSA using the FinnGen study (217 955 individuals) with 16 761 OSA patients identified using nationwide health registries.We estimated 0.08 (95% CI 0.06-0.11) heritability and identified five loci associated with OSA (p<5.0×10-8): rs4837016 near GAPVD1 (GTPase activating protein and VPS9 domains 1), rs10928560 near CXCR4 (C-X-C motif chemokine receptor type 4), rs185932673 near CAMK1D (calcium/calmodulin-dependent protein kinase ID) and rs9937053 near FTO (fat mass and obesity-associated protein; a variant previously associated with body mass index (BMI)). In a BMI-adjusted analysis, an association was observed for rs10507084 near RMST/NEDD1 (rhabdomyosarcoma 2 associated transcript/NEDD1 γ-tubulin ring complex targeting factor). We found high genetic correlations between OSA and BMI (rg=0.72 (95% CI 0.62-0.83)), and with comorbidities including hypertension, type 2 diabetes, coronary heart disease, stroke, depression, hypothyroidism, asthma and inflammatory rheumatic disease (rg>0.30). The polygenic risk score for BMI showed 1.98-fold increased OSA risk between the highest and the lowest quintile, and Mendelian randomisation supported a causal relationship between BMI and OSA.Our findings support the causal link between obesity and OSA, and the joint genetic basis between OSA and comorbidities.

Adult-onset diabetes in Middle Eastern immigrants to Sweden: Novel subgroups and diabetic complications-The All New Diabetes in Scania cohort diabetic complications and ethnicity.

BACKGROUND: Middle Eastern immigrants to Europe represent a high risk population for type 2 diabetes. We compared prevalence of novel subgroups and assessed risk of diabetic macro- and microvascular complications between diabetes patients of Middle Eastern and European origin. METHODS: This study included newly diagnosed diabetes patients born in Sweden (N = 10641) or Iraq (N = 286), previously included in the All New Diabetes in Scania cohort. The study was conducted between January 2008 and August 2016. Patients were followed to April 2017. Incidence rates in diabetic macro- and microvascular complications were assessed using cox-regression adjusting for the confounding effect of age at onset, sex, anthropometrics, glomerular filtration rate (eGFR) and HbA1c. FINDINGS: In Iraqi immigrants versus native Swedes, severe insulin-deficient diabetes was almost twice as common (27.9 vs. 16.2% p < 0.001) but severe insulin-resistant diabetes was less prevalent. Patients born in Iraq had higher risk of coronary events (hazard ratio [HR] 1.84, 95% CI 1.06-3.12) but considerably lower risk of chronic kidney disease (CKD) than Swedes (HR 0.19; 0.05-0.76). The lower risk in Iraqi immigrants was partially attributed to better eGFR. Genetic risk scores (GRS) showed more genetic variants associated with poor insulin secretion but lower risk of insulin resistance in the Iraqi than native Swedish group. INTERPRETATION: People with diabetes, born in the Middle East present with a more insulin-deficient phenotype and genotype than native Swedes. They have a higher risk of coronary events but lower risk of CKD. Ethnic differences should be considered in the preventive work towards diabetes and its complications.