Skip to main content

Diabetes risk loci-associated pathways are shared across metabolic tissues

Abstract

Aims/hypothesis

Numerous genome-wide association studies have been performed to understand the influence of genetic variation on type 2 diabetes etiology. Many identified risk variants are located in non-coding and intergenic regions, which complicates understanding of how genes and their downstream pathways are influenced. An integrative data approach will help to understand the mechanism and consequences of identified risk variants.

Methods

In the current study we use our previously developed method CONQUER to overlap 403 type 2 diabetes risk variants with regulatory, expression and protein data to identify tissue-shared disease-relevant mechanisms.

Results

One SNP rs474513 was found to be an expression-, protein- and metabolite QTL. Rs474513 influenced LPA mRNA and protein levels in the pancreas and plasma, respectively. On the pathway level, in investigated tissues most SNPs linked to metabolism. However, in eleven of the twelve tissues investigated nine SNPs were linked to differential expression of the ribosome pathway. Furthermore, seven SNPs were linked to altered expression of genes linked to the immune system. Among them, rs601945 was found to influence multiple HLA genes, including HLA-DQA2, in all twelve tissues investigated.

Conclusion

Our results show that in addition to the classical metabolism pathways, other pathways may be important to type 2 diabetes that show a potential overlap with type 1 diabetes.

Peer Review reports

Introduction

Several large genome-wide association studies (GWASs) have been performed to understand the genetic drivers of type 2 diabetes (T2D). The most recent GWAS in almost one million people identified 403 variants [13]. While for some risk variants the underlying mechanisms are relatively well understood, for most the mechanisms are largely unclear. In an attempt to elucidate such mechanisms, previous studies have undertaken efforts to integrate data from public repositories or perform functional follow-up of loci [1, 4]. For example, a recent study used public data to assign scores to genes near 101 T2D risk variants to identify the causal genes [4]. Most studies have used the precalculated cis expression quantitative trait loci (eQTLs) from the GTEx project [4, 5]. GTEx provides eQTLs in a one Mb region around transcription start sites of genes, which could lead to missed eQTLs of variants in intergenic regions or with more distant genes. In addition, most studies, focused on a single tissue rather than considering all diabetes-relevant tissues. Finally, only few studies used other layers of information available in public data repositories such as DNA methylation-, miRNA -, protein-QTLs, chromatin interactions, chromatin state segmentations and transcription factor binding sites and expression data.

In the current study, we investigate 403 previously identified T2D-associated single nucleotide polymorphisms (SNPs) to gain insight into the pathways under influence of T2D risk variants. We overlap SNPs in multiple tissues and molecular levels, including chromatin state segmentations, multiple QTL modalities, including cis and trans eQTLs. Instead of focusing only on the eQTLs we investigate co-expression networks of eQTLs and look at the effects of eQTL shared by tissues. Our results confirm pathways known in T2D including metabolic pathways, but also identify potential other pathways including the ribosome and auto-immunity pathways.

Results

Four hundred three previously published T2D-associated SNPs [1] were investigated both individually and together (Fig. 1).

Fig. 1
figure 1

Setup of the current study

The SLC22A3 locus influences SLC22A3 and LPA gene expression and plasma LPA protein levels

SNPs were compared to expression QTLs (eQTLs), protein QTLs (cis and trans pQTLs) and metabolite QTLs (mQTLs). Among the 403 SNPs, there were 189 eQTLs across the 12 metabolic tissues investigated (Table S2), 1 cis-pQTL, 2 trans-pQTLs and 6 mQTLs. One SNP rs1260326 was an eQTL, pQTL and mQTL (Fig. 2a, Table S3). There were multiple eQTLs in multiple tissues. Rs1260326-C was associated with lower SNX17 mRNA levels in the muscle (β = -0.32, P = 1.29·10–32, Table S3) and lower mRNA levels of NRBP1 in subcutaneous and visceral fat, whole blood, intestine, transverse colon. Of note, based on the H4 posterior probability (H4 PP), type 2 diabetes and NRBP1 (H4 PP ≥ 0.9) shared causal variants in five tissues, while the H4 PP for SNX17 was much lower (0.483). Rs1260326 has been associated with Glucokinase Regulatory Protein, but only in thyroid, this SNP was an eQTL with GCKR (β = 0.22, P = 1.84·10–8). In plasma, rs1260326-T was associated with upregulated alanine and multiple lipid levels, including very-low-density lipoproteins (VLDL) and triacylglycerols levels. Finally, in trans rs1260326-T was associated with higher levels of plasma Insulin-like growth factor-binding protein 1 (IGFBP1, P = 2.42·10–13, H4 PP = 0.873), Kallikrein B1 (KLKB1, P = 2.13·10–10), although the latter showed a very low H4 PP (4.23·10–46). KLKB1 is involved in blood coagulation and IGFBP1 is involved in metabolism.

Fig. 2
figure 2

Tissue shared eQTLs for diabetes-associated SNPs. a Overlap between expression, protein (cis), protein (trans) and metabolite QTLs. b QTL effects of rs474513 on SLC22A3 and LPA c Effect of rs474513 on LPA mRNA expression in the liver. X-axis, genotype; y-axis, normalized expression levels. d Effect of rs474513 on SLC22A3 mRNA expression in the liver. X-axis, genotype; y-axis, normalized expression levels. e Volcano plot of all the 1619 calculated eQTL—eGene pairs, with the tissue-shared eQTL highlighted. X-axis, normalized effect size; y-axis, P-value. Dots with similar colors represent the same eQTLs in multiple tissues

One SNP rs474513 was an eQTL and a cis pQTL (Fig. 2a, Table S4). In the liver, rs474513-A was an eQTL with apolipoprotein(a) (LPA, β = 0.29, P = 1.29·10–5, Figs. 2a and c). The A-allele of rs474513 gives an increased type 2 diabetes risk [1]. In plasma, rs474513-A influenced LPA protein levels (β = 0.23, P = 8.27·10–37). In multiple tissues, rs474513-A was an eQTL for SLC22A3 with an increase in expression, interestingly the liver was the sole tissue with decreased expression (β = -0.53, P = 2.23·10–12, Fig. 2d). Interestingly SLC22A3 encodes the organic cation transporter 3 (OCT3) which is involved in metformin transport [6]. Of note, the H4 PP do however suggest that these QTLs do not share causal variants with type 2 diabetes, limiting the implications and interpretation of these results.

The 189 SNPs that were identified as eQTL influenced the expression of 556 genes with in total 1619 eQTL-eGene pairs across investigated tissues (Fig. 2e, Table S2). Several diabetes risk variants increased the expression of genes in multiple tissues, including AP3S2 (rs4932265-T), CCDC92 (rs7978610-G), HLA-DQA2 (rs601945-G) and a lncRNA RP11-252K23.2 (Table S2, Fig. S1, rs3115960-G). Genes that showed decreased expression with diabetes-risk alleles included CEP68 (rs2052261-A, rs2249105-A), HSD17B12 (rs1061810-A) and the long non-coding RNA (lncRNA) RP11-613D13.10 (rs1061810-A). HLA-DQA2 is involved in multiple disease- and immune response-related pathways [7]. The two SNPs that influenced the expression of Centrosomal Protein 68 (CEP68) were in LD (r2 = 0.798). CEP68 is involved in centrosome cohesion [8]. AP3S2 is part of the AP-3 complex which is associated with the Golgi region and involved in vesicle transport [9]. HSD17B12 encodes Hydroxysteroid 17-Beta Dehydrogenase 12, which is involved in synthesis of fatty acids [7]. Rs1061810 influenced both the expression of HSD17B12 and RP11-613D13.10, where the latter is likely the HSD17B12 antisense, given that the expression of both transcripts was strongly correlated (Spearman’s ρ = 0.65–0.90, Fig. S2). CCDC92 has been associated with impaired adipogenesis [10]. Of note, based on the posterior probability, only HSD17B12, R11-613D13.10 shared causal variants with type 2 diabetes (Table S2).

Type 2 diabetes-associated eQTLs link to ribosome and autoimmunity pathways

Next, eQTLs were investigated in more detail to find common and tissue-specific pathways across tissues. The strongest enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were not those linked to metabolism (Table S5). Instead, the strongest enriched pathway was the Ribosome, which was enriched in six of the twelve tissues (Table S5, Fig. S3). Also when the enrichment was based on REACTOME instead of KEGG, the top enriched pathway was Eukaryotic Translation Elongation, in which the ribosomes play a key role (Table S6). The ribosome pathway was followed by pathways that are related to immunity and all these pathways were driven by the same set of genes as indicated by the lines connecting pathways (Fig. 3a-b, Fig. S4). Immune-related pathways were identified in 10 of the 12 tissues. For the ribosome pathway the number of eQTLs varied across tissues from a single in pituitary up to nine SNPs in subcutaneous fat. Across the twelve tissues investigated, two eQTL-eGene pairs were most consistently observed, that is rs12719778-T/RPL8 and rs12920022-A/RPL13 (Fig. 3c). The biggest effect size of rs12719778-T/RPL8 was observed in whole blood (NES = -0.10, P = 8.46·10–13, Fig. 3d), while for rs12920022-A/RPL13 the largest normalized effect size was observed in skeletal muscle (NES = -0.26, P = 3.36·10–25, Fig. 3e). The pathways that were identified in less tissues were those related to metabolism, for example in seven tissues Biosynthesis of unsaturated fatty acids which was driven by the effect of rs1061810 on HSD17B12 (Table S5).

Fig. 3
figure 3

Diabetes eQTLs are linked to metabolic pathways, including the ribosome pathway a-b. Enriched pathways in the pancreas (a) and the muscle (b). Grey lines between nodes indicate overlap between eGenes. c eQTLs associated with the Ribosome pathway. X-axis tissues investigated, y-axis eGenes and eQTLs. Circles indicate significant P-values. d Effect of rs12719778 on RPL8 mRNA expression in whole blood. X-axis, genotype; y-axis, normalized expression levels. e Effect of rs12920022 on RPL13 mRNA expression in skeletal muscle. X-axis, genotype; y-axis, normalized expression levels

Multiple diabetes-associated SNPs link to genes linked to autoimmunity

The immune-related pathways were generally enriched based on the same gene set in all tissues investigated (Fig. 3a-b, Fig. S4). Across tissues, the Antigen processing and presentation pathway was the strongest enriched pathway of the immune-related pathways (Table S5) which involved in total seven eQTLs and 18 genes. Comparison of the KEGG enrichment with that based on REACTOME resulted also in a strong enrichment for antigen presentation pathways, including PD-1 signaling and MHC class II antigen presentation (Table S6). Rs601945 was the eQTL that influenced the most genes involved (Fig. 4a). Rs601945 is a QTL in the HLA region, with the strongest positive effect on HLA-DQA2 in the skeletal muscle (NES = 1.19, P = 2.50·10–78, Fig. 4b), while rs601945 had a negative effect on HLA-DQB1 with the strongest effect, again, in skeletal muscle (NES = -0.50, P = 1.18·10–16, Fig. 4c). Based on histone modifications across tissues, the HLA locus was generally quiescent, except for blood cells where in multiple blood cell types enhancers were found (Fig. S5). In line with this, multiple chromatin interactions were observed in blood cells (CD34+, CD4+ memory, CD4+ naïve and CD4+ T-cells, Fig. 4d), 40 of which are interactions with loci located in HLA genes, with the top genes HLA-DQA1 (19 interactions) and HLA-DQB1 (13 interactions, Fig. 4d). Rs601945 has previously also been identified as a risk factor for autoimmunity diseases, including type 1 diabetes (P = 5.72·10–80), ulcerative colitis (P = 1.05·10–28) and inflammatory bowel disease (P = 2.34·10–24) based on data from the T1D portal (https://t1d.hugeamp.org).

Fig. 4
figure 4

Diabetes eQTLs linked to immune system-related pathways. a eQTLs associated with the Antigen processing and presentation pathway driven by rs601945. X-axis tissues investigated, y-axis eGenes and eQTLs. b Effect of rs601945 on HLA-DQA2 mRNA expression in the skeletal muscle. X-axis, genotype; y-axis, normalized expression levels. c Effect of rs12719778 on HLA-DQB1 mRNA expression in the skeletal muscle. X-axis, genotype; y-axis, normalized expression levels. d Circos plot of the rs601945 locus. Chromatin interactions in blood cells (CD34+, CD4+ memory, CD4+ naïve and CD4+ T) of the region that is in LD (r2 ≥ 0.8, lead SNP is shown in red) with rs601945 are shown as colored arcs. Chromatin states are shown for CD4+ memory (inner-ring) and CD4+ naïve cells respectively (middle-ring). Genes are shown in the outer ring in gray and red

Moreover, rs601945 was associated with blood cell counts of neutrophils (P = 2.53·10–112), eosinophils (P = 1.36·10–37) and monocytes (P = 1.70·10–15).

Discussion

In this study, we investigated 403 T2D-associated SNPs in more detail. We show that variants can have effects across multiple molecular layers, including gene, protein and metabolite levels. Moreover, risk variants are associated with altered expression of the same genes in multiple tissues, including HSD17B12, AP3S2, HLA-DQA2. This was also reflected in observed enriched pathways, where the same pathways were influenced by genetic risk for diabetes across tissues, including antigen processing, ribosome, proteasome and protein processing in endoplasmic reticulum.

One T2D-associated variant, intronic of SLC22A3 was both an expression- and protein-QTL. The T2D risk allele A of rs474513 was associated with higher expression of LPA expression in the liver and higher protein levels in plasma. High LPA levels have been associated with higher risk on cardiovascular disease events in people with diabetes [11]. In addition to LPA, rs474513-A was associated with higher expression of SLC22A3 (OCT3).

Several T2D-associated risk-alleles showed tissue-shared eQTL effects, including AP3S2, CCDC92, HLA-DQA2, CEP68, HSD17B12 and two lncRNAs RP11-613D13.10 and RP11-252K23.2. Both AP3S2 and HSD17B12 have previously been found in relation to T2D, but in limited number of tissues. AP3S2 in human pancreatic islets [12] and HSD17B12 in adipose, liver, muscle tissue and whole blood [3], which are relevant for T2D [13]. One of the strongest observed effects was between rs4932265 and AP3S2. AP3S2 is a subunit of the AP-3 complex which is involved in budding of vesicles from the Golgi membrane [9]. AP3S2 has been linked to T2D in six different GWASs [1419] investigating various populations (South Asian, Japanese and European ancestry) and with four different SNPs, three of which (rs12912009, rs2028299, rs8031576) are in LD (r2 ≥ 0.80) with rs4932265. In the current study we observed that AP3S2 has a higher expression in the twelve tissues in individuals carrying the risk allele of rs4932265. Despite increasing evidence for the role AP3S2 in T2D susceptibility it remains unclear how AP3S2 is involved, although there is some evidence pointing at a beta-cell defect [20]. Nonetheless, our data suggest that the effect may be more tissue-shared than limited to one cell type.

For HSD17B12, we show that the diabetes-risk risk allele of rs1061810 was associated with a lower expression in all twelve tissues investigated. The HSD17B12 gene encodes a bifunctional enzyme involved in the biosynthesis of estradiol and the elongation of very long chain fatty acids [21]. This result corroborates the finding that HSD17B12 expression is downregulated in the adipose tissue of insulin-resistant subjects [22] and that HSD17B12 plays a role in adipogenesis [23]. Rs1061810 was also associated with altered expression of RP11-613D13.10 and which is the antisense of HSD17B12. The role of the lncRNAs RP11-613D13.10 and RP11-252K23.2 is not clear. Furthermore, we observe in the current study that based on the H4 PP causal variants were shared between type 2 diabetes and HSD17B12 / R11-613D13.10.

On the pathway level as expected most SNPs linked to metabolic pathways. The metabolic pathways as curated by KEGG [7] consists of 1489 genes and is an encompassing term for all pathways that are involved in metabolism. Our results show that SNPs that are directly linked to metabolism do not influence a single metabolic process but are scattered among various metabolic pathways. Due to this dispersion of SNPs between numerous pathways it remains difficult to assign groups of SNPs to specific processes in specific tissues. This together with the variety of pathways to which SNPs are mapped shows that T2D has a lot of different points of engagement through which it can originate and progress, which is accordance with the heterogeneous nature of T2D [24].

In eleven of the twelve tissues, eQTLs were enriched for the ribosome pathway. Genetic susceptibility to T2D has previously not been linked to a decreased expression of ribosomal genes, although the association between ribosomal content and T2D has extensively been studied [2527]. Insulin and ribosomal content are tightly connected, where insulin stimulates the synthesis of ribosomal proteins in various tissues [28, 29] and a loss of ribosomal proteins is associated with an inhibition of AKT phosphorylation activity/insulin pathway [30].

Multiple enriched pathways linked to immunity and the eGenes in these pathways were mainly from the HLA class. Rs601945 was the key SNP in these pathways as it influences the expression of multiple HLA genes across tissues. HLA genes have previously been associated with T2D [3, 17], however, our results reveal that the effects are widespread as its association with altered expression of various HLA genes was observed in all investigated tissues. Interestingly, while the HLA region represents the highest risk for T1D [31], our results are pointing to a connection between HLA-DQA2 and T2D. Rs601945 has also been associated with other diseases where autoimmunity plays a key role, including type 1 diabetes, ulcerative colitis, inflammatory bowel disease. Moreover, rs601945 is strongly associated with blood cell counts of neutrophils, eosinophils and monocytes. These results suggest that the r601945 is associated with autoimmunity. These results also support previous results that diabetes could also be considered as a continuum rather than two separate diseases with overlap [32, 33]. Our data support that T2D has an immuno-metabolic component involving, like T1D, members of both innate and adaptive immune response. In addition, our pQTL analyses also highlighted immune response pathways. Also, in a previous study in blood of persons with T2D we found that HbA1c is associated with altered expression of immune response-related genes [34].

Strength of our study is the hypothesis-free approach in multiple tissues, which allowed us to investigate tissue-shared effects of T2D-associated SNPs. A limitation of our study is that we mostly rely on eQTLs where we do not know whether the observed changes in expression also translate to changes in protein levels. A second limitation is that we use QTL data from different sets of individuals, while ideally one would use regulatory, expression and protein data from the same set of individuals. Third, we rely on data from European descent, which limits the generalizability to other ethnicities. Fourth, a limitation is that in the current study we do not validate our findings in vitro, which is required to fully understand the observed results. Finally, the pQTLs included in the current study were only measured in plasma and are therefore not necessarily representative of the pQTLs within tissues.

Conclusion

Altogether, our data show biological processes that are subject to genetic influences. We show that they are not necessarily limited to single tissues but are shared across diabetes-relevant tissues. Our findings highlight the importance of an integrative tissue-wide approach where risk loci for T2D are not only seen as individual risk factors but also as a network of risk factors that may play a role across tissues.

Methods

Four hundred three SNPs from the GWAS summary statistics reported by Mahajan et al. [1] were extracted. SNPs plus those in LD were used in subsequent analyses, where LD was defined as r2 ≥ 0.80 around the respective variants. The SNPs were analyzed using our previously published package CONQUER [35]. Both the single SNP mode and multi-SNP mode were used to analyze the SNPs of interest. For the single SNP analysis we extracted data from relevant databases using our R-package CONQUER, including expression quantitative trait loci (eQTLs) in cis and trans, DNA methylation-, protein-QTLs, chromatin interactions, chromatin state segmentations, transcription factor binding sites and expression data [35].

QTL data

eQTLs were obtained from GTEx based on data from mainly European descent. All lead SNPs were tested against all genes (eGenes) in cis and trans on the GTEx API. eGenes are genes under influence of an eQTL. Trans-eQTLs were those more distant from the SNP but on the same chromosome and the tested region was defined by the range of the predicted chromosomal interactions. SNPs were considered significant if the P-value was below the GTEx threshold. Effect sizes in figures and tables are those for the effect allele of T2D relative to the alternative allele based on an additive model. The latter threshold is defined as the empirical P-value of the gene closest to the 0.05 FDR threshold. Significant protein-QTLs (pQTLs, P < 5·10–8) obtained from Yao et al. [36] and filtered for the 403 lead SNPs plus SNPs in LD with those SNPs. QTLs in Yao et al. [36] have been identified in people from European descent. Significant metabolite QTLs (mQTLs, P < 5·10–8) were obtained Gallois et al. [37] based on a Finnish population.

GTEx gene expression data

Expression data was obtained from GTEx V8 for tissues relevant for the etiology of T2D, including subcutaneous and visceral fat, sigmoid- and transverse colon, liver, skeletal muscle, pancreas, pituitary, terminal ileum of the small intestine, stomach, thyroid and whole blood. Of the included tissues, sample sizes range from N = 187 (terminal ileum) to N = 803 (skeletal muscle). The percentage males was relatively higher (63.1%—72.1%, Table S1) with the majority middle-aged (50–69 years, Table S1).

Colocalization analysis

To compare whether two traits shared a common causal variant, a colocalization analysis was performed using the R-package coloc. The measure for a single shared causal variant is the H4 posterior probability, where a H4 PP larger than 0.9 indicates a single causal variant shared by two traits. In the current study, we use the colocalization to compare the type 2 diabetes SNPs to the protein and expression QTLs.

Multi-SNP analysis

Modularization and pathway enrichment for included tissues was performed on all significant eQTLs (see above). Genes co-expressed with eGenes were identified (ρ ≥|0.90|, Spearman's rank correlation coefficient). The latter was done to identify entire pathways that are under influence of one or more SNPs. The eGenes and co-expressed genes were hierarchical clustered [38, 39] based on pairwise distance between genes (1 – ρ). The number of modules within the clustered data was optimized by maximizing the globalSEmax of the gap statistic [40] using the cluster R package [41]. Modules of co-expressed genes and eGenes were tested for pathway enrichment based on KEGG [7] pathways. For each module, association with pathways was determined with Fisher’s exact test [42], which resulted in odds ratios and accompanying P-values on the association. If a module did not contain an eQTL or was not enriched for a pathway, it was omitted from the results. Figures were produced in and obtained from CONQUER or additionally made using the R-package ggplot2 version 3.2.1. [43].

Availability of data and materials

All data presented in the current manuscript are publicly available and implemented in the R-package CONQUER.db (https://github.com/roderickslieker/CONQUER.db). eQTLs and expression data were obtained from GTEx. Pathways were obtained from KEGG (https://www.genome.jp/kegg/).

References

  1. Mahajan A, Taliun D, Thurner M, Robertson NR, Torres JM, Rayner NW, Steinthorsdottir V, Scott RA, Grarup N, Cook JP, et al. Fine-mapping of an expanded set of type 2 diabetes loci to single-variant resolution using high-density imputation and islet-specific epigenome maps. Nat Genet. 2018;50(11):1505–13.

    Article  CAS  Google Scholar 

  2. Xue A, Wu Y, Zhu Z, Zhang F, Kemper KE, Zheng Z, Yengo L, Lloyd-Jones LR, Sidorenko J, Wu Y, et al. Genome-wide association analyses identify 143 risk variants and putative regulatory mechanisms for type 2 diabetes. Nat Commun. 2018;9(1):2941.

    Article  Google Scholar 

  3. Scott RA, Scott LJ, Mägi R, Marullo L, Gaulton KJ, Kaakinen M, Pervjakova N, Pers TH, Johnson AD, Eicher JD, et al. An expanded genome-wide association study of type 2 diabetes in Europeans. Diabetes. 2017;66(11):2888–902.

    Article  CAS  Google Scholar 

  4. Fernández-Tajes J, Gaulton KJ, van de Bunt M, Torres J, Thurner M, Mahajan A, Gloyn AL, Lage K, McCarthy MI. Developing a network view of type 2 diabetes risk pathways through integration of genetic, genomic and functional data. Genome Med. 2019;11(1):19.

    Article  Google Scholar 

  5. Cirillo E, Kutmon M, Gonzalez Hernandez M, Hooimeijer T, Adriaens ME, Eijssen LMT, Parnell LD, Coort SL, Evelo CT. From SNPs to pathways: Biological interpretation of type 2 diabetes (T2DM) genome wide association study (GWAS) results. PLoS ONE. 2018;13(4):e0193515.

    Article  Google Scholar 

  6. Nies AT, Hofmann U, Resch C, Schaeffeler E, Rius M, Schwab M. Proton pump inhibitors inhibit metformin uptake by organic cation transporters (OCTs). PLoS One. 2011;6(7):e22163.

    Article  CAS  Google Scholar 

  7. Kanehisa M, Goto S. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 2000;28(1):27–30.

    Article  CAS  Google Scholar 

  8. Graser S, Stierhof YD, Nigg EA. Cep68 and Cep215 (Cdk5rap2) are required for centrosome cohesion. J Cell Sci. 2007;120(Pt 24):4321–31.

    Article  CAS  Google Scholar 

  9. Odorizzi G, Cowles CR, Emr SD. The AP-3 complex: a coat of many colours. Trends Cell Biol. 1998;8(7):282–8.

    Article  CAS  Google Scholar 

  10. Lotta LA, Gulati P, Day FR, Payne F, Ongen H, van de Bunt M, Gaulton KJ, Eicher JD, Sharp SJ, Luan J, et al. Integrative genomic analysis implicates limited peripheral adipose storage capacity in the pathogenesis of human insulin resistance. Nat Genet. 2017;49(1):17–26.

    Article  CAS  Google Scholar 

  11. Saeed A, Sun W, Agarwala A, Virani SS, Nambi V, Coresh J, Selvin E, Boerwinkle E, Jones PH, Ballantyne CM, et al. Lipoprotein(a) levels and risk of cardiovascular disease events in individuals with diabetes mellitus or prediabetes: the Atherosclerosis Risk in Communities study. Atherosclerosis. 2019;282:52–6.

    Article  CAS  Google Scholar 

  12. Khamis A, Canouil M, Siddiq A, Crouch H, Falchi M, Bulow Mv, Ehehalt F, Marselli L, Distler M, Richter D, et al. Laser capture microdissection of human pancreatic islets reveals novel eQTLs associated with type 2 diabetes. Mol Metab. 2019;24:98–107.

    Article  CAS  Google Scholar 

  13. DeFronzo RA. From the triumvirate to the ominous octet: a new paradigm for the treatment of type 2 diabetes mellitus. Diabetes. 2009;58(4):773.

    Article  CAS  Google Scholar 

  14. Kooner JS, Saleheen D, Sim X, Sehmi J, Zhang W, Frossard P, Been LF, Chia K-S, Dimas AS, Hassanali N, et al. Genome-wide association study in individuals of South Asian ancestry identifies six new type 2 diabetes susceptibility loci. Nat Genet. 2011;43(10):984–9.

    Article  CAS  Google Scholar 

  15. Hara K, Fujita H, Johnson TA, Yamauchi T, Yasuda K, Horikoshi M, Peng C, Hu C, Ma RCW, Imamura M, et al. Genome-wide association study identifies three novel loci for type 2 diabetes. Hum Mol Genet. 2013;23(1):239–46.

    Article  Google Scholar 

  16. Bonàs-Guarch S, Guindo-Martínez M, Miguel-Escalada I, Grarup N, Sebastian D, Rodriguez-Fos E, Sánchez F, Planas-Fèlix M, Cortes-Sánchez P, González S, et al. Re-analysis of public genetic data reveals a rare X-chromosomal variant associated with type 2 diabetes. Nat Commun. 2018;9(1):321.

    Article  Google Scholar 

  17. Zhao W, Rasheed A, Tikkanen E, Lee J-J, Butterworth AS, Howson JMM, Assimes TL, Chowdhury R, Orho-Melander M, Damrauer S, et al. Identification of new susceptibility loci for type 2 diabetes and shared etiological pathways with coronary heart disease. Nat Genet. 2017;49(10):1450–7.

    Article  CAS  Google Scholar 

  18. Suzuki K, Akiyama M, Ishigaki K, Kanai M, Hosoe J, Shojima N, Hozawa A, Kadota A, Kuriki K, Naito M, et al. Identification of 28 new susceptibility loci for type 2 diabetes in the Japanese population. Nat Genet. 2019;51(3):379–86.

    Article  CAS  Google Scholar 

  19. Kichaev G, Bhatia G, Loh P-R, Gazal S, Burch K, Freund MK, Schoech A, Pasaniuc B, Price AL. Leveraging polygenic functional enrichment to improve GWAS power. Am J Hum Genet. 2019;104(1):65–75.

    Article  CAS  Google Scholar 

  20. Wood AR, Jonsson A, Jackson AU, Wang N, van Leewen N, Palmer ND, Kobes S, Deelen J, Boquete-Vilarino L, Paananen J, et al. A genome-wide association study of IVGTT-based measures of first-phase insulin secretion refines the underlying physiology of type 2 diabetes variants. Diabetes. 2017;66(8):2296.

    Article  CAS  Google Scholar 

  21. Hiltunen JK, Kastaniotis AJ, Autio KJ, Jiang G, Chen Z, Glumoff T. 17B-hydroxysteroid dehydrogenases as acyl thioester metabolizing enzymes. Mol Cell Endocrinol. 2019;489:107–18.

    Article  CAS  Google Scholar 

  22. Elbein SC, Kern PA, Rasouli N, Yao-Borengasser A, Sharma NK, Das SK. Global gene expression profiles of subcutaneous adipose and muscle from glucose-tolerant, insulin-sensitive, and insulin-resistant individuals matched for BMI. Diabetes. 2011;60(3):1019–29.

    Article  CAS  Google Scholar 

  23. Söhle J, Machuy N, Smailbegovic E, Holtzmann U, Grönniger E, Wenck H, Stäb F, Winnefeld M. Identification of new genes involved in human adipogenesis and fat storage. PLoS One. 2012;7(2):e31193–e31193.

    Article  Google Scholar 

  24. DeFronzo RA, Ferrannini E, Groop L, Henry RR, Herman WH, Holst JJ, Hu FB, Kahn CR, Raz I, Shulman GI, et al. Type 2 diabetes mellitus. Nat Rev Dis Primers. 2015;1(1):15019.

    Article  Google Scholar 

  25. Stirewalt WS, Wool IG, Cavicchi P. The relation of RNA and protein synthesis to the sedimentation of muscle ribosomes: effect of diabetes and insulin. Proc Natl Acad Sci U S A. 1967;57(6):1885–92.

    Article  CAS  Google Scholar 

  26. Wool IG, Kurihara K. Determination of the number of active muscle ribosomes: effect of diabetes and insulin. Proc Natl Acad Sci U S A. 1967;58(6):2401–7.

    Article  CAS  Google Scholar 

  27. Ozoe A, Sone M, Fukushima T, Kataoka N, Chida K, Asano T, Hakuno F, Takahashi S-I. Insulin receptor substrate-1 associates with small nucleolar RNA which contributes to ribosome biogenesis. Front Endocrinol. 2014;5:24.

    Article  Google Scholar 

  28. Proud CG. Regulation of protein synthesis by insulin. Biochem Soc Trans. 2006;34(Pt 2):213–6.

    Article  CAS  Google Scholar 

  29. Proud CG, Denton RM. Molecular mechanisms for the control of translation by insulin. Biochem J. 1997;328((Pt 2)(Pt 2)):329–41.

    Article  CAS  Google Scholar 

  30. Heijnen HF, van Wijk R, Pereboom TC, Goos YJ, Seinen CW, van Oirschot BA, van Dooren R, Gastou M, Giles RH, van Solinge W, et al. Ribosomal protein mutations induce autophagy through S6 kinase inhibition of the insulin pathway. PLoS Genet. 2014;10(5):e1004371.

    Article  Google Scholar 

  31. Bradfield JP, Qu H-Q, Wang K, Zhang H, Sleiman PM, Kim CE, Mentch FD, Qiu H, Glessner JT, Thomas KA, et al. A genome-wide meta-analysis of six type 1 diabetes cohorts identifies multiple associated loci. PLoS Genet. 2011;7(9):e1002293.

    Article  CAS  Google Scholar 

  32. Candia J, Cheung F, Kotliarov Y, Fantoni G, Sellers B, Griesman T, Huang J, Stuccio S, Zingone A, Ryan BM, et al. Assessment of variability in the SOMAscan assay. Sci Rep. 2017;7(1):14248.

    Article  Google Scholar 

  33. Brooks-Worrell B, Palmer JP. Is diabetes mellitus a continuous spectrum? Clin Chem. 2011;57(2):158–61.

    Article  CAS  Google Scholar 

  34. Slieker RC, van der Heijden AAWA, van Leeuwen N, Mei H, Nijpels G, Beulens JWJ. ’t Hart LM: HbA1c is associated with altered expression in blood of cell cycle- and immune response-related genes. Diabetologia. 2018;61(1):138–46.

    Article  CAS  Google Scholar 

  35. Bouland GA, Beulens JWJ, Nap J, van der Slik AR, Zaldumbide A, 't Hart LM, Slieker RC. CONQUER: an interactive toolbox to understand functional consequences of GWAS hits. NAR Genom Bioinform. 2020;2(4):lqaa085.

  36. Yao C, Chen G, Song C, Keefe J, Mendelson M, Huan T, Sun BB, Laser A, Maranville JC, Wu H, et al. Genome-wide mapping of plasma protein QTLs identifies putatively causal genes and pathways for cardiovascular disease. Nat Commun. 2018;9(1):3268.

    Article  Google Scholar 

  37. Gallois A, Mefford J, Ko A, Vaysse A, Julienne H, Ala-Korpela M, Laakso M, Zaitlen N, Pajukanta P, Aschard H. A comprehensive study of metabolite genetics reveals strong pleiotropy and heterogeneity across time and context. Nat Commun. 2019;10(1):4788.

    Article  Google Scholar 

  38. Leonard Kaufman PJR. Finding groups in data: an introduction to cluster analysis. New York: John Wiley & Sons, Inc.; 1990.

    Book  Google Scholar 

  39. Belbin L, Faith DP, Milligan GW. A comparison of two approaches to beta-flexible clustering. Multivar Behav Res. 1992;27(3):417–33.

    Article  CAS  Google Scholar 

  40. Dudoit S, Fridlyand J. A prediction-based resampling method for estimating the number of clusters in a dataset. Genome Biol. 2002;3(7):research0036.0031.

    Article  Google Scholar 

  41. Maechler M: Cluster R-Package. Cran 2019.

  42. Fisher RA. The logic of inductive inference. J Roy Stat Soc. 1935;98(1):39–82.

    Article  Google Scholar 

  43. Wickham H: ggplot2: Elegant Graphics for Data Analysis: Springer-Verlag New York; 2009.

Download references

Acknowledgements

None.

Funding

This project has received funding from the Innovative Medicines Initiative 2 Joint Undertaking under grant agreement No 115881 (RHAPSODY). This Joint Undertaking receives support from the European Union’s Horizon 2020 research and innovation programme and EFPIA. This work is supported by the Swiss State Secretariat for Education‚ Research and Innovation (SERI) under contract number 16.0097–2. The opinions expressed and arguments employed herein do not necessarily reflect the official views of these funding bodies.

Author information

Authors and Affiliations

Authors

Contributions

GAB, RCS, LMTH designed the study. GAB and RCS analyzed the data and wrote the draft manuscript. GAB and RCS are the guarantors of this work and, as such, had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis. All authors critically read and revised the manuscript and approved the final version of the manuscript.

Corresponding author

Correspondence to Roderick C. Slieker.

Ethics declarations

Ethics approval and consent to participate

Ethical approval and consent to participate was obtained by studies from which the data was used in the current study. The study and obtained data were carried out in accordance with the relevant guidelines and regulations.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1. a

Frequency of tissues in which a certain eGene was identified. X-axis, genes under influence of a diabetes SNP, y-axis, frequency of tissues in which the QTL was identified. Figure S2. Scatterplot of expression of RP11-613D13.10 versus the expression of HSD17B12. X-axis, expression of RP11-613D13.10; y-axis expression of HSD17B12. Figure S3. a Frequency of the number of times an enriched pathway was found across the twelve tested tissues. X-axis, pathway investigated; y-axis, frequency of tissues. Blue bars indicate pathways associated with HLA-genes. Figure S4. Relation between KEGG pathways identified in each of the tissues. Tissues include subcutaneous fat (a), visceral fat (b), sigmoid colon (c), transverse colon (d), pituitary (e), small intestine (f), stomach (g), thyroid (h), whole blood (i). Figure S5. Chromatin state segmentations for rs601945 in various tissues and cell lines (data from Epigenomics Roadmap). Plus-symbol indicates a SNP in LD and blood cell types are colored in blue. X-axis, location on the genome; y-axis cell type.

Additional file 2: Table S1.

Characteristics of the individuals in the GTEX data.

Additional file 3: Table S2.

eQTLs associated with T2D SNPs.

Additional file 4: Table S3. 

Protein, metabolite and expression QTLs associated with rs1260326.

Additional file 5: Table S4. 

Protein and expression QTLs associated with rs474513.

Additional file 6: Table S5. 

Enriched KEGG pathways based on T2D SNPs.

Additional file 7: Table S6. 

Enriched REACTOME pathways based on T2D SNPs.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Bouland, G.A., Beulens, J.W.J., Nap, J. et al. Diabetes risk loci-associated pathways are shared across metabolic tissues. BMC Genomics 23, 368 (2022). https://doi.org/10.1186/s12864-022-08587-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-022-08587-5

Keywords