👤 M Laakso

A key methodological challenge for genome-wide association studies is how to leverage haplotype diversity and allelic heterogeneity to improve trait association power, especially in noncoding regions where it is difficult to predict variant impacts and define functional units for variant aggregation. Genealogy-based association methods have the potential to bridge this gap by testing combinations of common and rare haplotypes based purely on their ancestral relationships. In parallel work, we have developed an efficient local ancestry inference engine and a novel statistical method (LOCATER) for combining signals present on different branches of a locus-specific haplotype tree. Here, we develop a genome-wide LOCATER analysis pipeline and apply it to a genome sequencing study of 6795 Finnish individuals with 101 cardiometabolic traits and 18.9 million autosomal variants. We identify 351 significant trait associations at 47 distinct genomic loci and find that LOCATER boosts the single marker test (SMT) association signal at five loci by combining independent signals from distinct alleles. LOCATER successfully recovers known quantitative trait loci not found by SMT, including Show less

To delineate organ-specific and systemic drivers of metabolic dysfunction-associated steatotic liver disease (MASLD), we applied integrative causal inference across clinical, imaging, and proteomic domains in individuals with and without type 2 diabetes (T2D). Bayesian network analyses and complementary two-sample Mendelian randomization were used to quantify causal pathways linking adipose distribution, glycemia, and insulin dynamics with liver fat in the IMI-DIRECT prospective cohort study. Data included frequently sampled metabolic challenge tests, MRI-derived abdominal and hepatic fat content, serological biomarkers, and Olink plasma proteomics from 331 adults with new-onset T2D and 964 adults without diabetes, with harmonized protocols enabling replication. High basal insulin secretion rate (BasalISR), estimated via C-peptide deconvolution, emerged as the primary potential causal driver of liver fat accumulation in both cohorts. BasalISR, a clearance-independent measure of β-cell insulin output distinct from peripheral insulin levels, was independently linked to hepatic steatosis. Visceral adipose tissue exhibited bidirectional associations with liver fat, suggesting a self-reinforcing metabolic loop. Of 446 analyzed proteins, 34 mapped to these metabolic networks (27 in the non-diabetes network, 18 in the T2D network, and 11 shared). Key proteins directly associated with liver fat included GUSB, ALDH1A1, LPL, IGFBP1/2, CTSD, HMOX1, FGF21, AGRP, and ACE2. Sex-stratified analyses identified GUSB in females and LEP in males as the strongest protein predictors of liver fat. BasalISR may better capture early β-cell-driven disturbances contributing to MASLD. These findings outline a multifactorial, sex- and disease stage-specific proteo-metabolic architecture of hepatic steatosis and identify potential biomarkers or therapeutic targets. Show less

Hyperlipidemia is known to impair endothelial function. We have recently shown that hyperlipidemia also blunts native post-ischemic capillary enlargement that is important for efficient skeletal muscle recovery from ischemia as it supports the recovery of arterial driving pressure and through intussusception increases capillary density. The correction of capillary reactivity under hyperlipidemia could, therefore, improve post-ischemic skeletal muscle recovery. This study tested the ability of adenoviral (Ad) vascular endothelial growth factor (VEGF) gene therapy to rescue capillary enlargement and improve post-ischemic muscle repair in hyperlipidemic mice. AdVEGF or AdLacZ-control vector were delivered into the calf muscles of aged, hyperlipidemic LDLR It was found that AdVEGF gene therapy was able to promote capillary enlargement (P < 0.05) that led to recovery of arterial driving pressure in ischemic LDLR Hyperlipidemia or old age did not seem to impair AdVEGF-induced capillary enlargement. However, regarding the side-effects of capillary enlargement, therapies trying to promote post-ischemic skeletal muscle recovery through angiogenesis should consider not only capillary size or density but also timing and dynamics of the capillary changes. Show less

A key methodological challenge for genome-wide association studies is how to leverage haplotype diversity and allelic heterogeneity to improve trait association power, especially in noncoding regions where it is difficult to predict variant impacts and define functional units for variant aggregation. Genealogy-based association methods have the potential to bridge this gap by testing combinations of common and rare haplotypes based purely on their ancestral relationships. In parallel work, we have developed an efficient local ancestry inference engine and a novel statistical method (LOCATER) for combining signals present on different branches of a locus specific haplotype tree. Here, we developed a genome-wide LOCATER analysis pipeline and applied it to a genome sequencing study of 6,795 Finnish individuals with 101 cardiometabolic traits and 18.9 million autosomal variants. We identify 351 significant trait associations at 47 distinct genomic loci and find that LOCATER boosts single marker test (SMT) association signal at 5 loci by combining independent signals from distinct alleles. LOCATER successfully recovers known quantitative trait loci not found by SMT, including Show less

To delineate organ-specific and systemic drivers of metabolic dysfunction-associated steatotic liver disease (MASLD), we applied integrative causal inference across clinical, imaging, and proteomic domains in individuals with and without type 2 diabetes (T2D). We used Bayesian network analyses to quantify causal pathways linking adipose distribution, glycemia, and insulin dynamics with fatty liver using data from the IMI-DIRECT prospective cohort study. Measurements were made of glucose and insulin dynamics (using frequently-sampled metabolic challenge tests), MRI-derived abdominal and liver fat content, serological biomarkers, and Olink plasma proteomics from 331 adults with new-onset T2D and 964 adults free from diabetes at enrolment. The common protocols used in these two cohorts provided the opportunity for replication analyses to be performed. When the direction of the effect could not be determined with high probability through Bayesian networks, complementary two-sample Mendelian randomization (MR) was employed. High basal insulin secretion rate (BasalISR) was identified as the primary causal driver of liver fat accumulation in both diabetes and non-diabetes. Excess visceral adipose tissue (VAT) was bidirectionally associated with liver fat, indicating a self-reinforcing metabolic loop. Basal insulin clearance (Clinsb) worsened as a consequence of liver fat accumulation to a greater degree before the onset of T2D. Out of 446 analysed proteins, 34 mapped to these metabolic networks and 27 were identified in the non-diabetes network, 18 in the diabetes network, and 11 were common between the two networks. Key proteins directly associated with liver fat included GUSB, ALDH1A1, LPL, IGFBP1/2, CTSD, HMOX1, FGF21, AGRP, and ACE2. Sex-stratified analyses revealed distinct proteomic drivers: GUSB and LEP were most predictive of liver fat in females and males, respectively. Basal insulin hypersecretion is a modifiable, causal driver of MASLD, particularly prior to glycaemic decompensation. Our findings highlight a multifactorial, sex- and disease-stage-specific proteo-metabolic architecture of hepatic steatosis. Proteins such as GUSB, ALDH1A1, LPL, and IGFBPs warrant further investigation as potential biomarkers or therapeutic targets for MASLD prevention and treatment. Show less

Fatty acid desaturase (FADS1) variant-rs174550 strongly regulates polyunsaturated fatty acid (PUFA) biosynthesis. Additionally, the FADS1 is related to mitochondrial function. Thus, we investigated whether changes in mitochondrial function are associated with the genetic variation in FADS1 (rs174550) in human adipocytes isolated from individuals consuming diets enriched with either dietary alpha-linolenic (ALA) or linoleic acid (LA). Two cohorts of men homozygous for the genotype of FADS1 (rs174550) were studied: FADSDIET2 dietary intervention study with ALA- and LA-enriched diets and Kuopio Obesity Surgery study (KOBS), respectively. We could demonstrate that differentiated human adipose-derived stromal cells from subjects with the TT genotype had higher mitochondrial metabolism compared with subjects with the CC genotype of FADS1-rs174550 in the FADSDIET2. Responses to PUFA-enriched diets differed between the genotypes of FADS1-rs174550, showing that ALA, but not LA, -enriched diet stimulated mitochondrial metabolism more in subjects with the CC genotype when compared with subjects with the TT genotype. ALA, but not LA, proportion in plasma phospholipid fraction correlated positively with adipose tissue mitochondrial-DNA amount in subjects with the CC genotype of FADS1-rs174550 in the KOBS. These findings demonstrate that the FADS1-rs174550 is associated with modification in mitochondrial function in human adipocytes. Additionally, subjects with the CC genotype, when compared with the TT genotype, benefit more from the ALA-enriched diet, leading to enhanced energy metabolism in human adipocytes. Altogether, the FADS1-rs174550 could be a genetic marker to identify subjects who are most suitable to receive dietary PUFA supplementation, establishing also a personalized therapeutic strategy to improve mitochondrial function in metabolic diseases. Show less

Lipoprotein(a) [Lp(a)] is a causal, genetically determined cardiovascular risk factor. Limited evidence suggests that dietary unsaturated fat may increase serum Lp(a) concentration by 10-15 %. Linoleic acid may increase Lp(a) concentration through its endogenous conversion to arachidonic acid, a process regulated by the fatty acid desaturase (FADS) gene cluster. We aimed to compare the Lp(a) and other lipoprotein trait-modulating effects of dietary alpha-linolenic (ALA) and linoleic acids (LA). Additionally, we examined whether FADS1 rs174550 genotype modifies Lp(a) responses. A genotype-based randomized trial was performed in 118 men homozygous for FADS1 rs174550 SNP (TT or CC). After a 4-week run-in period, the participants were randomized to 8-week intervention diets enriched with either Camelina sativa oil (ALA diet) or sunflower oil (LA diet) 30-50 mL/day based on their BMI. Serum lipid profile was measured at baseline and at the end of the intervention. ALA diet lowered serum Lp(a) concentration by 7.3 % (p = 0.003) and LA diet by 9.5 % (p < 0.001) (p = 0.089 for between-diet difference). Both diets led to greater absolute decreases in individuals with higher baseline Lp(a) concentration (p < 0.001). Concentrations of LDL cholesterol (LDL-C), non-HDL-C, remnant-C, and apolipoprotein B were lowered more by the ALA diet (p < 0.01). Lipid or lipoprotein responses were not modified by the FADS1 rs174550 genotype. A considerable increase in either dietary ALA or LA from vegetable oils has a similar Lp(a)-lowering effect, whereas ALA may lower other major atherogenic lipids and lipoproteins to a greater extent than LA. Genetic differences in endogenous PUFA conversion may not influence serum Lp(a) concentration. Show less

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

Previous genome-wide association studies (GWASs) for adiponectin, a complex trait linked to type 2 diabetes and obesity, identified >20 associated loci. However, most loci were identified in populations of European ancestry, and many of the target genes underlying the associations remain unknown. We conducted a cross-ancestry adiponectin GWAS meta-analysis in ≤46,434 individuals from the Metabolic Syndrome in Men (METSIM) cohort and the ADIPOGen and AGEN consortiums. We combined study-specific association summary statistics using a fixed-effects, inverse variance-weighted approach. We identified 22 loci associated with adiponectin (p < 5×10 Show less

BackgroundApolipoprotein C-III (apoC-III) is a regulator of triglyceride (TG) metabolism, and due to its association with risk of cardiovascular disease, is an emergent target for pharmacological intervention. The impact of substantially lowering apoC-III on lipoprotein metabolism is not clear.MethodsWe investigated the kinetics of apolipoproteins B48 and B100 (apoB48 and apoB100) in chylomicrons, VLDL1, VLDL2, IDL, and LDL in patients heterozygous for a loss-of-function (LOF) mutation in the APOC3 gene. Studies were conducted in the postprandial state to provide a more comprehensive view of the influence of this protein on TG transport.ResultsCompared with non-LOF variant participants, a genetically determined decrease in apoC-III resulted in marked acceleration of lipolysis of TG-rich lipoproteins (TRLs), increased removal of VLDL remnants from the bloodstream, and substantial decrease in circulating levels of VLDL1, VLDL2, and IDL particles. Production rates for apoB48-containing chylomicrons and apoB100-containing VLDL1 and VLDL2 were not different between LOF carriers and noncarriers. Likewise, the rate of production of LDL was not affected by the lower apoC-III level, nor were the concentration and clearance rate of LDL-apoB100.ConclusionThese findings indicate that apoC-III lowering will have a marked effect on TRL and remnant metabolism, with possibly significant consequences for cardiovascular disease prevention.Trial registrationClinicalTrials.gov NCT04209816 and NCT01445730.FundingSwedish Heart-Lung Foundation, Swedish Research Council, ALF grant from the Sahlgrenska University Hospital, Novo Nordisk Foundation, Sigrid Juselius Foundation, Helsinki University Hospital Government Research funds, Finnish Heart Foundation, and Finnish Diabetes Research Foundation. Show less

The fatty acid composition of plasma lipids, which is associated with biomarkers and risk of non-communicable diseases, is regulated by dietary polyunsaturated fatty acids (PUFAs) and variants of fatty acid desaturase (FADS). We investigated the interactions between dietary PUFAs and FADS1 rs174550 variant. Participants (n = 118), homozygous for FADS1 rs174550 variant (TT and CC) followed a high alpha-linolenic acid (ALA, 5 percent of energy (E-%)) or a high linoleic acid (LA, 10 E-%) diet during an 8-week randomized controlled intervention. Fatty acid composition of plasma lipids and PUFA-derived lipid mediators were quantified by gas and liquid chromatography mass spectrometry, respectively. The high-LA diet increased the concentration of plasma LA, but not its lipid mediators. The concentration of plasma arachidonic acid decreased in carriers of CC and remained unchanged in the TT genotype. The high-ALA diet increased the concentration of plasma ALA and its cytochrome P450-derived epoxides and dihydroxys, and cyclooxygenase-derived monohydroxys. Concentrations of plasma eicosapentaenoic acid and its mono- and dihydroxys increased only in TT genotype carriers. These findings suggest the potential for genotype-based recommendations for PUFA consumption, resulting in modulation of bioactive lipid mediators which can exert beneficial effects in maintaining health. Show less

Fatty acid desaturase (FADS) variants associate with fatty acid (FA) and adipose tissue (AT) metabolism and inflammation. Thus, the role of FADS1 variants in the regulation of dietary linoleic acid (LA)-induced effects on AT inflammation was investigated. Subjects homozygotes for the TT and CC genotypes of the FADS1-rs174550 (TT, n = 25 and CC, n = 28) or -rs174547 (TT, n = 42 and CC, n = 28), were either recruited from the METabolic Syndrome In Men cohort to participate in an intervention with LA-enriched diet (FADSDIET) or from the Kuopio Obesity Surgery (KOBS) study. GC and LC-MS for plasma FA proportions and eicosanoid concentrations and AT gene expression for AT inflammatory score (AT-InSc) was determined. We observed a diet-genotype interaction between LA-enriched diet and AT-InSc in the FADSDIET. In the KOBS study, interleukin (IL)1 beta mRNA expression in AT was increased in subjects with the TT genotype and highest LA proportion. In the FADSDIET, n-6/LA proportions correlated positively with AT-InSc in those with the TT genotype but not with the CC genotype after LA-enriched diet. Specifically, LA- and AA-derived pro-inflammatory eicosanoids related to CYP450/sEH-pathways correlated positively with AT-InSc in those with the TT genotype, whereas in those with the CC genotype, the negative correlations between pro-inflammatory eicosanoids and AT-InSc related to COX/LOX-pathways. LA-enriched diet increases inflammatory AT gene expression in subjects with the TT genotype, while CC genotype could play a protective role against LA-induced AT inflammation. Overall, the FADS1 variant could modify the dietary LA-induced effects on AT inflammation through the differential biosynthesis of AA-derived eicosanoids. Show less

Fatty acid desaturase 1 (FADS1) gene encodes for delta-5 desaturase enzyme which is needed in conversion of linoleic acid (LA) to arachidonic acid (AA). Recent studies have shown that response to dietary PUFAs differs between the genotypes in circulating fatty acids. However, interactions between the FADS1 genotype and dietary LA on overall metabolism have not been studied. We aimed to examine the interactions of FADS1 rs174550 genotypes (TT and CC) and high-LA diet to identify plasma metabolites that respond differentially to dietary LA according to the FADS1 genotype. A total of 59 men (TT n = 26, CC n = 33) consumed a sunflower oil supplemented diet for 4 weeks. Daily dose of 30, 40, or 50 ml was calculated based on body mass index. It resulted in 17-28 g of LA on top of the usual daily intake. Fasting plasma samples at the beginning and at the end of the intervention were analyzed with LC-MS/MS non-targeted metabolomics method. At the baseline, the carriers of FADS1 rs174550-TT genotype had higher abundance of long-chain PUFA phospholipids compared to the FADS1 rs174550-CC one. In response to the high-LA diet, LA phospholipids and long-chain acylcarnitines increased and lysophospholipids decreased in fasting plasma similarly in both genotypes. LysoPE (20:4), LysoPC (20:4), and PC (16:0₂₀:4) decreased and cortisol increased in the carriers of rs174550-CC genotype; however, these genotype-diet interactions were not significant after correction for multiple testing. Our findings show that both FADS1 rs174550 genotype and high-LA diet modify plasma phospholipid composition. The study was registered to ClinicalTrials: NCT02543216, September 7, 2015 (retrospectively registered). Show less

Hypertrophic (HCM) and dilated (DCM) cardiomyopathies are among the leading causes of sudden cardiac death. We identified 38 pathogenic or likely pathogenic variant carriers for HCM in three sarcomere genes ( Show less

Alternate splicing events can create isoforms that alter gene function, and genetic variants associated with alternate gene isoforms may reveal molecular mechanisms of disease. We used subcutaneous adipose tissue of 426 Finnish men from the METSIM study and identified splice junction quantitative trait loci (sQTLs) for 6,077 splice junctions (FDR < 1%). In the same individuals, we detected expression QTLs (eQTLs) for 59,443 exons and 15,397 genes (FDR < 1%). We identified 595 genes with an sQTL and exon eQTL but no gene eQTL, which could indicate potential isoform differences. Of the significant sQTL signals, 2,114 (39.8%) included at least one proxy variant (linkage disequilibrium r 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

The article investigates the FADS1 rs174550 genotype interaction with dietary intakes of high linoleic acid (LA) and high alpha-linolenic acid (ALA) on the response of fatty acid composition of plasma phospholipids (PLs), and of markers of low-grade inflammation and glucose-insulin homeostasis. One-hundred thirty homozygotes men for FADS1 rs174550 SNP (TT and CC genotypes) were randomized to an 8-week intervention with either LA- or ALA-enriched diet (13 E% PUFA). The source of LA and ALA are 30-50 mL of sunflower oil (SFO, 62-63% LA) and Camelina sativa oil (CSO, 30- are randomized to an 35% ALA), respectively. In the SFO arm, there is a significant genotype x diet interaction for the proportion of arachidonic acid in plasma phospholipids (p < 0.001), disposition index (DI The FADS1 genotype modifies the response to high PUFA diets, especially to high-LA diet. These findings suggest that approaches considering FADS variation may be useful in personalized dietary counseling. Show less

Lean body mass (LM) plays an important role in mobility and metabolic function. We previously identified five loci associated with LM adjusted for fat mass in kilograms. Such an adjustment may reduce the power to identify genetic signals having an association with both lean mass and fat mass. To determine the impact of different fat mass adjustments on genetic architecture of LM and identify additional LM loci. We performed genome-wide association analyses for whole-body LM (20 cohorts of European ancestry with n = 38,292) measured using dual-energy X-ray absorptiometry) or bioelectrical impedance analysis, adjusted for sex, age, age2, and height with or without fat mass adjustments (Model 1 no fat adjustment; Model 2 adjustment for fat mass as a percentage of body mass; Model 3 adjustment for fat mass in kilograms). Seven single-nucleotide polymorphisms (SNPs) in separate loci, including one novel LM locus (TNRC6B), were successfully replicated in an additional 47,227 individuals from 29 cohorts. Based on the strengths of the associations in Model 1 vs Model 3, we divided the LM loci into those with an effect on both lean mass and fat mass in the same direction and refer to those as "sumo wrestler" loci (FTO and MC4R). In contrast, loci with an impact specifically on LM were termed "body builder" loci (VCAN and ADAMTSL3). Using existing available genome-wide association study databases, LM increasing alleles of SNPs in sumo wrestler loci were associated with an adverse metabolic profile, whereas LM increasing alleles of SNPs in "body builder" loci were associated with metabolic protection. In conclusion, we identified one novel LM locus (TNRC6B). Our results suggest that a genetically determined increase in lean mass might exert either harmful or protective effects on metabolic traits, depending on its relation to fat mass. Show less

The health benefits of substituting dietary polyunsaturated fatty acids (PUFAs) for saturated fatty acids are well known. However, limited information exists on how the response to dietary intake of linoleic acid (LA; 18:2n-6) is modified by polymorphisms in the fatty acid desaturase (FADS) gene cluster. The aim of the current study was to test the hypothesis that the FADS1 rs174550 genotype modifies the effect of dietary LA intake on the fatty acid composition of plasma lipids, fasting glucose, and high-sensitivity C-reactive protein (hsCRP). Associations were investigated between genotype, plasma PUFAs, fasting glucose, and hsCRP concentrations in the cross-sectional, population-based Metabolic Syndrome in Men cohort (n = 1337). In addition, 62 healthy men from the cohort who were homozygotes for the TT or CC genotype of the FADS1 rs174550 were recruited to a 4-wk intervention (FADSDIET) with an LA-enriched diet. The fatty acid composition of plasma PUFAs and concentrations of plasma fasting glucose, serum hsCRP, and plasma lipid mediators (eicosanoids and related analogs) were measured at the beginning and end of the 4-wk intervention period. In the FADSDIET trial, the plasma LA proportion increased in both genotype groups in response to an LA-enriched diet. Responses in concentrations of serum hsCRP and plasma fasting glucose and the proportion of arachidonic acid (20:4n-6) in plasma phospholipids and cholesteryl esters differed between genotype groups (interaction of diet × genotype, P < 0.05). In TT homozygous subjects, plasma eicosanoid concentrations correlated with the arachidonic acid proportion in plasma and with hsCRP (r = 0.4-0.7, P < 0.05), whereas in the CC genotype there were no correlations. Our findings show that the FADS1 genotype modifies metabolic responses to dietary LA. The emerging concept that personalized dietary counseling should be modified by the FADS1 genotype needs to be tested in larger randomized trials. The study was registered at clinicaltrials.gov as NCT02543216. Show less

Nationwide large-scale genetic and outcome studies in cohorts with hypertrophic cardiomyopathy (HCM) have not been previously published. We sequenced 59 cardiomyopathy-associated genes in 382 unrelated Finnish patients with HCM and found 24 pathogenic or likely pathogenic mutations in six genes in 38.2% of patients. Most mutations were located in sarcomere genes (MYBPC3, MYH7, TPM1, and MYL2). Previously reported mutations by our study group (MYBPC3-Gln1061Ter, MYH7-Arg1053Gln, and TPM1-Asp175Asn) and a fourth major mutation MYH7-Val606Met accounted for 28.0% of cases. Mutations in GLA and PRKAG2 were found in three patients. Furthermore, we found 49 variants of unknown significance in 31 genes in 20.4% of cases. During a 6.7 ± 4.2 year follow-up, annual all-cause mortality in 482 index patients and their relatives with HCM was higher than that in the matched Finnish population (1.70 vs. 0.87%; P < 0.001). Sudden cardiac deaths were rare (n = 8). Systolic heart failure (hazard ratio 17.256, 95% confidence interval 3.266-91.170, P = 0.001) and maximal left ventricular wall thickness (hazard ratio 1.223, 95% confidence interval 1.098-1.363, P < 0.001) were independent predictors of HCM-related mortality and life-threatening cardiac events. The patients with a pathogenic or likely pathogenic mutation underwent an implantable cardioverter defibrillator implantation more often than patients without a pathogenic or likely pathogenic mutation (12.9 vs. 3.5%, P < 0.001), but there was no difference in all-cause or HCM-related mortality between the two groups. Mortality due to HCM during 10 year follow-up among the 5.2 million population of Finland was studied from death certificates of the National Registry, showing 269 HCM-related deaths, of which 32% were sudden. We identified pathogenic and likely pathogenic mutations in 38% of Finnish patients with HCM. Four major sarcomere mutations accounted for 28% of HCM cases, whereas HCM-related mutations in non-sarcomeric genes were rare. Mortality in patients with HCM exceeded that of the general population. Finally, among 5.2 million Finns, there were at least 27 HCM-related deaths annually. 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

Increased adiposity is a hallmark of obesity and overweight, which affect 2.2 billion people world-wide. Understanding the genetic and molecular mechanisms that underlie obesity-related phenotypes can help to improve treatment options and drug development. Here we perform promoter Capture Hi-C in human adipocytes to investigate interactions between gene promoters and distal elements as a transcription-regulating mechanism contributing to these phenotypes. We find that promoter-interacting elements in human adipocytes are enriched for adipose-related transcription factor motifs, such as PPARG and CEBPB, and contribute to heritability of cis-regulated gene expression. We further intersect these data with published genome-wide association studies for BMI and BMI-related metabolic traits to identify the genes that are under genetic cis regulation in human adipocytes via chromosomal interactions. This integrative genomics approach identifies four cis-eQTL-eGene relationships associated with BMI or obesity-related traits, including rs4776984 and MAP2K5, which we further confirm by EMSA, and highlights 38 additional candidate genes. Show less

The sensitivity and specificity of the conventional 12-lead ECG to identify carriers of hypertrophic cardiomyopathy (HCM) - causing mutations without left ventricular hypertrophy (LVH) has been limited. We assessed the ability of novel electrocardiographic parameters to improve the detection of HCM mutation carriers. We studied 140 carriers (G+) of the TPM1-Asp175Asn or MYBPC3-Gln1061X pathogenic variants for HCM: The G+/LVH+ group (n = 98) consisted of mutation carriers with LVH and the G+/LVH- group (n = 42) without LVH. The control group consisted of 30 subjects. The standard 12-lead ECG was comprehensively analyzed and two novel ECG variables were introduced: RV1RV3 and septal remodeling. A subset of 65 individuals underwent cardiac magnetic resonance imaging and 2D strain echocardiography. Conventional major ECG criteria were sensitive (90%) and specific (97%) in identifying G+/LVH+ subjects. RV1RV3 and septal remodeling were more prevalent in the G+/LVH- subjects compared to the control group (33% vs 3%, p = 0.005 and 45% vs 3%, p < 0.001, respectively). The combination of RV1RV3 and Q waves and repolarization abnormalities (QR) differentiated between the G+/LVH- subjects and the control group with a sensitivity of 52% and specificity of 97%. The combination of septal remodeling and QR differentiated between G+/LVH- subjects and the control group with a sensitivity of 64% and specificity of 97%. The novel ECG-parameters RV1RV3 and septal remodeling were effective in identifying G+/LVH- subjects and could be useful in the diagnostics of new suspected HCM patients and in the screening and follow-up of HCM families. 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

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

Glucose-dependent insulinotropic polypeptide (GIP) is an incretin hormone with extrapancreatic effects beyond glycemic control. Here we demonstrate unexpected effects of GIP signaling in the vasculature. GIP induces the expression of the proatherogenic cytokine osteopontin (OPN) in mouse arteries via local release of endothelin-1 and activation of CREB. Infusion of GIP increases plasma OPN concentrations in healthy individuals. Plasma endothelin-1 and OPN concentrations are positively correlated in patients with critical limb ischemia. Fasting GIP concentrations are higher in individuals with a history of cardiovascular disease (myocardial infarction, stroke) when compared with control subjects. GIP receptor (GIPR) and OPN mRNA levels are higher in carotid endarterectomies from patients with symptoms (stroke, transient ischemic attacks, amaurosis fugax) than in asymptomatic patients, and expression associates with parameters that are characteristic of unstable and inflammatory plaques (increased lipid accumulation, macrophage infiltration, and reduced smooth muscle cell content). While GIPR expression is predominantly endothelial in healthy arteries from humans, mice, rats, and pigs, remarkable upregulation is observed in endothelial and smooth muscle cells upon culture conditions, yielding a "vascular disease-like" phenotype. Moreover, the common variant rs10423928 in the GIPR gene is associated with increased risk of stroke in patients with type 2 diabetes. 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

Limited information exists on how the relationship between dietary intake of fat and fatty acids in erythrocytes and plasma is modulated by polymorphisms in the FADS gene cluster. We examined gene-diet interaction of total marine PUFA intake with a known gene encoding Δ-5 desaturase enzyme (FADS1) variant (rs174550) for fatty acids in erythrocyte membranes and plasma phospholipids (PL), cholesteryl esters (CE), and triglycerides (TG). In this cross-sectional study, fatty acid compositions were measured using GC, and total intake of polyunsaturated fat from fish and fish oil was estimated using a food frequency questionnaire in a subsample (n = 962) of the Metabolic Syndrome in Men Study. We found nominally significant gene-diet interactions for eicosapentaenoic acid (EPA, 20:5n-3) in erythrocytes (pinteraction = 0.032) and for EPA in plasma PL (pinteraction = 0.062), CE (pinteraction = 0.035), and TG (pinteraction = 0.035), as well as for docosapentaenoic acid (22:5n-3) in PL (pinteraction = 0.007). After excluding omega-3 supplement users, we found a significant gene-diet interaction for EPA in erythrocytes (pinteraction < 0.003). In a separate cohort of the Kuopio Obesity Surgery Study, the same locus was strongly associated with hepatic mRNA expression of FADS1 (p = 1.5 × 10(-10) ). FADS1 variants may modulate the relationship between marine fatty acid intake and circulating levels of long-chain omega-3 fatty acids. Show less

We assessed the value of speckle tracking two-dimensional (2D) strain echocardiography (2DSE) measured mechanical dispersion (MD) with other imaging and electrocardiographic parameters in differentiating hypertrophic cardiomyopathy (HCM) patients with and without nonsustained ventricular tachycardia (NSVT) on 24-h ambulatory ECG monitoring. We studied 31 patients with HCM caused by the Finnish founder mutation MYBPC3-Q1061X and 20 control subjects with comprehensive 2DSE echocardiography and cardiac magnetic resonance imaging (CMRI). The presence of NSVT was assessed from ambulatory 24-h ECG monitoring. NSVT episodes were recorded in 11 (35%) patients with HCM. MD was significantly higher in HCM patients with NSVT (93 ± 41 ms) compared to HCM patients without NSVT (50 ± 18 ms, p = 0.012) and control subjects (41 ± 16 ms, p < 0.001). MD was the only variable independently associated with the presence of NSVT (OR: 1.60, 95% CI: 1.05-2.45, p = 0.030). Assessed by ROC curves, MD performed best in differentiating between HCM patients with and without NSVT (AUC = 0.81). Increased mechanical dispersion was associated with NSVT in HCM patients on 24-h ambulatory ECG monitoring. Key messages The prediction of sudden cardiac death in hypertrophic cardiomyopathy remains a challenge and novel imaging methods are required to identify individuals at risk of malignant ventricular arrhythmias. Mechanical dispersion by speckle tracking echocardiography is associated with NSVT on 24-h ambulatory ECG monitoring in patients with hypertrophic cardiomyopathy. Show less

Mutations in the cardiac myosin-binding protein C gene (MYBPC3) are the most common genetic cause of hypertrophic cardiomyopathy (HCM) worldwide. The molecular mechanisms leading to HCM are poorly understood. We investigated the metabolic profiles of mutation carriers with the HCM-causing MYBPC3-Q1061X mutation with and without left ventricular hypertrophy (LVH) and non-affected relatives, and the association of the metabolome to the echocardiographic parameters. 34 hypertrophic subjects carrying the MYBPC3-Q1061X mutation, 19 non-hypertrophic mutation carriers and 20 relatives with neither mutation nor hypertrophy were examined using comprehensive echocardiography. Plasma was analyzed for molecular lipids and polar metabolites using two metabolomics platforms. Concentrations of branched chain amino acids, triglycerides and ether phospholipids were increased in mutation carriers with hypertrophy as compared to controls and non-hypertrophic mutation carriers, and correlated with echocardiographic LVH and signs of diastolic and systolic dysfunction in subjects with the MYBPC3-Q1061X mutation. Our study implicates the potential role of branched chain amino acids, triglycerides and ether phospholipids in HCM, as well as suggests an association of these metabolites with remodeling and dysfunction of the left ventricle. Show less