ABSTRACT
Transcriptomics data have been integrated with genome-wide association studies (GWASs) to help understand disease/trait molecular mechanisms. The utility of metabolomics, integrated with transcriptomics and disease GWASs, to understand molecular mechanisms for metabolite levels or diseases has not been thoroughly evaluated. We performed probabilistic transcriptome-wide association and locus-level colocalization analyses to integrate transcriptomics results for 49 tissues in 706 individuals from the GTEx project, metabolomics results for 1,391 plasma metabolites in 6,136 Finnish men from the METSIM study, and GWAS results for 2,861 disease traits in 260,405 Finnish individuals from the FinnGen study. We found that genetic variants that regulate metabolite levels were more likely to influence gene expression and disease risk compared to the ones that do not. Integrating transcriptomics with metabolomics results prioritized 397 genes for 521 metabolites, including 496 previously identified gene-metabolite pairs with strong functional connections and suggested 33.3% of such gene-metabolite pairs shared the same causal variants with genetic associations of gene expression. Integrating transcriptomics and metabolomics individually with FinnGen GWAS results identified 1,597 genes for 790 disease traits. Integrating transcriptomics and metabolomics jointly with FinnGen GWAS results helped pinpoint metabolic pathways from genes to diseases. We identified putative causal effects of UGT1A1/UGT1A4 expression on gallbladder disorders through regulating plasma (E,E)-bilirubin levels, of SLC22A5 expression on nasal polyps and plasma carnitine levels through distinct pathways, and of LIPC expression on age-related macular degeneration through glycerophospholipid metabolic pathways. Our study highlights the power of integrating multiple sets of molecular traits and GWAS results to deepen understanding of disease pathophysiology.
Subject(s)
Genome-Wide Association Study , Transcriptome , Bilirubin , Carnitine , Glycerophospholipids , Humans , Male , Metabolomics , Quantitative Trait Loci/genetics , Solute Carrier Family 22 Member 5/genetics , Transcriptome/geneticsABSTRACT
CONTEXT: Lifestyle intervention prevents or delays type 2 diabetes (T2D) in subjects at a high risk of T2D. However, it is not known whether genetic variants modify the effect on incident T2D during lifestyle intervention. OBJECTIVE: To investigate whether a low or high genetic risk has effects on incident T2D in a group-based lifestyle intervention study. METHODS: The T2D-GENE trial involved 973 men from the Metabolic Syndrome in Men (METSIM) cohort, aged 50-75 years, body mass index ≥25 kg/m2, fasting plasma glucose 5.6-6.9 mmol/l, and HbA1c <48 mmol/mol and either a low or high genetic risk score for T2D. There were two intervention groups, a low (n=315) and high genetic risk for T2D (n=313). They were provided with a 3-year group-based intervention with access to a web portal focused on healthy diet and physical activity. There were also corresponding population-based control groups at low (n=196), and high genetic risk for T2D (n=149) who had two laboratory visits (0 and 3 years) and general health advice as a part of their METSIM cohort protocol. The primary outcome was incident T2D, and a secondary outcome glycemia. RESULTS: The intervention significantly lowered the risk of T2D among the participants with a high genetic risk for T2D (HR 0.30, 95% CI 0.16-0.56, p<0.001) whereas in the low genetic risk group the effect was not significant (HR 0.69, 95% CI 0.36-1.32, p=0.262). The intervention effect was not significantly different between the high and low genetic risk groups (p=0.135). The intervention significantly ameliorated the worsening of glycemia and decreased weight both in the low and high genetic risk groups. CONCLUSIONS: Our results showed that individuals with a high genetic risk for T2D benefited from a low-cost group-based intervention focusing on healthy diet and physical activity. Therefore, all individuals at risk of T2D should be encouraged to make lifestyle changes regardless of genetic risk.
ABSTRACT
Both genetic and non-genetic factors are important in the pathophysiology of non-alcoholic fatty liver disease (NAFLD). The aim of our study was to identify novel metabolites and pathways associated with NAFLD by including both genetic and non-genetic factors in statistical analyses. We genotyped six genetic variants in the PNPLA3, TM6SF2, MBOAT7, GCKR, PPP1R3B, and HSD17B13 genes reported to be associated with NAFLD. Non-targeted metabolomic profiling was performed from plasma samples. We applied a previously validated fatty liver index to identify participants with NAFLD. First, we associated the six genetic variants with 1098 metabolites in 2 339 men without NAFLD to determine the effects of the genetic variants on metabolites, and then in 2 535 men with NAFLD to determine the joint effects of genetic variants and non-genetic factors on metabolites. We identified several novel metabolites and metabolic pathways, especially for PNPLA3, GCKR, and PPP1R38 variants relevant to the pathophysiology of NAFLD. Importantly, we showed that each genetic variant for NAFLD had a specific metabolite signature. The plasma metabolite signature was unique for each genetic variant, suggesting that several metabolites and different pathways are involved in the risk of NAFLD. The FLI index reliably identifies metabolites for NAFLD in large population-based studies.
ABSTRACT
CONTEXT: Diabetic retinopathy (DR) is a specific microvascular complication in patients with diabetes and the leading cause of blindness. Recent advances in omics, especially metabolomics, offer the possibility identifying novel potential biomarkers for DR. OBJECTIVE: The aim was to identify metabolites associated with DR. METHODS: We performed a 12-year follow-up study including 1349 participants with type 2 diabetes (1021 without DR, 328 with DR) selected from the METSIM cohort. Individuals who had retinopathy before the baseline study were excluded (n = 63). The diagnosis of retinopathy was based on fundus photography examination. We performed nontargeted metabolomics profiling to identify metabolites. RESULTS: We found 17 metabolites significantly associated with incident DR after adjustment for confounding factors. Among amino acids, N-lactoyl isoleucine, N-lactoyl valine, N-lactoyl tyrosine, N-lactoyl phenylalanine, N-(2-furoyl) glycine, and 5-hydroxylysine were associated with an increased risk of DR, and citrulline with a decreased risk of DR. Among the fatty acids N,N,N-trimethyl-5-aminovalerate was associated with an increased risk of DR, and myristoleate (14:1n5), palmitoleate (16:1n7), and 5-dodecenoate (12:1n7) with a decreased risk of DR. Sphingomyelin (d18:2/24:2), a sphingolipid, was significantly associated with a decreased risk of DR. Carboxylic acid maleate and organic compounds 3-hydroxypyridine sulfate, 4-vinylphenol sulfate, 4-ethylcatechol sulfate, and dimethyl sulfone were significantly associated with an increased risk of DR. CONCLUSION: Our study is the first large population-based longitudinal study to identify metabolites for DR. We found multiple metabolites associated with an increased and decreased risk for DR from several different metabolic pathways.
Subject(s)
Diabetes Mellitus, Type 2 , Diabetic Retinopathy , Humans , Diabetes Mellitus, Type 2/complications , Diabetes Mellitus, Type 2/epidemiology , Diabetes Mellitus, Type 2/diagnosis , Diabetic Retinopathy/diagnosis , Diabetic Retinopathy/epidemiology , Diabetic Retinopathy/etiology , Follow-Up Studies , Longitudinal Studies , Risk Factors , SulfatesABSTRACT
Metabolites are small molecules that are useful for estimating disease risk and elucidating disease biology. Nevertheless, their causal effects on human diseases have not been evaluated comprehensively. We performed two-sample Mendelian randomization to systematically infer the causal effects of 1,099 plasma metabolites measured in 6,136 Finnish men from the METSIM study on risk of 2,099 binary disease endpoints measured in 309,154 Finnish individuals from FinnGen. We identified evidence for 282 causal effects of 70 metabolites on 183 disease endpoints (FDR<1%). We found 25 metabolites with potential causal effects across multiple disease domains, including ascorbic acid 2-sulfate affecting 26 disease endpoints in 12 disease domains. Our study suggests that N-acetyl-2-aminooctanoate and glycocholenate sulfate affect risk of atrial fibrillation through two distinct metabolic pathways and that N-methylpipecolate may mediate the causal effect of N6, N6-dimethyllysine on anxious personality disorder. This study highlights the broad causal impact of plasma metabolites and widespread metabolic connections across diseases.
ABSTRACT
Understanding differences in adipose gene expression between individuals with different levels of clinical traits may reveal the genes and mechanisms leading to cardiometabolic diseases. However, adipose is a heterogeneous tissue. To account for cell-type heterogeneity, we estimated cell-type proportions in 859 subcutaneous adipose tissue samples with bulk RNA sequencing (RNA-seq) using a reference single-nuclear RNA-seq data set. Cell-type proportions were associated with cardiometabolic traits; for example, higher macrophage and adipocyte proportions were associated with higher and lower BMI, respectively. We evaluated cell-type proportions and BMI as covariates in tests of association between >25,000 gene expression levels and 22 cardiometabolic traits. For >95% of genes, the optimal, or best-fit, models included BMI as a covariate, and for 79% of associations, the optimal models also included cell type. After adjusting for the optimal covariates, we identified 2,664 significant associations (P ≤ 2e-6) for 1,252 genes and 14 traits. Among genes proposed to affect cardiometabolic traits based on colocalized genome-wide association study and adipose expression quantitative trait locus signals, 25 showed a corresponding association between trait and gene expression levels. Overall, these results suggest the importance of modeling cell-type proportion when identifying gene expression associations with cardiometabolic traits.
Subject(s)
Cardiovascular Diseases , Genome-Wide Association Study , Humans , Body Mass Index , Obesity/genetics , Gene Expression , Cardiovascular Diseases/geneticsABSTRACT
Complete characterization of the genetic effects on gene expression is needed to elucidate tissue biology and the etiology of complex traits. Here, we analyzed 2,344 subcutaneous adipose tissue samples and identified 34K conditionally distinct expression quantitative trait locus (eQTL) signals in 18K genes. Over half of eQTL genes exhibited at least two eQTL signals. Compared to primary signals, non-primary signals had lower effect sizes, lower minor allele frequencies, and less promoter enrichment; they corresponded to genes with higher heritability and higher tolerance for loss of function. Colocalization of eQTL with conditionally distinct genome-wide association study signals for 28 cardiometabolic traits identified 3,605 eQTL signals for 1,861 genes. Inclusion of non-primary eQTL signals increased colocalized signals by 46%. Among 30 genes with ≥2 pairs of colocalized signals, 21 showed a mediating gene dosage effect on the trait. Thus, expanded eQTL identification reveals more mechanisms underlying complex traits and improves understanding of the complexity of gene expression regulation.
ABSTRACT
Organic-anion-transporting polypeptide 1B1 (OATP1B1), encoded by the solute carrier organic anion transporter family member 1B1 gene (SLCO1B1), is highly expressed in the liver and transports several endogenous metabolites into the liver, including statins. Previous studies have not investigated the association of SLCO1B1 rs4149056 variant with the risk of type 2 diabetes (T2D) or determined the metabolite signature of the C allele of SLCO1B1 rs4149056 (SLCO1B1 rs4149056-C allele) in a large randomly selected population. SLCO1B1 rs4149056-C inhibits OATP1B1 transporter and is associated with increased levels of blood simvastatin concentrations. Our study is to first to show that SLCO1B1 rs4149056 variant is not significantly associated with the risk of T2D, suggesting that simvastatin has a direct effect on the risk of T2D. Additionally, we investigated the effects of SLCO1B1 rs4149056-C on plasma metabolite concentrations in 1373 participants on simvastatin treatment and in 1368 age- and body-mass index (BMI)-matched participants without any statin treatment. We found 31 novel metabolites significantly associated with SLCO1B1 rs4149056-C in the participants on simvastatin treatment and in the participants without statin treatment. Simvastatin decreased concentrations of dicarboxylic acids, such as docosadioate and dodecanedioate, that may increase beta- and peroxisomal oxidation and increased the turnover of cholesterol into bile acids, resulting in a decrease in steroidogenesis due to limited availability of cholesterol for steroid synthesis. Our findings suggest that simvastatin exerts its effects on the lowering of low-density lipoprotein (LDL) cholesterol concentrations through several distinct pathways in the carriers of SLCO1B1 rs4149056-C, including dicarboxylic acids, bile acids, steroids, and glycerophospholipids.