Leo-Pekka Lyytikäinen

miR-21, miR-210, miR-34a, and miR-146a/b are up-regulated in human atherosclerotic plaques in the Tampere Vascular Study

OBJECTIVE MicroRNAs are small non-coding RNAs that inversely regulate their target gene expression. The whole miRNA profile of human atherosclerotic plaques has not been studied previously. The aim of this study was to investigate the miRNA expression profile in human atherosclerotic plaques as compared to non-atherosclerotic left internal thoracic arteries (LITA), and to connect this expression to the processes in atherosclerosis. METHODS The miRNA expression profiles of six LITAs and 12 atherosclerotic plaques obtained from aortic, carotid, and femoral atherosclerotic arteries from Tampere Vascular Study were analyzed. The analyses were performed with Agilents miRNA Microarray. The expression levels of over 4-fold up-regulated miRNAs were verified with qRT-PCR from a larger population (n=50). Messenger RNA levels were analyzed with Illuminas Expression BeadChip to study miRNA target expression. RESULTS Ten miRNAs were found to be differently expressed in atherosclerotic plaques when compared to controls (p<0.05). The expression of miR-21, -34a, -146a, -146b-5p, and -210 was verified and found to be significantly up-regulated in atherosclerotic arteries versus LITAs (p<0.001, fold changes 4.61, 2.55, 2.87, 2.82, and 3.92, respectively). Several predicted targets of these miRNAs were down-regulated, and gene set enrichment analysis showed several pathways which could be differently expressed due to this miRNA profile. CONCLUSIONS The microRNA expression profile differs significantly between atherosclerotic plaques and control arteries. The most up-regulated miRNAs are involved in processes known to be connected to atherosclerosis. Interfering with the miRNA expression in the artery wall is a potential way to affect atherosclerotic plaque and cardiovascular disease development.

A Genome-Wide Association Meta-Analysis of Circulating Sex Hormone–Binding Globulin Reveals Multiple Loci Implicated in Sex Steroid Hormone Regulation

Sex hormone-binding globulin (SHBG) is a glycoprotein responsible for the transport and biologic availability of sex steroid hormones, primarily testosterone and estradiol. SHBG has been associated with chronic diseases including type 2 diabetes (T2D) and with hormone-sensitive cancers such as breast and prostate cancer. We performed a genome-wide association study (GWAS) meta-analysis of 21,791 individuals from 10 epidemiologic studies and validated these findings in 7,046 individuals in an additional six studies. We identified twelve genomic regions (SNPs) associated with circulating SHBG concentrations. Loci near the identified SNPs included SHBG (rs12150660, 17p13.1, p = 1.8×10−106), PRMT6 (rs17496332, 1p13.3, p = 1.4×10−11), GCKR (rs780093, 2p23.3, p = 2.2×10−16), ZBTB10 (rs440837, 8q21.13, p = 3.4×10−09), JMJD1C (rs7910927, 10q21.3, p = 6.1×10−35), SLCO1B1 (rs4149056, 12p12.1, p = 1.9×10−08), NR2F2 (rs8023580, 15q26.2, p = 8.3×10−12), ZNF652 (rs2411984, 17q21.32, p = 3.5×10−14), TDGF3 (rs1573036, Xq22.3, p = 4.1×10−14), LHCGR (rs10454142, 2p16.3, p = 1.3×10−07), BAIAP2L1 (rs3779195, 7q21.3, p = 2.7×10−08), and UGT2B15 (rs293428, 4q13.2, p = 5.5×10−06). These genes encompass multiple biologic pathways, including hepatic function, lipid metabolism, carbohydrate metabolism and T2D, androgen and estrogen receptor function, epigenetic effects, and the biology of sex steroid hormone-responsive cancers including breast and prostate cancer. We found evidence of sex-differentiated genetic influences on SHBG. In a sex-specific GWAS, the loci 4q13.2-UGT2B15 was significant in men only (men p = 2.5×10−08, women p = 0.66, heterogeneity p = 0.003). Additionally, three loci showed strong sex-differentiated effects: 17p13.1-SHBG and Xq22.3-TDGF3 were stronger in men, whereas 8q21.12-ZBTB10 was stronger in women. Conditional analyses identified additional signals at the SHBG gene that together almost double the proportion of variance explained at the locus. Using an independent study of 1,129 individuals, all SNPs identified in the overall or sex-differentiated or conditional analyses explained ∼15.6% and ∼8.4% of the genetic variation of SHBG concentrations in men and women, respectively. The evidence for sex-differentiated effects and allelic heterogeneity highlight the importance of considering these features when estimating complex trait variance.

WNT16 influences bone mineral density, cortical bone thickness, bone strength, and osteoporotic fracture risk.

We aimed to identify genetic variants associated with cortical bone thickness (CBT) and bone mineral density (BMD) by performing two separate genome-wide association study (GWAS) meta-analyses for CBT in 3 cohorts comprising 5,878 European subjects and for BMD in 5 cohorts comprising 5,672 individuals. We then assessed selected single-nucleotide polymorphisms (SNPs) for osteoporotic fracture in 2,023 cases and 3,740 controls. Association with CBT and forearm BMD was tested for ∼2.5 million SNPs in each cohort separately, and results were meta-analyzed using fixed effect meta-analysis. We identified a missense SNP (Thr>Ile; rs2707466) located in the WNT16 gene (7q31), associated with CBT (effect size of −0.11 standard deviations [SD] per C allele, P = 6.2×10−9). This SNP, as well as another nonsynonymous SNP rs2908004 (Gly>Arg), also had genome-wide significant association with forearm BMD (−0.14 SD per C allele, P = 2.3×10−12, and −0.16 SD per G allele, P = 1.2×10−15, respectively). Four genome-wide significant SNPs arising from BMD meta-analysis were tested for association with forearm fracture. SNP rs7776725 in FAM3C, a gene adjacent to WNT16, was associated with a genome-wide significant increased risk of forearm fracture (OR = 1.33, P = 7.3×10−9), with genome-wide suggestive signals from the two missense variants in WNT16 (rs2908004: OR = 1.22, P = 4.9×10−6 and rs2707466: OR = 1.22, P = 7.2×10−6). We next generated a homozygous mouse with targeted disruption of Wnt16. Female Wnt16−/− mice had 27% (P<0.001) thinner cortical bones at the femur midshaft, and bone strength measures were reduced between 43%–61% (6.5×10−13<P<5.9×10−4) at both femur and tibia, compared with their wild-type littermates. Natural variation in humans and targeted disruption in mice demonstrate that WNT16 is an important determinant of CBT, BMD, bone strength, and risk of fracture.

Genetic Determinants of Serum Testosterone Concentrations in Men

Testosterone concentrations in men are associated with cardiovascular morbidity, osteoporosis, and mortality and are affected by age, smoking, and obesity. Because of serum testosterones high heritability, we performed a meta-analysis of genome-wide association data in 8,938 men from seven cohorts and followed up the genome-wide significant findings in one in silico (n = 871) and two de novo replication cohorts (n = 4,620) to identify genetic loci significantly associated with serum testosterone concentration in men. All these loci were also associated with low serum testosterone concentration defined as <300 ng/dl. Two single-nucleotide polymorphisms at the sex hormone-binding globulin (SHBG) locus (17p13-p12) were identified as independently associated with serum testosterone concentration (rs12150660, p = 1.2×10−41 and rs6258, p = 2.3×10−22). Subjects with ≥3 risk alleles of these variants had 6.5-fold higher risk of having low serum testosterone than subjects with no risk allele. The rs5934505 polymorphism near FAM9B on the X chromosome was also associated with testosterone concentrations (p = 5.6×10−16). The rs6258 polymorphism in exon 4 of SHBG affected SHBGs affinity for binding testosterone and the measured free testosterone fraction (p<0.01). Genetic variants in the SHBG locus and on the X chromosome are associated with a substantial variation in testosterone concentrations and increased risk of low testosterone. rs6258 is the first reported SHBG polymorphism, which affects testosterone binding to SHBG and the free testosterone fraction and could therefore influence the calculation of free testosterone using law-of-mass-action equation.

Genome-wide meta-analysis of observational studies shows common genetic variants associated with macronutrient intake

Background: Macronutrient intake varies substantially between individuals, and there is evidence that this variation is partly accounted for by genetic variants. Objective: The objective of the study was to identify common genetic variants that are associated with macronutrient intake. Design: We performed 2-stage genome-wide association (GWA) meta-analysis of macronutrient intake in populations of European descent. Macronutrients were assessed by using food-frequency questionnaires and analyzed as percentages of total energy consumption from total fat, protein, and carbohydrate. From the discovery GWA (n = 38,360), 35 independent loci associated with macronutrient intake at P < 5 × 10−6 were identified and taken forward to replication in 3 additional cohorts (n = 33,533) from the DietGen Consortium. For one locus, fat mass obesity-associated protein (FTO), cohorts with Illumina MetaboChip genotype data (n = 7724) provided additional replication data. Results: A variant in the chromosome 19 locus (rs838145) was associated with higher carbohydrate (β ± SE: 0.25 ± 0.04%; P = 1.68 × 10−8) and lower fat (β ± SE: −0.21 ± 0.04%; P = 1.57 × 10−9) consumption. A candidate gene in this region, fibroblast growth factor 21 (FGF21), encodes a fibroblast growth factor involved in glucose and lipid metabolism. The variants in this locus were associated with circulating FGF21 protein concentrations (P < 0.05) but not mRNA concentrations in blood or brain. The body mass index (BMI)–increasing allele of the FTO variant (rs1421085) was associated with higher protein intake (β ± SE: 0.10 ± 0.02%; P = 9.96 × 10−10), independent of BMI (after adjustment for BMI, β ± SE: 0.08 ± 0.02%; P = 3.15 × 10−7). Conclusion: Our results indicate that variants in genes involved in nutrient metabolism and obesity are associated with macronutrient consumption in humans. Trials related to this study were registered at clinicaltrials.gov as NCT00005131 (Atherosclerosis Risk in Communities), NCT00005133 (Cardiovascular Health Study), NCT00005136 (Family Heart Study), NCT00005121 (Framingham Heart Study), NCT00083369 (Genetic and Environmental Determinants of Triglycerides), NCT01331512 (InCHIANTI Study), and NCT00005487 (Multi-Ethnic Study of Atherosclerosis).

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.

Genome-wide study for circulating metabolites identifies 62 loci and reveals novel systemic effects of LPA

Genome-wide association studies have identified numerous loci linked with complex diseases, for which the molecular mechanisms remain largely unclear. Comprehensive molecular profiling of circulating metabolites captures highly heritable traits, which can help to uncover metabolic pathophysiology underlying established disease variants. We conduct an extended genome-wide association study of genetic influences on 123 circulating metabolic traits quantified by nuclear magnetic resonance metabolomics from up to 24,925 individuals and identify eight novel loci for amino acids, pyruvate and fatty acids. The LPA locus link with cardiovascular risk exemplifies how detailed metabolic profiling may inform underlying aetiology via extensive associations with very-low-density lipoprotein and triglyceride metabolism. Genetic fine mapping and Mendelian randomization uncover wide-spread causal effects of lipoprotein(a) on overall lipoprotein metabolism and we assess potential pleiotropic consequences of genetically elevated lipoprotein(a) on diverse morbidities via electronic health-care records. Our findings strengthen the argument for safe LPA-targeted intervention to reduce cardiovascular risk.

Mendelian Randomization Studies Do Not Support a Causal Role for Reduced Circulating Adiponectin Levels in Insulin Resistance and Type 2 Diabetes

Adiponectin is strongly inversely associated with insulin resistance and type 2 diabetes, but its causal role remains controversial. We used a Mendelian randomization approach to test the hypothesis that adiponectin causally influences insulin resistance and type 2 diabetes. We used genetic variants at the ADIPOQ gene as instruments to calculate a regression slope between adiponectin levels and metabolic traits (up to 31,000 individuals) and a combination of instrumental variables and summary statistics–based genetic risk scores to test the associations with gold-standard measures of insulin sensitivity (2,969 individuals) and type 2 diabetes (15,960 case subjects and 64,731 control subjects). In conventional regression analyses, a 1-SD decrease in adiponectin levels was correlated with a 0.31-SD (95% CI 0.26–0.35) increase in fasting insulin, a 0.34-SD (0.30–0.38) decrease in insulin sensitivity, and a type 2 diabetes odds ratio (OR) of 1.75 (1.47–2.13). The instrumental variable analysis revealed no evidence of a causal association between genetically lower circulating adiponectin and higher fasting insulin (0.02 SD; 95% CI −0.07 to 0.11; N = 29,771), nominal evidence of a causal relationship with lower insulin sensitivity (−0.20 SD; 95% CI −0.38 to −0.02; N = 1,860), and no evidence of a relationship with type 2 diabetes (OR 0.94; 95% CI 0.75–1.19; N = 2,777 case subjects and 13,011 control subjects). Using the ADIPOQ summary statistics genetic risk scores, we found no evidence of an association between adiponectin-lowering alleles and insulin sensitivity (effect per weighted adiponectin-lowering allele: −0.03 SD; 95% CI −0.07 to 0.01; N = 2,969) or type 2 diabetes (OR per weighted adiponectin-lowering allele: 0.99; 95% CI 0.95–1.04; 15,960 case subjects vs. 64,731 control subjects). These results do not provide any consistent evidence that interventions aimed at increasing adiponectin levels will improve insulin sensitivity or risk of type 2 diabetes.

Genetic determinants of trabecular and cortical volumetric bone mineral densities and bone microstructure

Most previous genetic epidemiology studies within the field of osteoporosis have focused on the genetics of the complex trait areal bone mineral density (aBMD), not being able to differentiate genetic determinants of cortical volumetric BMD (vBMD), trabecular vBMD, and bone microstructural traits. The objective of this study was to separately identify genetic determinants of these bone traits as analysed by peripheral quantitative computed tomography (pQCT). Separate GWA meta-analyses for cortical and trabecular vBMDs were performed. The cortical vBMD GWA meta-analysis (n = 5,878) followed by replication (n = 1,052) identified genetic variants in four separate loci reaching genome-wide significance (RANKL, rs1021188, p = 3.6×10−14; LOC285735, rs271170, p = 2.7×10−12; OPG, rs7839059, p = 1.2×10−10; and ESR1/C6orf97, rs6909279, p = 1.1×10−9). The trabecular vBMD GWA meta-analysis (n = 2,500) followed by replication (n = 1,022) identified one locus reaching genome-wide significance (FMN2/GREM2, rs9287237, p = 1.9×10−9). High-resolution pQCT analyses, giving information about bone microstructure, were available in a subset of the GOOD cohort (n = 729). rs1021188 was significantly associated with cortical porosity while rs9287237 was significantly associated with trabecular bone fraction. The genetic variant in the FMN2/GREM2 locus was associated with fracture risk in the MrOS Sweden cohort (HR per extra T allele 0.75, 95% confidence interval 0.60–0.93) and GREM2 expression in human osteoblasts. In conclusion, five genetic loci associated with trabecular or cortical vBMD were identified. Two of these (FMN2/GREM2 and LOC285735) are novel bone-related loci, while the other three have previously been reported to be associated with aBMD. The genetic variants associated with cortical and trabecular bone parameters differed, underscoring the complexity of the genetics of bone parameters. We propose that a genetic variant in the RANKL locus influences cortical vBMD, at least partly, via effects on cortical porosity, and that a genetic variant in the FMN2/GREM2 locus influences GREM2 expression in osteoblasts and thereby trabecular number and thickness as well as fracture risk.

Detailed metabolic and genetic characterization reveals new associations for 30 known lipid loci

Almost 100 genetic loci are known to affect serum cholesterol and triglyceride levels. For many of these loci, the biological function and causal variants remain unknown. We performed an association analysis of the reported 95 lipid loci against 216 metabolite measures, including 95 measurements on lipids and lipoprotein subclasses, obtained via serum nuclear magnetic resonance metabolomics and four enzymatic lipid traits in 8330 individuals from Finland. The genetic variation in the loci was investigated using a dense set of 440 807 directly genotyped and imputed variants around the previously identified lead single nucleotide polymorphisms (SNPs). For 30 of the 95 loci, we identified new metabolic or genetic associations (P < 5 × 10(-8)). In the majority of the loci, the strongest association was to a more specific metabolite measure than the enzymatic lipids. In four loci, the smallest high-density lipoprotein measures showed effects opposite to the larger ones, and 14 loci had associations beyond the individual lipoprotein measures. In 27 loci, we identified SNPs with a stronger association than the previously reported markers and 12 loci harboured multiple, statistically independent variants. Our data show considerable diversity in association patterns between the loci originally identified through associations with enzymatic lipid measures and reveal association profiles of far greater detail than from routine clinical lipid measures. Additionally, a dense marker set and a homogeneous population allow for detailed characterization of the genetic association signals to a resolution exceeding that achieved so far. Further understanding of the rich variability in genetic effects on metabolites provides insights into the biological processes modifying lipid levels.