Identification of pathway deregulation--gene expression based analysis of consistent signal transduction

Transcriptomic Profiling for the Autophagy Pathway in Colorectal Cancer

The role of autophagy in colorectal cancer (CRC) pathogenesis appears to be crucial. Autophagy acts both as a tumor suppressor, by removing redundant cellular material, and a tumor-promoting factor, by providing access to components necessary for growth, metabolism, and proliferation. To date, little is known about the expression of genes that play a basal role in the autophagy in CRC. In this study, we aimed to compare the expression levels of 46 genes involved in the autophagy pathway between tumor-adjacent and tumor tissue, employing large RNA sequencing (RNA-seq) and microarray datasets. Additionally, we verified our results using data on 38 CRC cell lines. Gene set enrichment analysis revealed a significant deregulation of autophagy-related gene sets in CRC. The unsupervised clustering of tumors using the mRNA levels of autophagy-related genes revealed the existence of two major clusters: microsatellite instability (MSI)-enriched and -depleted. In cluster 1 (MSI-depleted), ATG9B and LAMP1 genes were the most prominently expressed, whereas cluster 2 (MSI-enriched) was characterized by DRAM1 upregulation. CRC cell lines were also clustered according to MSI-enriched/-depleted subgroups. The moderate deregulation of autophagy-related genes in cancer tissue, as compared to adjacent tissue, suggests a prominent field cancerization or early disruption of autophagy. Genes differentiating these clusters are promising candidates for CRC targeting therapy worthy of further investigation.

Search for novel STAT3-dependent genes reveals SERPINA3 as a new STAT3 target that regulates invasion of human melanoma cells

Transcription factor signal transducer and activator of transcription 3 (STAT3) is constitutively activated in many cancers and promotes uncontrolled tumor growth and progression through multiple mechanisms. Compelling evidence shows tissue and cell-specific sets of STAT3 targets. Transcriptional targets of STAT3 in melanoma cells are largely unknown. Malignant melanoma is a deadly disease with highly aggressive and drug-resistant behavior. Less than 10% of patients with advanced melanomas reach the 5-year survival, partly due to the aggressive character of the tumor and ineffectiveness of current therapeutics for treating metastatic melanoma. STAT3 is constitutively activated in melanoma cells and plays important roles in its growth and angiogenesis in tumor xenograft studies. Moreover, highly metastatic melanoma cells have higher levels of active STAT3 than poorly metastatic ones. To identify genes that are driven by STAT3 in human melanoma cells, we performed JAK/STAT signaling specific and global gene expression profiling of human melanoma cells with silenced STAT3 expression. For selected genes, we performed computational identification of putative STAT3-binding sites and validated direct interactions STAT3 with defined promoters by using chromatin immunoprecipitation followed by qPCR. We found that STAT3 knockdown does not affect human melanoma cell viability, proliferation, or response to chemotherapeutics. We show that STAT3 regulates a discrete set of genes in melanoma cells, including SERPINA3, a novel STAT3 target gene, which is functionally involved in regulation of melanoma migration and invasion. Knockdown of STAT3 impaired cell migration and invasion, in part via regulation of its transcriptional target SERPINA3. Our results present novel targets and functions of STAT3 in melanoma cells.

Network-Based Approaches for Pathway Level Analysis

Identification of impacted pathways is an important problem because it allows us to gain insights into the underlying biology beyond the detection of differentially expressed genes. In the past decade, a plethora of methods have been developed for this purpose. The last generation of pathway analysis methods are designed to take into account various aspects of pathway topology in order to increase the accuracy of the findings. Here, we cover 34 such topology-based pathway analysis methods published in the past 13 years. We compare these methods on categories related to implementation, availability, input format, graph models, and statistical approaches used to compute pathway level statistics and statistical significance. We also discuss a number of critical challenges that need to be addressed, arising both in methodology and pathway representation, including inconsistent terminology, data format, lack of meaningful benchmarks, and, more importantly, a systematic bias that is present in most existing methods. © 2018 by John Wiley & Sons, Inc.

Injury-induced perivascular niche supports alternative differentiation of adult rodent CNS progenitor cells

Following CNS demyelination, oligodendrocyte progenitor cells (OPCs) are able to differentiate into either remyelinating oligodendrocytes (OLs) or remyelinating Schwann cells (SCs). However, the signals that determine which type of remyelinating cell is generated and the underlying mechanisms involved have not been identified. Here, we show that distinctive microenvironments created in discrete niches within demyelinated white matter determine fate decisions of adult OPCs. By comparative transcriptome profiling we demonstrate that an ectopic, injury-induced perivascular niche is enriched with secreted ligands of the BMP and Wnt signalling pathways, produced by activated OPCs and endothelium, whereas reactive astrocyte within non-vascular area express the dual BMP/Wnt antagonist Sostdc1. The balance of BMP/Wnt signalling network is instructive for OPCs to undertake fate decision shortly after their activation: disruption of the OPCs homeostasis during demyelination results in BMP4 upregulation, which, in the absence of Socstdc1, favours SCs differentiation.

Identification of perturbed signaling pathways from gene expression data using information divergence

Abnormal regulation of signaling pathways is the key causative factor in several diseases. Although many methods have been proposed to identify significantly differential pathways between two conditions via microarray gene expression datasets, most of them concentrate on differences in the pathway components-either the differential expression or the correlation of genes in a given pathway. However, as biological functional units, signaling pathways may have diverse activity patterns across different biological contexts. In order to detect overall changes in pathways, we propose an analysis model called SPAID (Signaling Pathway Analysis based on Information Divergence). SPAID is based on the concept of information divergence, which can be used to compare two conditions by computing the differential probability distribution of the regulation capacity. We compared SPAID with several classical algorithms using different datasets, and the results indicate that SPAID produces higher repeatability, has better performance and universality, and extracts more comprehensive information regarding the underlying biological processes. In conclusion, by introducing the idea of information divergence, our study measures differences in pathways from an overall perspective and will provide a complementary analysis framework for pathway analysis.

Prioritizing biological pathways by recognizing context in time-series gene expression data

Background

The primary goal of pathway analysis using transcriptome data is to find significantly perturbed pathways. However, pathway analysis is not always successful in identifying pathways that are truly relevant to the context under study. A major reason for this difficulty is that a single gene is involved in multiple pathways. In the KEGG pathway database, there are 146 genes, each of which is involved in more than 20 pathways. Thus activation of even a single gene will result in activation of many pathways. This complex relationship often makes the pathway analysis very difficult. While we need much more powerful pathway analysis methods, a readily available alternative way is to incorporate the literature information.

Results

In this study, we propose a novel approach for prioritizing pathways by combining results from both pathway analysis tools and literature information. The basic idea is as follows. Whenever there are enough articles that provide evidence on which pathways are relevant to the context, we can be assured that the pathways are indeed related to the context, which is termed as relevance in this paper. However, if there are few or no articles reported, then we should rely on the results from the pathway analysis tools, which is termed as significance in this paper. We realized this concept as an algorithm by introducing Context Score and Impact Score and then combining the two into a single score. Our method ranked truly relevant pathways significantly higher than existing pathway analysis tools in experiments with two data sets.

Conclusions

Our novel framework was implemented as ContextTRAP by utilizing two existing tools, TRAP and BEST. ContextTRAP will be a useful tool for the pathway based analysis of gene expression data since the user can specify the context of the biological experiment in a set of keywords. The web version of ContextTRAP is available at http://biohealth.snu.ac.kr/software/contextTRAP.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1335-8) contains supplementary material, which is available to authorized users.

Application of Systems Biology to Neuroproteomics: The Path to Enhanced Theranostics in Traumatic Brain Injury

The application of systems biology tools in analyzing heterogeneous data from multiple sources has become a necessity, especially in biomarker discovery. Such tools were developed with several approaches to address different types of research questions and hypotheses. In the field of neurotrauma and traumatic brain injury (TBI), three distinct approaches have been used so far as systems biology tools, namely functional group categorization, pathway analysis, and protein-protein interaction (PPI) networks. The databases allow for query of the system to identify candidate targets which can be further studied to elucidate potential downstream biomarkers indicative of disease progression, severity, and improvement. The various systems biology tools, databases, and strategies that can be implemented on available TBI data in neuroproteomic studies are discussed in this chapter.

A MATLAB tool for pathway enrichment using a topology-based pathway regulation score

BACKGROUND

Handling the vast amount of gene expression data generated by genome-wide transcriptional profiling techniques is a challenging task, demanding an informed combination of pre-processing, filtering and analysis methods if meaningful biological conclusions are to be drawn. For example, a range of traditional statistical and computational pathway analysis approaches have been used to identify over-represented processes in microarray data derived from various disease states. However, most of these approaches tend not to exploit the full spectrum of gene expression data, or the various relationships and dependencies. Previously, we described a pathway enrichment analysis tool created in MATLAB that yields a Pathway Regulation Score (PRS) by considering signalling pathway topology, and the overrepresentation and magnitude of differentially-expressed genes (J Comput Biol 19:563-573, 2012). Herein, we extended this approach to include metabolic pathways, and described the use of a graphical user interface (GUI).

RESULTS

Using input from a variety of microarray platforms and species, users are able to calculate PRS scores, along with a corresponding z-score for comparison. Further pathway significance assessment may be performed to increase confidence in the pathways obtained, and users can view Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway diagrams marked-up to highlight impacted genes.

CONCLUSIONS

The PRS tool provides a filter in the isolation of biologically-relevant insights from complex transcriptomic data.

Multivariate inference of pathway activity in host immunity and response to therapeutics

Developing a quantitative view of how biological pathways are regulated in response to environmental factors is central for understanding of disease phenotypes. We present a computational framework, named Multivariate Inference of Pathway Activity (MIPA), which quantifies degree of activity induced in a biological pathway by computing five distinct measures from transcriptomic profiles of its member genes. Statistical significance of inferred activity is examined using multiple independent self-contained tests followed by a competitive analysis. The method incorporates a new algorithm to identify a subset of genes that may regulate the extent of activity induced in a pathway. We present an in-depth evaluation of specificity, robustness, and reproducibility of our method. We benchmarked MIPA's false positive rate at less than 1%. Using transcriptomic profiles representing distinct physiological and disease states, we illustrate applicability of our method in (i) identifying gene–gene interactions in autophagy-dependent response to Salmonella infection, (ii) uncovering gene–environment interactions in host response to bacterial and viral pathogens and (iii) identifying driver genes and processes that contribute to wound healing and response to anti-TNFα therapy. We provide relevant experimental validation that corroborates the accuracy and advantage of our method.

Bayesian pathway analysis of cancer microarray data

High Throughput Biological Data (HTBD) requires detailed analysis methods and from a life science perspective, these analysis results make most sense when interpreted within the context of biological pathways. Bayesian Networks (BNs) capture both linear and nonlinear interactions and handle stochastic events in a probabilistic framework accounting for noise making them viable candidates for HTBD analysis. We have recently proposed an approach, called Bayesian Pathway Analysis (BPA), for analyzing HTBD using BNs in which known biological pathways are modeled as BNs and pathways that best explain the given HTBD are found. BPA uses the fold change information to obtain an input matrix to score each pathway modeled as a BN. Scoring is achieved using the Bayesian-Dirichlet Equivalent method and significance is assessed by randomization via bootstrapping of the columns of the input matrix. In this study, we improve on the BPA system by optimizing the steps involved in "Data Preprocessing and Discretization", "Scoring", "Significance Assessment", and "Software and Web Application". We tested the improved system on synthetic data sets and achieved over 98% accuracy in identifying the active pathways. The overall approach was applied on real cancer microarray data sets in order to investigate the pathways that are commonly active in different cancer types. We compared our findings on the real data sets with a relevant approach called the Signaling Pathway Impact Analysis (SPIA).

Stratification and delineation of gastric cancer signaling by in vitro transcription factor activity profiling and integrative genomics

Integrative functional genomic approaches are helpful in delineating the complex dysregulations in cancers. In the present study, in vitro activity profiling of 45 signaling pathway driven transcription factors in eight gastric cancer cell lines and direct comparison with genome-wide profiles of gastric tumors were performed and the integration resulted in the identification of three categories of factors/pathways: i) highly activated signaling pathways that stem from mutations are the critical oncogenic drivers, ii) constitutively activated stress responsive pathways which are activated not due to genetic alterations, and iii) consistently down-regulated nuclear receptor responsive factors. This functional profiling helps in discriminating therapeutic targets and signaling interactions.

Methods and approaches in the topology-based analysis of biological pathways

The goal of pathway analysis is to identify the pathways significantly impacted in a given phenotype. Many current methods are based on algorithms that consider pathways as simple gene lists, dramatically under-utilizing the knowledge that such pathways are meant to capture. During the past few years, a plethora of methods claiming to incorporate various aspects of the pathway topology have been proposed. These topology-based methods, sometimes referred to as “third generation,” have the potential to better model the phenomena described by pathways. Although there is now a large variety of approaches used for this purpose, no review is currently available to offer guidance for potential users and developers. This review covers 22 such topology-based pathway analysis methods published in the last decade. We compare these methods based on: type of pathways analyzed (e.g., signaling or metabolic), input (subset of genes, all genes, fold changes, gene p-values, etc.), mathematical models, pathway scoring approaches, output (one or more pathway scores, p-values, etc.) and implementation (web-based, standalone, etc.). We identify and discuss challenges, arising both in methodology and in pathway representation, including inconsistent terminology, different data formats, lack of meaningful benchmarks, and the lack of tissue and condition specificity.