Search

The gene regulatory network (GRN) of human cells encodes mechanisms to ensure proper functioning. However, if this GRN is dysregulated, the cell may enter into a disease state such as cancer. Understanding the GRN as a system can therefore help identify novel mechanisms underlying disease, which can lead to new therapies. To deduce regulatory interactions relevant to cancer, we applied a recent computational inference framework to data from perturbation experiments in squamous carcinoma cell line A431. GRNs were inferred using several methods, and the false discovery rate was controlled by the NestBoot framework. We developed a novel approach to assess the predictiveness of inferred GRNs against validation data, despite the lack of a gold standard. The best GRN was significantly more predictive than the null model, both in cross-validated benchmarks and for an independent dataset of the same genes under a different perturbation design. The inferred GRN captures many known regulatory interactions central to cancer-relevant processes in addition to predicting many novel interactions, some of which were experimentally validated, thus providing mechanistic insights that are useful for future cancer research.

Motivation Inference of gene regulatory networks (GRNs) from perturbation data can give detailed mechanistic insights of a biological system. Many inference methods exist, but the resulting GRN is generally sensitive to the choice of method-specific parameters. Even though the inferred GRN is optimal given the parameters, many links may be wrong or missing if the data is not informative. To make GRN inference reliable, a method is needed to estimate the support of each predicted link as the method parameters are varied. Results To achieve this we have developed a method called nested bootstrapping, which applies a bootstrapping protocol to GRN inference, and by repeated bootstrap runs assesses the stability of the estimated support values. To translate bootstrap support values to false discovery rates we run the same pipeline with shuffled data as input. This provides a general method to control the false discovery rate of GRN inference that can be applied to any setting of inference parameters, noise level, or data properties. We evaluated nested bootstrapping on a simulated dataset spanning a range of such properties, using the LASSO, Least Squares, RNI, GENIE3 and CLR inference methods. An improved inference accuracy was observed in almost all situations. Nested bootstrapping was incorporated into the GeneSPIDER package, which was also used for generating the simulated networks and data, as well as running and analyzing the inferences. Availability and implementation https://bitbucket.org/sonnhammergrni/genespider/src/NB/%2B Methods/NestBoot.m, Funding Agencies|National Cheng Kung University; Ministry of Science and Technology in Taiwan [105-2218-E-006-016-MY2]

Background: The progression of colorectal cancer (CRC) involves recurrent amplifications/mutations in the epidermal growth factor receptor (EGFR) and downstream signal transducers of the Ras pathway, KRAS and BRAF. Whether genetic events predicted to result in increased and constitutive signaling indeed lead to enhanced biological activity is often unclear and, due to technical challenges, unexplored. Here, we investigated proliferative signaling in CRC using a highly sensitive method for protein detection. The aim of the study was to determine whether multiple changes in proliferative signaling in CRC could be combined and exploited as a "complex biomarker" for diagnostic purposes. Methods: We used robotized capillary isoelectric focusing as well as conventional immunoblotting for the comprehensive analysis of epidermal growth factor receptor signaling pathways converging on extracellular regulated kinase 1/2 (ERK1/2), AKT, phospholipase C gamma 1 (PLC gamma 1) and c-SRC in normal mucosa compared with CRC stage II and IV. Computational analyses were used to test different activity patterns for the analyzed signal transducers. Results: Signaling pathways implicated in cell proliferation were differently dysregulated in CRC and, unexpectedly, several were downregulated in disease. Thus, levels of activated ERK1 (pERK1), but not pERK2, decreased in stage II and IV while total ERK1/2 expression remained unaffected. In addition, c-SRC expression was lower in CRC compared with normal tissues and phosphorylation on the activating residue Y418 was not detected. In contrast, PLC gamma 1 and AKT expression levels were elevated in disease. Immunoblotting of the different signal transducers, run in parallel to capillary isoelectric focusing, showed higher variability and lower sensitivity and resolution. Computational analyses showed that, while individual signaling changes lacked predictive power, using the combination of changes in three signaling components to create a "com

Gene regulatory network inference (that is, determination of the regulatory interactions between a set of genes) provides mechanistic insights of central importance to research in systems biology. Most contemporary network inference methods rely on a sparsity/regularization coefficient, which we call zeta (zeta), to determine the degree of sparsity of the network estimates, that is, the total number of links between the nodes. However, they offer little or no advice on how to select this sparsity coefficient, in particular, for biological data with few samples. We show that an empty network is more accurate than estimates obtained for a poor choice of zeta. In order to avoid such poor choices, we propose a method for optimization of zeta, which maximizes the accuracy of the inferred network for any sparsity-dependent inference method and data set. Our procedure is based on leave-one-out cross-optimization and selection of the zeta value that minimizes the prediction error. We also illustrate the adverse effects of noise, few samples, and uninformative experiments on network inference as well as our method for optimization of zeta. We demonstrate that our zeta optimization method for two widely used inference algorithms-Glmnet and NIR-gives accurate and informative estimates of the network structure, given that the data is informative enough., QC 20130627

DNA copy number aberrations (CNAs) are a hallmark of cancer genomes. However, little is known about how such changes affect global gene expression. We develop a modeling framework, EPoC (Endogenous Perturbation analysis of Cancer), to (1) detect disease-driving CNAs and their effect on target mRNA expression, and to (2) stratify cancer patients into long-and short-term survivors. Our method constructs causal network models of gene expression by combining genome-wide DNA-and RNA-level data. Prognostic scores are obtained from a singular value decomposition of the networks. By applying EPoC to glioblastoma data from The Cancer Genome Atlas consortium, we demonstrate that the resulting network models contain known disease-relevant hub genes, reveal interesting candidate hubs, and uncover predictors of patient survival. Targeted validations in four glioblastoma cell lines support selected predictions, and implicate the p53-interacting protein Necdin in suppressing glioblastoma cell growth. We conclude that large-scale network modeling of the effects of CNAs on gene expression may provide insights into the biology of human cancer. Free software in MATLAB and R is provided., QC 20110531