Identification of aberrant pathways and network activities from high-throughput data

Characteristics of genomic mutations and signaling pathway alterations in thymic epithelial tumors

BACKGROUND

To elucidate the mechanisms of thymic epithelial tumor (TET) canceration by characterizing genomic mutations and signaling pathway alterations.

METHODS

Primary tumor and blood samples were collected from 21 patients diagnosed with TETs (thymoma and thymic cancer), 15 of whom were screened by nucleic acid extraction and whole exon sequencing. Bioinformatics was used to comprehensively analyze the sequencing data for these samples, including gene mutation information and the difference of tumor mutation burden (TMB) between thymoma and thymic carcinoma groups. We performed signaling pathway and functional enrichment analysis using the WebGestalt 2017 toolkit.

RESULTS

ZNF429 (36%) was the gene with the highest mutation frequency in thymic carcinoma. Mutations in BAP1 (14%), ABI1 (7%), BCL9L (7%), and CHEK2 (7%) were exclusively detected in thymic carcinoma, whereas ZNF721 mutations (14%) and PABPC1 (14%) were found exclusively in thymoma. The mean TMB values for thymic carcinoma and thymoma were 0.722 and 0.663 mutations per megabase (Mb), respectively, and these differences were not statistically significant. The ErbB signaling pathway was enriched in the thymoma and intersection groups, and pathways of central carbon metabolism in cancer, longevity regulating and MAPK signaling were only found in the thymoma group, while pathways in cancer (hsa05200) was found in the thymoma and thymic carcinoma groups.

CONCLUSIONS

Multiple differences in somatic genes and pathways have been identified. Our findings provide insights into differences between thymoma and thymic carcinoma that could aid in designing personalized clinical therapeutic strategies.

System, Method and Software for Calculation of a Cannabis Drug Efficiency Index for the Reduction of Inflammation

There are many varieties of Cannabis sativa that differ from each other by composition of cannabinoids, terpenes and other molecules. The medicinal properties of these cultivars are often very different, with some being more efficient than others. This report describes the development of a method and software for the analysis of the efficiency of various cannabis extracts to detect the anti-inflammatory properties of the various cannabis extracts. The method uses high-throughput gene expression profiling data but can potentially use other omics data as well. According to the signaling pathway topology, the gene expression profiles are convoluted into the signaling pathway activities using a signaling pathway impact analysis (SPIA) method. The method was tested by inducing inflammation in human 3D epithelial tissues, including intestine, oral and skin, and then exposing these tissues to various extracts and then performing transcriptome analysis. The analysis showed a different efficiency of the various extracts in restoring the transcriptome changes to the pre-inflammation state, thus allowing to calculate a different cannabis drug efficiency index (CDEI).

Multi-omic Pathway and Network Analysis to Identify Biomarkers for Hepatocellular Carcinoma

The threat of Hepatocellular Carcinoma (HCC) is a growing problem, with incidence rates anticipated to near double over the next two decades. The increasing burden makes discovery of novel diagnostic, prognostic, and therapeutic biomarkers distinguishing HCC from underlying cirrhosis a significant focus. In this study, we analyzed tissue and serum samples from 40 HCC cases and 25 patients with liver cirrhosis (CIRR) to better understand the mechanistic differences between HCC and CIRR. Through pathway and network analysis, we are able to take a systems biology approach to conduct multi-omic analysis of transcriptomic, glycoproteomic, and metabolomic data acquired through various platforms. As a result, we are able to identify the FXR/RXR Activation pathway as being represented by molecules spanning multiple molecular compartments in these samples. Specifically, serum metabolites deoxycholate and chenodeoxycholic acid and serum glycoproteins C4A/C4B, KNG1, and HPX are biomarker candidates identified from this analysis that are of interest for future targeted studies. These results demonstrate the integrative power of multi-omic analysis to prioritize clinically and biologically relevant biomarker candidates that can increase understanding of molecular mechanisms driving HCC and make an impact in patient care.

Network-Based Protein Biomarker Discovery Platforms

The advances in mass spectrometry-based proteomics technologies have enabled the generation of global proteome data from tissue or body fluid samples collected from a broad spectrum of human diseases. Comparative proteomic analysis of global proteome data identifies and prioritizes the proteins showing altered abundances, called differentially expressed proteins (DEPs), in disease samples, compared to control samples. Protein biomarker candidates that can serve as indicators of disease states are then selected as key molecules among these proteins. Recently, it has been addressed that cellular pathways can provide better indications of disease states than individual molecules and also network analysis of the DEPs enables effective identification of cellular pathways altered in disease conditions and key molecules representing the altered cellular pathways. Accordingly, a number of network-based approaches to identify disease-related pathways and representative molecules of such pathways have been developed. In this review, we summarize analytical platforms for network-based protein biomarker discovery and key components in the platforms.

Using activation status of signaling pathways as mechanism-based biomarkers to predict drug sensitivity

Many complex traits, as drug response, are associated with changes in biological pathways rather than being caused by single gene alterations. Here, a predictive framework is presented in which gene expression data are recoded into activity statuses of signal transduction circuits (sub-pathways within signaling pathways that connect receptor proteins to final effector proteins that trigger cell actions). Such activity values are used as features by a prediction algorithm which can efficiently predict a continuous variable such as the IC50 value. The main advantage of this prediction method is that the features selected by the predictor, the signaling circuits, are themselves rich-informative, mechanism-based biomarkers which provide insight into or drug molecular mechanisms of action (MoA).

Data integration in biological research: an overview

Data sharing, integration and annotation are essential to ensure the reproducibility of the analysis and interpretation of the experimental findings. Often these activities are perceived as a role that bioinformaticians and computer scientists have to take with no or little input from the experimental biologist. On the contrary, biological researchers, being the producers and often the end users of such data, have a big role in enabling biological data integration. The quality and usefulness of data integration depend on the existence and adoption of standards, shared formats, and mechanisms that are suitable for biological researchers to submit and annotate the data, so it can be easily searchable, conveniently linked and consequently used for further biological analysis and discovery. Here, we provide background on what is data integration from a computational science point of view, how it has been applied to biological research, which key aspects contributed to its success and future directions.

A proteomics approach to identifying key protein targets involved in VEGF inhibitor mediated attenuation of bleomycin-induced pulmonary fibrosis

Idiopathic pulmonary fibrosis (IPF) is a progressive lung disease with a life expectancy of less than 5 years post diagnosis for most patients. Poor molecular characterization of IPF has led to insufficient understanding of the pathogenesis of the disease, resulting in lack of effective therapies. In this study, we have integrated a label-free LC-MS based approach with systems biology to identify signaling pathways and regulatory nodes within protein interaction networks that govern phenotypic changes that may lead to IPF. Ingenuity Pathway Analysis of proteins modulated in response to bleomycin treatment identified PI3K/Akt and Wnt signaling as the most significant profibrotic pathways. Similar analysis of proteins modulated in response to vascular endothelial growth factor (VEGF) inhibitor (CBO-P11) treatment identified natural killer cell signaling and PTEN signaling as the most significant antifibrotic pathways. Mechanistic/mammalian target of rapamycin (mTOR) and extracellular signal-regulated kinase (ERK) were identified to be key mediators of pro- and antifibrotic response, where bleomycin (BLM) treatment resulted in increased expression and VEGF inhibitor treatment attenuated expression of mTOR and ERK. Using a BLM mouse model of pulmonary fibrosis and VEGF inhibitor CBO-P11 as a therapeutic measure, we identified a comprehensive set of signaling pathways and proteins that contribute to the pathogenesis of pulmonary fibrosis that can be targeted for therapy against this fatal disease.

Identification of functional modules by integration of multiple data sources using a Bayesian network classifier

BACKGROUND

Prediction of functional modules is indispensable for detecting protein deregulation in human complex diseases such as cancer. Bayesian network is one of the most commonly used models to integrate heterogeneous data from multiple sources such as protein domain, interactome, functional annotation, genome-wide gene expression, and the literature.

METHODS AND RESULTS

In this article, we present a Bayesian network classifier that is customized to (1) increase the ability to integrate diverse information from different sources, (2) effectively predict protein-protein interactions, (3) infer aberrant networks with scale-free and small-world properties, and (4) group molecules into functional modules or pathways based on the primary function and biological features. Application of this model in discovering protein biomarkers of hepatocellular carcinoma leads to the identification of functional modules that provide insights into the mechanism of the development and progression of hepatocellular carcinoma. These functional modules include cell cycle deregulation, increased angiogenesis (eg, vascular endothelial growth factor, blood vessel morphogenesis), oxidative metabolic alterations, and aberrant activation of signaling pathways involved in cellular proliferation, survival, and differentiation.

CONCLUSIONS

The discoveries and conclusions derived from our customized Bayesian network classifier are consistent with previously published results. The proposed approach for determining Bayesian network structure facilitates the integration of heterogeneous data from multiple sources to elucidate the mechanisms of complex diseases.

Data- and knowledge-based modeling of gene regulatory networks: an update

Gene regulatory network inference is a systems biology approach which predicts interactions between genes with the help of high-throughput data. In this review, we present current and updated network inference methods focusing on novel techniques for data acquisition, network inference assessment, network inference for interacting species and the integration of prior knowledge. After the advance of Next-Generation-Sequencing of cDNAs derived from RNA samples (RNA-Seq) we discuss in detail its application to network inference. Furthermore, we present progress for large-scale or even full-genomic network inference as well as for small-scale condensed network inference and review advances in the evaluation of network inference methods by crowdsourcing. Finally, we reflect the current availability of data and prior knowledge sources and give an outlook for the inference of gene regulatory networks that reflect interacting species, in particular pathogen-host interactions.

Applications of systems approaches in the study of rheumatic diseases

The complex interaction of molecules within a biological system constitutes a functional module. These modules are then acted upon by both internal and external factors, such as genetic and environmental stresses, which under certain conditions can manifest as complex disease phenotypes. Recent advances in high-throughput biological analyses, in combination with improved computational methods for data enrichment, functional annotation, and network visualization, have enabled a much deeper understanding of the mechanisms underlying important biological processes by identifying functional modules that are temporally and spatially perturbed in the context of disease development. Systems biology approaches such as these have produced compelling observations that would be impossible to replicate using classical methodologies, with greater insights expected as both the technology and methods improve in the coming years. Here, we examine the use of systems biology and network analysis in the study of a wide range of rheumatic diseases to better understand the underlying molecular and clinical features.

Pathway and network approaches for identification of cancer signature markers from omics data

The advancement of high throughput omic technologies during the past few years has made it possible to perform many complex assays in a much shorter time than the traditional approaches. The rapid accumulation and wide availability of omic data generated by these technologies offer great opportunities to unravel disease mechanisms, but also presents significant challenges to extract knowledge from such massive data and to evaluate the findings. To address these challenges, a number of pathway and network based approaches have been introduced. This review article evaluates these methods and discusses their application in cancer biomarker discovery using hepatocellular carcinoma (HCC) as an example.

Inferring pathway crosstalk networks using gene set co-expression signatures

Constructing molecular interaction networks in cells is important for understanding the underlying mechanisms of biological processes. Except for single gene analysis, several gene set-based methods have been proposed to infer pathway crosstalk by analyzing large-scale gene expression data. But most of them take all pathway genes as a whole to infer the crosstalk. Biological evidence suggests that the pathway crosstalk usually occurs between some subsets rather than the whole sets of pathway genes. In this study, we propose a novel method, sGSCA (signature-based gene set co-expression analysis) which can use the co-expression correlations between subsets of pathway genes to infer the pathway crosstalk networks. The method applies sparse canonical correlation analysis (sCCA) to measure the pathway level co-expression and simultaneously obtain the subsets or signature genes that contribute to the co-expression of pathways. On simulated datasets, sGSCA can efficiently detect pathway crosstalk and the corresponding highly correlated signature genes. We applied sGSCA to two cancer gene expression datasets (one for hepatocellular cancer and the other for lung cancer). In the inferred networks, we found several important pathway crosstalks related to the cancers. The identified signature genes also show high enrichment for the cancer related genes. sGSCA can infer pathway crosstalk networks using large-scale gene expression data, and should be a useful tool for systematically studying the molecular mechanisms of complex diseases on both pathway and gene levels at the same time.

Multiscale models of breast cancer progression

Breast cancer initiation, invasion and metastasis span multiple length and time scales. Molecular events at short length scales lead to an initial tumorigenic population, which left unchecked by immune action, acts at increasingly longer length scales until eventually the cancer cells escape from the primary tumor site. This series of events is highly complex, involving multiple cell types interacting with (and shaping) the microenvironment. Multiscale mathematical models have emerged as a powerful tool to quantitatively integrate the convective-diffusion-reaction processes occurring on the systemic scale, with the molecular signaling processes occurring on the cellular and subcellular scales. In this study, we reviewed the current state of the art in cancer modeling across multiple length scales, with an emphasis on the integration of intracellular signal transduction models with pro-tumorigenic chemical and mechanical microenvironmental cues. First, we reviewed the underlying biomolecular origin of breast cancer, with a special emphasis on angiogenesis. Then, we summarized the development of tissue engineering platforms which could provide high-fidelity ex vivo experimental models to identify and validate multiscale simulations. Lastly, we reviewed top-down and bottom-up multiscale strategies that integrate subcellular networks with the microenvironment. We present models of a variety of cancers, in addition to breast cancer specific models. Taken together, we expect as the sophistication of the simulations increase, that multiscale modeling and bottom-up agent-based models in particular will become an increasingly important platform technology for basic scientific discovery, as well as the identification and validation of potentially novel therapeutic targets.

Identifying aberrant pathways through integrated analysis of knowledge in pharmacogenomics

MOTIVATION

Many complex diseases are the result of abnormal pathway functions instead of single abnormalities. Disease diagnosis and intervention strategies must target these pathways while minimizing the interference with normal physiological processes. Large-scale identification of disease pathways and chemicals that may be used to perturb them requires the integration of information about drugs, genes, diseases and pathways. This information is currently distributed over several pharmacogenomics databases. An integrated analysis of the information in these databases can reveal disease pathways and facilitate novel biomedical analyses.

RESULTS

We demonstrate how to integrate pharmacogenomics databases through integration of the biomedical ontologies that are used as meta-data in these databases. The additional background knowledge in these ontologies can then be used to enable novel analyses. We identify disease pathways using a novel multi-ontology enrichment analysis over the Human Disease Ontology, and we identify significant associations between chemicals and pathways using an enrichment analysis over a chemical ontology. The drug-pathway and disease-pathway associations are a valuable resource for research in disease and drug mechanisms and can be used to improve computational drug repurposing.

AVAILABILITY

http://pharmgkb-owl.googlecode.com

Genomics and bioinformatics of Parkinson's disease

Within the last two decades, genomics and bioinformatics have profoundly impacted our understanding of the molecular mechanisms of Parkinson's disease (PD). From the description of the first PD gene in 1997 until today, we have witnessed the emergence of new technologies that have revolutionized our concepts to identify genetic mechanisms implicated in human health and disease. Driven by the publication of the human genome sequence and followed by the description of detailed maps for common genetic variability, novel applications to rapidly scrutinize the entire genome in a systematic, cost-effective manner have become a reality. As a consequence, about 30 genetic loci have been unequivocally linked to the pathogenesis of PD highlighting essential molecular pathways underlying this common disorder. Herein we discuss how neurogenomics and bioinformatics are applied to dissect the nature of this complex disease with the overall aim of developing rational therapeutic interventions.