Publications

In the Greenlandic Inuit population, 4% are homozygous carriers of a genetic nonsense TBC1D4 p.Arg684Ter variant leading to loss of the muscle-specific isoform of TBC1D4 and an approximately tenfold increased risk of type 2 diabetes1. Here we show the metabolic consequences of this variant in four female and four male homozygous carriers and matched controls. An extended glucose tolerance test reveals prolonged hyperglycaemia followed by reactive hypoglycaemia in the carriers. Whole-body glucose disposal is impaired during euglycaemic-hyperinsulinaemic clamp conditions and associates with severe insulin resistance in skeletal muscle only. Notably, a marked reduction in muscle glucose transporter GLUT4 and associated proteins is observed. While metabolic regulation during exercise remains normal, the insulin-sensitizing effect of a single exercise bout is compromised. Thus, loss of the muscle-specific isoform of TBC1D4 causes severe skeletal muscle insulin resistance without baseline hyperinsulinaemia. However, physical activity can ameliorate this condition. These observations offer avenues for personalized interventions and targeted preventive strategies.

BACKGROUND: Preserved ratio impaired spirometry (PRISm) is defined as a forced expiratory volume in 1 s (FEV1) <80% predicted and FEV1/forced vital capacity ≥0.70. PRISm is associated with respiratory symptoms and comorbidities. Our objective was to discover novel genetic signals for PRISm and see if they provide insight into the pathogenesis of PRISm and associated comorbidities.

METHODS: We undertook a genome-wide association study (GWAS) of PRISm in UK Biobank participants (Stage 1), and selected single nucleotide polymorphisms (SNPs) reaching genome-wide significance for replication in 13 cohorts (Stage 2). A combined meta-analysis of Stage 1 and Stage 2 was done to determine top SNPs. We used cross-trait linkage disequilibrium score regression to estimate genome-wide genetic correlation between PRISm and pulmonary and extrapulmonary traits. Phenome-wide association studies of top SNPs were performed.

RESULTS: 22 signals reached significance in the joint meta-analysis, including four signals novel for lung function. A strong genome-wide genetic correlation (rg) between PRISm and spirometric COPD (rg=0.62, p<0.001) was observed, and genetic correlation with type 2 diabetes (rg=0.12, p=0.007). Phenome-wide association studies showed that 18 of 22 signals were associated with diabetic traits and seven with blood pressure traits.

CONCLUSION: This is the first GWAS to successfully identify SNPs associated with PRISm. Four of the signals, rs7652391 (nearest gene MECOM), rs9431040 (HLX), rs62018863 (TMEM114) and rs185937162 (HLA-B), have not been described in association with lung function before, demonstrating the utility of using different lung function phenotypes in GWAS. Genetic factors associated with PRISm are strongly correlated with risk of both other lung diseases and extrapulmonary comorbidity.

Incretin-based therapies are highly successful in combatting obesity and type 2 diabetes1. Yet both activation and inhibition of the glucose-dependent insulinotropic polypeptide (GIP) receptor (GIPR) in combination with glucagon-like peptide-1 (GLP-1) receptor (GLP-1R) activation have resulted in similar clinical outcomes, as demonstrated by the GIPR-GLP-1R co-agonist tirzepatide2 and AMG-133 (ref. 3) combining GIPR antagonism with GLP-1R agonism. This underlines the importance of a better understanding of the GIP system. Here we show the necessity of β-arrestin recruitment for GIPR function, by combining in vitro pharmacological characterization of 47 GIPR variants with burden testing of clinical phenotypes and in vivo studies. Burden testing of variants with distinct ligand-binding capacity, Gs activation (cyclic adenosine monophosphate production) and β-arrestin 2 recruitment and internalization shows that unlike variants solely impaired in Gs signalling, variants impaired in both Gs and β-arrestin 2 recruitment contribute to lower adiposity-related traits. Endosomal Gs-mediated signalling of the variants shows a β-arrestin dependency and genetic ablation of β-arrestin 2 impairs cyclic adenosine monophosphate production and decreases GIP efficacy on glucose control in male mice. This study highlights a crucial impact of β-arrestins in regulating GIPR signalling and overall preservation of biological activity that may facilitate new developments in therapeutic targeting of the GIPR system.

The microbial communities of the oral cavity are important elements of oral and systemic health. With emerging evidence highlighting the heritability of oral bacterial microbiota, this study aimed to identify host genome variants that influence oral microbial traits. Using data from 16S rRNA gene amplicon sequencing, we performed genome-wide association studies with univariate and multivariate traits of the salivary microbiota from 610 unrelated adults from the Danish ADDITION-PRO cohort. We identified six single nucleotide polymorphisms (SNPs) in human genomes that showed associations with abundance of bacterial taxa at different taxonomical tiers (P < 5 × 10-8). Notably, SNP rs17793860 surpassed our study-wide significance threshold (P < 1.19 × 10-9). Additionally, rs4530093 was linked to bacterial beta diversity (P < 5 × 10-8). Out of these seven SNPs identified, six exerted effects on metabolic traits, including glycated hemoglobin A1c, triglyceride and high-density lipoprotein cholesterol levels, the risk of type 2 diabetes and stroke. Our findings highlight the impact of specific host SNPs on the composition and diversity of the oral bacterial community. Importantly, our results indicate an intricate interplay between host genetics, the oral microbiota, and metabolic health. We emphasize the need for integrative approaches considering genetic, microbial, and metabolic factors.

BACKGROUND: Birth weight (BW) is associated with risk of cardiometabolic disease (CMD) in adulthood, which may depend on the state of obesity, in particular if developed at a young age. We hypothesised that BW and a polygenic score (PGS) for BW were associated with cardiometabolic risk and related plasma protein levels in children and adolescents. We aimed to determine the modifying effect of childhood obesity on these associations.

METHODS: We used data from The cross-sectional HOLBAEK Study with 4263 participants (median [IQR] age, 11.7 [9.2, 14.3] years; 57.1% girls and 42.9% boys; 48.6% from an obesity clinic and 51.4% from a population-based group). We gathered information on BW and gestational age, anthropometrics, cardiometabolic risk factors, calculated a PGS for BW, and measured plasma proteins using Olink Inflammation and Cardiovascular II panels. We employed multiple linear regression to examine the associations with BW as a continuous variable and performed interaction analyses to assess the effect of childhood obesity on cardiometabolic risk and plasma protein levels.

FINDINGS: BW and a PGS for BW associated with cardiometabolic risk and plasma protein levels in childhood and adolescence. Childhood obesity modified the associations between BW and measures of insulin resistance, including HOMA-IR (βadj [95% CI per SD] for obesity: -0.12 [-0.15, -0.08]; normal weight: -0.04 [-0.08, 0.00]; Pinteraction = 0.004), c-peptide (obesity: -0.11 [-0.14, -0.08]; normal weight: -0.02 [-0.06, 0.02]; Pinteraction = 5.05E-04), and SBP SDS (obesity: -0.12 [-0.16, -0.08]; normal weight: -0.06 [-0.11, -0.01]; Pinteraction = 0.0479). Childhood obesity also modified the associations between BW and plasma levels of 14 proteins (e.g., IL15RA, MCP1, and XCL1; Pinteraction < 0.05).

INTERPRETATION: We identified associations between lower BW and adverse metabolic phenotypes, particularly insulin resistance, blood pressure, and altered plasma protein levels, which were more pronounced in children with obesity. Developing effective prevention and treatment strategies for this group is needed to reduce the risk of future CMD.

FUNDING: Novo Nordisk Foundation (NNF15OC0016544, NNF0064142 to T.H., NNF15OC0016692 to T.H. and A.K., NNF18CC0033668 to S.E.S, NNF18SA0034956 to C.E.F., NNF20SA0067242 to DCA, NNF18CC0034900 to NNF CBMR), The Innovation Fund Denmark (0603-00484B to T.H.), The Danish Cardiovascular Academy (DCA) and the Danish Heart Foundation (HF) (PhD2021007-DCA to P.K.R, 18-R125-A8447-22088 (HF) and 21-R149-A10071-22193 (HF) to M.A.V.L., PhD2023009-HF to L.A.H), EU Horizon (668031, 847989, 825694, 964590 to A.K.), Innovative Health Initiative (101132901 for A.K.), A.P. Møller Foundation (19-L-0366 to T.H.), The Danish National Research Foundation, Steno Diabetes Center Sjælland, and The Region Zealand and Southern Denmark Health Scientific Research Foundation.

BACKGROUND: Disease prevalence and mean phenotype values differ between many populations, including Inuit and Europeans. Whether these differences are partly explained by genetic differences or solely due to differences in environmental exposures is still unknown, because estimates of the genetic contribution to these means, which we will here refer to as mean genotypic values, are easily confounded, and because studies across genetically diverse populations are lacking.

METHODS: Leveraging the unique genetic properties of the small, admixed and historically isolated Greenlandic population, we estimated the differences in mean genotypic value between Inuit and European genetic ancestry using an admixed sibling design. Analyses were performed across 26 metabolic phenotypes, in 1474 admixed sibling pairs present in a cohort of 5996 Greenlanders.

RESULTS: After FDR correction for multiple testing, we found significantly lower mean genotypic values in Inuit genetic ancestry compared to European genetic ancestry for body weight (effect size per percentage of Inuit genetic ancestry (se), -0.51 (0.16) kg/%), body mass index (-0.20 (0.06) kg/m2/%), fat percentage (-0.38 (0.13) %/%), waist circumference (-0.42 (0.16) cm/%), hip circumference (-0.38 (0.11) cm/%) and fasting serum insulin levels (-1.07 (0.51) pmol/l/%). The direction of the effects was consistent with the observed mean phenotype differences between Inuit and European genetic ancestry. No difference in mean genotypic value was observed for height, markers of glucose homeostasis, or circulating lipid levels.

CONCLUSIONS: We show that mean genotypic values for some metabolic phenotypes differ between two human populations using a method not easily confounded by possible differences in environmental exposures. Our study illustrates the importance of performing genetic studies in diverse populations.

Perturbation of lipid homoeostasis is a major risk factor for cardiovascular disease (CVD), the leading cause of death worldwide. We aimed to identify genetic variants affecting lipid levels, and thereby risk of CVD, in Greenlanders. Genome-wide association studies (GWAS) of six blood lipids, triglycerides, LDL-cholesterol, HDL-cholesterol, total cholesterol, as well as apolipoproteins A1 and B, were performed in up to 4473 Greenlanders. For genome-wide significant variants, we also tested for associations with additional traits, including CVD events. We identified 11 genome-wide significant loci associated with lipid traits. Most of these loci were already known in Europeans, however, we found a potential causal variant near PCSK9 (rs12117661), which was independent of the known PCSK9 loss-of-function variant (rs11491147). rs12117661 was associated with lower LDL-cholesterol (βSD(SE) = -0.22 (0.03), p = 6.5 × 10-12) and total cholesterol (-0.17 (0.03), p = 1.1 × 10-8) in the Greenlandic study population. Similar associations were observed in Europeans from the UK Biobank, where the variant was also associated with a lower risk of CVD outcomes. Moreover, rs12117661 was a top eQTL for PCSK9 across tissues in European data from the GTEx portal, and was located in a predicted regulatory element, supporting a possible causal impact on PCSK9 expression. Combined, the 11 GWAS signals explained up to 16.3% of the variance of the lipid traits. This suggests that the genetic architecture of lipid levels in Greenlanders is different from Europeans, with fewer variants explaining the variance.

BACKGROUND: The emerging use of biomarkers in research and tailored care introduces a need for information about the association between biomarkers and basic demographics and lifestyle factors revealing expectable concentrations in healthy individuals while considering general demographic differences.

METHODS: A selection of 47 biomarkers, including markers of inflammation and vascular stress, were measured in plasma samples from 9876 Danish Blood Donor Study participants. Using regression models, we examined the association between biomarkers and sex, age, Body Mass Index (BMI), and smoking.

RESULTS: Here we show that concentrations of inflammation and vascular stress biomarkers generally increase with higher age, BMI, and smoking. Sex-specific effects are observed for multiple biomarkers.

CONCLUSION: This study provides comprehensive information on concentrations of 47 plasma biomarkers in healthy individuals. The study emphasizes that knowledge about biomarker concentrations in healthy individuals is critical for improved understanding of disease pathology and for tailored care and decision support tools.

BACKGROUND: Ventricular repolarization time (ECG QT and JT intervals) is associated with malignant arrhythmia. Genome-wide association studies have identified 230 independent loci for QT and JT; however, 50% of their heritability remains unexplained. Previous work supports a causal effect of lower serum calcium concentrations on longer ventricular repolarization time. We hypothesized calcium interactions with QT and JT variant associations could explain a proportion of the missing heritability.

METHODS AND RESULTS: We performed genome-wide calcium interaction analyses for QT and JT intervals. Participants were stratified by their calcium level relative to the study distribution (top or bottom 20%). We performed a 2-stage analysis (genome-wide discovery [N=62 532] and replication [N=59 861] of lead variants) and a single-stage genome-wide meta-analysis (N=122 393, [European ancestry N=117 581, African ancestry N=4812]). We also calculated 2-degrees of freedom joint main and interaction and 1-degree of freedom interaction P values. In 2-stage and single-stage analyses, 50 and 98 independent loci, respectively, were associated with either QT or JT intervals (2-degrees of freedom joint main and interaction P value <5×10-8). No lead variant had a significant interaction result after correcting for multiple testing and sensitivity analyses provided similar findings. Two loci in the single-stage meta-analysis were not reported previously (SPPL2B and RFX6).

CONCLUSIONS: We have found limited support for an interaction effect of serum calcium on QT and JT variant associations despite sample sizes with suitable power to detect relevant effects. Therefore, such effects are unlikely to explain a meaningful proportion of the heritability of QT and JT, and factors including rare variation and other environmental interactions need to be considered.

There are multiple independent genetic signals at the Ras-responsive element binding protein 1 (RREB1) locus associated with type 2 diabetes risk, fasting glucose, ectopic fat, height, and bone mineral density. We have previously shown that loss of RREB1 in pancreatic beta cells reduces insulin content and impairs islet cell development and function. However, RREB1 is a widely expressed transcription factor and the metabolic impact of RREB1 loss in vivo remains unknown. Here, we show that male and female global heterozygous knockout (Rreb1 +/-) mice have reduced body length, weight, and fat mass on high-fat diet. Rreb1+/- mice have sex- and diet-specific decreases in adipose tissue and adipocyte size; male mice on high-fat diet had larger gonadal adipocytes, while males on standard chow and females on high-fat diet had smaller, more insulin sensitive subcutaneous adipocytes. Mouse and human precursor cells lacking RREB1 have decreased adipogenic gene expression and activated transcription of genes associated with osteoblast differentiation, which was associated with Rreb1 +/- mice having increased bone mineral density in vivo. Finally, human carriers of RREB1 T2D protective alleles have smaller adipocytes, consistent with RREB1 loss-of-function reducing diabetes risk.