Publications

OBJECTIVES: Experimental and genetic studies suggest that fibroblast growth factor 21 (FGF21) modulates macronutrient and alcohol preferences, but evidence of such regulation in humans remains scarce. To address this gap in translation, we aimed to map the relationships between plasma FGF21 levels, FGF21 genetic variation and habitual macronutrient intake in a large human population.

METHODS: We fine-mapped and performed colocalization of the FGF21 genetic region in GWAS summary statistics of plasma FGF21 levels and macronutrient intake. UK Biobank data were used to investigate the associations between FGF21 genetic variants, plasma FGF21 protein levels, and macronutrient intake (including alcohol) assessed with repeated 24-hour recalls. One- and two-sample mendelian randomization were performed to estimate the effects of plasma FGF21 on macronutrient intake.

RESULTS: We show that the main macronutrient-associated variant rs838133 and the FGF21 cis-pQTL rs838131, both in the FGF21 gene, are distinct genetic signals. Effect directions also suggest that the influence of FGF21 variation on macronutrient intake appear more complex than by direct mediation through plasma FGF21. Only when considering this complexity at FGF21, is plasma FGF21 estimated to reduce alcohol and increase protein and fat intake using mendelian randomization. Importantly, plasma FGF21 levels also appear markedly elevated by primarily high alcohol and low protein intake.

CONCLUSIONS: These findings support the feedback diet-regulatory mechanism of FGF21 in humans, but highlights the need for mechanistic characterization of the complex FGF21 genetic region.

To date only a fraction of the genetic footprint of thyroid function has been clarified. We report a genome-wide association study meta-analysis of thyroid function in up to 271,040 individuals of European ancestry, including reference range thyrotropin (TSH), free thyroxine (FT4), free and total triiodothyronine (T3), proxies for metabolism (T3/FT4 ratio) as well as dichotomized high and low TSH levels. We revealed 259 independent significant associations for TSH (61% novel), 85 for FT4 (67% novel), and 62 novel signals for the T3 related traits. The loci explained 14.1%, 6.0%, 9.5% and 1.1% of the total variation in TSH, FT4, total T3 and free T3 concentrations, respectively. Genetic correlations indicate that TSH associated loci reflect the thyroid function determined by free T3, whereas the FT4 associations represent the thyroid hormone metabolism. Polygenic risk score and Mendelian randomization analyses showed the effects of genetically determined variation in thyroid function on various clinical outcomes, including cardiovascular risk factors and diseases, autoimmune diseases, and cancer. In conclusion, our results improve the understanding of thyroid hormone physiology and highlight the pleiotropic effects of thyroid function on various diseases.

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: 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.

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.

BACKGROUND: The relationship between combined genetic predisposition and lifestyle and the risk of incident atrial fibrillation (AF) is unclear. Therefore, we aimed to assess a possible interaction between lifestyle and genetics on AF risk.

METHODS: We included AF cases and a randomly drawn subcohort of 4040 participants from the Danish Diet, Cancer and Health cohort. Lifestyle risk factors were assessed, a score was calculated, and participants were categorised as having a poor, intermediate, or ideal lifestyle. We calculated a genetic risk score comprising 142 variants, and categorised participants into low (quintile 1), intermediate (quintiles 2-4) or high (quintile 5) genetic risk of AF.

RESULTS: 3094 AF cases occurred during a median follow-up of 12.9 years. Regardless of genetic risk, incidence rates per 1000 person-years were gradually higher with worse lifestyle. For participants with high genetic risk, the incidence rates of AF per 1000 person-years were 5.0 (95% CI 3.4 to 7.3) among individuals with ideal lifestyle, 6.6 (95% CI 5.4 to 8.1) among those with intermediate lifestyle and 10.4 (95% CI 9.2 to 11.8) among participants with poor lifestyle. On an additive scale, there was a positive statistically significant interaction between genetic risk and lifestyle (relative excess risk due to interaction=0.86, 95% CI 0.68 to 1.03, p<0.001).

CONCLUSIONS: The rates of AF increased gradually with worse lifestyle within each category of genetic risk. We found a positive interaction on an additive scale between genetic risk and lifestyle, suggesting that risk factor modification is especially important in individuals with a high genetic risk of AF.

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.

Mitophagy degrades damaged mitochondria, but we show here that it can also target functional mitochondria. This latter scenario occurs during programmed mitophagy and involves the mitophagy receptors NIX and BNIP3. Although AMP-activated protein kinase (AMPK), the energy-sensing protein kinase, can influence damaged-induced mitophagy, its role in programmed mitophagy is unclear. We found that AMPK directly inhibits NIX-dependent mitophagy by triggering 14-3-3-mediated sequestration of ULK1, via ULK1 phosphorylation at two sites: Ser556 and an additional identified site, Ser694. By contrast, AMPK activation increases Parkin phosphorylation and enhances the rate of depolarization-induced mitophagy, independently of ULK1. We show that this happens both in cultured cells and tissues in vivo, using the mito-QC mouse model. Our work unveils a mechanism whereby AMPK activation downregulates mitophagy of functional mitochondria but enhances that of dysfunctional/damaged ones.