Comparative analysis of protein interactome networks prioritizes candidate genes with cancer signatures

Open a new epoch of arsenic trioxide investigation: ATOdb

Arsenic trioxide (ATO) is a great discovery in the treatment of acute promyelocytic leukemia (APL), which has been used in an increasing number of malignant diseases. Systematic integrative analysis will help to precisely understand the mechanism of ATO and find new combined drugs. Therefore, we developed a one-stop comprehensive database of ATO named ATOdb by manually compiling a wealth of experimentally supported ATO-related data from 3479 articles, and integrated analysis tools. The current version of ATOdb contains 8373 associations among 2300 ATO targets, 80 conditions and 262 combined drugs. Each entry in ATOdb contains detailed information on ATO targets, therapeutic/side effects, systems, cell names, cell types, regulations, detection methods, brief descriptions, references, etc. Furthermore, ATOdb also provides data visualization and analysis results such as the drug similarities, protein-protein interactions, and miRNA-mRNA relationships. An easy-to-use web interface was deployed in ATOdb for users to easily browse, search and download the data. In conclusion, ATOdb will serve as a valuable resource for in-depth study of the mechanism of ATO, discovery of new drug combination strategies, promotion of rational drug use and individualized treatments. ATOdb is freely available at http://bio-bigdata.hrbmu.edu.cn/ATOdb/index.jsp.

Polygenic and Network-Based Studies in Risk Identification and Demystification of cancer

INTRODUCTION

Diseases were initially thought to be the consequence of a single gene mutation. Advances in DNA sequencing tools and our understanding of gene behavior have revealed that complex diseases, such as cancer, are the product of genes cooperating with each other and with their environment in orchestrated communication networks. Seeing that the function of individual genes is still used to analyze cancer, the shift to using functionally interacting groups of genes as a new unit of study holds promise for demystifying cancer.

AREAS COVERED

The literature search focused on three types of cancer, namely breast, lung, and prostate, but arguments from other cancers were also included. The aim was to prove that multigene analyses can accurately predict and prognosticate cancer risk, subtype cancer for more personalized and effective treatments, and discover anti-cancer therapies. Computational intelligence is being harnessed to analyze this type of data and is proving indispensable to scientific progress.

EXPERT OPINION

In the future, comprehensive profiling of all kinds of patient data (e.g., serum molecules, environmental exposures) can be used to build universal networks that should help us elucidate the molecular mechanisms underlying diseases and provide appropriate preventive measures, ensuring lifelong health and longevity.

Comprehensive Interactome Mapping of Nuclear Receptors Using Proximity Biotinylation

Nuclear receptors, including hormone receptors, perform their cellular activities by modulating their protein-protein interactions. They engage with specific ligands and translocate to the nucleus, where they bind the DNA and activate extensive transcriptional programs. Therefore, gaining a comprehensive overview of the protein-protein interactions they establish requires methods that function effectively throughout the cell with fast dynamics and high reproducibility. Focusing on estrogen receptor alpha (ESR1), the founding member of the nuclear receptor family, this chapter describes a new lentiviral system that allows the expression of TurboID-hemagglutinin (HA)-2 × Strep tagged proteins in mammalian cells to perform fast proximity biotinylation assays. Key validation steps for these reagents and their use in interactome mapping experiments in two distinct breast cancer cell lines are described. Our protocol enabled the quantification of ESR1 interactome generated by cellular contexts that were hormone-sensitive or not.

A full-proteome, interaction-specific characterization of mutational hotspots across human cancers

Rapid accumulation of cancer genomic data has led to the identification of an increasing number of mutational hotspots with uncharacterized significance. Here we present a biologically informed computational framework that characterizes the functional relevance of all 1107 published mutational hotspots identified in approximately 25,000 tumor samples across 41 cancer types in the context of a human 3D interactome network, in which the interface of each interaction is mapped at residue resolution. Hotspots reside in network hub proteins and are enriched on protein interaction interfaces, suggesting that alteration of specific protein-protein interactions is critical for the oncogenicity of many hotspot mutations. Our framework enables, for the first time, systematic identification of specific protein interactions affected by hotspot mutations at the full proteome scale. Furthermore, by constructing a hotspot-affected network that connects all hotspot-affected interactions throughout the whole-human interactome, we uncover genome-wide relationships among hotspots and implicate novel cancer proteins that do not harbor hotspot mutations themselves. Moreover, applying our network-based framework to specific cancer types identifies clinically significant hotspots that can be used for prognosis and therapy targets. Overall, we show that our framework bridges the gap between the statistical significance of mutational hotspots and their biological and clinical significance in human cancers.

High expression of MCM10 is predictive of poor outcomes in lung adenocarcinoma

Backgrounds

Lung adenocarcinoma is a complex disease that results in over 1.8 million deaths a year. Recent advancements in treating and managing lung adenocarcinoma have led to modest decreases in associated mortality rates, owing in part to the multifactorial etiology of the disease. Novel prognostic biomarkers are needed to accurately stage the disease and act as the basis of adjuvant treatments.

Material and Methods

The microarray datasets GSE75037, GSE31210 and GSE32863 were downloaded from the Gene Expression Omnibus (GEO) database to identify prognostic biomarkers for lung adenocarcinoma and therapy. The differentially expressed genes (DEGs) were identified by GEO2R. Functional and pathway enrichment analysis were performed by Kyoto Encyclopedia of Genes and Genomes and Gene Ontology (GO). Validation was performed based on 72 pairs of lung adenocarcinoma and adjacent normal lung tissues.

Results

Results showed that the DEGs were mainly focused on cell cycle and DNA replication initiation. Forty-one hub genes were identified and further analyzed by CytoScape. Here, we provide evidence which suggests MCM10 is a potential target with prognostic, diagnostic and therapeutic value. We base this on an integrated approach of comprehensive bioinformatics analysis and in vitro validation using the A549 lung adenocarcinoma cell line. We show that MCM10 overexpression correlates with a poor prognosis, while silencing of this gene decreases aberrant growth by 2-fold. Finally, evaluation of 72 clinical biopsy samples suggests that overexpression of MCM10 in the lung adenocarcinoma highly correlates with larger tumor size. Together, this work suggests that MCM10 may be a clinically relevant gene with both predictive and therapeutic value in lung adenocarcinoma.

e-MutPath: computational modeling reveals the functional landscape of genetic mutations rewiring interactome networks

Understanding the functional impact of cancer somatic mutations represents a critical knowledge gap for implementing precision oncology. It has been increasingly appreciated that the interaction profile mediated by a genomic mutation provides a fundamental link between genotype and phenotype. However, specific effects on biological signaling networks for the majority of mutations are largely unknown by experimental approaches. To resolve this challenge, we developed e-MutPath (edgetic Mutation-mediated Pathway perturbations), a network-based computational method to identify candidate ‘edgetic’ mutations that perturb functional pathways. e-MutPath identifies informative paths that could be used to distinguish disease risk factors from neutral elements and to stratify disease subtypes with clinical relevance. The predicted targets are enriched in cancer vulnerability genes, known drug targets but depleted for proteins associated with side effects, demonstrating the power of network-based strategies to investigate the functional impact and perturbation profiles of genomic mutations. Together, e-MutPath represents a robust computational tool to systematically assign functions to genetic mutations, especially in the context of their specific pathway perturbation effect.

Network-based integration method for potential breast cancer gene identification

Breast cancer is the most common female death-causing cancer worldwide. A network-based integration method was proposed to identify potential breast cancer genes. First, genes were prioritized using a gene prioritization algorithm by the strategy of disease risks transferred between genes in a network with weighted vertexes and edges. Our prioritization algorithm was effectives and robust for top-ranked seed gene number and higher area under the curve values compared to ToppGene and ToppNet. Then, 20 potential breast cancer genes were identified as common genes of the top 50 candidate genes for their robustness in multiple prioritizations. These genes could accurately classify tumor and normal samples of all and paired sample sets and three independent datasets. Of potential breast cancer genes, 18 were verified by literature and 2 were novel genes that need further study. This study would contribute to the understanding of the genetic architecture for the diagnosis and treatment of breast cancer.

Context-specific interactions in literature-curated protein interaction databases

BACKGROUND

Databases of literature-curated protein-protein interactions (PPIs) are often used to interpret high-throughput interactome mapping studies and estimate error rates. These databases combine interactions across thousands of published studies and experimental techniques. Because the tendency for two proteins to interact depends on the local conditions, this heterogeneity of conditions means that only a subset of database PPIs are interacting during any given experiment. A typical use of these databases as gold standards in interactome mapping projects, however, assumes that PPIs included in the database are indeed interacting under the experimental conditions of the study.

RESULTS

Using raw data from 20 co-fractionation experiments and six published interactomes, we demonstrate that this assumption is often false, with up to 55% of purported gold standard interactions showing no evidence of interaction, on average. We identify a subset of CORUM database complexes that do show consistent evidence of interaction in co-fractionation studies, and we use this subset as gold standards to dramatically improve interactome mapping as judged by the number of predicted interactions at a given error rate.

CONCLUSIONS

We recommend using this CORUM subset as the gold standard set in future co-fractionation studies. More generally, we recommend using the subset of literature-curated PPIs that are specific to the experimental context whenever possible.

Integrative gene ontology and network analysis of coronary artery disease associated genes suggests potential role of ErbB pathway gene EGFR

Coronary artery disease (CAD) is a major cause of mortality in India, more importantly the young Indians. Combinatorial and integrative approaches to evaluate pathways and genes to gain an improved understanding and potential biomarkers for risk assessment are required. Therefore, 608 genes from the CADgene database version 2.0, classified into 12 functional classes representing the atherosclerotic disease process, were analyzed. Homology analysis of the unique list of gene ontologies (GO) from each functional class gave 8 GO terms represented in 11 and 10 functional classes. Using disease ontology analysis 80 genes belonging to 8 GO terms, using FunDO suggested that 29 of them were identified to be associated with CAD. Extended network analysis of these genes using STRING version 9.1 gave 328 nodes and 4,525 interactions of which the top 5% had a node degree of ≥75 associated with pathways including the ErbB signaling pathway with epidermal growth factor receptor (EGFR) gene as the central hub. Evaluation of EFGR protein levels in age and gender-matched 342 CAD patients vs. 342 control subjects demonstrated significant differences [controls=149.76±2.47 pg/ml and CAD patients stratified into stable angina (SA)=161.65±3.40 pg/ml and myocardial infarction (MI)=171.51±4.26 pg/ml]. Logistic regression analysis suggested that increased EGFR levels exhibit 3-fold higher risk of CAD [odds ratio (OR) 3.51, 95% confidence interval [CI] 1.96–6.28, P≤0.001], upon adjustment for hypertension, diabetes and smoking. A unit increase in EGFR levels increased the risk by 2-fold for SA (OR 2.58, 95% CI 1.25–5.33, P=0.01) and 3.8-fold for MI (OR 3.82, 95% CI 1.94–7.52, P≤0.001) following adjustment. Thus, the use of ontology mapping and network analysis in an integrative manner aids in the prioritization of biomarkers of complex disease.

Base-resolution stratification of cancer mutations using functional variomics

A complete understanding of human cancer variants requires new methods to systematically and efficiently assess the functional effects of genomic mutations at a large scale. Here, we describe a set of tools to rapidly clone and stratify thousands of cancer mutations at base resolution. This protocol provides a massively parallel pipeline to achieve high stringency and throughput. The approach includes high-throughput generation of mutant clones by Gateway, confirmation of variant identity by barcoding and next-generation sequencing, and stratification of cancer variants by multiplexed interaction profiling. Compared with alternative site-directed mutagenesis methods, our protocol requires less sequencing effort and enables robust statistical calling of allele-specific effects. To ensure the precision of variant interaction profiling, we further describe two complementary methods-a high-throughput enhanced yeast two-hybrid (HT-eY2H) assay and a mammalian-cell-based Gaussia princeps luciferase protein-fragment complementation assay (GPCA). These independent assays with standard controls validate mutational interaction profiles with high quality. This protocol provides experimentally derived guidelines for classifying candidate cancer alleles emerging from whole-genome or whole-exome sequencing projects as 'drivers' or 'passengers'. For ∼100 genomic mutations, the protocol-including target primer design, variant library construction, and sequence verification-can be completed within as little as 2-3 weeks, and cancer variant stratification can be completed within 2 weeks.