Topological analysis of metabolic networks integrating co-segregating transcriptomes and metabolomes in type 2 diabetic rat congenic series

Transcriptomics reveals dynamic changes in the "gene profiles" of rat supraspinatus tendon at three different time points after diabetes induction

OBJECTIVE

There is increasing evidence that type 2 diabetes mellitus (T2DM) is an independent risk factor for the occur of tendinopathy. Therefore, this study is the first to explore the dynamic changes of the "gene profile" of supraspinatus tendon in rats at different time points after T2DM induction through transcriptomics, providing potential molecular markers for exploring the pathogenesis of diabetic tendinopathy.

METHODS

A total of 40 Sprague-Dawley rats were randomly divided into normal (NG, n = 10) and T2DM groups (T2DM, n = 30) and subdivided into three groups according to the duration of diabetes: T2DM-4w, T2DM-8w, and T2DM-12w groups; the duration was calculated from the time point of T2DM rat model establishment. The three comparison groups were set up in this study, T2DM-4w group vs. NG, T2DM-8w group vs. NG, and T2DM-12w group vs. NG. Differentially expressed genes (DEGs) in 3 comparison groups were screened. The intersection of the three comparison groups' DEGs was defined as key genes that changed consistently in the supraspinatus tendon after diabetes induction. Cluster analysis, gene ontology (GO) functional annotation analysis and Kyoto encyclopedia of genes and genomes (KEGG) functional annotation and enrichment analysis were performed for DEGs.

RESULTS

T2DM-4w group vs. NG, T2DM-8w group vs. NG, and T2DM-12w group vs. NG detected 519 (251 up-regulated and 268 down-regulated), 459 (342 up-regulated and 117 down-regulated) and 328 (255 up-regulated and 73 down-regulated) DEGs, respectively. 103 key genes of sustained changes in the supraspinatus tendon following induction of diabetes, which are the first identified biomarkers of the supraspinatus tendon as it progresses through the course of diabetes.The GO analysis results showed that the most significant enrichment in biological processes was calcium ion transmembrane import into cytosol (3 DEGs). The most significant enrichment in cellular component was extracellular matrix (9 DEGs). The most significant enrichment in molecular function was glutamate-gated calcium ion channel activity (3 DEGs). The results of KEGG pathway enrichment analysis showed that there were 17 major pathways (p < 0.05) that diabetes affected supratinusculus tendinopathy, including cAMP signaling pathway and Calcium signaling pathway.

CONCLUSIONS

Transcriptomics reveals dynamic changes in the"gene profiles"of rat supraspinatus tendon at three different time points after diabetes induction. The 103 DEGs identified in this study may provide potential molecular markers for exploring the pathogenesis of diabetic tendinopathy, and the 17 major pathways enriched in KEGG may provide new ideas for exploring the pathogenesis of diabetic tendinopathy.

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.

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

Studies combining metabolomics and genetics, known as metabolite genome-wide association studies (mGWAS), have provided valuable insights into our understanding of the genetic control of metabolite levels. However, the biological interpretation of these associations remains challenging due to a lack of existing tools to annotate mGWAS gene-metabolite pairs beyond the use of conservative statistical significance threshold. Here, we computed the shortest reactional distance (SRD) based on the curated knowledge of the KEGG database to explore its utility in enhancing the biological interpretation of results from three independent mGWAS, including a case study on sickle cell disease patients. Results show that, in reported mGWAS pairs, there is an excess of small SRD values and that SRD values and p-values significantly correlate, even beyond the standard conservative thresholds. The added-value of SRD annotation is shown for identification of potential false negative hits, exemplified by the finding of gene-metabolite associations with SRD ≤1 that did not reach standard genome-wide significance cut-off. The wider use of this statistic as an mGWAS annotation would prevent the exclusion of biologically relevant associations and can also identify errors or gaps in current metabolic pathway databases. Our findings highlight the SRD metric as an objective, quantitative and easy-to-compute annotation for gene-metabolite pairs that can be used to integrate statistical evidence to biological networks.

The regulatory landscape of multiple brain regions in outbred heterogeneous stock rats

Heterogeneous Stock (HS) rats are a genetically diverse outbred rat population that is widely used for studying genetics of behavioral and physiological traits. Mapping Quantitative Trait Loci (QTL) associated with transcriptional changes would help to identify mechanisms underlying these traits. We generated genotype and transcriptome data for five brain regions from 88 HS rats. We identified 21 392 cis-QTLs associated with expression and splicing changes across all five brain regions and validated their effects using allele specific expression data. We identified 80 cases where eQTLs were colocalized with genome-wide association study (GWAS) results from nine physiological traits. Comparing our dataset to human data from the Genotype-Tissue Expression (GTEx) project, we found that the HS rat data yields twice as many significant eQTLs as a similarly sized human dataset. We also identified a modest but highly significant correlation between genetic regulatory variation among orthologous genes. Surprisingly, we found less genetic variation in gene regulation in HS rats relative to humans, though we still found eQTLs for the orthologs of many human genes for which eQTLs had not been found. These data are available from the RatGTEx data portal (RatGTEx.org) and will enable new discoveries of the genetic influences of complex traits.

Multi-Omic Approaches to Identify Genetic Factors in Metabolic Syndrome

Metabolic syndrome (MetS) is a highly heritable disease and a major public health burden worldwide. MetS diagnosis criteria are met by the simultaneous presence of any three of the following: high triglycerides, low HDL/high LDL cholesterol, insulin resistance, hypertension, and central obesity. These diseases act synergistically in people suffering from MetS and dramatically increase risk of morbidity and mortality due to stroke and cardiovascular disease, as well as certain cancers. Each of these component features is itself a complex disease, as is MetS. As a genetically complex disease, genetic risk factors for MetS are numerous, but not very powerful individually, often requiring specific environmental stressors for the disease to manifest. When taken together, all sequence variants that contribute to MetS disease risk explain only a fraction of the heritable variance, suggesting additional, novel loci have yet to be discovered. In this article, we will give a brief overview on the genetic concepts needed to interpret genome-wide association studies (GWAS) and quantitative trait locus (QTL) data, summarize the state of the field of MetS physiological genomics, and to introduce tools and resources that can be used by the physiologist to integrate genomics into their own research on MetS and any of its component features. There is a wealth of phenotypic and molecular data in animal models and humans that can be leveraged as outlined in this article. Integrating these multi-omic QTL data for complex diseases such as MetS provides a means to unravel the pathways and mechanisms leading to complex disease and promise for novel treatments. © 2022 American Physiological Society. Compr Physiol 12:1-40, 2022.

Fatty Acid Synthase Correlates With Prognosis-Related Abdominal Adipose Distribution and Metabolic Disorders of Clear Cell Renal Cell Carcinoma

Purpose: Lipid metabolism reprogramming is a major pathway in tumor evolution. This study investigated fatty acid synthase (FASN) mRNA expression in anthropometric adipose tissue and elucidated the prognostic value and potential mechanism of clear cell renal cell carcinoma (ccRCC).

Materials and Methods: Transcription profiles were obtained from 533 ccRCC samples in The Cancer Genome Atlas (TCGA) cohorts. Real-time quantitative PCR (RT-qPCR) and immunohistochemistry were performed to detect FASN expression in 380 paired ccRCC and normal tissues from the Fudan University Shanghai Cancer Center (FUSCC). Visceral adipose tissue (VAT) and subcutaneous adipose tissue were at the level of the umbilicus as measured by magnetic resonance imaging (MRI). Non-targeted metabolomics and in vitro experiments were used to reveal the biological functions of FASN.

Results: Increased FASN expression was significantly relevant to advanced T, N, and American Joint Committee on Cancer (AJCC) stages (p < 0.01) and significantly correlated to poor progression-free survival (PFS) and overall survival (OS) of 913 ccRCC patients in FUSCC and TCGA cohorts. Pearson's correlation coefficient indicated that FASN amplification was positively correlated to VAT% (r = 0.772, p < 0.001), which significantly correlated to poor PFS (HR = 2.066, p = 0.028) and OS (HR = 2.773, p = 0.023) in the FUSCC cohort. Transient inhibition or overexpression of FASN significantly regulated A498 and 786O cell proliferation and migration by regulating epithelial–mesenchymal transition. Inhibition of FASN led to a higher apoptotic rate and decreased lipid droplet formation compared with normal control in ccRCC cells. Non-targeted metabolomics showed that decreased de novo lipogenesis might be required to sustain an elevation of glycolytic activity in 786O cells by regulating galactinol, dl-lactate, N-acetylaspartylglutamate, and sucrose, thereby participating in carcinogenesis and progression of ccRCC.

Conclusion: This study demonstrated that FASN expression is positively related to aggressive cell proliferation, migration, apoptosis, and lipid droplet formation and regulates metabolic disorders of the ccRCC microenvironment. Additionally, elevated FASN mRNA expression is significantly correlated to the abdominal obesity distribution, especially VAT%, which is a significant predictor of a poor prognosis for ccRCC patients.

A large-scale genome-lipid association map guides lipid identification

Despite the crucial roles of lipids in metabolism, we are still at the early stages of comprehensively annotating lipid species and their genetic basis. Mass spectrometry-based discovery lipidomics offers the potential to globally survey lipids and their relative abundances in various biological samples. To discover the genetics of lipid features obtained through high-resolution liquid chromatography-tandem mass spectrometry, we analysed liver and plasma from 384 diversity outbred mice, and quantified 3,283 molecular features. These features were mapped to 5,622 lipid quantitative trait loci and compiled into a public web resource termed LipidGenie. The data are cross-referenced to the human genome and offer a bridge between genetic associations in humans and mice. Harnessing this resource, we used genome-lipid association data as an additional aid to identify a number of lipids, for example gangliosides through their association with B4galnt1, and found evidence for a group of sex-specific phosphatidylcholines through their shared locus. Finally, LipidGenie's ability to query either mass or gene-centric terms suggests acyl-chain-specific functions for proteins of the ABHD family.

CMap analysis identifies Atractyloside as a potential drug candidate for type 2 diabetes based on integration of metabolomics and transcriptomics

Background

This research aimed at exploring the mechanisms of alterations of metabolites and pathways in T2D from the perspective of metabolomics and transcriptomics, as well as uncovering novel drug candidate for T2D treatment.

Methods

Metabolites in human plasma from 42 T2D patients and 45 non‐diabetic volunteers were detected by liquid chromatography‐mass spectrometer (LC‐MS). Microarray dataset of the transcriptome was obtained from Gene Expression Omnibus (GEO) database. Kyoto Encyclopedia of Genes and Genomes (KEGG) database was used to conduct pathway enrichment analysis. Connectivity Map (CMap) was employed to select potential drugs for T2D therapy. In vivo assay was performed to verify above findings. The protein expression levels of ME1, ME2 and MDH1 were detected by Western blot to determine the status of NAD/NADH cofactor system.

Results

In our study, differentially expressed metabolites were selected out between healthy samples and T2D samples with selection criteria P value < .05, |Fold Change| > 2, including N‐acetylglutamate and Malate. Genes set enrichment analysis (GSEA) revealed that 34 pathways were significantly enriched in T2D. Based on CMap analysis and animal experiments, Atractyloside was identified as a potential novel drug for T2D treatment via targeting ME1, ME2 and MDH1 and regulating the NAD/NADH cofactor system.

Conclusion

The present research revealed differentially expressed metabolites and genes, as well as significantly altered pathways in T2D via an integration of metabolomics, transcriptomics and CMap analysis. It was also demonstrated that comprehensive analysis based on metabolomics and transcriptomics was an effective approach for identification and verification of metabolic biomarkers and alternated pathways.

Systems genetics of the Drosophila metabolome

How effects of DNA sequence variants are transmitted through intermediate endophenotypes to modulate organismal traits remains a central question in quantitative genetics. This problem can be addressed through a systems approach in a population in which genetic polymorphisms, gene expression traits, metabolites, and complex phenotypes can be evaluated on the same genotypes. Here, we focused on the metabolome, which represents the most proximal link between genetic variation and organismal phenotype, and quantified metabolite levels in 40 lines of the Drosophila melanogaster Genetic Reference Panel. We identified sex-specific modules of genetically correlated metabolites and constructed networks that integrate DNA sequence variation and variation in gene expression with variation in metabolites and organismal traits, including starvation stress resistance and male aggression. Finally, we asked to what extent SNPs and metabolites can predict trait phenotypes and generated trait- and sex-specific prediction models that provide novel insights about the metabolomic underpinnings of complex phenotypes.

Metabolomics technology and bioinformatics for precision medicine

Precision medicine is rapidly emerging as a strategy to tailor medical treatment to a small group or even individual patients based on their genetics, environment and lifestyle. Precision medicine relies heavily on developments in systems biology and omics disciplines, including metabolomics. Combination of metabolomics with sophisticated bioinformatics analysis and mathematical modeling has an extreme power to provide a metabolic snapshot of the patient over the course of disease and treatment or classifying patients into subpopulations and subgroups requiring individual medical treatment. Although a powerful approach, metabolomics have certain limitations in technology and bioinformatics. We will review various aspects of metabolomics technology and bioinformatics, from data generation, bioinformatics analysis, data fusion and mathematical modeling to data management, in the context of precision medicine.

STATegra, a comprehensive multi-omics dataset of B-cell differentiation in mouse

Multi-omics approaches use a diversity of high-throughput technologies to profile the different molecular layers of living cells. Ideally, the integration of this information should result in comprehensive systems models of cellular physiology and regulation. However, most multi-omics projects still include a limited number of molecular assays and there have been very few multi-omic studies that evaluate dynamic processes such as cellular growth, development and adaptation. Hence, we lack formal analysis methods and comprehensive multi-omics datasets that can be leveraged to develop true multi-layered models for dynamic cellular systems. Here we present the STATegra multi-omics dataset that combines measurements from up to 10 different omics technologies applied to the same biological system, namely the well-studied mouse pre-B-cell differentiation. STATegra includes high-throughput measurements of chromatin structure, gene expression, proteomics and metabolomics, and it is complemented with single-cell data. To our knowledge, the STATegra collection is the most diverse multi-omics dataset describing a dynamic biological system.

Conserved properties of genetic architecture of renal and fat transcriptomes in rat models of insulin resistance

To define renal molecular mechanisms that are affected by permanent hyperglycaemia and might promote phenotypes relevant to diabetic nephropathy, we carried out linkage analysis of genome-wide gene transcription in the kidneys of F2 offspring from the Goto-Kakizaki (GK) rat model of type 2 diabetes and normoglycaemic Brown Norway (BN) rats. We mapped 2526 statistically significant expression quantitative trait loci (eQTLs) in the cross. More than 40% of eQTLs mapped in the close vicinity of the linked transcripts, underlying possible cis-regulatory mechanisms of gene expression. We identified eQTL hotspots on chromosomes 5 and 9 regulating the expression of 80-165 genes, sex or cross direction effects, and enriched metabolic and immunological processes by segregating GK alleles. Comparative analysis with adipose tissue eQTLs in the same cross showed that 496 eQTLs, in addition to the top enriched biological pathways, are conserved in the two tissues. Extensive similarities in eQTLs mapped in the GK rat and in the spontaneously hypertensive rat (SHR) suggest a common aetiology of disease phenotypes common to the two strains, including insulin resistance, which is a prominent pathophysiological feature in both GK rats and SHRs. Our data shed light on shared and tissue-specific molecular mechanisms that might underlie aetiological aspects of insulin resistance in the context of spontaneously occurring hyperglycaemia and hypertension.

Summary: Kidney and fat expression QTL mapping in rat models of spontaneously occurring insulin resistance associated with either diabetes or hypertension reveals conserved gene expression regulation, suggesting shared aetiology of disease phenotypes.

Simultaneous Extraction of RNA and Metabolites from Single Kidney Tissue Specimens for Combined Transcriptomic and Metabolomic Profiling

Tissue analysis represents a powerful tool for the investigation of disease pathophysiology. However, the heterogeneous nature of tissue samples, in particular of neoplastic, may affect the outcome of such analysis and hence obscure interpretation of results. Thus, comprehensive isolation and extraction of transcripts and metabolites from an identical tissue specimen would minimize variations and enable the economic use of biopsy material which is usually available in limited amounts. Here we demonstrate a fast and simple protocol for combined transcriptomics and metabolomics analysis in homogenates prepared from one single tissue sample. Metabolites were recovered by protein precipitation from lysates originally prepared for RNA isolation and were analyzed by LC-QTOF-MS after HILIC and RPLC separation, respectively. Strikingly, although ion suppression was observed, over 80% of the 2885 detected metabolic features could be extracted and analyzed with high reproducibility (CV ≤ 20%). Moreover fold changes of different tumor and nontumor kidney tissues were correlated to an established metabolomics protocol and revealed a strong correlation ( r_p ≥ 0.75). In order to demonstrate the feasibility of the combined analysis of RNA and metabolites, the protocol was applied to kidney tissue of metformin treated mice to investigate drug induced alterations.

MetaboSignal: a network-based approach for topological analysis of metabotype regulation via metabolic and signaling pathways

Summary

MetaboSignal is an R package that allows merging metabolic and signaling pathways reported in the Kyoto Encyclopaedia of Genes and Genomes (KEGG). It is a network-based approach designed to navigate through topological relationships between genes (signaling- or metabolic-genes) and metabolites, representing a powerful tool to investigate the genetic landscape of metabolic phenotypes.

Availability and Implementation

MetaboSignal is available from Bioconductor: https://bioconductor.org/packages/MetaboSignal/

Supplementary information

Supplementary data are available at Bioinformatics online.