Gene Expression Integration into Pathway Modules Reveals a Pan-Cancer Metabolic Landscape

BIO2014-57291-R and SAF2017-88908-R from the Spanish Ministry of Economy and Competitivenessgrant PI15/00854 from the FIS“Plataforma de Recursos Biomoleculares y Bioinformáticos” PT17/0009/0006 from the ISCIII, cofunded with European Regional Development FundsFP7-PEOPLE-2012-ITN MLPM2012EU H2020-INFRADEV-1-2015-1 ELIXIR-EXCELERAT

A dynamic network approach for the study of human phenotypes

The use of networks to integrate different genetic, proteomic, and metabolic datasets has been proposed as a viable path toward elucidating the origins of specific diseases. Here we introduce a new phenotypic database summarizing correlations obtained from the disease history of more than 30 million patients in a Phenotypic Disease Network (PDN). We present evidence that the structure of the PDN is relevant to the understanding of illness progression by showing that (1) patients develop diseases close in the network to those they already have; (2) the progression of disease along the links of the network is different for patients of different genders and ethnicities; (3) patients diagnosed with diseases which are more highly connected in the PDN tend to die sooner than those affected by less connected diseases; and (4) diseases that tend to be preceded by others in the PDN tend to be more connected than diseases that precede other illnesses, and are associated with higher degrees of mortality. Our findings show that disease progression can be represented and studied using network methods, offering the potential to enhance our understanding of the origin and evolution of human diseases. The dataset introduced here, released concurrently with this publication, represents the largest relational phenotypic resource publicly available to the research community.Comment: 28 pages (double space), 6 figure

Increased entropy of signal transduction in the cancer metastasis phenotype

Studies into the statistical properties of biological networks have led to important biological insights, such as the presence of hubs and hierarchical modularity. There is also a growing interest in studying the statistical properties of networks in the context of cancer genomics. However, relatively little is known as to what network features differ between the cancer and normal cell physiologies, or between different cancer cell phenotypes. Based on the observation that frequent genomic alterations underlie a more aggressive cancer phenotype, we asked if such an effect could be detectable as an increase in the randomness of local gene expression patterns. Using a breast cancer gene expression data set and a model network of protein interactions we derive constrained weighted networks defined by a stochastic information flux matrix reflecting expression correlations between interacting proteins. Based on this stochastic matrix we propose and compute an entropy measure that quantifies the degree of randomness in the local pattern of information flux around single genes. By comparing the local entropies in the non-metastatic versus metastatic breast cancer networks, we here show that breast cancers that metastasize are characterised by a small yet significant increase in the degree of randomness of local expression patterns. We validate this result in three additional breast cancer expression data sets and demonstrate that local entropy better characterises the metastatic phenotype than other non-entropy based measures. We show that increases in entropy can be used to identify genes and signalling pathways implicated in breast cancer metastasis. Further exploration of such integrated cancer expression and protein interaction networks will therefore be a fruitful endeavour.Comment: 5 figures, 2 Supplementary Figures and Table

Predicting glioblastoma prognosis networks using weighted gene co-expression network analysis on TCGA data

<p>Abstract</p> <p>Background</p> <p>Using gene co-expression analysis, researchers were able to predict clusters of genes with consistent functions that are relevant to cancer development and prognosis. We applied a weighted gene co-expression network (WGCN) analysis algorithm on glioblastoma multiforme (GBM) data obtained from the TCGA project and predicted a set of gene co-expression networks which are related to GBM prognosis.</p> <p>Methods</p> <p>We modified the Quasi-Clique Merger algorithm (QCM algorithm) into edge-covering Quasi-Clique Merger algorithm (eQCM) for mining weighted sub-network in WGCN. Each sub-network is considered a set of features to separate patients into two groups using K-means algorithm. Survival times of the two groups are compared using log-rank test and Kaplan-Meier curves. Simulations using random sets of genes are carried out to determine the thresholds for log-rank test p-values for network selection. Sub-networks with p-values less than their corresponding thresholds were further merged into clusters based on overlap ratios (>50%). The functions for each cluster are analyzed using gene ontology enrichment analysis.</p> <p>Results</p> <p>Using the eQCM algorithm, we identified 8,124 sub-networks in the WGCN, out of which 170 sub-networks show p-values less than their corresponding thresholds. They were then merged into 16 clusters.</p> <p>Conclusions</p> <p>We identified 16 gene clusters associated with GBM prognosis using the eQCM algorithm. Our results not only confirmed previous findings including the importance of cell cycle and immune response in GBM, but also suggested important epigenetic events in GBM development and prognosis.</p

On dynamic network entropy in cancer

The cellular phenotype is described by a complex network of molecular interactions. Elucidating network properties that distinguish disease from the healthy cellular state is therefore of critical importance for gaining systems-level insights into disease mechanisms and ultimately for developing improved therapies. By integrating gene expression data with a protein interaction network to induce a stochastic dynamics on the network, we here demonstrate that cancer cells are characterised by an increase in the dynamic network entropy, compared to cells of normal physiology. Using a fundamental relation between the macroscopic resilience of a dynamical system and the uncertainty (entropy) in the underlying microscopic processes, we argue that cancer cells will be more robust to random gene perturbations. In addition, we formally demonstrate that gene expression differences between normal and cancer tissue are anticorrelated with local dynamic entropy changes, thus providing a systemic link between gene expression changes at the nodes and their local network dynamics. In particular, we also find that genes which drive cell-proliferation in cancer cells and which often encode oncogenes are associated with reductions in the dynamic network entropy. In summary, our results support the view that the observed increased robustness of cancer cells to perturbation and therapy may be due to an increase in the dynamic network entropy that allows cells to adapt to the new cellular stresses. Conversely, genes that exhibit local flux entropy decreases in cancer may render cancer cells more susceptible to targeted intervention and may therefore represent promising drug targets.Comment: 10 pages, 3 figures, 4 tables. Submitte

ProteoLens: a visual analytic tool for multi-scale database-driven biological network data mining

Background New systems biology studies require researchers to understand how interplay among myriads of biomolecular entities is orchestrated in order to achieve high-level cellular and physiological functions. Many software tools have been developed in the past decade to help researchers visually navigate large networks of biomolecular interactions with built-in template-based query capabilities. To further advance researchers' ability to interrogate global physiological states of cells through multi-scale visual network explorations, new visualization software tools still need to be developed to empower the analysis. A robust visual data analysis platform driven by database management systems to perform bi-directional data processing-to-visualizations with declarative querying capabilities is needed. Results We developed ProteoLens as a JAVA-based visual analytic software tool for creating, annotating and exploring multi-scale biological networks. It supports direct database connectivity to either Oracle or PostgreSQL database tables/views, on which SQL statements using both Data Definition Languages (DDL) and Data Manipulation languages (DML) may be specified. The robust query languages embedded directly within the visualization software help users to bring their network data into a visualization context for annotation and exploration. ProteoLens supports graph/network represented data in standard Graph Modeling Language (GML) formats, and this enables interoperation with a wide range of other visual layout tools. The architectural design of ProteoLens enables the de-coupling of complex network data visualization tasks into two distinct phases: 1) creating network data association rules, which are mapping rules between network node IDs or edge IDs and data attributes such as functional annotations, expression levels, scores, synonyms, descriptions etc; 2) applying network data association rules to build the network and perform the visual annotation of graph nodes and edges according to associated data values. We demonstrated the advantages of these new capabilities through three biological network visualization case studies: human disease association network, drug-target interaction network and protein-peptide mapping network. Conclusion The architectural design of ProteoLens makes it suitable for bioinformatics expert data analysts who are experienced with relational database management to perform large-scale integrated network visual explorations. ProteoLens is a promising visual analytic platform that will facilitate knowledge discoveries in future network and systems biology studies

The Evolving Transcriptome of Head and Neck Squamous Cell Carcinoma: A Systematic Review

BACKGROUND: Numerous studies were performed to illuminate mechanisms of tumorigenesis and metastases from gene expression profiles of Head and Neck Squamous Cell Carcinoma (HNSCC). The objective of this review is to conduct a network-based meta-analysis to identify the underlying biological signatures of the HNSCC transcriptome. METHODS AND FINDINGS: We included 63 HNSCC transcriptomic studies into three specific categories of comparisons: Pre, premalignant lesions v.s. normal; TvN, primary tumors v.s. normal; and Meta, metastatic or invasive v.s. primary tumors. Reported genes extracted from the literature were systematically analyzed. Participation of differential gene activities across three progressive stages deciphered the evolving nature of HNSCC. In total, 1442 genes were verified, i.e. reported at least twice, with ECM1, EMP1, CXCL10 and POSTN shown to be highly reported across all three stages. Knowledge-based networks of the HNSCC transcriptome were constructed, demonstrating integrin signaling and antigen presentation pathways as highly enriched. Notably, functional estimates derived from topological characteristics of integrin signaling networks identified such important genes as ITGA3 and ITGA5, which were supported by findings of invasiveness in vitro. Moreover, we computed genome-wide probabilities of reporting differential gene activities for the Pre, TvN, and Meta stages, respectively. Results highlighted chromosomal regions of 6p21, 19p13 and 19q13, where genomic alterations were shown to be correlated with the nodal status of HNSCC. CONCLUSIONS: By means of a systems-biology approach via network-based meta-analyses, we provided a deeper insight into the evolving nature of the HNSCC transcriptome. Enriched canonical signaling pathways, hot-spots of transcriptional profiles across the genome, as well as topologically significant genes derived from network analyses were highlighted for each of the three progressive stages, Pre, TvN, and Meta, respectively

Dynamic Regulation of Oct1 during Mitosis by Phosphorylation and Ubiquitination

Transcription factor Oct1 regulates multiple cellular processes. It is known to be phosphorylated during the cell cycle and by stress, however the upstream kinases and downstream consequences are not well understood. One of these modified forms, phosphorylated at S335, lacks the ability to bind DNA. Other modification states besides phosphorylation have not been described.We show that Oct1 is phosphorylated at S335 in the Oct1 DNA binding domain during M-phase by the NIMA-related kinase Nek6. Phospho-Oct1 is also ubiquitinated. Phosphorylation excludes Oct1 from mitotic chromatin. Instead, Oct1(pS335) concentrates at centrosomes, mitotic spindle poles, kinetochores and the midbody. Oct1 siRNA knockdown diminishes the signal at these locations. Both Oct1 ablation and overexpression result in abnormal mitoses. S335 is important for the overexpression phenotype, implicating this residue in mitotic regulation. Oct1 depletion causes defects in spindle morphogenesis in Xenopus egg extracts, establishing a mitosis-specific function of Oct1. Oct1 colocalizes with lamin B1 at the spindle poles and midbody. At the midbody, both proteins are mutually required to correctly localize the other. We show that phospho-Oct1 is modified late in mitosis by non-canonical K11-linked polyubiquitin chains. Ubiquitination requires the anaphase-promoting complex, and we further show that the anaphase-promoting complex large subunit APC1 and Oct1(pS335) interact.These findings reveal mechanistic coupling between Oct1 phosphorylation and ubquitination during mitotic progression, and a role for Oct1 in mitosis