Background

The gut microbiome is a complex and metabolically active community that directly influences host phenotypes. In this study, we profile gut microbiota using 16S rRNA gene sequencing in 531 well-phenotyped Finnish men from the Metabolic Syndrome In Men (METSIM) study.

Results

We investigate gut microbiota relationships with a variety of factors that have an impact on the development of metabolic and cardiovascular traits. We identify novel associations between gut microbiota and fasting serum levels of a number of metabolites, including fatty acids, amino acids, lipids, and glucose. In particular, we detect associations with fasting plasma trimethylamine N-oxide (TMAO) levels, a gut microbiota-dependent metabolite associated with coronary artery disease and stroke. We further investigate the gut microbiota composition and microbiota-metabolite relationships in subjects with different body mass index and individuals with normal or altered oral glucose tolerance. Finally, we perform microbiota co-occurrence network analysis, which shows that certain metabolites strongly correlate with microbial community structure and that some of these correlations are specific for the pre-diabetic state.

Conclusions

Our study identifies novel relationships between the composition of the gut microbiota and circulating metabolites and provides a resource for future studies to understand host-gut microbiota relationships.

Genome-wide association studies (GWASs) have identified more than 150 loci associated with blood lipid and cholesterol levels; however, the functional and molecular mechanisms for many associations are unknown. We examined the functional regulatory effects of candidate variants at the GALNT2 locus associated with high-density lipoprotein cholesterol (HDL-C). Fine-mapping and conditional analyses in the METSIM study identified a single locus harboring 25 noncoding variants (r(2) > 0.7 with the lead GWAS variants) strongly associated with total cholesterol in medium-sized HDL (e.g., rs17315646, p = 3.5 × 10(-12)). We used luciferase reporter assays in HepG2 cells to test all 25 variants for allelic differences in regulatory enhancer activity. rs2281721 showed allelic differences in transcriptional activity (75-fold [T] versus 27-fold [C] more than the empty-vector control), as did a separate 780-bp segment containing rs4846913, rs2144300, and rs6143660 (49-fold [AT(-) haplotype] versus 16-fold [CC(+) haplotype] more). Using electrophoretic mobility shift assays, we observed differential CEBPB binding to rs4846913, and we confirmed this binding in a native chromatin context by performing chromatin-immunoprecipitation (ChIP) assays in HepG2 and Huh-7 cell lines of differing genotypes. Additionally, sequence reads in HepG2 DNase-I-hypersensitivity and CEBPB ChIP-seq signals spanning rs4846913 showed significant allelic imbalance. Allelic-expression-imbalance assays performed with RNA from primary human hepatocyte samples and expression-quantitative-trait-locus (eQTL) data in human subcutaneous adipose tissue samples confirmed that alleles associated with increased HDL-C are associated with a modest increase in GALNT2 expression. Together, these data suggest that at least rs4846913 and rs2281721 play key roles in influencing GALNT2 expression at this HDL-C locus.

We investigated the association of the levels of ketone bodies (KBs) with hyperglycemia and with 62 genetic risk variants regulating glucose levels or type 2 diabetes in the population-based Metabolic Syndrome in Men (METSIM) study, including 9,398 Finnish men without diabetes or newly diagnosed type 2 diabetes. Increasing fasting and 2-h plasma glucose levels were associated with elevated levels of acetoacetate (AcAc) and β-hydroxybutyrate (BHB). AcAc and BHB predicted an increase in the glucose area under the curve in an oral glucose tolerance test, and AcAc predicted the conversion to type 2 diabetes in a 5-year follow-up of the METSIM cohort. Impaired insulin secretion, but not insulin resistance, explained these findings. Of the 62 single nucleotide polymorphisms associated with the risk of type 2 diabetes or hyperglycemia, the glucose-increasing C allele of GCKR significantly associated with elevated levels of fasting BHB levels. Adipose tissue mRNA expression levels of genes involved in ketolysis were significantly associated with insulin sensitivity (Matsuda index). In conclusion, high levels of KBs predicted subsequent worsening of hyperglycemia, and a common variant of GCKR was significantly associated with BHB levels.

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.

The prevalence of hypertension in African Americans (AAs) is higher than in other US groups; yet, few have performed genome-wide association studies (GWASs) in AA. Among people of European descent, GWASs have identified genetic variants at 13 loci that are associated with blood pressure. It is unknown if these variants confer susceptibility in people of African ancestry. Here, we examined genome-wide and candidate gene associations with systolic blood pressure (SBP) and diastolic blood pressure (DBP) using the Candidate Gene Association Resource (CARe) consortium consisting of 8591 AAs. Genotypes included genome-wide single-nucleotide polymorphism (SNP) data utilizing the Affymetrix 6.0 array with imputation to 2.5 million HapMap SNPs and candidate gene SNP data utilizing a 50K cardiovascular gene-centric array (ITMAT-Broad-CARe [IBC] array). For Affymetrix data, the strongest signal for DBP was rs10474346 (P= 3.6 × 10(-8)) located near GPR98 and ARRDC3. For SBP, the strongest signal was rs2258119 in C21orf91 (P= 4.7 × 10(-8)). The top IBC association for SBP was rs2012318 (P= 6.4 × 10(-6)) near SLC25A42 and for DBP was rs2523586 (P= 1.3 × 10(-6)) near HLA-B. None of the top variants replicated in additional AA (n = 11 882) or European-American (n = 69 899) cohorts. We replicated previously reported European-American blood pressure SNPs in our AA samples (SH2B3, P= 0.009; TBX3-TBX5, P= 0.03; and CSK-ULK3, P= 0.0004). These genetic loci represent the best evidence of genetic influences on SBP and DBP in AAs to date. More broadly, this work supports that notion that blood pressure among AAs is a trait with genetic underpinnings but also with significant complexity.

Plasminogen activator inhibitor type 1 (PAI-1) plays an essential role in the fibrinolysis system and thrombosis. Population studies have reported that blood PAI-1 levels are associated with increased risk of coronary heart disease (CHD). However, it is unclear whether the association reflects a causal influence of PAI-1 on CHD risk. To evaluate the association between PAI-1 and CHD, we applied a 3-step strategy. First, we investigated the observational association between PAI-1 and CHD incidence using a systematic review based on a literature search for PAI-1 and CHD studies. Second, we explored the causal association between PAI-1 and CHD using a Mendelian randomization approach using summary statistics from large genome-wide association studies. Finally, we explored the causal effect of PAI-1 on cardiovascular risk factors including metabolic and subclinical atherosclerosis measures. In the systematic meta-analysis, the highest quantile of blood PAI-1 level was associated with higher CHD risk comparing with the lowest quantile (odds ratio=2.17; 95% CI: 1.53, 3.07) in an age- and sex-adjusted model. The effect size was reduced in studies using a multivariable-adjusted model (odds ratio=1.46; 95% CI: 1.13, 1.88). The Mendelian randomization analyses suggested a causal effect of increased PAI-1 level on CHD risk (odds ratio=1.22 per unit increase of log-transformed PAI-1; 95% CI: 1.01, 1.47). In addition, we also detected a causal effect of PAI-1 on elevating blood glucose and high-density lipoprotein cholesterol. Our study indicates a causal effect of elevated PAI-1 level on CHD risk, which may be mediated by glucose dysfunction.