Pirkka-Pekka Laurila

Novel Loci for Metabolic Networks and Multi-Tissue Expression Studies Reveal Genes for Atherosclerosis

Association testing of multiple correlated phenotypes offers better power than univariate analysis of single traits. We analyzed 6,600 individuals from two population-based cohorts with both genome-wide SNP data and serum metabolomic profiles. From the observed correlation structure of 130 metabolites measured by nuclear magnetic resonance, we identified 11 metabolic networks and performed a multivariate genome-wide association analysis. We identified 34 genomic loci at genome-wide significance, of which 7 are novel. In comparison to univariate tests, multivariate association analysis identified nearly twice as many significant associations in total. Multi-tissue gene expression studies identified variants in our top loci, SERPINA1 and AQP9, as eQTLs and showed that SERPINA1 and AQP9 expression in human blood was associated with metabolites from their corresponding metabolic networks. Finally, liver expression of AQP9 was associated with atherosclerotic lesion area in mice, and in human arterial tissue both SERPINA1 and AQP9 were shown to be upregulated (6.3-fold and 4.6-fold, respectively) in atherosclerotic plaques. Our study illustrates the power of multi-phenotype GWAS and highlights candidate genes for atherosclerosis.

Genetic Association and Interaction Analysis of USF1 and APOA5 on Lipid Levels and Atherosclerosis

Objective—USF1 is a ubiquitous transcription factor governing the expression of numerous genes of lipid and glucose metabolism. APOA5 is a well-established candidate gene regulating triglyceride (TG) levels and has been identified as a downstream target of upstream stimulatory factor. No detailed studies about the effect of APOA5 on atherosclerotic lesion formation have been conducted, nor has its potential interaction with USF1 been examined. Methods and Results—We analyzed allelic variants of USF1 and APOA5 in families (n=516) ascertained for atherogenic dyslipidemia and in an autopsy series of middle-aged men (n=300) with precise quantitative measurements of atherosclerotic lesions. The impact of previously associated APOA5 variants on TGs was observed in the dyslipidemic families, and variant rs3135506 was associated with size of fibrotic aortic lesions in the autopsy series. The USF1 variant rs2516839, associated previously with atherosclerotic lesions, showed an effect on TGs in members of the dyslipidemic families with documented coronary artery disease. We provide preliminary evidence of gene-gene interaction between these variants in an autopsy series with a fibrotic lesion area in the abdominal aorta (P=0.0028), with TGs in dyslipidemic coronary artery disease subjects (P=0.03), and with high-density lipoprotein cholesterol (P=0.008) in a large population cohort of coronary artery disease patients (n=1065) in which the interaction for TGs was not replicated. Conclusion—Our findings in these unique samples reinforce the roles of APOA5 and USF1 variants on cardiovascular phenotypes and suggest that both genes contribute to lipid levels and aortic atherosclerosis individually and possibly through epistatic effects.

Adipose tissue gene expression analysis reveals changes in inflammatory, mitochondrial respiratory and lipid metabolic pathways in obese insulin-resistant subjects

BackgroundTo get insight into molecular mechanisms underlying insulin resistance, we compared acute in vivo effects of insulin on adipose tissue transcriptional profiles between obese insulin-resistant and lean insulin-sensitive women.MethodsSubcutaneous adipose tissue biopsies were obtained before and after 3 and 6 hours of intravenously maintained euglycemic hyperinsulinemia from 9 insulin-resistant and 11 insulin-sensitive females. Gene expression was measured using Affymetrix HG U133 Plus 2 microarrays and qRT-PCR. Microarray data and pathway analyses were performed with Chipster v1.4.2 and by using in-house developed nonparametric pathway analysis software.ResultsThe most prominent difference in gene expression of the insulin-resistant group during hyperinsulinemia was reduced transcription of nuclear genes involved in mitochondrial respiration (mitochondrial respiratory chain, GO:0001934). Inflammatory pathways with complement components (inflammatory response, GO:0006954) and cytokines (chemotaxis, GO:0042330) were strongly up-regulated in insulin-resistant as compared to insulin-sensitive subjects both before and during hyperinsulinemia. Furthermore, differences were observed in genes contributing to fatty acid, cholesterol and triglyceride metabolism (FATP2, ELOVL6, PNPLA3, SREBF1) and in genes involved in regulating lipolysis (ANGPTL4) between the insulin-resistant and -sensitive subjects especially during hyperinsulinemia.ConclusionsThe major finding of this study was lower expression of mitochondrial respiratory pathway and defective induction of lipid metabolism pathways by insulin in insulin-resistant subjects. Moreover, the study reveals several novel genes whose aberrant regulation is associated with the obese insulin-resistant phenotype.

Genomic, Transcriptomic, and Lipidomic Profiling Highlights the Role of Inflammation in Individuals With Low High-density Lipoprotein Cholesterol

Objective—Low high-density lipoprotein cholesterol (HDL-C) is associated with cardiometabolic pathologies. In this study, we investigate the biological pathways and individual genes behind low HDL-C by integrating results from 3 high-throughput data sources: adipose tissue transcriptomics, HDL lipidomics, and dense marker genotypes from Finnish individuals with low or high HDL-C (n=450). Approach and Results—In the pathway analysis of genetic data, we demonstrate that genetic variants within inflammatory pathways were enriched among low HDL-C associated single-nucleotide polymorphisms, and the expression of these pathways upregulated in the adipose tissue of low HDL-C subjects. The lipidomic analysis highlighted the change in HDL particle quality toward putatively more inflammatory and less vasoprotective state in subjects with low HDL-C, as evidenced by their decreased antioxidative plasmalogen contents. We show that the focal point of these inflammatory pathways seems to be the HLA region with its low HDL-associated alleles also associating with more abundant local transcript levels in adipose tissue, increased plasma vascular cell adhesion molecule 1 (VCAM1) levels, and decreased HDL particle plasmalogen contents, markers of adipose tissue inflammation, vascular inflammation, and HDL antioxidative potential, respectively. In a population-based look-up of the inflammatory pathway single-nucleotide polymorphisms in a large Finnish cohorts (n=11 211), no association of the HLA region was detected for HDL-C as quantitative trait, but with extreme HDL-C phenotypes, implying the presence of low or high HDL genes in addition to the population-genomewide association studies–identified HDL genes. Conclusions—Our study highlights the role of inflammation with a genetic component in subjects with low HDL-C and identifies novel cis-expression quantitative trait loci (cis-eQTL) variants in HLA region to be associated with low HDL-C.

A Genome-Wide Screen for Interactions Reveals a New Locus on 4p15 Modifying the Effect of Waist-to-Hip Ratio on Total Cholesterol

Recent genome-wide association (GWA) studies described 95 loci controlling serum lipid levels. These common variants explain ∼25% of the heritability of the phenotypes. To date, no unbiased screen for gene–environment interactions for circulating lipids has been reported. We screened for variants that modify the relationship between known epidemiological risk factors and circulating lipid levels in a meta-analysis of genome-wide association (GWA) data from 18 population-based cohorts with European ancestry (maximum N = 32,225). We collected 8 further cohorts (N = 17,102) for replication, and rs6448771 on 4p15 demonstrated genome-wide significant interaction with waist-to-hip-ratio (WHR) on total cholesterol (TC) with a combined P-value of 4.79×10−9. There were two potential candidate genes in the region, PCDH7 and CCKAR, with differential expression levels for rs6448771 genotypes in adipose tissue. The effect of WHR on TC was strongest for individuals carrying two copies of G allele, for whom a one standard deviation (sd) difference in WHR corresponds to 0.19 sd difference in TC concentration, while for A allele homozygous the difference was 0.12 sd. Our findings may open up possibilities for targeted intervention strategies for people characterized by specific genomic profiles. However, more refined measures of both body-fat distribution and metabolic measures are needed to understand how their joint dynamics are modified by the newly found locus.

USF1 deficiency activates brown adipose tissue and improves cardiometabolic health

Deficiency of USF1 protects against obesity, insulin resistance, and cardiovascular disease in mice and humans, and induces brown adipose tissue to burn triglycerides and glucose. Boosting metabolism with low USF1 Laurila et al. show that a deficiency of the transcription factor USF1 protects against obesity, insulin resistance, and cardiovascular disease. Even when fed a high-fat diet, USF1-deficient mice stayed lean and maintained a beneficial lipid profile with low triglycerides and high high-density lipoprotein (HDL) cholesterol. The mice had elevated energy expenditure because their brown adipose tissue was more active. In USF1-deficient mice, glucose and lipids were rapidly cleared from the circulation to be burned by brown fat. In humans, individuals with reduced USF1 expression also had higher plasma HDL cholesterol and lower triglycerides, and were more insulin-sensitive and less prone to hardening of the arteries. These findings identify USF1 as a potential therapeutic target for treating metabolic and cardiac diseases. USF1 (upstream stimulatory factor 1) is a transcription factor associated with familial combined hyperlipidemia and coronary artery disease in humans. However, whether USF1 is beneficial or detrimental to cardiometabolic health has not been addressed. By inactivating USF1 in mice, we demonstrate protection against diet-induced dyslipidemia, obesity, insulin resistance, hepatic steatosis, and atherosclerosis. The favorable plasma lipid profile, including increased high-density lipoprotein cholesterol and decreased triglycerides, was coupled with increased energy expenditure due to activation of brown adipose tissue (BAT). Usf1 inactivation directs triglycerides from the circulation to BAT for combustion via a lipoprotein lipase–dependent mechanism, thus enhancing plasma triglyceride clearance. Mice lacking Usf1 displayed increased BAT-facilitated, diet-induced thermogenesis with up-regulation of mitochondrial respiratory chain complexes, as well as increased BAT activity even at thermoneutrality and after BAT sympathectomy. A direct effect of USF1 on BAT activation was demonstrated by an amplified adrenergic response in brown adipocytes after Usf1 silencing, and by augmented norepinephrine-induced thermogenesis in mice lacking Usf1. In humans, individuals carrying SNP (single-nucleotide polymorphism) alleles that reduced USF1 mRNA expression also displayed a beneficial cardiometabolic profile, featuring improved insulin sensitivity, a favorable lipid profile, and reduced atherosclerosis. Our findings identify a new molecular link between lipid metabolism and energy expenditure, and point to the potential of USF1 as a therapeutic target for cardiometabolic disease.

The Contribution of GWAS Loci in Familial Dyslipidemias

Familial combined hyperlipidemia (FCH) is a complex and common familial dyslipidemia characterized by elevated total cholesterol and/or triglyceride levels with over five-fold risk of coronary heart disease. The genetic architecture and contribution of rare Mendelian and common variants to FCH susceptibility is unknown. In 53 Finnish FCH families, we genotyped and imputed nine million variants in 715 family members with DNA available. We studied the enrichment of variants previously implicated with monogenic dyslipidemias and/or lipid levels in the general population by comparing allele frequencies between the FCH families and population samples. We also constructed weighted polygenic scores using 212 lipid-associated SNPs and estimated the relative contributions of Mendelian variants and polygenic scores to the risk of FCH in the families. We identified, across the whole allele frequency spectrum, an enrichment of variants known to elevate, and a deficiency of variants known to lower LDL-C and/or TG levels among both probands and affected FCH individuals. The score based on TG associated SNPs was particularly high among affected individuals compared to non-affected family members. Out of 234 affected FCH individuals across the families, seven (3%) carried Mendelian variants and 83 (35%) showed high accumulation of either known LDL-C or TG elevating variants by having either polygenic score over the 90th percentile in the population. The positive predictive value of high score was much higher for affected FCH individuals than for similar sporadic cases in the population. FCH is highly polygenic, supporting the hypothesis that variants across the whole allele frequency spectrum contribute to this complex familial trait. Polygenic SNP panels improve identification of individuals affected with FCH, but their clinical utility remains to be defined.