Transcriptional and metabolic data integration and modeling for identification of active pathways

Transcriptomics- and metabolomics-based integration analyses revealed the potential pharmacological effects and functional pattern of in vivo Radix Paeoniae Alba administration

Background

Radix Paeoniae Alba (RPA) and other natural medicines have remarkable curative effects and are widely used in traditional Chinese Medicine (TCM). However, due to their multi-component and multi-target characteristics, it is difficult to study the detailed pharmacological mechanisms for those natural medicines in vivo. Therefore, their real effects on organisms is still uncertain.

Methods

RPA was selected as research object, the present study was designed to study the complex mechanisms of RPA in vivo by integrating and interpreting the transcriptomic based RNA-seq and metabolomic based NMR spectrum after RPA administration in mice. A variety of dimension-reduction algorithms and classifier models were applied to the processing of high-throughput data.

Results

Among serum metabolites, the contents of PC and glucose were significantly increased, while the contents of various amino acids, lipids and their metabolites were significantly decreased in mice after RPA administration. Based on the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, differential analysis showed that the liver was the site where RPA exerted a significant effect, which confirmed the rationality of “meridian tropism” in the theory in TCM. In addition, RPA played a role in lipid metabolism by regulating genes encoding enzymes of the glycerolipid metabolism pathway, such as 1-acyl-sn-glycerol-3-phosphate acyltransferase (Agpat), phosphatidate phosphatase (Lpin), phospholipid phosphatase (Plpp) and endothelial lipase (Lipg). We also found that RPA regulates several substance addiction pathways in the brain, such as the cocaine addiction pathway, and the related targets were predicted based on the sequencing data from pathological model in the GEO database. The overall effective pattern of RPA was intuitively presented with a multidimensional radar map through a self-designed model which found that liver and brain were mainly regulated by RPA compared with the traditional meridian tropism theory.

Conclusions

Overall this study expanded the potential application of RPA and provided possible targets and directions for further mechanism study, meanwhile, it also established a multi-dimensional evaluation model to represent the overall effective pattern of TCM for the first time. In the future, such study based on the high-throughput data sets can be used to interpret the theory of TCM and to provide a valuable research model and clinical medication reference for the TCM researchers and doctors.

Review of Statistical Learning Methods in Integrated Omics Studies (An Integrated Information Science)

Integrated omics is becoming a new channel for investigating the complex molecular system in modern biological science and sets a foundation for systematic learning for precision medicine. The statistical/machine learning methods that have emerged in the past decade for integrated omics are not only innovative but also multidisciplinary with integrated knowledge in biology, medicine, statistics, machine learning, and artificial intelligence. Here, we review the nontrivial classes of learning methods from the statistical aspects and streamline these learning methods within the statistical learning framework. The intriguing findings from the review are that the methods used are generalizable to other disciplines with complex systematic structure, and the integrated omics is part of an integrated information science which has collated and integrated different types of information for inferences and decision making. We review the statistical learning methods of exploratory and supervised learning from 42 publications. We also discuss the strengths and limitations of the extended principal component analysis, cluster analysis, network analysis, and regression methods. Statistical techniques such as penalization for sparsity induction when there are fewer observations than the number of features and using Bayesian approach when there are prior knowledge to be integrated are also included in the commentary. For the completeness of the review, a table of currently available software and packages from 23 publications for omics are summarized in the appendix.

Hierarchical Feature Selection Incorporating Known and Novel Biological Information: Identifying Genomic Features Related to Prostate Cancer Recurrence

Our work is motivated by a prostate cancer study aimed at identifying mRNA and miRNA biomarkers that are predictive of cancer recurrence after prostatectomy. It has been shown in the literature that incorporating known biological information on pathway memberships and interactions among biomarkers improves feature selection of high-dimensional biomarkers in relation to disease risk. Biological information is often represented by graphs or networks, in which biomarkers are represented by nodes and interactions among them are represented by edges; however, biological information is often not fully known. For example, the role of microRNAs (miRNAs) in regulating gene expression is not fully understood and the miRNA regulatory network is not fully established, in which case new strategies are needed for feature selection. To this end, we treat unknown biological information as missing data (i.e., missing edges in graphs), different from commonly encountered missing data problems where variable values are missing. We propose a new concept of imputing unknown biological information based on observed data and define the imputed information as the novel biological information. In addition, we propose a hierarchical group penalty to encourage sparsity and feature selection at both the pathway level and the within-pathway level, which, combined with the imputation step, allows for incorporation of known and novel biological information. While it is applicable to general regression settings, we develop and investigate the proposed approach in the context of semiparametric accelerated failure time models motivated by our data example. Data application and simulation studies show that incorporation of novel biological information improves performance in risk prediction and feature selection and the proposed penalty outperforms the extensions of several existing penalties.

Pathway crosstalk effects: Shrinkage and disentanglement using a Bayesian hierarchical model

Identifying the biological pathways that are related to various clinical phenotypes is an important concern in biomedical research. Based on estimated expression levels and/or p-values, over-representation analysis (ORA) methods provide rankings of pathways, but they are tainted because pathways overlap. This crosstalk phenomenon has not been rigorously studied and classical ORA does not take into consideration: (i) that crosstalk effects in cases of overlapping pathways can cause incorrect rankings of pathways, (ii) that crosstalk effects can cause both excess type I errors and type II errors, (iii) that rankings of small pathways are unreliable and (iv) that type I error rates can be inflated due to multiple comparisons of pathways. We develop a Bayesian hierarchical model that addresses these problems, providing sensible estimates and rankings, and reducing error rates. We show, on both real and simulated data, that the results of our method are more accurate than the results produced by the classical over-representation analysis, providing a better understanding of the underlying biological phenomena involved in the phenotypes under study. The R code and the binary datasets for implementing the analyses described in this article are available online at: http://www.eng.wayne.edu/page.php?id=6402.

A non-negative matrix factorization method for detecting modules in heterogeneous omics multi-modal data

MOTIVATION

Recent advances in high-throughput omics technologies have enabled biomedical researchers to collect large-scale genomic data. As a consequence, there has been growing interest in developing methods to integrate such data to obtain deeper insights regarding the underlying biological system. A key challenge for integrative studies is the heterogeneity present in the different omics data sources, which makes it difficult to discern the coordinated signal of interest from source-specific noise or extraneous effects.

RESULTS

We introduce a novel method of multi-modal data analysis that is designed for heterogeneous data based on non-negative matrix factorization. We provide an algorithm for jointly decomposing the data matrices involved that also includes a sparsity option for high-dimensional settings. The performance of the proposed method is evaluated on synthetic data and on real DNA methylation, gene expression and miRNA expression data from ovarian cancer samples obtained from The Cancer Genome Atlas. The results show the presence of common modules across patient samples linked to cancer-related pathways, as well as previously established ovarian cancer subtypes.

AVAILABILITY AND IMPLEMENTATION

The source code repository is publicly available at https://github.com/yangzi4/iNMF.

CONTACT

gmichail@umich.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The discordant method: a novel approach for differential correlation

MOTIVATION

Current differential correlation methods are designed to determine molecular feature pairs that have the largest magnitude of difference between correlation coefficients. These methods do not easily capture molecular feature pairs that experience no correlation in one group but correlation in another, which may reflect certain types of biological interactions. We have developed a tool, the Discordant method, which categorizes the correlation types for each group to make this possible.

RESULTS

We compare the Discordant method to existing approaches using simulations and two biological datasets with different types of -omics data. In contrast to other methods, Discordant identifies phenotype-related features at a similar or higher rate while maintaining reasonable computational tractability and usability.

AVAILABILITY AND IMPLEMENTATION

R code and sample data are available at https://github.com/siskac/discordant

CONTACT

katerina.kechris@ucdenver.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Transcriptomic and metabolomic data integration

Many studies now produce parallel data sets from different omics technologies; however, the task of interpreting the acquired data in an integrated fashion is not trivial. This review covers those methods that have been used over the past decade to statistically integrate and interpret metabolomics and transcriptomic data sets. It defines four categories of approaches, correlation-based integration, concatenation-based integration, multivariate-based integration and pathway-based integration, into which all existing statistical methods fit. It also explores the choices in study design for generating samples for analysis by these omics technologies and the impact that these technical decisions have on the subsequent data analysis options.

Modern plant metabolomics: advanced natural product gene discoveries, improved technologies, and future prospects

Plant metabolomics has matured and modern plant metabolomics has accelerated gene discoveries and the elucidation of a variety of plant natural product biosynthetic pathways. This review covers the approximate period of 2000 to 2014, and highlights specific examples of the discovery and characterization of novel genes and enzymes associated with the biosynthesis of natural products such as flavonoids, glucosinolates, terpenoids, and alkaloids. Additional examples of the integration of metabolomics with genome-based functional characterizations of plant natural products that are important to modern pharmaceutical technology are also reviewed. This article also provides a substantial review of recent technical advances in mass spectrometry imaging, nuclear magnetic resonance imaging, integrated LC-MS-SPE-NMR for metabolite identifications, and X-ray crystallography of microgram quantities for structural determinations. The review closes with a discussion on the future prospects of metabolomics related to crop species and herbal medicine.

Comparative analysis and modeling of the severity of steatohepatitis in DDC-treated mouse strains

Background

Non-alcoholic fatty liver disease (NAFLD) has a broad spectrum of disease states ranging from mild steatosis characterized by an abnormal retention of lipids within liver cells to steatohepatitis (NASH) showing fat accumulation, inflammation, ballooning and degradation of hepatocytes, and fibrosis. Ultimately, steatohepatitis can result in liver cirrhosis and hepatocellular carcinoma.

Methodology and Results

In this study we have analyzed three different mouse strains, A/J, C57BL/6J, and PWD/PhJ, that show different degrees of steatohepatitis when administered a 3,5-diethoxycarbonyl-1,4-dihydrocollidine (DDC) containing diet. RNA-Seq gene expression analysis, protein analysis and metabolic profiling were applied to identify differentially expressed genes/proteins and perturbed metabolite levels of mouse liver samples upon DDC-treatment. Pathway analysis revealed alteration of arachidonic acid (AA) and S-adenosylmethionine (SAMe) metabolism upon other pathways. To understand metabolic changes of arachidonic acid metabolism in the light of disease expression profiles a kinetic model of this pathway was developed and optimized according to metabolite levels. Subsequently, the model was used to study in silico effects of potential drug targets for steatohepatitis.

Conclusions

We identified AA/eicosanoid metabolism as highly perturbed in DDC-induced mice using a combination of an experimental and in silico approach. Our analysis of the AA/eicosanoid metabolic pathway suggests that 5-hydroxyeicosatetraenoic acid (5-HETE), 15-hydroxyeicosatetraenoic acid (15-HETE) and prostaglandin D2 (PGD2) are perturbed in DDC mice. We further demonstrate that a dynamic model can be used for qualitative prediction of metabolic changes based on transcriptomics data in a disease-related context. Furthermore, SAMe metabolism was identified as being perturbed due to DDC treatment. Several genes as well as some metabolites of this module show differences between A/J and C57BL/6J on the one hand and PWD/PhJ on the other.

A global approach to analysis and interpretation of metabolic data for plant natural product discovery

Discovering molecular components and their functionality is key to the development of hypotheses concerning the organization and regulation of metabolic networks. The iterative experimental testing of such hypotheses is the trajectory that can ultimately enable accurate computational modelling and prediction of metabolic outcomes. This information can be particularly important for understanding the biology of natural products, whose metabolism itself is often only poorly defined. Here, we describe factors that must be in place to optimize the use of metabolomics in predictive biology. A key to achieving this vision is a collection of accurate time-resolved and spatially defined metabolite abundance data and associated metadata. One formidable challenge associated with metabolite profiling is the complexity and analytical limits associated with comprehensively determining the metabolome of an organism. Further, for metabolomics data to be efficiently used by the research community, it must be curated in publicly available metabolomics databases. Such databases require clear, consistent formats, easy access to data and metadata, data download, and accessible computational tools to integrate genome system-scale datasets. Although transcriptomics and proteomics integrate the linear predictive power of the genome, the metabolome represents the nonlinear, final biochemical products of the genome, which results from the intricate system(s) that regulate genome expression. For example, the relationship of metabolomics data to the metabolic network is confounded by redundant connections between metabolites and gene-products. However, connections among metabolites are predictable through the rules of chemistry. Therefore, enhancing the ability to integrate the metabolome with anchor-points in the transcriptome and proteome will enhance the predictive power of genomics data. We detail a public database repository for metabolomics, tools and approaches for statistical analysis of metabolomics data, and methods for integrating these datasets with transcriptomic data to create hypotheses concerning specialized metabolisms that generate the diversity in natural product chemistry. We discuss the importance of close collaborations among biologists, chemists, computer scientists and statisticians throughout the development of such integrated metabolism-centric databases and software.

Inter-individual differences in response to dietary intervention: integrating omics platforms towards personalised dietary recommendations

Technologic advances now make it possible to collect large amounts of genetic, epigenetic, metabolomic and gut microbiome data. These data have the potential to transform approaches towards nutrition counselling by allowing us to recognise and embrace the metabolic, physiologic and genetic differences among individuals. The ultimate goal is to be able to integrate these multi-dimensional data so as to characterise the health status and disease risk of an individual and to provide personalised recommendations to maximise health. To this end, accurate and predictive systems-based measures of health are needed that incorporate molecular signatures of genes, transcripts, proteins, metabolites and microbes. Although we are making progress within each of these omics arenas, we have yet to integrate effectively multiple sources of biologic data so as to provide comprehensive phenotypic profiles. Observational studies have provided some insights into associative interactions between genetic or phenotypic variation and diet and their impact on health; however, very few human experimental studies have addressed these relationships. Dietary interventions that test prescribed diets in well-characterised study populations and that monitor system-wide responses (ideally using several omics platforms) are needed to make correlation-causation connections and to characterise phenotypes under controlled conditions. Given the growth in our knowledge, there is the potential to develop personalised dietary recommendations. However, developing these recommendations assumes that an improved understanding of the phenotypic complexities of individuals and their responses to the complexities of their diets will lead to a sustainable, effective approach to promote health and prevent disease - therein lies our challenge.