Reporter pathway analysis from transcriptome data: Metabolite-centric versus Reaction-centric approach

Multilayer Analysis of RNA Sequencing Data in Alzheimer's Disease to Unravel Molecular Mysteries

Alzheimer's disease (AD) is a complex disease, and numerous cellular events may be involved in etiology. RNAseq-based transcriptome data hold multilayer information content, which could be crucial in unraveling molecular mysteries of AD. It enables quantification of gene expression levels, identification of genomic variants, and elucidation of splicing anomalies such as exon skipping and intron retention. Additional integration of this information into protein-protein interaction networks and genome-scale metabolic models from the literature has potential to decipher functional modules and affected mechanisms for complex scenarios such as AD. In this chapter, we review the application areas of the multilayer content of RNAseq and associated integrative approaches available, with a special focus on AD.

Integrated small RNA, mRNA and protein omics reveal a miRNA network orchestrating metabolic maturation of the developing human heart

BACKGROUND

As the fetal heart develops, cardiomyocyte proliferation potential decreases while fatty acid oxidative capacity increases in a highly regulated transition known as cardiac maturation. Small noncoding RNAs, such as microRNAs (miRNAs), contribute to the establishment and control of tissue-specific transcriptional programs. However, small RNA expression dynamics and genome-wide miRNA regulatory networks controlling maturation of the human fetal heart remain poorly understood.

RESULTS

Transcriptome profiling of small RNAs revealed the temporal expression patterns of miRNA, piRNA, circRNA, snoRNA, snRNA and tRNA in the developing human heart between 8 and 19 weeks of gestation. Our analysis demonstrated that miRNAs were the most dynamically expressed small RNA species throughout mid-gestation. Cross-referencing differentially expressed miRNAs and mRNAs predicted 6200 mRNA targets, 2134 of which were upregulated and 4066 downregulated as gestation progressed. Moreover, we found that downregulated targets of upregulated miRNAs, including hsa-let-7b, miR-1-3p, miR-133a-3p, miR-143-3p, miR-499a-5p, and miR-30a-5p predominantly control cell cycle progression. In contrast, upregulated targets of downregulated miRNAs, including hsa-miR-1276, miR-183-5p, miR-1229-3p, miR-615-3p, miR-421, miR-200b-3p and miR-18a-3p, are linked to energy sensing and oxidative metabolism. Furthermore, integrating miRNA and mRNA profiles with proteomes and reporter metabolites revealed that proteins encoded in mRNA targets and their associated metabolites mediate fatty acid oxidation and are enriched as the heart develops.

CONCLUSIONS

This study presents the first comprehensive analysis of the small RNAome of the maturing human fetal heart. Our findings suggest that coordinated activation and repression of miRNA expression throughout mid-gestation is essential to establish a dynamic miRNA-mRNA-protein network that decreases cardiomyocyte proliferation potential while increasing the oxidative capacity of the maturing human fetal heart. Our results provide novel insights into the molecular control of metabolic maturation of the human fetal heart.

Transcriptome data analysis of primary cardiomyopathies reveals perturbations in arachidonic acid metabolism

Introduction

Cardiomyopathies are complex heart diseases with significant prevalence around the world. Among these, primary forms are the major contributors to heart failure and sudden cardiac death. As a high-energy demanding engine, the heart utilizes fatty acids, glucose, amino acid, lactate and ketone bodies for energy to meet its requirement. However, continuous myocardial stress and cardiomyopathies drive towards metabolic impairment that advances heart failure (HF) pathogenesis. So far, metabolic profile correlation across different cardiomyopathies remains poorly understood.

Methods

In this study, we systematically explore metabolic differences amongst primary cardiomyopathies. By assessing the metabolic gene expression of all primary cardiomyopathies, we highlight the significantly shared and distinct metabolic pathways that may represent specialized adaptations to unique cellular demands. We utilized publicly available RNA-seq datasets to profile global changes in the above diseases (|log2FC| ≥ 0.28 and BH adjusted p-val 0.1) and performed gene set analysis (GSA) using the PAGE statistics on KEGG pathways.

Results

Our analysis demonstrates that genes in arachidonic acid metabolism (AA) are significantly perturbed across cardiomyopathies. In particular, the arachidonic acid metabolism gene PLA2G2A interacts with fibroblast marker genes and can potentially influence fibrosis during cardiomyopathy.

Conclusion

The profound significance of AA metabolism within the cardiovascular system renders it a key player in modulating the phenotypes of cardiomyopathies.

Revealing holistic metabolic responses associated with lipid and docosahexaenoic acid (DHA) production in Aurantiochytrium sp. SW1

Aurantiochytrium sp. SW1, a marine thraustochytrid, has been regarded as a potential candidate as a docosahexaenoic acid (DHA) producer. Even though the genomics of Aurantiochytrium sp. are available, the metabolic responses at a systems level are largely unknown. Therefore, this study aimed to investigate the global metabolic responses to DHA production in Aurantiochytrium sp. through transcriptome and genome-scale network-driven analysis. Of a total of 13,505 genes, 2527 differentially expressed genes (DEGs) were identified in Aurantiochytrium sp., unravelling the transcriptional regulations behinds lipid and DHA accumulation. The highest number of DEG were found for pairwise comparison between growth phase and lipid accumulating phase where a total of 1435 genes were down-regulated with 869 genes being up-regulated. These uncovered several metabolic pathways that contributing in DHA and lipid accumulation including amino acid and acetate metabolism which involve in the generation of crucial precursors. Upon applying network-driven analysis, hydrogen sulphide was found as potential reporter metabolite that could be associated with the genes related to acetyl-CoA synthesis for DHA production. Our findings suggest that the transcriptional regulation of these pathways is a ubiquitous feature in response to specific cultivation phases during DHA overproduction in Aurantiochytrium sp. SW1.

In silico analysis of metabolic effects of bipolar disorder on prefrontal cortex identified altered GABA, glutamate-glutamine cycle, energy metabolism and amino acid synthesis pathways

Bipolar disorder (BP) is a lifelong psychiatric condition, which often disrupts the daily life of the patients. It is characterized by unstable and periodic mood changes, which cause patients to display unusual shifts in mood, energy, activity levels, concentration and the ability to carry out day-to-day tasks. BP is a major psychiatric condition, and it is still undertreated. The causes and neural mechanisms of bipolar disorder are unclear, and diagnosis is still mostly based on psychiatric examination, furthermore the unstable character of the disorder makes diagnosis challenging. Identification of the molecular mechanisms underlying the disease may improve the diagnosis and treatment rates. Single nucleotide polymorphisms (SNP) and transcriptome profiles of patients were studied along with signalling pathways that are thought to be associated with bipolar disorder. Here, we present a computational approach that uses publicly available transcriptome data from bipolar disorder patients and healthy controls. Along with statistical analyses, data are integrated with a genome-scale metabolic model and protein-protein interaction network. Healthy individuals and bipolar disorder patients are compared based on their metabolic profiles. We hypothesize that energy metabolism alterations in bipolar disorder relate to perturbations in amino-acid metabolism and neuron-astrocyte exchange reactions. Due to changes in amino acid metabolism, neurotransmitters and their secretion from neurons and metabolic exchange pathways between neurons and astrocytes such as the glutamine-glutamate cycle are also altered. Changes in negatively charged (-1) KIV and KMV molecules are also observed, and it indicates that charge balance in the brain is highly altered in bipolar disorder. Due to this fact, we also hypothesize that positively charged lithium ions may stabilize the disturbed charge balance in neurons in addition to its effects on neurotransmission. To the best of our knowledge, our approach is unique as it is the first study using genome-scale metabolic models in neuropsychiatric research.

Cancer Recurrence and Omics: Metabolic Signatures of Cancer Dormancy Revealed by Transcriptome Mapping of Genome-Scale Networks

A major problem in medicine and oncology is cancer recurrence through the activation of dormant cancer cells. A system scale examination of metabolic dysregulations associated with the cancer dormancy offers promise for the discovery of novel molecular targets for cancer precision medicine, and importantly, for the prevention of cancer recurrence. In this study, we mapped the total mRNA sequencing-based transcriptomic data from dormant cancer cell lines and nondormant cancer controls onto a human genome-scale metabolic network by using a graph-based approach, and two mass balance-based approaches with one based on reaction activity/inactivity and the other one on flux changes. The gene expression datasets were accessed from Gene Expression Omnibus (GSE83142 and GSE114012). This analysis included two diverse cancer types, a liquid and a solid cancer, namely, acute lymphoblastic leukemia and colorectal cancer. For the dormant cancer state, we observed changes in major adenosine triphosphate-producing pathways, including the citric acid cycle, oxidative phosphorylation, and glycolysis/gluconeogenesis, indicating a reprogramming in the metabolism of dormant cells away from Warburg-based energy metabolism. All three computational approaches unanimously predicted that folate metabolism, pyruvate metabolism, and glutamate metabolism, as well as valine/leucine/isoleucine metabolism are likely dysregulated in cancer dormancy. These findings provide new insights and molecular pathway targets on cancer dormancy, comprehensively catalog dormancy-associated metabolic pathways, and inform future research aimed at prevention of cancer recurrence in particular.

Genome-scale metabolic modelling predicts biomarkers and therapeutic targets for neuropsychiatric disorders

Distinguishing neuropsychiatric disorders is challenging due to the overlap in symptoms and genetic risk factors. People suffering from these disorders face personal and professional challenges. Understanding the dysregulation of brain metabolism under disease condition can aid in effective diagnosis and in developing treatment strategies based on the metabolism. In this study, we reconstructed the metabolic network of three major neuropsychiatric disorders, schizophrenia (SCZ), bipolar disorder (BD) and major depressive disorder (MDD) using transcriptomic data and constrained based modelling approach. We integrated brain transcriptomic data from six independent studies with a recent comprehensive genome-scale metabolic model Recon3D. The analysis of the reconstructed network revealed the flux-level alterations in the peroxisome-mitochondria-golgi axis in neuropsychiatric disorders. We also extracted reporter metabolites and pathways that distinguish these three neuropsychiatric disorders. We found differences with respect to fatty acid oxidation, aromatic and branched chain amino acid metabolism, bile acid synthesis, glycosaminoglycans synthesis and modifications, and phospholipid metabolism. Further, we predicted network perturbations that transform the disease metabolic state to a healthy metabolic state for each disorder. These analyses provide local and global views of the metabolic changes in SCZ, BD and MDD, which may have clinical implications.

Integration of transcriptomic profile of SARS-CoV-2 infected normal human bronchial epithelial cells with metabolic and protein-protein interaction networks

A novel coronavirus (SARS-CoV-2, formerly known as nCoV-2019) that causes an acute respiratory disease has emerged in Wuhan, China and spread globally in early 2020. On January the 30th, the World Health Organization (WHO) declared spread of this virus as an epidemic and a public health emergency. With its highly contagious characteristic and long incubation time, confinement of SARS-CoV-2 requires drastic lock-down measures to be taken and therefore early diagnosis is crucial. We analysed transcriptome of SARS-CoV-2 infected human lung epithelial cells, compared it with mock-infected cells, used network-based reporter metabolite approach and integrated the transcriptome data with protein-protein interaction network to elucidate the early cellular response. Significantly affected metabolites have the potential to be used in diagnostics while pathways of protein clusters have the potential to be used as targets for supportive or novel therapeutic approaches. Our results are in accordance with the literature on response of IL6 family of cytokines and their importance, in addition, we find that matrix metalloproteinase 2 (MMP2) and matrix metalloproteinase 9 (MMP9) with keratan sulfate synthesis pathway may play a key role in the infection. We hypothesize that MMP9 inhibitors have potential to prevent "cytokine storm" in severely affected patients.

Network-based metabolic characterization of renal cell carcinoma

An emerging hallmark of cancer is metabolic reprogramming, which presents opportunities for cancer diagnosis and treatment based on metabolism. We performed a comprehensive metabolic network analysis of major renal cell carcinoma (RCC) subtypes including clear cell, papillary and chromophobe by integrating transcriptomic data with the human genome-scale metabolic model to understand the coordination of metabolic pathways in cancer cells. We identified metabolic alterations of each subtype with respect to tumor-adjacent normal samples and compared them to understand the differences between subtypes. We found that genes of amino acid metabolism and redox homeostasis are significantly altered in RCC subtypes. Chromophobe showed metabolic divergence compared to other subtypes with upregulation of genes involved in glutamine anaplerosis and aspartate biosynthesis. A difference in transcriptional regulation involving HIF1A is observed between subtypes. We identified E2F1 and FOXM1 as other major transcriptional activators of metabolic genes in RCC. Further, the co-expression pattern of metabolic genes in each patient showed the variations in metabolism within RCC subtypes. We also found that co-expression modules of each subtype have tumor stage-specific behavior, which may have clinical implications.

Current Status and Future Prospects of Genome-Scale Metabolic Modeling to Optimize the Use of Mesenchymal Stem Cells in Regenerative Medicine

Mesenchymal stem cells are a promising source for externally grown tissue replacements and patient-specific immunomodulatory treatments. This promise has not yet been fulfilled in part due to production scaling issues and the need to maintain the correct phenotype after re-implantation. One aspect of extracorporeal growth that may be manipulated to optimize cell growth and differentiation is metabolism. The metabolism of MSCs changes during and in response to differentiation and immunomodulatory changes. MSC metabolism may be linked to functional differences but how this occurs and influences MSC function remains unclear. Understanding how MSC metabolism relates to cell function is however important as metabolite availability and environmental circumstances in the body may affect the success of implantation. Genome-scale constraint based metabolic modeling can be used as a tool to fill gaps in knowledge of MSC metabolism, acting as a framework to integrate and understand various data types (e.g., genomic, transcriptomic and metabolomic). These approaches have long been used to optimize the growth and productivity of bacterial production systems and are being increasingly used to provide insights into human health research. Production of tissue for implantation using MSCs requires both optimized production of cell mass and the understanding of the patient and phenotype specific metabolic situation. This review considers the current knowledge of MSC metabolism and how it may be optimized along with the current and future uses of genome scale constraint based metabolic modeling to further this aim.

Metabolite-Centric Reporter Pathway and Tripartite Network Analysis of Arabidopsis Under Cold Stress

The study of plant resistance to cold stress and the metabolic processes underlying its molecular mechanisms benefit crop improvement programs. Here we investigate the effects of cold stress on the metabolic pathways of Arabidopsis when directly inferred at system level from transcriptome data. A metabolite-centric reporter pathway analysis approach enabled the computation of metabolites associated with transcripts at four time points of cold treatment. Tripartite networks of gene-metabolite-pathway connectivity outlined the response of metabolites and pathways to cold stress. Our metabolome-independent analysis revealed stress-associated metabolites in pathway routes of the cold stress response, including amino acid, carbohydrate, lipid, hormone, energy, photosynthesis, and signaling pathways. Cold stress first triggered the mobilization of energy from glycolysis and ethanol degradation to enhance TCA cycle activity via acetyl-CoA. Interestingly, tripartite networks lacked power law behavior and scale free connectivity, favoring modularity. Network rewiring explicitly involved energetics, signal, carbon and redox metabolisms and membrane remodeling.

Systems-level organization of non-alcoholic fatty liver disease progression network

Non-Alcoholic Fatty Liver Disease (NAFLD) is a complex spectrum of diseases ranging from simple steatosis to Non-Alcoholic Steatohepatitis (NASH) with fibrosis, which can progress to cirrhosis and hepatocellular carcinoma. The pathogenesis of NAFLD is complex, involving crosstalk between multiple organs, cell-types, and environmental and genetic factors. Dysfunction of the adipose tissue plays a central role in NAFLD progression. Here, we analysed transcriptomics data obtained from the Visceral Adipose Tissue (VAT) of NAFLD patients to understand how the VAT metabolism is altered at the genome scale and co-regulated with other cellular processes during the progression from obesity to NASH with fibrosis. For this purpose, we performed Weighted Gene Co-expression Network Analysis (WGCNA), a method that organizes the disease transcriptome into functional modules of cellular processes and pathways. Our analysis revealed the coordination of metabolic and inflammatory modules (termed "immunometabolism") in the VAT of NAFLD patients. We found that genes of arachidonic acid, sphingolipid and glycosphingolipid metabolism were upregulated and co-expressed with genes of proinflammatory signalling pathways and hypoxia in NASH/NASH with fibrosis. We hypothesize that these metabolic alterations might play a role in sustaining VAT inflammation. Furthermore, immunometabolism related genes were also co-expressed with genes involved in Extracellular Matrix (ECM) degradation. Our analysis indicates that upregulation of both ECM degrading enzymes and their inhibitors (incoherent feedforward loop) potentially leads to the ECM deposition in the VAT of NASH with fibrosis patients.

Genome-Scale Brain Metabolic Networks as Scaffolds for the Systems Biology of Neurodegenerative Diseases: Mapping Metabolic Alterations

Systems-based investigation of diseases requires integrated analysis of cellular networks and high-throughput data of gene products. The use of genome-scale metabolic networks for such integration has led to the elucidation of cellular mechanisms for several cell types from microorganisms to plants. It has become easier and cheaper to generate high-throughput data over years in the form of transcriptome, proteome and metabolome. This has tremendously improved the quality and quantity of information extracted from such data enabling the documentation of active pathways and reactions in cell metabolism. A number of omics-based datasets for several neurodegenerative diseases are now available in public repositories. This increases the potential of using genome-scale brain metabolic networks as a scaffold for this type of data to map metabolic alterations for the purpose of elucidating disease mechanisms and for the diagnosis and treatment of such disorders. This chapter first reviews omics data collected for neurodegenerative diseases to map their effect on metabolism. Later, the potential for genome-scale metabolic modeling of such data is reviewed and discussed in light of recently reconstructed brain metabolic networks at genome-scale.