Integration of gene expression data with prior knowledge for network analysis and validation

Background

Reconstruction of protein-protein interaction or metabolic networks based on expression data often involves in silico predictions, while on the other hand, there are unspecific networks of in vivo interactions derived from knowledge bases.We analyze networks designed to come as close as possible to data measured in vivo, both with respect to the set of nodes which were taken to be expressed in experiment as well as with respect to the interactions between them which were taken from manually curated databases

Results

A signaling network derived from the TRANSPATH database and a metabolic network derived from KEGG LIGAND are each filtered onto expression data from breast cancer (SAGE) considering different levels of restrictiveness in edge and vertex selection.We perform several validation steps, in particular we define pathway over-representation tests based on refined null models to recover functional modules. The prominent role of the spindle checkpoint-related pathways in breast cancer is exhibited. High-ranking key nodes cluster in functional groups retrieved from literature. Results are consistent between several functional and topological analyses and between signaling and metabolic aspects.

Conclusions

This construction involved as a crucial step the passage to a mammalian protein identifier format as well as to a reaction-based semantics of metabolism. This yielded good connectivity but also led to the need to perform benchmark tests to exclude loss of essential information. Such validation, albeit tedious due to limitations of existing methods, turned out to be informative, and in particular provided biological insights as well as information on the degrees of coherence of the networks despite fragmentation of experimental data.Key node analysis exploited the networks for potentially interesting proteins in view of drug target prediction.

     

Comprehensive annotation of microRNA expression profiles

Background

MicroRNAs (miRNAs) regulate many biological processes by post-translational gene silencing. Analysis of miRNA expression profiles is a reliable method for investigating particular biological processes due to the stability of miRNA and the development of advanced sequencing methods. However, this approach is limited by the broad specificity of miRNAs, which may target several mRNAs.

Result

In this study, we developed a method for comprehensive annotation of miRNA array or deep sequencing data for investigation of cellular biological effects. Using this method, the specific pathways and biological processes involved in Alzheimer’s disease were predicted with high correlation in four independent samples. Furthermore, this method was validated for evaluation of cadmium telluride (CdTe) nanomaterial cytotoxicity. As a result, apoptosis pathways were selected as the top pathways associated with CdTe nanoparticle exposure, which is consistent with previous studies.

Conclusions

Our findings contribute to the validation of miRNA microarray or deep sequencing results for early diagnosis of disease and evaluation of the biological safety of new materials and drugs.

     

Comparison,alignment, and synchronization of cell line information between CLO and EFO

Background

The Experimental Factor Ontology (EFO) is an application ontology driven by experimental variables including cell lines to organize and describe the diverse experimental variables and data resided in the EMBL-EBI resources. The Cell Line Ontology (CLO) is an OBO community-based ontology that contains information of immortalized cell lines and relevant experimental components. EFO integrates and extends ontologies from the bio-ontology community to drive a number of practical applications. It is desirable that the community shares design patterns and therefore that EFO reuses the cell line representation from the Cell Line Ontology (CLO). There are, however, challenges to be addressed when developing a common ontology design pattern for representing cell lines in both EFO and CLO.

Results

In this study, we developed a strategy to compare and map cell line terms between EFO and CLO. We examined Cellosaurus resources for EFO-CLO cross-references. Text labels of cell lines from both ontologies were verified by biological information axiomatized in each source. The study resulted in the identification 873 EFO-CLO aligned and 344 EFO unique immortalized permanent cell lines. All of these cell lines were updated to CLO and the cell line related information was merged. A design pattern that integrates EFO and CLO was also developed.

Conclusion

Our study compared, aligned, and synchronized the cell line information between CLO and EFO. The final updated CLO will be examined as the candidate ontology to import and replace eligible EFO cell line classes thereby supporting the interoperability in the bio-ontology domain. Our mapping pipeline illustrates the use of ontology in aiding biological data standardization and integration through the biological and semantics content of cell lines.

     

Use artificial neural network to align biological ontologies

Background

Being formal, declarative knowledge representation models, ontologies help to address the problem of imprecise terminologies in biological and biomedical research. However, ontologies constructed under the auspices of the Open Biomedical Ontologies (OBO) group have exhibited a great deal of variety, because different parties can design ontologies according to their own conceptual views of the world. It is therefore becoming critical to align ontologies from different parties. During automated/semi-automated alignment across biological ontologies, different semantic aspects, i.e., concept name, concept properties, and concept relationships, contribute in different degrees to alignment results. Therefore, a vector of weights must be assigned to these semantic aspects. It is not trivial to determine what those weights should be, and current methodologies depend a lot on human heuristics.

Results

In this paper, we take an artificial neural network approach to learn and adjust these weights, and thereby support a new ontology alignment algorithm, customized for biological ontologies, with the purpose of avoiding some disadvantages in both rule-based and learning-based aligning algorithms. This approach has been evaluated by aligning two real-world biological ontologies, whose features include huge file size, very few instances, concept names in numerical strings, and others.

Conclusion

The promising experiment results verify our proposed hypothesis, i.e., three weights for semantic aspects learned from a subset of concepts are representative of all concepts in the same ontology. Therefore, our method represents a large leap forward towards automating biological ontology alignment.

     

Pathway-based analyses

Background

New technologies for acquisition of genomic data, while offering unprecedented opportunities for genetic discovery, also impose severe burdens of interpretation andpenalties for multiple testing.

Methods

The Pathway-based Analyses Group of the Genetic Analysis Workshop 19 (GAW19) sought reduction of multiple-testing burden through various approaches to aggregation of highdimensional data in pathways informed by prior biological knowledge.

Results

Experimental methods testedincluded the use of "synthetic pathways" (random sets of genes) to estimate power and false-positive error rate of methods applied to simulated data; data reduction via independent components analysis, single-nucleotide polymorphism (SNP)-SNP interaction, and use of gene sets to estimate genetic similarity; and general assessment of the efficacy of prior biological knowledge to reduce the dimensionality of complex genomic data.

Conclusions

The work of this group explored several promising approaches to managing high-dimensional data, with the caveat that these methods are necessarily constrained by the quality of external bioinformatic annotation.

     

Identification of protein complexes from multi-relationship protein interaction networks

Background

Protein complexes play an important role in biological processes. Recent developments in experiments have resulted in the publication of many high-quality, large-scale protein-protein interaction (PPI) datasets, which provide abundant data for computational approaches to the prediction of protein complexes. However, the precision of protein complex prediction still needs to be improved due to the incompletion and noise in PPI networks.

Results

There exist complex and diverse relationships among proteins after integrating multiple sources of biological information. Considering that the influences of different types of interactions are not the same weight for protein complex prediction, we construct a multi-relationship protein interaction network (MPIN) by integrating PPI network topology with gene ontology annotation information. Then, we design a novel algorithm named MINE (identifying protein complexes based on Multi-relationship protein Interaction NEtwork) to predict protein complexes with high cohesion and low coupling from MPIN.

Conclusions

The experiments on yeast data show that MINE outperforms the current methods in terms of both accuracy and statistical significance.

     

Generalized logical model based on network topology to capture the dynamical trends of cellular signaling pathways

Background

Cellular responses to extracellular perturbations require signaling pathways to capture and transmit the signals. However, the underlying molecular mechanisms of signal transduction are not yet fully understood, thus detailed and comprehensive models may not be available for all the signaling pathways. In particular, insufficient knowledge of parameters, which is a long-standing hindrance for quantitative kinetic modeling necessitates the use of parameter-free methods for modeling and simulation to capture dynamic properties of signaling pathways.

Results

We present a computational model that is able to simulate the graded responses to degradations, the sigmoidal biological relationships between signaling molecules and the effects of scheduled perturbations to the cells. The simulation results are validated using experimental data of protein phosphorylation, demonstrating that the proposed model is capable of capturing the main trend of protein activities during the process of signal transduction. Compared with existing simulators, our model has better performance on predicting the state transitions of signaling networks.

Conclusion

The proposed simulation tool provides a valuable resource for modeling cellular signaling pathways using a knowledge-based method.

     

How the chemical features of molecules may have addressed the settlement of metabolic steps

Introduction

While the evolutionary adaptation of enzymes to their own substrates is a well assessed and rationalized field, how molecules have been originally selected in order to initiate and assemble convenient metabolic pathways is a fascinating, but still debated argument.

Objectives

Aim of the present study is to give a rationale for the preferential selection of specific molecules to generate metabolic pathways.

Methods

The comparison of structural features of molecules, through an inductive methodological approach, offer a reading key to cautiously propose a determining factor for their metabolic recruitment.

Results

Starting with some commonplaces occurring in the structural representation of relevant carbohydrates, such as glucose, fructose and ribose, arguments are presented in associating stable structural determinants of these molecules and their peculiar occurrence in metabolic pathways.

Conclusions

Among other possible factors, the reliability of the structural asset of a molecule may be relevant or its selection among structurally and, a priori, functionally similar molecules.

     

Untargeted and stable isotope-assisted metabolomic analysis of MDA-MB-231 cells under hypoxia

Introduction

Hypoxia commonly occurs in cancers and is highly related with the occurrence, development and metastasis of cancer. Treatment of triple negative breast cancer remains challenge. Knowledge about the metabolic status of triple negative breast cancer cell lines in hypoxia is valuable for the understanding of molecular mechanisms of this tumor subtype to develop effective therapeutics.

Objectives

Comprehensively characterize the metabolic profiles of triple negative breast cancer cell line MDA-MB-231 in normoxia and hypoxia and the pathways involved in metabolic changes in hypoxia.

Methods

Differences in metabolic profiles affected pathways of MDA-MB-231 cells in normoxia and hypoxia were characterized using GC–MS based untargeted and stable isotope assisted metabolomic techniques.

Results

Thirty-three metabolites were significantly changed in hypoxia and nine pathways were involved. Hypoxia increased glycolysis, inhibited TCA cycle, pentose phosphate pathway and pyruvate carboxylation, while increased glutaminolysis in MDA-MB-231 cells.

Conclusion

The current results provide metabolic differences of MDA-MB-231 cells in normoxia and hypoxia conditions as well as the involved metabolic pathways, demonstrating the power of combined use of untargeted and stable isotope-assisted metabolomic methods in comprehensive metabolomic analysis.

     

Prediction of heme binding residues from protein sequences with integrative sequence profiles

Background

The heme-protein interactions are essential for various biological processes such as electron transfer, catalysis, signal transduction and the control of gene expression. The knowledge of heme binding residues can provide crucial clues to understand these activities and aid in functional annotation, however, insufficient work has been done on the research of heme binding residues from protein sequence information.

Methods

We propose a sequence-based approach for accurate prediction of heme binding residues by a novel integrative sequence profile coupling position specific scoring matrices with heme specific physicochemical properties. In order to select the informative physicochemical properties, we design an intuitive feature selection scheme by combining a greedy strategy with correlation analysis.

Results

Our integrative sequence profile approach for prediction of heme binding residues outperforms the conventional methods using amino acid and evolutionary information on the 5-fold cross validation and the independent tests.

Conclusions

The novel feature of an integrative sequence profile achieves good performance using a reduced set of feature vector elements.

     

Untargeted metabolomics suffers from incomplete raw data processing

Introduction

Untargeted metabolomics is a powerful tool for biological discoveries. To analyze the complex raw data, significant advances in computational approaches have been made, yet it is not clear how exhaustive and reliable the data analysis results are.

Objectives

Assessment of the quality of raw data processing in untargeted metabolomics.

Methods

Five published untargeted metabolomics studies, were reanalyzed.

Results

Omissions of at least 50 relevant compounds from the original results as well as examples of representative mistakes were reported for each study.

Conclusion

Incomplete raw data processing shows unexplored potential of current and legacy data.

     

An environment for relation mining over richly annotated corpora: the case of GENIA

Background

The biomedical domain is witnessing a rapid growth of the amount of published scientific results, which makes it increasingly difficult to filter the core information. There is a real need for support tools that 'digest' the published results and extract the most important information.

Results

We describe and evaluate an environment supporting the extraction of domain-specific relations, such as protein-protein interactions, from a richly-annotated corpus. We use full, deep-linguistic parsing and manually created, versatile patterns, expressing a large set of syntactic alternations, plus semantic ontology information.

Conclusion

The experiments show that our approach described is capable of delivering high-precision results, while maintaining sufficient levels of recall. The high level of abstraction of the rules used by the system, which are considerably more powerful and versatile than finite-state approaches, allows speedy interactive development and validation.

     

A methodology for elucidating regulatory mechanisms leading to changes in lipid profiles

Introduction

It is difficult to elucidate the metabolic and regulatory factors causing lipidome perturbations.

Objectives

This work simplifies this process.

Methods

A method has been developed to query an online holistic lipid metabolic network (of 7923 metabolites) to extract the pathways that connect the input list of lipids.

Results

The output enables pathway visualisation and the querying of other databases to identify potential regulators. When used to a study a plasma lipidome dataset of polycystic ovary syndrome, 14 enzymes were identified, of which 3 are linked to ELAVL1—an mRNA stabiliser.

Conclusion

This method provides a simplified approach to identifying potential regulators causing lipid-profile perturbations.

     

From correlation to causation: analysis of metabolomics data using systems biology approaches

Introduction

Metabolomics is a well-established tool in systems biology, especially in the top–down approach. Metabolomics experiments often results in discovery studies that provide intriguing biological hypotheses but rarely offer mechanistic explanation of such findings. In this light, the interpretation of metabolomics data can be boosted by deploying systems biology approaches.

Objectives

This review aims to provide an overview of systems biology approaches that are relevant to metabolomics and to discuss some successful applications of these methods.

Methods

We review the most recent applications of systems biology tools in the field of metabolomics, such as network inference and analysis, metabolic modelling and pathways analysis.

Results

We offer an ample overview of systems biology tools that can be applied to address metabolomics problems. The characteristics and application results of these tools are discussed also in a comparative manner.

Conclusions

Systems biology-enhanced analysis of metabolomics data can provide insights into the molecular mechanisms originating the observed metabolic profiles and enhance the scientific impact of metabolomics studies.

     

Targeted metabolomic profiling of low and high grade serous epithelial ovarian cancer tissues: a pilot study

Introduction

Epithelial ovarian cancer (EOC) remains the leading cause of death from gynecologic malignancies and has an alarming global fatality rate. Besides the differences in underlying pathogenesis, distinguishing between high grade (HG) and low grade (LG) EOC is imperative for the prediction of disease progression and responsiveness to chemotherapy.

Objectives

The aim of this study was to investigate, the tissue metabolome associated with HG and LG serous epithelial ovarian cancer.

Methods

A combination of one dimensional proton nuclear magnetic resonance (1D H NMR) spectroscopy and targeted mass spectrometry (MS) was employed to profile the tissue metabolome of HG, LG serous EOCs, and controls.

Results

Using partial least squares-discriminant analysis, we observed significant separation between all groups (p?<?0.05) following cross validation. We identified which metabolites were significantly perturbed in each EOC grade as compared with controls and report the biochemical pathways which were perturbed due to the disease. Among these metabolic pathways, ascorbate and aldarate metabolism was identified, for the first time, as being significantly altered in both LG and HG serous cancers. Further, we have identified potential biomarkers of EOC and generated predictive algorithms with AUC (CI)?=?0.940 and 0.929 for HG and LG, respectively.

Conclusion

These previously unreported biochemical changes provide a framework for future metabolomic studies for the development of EOC biomarkers. Finally, pharmacologic targeting of the key metabolic pathways identified herein could lead to novel and effective treatments of EOC.

     

Profiling the interactome of protein kinase C ζ by proteomics and bioinformatics

Background

Protein kinase C ζ (PKCζ), an isoform of the atypical protein kinase C, is a pivotal regulator in cancer. However, the molecular and cellular mechanisms whereby PKCζ regulates tumorigenesis and metastasis are still not fully understood. In this study, proteomics and bioinformatics analyses were performed to establish a protein-protein interaction (PPI) network associated with PKCζ, laying a stepping stone to further understand the diverse biological roles of PKCζ.

Methods

Protein complexes associated with PKCζ were purified by co-immunoprecipitation from breast cancer cell MDA-MB-231 and identified by LC-MS/MS. Two biological replicates and two technical replicates were analyzed. The observed proteins were filtered using the CRAPome database to eliminate the potential false positives. The proteomics identification results were combined with PPI database search to construct the interactome network. Gene ontology (GO) and pathway analysis were performed by PANTHER database and DAVID. Next, the interaction between PKCζ and protein phosphatase 2 catalytic subunit alpha (PPP2CA) was validated by co-immunoprecipitation, Western blotting and immunofluorescence. Furthermore, the TCGA database and the COSMIC database were used to analyze the expressions of these two proteins in clinical samples.

Results

The PKCζ centered PPI network containing 178 nodes and 1225 connections was built. Network analysis showed that the identified proteins were significantly associated with several key signaling pathways regulating cancer related cellular processes.

Conclusions

Through combining the proteomics and bioinformatics analyses, a PKCζ centered PPI network was constructed, providing a more complete picture regarding the biological roles of PKCζ in both cancer regulation and other aspects of cellular biology.

     

Genome-wide discovery of missing genes in biological pathways of prokaryotes

Background

Reconstruction of biological pathways is typically done through mapping well-characterized pathways of model organisms to a target genome, through orthologous gene mapping. A limitation of such pathway-mapping approaches is that the mapped pathway models are constrained by the composition of the template pathways, e.g., some genes in a target pathway may not have corresponding genes in the template pathways, the so-called “missing gene” problem.

Methods

We present a novel pathway-expansion method for identifying additional genes that are possibly involved in a target pathway after pathway mapping, to fill holes caused by missing genes as well as to expand the mapped pathway model. The basic idea of the algorithm is to identify genes in the target genome whose homologous genes share common operons with homologs of any mapped pathway genes in some reference genome, and to add such genes to the target pathway if their functions are consistent with the cellular function of the target pathway.

Results

We have implemented this idea using a graph-theoretic approach and demonstrated the effectiveness of the algorithm on known pathways of E. coli in the KEGG database. On all KEGG pathways containing at least 5 genes, our method achieves an average of 60% positive predictive value (PPV) and the performance is increased with more seed genes added. Analysis shows that our method is highly robust.

Conclusions

An effective method is presented to find missing genes in biological pathways of prokaryotes, which achieves high prediction reliability on E. coli at a genome level. Numerous missing genes are found to be related to knwon E. coli pathways, which can be further validated through biological experiments. Overall this method is robust and can be used for functional inference.

     

Prediction of intravenous busulfan clearance by endogenous plasma biomarkers using global pharmacometabolomics

Introduction

High-dose busulfan (busulfan) is an integral part of the majority of hematopoietic cell transplantation conditioning regimens. Intravenous (IV) busulfan doses are personalized using pharmacokinetics (PK)-guided dosing where the patient’s IV busulfan clearance is calculated after the first dose and is used to personalize subsequent doses to a target plasma exposure. PK-guided dosing has improved patient outcomes and is clinically accepted but highly resource-intensive.

Objective

We sought to discover endogenous plasma biomarkers predictive of IV busulfan clearance using a global pharmacometabolomics-based approach

Methods

Using LC-QTOF, we analyzed 59 (discovery) and 88 (validation) plasma samples obtained before IV busulfan administration.

Results

In the discovery dataset, we evaluated the association of the relative abundance of 1885 ions with IV busulfan clearance and found 21 ions that were associated with IV busulfan clearance tertiles (r² ≥ 0.3). Identified compounds were deoxycholic acid and/or chenodeoxycholic acid, and linoleic acid. We used these 21 ions to develop a parsimonious seven-ion linear predictive model that accurately predicted IV busulfan clearance in 93 % (discovery) and 78 % (validation) of samples.

Conclusion

IV busulfan clearance was significantly correlated with the relative abundance of 21 ions, seven of which were included in a predictive model that accurately predicted IV busulfan clearance in the majority of the validation samples. These results reinforce the potential of pharmacometabolomics as a critical tool in personalized medicine, with the potential to improve the personalized dosing of drugs with a narrow therapeutic index such as busulfan.

     

Tracing conformational changes in proteins

Background

Many proteins undergo extensive conformational changes as part of their functionality. Tracing these changes is important for understanding the way these proteins function. Traditional biophysics-based conformational search methods require a large number of calculations and are hard to apply to large-scale conformational motions.

Results

In this work we investigate the application of a robotics-inspired method, using backbone and limited side chain representation and a coarse grained energy function to trace large-scale conformational motions. We tested the algorithm on four well known medium to large proteins and we show that even with relatively little information we are able to trace low-energy conformational pathways efficiently. The conformational pathways produced by our methods can be further filtered and refined to produce more useful information on the way proteins function under physiological conditions.

Conclusions

The proposed method effectively captures large-scale conformational changes and produces pathways that are consistent with experimental data and other computational studies. The method represents an important first step towards a larger scale modeling of more complex biological systems.

     

Metabolites predict cardiovascular disease events in persons living with HIV: a pilot case–control study

Introduction

Persons living with HIV (PLWH) are at higher risk for cardiovascular disease (CVD) events than uninfected persons. Current risk-stratification methods to define PLWH at highest risk for CVD events are lacking.

Methods

Using tandem flow injection mass spectrometry, we quantified plasma levels of 60 metabolites in 24 matched pairs of PLWH [1:1 with and without known coronary artery disease (CAD)]. Metabolite levels were reduced to interpretable factors using principal components analysis.

Results

Factors derived from short-chain dicarboxylacylcarnitines (SCDA) (p?=?0.08) and glutamine/valine (p?=?0.003) were elevated in CAD cases compared to controls.

Conclusion

SCDAs and glutamine/valine may be valuable markers of cardiovascular risk among persons living with HIV in the future, pending validation in larger cohorts.