👤 Ville Männistö

Abdominal obesity increases the risk for non-alcoholic fatty liver disease (NAFLD), now known as metabolic dysfunction-associated steatotic liver disease (MASLD). To elucidate the directional cell-type level biological mechanisms underlying the association between abdominal obesity and MASLD, we integrated adipose and liver single nucleus RNA-sequencing and bulk cis-expression quantitative trait locus (eQTL) data with the UK Biobank genome-wide association study (GWAS) data using colocalization. Then we used colocalized cis-eQTL variants as instrumental variables in Mendelian randomization (MR) analyses, followed by functional validation experiments on the target genes of the cis-eQTL variants. We identified 17 colocalized abdominal obesity GWAS variants, regulating 17 adipose cell-type marker genes. Incorporating these 17 variants into MR discovers a putative tissue-of-origin, cell-type-aware causal effect of abdominal obesity on MASLD consistently with multiple MR methods without significant evidence for pleiotropy or heterogeneity. Single cell data confirm the adipocyte-enriched mean expression of the 17 genes. Our cellular experiments across human adipogenesis identify risk variant -specific epigenetic and transcriptional mechanisms. Knocking down two of the 17 genes, PPP2R5A and SH3PXD2B, shows a marked decrease in adipocyte lipidation and significantly alters adipocyte function and adipogenesis regulator genes, including DGAT2, LPL, ADIPOQ, PPARG, and SREBF1. Furthermore, the 17 genes capture a characteristic MASLD expression signature in subcutaneous adipose tissue. Overall, we discover a significant cell-type level effect of abdominal obesity on MASLD and trace its biological effect to adipogenesis. NIH grants R01HG010505, R01DK132775, and R01HL170604; the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation program (Grant No. 802825), Academy of Finland (Grants Nos. 333021), the Finnish Foundation for Cardiovascular Research the Sigrid Jusélius Foundation and the Jane and Aatos Erkko Foundation; American Association for the Study of Liver Diseases (AASLD) Advanced Transplant Hepatology award and NIH/NIDDK (P30DK41301) Pilot and Feasibility award; NIH/NIEHS F32 award (F32ES034668); Finnish Diabetes Research Foundation, Kuopio University Hospital Project grant (EVO/VTR grants 2005-2021), the Academy of Finland grant (Contract no. 138006); Academy of Finland (Grant Nos 335443, 314383, 272376 and 266286), Sigrid Jusélius Foundation, Finnish Medical Foundation, Finnish Diabetes Research Foundation, Novo Nordisk Foundation (#NNF20OC0060547, NNF17OC0027232, NNF10OC1013354) and Government Research Funds to Helsinki University Hospital; Orion Research Foundation, Maud Kuistila Foundation, Finish Medical Foundation, and University of Helsinki. Show less

Non-alcoholic fatty liver disease has been associated with increased mRNA expression of FADS2 in the liver and estimated activity of delta-6 desaturase in serum, encoded by the FADS2 gene. Since DNA methylation in the FADS1/2/3 gene cluster has been previously linked with genetic variants and desaturase activities, we now aimed to discover factors regulating DNA methylation of the CpG sites annotated to FADS1/2 genes. DNA methylation levels in the CpG sites annotated to FADS2 and FADS1 were analyzed from liver samples of 95 obese participants of the Kuopio Obesity Surgery Study (34 men and 61 women, age 49.5 ± 7.7 years, BMI 43.0 ± 5.7 kg/m We found a negative correlation between DNA methylation levels of cg06781209 and cg07999042 and hepatic FADS2 mRNA expression (both p < 0.05), and with estimated delta-6 desaturase activity based on both liver and serum fatty acids (all p < 0.05). Interestingly, the methylation level of cg07999042 (p = 0.001) but not of cg06781209 (p = 0.874) was associated with FADS2 variant rs174616. Genetic variants of FADS2 may contribute to the pathogenesis of non-alcoholic fatty liver disease by modifying DNA methylation. 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

Non-alcoholic steatohepatitis (NASH) is associated with changes in fatty acid (FA) metabolism. However, specific changes in metabolism and hepatic mRNA expression related to NASH independent of simple steatosis, obesity and diet are unknown. Liver histology, serum and liver FA composition and estimated enzyme activities based on the FA ratios in cholesteryl esters and triglycerides were assessed in 92 obese participants of the Kuopio Obesity Surgery Study (KOBS) divided to those with normal liver, steatosis or NASH (30 men and 62 women, age 46.8±9.5years (mean±SD), BMI 44.2±6.2kg/m(2)). Plasma FA composition was also investigated in the Metabolic Syndrome in Men (METSIM) Study (n=769), in which serum alanine aminotransferase (ALT) was used as a marker of liver disease. Obese individuals with NASH had higher activity of estimated activities of delta-6 desaturase (D6D, p<0.002) and stearoyl-CoA desaturase 1 (SCD1, p<0.002) and lower activity of delta-5 desaturase (D5D, p<0.002) when compared to individuals with normal liver. Estimated activities of D5D, D6D and SCD1 correlated positively between liver and serum indicating that serum estimates reflected liver metabolism. Accordingly, NASH was associated with higher hepatic mRNA expression of corresponding genes FADS1, FADS2 and SCD. Finally, differences in FA metabolism that associated with NASH in obese individuals were also associated with high ALT in the METSIM Study. We demonstrated alterations in FA metabolism and endogenous desaturase activities that associate with NASH, independent of obesity and diet. This suggests that changes in endogenous FA metabolism are related to NASH and that they may contribute to the progression of the disease. Show less

Obesity is associated with disturbed lipid metabolism and low-grade inflammation in tissues. The aim of this study was to investigate the association between FA metabolism and adipose tissue (AT) inflammation in the Kuopio Obesity Surgery study. We investigated the association of surgery-induced weight loss and FA desaturase (FADS)1/2 genotypes with serum and AT FA profile and with AT inflammation, measured as interleukin (IL)-1β and NFκB pathway gene expression, in order to find potential gene-environment interactions. We demonstrated an association between serum levels of saturated and polyunsaturated n-6 FAs, and estimated enzyme activities of FADS1/2 genes with IL-1β expression in AT both at baseline and at follow-up. Variation in the FADS1/2 genes associated with IL-1β and NFκB pathway gene expression in SAT after weight reduction, but not at baseline. In addition, the FA composition in subcutaneous and visceral fat correlated with serum FAs, and the associations between serum PUFAs and estimated D6D enzyme activity with AT inflammation were also replicated with corresponding AT FAs and AT inflammation. We conclude that the polymorphism in FADS1/2 genes associates with FA metabolism and AT inflammation, leading to an interaction between weight loss and FADS1/2 genes in the regulation of AT inflammation. Show less

Coffee, a major dietary source of caffeine, is among the most widely consumed beverages in the world and has received considerable attention regarding health risks and benefits. We conducted a genome-wide (GW) meta-analysis of predominately regular-type coffee consumption (cups per day) among up to 91,462 coffee consumers of European ancestry with top single-nucleotide polymorphisms (SNPs) followed-up in ~30 062 and 7964 coffee consumers of European and African-American ancestry, respectively. Studies from both stages were combined in a trans-ethnic meta-analysis. Confirmed loci were examined for putative functional and biological relevance. Eight loci, including six novel loci, met GW significance (log10Bayes factor (BF)>5.64) with per-allele effect sizes of 0.03-0.14 cups per day. Six are located in or near genes potentially involved in pharmacokinetics (ABCG2, AHR, POR and CYP1A2) and pharmacodynamics (BDNF and SLC6A4) of caffeine. Two map to GCKR and MLXIPL genes related to metabolic traits but lacking known roles in coffee consumption. Enhancer and promoter histone marks populate the regions of many confirmed loci and several potential regulatory SNPs are highly correlated with the lead SNP of each. SNP alleles near GCKR, MLXIPL, BDNF and CYP1A2 that were associated with higher coffee consumption have previously been associated with smoking initiation, higher adiposity and fasting insulin and glucose but lower blood pressure and favorable lipid, inflammatory and liver enzyme profiles (P<5 × 10(-8)).Our genetic findings among European and African-American adults reinforce the role of caffeine in mediating habitual coffee consumption and may point to molecular mechanisms underlying inter-individual variability in pharmacological and health effects of coffee. Show less