Niels Grarup

Meta-analysis of genome-wide association data and large-scale replication identifies additional susceptibility loci for type 2 diabetes

Genome-wide association (GWA) studies have identified multiple loci at which common variants modestly but reproducibly influence risk of type 2 diabetes (T2D). Established associations to common and rare variants explain only a small proportion of the heritability of T2D. As previously published analyses had limited power to identify variants with modest effects, we carried out meta-analysis of three T2D GWA scans comprising 10,128 individuals of European descent and ∼2.2 million SNPs (directly genotyped and imputed), followed by replication testing in an independent sample with an effective sample size of up to 53,975. We detected at least six previously unknown loci with robust evidence for association, including the JAZF1 (P = 5.0 × 10−14), CDC123-CAMK1D (P = 1.2 × 10−10), TSPAN8-LGR5 (P = 1.1 × 10−9), THADA (P = 1.1 × 10−9), ADAMTS9 (P = 1.2 × 10−8) and NOTCH2 (P = 4.1 × 10−8) gene regions. Our results illustrate the value of large discovery and follow-up samples for gaining further insights into the inherited basis of T2D.

Richness of human gut microbiome correlates with metabolic markers

We are facing a global metabolic health crisis provoked by an obesity epidemic. Here we report the human gut microbial composition in a population sample of 123 non-obese and 169 obese Danish individuals. We find two groups of individuals that differ by the number of gut microbial genes and thus gut bacterial richness. They contain known and previously unknown bacterial species at different proportions; individuals with a low bacterial richness (23% of the population) are characterized by more marked overall adiposity, insulin resistance and dyslipidaemia and a more pronounced inflammatory phenotype when compared with high bacterial richness individuals. The obese individuals among the lower bacterial richness group also gain more weight over time. Only a few bacterial species are sufficient to distinguish between individuals with high and low bacterial richness, and even between lean and obese participants. Our classifications based on variation in the gut microbiome identify subsets of individuals in the general white adult population who may be at increased risk of progressing to adiposity-associated co-morbidities.

Resequencing of 200 human exomes identifies an excess of low-frequency non-synonymous coding variants

Targeted capture combined with massively parallel exome sequencing is a promising approach to identify genetic variants implicated in human traits. We report exome sequencing of 200 individuals from Denmark with targeted capture of 18,654 coding genes and sequence coverage of each individual exome at an average depth of 12-fold. On average, about 95% of the target regions were covered by at least one read. We identified 121,870 SNPs in the sample population, including 53,081 coding SNPs (cSNPs). Using a statistical method for SNP calling and an estimation of allelic frequencies based on our population data, we derived the allele frequency spectrum of cSNPs with a minor allele frequency greater than 0.02. We identified a 1.8-fold excess of deleterious, non-syonomyous cSNPs over synonymous cSNPs in the low-frequency range (minor allele frequencies between 2% and 5%). This excess was more pronounced for X-linked SNPs, suggesting that deleterious substitutions are primarily recessive.

Loss-of-function mutations in SLC30A8 protect against type 2 diabetes

Loss-of-function mutations protective against human disease provide in vivo validation of therapeutic targets, but none have yet been described for type 2 diabetes (T2D). Through sequencing or genotyping of ∼150,000 individuals across 5 ancestry groups, we identified 12 rare protein-truncating variants in SLC30A8, which encodes an islet zinc transporter (ZnT8) and harbors a common variant (p.Trp325Arg) associated with T2D risk and glucose and proinsulin levels. Collectively, carriers of protein-truncating variants had 65% reduced T2D risk (P = 1.7 × 10−6), and non-diabetic Icelandic carriers of a frameshift variant (p.Lys34Serfs*50) demonstrated reduced glucose levels (−0.17 s.d., P = 4.6 × 10−4). The two most common protein-truncating variants (p.Arg138* and p.Lys34Serfs*50) individually associate with T2D protection and encode unstable ZnT8 proteins. Previous functional study of SLC30A8 suggested that reduced zinc transport increases T2D risk, and phenotypic heterogeneity was observed in mouse Slc30a8 knockouts. In contrast, loss-of-function mutations in humans provide strong evidence that SLC30A8 haploinsufficiency protects against T2D, suggesting ZnT8 inhibition as a therapeutic strategy in T2D prevention.

Studies of association of variants near the HHEX, CDKN2A/B and IGF2BP2 genes with type 2 diabetes and impaired insulin release in 10,705 Danish subjects – validation and extension of genome-wide association studies

OBJECTIVE— In the present study, we aimed to validate the type 2 diabetes susceptibility alleles identified in six recent genome-wide association studies in the HHEX/KIF11/IDE (rs1111875), CDKN2A/B (rs10811661), and IGF2BP2 (rs4402960) loci, as well as the intergenic rs9300039 variant. Furthermore, we aimed to characterize quantitative metabolic risk phenotypes of the four variants. RESEARCH DESIGN AND METHODS— The variants were genotyped in the population-based Inter99 cohort (n = 5,970), the ADDITION Study (n = 1,626), a population-based sample of young healthy subjects (n = 377), and in additional type 2 diabetic case (n = 2,111) and glucose-tolerant (n = 521) subjects. The case-control studies involved a total of 4,089 type 2 diabetic patients and 5,043 glucose-tolerant control subjects. RESULTS— We validated association of variants near HHEX/KIF11/IDE, CDKN2A/B, and IGF2BP2 with type 2 diabetes. Interestingly, in middle-aged people, the rs1111875 C-allele of HHEX/KIF11/IDE strongly associated with lower acute insulin response during an oral glucose tolerance test (P = 6 × 10−7). In addition, decreased insulin release following intravenous tolbutamide injection was observed in young healthy subjects (P = 0.02). Also, a reduced insulin release was observed for the CDKN2A/B rs10811661 T-allele after both oral and intravenous glucose challenges (P = 0.001 and P = 0.009, respectively). CONCLUSIONS— We validate that variants in the proximity of the HHEX/KIF11/IDE, CDKN2A/B, and IFG2BP2 loci associate with type 2 diabetes. Importantly, variations within the HHEX/KIF11/IDE and CDKN2A/B loci confer impaired glucose- and tolbutamide-induced insulin release in middle-aged and young healthy subjects, suggesting a role for these variants in the pathogenesis of pancreatic β-cell dysfunction.

A genome-wide association study in the Japanese population identifies susceptibility loci for type 2 diabetes at UBE2E2 and C2CD4A-C2CD4B

We conducted a genome-wide association study of type 2 diabetes (T2D) using 459,359 SNPs in a Japanese population with a three-stage study design (stage 1, 4,470 cases and 3,071 controls; stage 2, 2,886 cases and 3,087 controls; stage 3, 3,622 cases and 2,356 controls). We identified new associations in UBE2E2 on chromosome 3 and in C2CD4A-C2CD4B on chromosome 15 at genome-wide significant levels (rs7612463 in UBE2E2, combined P = 2.27 × 10−9; rs7172432 in C2CD4A-C2CD4B, combined P = 3.66 × 10−9). The association of these two loci with T2D was replicated in other east Asian populations. In the European populations, the C2CD4A-C2CD4B locus was significantly associated with T2D, and a combined analysis of all populations gave P = 8.78 × 10−14, whereas the UBE2E2 locus did not show association to T2D. In conclusion, we identified two new loci at UBE2E2 and C2CD4A-C2CD4B associated with susceptibility to T2D.

Impact of type 2 diabetes susceptibility variants on quantitative glycemic traits reveals mechanistic heterogeneity

Patients with established type 2 diabetes display both β-cell dysfunction and insulin resistance. To define fundamental processes leading to the diabetic state, we examined the relationship between type 2 diabetes risk variants at 37 established susceptibility loci, and indices of proinsulin processing, insulin secretion, and insulin sensitivity. We included data from up to 58,614 nondiabetic subjects with basal measures and 17,327 with dynamic measures. We used additive genetic models with adjustment for sex, age, and BMI, followed by fixed-effects, inverse-variance meta-analyses. Cluster analyses grouped risk loci into five major categories based on their relationship to these continuous glycemic phenotypes. The first cluster (PPARG, KLF14, IRS1, GCKR) was characterized by primary effects on insulin sensitivity. The second cluster (MTNR1B, GCK) featured risk alleles associated with reduced insulin secretion and fasting hyperglycemia. ARAP1 constituted a third cluster characterized by defects in insulin processing. A fourth cluster (TCF7L2, SLC30A8, HHEX/IDE, CDKAL1, CDKN2A/2B) was defined by loci influencing insulin processing and secretion without a detectable change in fasting glucose levels. The final group contained 20 risk loci with no clear-cut associations to continuous glycemic traits. By assembling extensive data on continuous glycemic traits, we have exposed the diverse mechanisms whereby type 2 diabetes risk variants impact disease predisposition.

Identification of low-frequency and rare sequence variants associated with elevated or reduced risk of type 2 diabetes.

Through whole-genome sequencing of 2,630 Icelanders and imputation into 11,114 Icelandic cases and 267,140 controls followed by testing in Danish and Iranian samples, we discovered 4 previously unreported variants affecting risk of type 2 diabetes (T2D). A low-frequency (1.47%) variant in intron 1 of CCND2, rs76895963[G], reduces risk of T2D by half (odds ratio (OR) = 0.53, P = 5.0 × 10−21) and is correlated with increased CCND2 expression. Notably, this variant is also associated with both greater height and higher body mass index (1.17 cm per allele, P = 5.5 × 10−12 and 0.56 kg/m2 per allele, P = 6.5 × 10−7, respectively). In addition, two missense variants in PAM, encoding p.Asp563Gly (frequency of 4.98%) and p.Ser539Trp (frequency of 0.65%), confer moderately higher risk of T2D (OR = 1.23, P = 3.9 × 10−10 and OR = 1.47, P = 1.7 × 10−5, respectively), and a rare (0.20%) frameshift variant in PDX1, encoding p.Gly218Alafs*12, associates with high risk of T2D (OR = 2.27, P = 7.3 × 10−7).

A common Greenlandic TBC1D4 variant confers muscle insulin resistance and type 2 diabetes

The Greenlandic population, a small and historically isolated founder population comprising about 57,000 inhabitants, has experienced a dramatic increase in type 2 diabetes (T2D) prevalence during the past 25 years. Motivated by this, we performed association mapping of T2D-related quantitative traits in up to 2,575 Greenlandic individuals without known diabetes. Using array-based genotyping and exome sequencing, we discovered a nonsense p.Arg684Ter variant (in which arginine is replaced by a termination codon) in the gene TBC1D4 with an allele frequency of 17%. Here we show that homozygous carriers of this variant have markedly higher concentrations of plasma glucose (β = 3.8 mmol l−1, P = 2.5 × 10−35) and serum insulin (β = 165 pmol l−1, P = 1.5 × 10−20) 2 hours after an oral glucose load compared with individuals with other genotypes (both non-carriers and heterozygous carriers). Furthermore, homozygous carriers have marginally lower concentrations of fasting plasma glucose (β = −0.18 mmol l−1, P = 1.1 × 10−6) and fasting serum insulin (β = −8.3 pmol l−1, P = 0.0014), and their T2D risk is markedly increased (odds ratio (OR) = 10.3, P = 1.6 × 10−24). Heterozygous carriers have a moderately higher plasma glucose concentration 2 hours after an oral glucose load than non-carriers (β = 0.43 mmol l−1, P = 5.3 × 10−5). Analyses of skeletal muscle biopsies showed lower messenger RNA and protein levels of the long isoform of TBC1D4, and lower muscle protein levels of the glucose transporter GLUT4, with increasing number of p.Arg684Ter alleles. These findings are concomitant with a severely decreased insulin-stimulated glucose uptake in muscle, leading to postprandial hyperglycaemia, impaired glucose tolerance and T2D. The observed effect sizes are several times larger than any previous findings in large-scale genome-wide association studies of these traits and constitute further proof of the value of conducting genetic association studies outside the traditional setting of large homogeneous populations.

Common type 2 diabetes risk gene variants associate with gestational diabetes.

OBJECTIVE We aimed to examine the association between gestational diabetes mellitus (GDM) and 11 recently identified type 2 diabetes susceptibility loci. RESEARCH DESIGN AND METHODS Type 2 diabetes risk variants in TCF7L2, CDKAL1, SLC30A8, HHEX/IDE, CDKN2A/2B, IGF2BP2, FTO, TCF2, PPARG, KCNJ11, and WFS1 loci were genotyped in a cohort of women with a history of GDM (n = 283) and glucose-tolerant women of the population-based Inter99 cohort (n = 2446). RESULTS All the risk alleles in the 11 examined type 2 diabetes risk variants showed an odds ratio (OR) greater than 1 for the GDM group compared with the control group ranging from 1.13 [95% confidence interval (CI) 0.88-1.46] to 1.44 (95% CI 1.19-1.74) except for the WFS1 rs10010131 variant with OR 0.87 (95% CI 0.73-1.05). Combined analysis of all 11 variants showed a highly significant additive effect of multiple risk alleles on risk of GDM [OR 1.18 (95% CI 1.10-1.27)] per risk allele, P = 3.2 x 10(-6)). Applying receiver-operating characteristic showed an area under the receiver-operating characteristic curve of 0.62 for the genetic test alone and 0.73 when combining information on age, body mass index, and genotypes of the 11 gene variants. CONCLUSIONS The prevalence in a prior GDM group of several previously proven type 2 diabetes risk alleles equals the findings from association studies on type 2 diabetes. This supports the hypothesis that GDM and type 2 diabetes are two of the same entity.