👤 Niels Grarup

Insulin sensitivity (IS) is a key determinant of metabolic health and may share genetic factors with obesity-related traits. Previous large-scale genetic studies have identified variants associated with IS as well as obesity related traits like body mass index (BMI) and waist-to-hip ratio (WHR). Notably, many of these associations are shared across traits, indicating a potential genetic overlap. However, the genetic intersection between IS and obesity-related traits remains underexplored. To explore this gap, we investigated associations between six IS indices, including fasting and post-glucose load measures, and genetic variants linked to BMI and WHR to determine their influence on IS and related cardiometabolic traits. To achieve this, we calculated six IS indices using fasting and oral glucose tolerance test (OGTT) data from 5,007 non-diabetic individuals, grouping them into fasting, OGTT Show less

Glucose-dependent insulinotropic polypeptide (GIP) and its receptor (GIPR) are involved in multiple physiological systems related to glucose metabolism, bone homeostasis and fat deposition. Recent research has surprisingly indicated that both agonists and antagonists of GIPR may be useful in the treatment of obesity and type 2 diabetes, as both result in weight loss when combined with GLP-1 receptor activation. To understand the receptor signaling related with weight loss, we examined the pharmacological properties of two rare missense Show less

Pancreatic beta-cell glucose sensitivity is the slope of the plasma glucose-insulin secretion relationship and is a key predictor of deteriorating glucose tolerance and development of type 2 diabetes. However, there are no large-scale studies looking at the genetic determinants of beta-cell glucose sensitivity. To understand the genetic determinants of pancreatic beta-cell glucose sensitivity using genome-wide meta-analysis and candidate gene studies. We performed a genome-wide meta-analysis for beta-cell glucose sensitivity in subjects with type 2 diabetes and nondiabetic subjects from 6 independent cohorts (n = 5706). Beta-cell glucose sensitivity was calculated from mixed meal and oral glucose tolerance tests, and its associations between known glycemia-related single nucleotide polymorphisms (SNPs) and genome-wide association study (GWAS) SNPs were estimated using linear regression models. Beta-cell glucose sensitivity was moderately heritable (h2 ranged from 34% to 55%) using SNP and family-based analyses. GWAS meta-analysis identified multiple correlated SNPs in the CDKAL1 gene and GIPR-QPCTL gene loci that reached genome-wide significance, with SNP rs2238691 in GIPR-QPCTL (P value = 2.64 × 10-9) and rs9368219 in the CDKAL1 (P value = 3.15 × 10-9) showing the strongest association with beta-cell glucose sensitivity. These loci surpassed genome-wide significance when the GWAS meta-analysis was repeated after exclusion of the diabetic subjects. After correction for multiple testing, glycemia-associated SNPs in or near the HHEX and IGF2B2 loci were also associated with beta-cell glucose sensitivity. We show that, variation at the GIPR-QPCTL and CDKAL1 loci are key determinants of pancreatic beta-cell glucose sensitivity. Show less

Although hundreds of genome-wide association studies-implicated loci have been reported for adult obesity-related traits, less is known about the genetics specific for early-onset obesity and with only a few studies conducted in non-European populations to date. Searching for additional genetic variants associated with childhood obesity, we performed a trans-ancestral meta-analysis of 30 studies consisting of up to 13 005 cases (≥95th percentile of body mass index (BMI) achieved 2-18 years old) and 15 599 controls (consistently <50th percentile of BMI) of European, African, North/South American and East Asian ancestry. Suggestive loci were taken forward for replication in a sample of 1888 cases and 4689 controls from seven cohorts of European and North/South American ancestry. In addition to observing 18 previously implicated BMI or obesity loci, for both early and late onset, we uncovered one completely novel locus in this trans-ancestral analysis (nearest gene, METTL15). The variant was nominally associated with only the European subgroup analysis but had a consistent direction of effect in other ethnicities. We then utilized trans-ancestral Bayesian analysis to narrow down the location of the probable causal variant at each genome-wide significant signal. Of all the fine-mapped loci, we were able to narrow down the causative variant at four known loci to fewer than 10 single nucleotide polymorphisms (SNPs) (FAIM2, GNPDA2, MC4R and SEC16B loci). In conclusion, an ethnically diverse setting has enabled us to both identify an additional pediatric obesity locus and further fine-map existing loci. Show less

Consanguine families display a high degree of homozygosity which increases the risk of family members suffering from autosomal recessive disorders. Thus, homozygous mutations in monogenic obesity genes may be a more frequent cause of childhood obesity in a consanguineous population. We identified 23 probands from 23 Pakistani families displaying autosomal recessive obesity. We have previously excluded mutations in MC4R, LEP and LEPR in all probands. Using a chip-based, target-region capture array, 31 genes involved in monogenic forms of obesity, were screened in all probands. We identified 31 rare non-synonymous possibly pathogenic variants (28 missense and three nonsense) within the 31 selected genes. All variants were heterozygous, thus no homozygous pathogenic variants were found. Two of the rare heterozygous nonsense variants identified (p.R75X and p.R481X) were found in BBS9 within one proband, suggesting that obesity is caused by compound heterozygosity. Sequencing of the parents supported the compound heterozygous nature of obesity as each parent was carrying one of the variants. Subsequent clinical investigation strongly indicated that the proband had Bardet-Biedl syndrome. Mutation screening in 31 genes among probands with severe early-onset obesity from Pakistani families did not reveal the presence of homozygous obesity causing variants. However, a compound heterozygote carrier of BBS9 mutations was identified, indicating that compound heterozygosity must not be overlooked when investigating the genetic etiology of severe childhood obesity in populations with a high degree of consanguinity. Show less

Genome-wide association studies (GWAS) have identified >250 loci for body mass index (BMI), implicating pathways related to neuronal biology. Most GWAS loci represent clusters of common, noncoding variants from which pinpointing causal genes remains challenging. Here we combined data from 718,734 individuals to discover rare and low-frequency (minor allele frequency (MAF) < 5%) coding variants associated with BMI. We identified 14 coding variants in 13 genes, of which 8 variants were in genes (ZBTB7B, ACHE, RAPGEF3, RAB21, ZFHX3, ENTPD6, ZFR2 and ZNF169) newly implicated in human obesity, 2 variants were in genes (MC4R and KSR2) previously observed to be mutated in extreme obesity and 2 variants were in GIPR. The effect sizes of rare variants are ~10 times larger than those of common variants, with the largest effect observed in carriers of an MC4R mutation introducing a stop codon (p.Tyr35Ter, MAF = 0.01%), who weighed ~7 kg more than non-carriers. Pathway analyses based on the variants associated with BMI confirm enrichment of neuronal genes and provide new evidence for adipocyte and energy expenditure biology, widening the potential of genetically supported therapeutic targets in obesity. Show less

Mutations in the genes encoding leptin (LEP), the leptin receptor (LEPR), and the melanocortin 4 receptor (MC4R) are known to cause severe early-onset childhood obesity. The aim of the current study was to examine the prevalence of damaging LEP, LEPR, and MC4R mutations in Pakistani families having a recessive heritance of early-onset obesity. Using targeted resequencing, the presence of rare mutations in LEP, LEPR, and MC4R, was investigated in individuals from 25 families suspected of having autosomal recessive early-onset obesity. Segregation patterns of variants were assessed based on chip-based genotyping. Homozygous LEPR variants were identified in two probands. One carried a deletion (c.3260AG) resulting in the frameshift mutation p.Ser1090Trpfs*6, and the second carried a substitution (c.2675C > G) resulting in the missense mutation p.Pro892Arg. Both mutations were located within regions of homozygosity shared only among affected individuals. Both probands displayed early-onset obesity, hyperphagia and diabetes. No mutations were found in LEP and MC4R. The current study highlights the implication of LEPR mutations in cases of severe early-onset obesity in consanguineous Pakistani families. Through targeted resequencing, we identified novel damaging mutations, and our approach may therefore be utilized in clinical testing or diagnosis of known forms of monogenic obesity with the aim of optimizing obesity treatment. Show less

We have identified a variant in ADCY3 (encoding adenylate cyclase 3) associated with markedly increased risk of obesity and type 2 diabetes in the Greenlandic population. The variant disrupts a splice acceptor site, and carriers have decreased ADCY3 RNA expression. Additionally, we observe an enrichment of rare ADCY3 loss-of-function variants among individuals with type 2 diabetes in trans-ancestry cohorts. These findings provide new information on disease etiology relevant for future treatment strategies. Show less

We screened variants on an exome-focused genotyping array in >300,000 participants (replication in >280,000 participants) and identified 444 independent variants in 250 loci significantly associated with total cholesterol (TC), high-density-lipoprotein cholesterol (HDL-C), low-density-lipoprotein cholesterol (LDL-C), and/or triglycerides (TG). At two loci (JAK2 and A1CF), experimental analysis in mice showed lipid changes consistent with the human data. We also found that: (i) beta-thalassemia trait carriers displayed lower TC and were protected from coronary artery disease (CAD); (ii) excluding the CETP locus, there was not a predictable relationship between plasma HDL-C and risk for age-related macular degeneration; (iii) only some mechanisms of lowering LDL-C appeared to increase risk for type 2 diabetes (T2D); and (iv) TG-lowering alleles involved in hepatic production of TG-rich lipoproteins (TM6SF2 and PNPLA3) tracked with higher liver fat, higher risk for T2D, and lower risk for CAD, whereas TG-lowering alleles involved in peripheral lipolysis (LPL and ANGPTL4) had no effect on liver fat but decreased risks for both T2D and CAD. Show less

Hand grip strength is a widely used proxy of muscular fitness, a marker of frailty, and predictor of a range of morbidities and all-cause mortality. To investigate the genetic determinants of variation in grip strength, we perform a large-scale genetic discovery analysis in a combined sample of 195,180 individuals and identify 16 loci associated with grip strength (P<5 × 10 Show less

Understanding the physiological mechanisms by which common variants predispose to type 2 diabetes requires large studies with detailed measures of insulin secretion and sensitivity. Here we performed the largest genome-wide association study of first-phase insulin secretion, as measured by intravenous glucose tolerance tests, using up to 5,567 individuals without diabetes from 10 studies. We aimed to refine the mechanisms of 178 known associations between common variants and glycemic traits and identify new loci. Thirty type 2 diabetes or fasting glucose-raising alleles were associated with a measure of first-phase insulin secretion at Show less

A large number of genetic loci are associated with adult body mass index. However, the genetics of childhood body mass index are largely unknown. We performed a meta-analysis of genome-wide association studies of childhood body mass index, using sex- and age-adjusted standard deviation scores. We included 35 668 children from 20 studies in the discovery phase and 11 873 children from 13 studies in the replication phase. In total, 15 loci reached genome-wide significance (P-value < 5 × 10(-8)) in the joint discovery and replication analysis, of which 12 are previously identified loci in or close to ADCY3, GNPDA2, TMEM18, SEC16B, FAIM2, FTO, TFAP2B, TNNI3K, MC4R, GPR61, LMX1B and OLFM4 associated with adult body mass index or childhood obesity. We identified three novel loci: rs13253111 near ELP3, rs8092503 near RAB27B and rs13387838 near ADAM23. Per additional risk allele, body mass index increased 0.04 Standard Deviation Score (SDS) [Standard Error (SE) 0.007], 0.05 SDS (SE 0.008) and 0.14 SDS (SE 0.025), for rs13253111, rs8092503 and rs13387838, respectively. A genetic risk score combining all 15 SNPs showed that each additional average risk allele was associated with a 0.073 SDS (SE 0.011, P-value = 3.12 × 10(-10)) increase in childhood body mass index in a population of 1955 children. This risk score explained 2% of the variance in childhood body mass index. This study highlights the shared genetic background between childhood and adult body mass index and adds three novel loci. These loci likely represent age-related differences in strength of the associations with body mass index. Show less

Human complex metabolic traits are in part regulated by genetic determinants. Here we applied exome sequencing to identify novel associations of coding polymorphisms at minor allele frequencies (MAFs) >1% with common metabolic phenotypes. The study comprised three stages. We performed medium-depth (8×) whole exome sequencing in 1,000 cases with type 2 diabetes, BMI >27.5 kg/m(2) and hypertension and in 1,000 controls (stage 1). We selected 16,192 polymorphisms nominally associated (p < 0.05) with case-control status, from four selected annotation categories or from loci reported to associate with metabolic traits. These variants were genotyped in 15,989 Danes to search for association with 12 metabolic phenotypes (stage 2). In stage 3, polymorphisms showing potential associations were genotyped in a further 63,896 Europeans. Exome sequencing identified 70,182 polymorphisms with MAF >1%. In stage 2 we identified 51 potential associations with one or more of eight metabolic phenotypes covered by 45 unique polymorphisms. In meta-analyses of stage 2 and stage 3 results, we demonstrated robust associations for coding polymorphisms in CD300LG (fasting HDL-cholesterol: MAF 3.5%, p = 8.5 × 10(-14)), COBLL1 (type 2 diabetes: MAF 12.5%, OR 0.88, p = 1.2 × 10(-11)) and MACF1 (type 2 diabetes: MAF 23.4%, OR 1.10, p = 8.2 × 10(-10)). We applied exome sequencing as a basis for finding genetic determinants of metabolic traits and show the existence of low-frequency and common coding polymorphisms with impact on common metabolic traits. Based on our study, coding polymorphisms with MAF above 1% do not seem to have particularly high effect sizes on the measured metabolic traits. Show less

Numerous obesity loci have been identified using genome-wide association studies. A UK study indicated that physical activity may attenuate the cumulative effect of 12 of these loci, but replication studies are lacking. Therefore, we tested whether the aggregate effect of these loci is diminished in adults of European ancestry reporting high levels of physical activity. Twelve obesity-susceptibility loci were genotyped or imputed in 111,421 participants. A genetic risk score (GRS) was calculated by summing the BMI-associated alleles of each genetic variant. Physical activity was assessed using self-administered questionnaires. Multiplicative interactions between the GRS and physical activity on BMI were tested in linear and logistic regression models in each cohort, with adjustment for age, age(2), sex, study center (for multicenter studies), and the marginal terms for physical activity and the GRS. These results were combined using meta-analysis weighted by cohort sample size. The meta-analysis yielded a statistically significant GRS × physical activity interaction effect estimate (Pinteraction  = 0.015). However, a statistically significant interaction effect was only apparent in North American cohorts (n = 39,810, Pinteraction  = 0.014 vs. n = 71,611, Pinteraction  = 0.275 for Europeans). In secondary analyses, both the FTO rs1121980 (Pinteraction  = 0.003) and the SEC16B rs10913469 (Pinteraction  = 0.025) variants showed evidence of SNP × physical activity interactions. This meta-analysis of 111,421 individuals provides further support for an interaction between physical activity and a GRS in obesity disposition, although these findings hinge on the inclusion of cohorts from North America, indicating that these results are either population-specific or non-causal. Show less

Thirty-two common variants associated with body mass index (BMI) have been identified in genome-wide association studies, explaining ∼1.45% of BMI variation in general population cohorts. We performed a genome-wide association study in a sample of young adults enriched for extremely overweight individuals. We aimed to identify new loci associated with BMI and to ascertain whether using an extreme sampling design would identify the variants known to be associated with BMI in general populations. From two large Danish cohorts we selected all extremely overweight young men and women (n = 2,633), and equal numbers of population-based controls (n = 2,740, drawn randomly from the same populations as the extremes, representing ∼212,000 individuals). We followed up novel (at the time of the study) association signals (p<0.001) from the discovery cohort in a genome-wide study of 5,846 Europeans, before attempting to replicate the most strongly associated 28 SNPs in an independent sample of Danish individuals (n = 20,917) and a population-based cohort of 15-year-old British adolescents (n = 2,418). Our discovery analysis identified SNPs at three loci known to be associated with BMI with genome-wide confidence (P<5×10(-8); FTO, MC4R and FAIM2). We also found strong evidence of association at the known TMEM18, GNPDA2, SEC16B, TFAP2B, SH2B1 and KCTD15 loci (p<0.001), and nominal association (p<0.05) at a further 8 loci known to be associated with BMI. However, meta-analyses of our discovery and replication cohorts identified no novel associations. Our results indicate that the detectable genetic variation associated with extreme overweight is very similar to that previously found for general BMI. This suggests that population-based study designs with enriched sampling of individuals with the extreme phenotype may be an efficient method for identifying common variants that influence quantitative traits and a valid alternative to genotyping all individuals in large population-based studies, which may require tens of thousands of subjects to achieve similar power. Show less

A genome-wide association study in the Japanese population reported two genome-wide significant loci associated with type 2 diabetes of which the VPS13C/C2CD4A/C2CD4B locus was replicated in Europeans. We looked for potential associations between the diabetogenic VPS13C/C2CD4A/C2CD4B rs7172432 variant and diabetes-related intermediary traits. We genotyped the rs7172432 variant in the population-based Inter99 cohort (n = 6,784) and analysed quantitative diabetes-related traits in 5,722 non-diabetic participants who all were examined by an OGTT. The diabetes-associated A allele was associated with 0.60 cm higher waist circumference (p = 0.004), 0.037 mmol/l higher fasting plasma glucose (p = 4 × 10(-5)) and 0.11 mmol/l higher plasma glucose at 30 min during an OGTT (p = 4 × 10(-4)). In analyses adjusted for concomitant insulin sensitivity levels the diabetogenic allele was associated with a lower acute glucose-stimulated insulin response (GSIR) as estimated by 30 min serum insulin (β = -0.039, p = 2 × 10(-7)), insulinogenic index (β = -0.057, p = 1 × 10(-8)) and BIGTT-acute insulin release (β = -0.041, p = 9 × 10(-9)). As rs7172432 is situated in a region previously associated with glycaemic traits, we tested linkage disequilibrium (LD) with the reported regional lead single-nucleotide polymorphisms for fasting (rs11071657) and 2 h plasma glucose (rs17271305), and performed conditional analyses of rs7172432. Rs7172432 showed moderate LD with rs11071657 and rs17271305 (R (2) < 0.34) and we found strong association by almost unchanged effect sizes of rs7172432 with plasma glucose and estimates of GSIR in analyses conditional on rs11071657 and rs17271305. The diabetogenic VPS13C/C2CD4A/C2CD4B rs7172432 A allele associates with GSIR in non-diabetic individuals from the general population, suggesting an impaired beta cell function as an intermediary diabetes-related trait. Show less

Glucose levels 2 h after an oral glucose challenge are a clinical measure of glucose tolerance used in the diagnosis of type 2 diabetes. We report a meta-analysis of nine genome-wide association studies (n = 15,234 nondiabetic individuals) and a follow-up of 29 independent loci (n = 6,958-30,620). We identify variants at the GIPR locus associated with 2-h glucose level (rs10423928, beta (s.e.m.) = 0.09 (0.01) mmol/l per A allele, P = 2.0 x 10(-15)). The GIPR A-allele carriers also showed decreased insulin secretion (n = 22,492; insulinogenic index, P = 1.0 x 10(-17); ratio of insulin to glucose area under the curve, P = 1.3 x 10(-16)) and diminished incretin effect (n = 804; P = 4.3 x 10(-4)). We also identified variants at ADCY5 (rs2877716, P = 4.2 x 10(-16)), VPS13C (rs17271305, P = 4.1 x 10(-8)), GCKR (rs1260326, P = 7.1 x 10(-11)) and TCF7L2 (rs7903146, P = 4.2 x 10(-10)) associated with 2-h glucose. Of the three newly implicated loci (GIPR, ADCY5 and VPS13C), only ADCY5 was found to be associated with type 2 diabetes in collaborating studies (n = 35,869 cases, 89,798 controls, OR = 1.12, 95% CI 1.09-1.15, P = 4.8 x 10(-18)). Show less

A meta-analysis of 21 genome-wide association studies identified 11 novel genetic loci implicated in fasting glucose homeostasis. We aimed to evaluate the impact of these variants on insulin release and insulin sensitivity estimated from OGTTs. Eleven variants in or near DGKB/TMEM195, ADCY5, MADD, ADRA2A, FADS1, CRY2, SLC2A2, GLIS3, PROX1, C2CD4B and IGF1 were genotyped in 6,784 middle-aged participants of the population-based Inter99 cohort. Association studies of quantitative estimates of insulin release and insulin sensitivity were performed in 5,722 non-diabetic Danish participants on whom an OGTT was performed. Assuming an additive genetic model, carriers of the alleles increasing fasting glucose in DGKB/TMEM195, ADRA2A, GLIS3 and C2CD4B showed decreased glucose-stimulated insulin release as assessed by the BIGTT-acute insulin response index (2.7-3.5%; p < 0.005 for all) and by corrected insulin response (2.8-5.9%; p < 0.03 for all). In addition, the PROX1 glucose-raising allele showed a 2.9% decreased corrected insulin response (p = 0.03), while the hyperglycaemic allele of variants in or near ADRA2A, FADS1, CRY2 and C2CD4B were associated with a 2.6% to 9.3% decrease in one or both of two different OGTT-based disposition indices (p < 0.02 for all). After correction for multiple testing, variants in the DGKB/TMEM195, ADRA2A, GLIS3 and C2CD4B loci were associated with estimates of beta cell function. We found that the lead variants at the DGKB/TMEM195, ADRA2A, GLIS3 and C2CD4B loci were associated with decreased glucose-stimulated insulin response. This association underlines the importance of pancreatic beta cell dysfunction in the genetic predisposition to hyperglycaemia and type 2 diabetes. Show less

Levels of circulating glucose are tightly regulated. To identify new loci influencing glycemic traits, we performed meta-analyses of 21 genome-wide association studies informative for fasting glucose, fasting insulin and indices of beta-cell function (HOMA-B) and insulin resistance (HOMA-IR) in up to 46,186 nondiabetic participants. Follow-up of 25 loci in up to 76,558 additional subjects identified 16 loci associated with fasting glucose and HOMA-B and two loci associated with fasting insulin and HOMA-IR. These include nine loci newly associated with fasting glucose (in or near ADCY5, MADD, ADRA2A, CRY2, FADS1, GLIS3, SLC2A2, PROX1 and C2CD4B) and one influencing fasting insulin and HOMA-IR (near IGF1). We also demonstrated association of ADCY5, PROX1, GCK, GCKR and DGKB-TMEM195 with type 2 diabetes. Within these loci, likely biological candidate genes influence signal transduction, cell proliferation, development, glucose-sensing and circadian regulation. Our results demonstrate that genetic studies of glycemic traits can identify type 2 diabetes risk loci, as well as loci containing gene variants that are associated with a modest elevation in glucose levels but are not associated with overt diabetes. Show less

Gene-environment interaction can be viewed as a departure from an otherwise expected additivity of genetic and environmental factors on a given outcome measure. Important genetic and environmental factors contribute to the pathogenesis of type 2 diabetes and intermediary traits, probably modulated by their complex interaction. This paper provides an update on the current literature investigating gene-environment interactions of type 2 diabetes and metabolic phenotypes, and discusses the future perspectives of this research. Recent advances in gene-environment interaction studies of metabolism have involved LIPC, APOA5 and PPARG variation, and nutrition and physical activity, of which the most consistently replicated observations have been obtained for APOA5. Also, intervention studies of the promising TCF7L2 type 2 diabetes gene and possible future strategies are discussed. Possibly as a result of the complexity of these multifactorial diseases, recent years have seen only limited success in unravelling significant gene-environment interactions, but important insights have been gained and they hold promise for implementation in lifestyle intervention strategies. We need to evolve to more complex, but realistic, scenarios involving several genes and environmental factors. Recent progress in statistical methods allowing for higher-order interactions may make this possible. Show less