Integrating transcriptomics, metabolomics, and GWAS helps reveal molecular mechanisms for metabolite levels and disease risk

Recent advances in omics and the integration of multi-omics in osteoarthritis research

Osteoarthritis (OA) is a complex disorder driven by the combination of environmental and genetic factors. Given its high global prevalence and heterogeneity, developing effective and personalized treatment methods is crucial. This requires identifying new disease mechanisms, drug targets, and biomarkers. Various omics approaches have been applied to identify OA-related genes, pathways, and biomarkers, including genomics, epigenomics, transcriptomics, proteomics, and metabolomics. These omics studies have generated vast datasets that are shaping the field of OA research. The emergence of high-resolution methodologies, such as single-cell and spatial omics techniques, further enhances our ability to dissect molecular complexities within the OA microenvironment. By integrating these multi-layered datasets, researchers can uncover central signaling hubs and disease mechanisms, ultimately facilitating the development of targeted therapies and precision medicine approaches for OA treatment.

Recent developments in population biobanks and the genetic architecture of complex disease

Population biobanks have radically transformed our understanding of complex disease genetics. Recent technological advances and the inclusion of diverse populations have accelerated the discovery and interpretation of variant associations. For instance, population-scale whole-genome sequencing now allows deep exploration of rare and structural variant associations, while multi-omics approaches integrating genome-wide association studies with proteomics, metabolomics, and advanced statistical methods like Mendelian randomization provide nuanced insights into genetic disease mechanisms. Additionally, cross-biobank collaborations and meta-analyses have been particularly impactful, dramatically increasing the statistical power for discovery. These efforts have identified novel genetic associations across numerous complex diseases, with significant contributions from non-European populations. However, data integration complexities, privacy concerns, and methodological limitations continue to constrain research. Here we review how recent advances have contributed to genetic discovery.

Advances in Metabolomics: A Comprehensive Review of Type 2 Diabetes and Cardiovascular Disease Interactions

Type 2 diabetes (T2D) and cardiovascular diseases (CVDs) are major public health challenges worldwide. Metabolomics, the exhaustive assessment of metabolites in biological systems, offers important insights regarding the metabolic disturbances related to these disorders. Recent advances toward the integration of metabolomics into clinical practice to facilitate the discovery of novel biomarkers that can improve the diagnosis, prognosis, and treatment of T2D and CVDs are discussed in this review. Metabolomics offers the potential to characterize the key metabolic alterations associated with disease pathophysiology and treatment. T2D is a heterogeneous disease that develops through diverse pathophysiological processes and molecular mechanisms; therefore, the disease-causing pathways of T2D are not completely understood. Recent studies have identified several robust clusters of T2D variants representing biologically meaningful, distinct pathways, such as the beta cell and proinsulin cluster related to pancreatic insulin secretion, obesity, lipodystrophy, the liver/lipid cluster, glycemia, and blood pressure, and metabolic syndrome clusters representing different pathways causing insulin resistance. Regarding CVDs, recent studies have allowed the metabolomic profile to delineate pathways that contribute to atherosclerosis and heart failure, as well as to the development of targeted therapy. This review also covers the role of metabolomics in integrated metabolic genomics and other omics platforms to better understand disease mechanisms, along with the transition toward precision medicine. This review further investigates the use of metabolomics in multi-metabolite modeling to enhance risk prediction models for predicting the first occurrence of major adverse cardiovascular events among individuals with T2D, highlighting the value of such approaches in optimizing the preventive and therapeutic models used in clinical practice.

Examining the link between 179 lipid species and 7 diseases using genetic predictors

BACKGROUND

Genome-wide association studies of lipid species have identified several loci shared with various diseases, however, the relationship between lipid species and disease risk remains poorly understood. Here we investigated whether the plasma levels of lipid species are causally linked to disease risk.

METHODS

We built genetic predictors of 179 lipid species, measured in 7174 Finnish individuals, by utilising either 11 high-impact genomic loci or genome-wide polygenic scores (PGS). We assessed the impact of the lipid species on seven diseases by performing disease association across FinnGen (n = 500,348), UK Biobank (n = 420,531), and Generation Scotland (n = 20,032). We performed univariable Mendelian randomisation (MR) and multivariable MR (MVMR) analyses to examine whether lipid species impact disease risk independently of standard lipids.

FINDINGS

PGS explained >4% of the variance for 34 lipid species but variants outside the high-impact loci had only a marginal contribution. Variants within the high-impact loci showed association with all seven diseases. MVMR supported a causal role of ApoB in ischaemic heart disease after accounting for lipid species. Phosphatidylethanolamine-increasing LIPC variants seemed to lower age-related macular degeneration risk independently of HDL-cholesterol. MVMR suggested a protective effect of four lipid species containing arachidonic acid on cholelithiasis risk independently of Total Cholesterol.

INTERPRETATION

Our study demonstrates how genetic predictors of lipid species can be utilised to gain insights into disease risk. We report potential links between lipid species and age-related macular degeneration and cholelithiasis risk, which can be explored for their utility in disease risk prediction and therapy.

FUNDING

The funders had no role in the study design, data analyses, interpretation, or writing of this article.

Transcriptomic and Metabolomic Analyses in Monozygotic and Dizygotic Twins

Monozygotic (MZ) and dizygotic (DZ) twins are studied to understand genetic and environmental influences on complex traits, however the mechanisms behind twinning are not completely understood. (Epi)genomic studies identified SNPs associated with DZ twinning and DNA methylation sites with MZ twinning. To find molecular biomarkers of twinning, we compared transcriptomics and metabolomics data from MZ and DZ twins. We analyzed 42,663 RNA transcripts in 1453 MZ twins and 1294 DZ twins from the Netherlands Twin Register (NTR), followed by sex-stratified analyses. The top 5% transcripts with lowest p-values were analyzed for replication in 217 MZ and 158 DZ twins from the older Finnish Twin cohort (FTC). In the NTR, one transcript (PURG) was significantly differentially expressed between MZ and DZ twins; but this did not replicate in FTC. Pathway analyses highlighted the WNT-pathway, previously associated with MZ twinning, and the TGF-B and SMAD pathway, previously associated with DZ twinning. Meta-analysis of 169 serum metabolites in 2797 MZ and 2040 DZ twins from the NTR, FTC and FinnTwin12, showed no metabolomic differences. Overall, we did not find replicable transcript-level expression differences in blood between MZ and DZ twins, but highlighted the TGF-B/SMAD pathway as a potential transcriptional biomarker for DZ twinning.

Assessing the potential causal effects of 1099 plasma metabolites on 2099 binary disease endpoints

Metabolites are small molecules that are useful for estimating disease risk and elucidating disease biology. Here, we perform two-sample Mendelian randomization to systematically infer the potential causal effects of 1099 plasma metabolites measured in 6136 Finnish men from the METSIM study on risk of 2099 binary disease endpoints measured in 309,154 Finnish individuals from FinnGen. We find evidence for 282 putative causal effects of 70 metabolites on 183 disease endpoints. We also identify 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 putative causal effect of N6,N6-dimethyllysine on anxious personality disorder.

Multi-INTACT: integrative analysis of the genome, transcriptome, and proteome identifies causal mechanisms of complex traits

We present multi-integration of transcriptome-wide association studies and colocalization (Multi-INTACT), an algorithm that models multiple "gene products" (e.g., encoded RNA transcript and protein levels) to implicate causal genes and relevant gene products. In simulations, Multi-INTACT achieves higher power than existing methods, maintains calibrated false discovery rates, and detects the true causal gene product(s). We apply Multi-INTACT to GWAS on 1408 metabolites, integrating the GTEx expression and UK Biobank protein QTL datasets. Multi-INTACT infers 52 to 109% more metabolite causal genes than protein-alone or expression-alone analyses and indicates both gene products are relevant for most gene nominations.

Mapping integral cell-type-specific interferon-induced gene regulatory networks (GRNs) involved in systemic lupus erythematosus using systems and computational analysis

Systemic lupus erythematosus (SLE) is a systemic autoimmune disorder characterized by the production of autoantibodies, resulting in inflammation and organ damage. Although extensive research has been conducted on SLE pathogenesis, a comprehensive understanding of its molecular landscape in different cell types has not been achieved. This study uncovers the molecular mechanisms of the disease by thoroughly examining gene regulatory networks within neutrophils, dendritic cells, T cells, and B cells. Firstly, we identified genes and ncRNAs with differential expression in SLE patients compared to controls for different cell types. Furthermore, the derived differentially expressed genes were curated based on immune functions using functional enrichment analysis to create a protein-protein interaction network. Topological network analysis of the formed genes revealed key hub genes associated with each of the cell types. To understand the regulatory mechanism through which these hub genes function in the diseased state, their associations with transcription factors, and non-coding RNAs in different immune cell types were investigated through correlation analysis and regression models. Finally, by integrating these findings, distinct gene regulatory networks were constructed, which provide a novel perspective on the molecular, cellular, and immunological landscapes of SLE. Importantly, we reveal the crucial role of IRF3, IRF7, and STAT1 in neutrophils, dendritic cells, and T cells, where their aberrant upregulation in disease states might enhance the production of type I IFN. Furthermore, we found MYB to be a crucial regulator that might activate T cells toward autoimmune responses in SLE. Similarly, in B-cell lymphocytes, we found FOXO1 to be a key player in autophagy and chemokine regulation. These findings were also validated using single-cell RNASeq analysis using an independent dataset. Genotype variations of these genes were also explored using the GWAS central database to ensure their targetability. These findings indicate that IRF3, IRF7, STAT1, MYB, and FOXO1 are promising targets for therapeutic interventions for SLE.

Rapid and accurate multi-phenotype imputation for millions of individuals

Deep phenotyping can enhance the power of genetic analysis, including genome-wide association studies (GWAS), but the occurrence of missing phenotypes compromises the potential of such resources. Although many phenotypic imputation methods have been developed, the accurate imputation of millions of individuals remains challenging. In the present study, we have developed a multi-phenotype imputation method based on mixed fast random forest (PIXANT) by leveraging efficient machine learning (ML)-based algorithms. We demonstrate by extensive simulations that PIXANT is reliable, robust and highly resource-efficient. We then apply PIXANT to the UKB data of 277,301 unrelated White British citizens and 425 traits, and GWAS is subsequently performed on the imputed phenotypes, 18.4% more GWAS loci are identified than before imputation (8710 vs 7355). The increased statistical power of GWAS identified some additional candidate genes affecting heart rate, such as RNF220, SCN10A, and RGS6, suggesting that the use of imputed phenotype data from a large cohort may lead to the discovery of additional candidate genes for complex traits.

Refinement of a Published Gene-Physical Activity Interaction Impacting HDL-Cholesterol: Role of Sex and Lipoprotein Subfractions

Large-scale gene-environment interaction (GxE) discovery efforts often involve analytical compromises for the sake of data harmonization and statistical power. Refinement of exposures, covariates, outcomes, and population subsets may be helpful to establish often-elusive replication and evaluate potential clinical utility. Here, we used additional datasets, an expanded set of statistical models, and interrogation of lipoprotein metabolism via nuclear magnetic resonance (NMR)-based lipoprotein subfractions to refine a previously discovered GxE modifying the relationship between physical activity (PA) and HDL-cholesterol (HDL-C). We explored this GxE in the Women's Genome Health Study (WGHS; N = 23,294; the strongest cohort-specific signal in the original meta-analysis), the UK Biobank (UKB; N = 281,380), and the Multi-Ethnic Study of Atherosclerosis (MESA; N = 4587), using self-reported PA (MET-min/wk) and genotypes at rs295849 (nearest gene: LHX1). As originally reported, minor allele carriers of rs295849 in WGHS had a stronger positive association between PA and HDL-C (pint = 0.002). When testing available NMR metabolites to refine the HDL-C outcome, we found a stronger interaction effect on medium-sized HDL particle concentrations (M-HDL-P; pint = 1.0 × 10-4) than HDL-C. Meta-regression revealed a systematically larger interaction effect in cohorts from the original meta-analysis with a greater fraction of women (p = 0.018). In the UKB, GxE effects were stronger in women and using M-HDL-P as the outcome. In MESA, the primary interaction for HDL-C showed nominal significance (pint = 0.013), but without clear sex differences and with a greater magnitude for large HDL-P. Our work provides additional insights into a known gene-PA interaction while illustrating the importance of phenotype and model refinement toward understanding and replicating GxEs.

Technological advancements in functional interpretation of genome-wide association studies (GWAS) findings: bridging the gap to clinical translation

Genome-wide association studies (GWAS) significantly advanced our understanding of the genetic underpinnings of diseases. However, challenges persist, particularly in interpreting non-coding variants in linkage disequilibrium that affect genes in disease-relevant cells. Addressing key obstacles-identifying causal variants, uncovering target genes, and understanding their network impact-is crucial. This graphical review navigates advanced techniques to fully leverage GWAS for future therapeutic breakthroughs.

Genetic architecture of cerebrospinal fluid and brain metabolite levels and the genetic colocalization of metabolites with human traits

Brain metabolism perturbation can contribute to traits and diseases. We conducted a genome-wide association study for cerebrospinal fluid (CSF) and brain metabolite levels, identifying 205 independent associations (47.3% new signals, containing 11 new loci) for 139 CSF metabolites, and 32 independent associations (43.8% new signals, containing 4 new loci) for 31 brain metabolites. Of these, 96.9% (CSF) and 71.4% (brain) of the new signals belonged to previously analyzed metabolites in blood or urine. We integrated the metabolite quantitative trait loci (MQTLs) with 23 neurological, psychiatric and common human traits and diseases through colocalization to identify metabolites and biological processes implicated in these phenotypes. Combining CSF and brain, we identified 71 metabolite-trait associations, such as glycerophosphocholines with Alzheimer's disease, O-sulfo-L-tyrosine with Parkinson's disease, glycine, xanthine with waist-to-hip ratio and ergothioneine with inflammatory bowel disease. Our study expanded the knowledge of MQTLs in the central nervous system, providing insights into human traits.

Leveraging large-scale multi-omics evidences to identify therapeutic targets from genome-wide association studies

BACKGROUND

Therapeutic targets supported by genetic evidence from genome-wide association studies (GWAS) show higher probability of success in clinical trials. GWAS is a powerful approach to identify links between genetic variants and phenotypic variation; however, identifying the genes driving associations identified in GWAS remains challenging. Integration of molecular quantitative trait loci (molQTL) such as expression QTL (eQTL) using mendelian randomization (MR) and colocalization analyses can help with the identification of causal genes. Careful interpretation remains warranted because eQTL can affect the expression of multiple genes within the same locus.

METHODS

We used a combination of genomic features that include variant annotation, activity-by-contact maps, MR, and colocalization with molQTL to prioritize causal genes across 4,611 disease GWAS and meta-analyses from biobank studies, namely FinnGen, Estonian Biobank and UK Biobank.

RESULTS

Genes identified using this approach are enriched for gold standard causal genes and capture known biological links between disease genetics and biology. In addition, we find that eQTL colocalizing with GWAS are statistically enriched for corresponding disease-relevant tissues. We show that predicted directionality from MR is generally consistent with matched drug mechanism of actions (> 85% for approved drugs). Compared to the nearest gene mapping method, genes supported by multi-omics evidences displayed higher enrichment in approved therapeutic targets (risk ratio 1.75 vs. 2.58 for genes with the highest level of support). Finally, using this approach, we detected anassociation between the IL6 receptor signal transduction gene IL6ST and polymyalgia rheumatica, an indication for which sarilumab, a monoclonal antibody against IL-6, has been recently approved.

CONCLUSIONS

Combining variant annotation, activity-by-contact maps, and molQTL increases performance to identify causal genes, while informing on directionality which can be translated to successful target identification and drug development.

SNP-Based and Kmer-Based eQTL Analysis Using Transcriptome Data

Traditional expression quantitative trait locus (eQTL) mapping associates single nucleotide polymorphisms (SNPs) with gene expression, where the SNPs are derived from large-scale whole-genome sequencing (WGS) data or transcriptome data. While WGS provides a high SNP density, it also incurs substantial sequencing costs. In contrast, RNA-seq data, which are more accessible and less expensive, can simultaneously yield gene expressions and SNPs. Thus, eQTL analysis based on RNA-seq offers significant potential applications. Two primary strategies were employed for eQTL in this study. The first involved analyzing expression levels in relation to variant sites detected between populations from RNA-seq data. The second approach utilized kmers, which are sequences of length k derived from RNA-seq reads, to represent variant sites and associated these kmer genotypes with gene expression. We discovered 87 significant association signals involving eGene on the basis of the SNP-based eQTL analysis. These genes include DYNLT1, NMNAT1, and MRLC2, which are closely related to neurological functions such as motor coordination and homeostasis, play a role in cellular energy metabolism, and function in regulating calcium-dependent signaling in muscle contraction, respectively. This study compared the results obtained from eQTL mapping using RNA-seq identified SNPs and gene expression with those derived from kmers. We found that the vast majority (23/30) of the association signals overlapping the two methods could be verified by haplotype block analysis. This comparison elucidates the strengths and limitations of each method, providing insights into their relative efficacy for eQTL identification.

Transcriptomics integrated with metabolomics reveals partial molecular mechanisms of nutritional risk and neurodevelopment in children with congenital heart disease

Purpose

To explore molecular mechanisms affecting nutritional risk and neurodevelopment in children with congenital heart disease (CHD) by combining transcriptome and metabolome analysis.

Methods

A total of 26 blood and serum samples from 3 groups of children with CHD low nutritional risk combined with normal neurodevelopment (group A), low nutritional risk combined with neurodevelopmental disorders (group B) and high nutritional risk combined with normal neurodevelopment (group C) were analyzed by transcriptome and metabolomics to search for differentially expressed genes (DEGs) and metabolites (DEMs). Functional analysis was conducted for DEGs and DEMs. Further, the joint pathway analysis and correlation analysis of DEGs and DEMs were performed.

Results

A total of 362 and 1,351 DEGs were detected in group B and C compared to A, respectively. A total of 6 and 7 DEMs were detected in group B and C compared to A in positive mode, respectively. There were 39 and 31 DEMs in group B and C compared to A in negative mode. Transcriptomic analysis indicated that neurodevelopment may be regulated by some genes such as NSUN7, SLC6A8, CXCL1 and LCN8, nutritional risk may be regulated by SLC1A3 and LCN8. Metabolome analysis and joint pathway analysis showed that tryptophan metabolism, linoleic and metabolism and glycerophospholipid metabolism may be related to neurodevelopment, and glycerophospholipid metabolism pathway may be related to nutritional risk.

Conclusion

By integrating transcriptome and metabolome analyses, this study revealed key genes and metabolites associated with nutritional risk and neurodevelopment in children with CHD, as well as significantly altered pathways. It has important clinical translational significance.

Metabolic gene function discovery platform GeneMAP identifies SLC25A48 as necessary for mitochondrial choline import

Organisms maintain metabolic homeostasis through the combined functions of small-molecule transporters and enzymes. While many metabolic components have been well established, a substantial number remains without identified physiological substrates. To bridge this gap, we have leveraged large-scale plasma metabolome genome-wide association studies (GWAS) to develop a multiomic Gene-Metabolite Association Prediction (GeneMAP) discovery platform. GeneMAP can generate accurate predictions and even pinpoint genes that are distant from the variants implicated by GWAS. In particular, our analysis identified solute carrier family 25 member 48 (SLC25A48) as a genetic determinant of plasma choline levels. Mechanistically, SLC25A48 loss strongly impairs mitochondrial choline import and synthesis of its downstream metabolite betaine. Integrative rare variant and polygenic score analyses in UK Biobank provide strong evidence that the SLC25A48 causal effects on human disease may in part be mediated by the effects of choline. Altogether, our study provides a discovery platform for metabolic gene function and proposes SLC25A48 as a mitochondrial choline transporter.

Refinement of a published gene-physical activity interaction impacting HDL-cholesterol: role of sex and lipoprotein subfractions

Large-scale gene-environment interaction (GxE) discovery efforts often involve compromises in the definition of outcomes and choice of covariates for the sake of data harmonization and statistical power. Consequently, refinement of exposures, covariates, outcomes, and population subsets may be helpful to establish often-elusive replication and evaluate potential clinical utility. Here, we used additional datasets, an expanded set of statistical models, and interrogation of lipoprotein metabolism via nuclear magnetic resonance (NMR)-based lipoprotein subfractions to refine a previously discovered GxE modifying the relationship between physical activity (PA) and HDL-cholesterol (HDL-C). This GxE was originally identified by Kilpeläinen et al., with the strongest cohort-specific signal coming from the Women's Genome Health Study (WGHS). We thus explored this GxE further in the WGHS (N = 23,294), with follow-up in the UK Biobank (UKB; N = 281,380), and the Multi-Ethnic Study of Atherosclerosis (MESA; N = 4,587). Self-reported PA (MET-hrs/wk), genotypes at rs295849 (nearest gene: LHX1), and NMR metabolomics data were available in all three cohorts. As originally reported, minor allele carriers of rs295849 in WGHS had a stronger positive association between PA and HDL-C (pint = 0.002). When testing a range of NMR metabolites (primarily lipoprotein and lipid subfractions) to refine the HDL-C outcome, we found a stronger interaction effect on medium-sized HDL particle concentrations (M-HDL-P; pint = 1.0×10-4) than HDL-C. Meta-regression revealed a systematically larger interaction effect in cohorts from the original meta-analysis with a greater fraction of women (p = 0.018). In the UKB, GxE effects were stronger both in women and using M-HDL-P as the outcome. In MESA, the primary interaction for HDL-C showed nominal significance (pint = 0.013), but without clear differences by sex and with a greater magnitude using large, rather than medium, HDL-P as an outcome. Towards reconciling these observations, further exploration leveraging NMR platform-specific HDL subfraction diameter annotations revealed modest agreement across all cohorts in the interaction affecting medium-to-large particles. Taken together, our work provides additional insights into a specific known gene-PA interaction while illustrating the importance of phenotype and model refinement towards understanding and replicating GxEs.

Gene-metabolite annotation with shortest reactional distance enhances metabolite genome-wide association studies results

Metabolite genome-wide association studies (mGWAS) have advanced our understanding of the genetic control of metabolite levels. However, interpreting these associations remains challenging due to a lack of tools to annotate gene-metabolite pairs beyond the use of conservative statistical significance threshold. Here, we introduce the shortest reactional distance (SRD) metric, drawing from the comprehensive KEGG database, to enhance the biological interpretation of mGWAS results. We applied this approach to three independent mGWAS, including a case study on sickle cell disease patients. Our analysis reveals an enrichment of small SRD values in reported mGWAS pairs, with SRD values significantly correlating with mGWAS p values, even beyond the standard conservative thresholds. We demonstrate the utility of SRD annotation in identifying potential false negatives and inaccuracies within current metabolic pathway databases. Our findings highlight the SRD metric as an objective, quantitative and easy-to-compute annotation for gene-metabolite pairs, suitable to integrate statistical evidence to biological networks.

GeneMAP: A discovery platform for metabolic gene function

Organisms maintain metabolic homeostasis through the combined functions of small molecule transporters and enzymes. While many of the metabolic components have been well-established, a substantial number remains without identified physiological substrates. To bridge this gap, we have leveraged large-scale plasma metabolome genome-wide association studies (GWAS) to develop a multiomic Gene-Metabolite Associations Prediction (GeneMAP) discovery platform. GeneMAP can generate accurate predictions, even pinpointing genes that are distant from the variants implicated by GWAS. In particular, our work identified SLC25A48 as a genetic determinant of plasma choline levels. Mechanistically, SLC25A48 loss strongly impairs mitochondrial choline import and synthesis of its downstream metabolite, betaine. Rare variant testing and polygenic risk score analyses have elucidated choline-relevant phenomic consequences of SLC25A48 dysfunction. Altogether, our study proposes SLC25A48 as a mitochondrial choline transporter and provides a discovery platform for metabolic gene function.

Clinical and genetic definition of serum bilirubin levels for the diagnosis of Gilbert syndrome and hypobilirubinemia

BACKGROUND AND AIMS

Gilbert syndrome (GS) is genotypically predetermined by UGT1A1*28 homozygosity in Europeans and is phenotypically defined by hyperbilirubinemia using total bilirubin (TB) cutoff ≥1mg/dL (17 μmol/L). The prevalence of illnesses associated with GS and hypobilirubinemia has never been studied prospectively. As TB varies with UGT1A1*28 genotyping, sex, and age, we propose stratified definitions of TB reference intervals and report the prevalence of illnesses and adjusted 15 years survival.

METHODS

UK Biobank with apparently healthy liver participants (middle-aged, n=138,125) were analyzed after the exclusion of of nonhealthy individuals. The stratified TB was classified as GS when TB >90th centile; <10th centile indicated hypobilirubinemia, and between the 10th and 90th centile was normobilirubinemia. We compared the prevalence and survival rates of 54 illnesses using odds ratio (OR), logistic regression, and Cox models adjusted for confounders, and causality by Mendelian randomizations.

RESULTS

In women, we identified 10% (7,741/76,809) of GS versus 3.7% (2,819/76,809) using the historical cutoff of ≥1 mg/dL (P<0.0001). When GS and hypobilirubinemia participants were compared with normobilirubinemia, after adjustment and Mendelian randomizations, only cholelithiasis prevalence was significantly higher (OR=1.50; 95% CI [1.3-1.7], P=0.001) in men with GS compared with normobilirubinemia and in causal association with bilirubin (P=0.04). No adjusted survival was significantly associated with GS or hypobilirubinemia.

CONCLUSIONS

In middle-aged Europeans, the stratified TB demonstrates a careless GS underestimation in women when using the standard unisex 1 mg/dL cutoff. The prevalence of illnesses is different in GS and hypobilirubinemia as well as survivals before adjusting for confounding factors. With the exception of cholelithiasis in men, these differences were no more significant after adjustment and Mendelian randomization.

AI-powered therapeutic target discovery

Disease modeling and target identification are the most crucial initial steps in drug discovery, and influence the probability of success at every step of drug development. Traditional target identification is a time-consuming process that takes years to decades and usually starts in an academic setting. Given its advantages of analyzing large datasets and intricate biological networks, artificial intelligence (AI) is playing a growing role in modern drug target identification. We review recent advances in target discovery, focusing on breakthroughs in AI-driven therapeutic target exploration. We also discuss the importance of striking a balance between novelty and confidence in target selection. An increasing number of AI-identified targets are being validated through experiments and several AI-derived drugs are entering clinical trials; we highlight current limitations and potential pathways for moving forward.

Approach strategies and application of metabolomics to biotechnology in plants

Metabolomics refers to the technology for the comprehensive analysis of metabolites and low-molecular-weight compounds in a biological system, such as cells or tissues. Metabolites play an important role in biological phenomena through their direct involvement in the regulation of physiological mechanisms, such as maintaining cell homeostasis or signal transmission through protein-protein interactions. The current review aims provide a framework for how the integrated analysis of metabolites, their functional actions and inherent biological information can be used to understand biological phenomena related to the regulation of metabolites and how this information can be applied to safety assessments of crops created using biotechnology. Advancement in technology and analytical instrumentation have led new ways to examine the convergence between biology and chemistry, which has yielded a deeper understanding of complex biological phenomena. Metabolomics can be utilized and applied to safety assessments of biotechnology products through a systematic approach using metabolite-level data processing algorithms, statistical techniques, and database development. The integration of metabolomics data with sequencing data is a key step towards improving additional phenotypical evidence to elucidate the degree of environmental affects for variants found in genome associated with metabolic processes. Moreover, information analysis technology such as big data, machine learning, and IT investment must be introduced to establish a system for data extraction, selection, and metabolomic data analysis for the interpretation of biological implications of biotechnology innovations. This review outlines the integrity of metabolomics assessments in determining the consequences of genetic engineering and biotechnology in plants.

Early hyperlipidemia triggers metabolomic reprogramming with increased SAH, increased acetyl-CoA-cholesterol synthesis, and decreased glycolysis

To identify metabolomic reprogramming in early hyperlipidemia, unbiased metabolome was screened in four tissues from ApoE-/- mice fed with high fat diet (HFD) for 3 weeks. 30, 122, 67, and 97 metabolites in the aorta, heart, liver, and plasma, respectively, were upregulated. 9 upregulated metabolites were uremic toxins, and 13 metabolites, including palmitate, promoted a trained immunity with increased syntheses of acetyl-CoA and cholesterol, increased S-adenosylhomocysteine (SAH) and hypomethylation and decreased glycolysis. The cross-omics analysis found upregulation of 11 metabolite synthetases in ApoE‾/‾ aorta, which promote ROS, cholesterol biosynthesis, and inflammation. Statistical correlation of 12 upregulated metabolites with 37 gene upregulations in ApoE‾/‾ aorta indicated 9 upregulated new metabolites to be proatherogenic. Antioxidant transcription factor NRF2-/- transcriptome analysis indicated that NRF2 suppresses trained immunity-metabolomic reprogramming. Our results have provided novel insights on metabolomic reprogramming in multiple tissues in early hyperlipidemia oriented toward three co-existed new types of trained immunity.

mGWAS-Explorer 2.0: Causal Analysis and Interpretation of Metabolite-Phenotype Associations

Metabolomics-based genome-wide association studies (mGWAS) are key to understanding the genetic regulations of metabolites in complex phenotypes. We previously developed mGWAS-Explorer 1.0 to link single-nucleotide polymorphisms (SNPs), metabolites, genes and phenotypes for hypothesis generation. It has become clear that identifying potential causal relationships between metabolites and phenotypes, as well as providing deep functional insights, are crucial for further downstream applications. Here, we introduce mGWAS-Explorer 2.0 to support the causal analysis between >4000 metabolites and various phenotypes. The results can be interpreted within the context of semantic triples and molecular quantitative trait loci (QTL) data. The underlying R package is released for reproducible analysis. Using two case studies, we demonstrate that mGWAS-Explorer 2.0 is able to detect potential causal relationships between arachidonic acid and Crohn's disease, as well as between glycine and coronary heart disease.

Metabolome-wide Mendelian randomization characterizes heterogeneous and shared causal effects of metabolites on human health

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.

Prioritization of nasal polyp-associated genes by integrating GWAS and eQTL summary data

Background: Nasal polyps (NP) are benign inflammatory growths of nasal and paranasal sinus mucosa that can substantially impair patients' quality of life by various symptoms such as nasal obstruction, insomnia, and anosmia. NP often relapse even after surgical treatment, and the curative therapy would be challenging without understanding the underlying mechanisms. Genome wide association studies (GWASs) on NP have been conducted; however, few genes that are causally associated with NP have been identified. Methods: We aimed to prioritize NP associated genes for functional follow-up studies using the summary data-based Mendelian Randomization (SMR) and Bayesian colocalization (COLOC) methods to integrate the summary-level data of the GWAS on NP and the expression quantitative trait locus (eQTL) study in blood. We utilized the GWAS data including 5,554 NP cases and 258,553 controls with 34 genome-wide significant loci from the FinnGen consortium (data freeze 8) and the eQTL data from 31,684 participants of predominantly European ancestry from the eQTLGen consortium. Results: The SMR analysis identified several genes including TNFRSF18, CTSK, and IRF1 that were associated with NP due to not linkage but pleiotropy or causality. The COLOC analysis strongly suggested that these genes and the trait of NP were affected by shared causal variants, and thus were colocalized. An enrichment analysis by Metascape suggested that these genes might be involved in the biological process of cellular response to cytokine stimulus. Conclusion: We could prioritize several NP associated genes including TNFRSF18, CTSK, and IRF1 for follow-up functional studies in future to elucidate the underlying disease mechanisms.

Plasma metabolomic response to high-carbohydrate meals of differing glycaemic load in overweight women

BACKGROUND

Metabolomic dysregulation following a meal in overweight individuals with the Metabolic Syndrome (MetS) involves multiple pathways of nutrient storage and oxidation.

OBJECTIVE

The aim of the current study was to perform an acute cross-over intervention to examine the interactive actions of meal glycaemic load (GL) on the dynamic responses of the plasma metabolome in overweight females.

METHODS

Postmenopausal women [63 ± 1.23y; Healthy (n = 20) and MetS (n = 20)] ingested two differing high-carbohydrate test meals (73 g carbohydrate; 51% energy) composed of either low glycemic index (LGI) or high (HGI) foods in a randomised sequence. Plasma metabolome was analysed using liquid chromatography-mass spectrometry (LC-MS).

RESULTS

In the overweight women with MetS, there were suppressed postprandial responses for several amino acids (AAs), including phenylalanine, leucine, valine, and tryptophan, p < 0.05), irrespective of the meal type. Meal GL exerted a limited impact on the overall metabolomic response, although the postprandial levels of alanine were higher with the low GL meal and uric acid was greater following the high GL meal (p < 0.05).

CONCLUSIONS

MetS participants exhibited reduced differences in the concentrations of a small set of AAs and a limited group of metabolites implicated in energy metabolism following the meals. However, the manipulation of meal GL had minimal impact on the postprandial metabolome. This study suggests that the GL of a meal is not a major determinant of postprandial response, with a greater impact exerted by the metabolic health of the individual. Trial registration Australia New Zealand Clinical Trials Registry: ACTRN12615001108505 (21/10/2015).

Exploiting the mediating role of the metabolome to unravel transcript-to-phenotype associations

Despite the success of genome-wide association studies (GWASs) in identifying genetic variants associated with complex traits, understanding the mechanisms behind these statistical associations remains challenging. Several methods that integrate methylation, gene expression, and protein quantitative trait loci (QTLs) with GWAS data to determine their causal role in the path from genotype to phenotype have been proposed. Here, we developed and applied a multi-omics Mendelian randomization (MR) framework to study how metabolites mediate the effect of gene expression on complex traits. We identified 216 transcript-metabolite-trait causal triplets involving 26 medically relevant phenotypes. Among these associations, 58% were missed by classical transcriptome-wide MR, which only uses gene expression and GWAS data. This allowed the identification of biologically relevant pathways, such as between ANKH and calcium levels mediated by citrate levels and SLC6A12 and serum creatinine through modulation of the levels of the renal osmolyte betaine. We show that the signals missed by transcriptome-wide MR are found, thanks to the increase in power conferred by integrating multiple omics layer. Simulation analyses show that with larger molecular QTL studies and in case of mediated effects, our multi-omics MR framework outperforms classical MR approaches designed to detect causal relationships between single molecular traits and complex phenotypes.

Prior metabolite extraction fully preserves RNAseq quality and enables integrative multi-'omics analysis of the liver metabolic response to viral infection

Here, we provide an in-depth analysis of the usefulness of single-sample metabolite/RNA extraction for multi-'omics readout. Using pulverized frozen livers of mice injected with lymphocytic choriomeningitis virus (LCMV) or vehicle (Veh), we isolated RNA prior (RNA) or following metabolite extraction (MetRNA). RNA sequencing (RNAseq) data were evaluated for differential expression analysis and dispersion, and differential metabolite abundance was determined. Both RNA and MetRNA clustered together by principal component analysis, indicating that inter-individual differences were the largest source of variance. Over 85% of LCMV versus Veh differentially expressed genes were shared between extraction methods, with the remaining 15% evenly and randomly divided between groups. Differentially expressed genes unique to the extraction method were attributed to randomness around the 0.05 FDR cut-off and stochastic changes in variance and mean expression. In addition, analysis using the mean absolute difference showed no difference in the dispersion of transcripts between extraction methods. Altogether, our data show that prior metabolite extraction preserves RNAseq data quality, which enables us to confidently perform integrated pathway enrichment analysis on metabolomics and RNAseq data from a single sample. This analysis revealed pyrimidine metabolism as the most LCMV-impacted pathway. Combined analysis of genes and metabolites in the pathway exposed a pattern in the degradation of pyrimidine nucleotides leading to uracil generation. In support of this, uracil was among the most differentially abundant metabolites in serum upon LCMV infection. Our data suggest that hepatic uracil export is a novel phenotypic feature of acute infection and highlight the usefulness of our integrated single-sample multi-'omics approach.