High-throughput Genetic Clustering of Type 2 Diabetes Loci Reveals Heterogeneous Mechanistic Pathways of Metabolic Disease
==========================================================================================================================

* Hyunkyung Kim
* Kenneth E. Westerman
* Kirk Smith
* Joshua Chiou
* Joanne B. Cole
* Timothy Majarian
* Marcin von Grotthuss
* Josep M. Mercader
* Soo Heon Kwak
* Jaegil Kim
* Jose C. Florez
* Kyle Gaulton
* Alisa K. Manning
* Miriam S. Udler

## Abstract

**Aims/hypothesis** Type 2 diabetes (T2D) is highly polygenic and influenced by multiple biological pathways. Rapid expansion in the number of T2D loci can be leveraged to identify such pathways, thus facilitating improved disease management.

**Methods** We developed a high-throughput pipeline to enable clustering of T2D loci based on variant-trait associations. Our pipeline extracted summary statistics from genome-wide association studies (GWAS) for T2D and related traits to generate a matrix of 324 variant x 64 trait associations and applied Bayesian Non-negative Factorization (bNMF) to identify genetic components of T2D. We generated cluster-specific polygenic scores and performed regression analysis in an independent cohort (N=25,419) to assess for clinical relevance.

**Results** We identified ten clusters, replicating the five from our prior analysis as well as novel clusters related to beta-cell dysfunction, pronounced insulin secretion, and levels of alkaline phosphatase, lipoprotein-A, and sex hormone-binding globulin. Four clusters related to mechanisms of insulin deficiency, five to insulin resistance, and one had an unclear mechanism. The clusters displayed tissue-specific epigenomic enrichment, notably with the two beta-cell clusters differentially enriched in functional and stressed pancreatic beta-cell states. Additionally, cluster-specific polygenic scores were differentially associated with patient clinical characteristics and outcomes. The pipeline was applied to coronary artery disease and chronic kidney disease, identifying multiple shared genetic pathways with T2D.

**Conclusions/interpretation** Our approach stratifies T2D loci into physiologically meaningful genetic clusters associated with distinct tissues and clinical outcomes. The pipeline allows for efficient updating as additional GWAS become available and can be readily applied to other conditions, facilitating clinical translation of GWAS findings. Software to perform this clustering pipeline is freely available.

## Introduction

Type 2 Diabetes (T2D) has variable contributions of insulin resistance and beta cell dysfunction, and influenced by multiple risk factors, including genetics [1]. Untangling the heterogeneity of T2D may be fundamental to improving patient management and facilitating precision medicine.

Hundreds of loci associated with T2D have been identified in large-scale genetic studies, including genome-wide association studies (GWAS), however, the translation of these established T2D genetic loci into improved understanding of disease pathophysiology has been challenging, due in large part to the variants being non-protein coding [2–5]. Recent studies have leveraged a growing number of available GWAS datasets to connect genetic loci to mechanistic pathways by clustering loci based on shared patterns of associations across multiple traits [2, 6–8]. In our previous work [6], clustering was performed on 94 T2D variants identified by manual curation of published T2D GWAS manuscripts. Soft clustering analysis with Bayesian Non-negative Matrix Factorization (bNMF) of the associations of these 94 T2D variants with 47 diabetes-related traits identified five distinct clusters, recognizable as biological pathways of T2D. A similar set of five clusters of T2D loci were independently identified by Mahajan *et al*., along with a sixth cluster of “undetermined” physiological impact[2]. Of these five shared clusters, two related to beta-cell dysfunction, and the other three clusters represented different mechanisms of insulin resistance: obesity-mediated, abnormal lipodystrophy-like fat distribution, and altered hepatic lipid metabolism [9].

With new T2D loci continuously being discovered and additional GWAS trait summary statistics becoming publicly available, we sought to expand our prior work which involved manual curation of a smaller set of T2D loci. We developed a high-throughput pipeline to enable extraction of hundreds of genetic variants and traits from multiple GWAS to be used for cluster analysis in order to identify new genetic pathways of disease.

## Material and methods

### Pipeline for input variant-trait association matrix for clustering

An overview of preprocessing steps for variants and traits used for generating the input matrix for variant-trait association clustering analysis is illustrated in **Figure S1** with additional details in the **Supplementary Methods**. To obtain a comprehensive set of independent genetic variants associated with T2D, we extracted 21,666 variants reaching genome-wide significance (*P* < 5×10−8) from multiple large-scale T2D studies [2, 3, 10–14] in the AMP-Common Metabolic Disease Knowledge Portal (CMDKP) [15] (**Table S1**) and performed stringent LD-pruning of variants at r2 < 0.1 (**Table S2**) as well as filtering and removal of multi-allelic, ambiguous (A/T or C/G), or poorly represented (<80%) of trait GWAS datasets.

For trait selection, we utilized summary statistics available for 75 GWAS of glycemic traits, anthropometric traits, vital signs, and additional laboratory measures in the AMP-CMDKP or UK Biobank [15, 16] (**Table S3, Supplementary Methods**). Our goal was to let the genetics guide which traits were included in the clustering analysis, and thus traits were used only if the minimum *P*-value across the final set of variants was lower than a Bonferroni *P*-value cutoff of 0.05/N_final_variants (N=324). We then removed highly correlated traits (with |r| ≥ 0.85) to reduce redundancy.

For the selected lists of variants and traits, we utilized the GWAS summary statistics to generate a matrix of standardized Z-scores, choosing the T2D risk-increasing allele for each variant and dividing the estimated regression coefficient beta by the standard error. To account for the differences in sample size across trait GWAS studies, we scaled the standardized Z-scores by dividing by the square root of the sample size for each variant in each trait GWAS and then multiplied by the mean of square root of median sample size across all SNPs in each GWAS.

This pipeline was also used for coronary artery disease (CAD) and chronic kidney disease (CKD) with six CAD GWAS [17–19] and 39 CKD-related GWAS [11, 16, 20–23] queried.

### bNMF clustering

The variant-trait association matrix Z (m by n, m: # of variants, n: # of traits) was constructed as above. We then generated a non-negative input matrix X (2m by n) by concatenating two separate modifications of the original Z matrix: one containing all positive standardized Z-scores (zero otherwise) and the other all negative standardized Z-scores multiplied by -1. The bNMF procedure factorizes X into two matrices, W (2m by K) and HT (n by K), as X ∼ WH with an optimal rank K, corresponding to the association matrix of variants and traits to the number of clusters, respectively (**Supplementary Methods**) [6]. The key features for each cluster are determined by the most strongly associated traits, a natural output of the bNMF approach. bNMF algorithm was performed in R Studio for 1,000 iterations with maximum number of cluster K set to 20. To define a set of strongest-weighted variants in each cluster, we employed a method to determine a cluster weight cutoff that maximized the signal to noise ratio of weights (**Figure S2**).

### Cluster associations with relevant phenotypes using GWAS summary statistics

We generated GWAS-partitioned polygenic scores (GWAS pPS) for each cluster utilizing inverse-variance weighted fixed effects meta-analysis of GWAS summary statistics including the set of strongest-weighted variants above the weight cutoff for each cluster using the dmetar package in R [24] (**Supplementary Methods**). For testing T2D cluster associations with cardiometabolic outcomes, the significance threshold was set to 0.05/(7×K), representing a Bonferroni correction for K clusters and 7 outcomes (**Table S4**).

### Functional annotation and enrichment analysis

At each locus, we calculated approximate Bayes Factors (aBF) for all variants 500 kb upstream and downstream with r2 ≥ 0.1, with the index variant (100% credible set) from effect size estimates and standard errors, using the approach of Wakefield [25]. We then calculated a posterior probability for each variant by dividing the aBF by the sum of all aBF in the credible set. We obtained previously published 13-state ChromHMM [26] chromatin state calls for 28 cell types, excluding cancer cell lines [27]. We also compiled candidate cis-regulatory elements (cCREs) for 14 cell types and subtypes from published single cell chromatin accessibility datasets [28, 29]. We assessed enrichment of annotations within clusters by overlapping 100% credible set variants for signals in each cluster with cell type epigenomic annotations (chromatin states and cCREs). We also assessed epigenomic enrichment in single cell pancreatic tissue using a second method. As previously described [30], we subset loci from the Beta-cell 1 and 2 clusters, annotated variants using cCREs from INShigh and INSlow beta cells, and applied fgwas [31] in the fine mapping mode (**Supplementary Methods**).

### Partitioned Polygenic Score (pPS) analysis in the Mass General Brigham Biobank

The Mass General Brigham (MGB) Biobank (formerly Partners Biobank) includes more than 120,000 consented patients across the MGB healthcare system [32, 33] (**Supplementary Methods**). Approval for data analysis was obtained by the MGB IRB, study 2016P001018. We employed algorithmically defined outcome phenotypes for T2D, CAD, and CKD developed by the MGB Biobank [34]. Additional phenotypic data (laboratory measures, vital signs, and anthropometric measures) were extracted from outpatient clinic encounters between 2015-2020, and median values were generated. Up to 36,000 samples were genotyped using three versions of the Biobank SNP array offered by Illumina, which underwent harmonization and rigorous quality control (**Supplementary Methods**). Phasing was performed with SHAPEIT [35] and then imputed with the Haplotype Consortium Reference Panel [36] using the Michigan Imputation Server [37]. All SNPs were genotyped or imputed with high quality (r2 values > 0.95).

We performed individual-level analyses on individuals of European ancestry based on self-reported ancestry and genetic PC’s, totaling 25,419 individuals. T2D partitioned polygenic scores (pPSs) for each cluster were generated by multiplying a variant’s genotype dosage by its cluster weight, with only the top-weighted variants included, as defined above. Logistic and linear regression were performed in R v3.6.2, adjusting for age, sex, and PC’s.

## Results

### Ten T2D genetic clusters identified by high-throughput approach

We employed a novel high-throughput pipeline to enable extraction of loci from GWAS summary statistics files and generate a variant-trait association input matrix for clustering analysis (**Figure S1**). Our pipeline started with summary statistics from 13 T2D GWAS studies available in the AMP-CMDKP [15], from which we extracted 21,666 variants associated with T2D reaching genome-wide significance and performed stringent LD-pruning and optimization to identify 324 T2D variants, representing independent T2D risk loci (see **Materials and methods**). We let the T2D genetics guide selection of 64 relevant traits to include in subsequent cluster analysis, such that each trait was significantly associated with at least one T2D variant. Soft clustering of the resulting 324 by 64 variant-trait association matrix was performed using bNMF.

The plurality of bNMF iteration results converged on ten clusters (36.3%), which captured the five clusters identified in our previous work [6] as well as five novel clusters (**Table S5, S6**). The remaining bNMF iterations converged on nested clusters, with 6 clusters in 0.3%, 7 clusters in 1.1%, 8 clusters in 8.3%, 9 clusters in 26.6%, 11 clusters in 22.6%, 12 clusters in 4.4% and 13 clusters in 0.4%. The same six clusters (Beta-cell 1, Beta-cell 2, Proinsulin, Obesity, Lipodystrophy, Liver/Lipid, as described below) were identified in all iterations, based on inspection of constituent variants and traits.

To interpret these clusters, we examined their most highly weighted loci and traits, as well as the aggregate associations of cluster loci with the traits via GWAS pPS (**Table 1, Table S7, Figure S3**). The clusters were named after their most defining traits. Four of the clusters (Beta-cell 1, Beta-cell 2, Proinsulin, and Lipoprotein A) appeared to represent mechanisms of insulin deficiency, with T2D risk-increasing alleles in each cluster collectively associated with reduced fasting insulin and reduced homeostatic assessment of beta cell function (HOMA-B; GWAS pPS *P*-values<0.05). Another five clusters (Obesity, Lipodystrophy, Liver/Lipid, ALP negative, Hyper Insulin Secretion) reflected mechanisms of impaired insulin action, with the T2D risk alleles in these clusters associated with increased fasting insulin and homeostatic assessment of insulin resistance (HOMA-IR; GWAS pPS *P-*values<0.05). The remaining cluster (SHBG) was driven by one T2D allele which was not significantly associated with fasting insulin, but trended toward increased levels (**Table 1, Figure 1, Figure S3, Table S7**).

View this table:
[Table 1.](http://medrxiv.org/content/early/2022/08/05/2022.07.11.22277436/T1)

Table 1. Overview of T2D clusters
Physiological impact, key traits, key loci, suspected mechanism for each cluster. The numbers in parentheses next to cluster names indicate the numbers of top-weighted variants in each of the clusters.

![Figure 1.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/08/05/2022.07.11.22277436/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2022/08/05/2022.07.11.22277436/F1)

Figure 1. Cluster associations with metabolic traits using GWAS
Standardized effect sizes with 95% confidence intervals of cluster pPS-trait associations derived from GWAS summary statistics shown in forest plot. Three metabolic traits (Fasting Insulin, Fasting proinsulin adjusted for fasting insulin, Disposition Index) that help discriminate clusters are displayed. The numbers in the parenthesis next to cluster names indicate the number of variants included in the analysis in each cluster. “All SNPs” include all the variants that are top-weighted in at least one cluster. Filled points indicate *P*-values less than 0.05.

![Figure 2.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/08/05/2022.07.11.22277436/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2022/08/05/2022.07.11.22277436/F2)

Figure 2. Clusters of Type 2 Diabetes Loci
Top-weighted loci and traits in each of the ten clusters are represented in circular plots. Green bars represent top-weighted loci, red bars represent increased traits, and blue bars represent reduced traits in each cluster. A maximum of 35 elements (loci and traits) based on highest weights are displayed in each cluster. Blue outline indicates clusters associated with decreased fasting insulin levels, and red outline indicates clusters associated with increased fasting insulin levels.

Of the four clusters related to insulin deficiency (Beta-cell 1, Beta-cell 2, Proinsulin, Lipoprotein A), Beta-cell 1 and Beta-cell 2 appeared to be a combination of the single Beta-cell cluster in our previous work [6], both including several well-known loci (**Table 1**). In Beta-cell 1, the top-weighted traits were decreased corrected insulin response (CIR) and decreased disposition index (DI), both indicators of reduced pancreatic beta-cell function; the most strongly weighted loci included known beta cell loci *MTNR1B, CDKAL1, HHEX, C2CD4A*, and *SLC30A8* [38–40] (**Table 1, Table S5**). Beta-cell 2 cluster was notable for having increased fasting proinsulin adjusted for fasting insulin, reduced HOMA-B and fasting insulin among the most strongly weighted traits; *TCF7L2, ADCY5, GCK, DGKB/TMEM195*, and *GLIS3* were among the top-weighted loci [41–43] (**Table S5, S6**).

The Beta-cell 1 and Beta-cell 2 clusters differed from each other with regard to the magnitude of glycemic trait effects, with Beta-cell 1 (N loci=63) having more marked association with reduced DI compared to Beta-cell 2 (beta=-0.05, *P*=3.69×10−61 vs beta=-0.03, *P*=9.02×10−9), while Beta-cell 2 (N loci=28) had a more marked association with increased fasting proinsulin adjusted for fasting insulin (beta=0.02, *P*=9.81×10−43 vs beta=0.006, *P*=9.81×10−7) and fasting glucose (beta=0.02, *P*=1.87×10−88 vs beta=0.008, *P*=8.78×10−45), compared to Beta-cell 1 cluster (**Figure 1, Table S7**). Proinsulin is a prohormone precursor to insulin, and elevated fasting proinsulin levels relative to fasting insulin levels indicates defective proinsulin processing, particularly related to beta cell stress [44, 45]. The stronger association with increased proinsulin levels for Beta-cell 2 vs Beta-cell 1 could therefore indicate that Beta-cell 2 relates more specifically to beta cell stress.

The Proinsulin cluster, also identified in our previous work, had top-weighted traits of decreased fasting proinsulin adjusted for fasting insulin and reduced HOMA-B (**Table S5, S6**). The top-weighted loci included distinct signals in the *ARAP1/STARD10* locus; beta cell-selective deletion of *StarD10* in mice has previously been shown to cause impaired insulin secretion [46]. In contrast to the other insulin deficiency clusters, the Proinsulin cluster (N loci=18) was significantly associated with decreased fasting proinsulin adjusted for fasting insulin (GWAS pPS *P*=3.51×10−36) (**Figure 1, Table S7**), potentially indicating a mechanism of lack of proinsulin substrate for insulin synthesis and secretion.

The Lipoprotein A cluster was novel to the present analysis and had a single highly weighted trait, increased serum lipoprotein A (Lp(a)), and a single highly weighted locus, *SLC22A3/LPA* (**Table S5, S6**). *SLC22A3/LPA* has been previously associated with serum Lp(a) levels [47] and contains the gene *LPA* encoding Lp(a). The T2D-risk-increasing allele of Lipoprotein A cluster variant (rs487152) was strongly associated with increased Lp(a) levels in the UK Biobank (GWAS pPS *P*=4.06×10−1586) (**Table S7**), but the underlying mechanism relating to insulin deficiency is unknown.

Of the five clusters related to mechanisms of insulin response (Obesity, Lipodystrophy, Liver/Lipid, Hyper Insulin Secretion, ALP negative), three (Obesity, Lipodystrophy, and Liver/Lipid) were also identified in our previous work, but gained additional loci (and traits) in this expanded analysis.

The Obesity cluster had most-strongly weighted traits of increased body mass index (BMI), waist circumference, percent body fat, and C-reactive protein (CRP), and key genetic signals included the well-known obesity loci *FTO* and *MC4R* [48] (**Table S5, S6**). GWAS pPS for the Obesity cluster (N loci=35) identified significant associations with increased fasting insulin (*P*=7.92×10−22), HOMA-IR (*P*=7.58×10−19), BMI (*P*=1.87×10−1398), percent body fat (*P*=6.94×10−83) and CRP (*P*=6.47×10−260), supporting a mechanism of obesity-mediated insulin resistance.

The Lipodystrophy cluster had top-weighted traits and loci suggestive of “lipodystrophy-like” or fat distribution-mediated insulin resistance as in our prior work [6]; these included decreased adiponectin, HDL cholesterol, and modified Stumvoll insulin sensitivity index (adjusted for age, sex and BMI), and increased triglycerides and waist-hip ratio, as well as top-weighted loci *IRS1, KLF14*, and *PPARG* [49, 50] (**Table S5, S6**). GWAS pPS for the Lipodystrophy cluster (N loci=54) was associated with increased fasting insulin (*P*=3.16×10−43), HOMA-IR (*P*=7.47×10−29) and triglycerides (*P*=1.18×10−612), decreased insulin sensitivity index (*P*=1.84×10−38) and HDL (*P*=5.19×10−535).

The Liver/Lipid cluster was defined by decreased triglycerides and gamma-glutamyl transferase levels, and multiple loci previously connected to hepatic lipid or glycogen metabolism, including *GCKR, HNF1A, PPP1R3B, TOMM40/APOE*, and *PNPLA3* (**Table S5, S6**) [51–56]. GWAS pPS for this cluster (N loci=11) were associated with reduced triglycerides (*P*=3.64×10−181) and interestingly also reduced CRP (*P*=7.75×10−106) and white blood cell count (*P*=1.42×10−49).

The two remaining insulin response clusters were novel (labeled ALP negative and Hyper Insulin Secretion). The ALP negative cluster had decreased alkaline phosphatase (ALP) level as its top-weighed trait, and the *ABO* locus as the top-weighted locus (**Table S5, S6**). GWAS pPS in this cluster (N loci=4) was associated with decreased serum ALP (*P*=1.97×10−1431) and triglycerides (*P*=4.49×10−247). The Hyper Insulin Secretion cluster included top-weighted traits of increased DI and CIR, and loci *PPP1R3B, CNTN2, DTNB, SREBF1*, and *TNF*. The Hyper Insulin Secretion GWAS pPS (N loci=32) was associated with increased CIR (*P*=1.16×10−14), DI (*P*=2.89×10−14), BMI (*P*=1.01×10−26), and reduced HDL (*P*=1.09×10−110) and SHBG (*P*=1.07×10−100).

The final cluster, labeled SHBG, was novel to the current work and not significantly associated with fasting insulin (GWAS pPS *P*=0.36). The cluster was driven by a single trait and locus: decreased SHBG levels and the *SHBG* locus (**Table S5, S6**). GWAS pPS in the SHBG cluster (N loci=1) was significantly associated with reduced SHBG (*P*=1.2×10−1784) and reduced IGF-1 (*P*=4.12×10−13).

### T2D Clusters differ in tissue enrichment including single cell islets

To acquire further evidence for the suspected mechanistic pathways represented by clusters and assess the biological difference between the clusters, we analyzed the top weighted loci in each cluster for enrichment of epigenomic annotations across 28 tissues. The T2D clusters displayed clearly different patterns of tissue enhancer/promoter enrichment (**Figure 3A, Table S8a**). In line with expected mechanisms, the Beta-cell 1, Beta-cell 2, and Proinsulin clusters were significantly enriched for pancreatic islet tissue (FDR<0.05). The Liver/Lipid and ALP negative clusters were significantly enriched in liver tissue (FDR<0.01). The Lipodystrophy cluster was strongly enriched for adipose tissue (FDR<0.01). Additionally, both Beta-cell 1 and 2 clusters had enrichment in adipose and the brain hippocampus (FDR<0.01). The Obesity cluster was most transcriptionally enriched in human epidermal keratinocytes (NHEK) and hASC-t3 pre-adipose cells, both at nominal significance (*P*<0.05, FDR=0.11); of note, a stigmata of insulin resistance commonly seen in obese patients with T2D is acanthosis nigricans, which is hyperpigmentation of skin driven by proliferation of epidermal keratinocytes [57].

![Figure 3.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/08/05/2022.07.11.22277436/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2022/08/05/2022.07.11.22277436/F3)

Figure 3. Enrichment for tissue-specific enhancers in T2D clusters
(A) Heatmap of tissue enhancer/promoter enrichment analysis result. (B) Heatmap of pancreatic islet cell enrichment analysis result. Significance was indicated as follows: \***| FDR < 0.001, ** FDR < 0.01, * FDR < 0.1, • P < 0.05. (C) Comparison of Beta-cell 1 and Beta-cell 2 clusters in fgwas enrichment analysis in functional and stressed beta-cell states shown in a forest plot.

We also interrogated newly available chromatin profiles from 14.3k pancreatic islet cells, which Chiou *et al*. subsetted based on their chromatin profiles [30]. In prior work, the islets were found to have two epigenomic subsets, labeled Beta INShigh and Beta INSlow, indicating high or low insulin gene (*INS)* promoter accessibility; the Beta INShigh islet cells were noted to be enriched for promoter accessibility for genes involved in insulin secretion, whereas the Beta INSlow was enriched for genes involved in stress-induced signaling response [30]. When assessing enrichment of our genetic clusters, we found that our Beta-cell 1 genetic cluster was enriched only in Beta INShigh cells (P=0.0001, FDR=0.0014), whereas our Beta-cell 2 genetic cluster was nominally enriched in both Beta INShigh and Beta INSlow cells (P=0.025, P=0.013, respectively, FDR=0.18 for both), (**Figure 3B, Table S8b**). Further supporting the delineation of the Beta cell loci into two separate sub-pathways, the same trend was observed in our fgwas enrichment analysis: Beta-cell 1 was significantly enriched only in INShigh (ln(enrichment) (95% CI) INShigh 2.32 (1.31-3.12); INSlow -0.36 (−1.79-0.55)) whereas Beta-cell 2 was significantly enriched in both single cell subsets (ln(enrichment) (95% CI), INShigh 1.61 (0.22 - 2.96); INSlow 2.11 (0.73 - 3.46)) (**Figure 3C**). Together these results supported that Beta-cell 1 and Beta-cell 2 clusters relate to distinct physiological mechanisms, with Beta-cell 2 again connected to a stress-induced pancreatic state.

Also of interest, within the pancreas single cell data, the Liver/Lipid cluster was most enriched for alpha cells, (P=0.007, FDR=0.099); alpha cells secrete glucagon, which acts to release glucose from the glycogen stores in the liver, providing further connection between these T2D genetic loci with liver function.

### T2D clusters are differentially associated with clinical traits and outcomes

To assess translation of the clusters to individuals, we generated cluster pPSs in the hospital-based MGB Biobank (N=25,419). We first confirmed that cluster pPSs were associated with expected traits in this study population both in all individuals and in just those with T2D (**Table S11**).

We next tested whether the cluster pPS were associated cardiometabolic clinical outcomes related to T2D using GWAS summary statistics: CAD, CKD, eGFR, hypertension, ischemic stroke, diabetic retinopathy, and diabetic neuropathy (**Table S4, Table S10, Figure 4a, Figure S4a**). All ten T2D clusters were associated with at least one outcome. The GWAS pPS results for eGFR highlighted the utility of cluster-specific scores, with individual clusters having more significant associations than the full set of combined T2D SNPs: increased pPSs for the Liver/Lipid, ALP negative, and SHBG clusters were associated with reduced eGFR (*P*<5×10−4), whereas all cluster T2D SNPs together did not reach Bonferroni-corrected statistical significance (*P*=0.03, **Table S10**). The most significant of these GWAS pPS were replicated using individual-data from MGB Biobank: increased Obesity cluster pPS had risk of hypertension, increased Lipodystrophy cluster pPS and risk of CAD, and increased Liver/Lipid cluster pPS and reduced eGFR, in all individuals with and without adjustment for T2D status (**Figure 4, Table S11**).

![Figure 4.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/08/05/2022.07.11.22277436/F4.medium.gif)

[Figure 4.](http://medrxiv.org/content/early/2022/08/05/2022.07.11.22277436/F4)

Figure 4. Forest plot of cluster associations with outcomes using (A) GWAS and (B) individual-level data from MGB Biobank
(A) Standardized effect sizes with 95% confidence intervals of cluster pPS-outcome associations derived from GWAS summary statistics shown in forest plot. Three metabolic outcomes (T2D, CAD and CKD, all T2D unadjusted) are displayed. The numbers in the parenthesis next to cluster names indicate the number of variants included in the analysis in each cluster. Filled points indicate *P*-values less than 0.05. (B) Associations of pPSs in individuals in the MGB Biobank with clinical outcomes are shown in forest plot. Three outcomes including T2D are displayed.

### Clusters from CAD and CKD share mechanistic pathways with T2D

To identify shared pathways among cardiometabolic outcomes and demonstrate portability of our high-throughput pipeline, we applied the same clustering approach to 219 CAD loci and 70 CKD loci. Five clusters were identified in clustering analysis of 219 CAD loci (ALP negative, Lipoprotein A, HDL negative, Cholesterol and Blood markers increased) and five clusters were identified for the 70 CKD loci (Blood markers increased, Urea increased, Reduced hematopoiesis, Beta-cell opposite and Lipoprotein A) (**Tables S12-15, Figure S5**). Based on inspection of constituent variants and traits in the clusters of T2D, CAD, and CKD, one cluster, Lipoprotein A, was shared by all three diseases. Similarly, the ALP negative cluster was shared between T2D and CAD, and the Blood markers increased cluster between CAD and CKD.

## Discussion

Novel approaches are needed to connect the currently identified hundreds of T2D genetic loci to disease pathophysiology and also accommodate the rapid pace of new loci discovery. Here, we describe expanded clustering of T2D variants, using a high-throughput pipeline for extracting and preprocessing variants from multiple GWAS datasets and generating a variant-trait association matrix. The resulting matrix consisted of 324 T2D genetic variants and 64 diabetes-related metabolic traits from publicly available GWAS datasets. By applying bNMF soft clustering to this matrix, we identified ten robust clusters of T2D variants, representing biologically meaningful mechanistic pathways.

Among the ten clusters, we replicated the five identified in our previous work of 94 T2D variants (Beta-cell, Proinsulin, Obesity, Lipodystrophy, Liver/Lipid) [6], with the Beta-cell cluster now subdivided into two distinct clusters, and also identified four additional novel clusters related to pronounced insulin secretion, levels of alkaline phosphatase, lipoprotein-A, and sex hormone-binding globulin. In contrast to our prior work, which involved manual curation of loci published in GWAS manuscripts to generate the input list of variants, the current approach allowed for use of uncurated GWAS summary statistics files and included additional newly available datasets, more than tripling the number of input genetic loci. Thus, replication of the previously identified five clusters provides strong validation of this high-throughput approach, with the newly identified clusters driven by traits or loci not available in the prior analysis.

Three of the ten T2D clusters identified in this work (Beta-cell 1, Beta-cell 2, and Proinsulin) clearly related to pancreatic beta-cell function, with the two Beta-cell clusters differing from the Proinsulin cluster with regard to the direction of association with fasting proinsulin adjusted for fasting insulin. All three clusters were enriched in pancreatic islet tissue enhancers and promoters in the epigenomics analysis. Additionally, loci in the Beta-cell 1 cluster were significantly enriched for a subset of single beta cells predicted independent of our work based on RNA transcript levels to represent a normal state of pancreatic beta cell function, whereas loci in Beta-cell 2 cluster had a unique signal of enrichment for single beta cells predicted to be in a stressed state [30]; these functional distinctions between Beta-cell 1 and 2 supported our independent approach of phenotypically informed clustering T2D loci.

Three other T2D genetic clusters (Obesity, Lipodystrophy, Liver/Lipid) replicated findings from our prior work related to pathways of insulin resistance, gaining additional loci and traits compared to the prior analysis. Loci in these three clusters were most enriched for enhancers in tissues for the suspected mechanisms: pre-adipocytes, adipocytes, and liver tissue, respectively. The distinction between fat accumulation in the Obesity cluster and abnormal fat compartmentalization in the Lipodystrophy cluster may be supported by the differential enhancer enrichment shown for different developmental stages of the same adipocyte lineage.

In addition to a second Beta-cell cluster, there were four newly identified T2D genetic clusters from this work: ALP negative (containing the *ABO* locus), Lipoprotein A, SHBG, and Hyper Insulin Secretion.

The ALP negative cluster seen for T2D and CAD was driven by reduced serum alkaline phosphatase levels and the *ABO* locus. Isoforms of alkaline phosphatase have been shown to vary in level by blood group [58]. The *ABO* locus and blood type have previously been connected to T2D [59–61] and CAD [62] risk, but the causal mechanisms are not fully understood.

The Lipoprotein A cluster seen for T2D, CAD, and CKD all included the locus (*SLC22A3*/*LPA* tagged by rs487152) and biomarker Lp(a), pointing to a genetic pathway leading to increased Lp(a) levels and also increased risk of T2D, CAD, and CKD. The relationship between Lp(a) and cardiometabolic disease is complex, and genetic interrogation of *LPA* has been complicated by the fact that plasma concentration of Lp(a) is influenced by kringle IV type 2 (KIV-2) repeats in addition to other genetic variation [63]. While epidemiological studies have connected elevated Lp(a) levels with increased risk of CAD and CKD [64, 65], an inverse association has been reported for T2D [66, 67]. Our genetic findings for T2D therefore indicate that there are likely multiple pathways impacting Lp(a) level that may have differential effects on T2D risk.

The SHBG cluster was also driven by a single locus and biomarker. Our results point to a genetic pathway whereby alteration of the *SHBG* locus leads to reduced SHBG levels and increased T2D risk, which was consistent with previous epidemiological and genetic studies indicating that low circulating levels of SHBG were causally related to increased risk of T2D in women and men [68, 69].

We assessed the impact of cluster pPSs in individuals, finding that individuals with increased cluster pPS had significant associations with clinical traits and outcomes (**Figure 4, Table S9, S11**), supporting prior findings for the original five clusters [70]. These results point to marked heterogeneity in T2D genetic associations that could be missed without delineation of loci into clusters and also suggest important physiological differences between clusters. At the same time, the effect sizes of the present pPSs on clinical outcomes were likely too small to be of clinical utility at the individual-level.

The strengths of this study include the high-throughput approach for preprocessing variants and traits from multiple GWAS datasets in a semi-automated way. This method can also be readily applied to other diseases beyond T2D to identify key pathways and code has been made freely available. We included here application of the pipeline to CAD and CKD, demonstrating transferability of the approach and potential shared pathways among the three cardiometabolic outcomes. Limitations include clustering of only available phenotypes from GWAS. It is possible that additional pathways exist that are not captured using the set of traits included in the analysis. Additionally, due to methodological limitations and data availability we have focused on GWAS datasets from populations of European ancestry, although we are actively pursuing application of this method in non-European populations through additional efforts. It is worth noting that bNMF generates weights for all included elements in the matrix, and it is not known how best determine a cut-off threshold for cluster membership; we have applied a reasonable strategy to maximize signal to noise. Finally, whether the associations of specific genetically derived clusters with metabolic traits and outcomes remain constant throughout the disease course has not been examined in this cross-sectional analysis.

In summary, we have identified ten robust genetic clusters pointing to mechanistic pathways of T2D using a high-throughput clustering pipeline of GWAS summary statistics. These clusters displayed tissue-specific enrichment patterns even within single cell pancreatic tissue and could be used to generate pPSs that stratify patients genetically with distinct associations with clinical features and cardiometabolic outcomes. We demonstrate that our approach can be applied to other complex diseases, with identification of shared genetic pathways between T2D, CAD, and CKD. Thus, we contribute to further delineation of cardiometabolic disease genetic pathways using a data-driven approach informed by physiology.

## Supporting information

Supplemental Tables [[supplements/277436_file02.xlsx]](pending:yes)

Supplemental Figures [[supplements/277436_file03.docx]](pending:yes)

Supplemental Methods [[supplements/277436_file04.docx]](pending:yes)

## Data Availability

All data produced in the present study are available upon reasonable request to the authors

## Description of Supplemental Data

Supplemental Data include five figures and fifteen tables.

## Declaration of Interests

The authors declare no competing interests.

## Data and Code Availability

Code for variant pre-processing, bNMF clustering, and basic visualizations is available at [https://github.com/gwas-partitioning/bnmf-clustering](https://github.com/gwas-partitioning/bnmf-clustering).

## Funding

This work was supported by FNIH RFP-13 and the MGH Transformative Scholars Award.

## Contribution Statement

MSU, JCF, MG, and JC conceived the research question. MSU, JK, MG, KS, KEW, JC, and KG conceived the methodology, which included implementation of the clustering computational pipeline. HK, JC, MG, TM, JMM and MSU curated the data. HK, KS and JC conducted the analysis and visualized the results. HK, JC, and MSU wrote the initial draft of the paper and incorporated co-author comments. KEW, KS, JBC, TM, MG, JMM, SK, JCF, KG and AKM provided feedback on the analysis, and critically reviewed the manuscript. All co-authors approved the final version of the paper. MSU is the guarantor of this work and, as such, had full access to all the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis.

## Web Resources

Interactive results are viewable on the Common Metabolic Disease Knowledge Portal ([https://hugeamp.org/](https://hugeamp.org/)).

## Footnotes

*   Methods section has been shortened and Supplementary methods has been added; Supplementary figures, Supplementary tables, and manuscript main text related to chronic kidney disease (CKD) clusters have been updated.

*   Received July 11, 2022.
*   Revision received August 5, 2022.
*   Accepted August 5, 2022.


*   © 2022, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial 4.0 International), CC BY-NC 4.0, as described at [http://creativecommons.org/licenses/by-nc/4.0/](http://creativecommons.org/licenses/by-nc/4.0/)

## References

1.  1.Redondo MJ, Hagopian WA, Oram R, et al. (2020) The clinical consequences of heterogeneity within and between different diabetes types. Diabetologia 63(10):2040–2048
    
    

2.  2.Mahajan A, Wessel J, Willems SM, et al. (2018) Refining the accuracy of validated target identification through coding variant fine-mapping in type 2 diabetes. Nat Genet 50(4):559–571
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0084-1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29632382&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 

3.  3.Mahajan A, Taliun D, Thurner M, et al. (2018) Fine-mapping type 2 diabetes loci to single-variant resolution using high-density imputation and islet-specific epigenome maps. Nat Genet 50(11):1505–1513
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0241-6&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30297969&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 

4.  4.Vujkovic M, Keaton JM, Lynch JA, et al. (2020) Discovery of 318 new risk loci for type 2 diabetes and related vascular outcomes among 1.4 million participants in a multi-ancestry meta-analysis. Nat Genet 52(7):680–691
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-020-0637-y&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 

5.  5.Mahajan A, Spracklen CN, Zhang W, et al. (2022) Multi-ancestry genetic study of type 2 diabetes highlights the power of diverse populations for discovery and translation. Nat Genet 54(5):560–572
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-022-01058-3&link_type=DOI) 

6.  6.Udler MS, Kim J, von Grotthuss M, et al. (2018) Type 2 diabetes genetic loci informed by multi-trait associations point to disease mechanisms and subtypes: A soft clustering analysis. PLoS Med 15(9):e1002654
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 

7.  7.Goodarzi MO, Palmer ND, Cui J, et al. (2020) Classification of type 2 diabetes genetic variants and a novel genetic risk score association with insulin clearance. J Clin Endocrinol Metab 105(4):1251–1260
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 

8.  8.Dimas AS, Lagou V, Barker A, et al. (2014) Impact of type 2 diabetes susceptibility variants on quantitative glycemic traits reveals mechanistic heterogeneity. Diabetes 63(6):2158–2171
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZGlhYmV0ZXMiO3M6NToicmVzaWQiO3M6OToiNjMvNi8yMTU4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDgvMDUvMjAyMi4wNy4xMS4yMjI3NzQzNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

9.  9.Udler MS, McCarthy MI, Florez JC, Mahajan A (2019) Genetic Risk Scores for Diabetes Diagnosis and Precision Medicine. Endocr Rev 40(6):1500–1520
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1210/er.2019-00088&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 

10. 10.Scott RA, Scott LJ, Mägi R, et al. (2017) An expanded genome-wide association study of type 2 diabetes in Europeans. Diabetes 66(11):2888–2902
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZGlhYmV0ZXMiO3M6NToicmVzaWQiO3M6MTA6IjY2LzExLzI4ODgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wOC8wNS8yMDIyLjA3LjExLjIyMjc3NDM2LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

11. 11.Kurki MI, Karjalainen J, Palta P, et al. (2022) FinnGen: Unique genetic insights from combining isolated population and national health register data. bioRxiv
    
    

12. 12.Bonàs-Guarch S, Guindo-Martínez M, Miguel-Escalada I, et al. (2018) Re-analysis of public genetic data reveals a rare X-chromosomal variant associated with type 2 diabetes. Nat Commun 9(1):321
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-017-02380-9&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 

13. 13.Flannick J, Mercader JM, Fuchsberger C, et al. (2019) Exome sequencing of 20,791 cases of type 2 diabetes and 24,440 controls. Nature 570(7759):71–76
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-019-1231-2&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31118516&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 

14. 14.Guindo-Martínez M, Amela R, Bonàs-Guarch S, et al. (2021) The impact of non-additive genetic associations on age-related complex diseases. Nat Commun 12(1):2436
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-021-21952-4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33893285&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 

15. 15.Human Genetics Knowledge Portal - Home. [https://hugeamp.org/](https://hugeamp.org/). Accessed 19 Mar 2021
    
    

16. 16.UK biobank —. In: Neale lab. [http://www.nealelab.is/uk-biobank](http://www.nealelab.is/uk-biobank). Accessed 18 Jul 2022
    
    

17. 17.van der Harst P, Verweij N (2018) Identification of 64 novel genetic loci provides an expanded view on the genetic architecture of coronary artery disease. Circ Res 122(3):433–443
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTA6ImNpcmNyZXNhaGEiO3M6NToicmVzaWQiO3M6OToiMTIyLzMvNDMzIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDgvMDUvMjAyMi4wNy4xMS4yMjI3NzQzNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

18. 18.Nikpay M, Goel A, Won H-H, et al. (2015) A comprehensive 1,000 Genomes-based genome-wide association meta-analysis of coronary artery disease. Nat Genet 47(10):1121–1130
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3396&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26343387&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 

19. 19. Myocardial Infarction Genetics and  CARDIoGRAM Exome Consortia Investigators, Stitziel NO, Stirrups KE, et al. (2016) Coding variation in ANGPTL4, LPL, and SVEP1 and the risk of coronary disease. N Engl J Med 374(12):1134–1144
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa1507652&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26934567&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 

20. 20.Wuttke M, Li Y, Li M, et al. (2019) A catalog of genetic loci associated with kidney function from analyses of a million individuals. Nat Genet 51(6):957–972
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-019-0407-x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31152163&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 

21. 21.Salem RM, Todd JN, Sandholm N, et al. (2019) Genome-wide association study of diabetic kidney disease highlights biology involved in glomerular basement membrane collagen. J Am Soc Nephrol 30(10):2000–2016
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiam5lcGhyb2wiO3M6NToicmVzaWQiO3M6MTA6IjMwLzEwLzIwMDAiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wOC8wNS8yMDIyLjA3LjExLjIyMjc3NDM2LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

22. 22.van Zuydam NR, Ahlqvist E, Sandholm N, et al. (2018) A genome-wide association study of diabetic kidney disease in subjects with type 2 diabetes. Diabetes 67(7):1414–1427
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2337/db18-1414-P&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29703844&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 

23. 23.Locke AE, Steinberg KM, Chiang CWK, et al. (2019) Exome sequencing of Finnish isolates enhances rare-variant association power. Nature 572(7769):323–328
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-019-1457-z&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 

24. 24.Companion R package for the guide Doing Meta-Analysis in R. [http://dmetar.protectlab.org/](http://dmetar.protectlab.org/). Accessed 18 Jul 2022
    
    

25. 25.Wakefield J (2007) A Bayesian measure of the probability of false discovery in genetic epidemiology studies. Am J Hum Genet 81(2):208–227
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/519024&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17668372&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000248540400002&link_type=ISI) 

26. 26.Ernst J, Kellis M (2012) ChromHMM: automating chromatin-state discovery and characterization. Nat Methods 9(3):215–216
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nmeth.1906&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22373907&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000300890400006&link_type=ISI) 

27. 27.Varshney A, Scott LJ, Welch RP, et al. (2017) Genetic regulatory signatures underlying islet gene expression and type 2 diabetes. Proc Natl Acad Sci U S A 114(9):2301–2306
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMDoiMTE0LzkvMjMwMSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzA4LzA1LzIwMjIuMDcuMTEuMjIyNzc0MzYuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

28. 28.Zhang K, Hocker JD, Miller M, et al. (2021) A cell atlas of chromatin accessibility across 25 adult human tissues. Cold Spring Harbor Laboratory 2021.02.17.431699
    
    

29. 29.Chiou J, Geusz RJ, Okino M-L, et al. (2021) Interpreting type 1 diabetes risk with genetics and single-cell epigenomics. Nature 1–5
    
    

30. 30.Chiou J, Zeng C, Cheng Z, et al. (2021) Single-cell chromatin accessibility identifies pancreatic islet cell type- and state-specific regulatory programs of diabetes risk. Nat Genet 53(4):455–466
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 

31. 31.Pickrell JK (2014) Joint Analysis of Functional Genomic Data and Genome-wide Association Studies of 18 Human Traits. Am J Hum Genet 94(4):559
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2014.03.004&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24702953&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 

32. 32.Smoller JW, Karlson EW, Green RC, et al. (2016) An eMERGE Clinical Center at Partners Personalized Medicine. Journal of Personalized Medicine 6(1):5
    
    

33. 33.Karlson EW, Boutin NT, Hoffnagle AG, Allen NL (2016) Building the Partners HealthCare Biobank at Partners Personalized Medicine: Informed Consent, Return of Research Results, Recruitment Lessons and Operational Considerations. J Pers Med 6(1). [https://doi.org/10.3390/jpm6010002](https://doi.org/10.3390/jpm6010002)
    
    

34. 34.Yu S, Liao KP, Shaw SY, et al. (2015) Toward high-throughput phenotyping: unbiased automated feature extraction and selection from knowledge sources. J Am Med Inform Assoc 22(5):993–1000
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/jamia/ocv034&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25929596&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 

35. 35.Delaneau O, Zagury J-F, Marchini J (2013) Improved whole-chromosome phasing for disease and population genetic studies. Nat Methods 10(1):5–6
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nchembio.1414&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23269371&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000312810100003&link_type=ISI) 

36. 36.McCarthy S, Das S, Kretzschmar W, et al. (2016) A reference panel of 64,976 haplotypes for genotype imputation. Nat Genet 48(10):1279–1283
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3643&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27548312&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 

37. 37.Das S, Forer L, Schönherr S, et al. (2016) Next-generation genotype imputation service and methods. Nat Genet 48(10):1284–1287
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3656&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27571263&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 

38. 38.Rosengren AH, Braun M, Mahdi T, et al. (2012) Reduced Insulin Exocytosis in Human Pancreatic β-Cells With Gene Variants Linked to Type 2 Diabetes. Diabetes 61(7):1726–1733
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZGlhYmV0ZXMiO3M6NToicmVzaWQiO3M6OToiNjEvNy8xNzI2IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDgvMDUvMjAyMi4wNy4xMS4yMjI3NzQzNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

39. 39.Harder MN, Ribel-Madsen R, Justesen JM, et al. (2013) Type 2 diabetes risk alleles near BCAR1 and in ANK1 associate with decreased β-cell function whereas risk alleles near ANKRD55 and GRB14 associate with decreased insulin sensitivity in the Danish Inter99 cohort. J Clin Endocrinol Metab 98(4):E801–6
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1210/jc.2012-4169&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23457408&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 

40. 40.Xu J, Wijesekara N, Regeenes R, et al. (2021) Pancreatic β cell-selective zinc transporter 8 insufficiency accelerates diabetes associated with islet amyloidosis. JCI Insight 6(10). [https://doi.org/10.1172/jci.insight.143037](https://doi.org/10.1172/jci.insight.143037)
    
    

41. 41.Zhou Y, Park S-Y, Su J, et al. (2014) TCF7L2 is a master regulator of insulin production and processing. Hum Mol Genet 23(24):6419–6431
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/hmg/ddu359&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25015099&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 

42. 42.Hodson DJ, Mitchell RK, Marselli L, et al. (2014) ADCY5 couples glucose to insulin secretion in human islets. Diabetes 63(9):3009–3021
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZGlhYmV0ZXMiO3M6NToicmVzaWQiO3M6OToiNjMvOS8zMDA5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDgvMDUvMjAyMi4wNy4xMS4yMjI3NzQzNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

43. 43.Boesgaard TW, Grarup N, Jørgensen T, et al. (2010) Variants at DGKB/TMEM195, ADRA2A, GLIS3 and C2CD4B loci are associated with reduced glucose-stimulated beta cell function in middle-aged Danish people. Diabetologia 53(8):1647–1655
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00125-010-1753-5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20419449&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000279193800016&link_type=ISI) 

44. 44.Mezza T, Ferraro PM, Sun VA, et al. (2018) Increased β-Cell Workload Modulates Proinsulin-to-Insulin Ratio in Humans. Diabetes 67(11):2389–2396
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2337/db18-2389-PUB&link_type=DOI) 

45. 45.Røder ME, Porte D Jr., Schwartz RS, Kahn SE (1998) Disproportionately elevated proinsulin levels reflect the degree of impaired B cell secretory capacity in patients with noninsulin-dependent diabetes mellitus. J Clin Endocrinol Metab 83(2):604–608
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1210/jc.83.2.604&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9467581&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000071823900051&link_type=ISI) 

46. 46.Carrat GR, Hu M, Nguyen-Tu M-S, et al. (2017) Decreased STARD10 Expression Is Associated with Defective Insulin Secretion in Humans and Mice. Am J Hum Genet 100(2):238–256
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2017.01.011&link_type=DOI) 

47. 47.Qi Q, Workalemahu T, Zhang C, Hu FB, Qi L (2011) Genetic variants, plasma lipoprotein(a) levels, and risk of cardiovascular morbidity and mortality among two prospective cohorts of type 2 diabetes. Eur Heart J 33(3):325–334
    
    

48. 48.Choquet H, Meyre D (2011) Genetics of Obesity: What have we Learned? Curr Genomics 12(3):169–179
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2174/138920211795677895&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22043165&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000290618200003&link_type=ISI) 

49. 49.Yaghootkar H, Scott RA, White CC, et al. (2014) Genetic Evidence for a Normal-Weight “Metabolically Obese” Phenotype Linking Insulin Resistance, Hypertension, Coronary Artery Disease, and Type 2 Diabetes. Diabetes 63(12):4369–4377
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZGlhYmV0ZXMiO3M6NToicmVzaWQiO3M6MTA6IjYzLzEyLzQzNjkiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wOC8wNS8yMDIyLjA3LjExLjIyMjc3NDM2LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

50. 50.Yaghootkar H, Lotta LA, Tyrrell J, et al. (2016) Genetic Evidence for a Link Between Favorable Adiposity and Lower Risk of Type 2 Diabetes, Hypertension, and Heart Disease. Diabetes 65(8):2448–2460
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZGlhYmV0ZXMiO3M6NToicmVzaWQiO3M6OToiNjUvOC8yNDQ4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDgvMDUvMjAyMi4wNy4xMS4yMjI3NzQzNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

51. 51.Speliotes EK, Yerges-Armstrong LM, Wu J, et al. (2011) Genome-wide association analysis identifies variants associated with nonalcoholic fatty liver disease that have distinct effects on metabolic traits. PLoS Genet 7(3):e1001324
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgen.1001324&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21423719&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 

52. 52.Mahdessian H, Taxiarchis A, Popov S, et al. (2014) TM6SF2 is a regulator of liver fat metabolism influencing triglyceride secretion and hepatic lipid droplet content. Proc Natl Acad Sci U S A 111(24):8913–8918
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTExLzI0Lzg5MTMiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wOC8wNS8yMDIyLjA3LjExLjIyMjc3NDM2LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

53. 53.Kozlitina J, Smagris E, Stender S, et al. (2014) Exome-wide association study identifies a TM6SF2 variant that confers susceptibility to nonalcoholic fatty liver disease. Nat Genet 46(4):352–356
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.2901&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24531328&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 

54. 54.Smagris E, Gilyard S, BasuRay S, Cohen JC, Hobbs HH (2016) Inactivation of Tm6sf2, a gene defective in fatty liver disease, impairs lipidation but not secretion of very low density lipoproteins. J Biol Chem 291(20):10659–10676
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamJjIjtzOjU6InJlc2lkIjtzOjEyOiIyOTEvMjAvMTA2NTkiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wOC8wNS8yMDIyLjA3LjExLjIyMjc3NDM2LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

55. 55.Raimondo A, Rees MG, Gloyn AL (2015) Glucokinase regulatory protein: complexity at the crossroads of triglyceride and glucose metabolism. Curr Opin Lipidol 26(2):88–95
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/MOL.0000000000000155&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25692341&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 

56. 56.Smagris E, BasuRay S, Li J, et al. (2015) Pnpla3I148M knockin mice accumulate PNPLA3 on lipid droplets and develop hepatic steatosis. Hepatology 61(1):108–118
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/hep.27242&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24917523&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 

57. 57.Higgins SP, Freemark M, Prose NS (2008) Acanthosis nigricans: a practical approach to evaluation and management. Dermatol Online J 14(9):2
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19061584&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 

58. 58.Domar U, Hirano K, Stigbrand T (1991) Serum levels of human alkaline phosphatase isozymes in relation to blood groups. Clin Chim Acta 203(2–3):305–313
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0009-8981(91)90303-T&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1777990&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 

59. 59.Li-Gao R, Carlotti F, de Mutsert R, et al. (2019) Genome-Wide Association Study on the Early-Phase Insulin Response to a Liquid Mixed Meal: Results From the NEO Study. Diabetes 68(12):2327–2336
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2337/db19-2327-PUB&link_type=DOI) 

60. 60.Qi L, Cornelis MC, Kraft P, et al. (2010) Genetic variants in ABO blood group region, plasma soluble E-selectin levels and risk of type 2 diabetes. Hum Mol Genet 19(9):1856–1862
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/hmg/ddq057&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20147318&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000276526700019&link_type=ISI) 

61. 61.Fagherazzi G, Gusto G, Clavel-Chapelon F, Balkau B, Bonnet F (2015) ABO and Rhesus blood groups and risk of type 2 diabetes: evidence from the large E3N cohort study. Diabetologia 58(3):519–522
    
    

62. 62.Chen Z, Yang S-H, Xu H, Li J-J (2016) ABO blood group system and the coronary artery disease: an updated systematic review and meta-analysis. Sci Rep 6(1):23250
    
    

63. 63.Tolbus A, Mortensen MB, Nielsen SF, Kamstrup PR, Bojesen SE, Nordestgaard BG (2017) Kringle IV Type 2, Not Low Lipoprotein(a), as a Cause of Diabetes: A Novel Genetic Approach Using SNPs Associated Selectively with Lipoprotein(a) Concentrations or with Kringle IV Type 2 Repeats. Clin Chem 63(12):1866–1876
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiY2xpbmNoZW0iO3M6NToicmVzaWQiO3M6MTA6IjYzLzEyLzE4NjYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wOC8wNS8yMDIyLjA3LjExLjIyMjc3NDM2LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

64. 64.Tipping RW, Ford CE, Simpson LM, et al. (2009) Lipoprotein(a) concentration and the risk of coronary heart disease, stroke, and nonvascular mortality. JAMA 302:412–423
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2009.1063&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19622820&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000268139200024&link_type=ISI) 

65. 65.Xuan L, Wang T, Dai H, et al. (2020) Serum lipoprotein (a) associates with a higher risk of reduced renal function: a prospective investigation. J Lipid Res 61(10):1320–1327
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamxyIjtzOjU6InJlc2lkIjtzOjEwOiI2MS8xMC8xMzIwIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDgvMDUvMjAyMi4wNy4xMS4yMjI3NzQzNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

66. 66.Mora S, Kamstrup PR, Rifai N, Nordestgaard BG, Buring JE, Ridker PM (2010) Lipoprotein(a) and risk of type 2 diabetes. Clin Chem 56(8):1252–1260222.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiY2xpbmNoZW0iO3M6NToicmVzaWQiO3M6OToiNTYvOC8xMjUyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDgvMDUvMjAyMi4wNy4xMS4yMjI3NzQzNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

67. 67.Habib SS, Aslam M, Shah SFA, Naveed AK (2009) Lipoproteína (a) está associada com níveis basais de insulina em pacientes com Diabetes Mellitus tipo 2. Arq Bras Cardiol 93(1). [https://doi.org/10.1590/s0066-782x2009000700006](https://doi.org/10.1590/s0066-782x2009000700006)
    
    

68. 68.Ding EL, Song Y, Manson JE, et al. (2009) Sex hormone-binding globulin and risk of type 2 diabetes in women and men. N Engl J Med 361(12):1152–1163
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa0804381&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19657112&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000269828000006&link_type=ISI) 

69. 69.Perry JRB, Weedon MN, Langenberg C, et al. (2010) Genetic evidence that raised sex hormone binding globulin (SHBG) levels reduce the risk of type 2 diabetes. Hum Mol Genet 19(3):535–544
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/hmg/ddp522&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19933169&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F08%2F05%2F2022.07.11.22277436.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000273227200013&link_type=ISI) 

70. 70.DiCorpo D, LeClair J, Cole JB, et al. (2022) Type 2 diabetes partitioned polygenic scores associate with disease outcomes in 454,193 individuals across 13 Cohorts. Diabetes Care 45(3):674–683