Dynamic modularity in protein interaction networks predicts breast cancer outcome

A large accessory protein interactome is rewired across environments

To characterize how protein-protein interaction (PPI) networks change, we quantified the relative PPI abundance of 1.6 million protein pairs in the yeast Saccharomyces cerevisiae across nine growth conditions, with replication, for a total of 44 million measurements. Our multi-condition screen identified 13,764 pairwise PPIs, a threefold increase over PPIs identified in one condition. A few 'immutable' PPIs are present across all conditions, while most 'mutable' PPIs are rarely observed. Immutable PPIs aggregate into highly connected 'core' network modules, with most network remodeling occurring within a loosely connected 'accessory' module. Mutable PPIs are less likely to co-express, co-localize, and be explained by simple mass action kinetics, and more likely to contain proteins with intrinsically disordered regions, implying that environment-dependent association and binding is critical to cellular adaptation. Our results show that protein interactomes are larger than previously thought and contain highly dynamic regions that reorganize to drive or respond to cellular changes.

Identification of key pathways and genes responsible for aggressive behavior

Aggression is a complex behavior, underpinned by cross talk between several biomolecules. To date a composite molecular network of the behavioral disorder has not been constructed. The present study aims to develop the same from the system network analyses recruiting genes with empirical evidence demonstrating their role in the incidence and progression of aggression. In short, 327 genes were recruited in the study after extensive literature survey and subsequent shortlisting by sieving out the comorbidities like cancer and other pathological and physiological ailments, other languages and repeated citations. Subsequent String network analysis coalesces 275 genes in a network with 2223 edges. The developed network was then subjected to delineate modules using MCODE which via gene clustering on the basis of gene ontology segregate all genes into 14 modules. Of these, as expected top 5 modules involved entailing of neuronal signaling pathways with redundant repetitions. Finally, 10 genes (known) were picked randomly, accounting average module size, and subjected to the network analysis with 100,000 bootstrap replicates. This results in the detection of certain novel genes that lacks empirical evidence for their association with the aggression. Amongst those, most notable are genes involved in protein turnover regulation like UBC, UBA, mitogenic proteins such as Rho and Myc, transcription factors like Tp53. The findings in turn fill caveats in the molecular resolution of cross talk that underscore the development of aggressive behavior and may then be exploited as screening biomarker and/or therapeutic intervention for aggression.

Integrative biomarker detection on high-dimensional gene expression data sets: a survey on prior knowledge approaches

Gene expression data provide the expression levels of tens of thousands of genes from several hundred samples. These data are analyzed to detect biomarkers that can be of prognostic or diagnostic use. Traditionally, biomarker detection for gene expression data is the task of gene selection. The vast number of genes is reduced to a few relevant ones that achieve the best performance for the respective use case. Traditional approaches select genes based on their statistical significance in the data set. This results in issues of robustness, redundancy and true biological relevance of the selected genes. Integrative analyses typically address these shortcomings by integrating multiple data artifacts from the same objects, e.g. gene expression and methylation data. When only gene expression data are available, integrative analyses instead use curated information on biological processes from public knowledge bases. With knowledge bases providing an ever-increasing amount of curated biological knowledge, such prior knowledge approaches become more powerful. This paper provides a thorough overview on the status quo of biomarker detection on gene expression data with prior biological knowledge. We discuss current shortcomings of traditional approaches, review recent external knowledge bases, provide a classification and qualitative comparison of existing prior knowledge approaches and discuss open challenges for this kind of gene selection.

Gene dysregulation analysis builds a mechanistic signature for prognosis and therapeutic benefit in colorectal cancer

The implementation of cancer precision medicine requires biomarkers or signatures for predicting prognosis and therapeutic benefits. Most of current efforts in this field are paying much more attention to predictive accuracy than to molecular mechanistic interpretability. Mechanism-driven strategy has recently emerged, aiming to build signatures with both predictive power and explanatory power. Driven by this strategy, we developed a robust gene dysregulation analysis framework with machine learning algorithms, which is capable of exploring gene dysregulations underlying carcinogenesis from high-dimensional data with cooperativity and synergy between regulators and several other transcriptional regulation rules taken into consideration. We then applied the framework to a colorectal cancer (CRC) cohort from The Cancer Genome Atlas. The identified CRC-related dysregulations significantly covered known carcinogenic processes and exhibited good prognostic effect. By choosing dysregulations with greedy strategy, we built a four-dysregulation (4-DysReg) signature, which has the capability of predicting prognosis and adjuvant chemotherapy benefit. 4-DysReg has the potential to explain carcinogenesis in terms of dysfunctional transcriptional regulation. These results demonstrate that our gene dysregulation analysis framework could be used to develop predictive signature with mechanistic interpretability for cancer precision medicine, and furthermore, elucidate the mechanisms of carcinogenesis.

(Differential) Co-Expression Analysis of Gene Expression: A Survey of Best Practices

Analysis of gene expression data is widely used in transcriptomic studies to understand functions of molecules inside a cell and interactions among molecules. Differential co-expression analysis studies diseases and phenotypic variations by finding modules of genes whose co-expression patterns vary across conditions. We review the best practices in gene expression data analysis in terms of analysis of (differential) co-expression, co-expression network, differential networking, and differential connectivity considering both microarray and RNA-seq data along with comparisons. We highlight hurdles in RNA-seq data analysis using methods developed for microarrays. We include discussion of necessary tools for gene expression analysis throughout the paper. In addition, we shed light on scRNA-seq data analysis by including preprocessing and scRNA-seq in co-expression analysis along with useful tools specific to scRNA-seq. To get insights, biological interpretation and functional profiling is included. Finally, we provide guidelines for the analyst, along with research issues and challenges which should be addressed.

An integrative framework of heterogeneous genomic data for cancer dynamic modules based on matrix decomposition

Cancer progression is dynamic, and tracking dynamic modules is promising for cancer diagnosis and therapy. Accumulated genomic data provide us an opportunity to investigate the underlying mechanisms of cancers. However, as far as we know, no algorithm has been designed for dynamic modules by integrating heterogeneous omics data. To address this issue, we propose an integrative framework for dynamic module detection based on regularized nonnegative matrix factorization method (DrNMF) by integrating the gene expression and protein interaction network. To remove the heterogeneity of genomic data, we divide the samples of expression profiles into groups to construct gene co-expression networks. To characterize the dynamics of modules, the temporal smoothness framework is adopted, in which the gene co-expression network at the previous stage and protein interaction network are incorporated into the objective function of DrNMF via regularization. The experimental results demonstrate that DrNMF is superior to state-of-the-art methods in terms of accuracy. For breast cancer data, the obtained dynamic modules are more enriched by the known pathways, and can be used to predict the stages of cancers and survival time of patients. The proposed model and algorithm provide an effective integrative analysis of heterogeneous genomic data for cancer progression.

Functional genomics identifies new synergistic therapies for retinoblastoma

Local intravitreal or intra-arterial chemotherapy has improved therapeutic success for the pediatric cancer retinoblastoma (RB), but toxicity remains a major caveat. RB initiates primarily with RB1 loss or, rarely, MYCN amplification, but the critical downstream networks are incompletely understood. We set out to uncover perturbed molecular hubs, identify synergistic drug combinations to target these vulnerabilities, and expose and overcome drug resistance. We applied dynamic transcriptomic analysis to identify network hubs perturbed in RB versus normal fetal retina, and performed in vivo RNAi screens in RB1^null and RB1^wt;MYCN^amp orthotopic xenografts to pinpoint essential hubs. We employed in vitro and in vivo studies to validate hits, define mechanism, develop new therapeutic modalities, and understand drug resistance. We identified BRCA1 and RAD51 as essential for RB cell survival. Their oncogenic activity was independent of BRCA1 functions in centrosome, heterochromatin, or ROS regulation, and instead linked to DNA repair. RAD51 depletion or inhibition with the small molecule inhibitor, B02, killed RB cells in a Chk1/Chk2/p53-dependent manner. B02 further synergized with clinically relevant topotecan (TPT) to engage this pathway, activating p53-BAX mediated killing of RB but not human retinal progenitor cells. Paradoxically, a B02/TPT-resistant tumor exhibited more DNA damage than sensitive RB cells. Resistance reflected dominance of the p53-p21 axis, which mediated cell cycle arrest instead of death. Deleting p21 or applying the BCL2/BCL2L1 inhibitor Navitoclax re-engaged the p53-BAX axis, and synergized with B02, TPT or both to override resistance. These data expose new synergistic therapies to trigger p53-induced killing in diverse RB subtypes.

Construction, comparison and evolution of networks in life sciences and other disciplines

Network approaches have become pervasive in many research fields. They allow for a more comprehensive understanding of complex relationships between entities as well as their group-level properties and dynamics. Many networks change over time, be it within seconds or millions of years, depending on the nature of the network. Our focus will be on comparative network analyses in life sciences, where deciphering temporal network changes is a core interest of molecular, ecological, neuropsychological and evolutionary biologists. Further, we will take a journey through different disciplines, such as social sciences, finance and computational gastronomy, to present commonalities and differences in how networks change and can be analysed. Finally, we envision how borrowing ideas from these disciplines could enrich the future of life science research.

Methods of Computational Interactomics for Investigating Interactions of Human Proteoforms

Human genome contains ca. 20,000 protein-coding genes that could be translated into millions of unique protein species (proteoforms). Proteoforms coded by a single gene often have different functions, which implies different protein partners. By interacting with each other, proteoforms create a network reflecting the dynamics of cellular processes in an organism. Perturbations of protein-protein interactions change the network topology, which often triggers pathological processes. Studying proteoforms is a relatively new research area in proteomics, and this is why there are comparatively few experimental studies on the interaction of proteoforms. Bioinformatics tools can facilitate such studies by providing valuable complementary information to the experimental data and, in particular, expanding the possibilities of the studies of proteoform interactions.

Comparative evaluation of network features for the prediction of breast cancer metastasis

Background

Discovering a highly accurate and robust gene signature for the prediction of breast cancer metastasis from gene expression profiling of primary tumors is one of the most challenging tasks to reduce the number of deaths in women. Due to the limited success of gene-based features in achieving satisfactory prediction accuracy, many methodologies have been proposed in recent years to develop network-based features by integrating network information with gene expression. However, evaluation results are inconsistent to confirm the effectiveness of network-based features, because of many confounding factors involved in classification model learning process, such as data normalization, dimension reduction, and feature selection. An unbiased comparative evaluation is essential for uncovering the strength of network-based features.

Methods

In this study, we compared several types of network-based features obtained using different mathematical operators (Mean, Maximum, Minimum, Median, Variance) on geneset (i.e., a gene and its’ neighbors in the network) in protein-protein interaction network and gene co-expression network for their ability in predicting breast cancer metastasis using gene expression data from more than 10 patient cohorts.

Results

While network-based features are usually statistically more significant than gene-based feature, a consistent improvement of prediction performance using network-based features requires a substantial number of patients in the dataset. In contrary to many previous reports, no evidence was found to support the robustness of network-based features and we argue some of the robustness may be due to the inherent bias associated with node degree in the network. In addition, different types of network features seem to cover different pathways and are complementary to each other. Consequently, an ensemble classifier combining different network features was proposed and was found to significantly outperform classifiers based on gene-based feature or any single type of network-based features.

Conclusions

Network-based features and their combination show promise for improving the prediction of breast cancer metastasis but may require a large amount of training data. Robustness claim of network-based features needs to be re-examined with network node degree and other confounding factors in consideration.

Identification of common key genes and pathways between type 1 diabetes and multiple sclerosis using transcriptome and interactome analysis

PURPOSE

Type 1 diabetes (T1D) and multiple sclerosis (MS) are classified as T cell-mediated autoimmune diseases. Although convergent evidence proposed common genetic architecture for autoimmune diseases, it remains a challenge to identify them. This study aimed to determine common gene signature and pathways in T1D and MS via systems biology approach.

METHODS

Gene expression profiles of peripheral blood mononuclear cells (PBMCs) and pancreatic-β cells in T1D as well as PBMCs and cerebrospinal fluid (CSF) in MS were analyzed in our previous published data, and differential expressed genes were integrated with protein-protein interactions data to construct Query-Query PPI (QQPPI) networks. In this study, QQPPI networks were further analyzed to investigate more central genes, functional modules and complexes shared in T1D and MS progression. Lastly, the interaction of common genes with drugs was also explored.

RESULTS

Several cytokines such as IL-23A, IL-32, IL-34, and IL-37 tend to be differentially expressed in both diseases. In addition, PSMA1, MYC, SRPK1, YBX1, HNRNPM, NF-κB2, IKBKE, RAC1, FN1, ARRB2, ESR1, HSP90AB1, and PPP1CA were common high central genes in QQPPI networks corresponding to each disease. Proteasome, spliceosome, immune responses, apoptosis, cellular communication/signaling transduction mechanism, interaction with environment, and activity of intercellular mediators were shared biological processes in T1D and MS. Finally, azathioprine, melatonin, resveratrol, and geldanamycin identified as prioritized drugs for the treatment of patients with T1D and MS.

CONCLUSIONS

This study represented novel key genes and pathways shared between T1D and MS, which may facilitate the identification of potential therapeutic targets in these diseases.

The self-organization model reveals systematic characteristics of aging

BACKGROUND

Aging is a fundamental biological process, where key bio-markers interact with each other and synergistically regulate the aging process. Thus aging dysfunction will induce many disorders. Finding aging markers and re-constructing networks based on multi-omics data (i.e. methylation, transcriptional and so on) are informative to study the aging process. However, optimizing the model to predict aging have not been performed systemically, although it is critical to identify potential molecular mechanism of aging related diseases.

METHODS

This paper aims to model the aging self-organization system using a series of supervised learning methods, and study complex molecular mechanisms of aging at system level: i.e. optimizing the aging network; summarizing interactions between aging markers; accumulating patterns of aging markers within module; finding order-parameters in the aging self-organization system.

RESULTS

In this work, the normal aging process is modeled based on multi-omics profiles across different tissues. In addition, the computational pipeline aims to model aging self-organizing systems and study the relationship between aging and related diseases (i.e. cancers), thus provide useful indicators of aging related diseases and could help to improve prediction abilities of diagnostics.

CONCLUSIONS

The aging process could be studied thoroughly by modelling the self-organization system, where key functions and the crosstalk between aging and cancers were identified.

A deep learning model based on sparse auto-encoder for prioritizing cancer-related genes and drug target combinations

Prioritization of cancer-related genes from gene expression profiles and proteomic data is vital to improve the targeted therapies research. Although computational approaches have been complementing high-throughput biological experiments on the understanding of human diseases, it still remains a big challenge to accurately discover cancer-related proteins/genes via automatic learning from large-scale protein/gene expression data and protein-protein interaction data. Most of the existing methods are based on network construction combined with gene expression profiles, which ignore the diversity between normal samples and disease cell lines. In this study, we introduced a deep learning model based on a sparse auto-encoder to learn the specific characteristics of protein interactions in cancer cell lines integrated with protein expression data. The model showed learning ability to identify cancer-related proteins/genes from the input of different protein expression profiles by extracting the characteristics of protein interaction information, which could also predict cancer-related protein combinations. Comparing with other reported methods including differential expression and network-based methods, our model got the highest area under the curve value (>0.8) in predicting cancer-related genes. Our study prioritized ~500 high-confidence cancer-related genes; among these genes, 211 already known cancer drug targets were found, which supported the accuracy of our method. The above results indicated that the proposed auto-encoder model could computationally prioritize candidate proteins/genes involved in cancer and improve the targeted therapies research.

Network Properties of Cancer Prognostic Gene Signatures in the Human Protein Interactome

Prognostic gene signatures are critical in cancer prognosis assessments and their pinpoint treatments. However, their network properties remain unclear. Here, we obtained nine prognostic gene sets including 1439 prognostic genes of different cancers from related publications. Four network centralities were used to examine the network properties of prognostic genes (PG) compared with other gene sets based on the Human Protein Reference Database (HPRD) and String networks. We also proposed three novel network measures for further investigating the network properties of prognostic gene sets (PGS) besides clustering coefficient. The results showed that PG did not occupy key positions in the human protein interaction network and were more similar to essential genes rather than cancer genes. However, PGS had significantly smaller intra-set distance (IAD) and inter-set distance (IED) in comparison with random sets (p-value < 0.001). Moreover, we also found that PGS tended to be distributed within network modules rather than between modules (p-value < 0.01), and the functional intersection of the modules enriched with PGS was closely related to cancer development and progression. Our research reveals the common network properties of cancer prognostic gene signatures in the human protein interactome. We argue that these are biologically meaningful and useful for understanding their molecular mechanism.

Biclustering-based association rule mining approach for predicting cancer-associated protein interactions

Protein-protein interactions (PPIs) have been widely used to understand different biological processes and cellular functions associated with several diseases like cancer. Although some cancer-related protein interaction databases are available, lack of experimental data and conflicting PPI data among different available databases have slowed down the cancer research. Therefore, in this study, the authors have focused on various proteins that are directly related to different types of cancer disease. They have prepared a PPI database between cancer-associated proteins with the rest of the human proteins. They have also incorporated the annotation type and direction of each interaction. Subsequently, a biclustering-based association rule mining algorithm is applied to predict new interactions with type and direction. This study shows the prediction power of association rule mining algorithm over the traditional classifier model without choosing a negative data set. The time complexity of the biclustering-based association rule mining is also analysed and compared to traditional association rule mining. The authors are able to discover 38 new PPIs which are not present in the cancer database. The biological relevance of these newly predicted interactions is analysed by published literature. Recognition of such interactions may accelerate a way of developing new drugs to prevent different cancer-related diseases.

Mechanism of tanshinones and phenolic acids from Danshen in the treatment of coronary heart disease based on co-expression network

Background

The tanshinones and phenolic acids in Salvia miltiorrhiza (also named Danshen) have been confirmed for the treatment of coronary heart disease (CHD), but the action mechanisms remain elusive.

Methods

In the current study, the co-expression protein interaction network (Ce-PIN) was used to illustrate the differences between the tanshinones and phenolic acids of Danshen in the treatment of CHD. By integrating the gene expression profile data and protein-protein interactions (PPIs) data, the Ce-PINs of tanshinones and phenolic acids were constructed. Then, the Ce-PINs were analyzed by gene ontology enrichment analyzed based on the optimal algorithm.

Results

It turned out that Danshen is able to treat CHD by regulating the blood circulation, immune response and lipid metabolism. However, phenolic acids may regulate the blood circulation by Extracellular calcium-sensing receptor (CaSR), Endothelin-1 receptor (EDNRA), Endothelin-1 receptor (EDNRB), Kininogen-1 (KNG1), tanshinones may regulate the blood circulation by Guanylate cyclase soluble subunit alpha-1 (GUCY1A3) and Guanylate cyclase soluble subunit beta-1 (GUCY1B3). In addition, both the phenolic acids and tanshinones may regulate the immune response or inflammation by T-cell surface glycoprotein CD4 (CD4), Receptor-type tyrosine-protein phosphatase C (PTPRC).

Conclusion

Through the same targets of the same biological process and different targets of the same biological process, the tanshinones and phenolic acids synergistically treat coronary heart disease.

Model-based clustering of time-evolving networks through temporal exponential-family random graph models

Dynamic networks are a general language for describing time-evolving complex systems, and discrete time network models provide an emerging statistical technique for various applications. It is a fundamental research question to detect a set of nodes sharing similar connectivity patterns in time-evolving networks. Our work is primarily motivated by detecting groups based on interesting features of the time-evolving networks (e.g., stability). In this work, we propose a model-based clustering framework for time-evolving networks based on discrete time exponential-family random graph models, which simultaneously allows both modeling and detecting group structure. To choose the number of groups, we use the conditional likelihood to construct an effective model selection criterion. Furthermore, we propose an efficient variational expectation-maximization (EM) algorithm to find approximate maximum likelihood estimates of network parameters and mixing proportions. The power of our method is demonstrated in simulation studies and empirical applications to international trade networks and the collaboration networks of a large research university.

DCARS: differential correlation across ranked samples

MOTIVATION

Genes act as a system and not in isolation. Thus, it is important to consider coordinated changes of gene expression rather than single genes when investigating biological phenomena such as the aetiology of cancer. We have developed an approach for quantifying how changes in the association between pairs of genes may inform the outcome of interest called Differential Correlation across Ranked Samples (DCARS). Modelling gene correlation across a continuous sample ranking does not require the dichotomisation of samples into two distinct classes and can identify differences in gene correlation across early, mid or late stages of the outcome of interest.

RESULTS

When we evaluated DCARS against the typical Fisher Z-transformation test for differential correlation, as well as a typical approach testing for interaction within a linear model, on real TCGA data, DCARS significantly ranked gene pairs containing known cancer genes more highly across several cancers. Similar results are found with our simulation study. DCARS was applied to 13 cancers datasets in TCGA, revealing several distinct relationships for which survival ranking was found to be associated with a change in correlation between genes. Furthermore, we demonstrated that DCARS can be used in conjunction with network analysis techniques to extract biological meaning from multi-layered and complex data.

AVAILABILITY AND IMPLEMENTATION

DCARS R package and sample data are available at https://github.com/shazanfar/DCARS. Publicly available data from The Cancer Genome Atlas (TCGA) was used using the TCGABiolinks R package. Supplementary Files and DCARS R package is available at https://github.com/shazanfar/DCARS.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

DNLC: differential network local consistency analysis

BACKGROUND

The biological network is highly dynamic. Functional relations between genes can be activated or deactivated depending on the biological conditions. On the genome-scale network, subnetworks that gain or lose local expression consistency may shed light on the regulatory mechanisms related to the changing biological conditions, such as disease status or tissue developmental stages.

RESULTS

In this study, we develop a new method to select genes and modules on the existing biological network, in which local expression consistency changes significantly between clinical conditions. The method is called DNLC: Differential Network Local Consistency. In simulations, our algorithm detected artificially created local consistency changes effectively. We applied the method on two publicly available datasets, and the method detected novel genes and network modules that were biologically plausible.

CONCLUSIONS

The new method is effective in finding modules in which the gene expression consistency change between clinical conditions. It is a useful tool that complements traditional differential expression analyses to make discoveries from gene expression data. The R package is available at https://cran.r-project.org/web/packages/DNLC.

Comparative analysis of differential proteome-wide protein-protein interaction network of Methanobrevibacter ruminantium M1

A proteome-wide protein-protein interaction (PPI) network of Methanobrevibacter ruminantium M1 (MRU), a predominant rumen methanogen, was constructed from its metabolic genes using a gene neighborhood algorithm and then compared with closely related rumen methanogens Using proteome-wide PPI approach, we constructed network encompassed 2194 edges and 637 nodes interacting with 634 genes. Network quality and robustness of functional modules were assessed with gene ontology terms. A structure-function-metabolism mapping for each protein has been carried out with efforts to extract experimental PPI concomitant information from the literature. The results of our study revealed that some topological properties of its network were robust for sharing homologous protein interactions across heterotrophic and hydrogenotrophic methanogens. MRU proteome has shown to establish many PPI sub-networks for associated metabolic subsystems required to survive in the rumen environment. MRU genome found to share interacting proteins from its PPI network involved in specific metabolic subsystems distinct to heterotrophic and hydrogenotrophic methanogens. Across these proteomes, the interacting proteins from differential PPI networks were shared in common for the biosynthesis of amino acids, nucleosides, and nucleotides and energy metabolism in which more fractions of protein pairs shared with Methanosarcina acetivorans. Our comparative study expedites our knowledge to understand a complex proteome network associated with typical metabolic subsystems of MRU and to improve its genome-scale reconstruction in the future.

Estrogen receptors promote NSCLC progression by modulating the membrane receptor signaling network: a systems biology perspective

Background

Estrogen receptors (ERs) are thought to play an important role in non-small cell lung cancer (NSCLC). However, the effect of ERs in NSCLC is still controversial and needs further investigation. A new consideration is that ERs may affect NSCLC progression through complicated molecular signaling networks rather than individual targets. Therefore, this study aims to explore the effect of ERs in NSCLC from the perspective of cancer systems biology.

Methods

The gene expression profile of NSCLC samples in TCGA dataset was analyzed by bioinformatics method. Variations of cell behaviors and protein expression were detected in vitro. The kinetic process of molecular signaling network was illustrated by a systemic computational model. At last, immunohistochemical (IHC) and survival analysis was applied to evaluate the clinical relevance and prognostic effect of key receptors in NSCLC.

Results

Bioinformatics analysis revealed that ERs might affect many cancer-related molecular events and pathways in NSCLC, particularly membrane receptor activation and signal transduction, which might ultimately lead to changes in cell behaviors. Experimental results confirmed that ERs could regulate cell behaviors including cell proliferation, apoptosis, invasion and migration; ERs also regulated the expression or activation of key members in membrane receptor signaling pathways such as epidermal growth factor receptor (EGFR), Notch1 and Glycogen synthase kinase-3β/β-Catenin (GSK3β/β-Catenin) pathways. Modeling results illustrated that the promotive effect of ERs in NSCLC was implemented by modulating the signaling network composed of EGFR, Notch1 and GSK3β/β-Catenin pathways; ERs maintained and enhanced the output of oncogenic signals by adding redundant and positive-feedback paths into the network. IHC results echoed that high expression of ERs, EGFR and Notch1 had a synergistic effect on poor prognosis of advanced NSCLC.

Conclusions

This study indicated that ERs were likely to promote NSCLC progression by modulating the integrated membrane receptor signaling network composed of EGFR, Notch1 and GSK3β/β-Catenin pathways and then affecting tumor cell behaviors. It also complemented the molecular mechanisms underlying the progression of NSCLC and provided new opportunities for optimizing therapeutic scheme of NSCLC.

GAIL: An interactive webserver for inference and dynamic visualization of gene-gene associations based on gene ontology guided mining of biomedical literature

In systems biology, inference of functional associations among genes is compelling because the construction of functional association networks facilitates biomarker discovery. Specifically, such gene associations in human can help identify putative biomarkers that can be used as diagnostic tools in treating patients. Although biomedical literature is considered a valuable data source for this task, currently only a limited number of webservers are available for mining gene-gene associations from the vast amount of biomedical literature using text mining techniques. Moreover, these webservers often have limited coverage of biomedical literature and also lack efficient and user-friendly tools to interpret and visualize mined relationships among genes. To address these limitations, we developed GAIL (Gene-gene Association Inference based on biomedical Literature), an interactive webserver that infers human gene-gene associations from Gene Ontology (GO) guided biomedical literature mining and provides dynamic visualization of the resulting association networks and various gene set enrichment analysis tools. We evaluate the utility and performance of GAIL with applications to gene signatures associated with systemic lupus erythematosus and breast cancer. Results show that GAIL allows effective interrogation and visualization of gene-gene networks and their subnetworks, which facilitates biological understanding of gene-gene associations. GAIL is available at http://chunglab.io/GAIL/.

Identification of LIFR, PIK3R1, and MMP12 as Novel Prognostic Signatures in Gallbladder Cancer Using Network-Based Module Analysis

Background: Gallbladder cancer (GBC) is a rare and aggressive malignancy of the biliary tract with a dismal survival rate. Effective biomarkers and therapeutic targets are urgently needed.

Methods: We analyzed gene expression profiles of GBC to identify differentially expressed genes (DEGs) and then used these DEGs to identify functional module biomarkers based on protein functional interaction (FI) networks. We further evaluated the module-gene protein expression and clinical significance with immunohistochemistry staining (IHC) in a tissue microarray (TMA) from 80 GBC samples.

Results: Five functional modules were identified. Module 0 included classical cancer signaling pathways, such as Ras and PI3K-Akt; and modules 1–4 included genes associated with muscle cells, fibrinogen, extracellular matrix, and integrins, respectively. We validated the expression of LIFR, PIK3R1, and MMP12, which were hubs or functional nodes in modules. Compared with paired peritumoural tissues, we found that the expression of LIFR (P = 0.002) and PIK3R1 (P = 0.046) proteins were significantly downregulated, and MMP12 (P = 0.006) was significantly upregulated. Further prognostic analysis showed that patients with low expression of LIFR had shorter overall survival than those with high expression (log-rank test P = 0.028), the same trend as for PIK3R1 (P = 0.053) and MMP12 (P = 0.006). Multivariate analysis indicated that expression of MMP12 protein (hazard ratio [HR] = 0.429; 95% confidence interval [CI] 0.198, 0.930; P = 0.032) was one of the significant independent prognostic factors for overall survival.

Conclusions: We found a highly reliable FI network, which revealed LIFR, PIK3R1, and MMP12 as novel prognostic biomarker candidates for GBC. These findings could accelerate biomarker discovery and therapeutic development in this cancer.

Network Medicine in the Age of Biomedical Big Data

Network medicine is an emerging area of research dealing with molecular and genetic interactions, network biomarkers of disease, and therapeutic target discovery. Large-scale biomedical data generation offers a unique opportunity to assess the effect and impact of cellular heterogeneity and environmental perturbations on the observed phenotype. Marrying the two, network medicine with biomedical data provides a framework to build meaningful models and extract impactful results at a network level. In this review, we survey existing network types and biomedical data sources. More importantly, we delve into ways in which the network medicine approach, aided by phenotype-specific biomedical data, can be gainfully applied. We provide three paradigms, mainly dealing with three major biological network archetypes: protein-protein interaction, expression-based, and gene regulatory networks. For each of these paradigms, we discuss a broad overview of philosophies under which various network methods work. We also provide a few examples in each paradigm as a test case of its successful application. Finally, we delineate several opportunities and challenges in the field of network medicine. We hope this review provides a lexicon for researchers from biological sciences and network theory to come on the same page to work on research areas that require interdisciplinary expertise. Taken together, the understanding gained from combining biomedical data with networks can be useful for characterizing disease etiologies and identifying therapeutic targets, which, in turn, will lead to better preventive medicine with translational impact on personalized healthcare.

PIMKL: Pathway-Induced Multiple Kernel Learning

Reliable identification of molecular biomarkers is essential for accurate patient stratification. While state-of-the-art machine learning approaches for sample classification continue to push boundaries in terms of performance, most of these methods are not able to integrate different data types and lack generalization power, limiting their application in a clinical setting. Furthermore, many methods behave as black boxes, and we have very little understanding about the mechanisms that lead to the prediction. While opaqueness concerning machine behavior might not be a problem in deterministic domains, in health care, providing explanations about the molecular factors and phenotypes that are driving the classification is crucial to build trust in the performance of the predictive system. We propose Pathway-Induced Multiple Kernel Learning (PIMKL), a methodology to reliably classify samples that can also help gain insights into the molecular mechanisms that underlie the classification. PIMKL exploits prior knowledge in the form of a molecular interaction network and annotated gene sets, by optimizing a mixture of pathway-induced kernels using a Multiple Kernel Learning (MKL) algorithm, an approach that has demonstrated excellent performance in different machine learning applications. After optimizing the combination of kernels to predict a specific phenotype, the model provides a stable molecular signature that can be interpreted in the light of the ingested prior knowledge and that can be used in transfer learning tasks.

The reliable classification of biomedical samples to predict phenotypic differences requires not only robust methods, but also interpretable approaches that can explain the reasons behind a prediction. A team led by María Rodríguez Martínez at IBM Research - Zürich has developed PIMKL, a methodology that exploits prior knowledge and enables the integration of multiple types of data with varying predictive power. Even when noisy datasets are simultaneously analyzed, PIMKL is able to discard uninformative data and achieve strong prediction power. Importantly, PIMKL produces a molecular signature that enables the interpretation of the results in terms of known biological functions. This signature can be transferred to other cohorts without loss of performance, demonstrating surprising robustness across cohorts. Interpretable algorithms can effectively help gain insights about disease mechanisms and build trust in a model.

Identifying cancer prognostic modules by module network analysis

Background

The identification of prognostic genes that can distinguish the prognostic risks of cancer patients remains a significant challenge. Previous works have proven that functional gene sets were more reliable for this task than the gene signature. However, few works have considered the cross-talk among functional gene sets, which may result in neglecting important prognostic gene sets for cancer.

Results

Here, we proposed a new method that considers both the interactions among modules and the prognostic correlation of the modules to identify prognostic modules in cancers. First, dense sub-networks in the gene co-expression network of cancer patients were detected. Second, cross-talk between every two modules was identified by a permutation test, thus generating the module network. Third, the prognostic correlation of each module was evaluated by the resampling method. Then, the GeneRank algorithm, which takes the module network and the prognostic correlations of all the modules as input, was applied to prioritize the prognostic modules. Finally, the selected modules were validated by survival analysis in various data sets. Our method was applied in three kinds of cancers, and the results show that our method succeeded in identifying prognostic modules in all the three cancers. In addition, our method outperformed state-of-the-art methods. Furthermore, the selected modules were significantly enriched with known cancer-related genes and drug targets of cancer, which may indicate that the genes involved in the modules may be drug targets for therapy.

Conclusions

We proposed a useful method to identify key modules in cancer prognosis and our prognostic genes may be good candidates for drug targets.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-2674-z) contains supplementary material, which is available to authorized users.

Identification-Method Research for Open-Source Software Ecosystems

In recent years, open-source software (OSS) development has grown, with many developers around the world working on different OSS projects. A variety of open-source software ecosystems have emerged, for instance, GitHub, StackOverflow, and SourceForge. One of the most typical social-programming and code-hosting sites, GitHub, has amassed numerous open-source-software projects and developers in the same virtual collaboration platform. Since GitHub itself is a large open-source community, it hosts a collection of software projects that are developed together and coevolve. The great challenge here is how to identify the relationship between these projects, i.e., project relevance. Software-ecosystem identification is the basis of other studies in the ecosystem. Therefore, how to extract useful information in GitHub and identify software ecosystems is particularly important, and it is also a research area in symmetry. In this paper, a Topic-based Project Knowledge Metrics Framework (TPKMF) is proposed. By collecting the multisource dataset of an open-source ecosystem, project-relevance analysis of the open-source software is carried out on the basis of software-ecosystem identification. Then, we used our Spectral Clustering algorithm based on Core Project (CP-SC) to identify software-ecosystem projects and further identify software ecosystems. We verified that most software ecosystems usually contain a core software project, and most other projects are associated with it. Furthermore, we analyzed the characteristics of the ecosystem, and we also found that interactive information has greater impact on project relevance. Finally, we summarize the Topic-based Project Knowledge Metrics Framework.

A data-driven interactome of synergistic genes improves network-based cancer outcome prediction

Robustly predicting outcome for cancer patients from gene expression is an important challenge on the road to better personalized treatment. Network-based outcome predictors (NOPs), which considers the cellular wiring diagram in the classification, hold much promise to improve performance, stability and interpretability of identified marker genes. Problematically, reports on the efficacy of NOPs are conflicting and for instance suggest that utilizing random networks performs on par to networks that describe biologically relevant interactions. In this paper we turn the prediction problem around: instead of using a given biological network in the NOP, we aim to identify the network of genes that truly improves outcome prediction. To this end, we propose SyNet, a gene network constructed ab initio from synergistic gene pairs derived from survival-labelled gene expression data. To obtain SyNet, we evaluate synergy for all 69 million pairwise combinations of genes resulting in a network that is specific to the dataset and phenotype under study and can be used to in a NOP model. We evaluated SyNet and 11 other networks on a compendium dataset of >4000 survival-labelled breast cancer samples. For this purpose, we used cross-study validation which more closely emulates real world application of these outcome predictors. We find that SyNet is the only network that truly improves performance, stability and interpretability in several existing NOPs. We show that SyNet overlaps significantly with existing gene networks, and can be confidently predicted (~85% AUC) from graph-topological descriptions of these networks, in particular the breast tissue-specific network. Due to its data-driven nature, SyNet is not biased to well-studied genes and thus facilitates post-hoc interpretation. We find that SyNet is highly enriched for known breast cancer genes and genes related to e.g. histological grade and tamoxifen resistance, suggestive of a role in determining breast cancer outcome.

Graph-Based Hub Gene Selection Technique Using Protein Interaction Information: Application to Sample Classification

Classification of samples of gene expression profile plays a significant role in prediction and diagnosis of diseases. In the task of sample classification, a robust feature selection algorithm is very much essential to identify the important genes from the high dimensional gene expression data. This paper explores the information of protein-protein interaction with a graph mining technique for finding a proper subset of features (genes), which further takes part in sample classification. Here, our contribution for feature selection is three-fold: first, all the genes are grouped into different clusters based on the integrated information of the gene expression values and their protein interactions using a multi-objective optimization based clustering approach. Second, the confidence scores of the protein interactions are incorporated in a popular graph mining algorithm namely Goldberg algorithm to find out the relevant features. These features are the topologically and functionally significant genes, named as hub genes. Finally, these hub genes are identified varying the degrees of the nodes, and those are utilized for the sample classification task. Different machine learning classifiers are exploited for this purpose, and the classification performance is measured with respect to various performance metrics namely accuracy, sensitivity, specificity, precision, F-measure, and Mathews coefficient correlation. Comparative analysis with respect to two baselines and several existing approaches proves the efficiency of the proposed approach. Furthermore, the robustness of the identified hub-gene modules is endorsed using some strong biological significance analysis.

PPICurator: A Tool for Extracting Comprehensive Protein-Protein Interaction Information

Protein-protein interaction extraction through biological literature curation is widely employed for proteome analysis. There is a strong need for a tool that can assist researchers in extracting comprehensive PPI information through literature curation, which is critical in research on protein, for example, construction of protein interaction network, identification of protein signaling pathway, and discovery of meaningful protein interaction. However, most of current tools can only extract PPI relations. None of them are capable of extracting other important PPI information, such as interaction directions, effects, and functional annotations. To address these issues, this paper proposes PPICurator, a novel tool for extracting comprehensive PPI information with a variety of logic and syntax features based on a new support vector machine classifier. PPICurator provides a friendly web-based user interface. It is a platform that automates the extraction of comprehensive PPI information through literature, including PPI relations, as well as their confidential scores, interaction directions, effects, and functional annotations. Thus, PPICurator is more comprehensive than state-of-the-art tools. Moreover, it outperforms state-of-the-art tools in the accuracy of PPI relation extraction measured by F-score and recall on the widely used open datasets. PPICurator is available at https://ppicurator.hupo.org.cn.

Integrating molecular networks with genetic variant interpretation for precision medicine

More reliable and cheaper sequencing technologies have revealed the vast mutational landscapes characteristic of many phenotypes. The analysis of such genetic variants has led to successful identification of altered proteins underlying many Mendelian disorders. Nevertheless the simple one‐variant one‐phenotype model valid for many monogenic diseases does not capture the complexity of polygenic traits and disorders. Although experimental and computational approaches have improved detection of functionally deleterious variants and important interactions between gene products, the development of comprehensive models relating genotype and phenotypes remains a challenge in the field of genomic medicine. In this context, a new view of the pathologic state as significant perturbation of the network of interactions between biomolecules is crucial for the identification of biochemical pathways associated with complex phenotypes. Seminal studies in systems biology combined the analysis of genetic variation with protein–protein interaction networks to demonstrate that even as biological systems evolve to be robust to genetic variation, their topologies create disease vulnerabilities. More recent analyses model the impact of genetic variants as changes to the “wiring” of the interactome to better capture heterogeneity in genotype–phenotype relationships. These studies lay the foundation for using networks to predict variant effects at scale using machine‐learning or algorithmic approaches. A wealth of databases and resources for the annotation of genotype–phenotype relationships have been developed to support developments in this area. This overview describes how study of the molecular interactome has generated insights linking the organization of biological systems to disease mechanism, and how this information can enable precision medicine.

This article is categorized under:

Translational, Genomic, and Systems Medicine > Translational Medicine

Biological Mechanisms > Cell Signaling

Models of Systems Properties and Processes > Mechanistic Models

Analytical and Computational Methods > Computational Methods

Interaction paths promote module integration and network-level robustness of spliceosome to cascading effects

The functionality of distinct types of protein networks depends on the patterns of protein-protein interactions. A problem to solve is understanding the fragility of protein networks to predict system malfunctioning due to mutations and other errors. Spectral graph theory provides tools to understand the structural and dynamical properties of a system based on the mathematical properties of matrices associated with the networks. We combined two of such tools to explore the fragility to cascading effects of the network describing protein interactions within a key macromolecular complex, the spliceosome. Using S. cerevisiae as a model system we show that the spliceosome network has more indirect paths connecting proteins than random networks. Such multiplicity of paths may promote routes to cascading effects to propagate across the network. However, the modular network structure concentrates paths within modules, thus constraining the propagation of such cascading effects, as indicated by analytical results from the spectral graph theory and by numerical simulations of a minimal mathematical model parameterized with the spliceosome network. We hypothesize that the concentration of paths within modules favors robustness of the spliceosome against failure, but may lead to a higher vulnerability of functional subunits, which may affect the temporal assembly of the spliceosome. Our results illustrate the utility of spectral graph theory for identifying fragile spots in biological systems and predicting their implications.

Butyrate-Induced In Vitro Colonocyte Differentiation Network Model Identifies ITGB1, SYK, CDKN2A, CHAF1A, and LRP1 as the Prognostic Markers for Colorectal Cancer Recurrence

Numerous mechanisms are believed to contribute to the role of dietary fiber-derived butyrate in the protection against the development of colorectal cancers (CRCs). To identify the most crucial butyrate-regulated genes, we exploited whole genome microarray of HT-29 cells differentiated in vitro by butyrate treatment. Butyrate differentiates HT-29 cells by relaxing the perturbation, caused by mutations of Adenomatous polyposis coli (APC) and TP53 genes, the most frequent mutations observed in CRC. We constructed protein-protein interaction network (PPIN) with the differentially expressed genes after butyrate treatment and extracted the hub genes from the PPIN, which also participated in the APC-TP53 network. The idea behind this approach was that the expression of these hub genes also regulated cell differentiation, and subsequently CRC prognosis by evading the APC-TP53 mutational effect. We used mRNA expression profile of these critical hub genes from seven large CRC cohorts. Logistic Regression showed strong evidence for association of these common hubs with CRC recurrence. In this study, we exploited PPIN to reduce the dimension of microarray biologically and identified five prognostic markers for the CRC recurrence, which were validated across different datasets. Moreover, these five biomarkers we identified increase the predictive value of the TNM staging for CRC recurrence.

Pathway-based subnetworks enable cross-disease biomarker discovery

Biomarkers lie at the heart of precision medicine. Surprisingly, while rapid genomic profiling is becoming ubiquitous, the development of biomarkers usually involves the application of bespoke techniques that cannot be directly applied to other datasets. There is an urgent need for a systematic methodology to create biologically-interpretable molecular models that robustly predict key phenotypes. Here we present SIMMS (Subnetwork Integration for Multi-Modal Signatures): an algorithm that fragments pathways into functional modules and uses these to predict phenotypes. We apply SIMMS to multiple data types across five diseases, and in each it reproducibly identifies known and novel subtypes, and makes superior predictions to the best bespoke approaches. To demonstrate its ability on a new dataset, we profile 33 genes/nodes of the PI3K pathway in 1734 FFPE breast tumors and create a four-subnetwork prediction model. This model out-performs a clinically-validated molecular test in an independent cohort of 1742 patients. SIMMS is generic and enables systematic data integration for robust biomarker discovery.

Analyzing Differential Regulatory Networks Modulated by Continuous-State Genomic Features in Glioblastoma Multiforme

Gene regulatory networks are a global representation of complex interactions between molecules that dictate cellular behavior. Study of a regulatory network modulated by single or multiple modulators' expression levels, including microRNAs (miRNAs) and transcription factors (TFs), in different conditions can further reveal the modulators' roles in diseases such as cancers. Existing computational methods for identifying such modulated regulatory networks are typically carried out by comparing groups of samples dichotomized with respect to the modulator status, ignoring the fact that most biological features are intrinsically continuous variables. Here, we devised a sliding window-based regression scheme and proposed the Regression-based Inference of Modulation (RIM) algorithm to infer the dynamic gene regulation modulated by continuous-state modulators. We demonstrated the improvement in performance as well as computation efficiency achieved by RIM. Applying RIM to genome-wide expression profiles of 520 glioblastoma multiforme (GBM) tumors, we investigated miRNA- and TF-modulated gene regulatory networks and showed their association with dynamic cellular processes and brain-related functions in GBM. Overall, the proposed algorithm provides an efficient and robust scheme for comprehensively studying modulated gene regulatory networks.

wTO: an R package for computing weighted topological overlap and a consensus network with integrated visualization tool

BACKGROUND

Network analyses, such as of gene co-expression networks, metabolic networks and ecological networks have become a central approach for the systems-level study of biological data. Several software packages exist for generating and analyzing such networks, either from correlation scores or the absolute value of a transformed score called weighted topological overlap (wTO). However, since gene regulatory processes can up- or down-regulate genes, it is of great interest to explicitly consider both positive and negative correlations when constructing a gene co-expression network.

RESULTS

Here, we present an R package for calculating the weighted topological overlap (wTO), that, in contrast to existing packages, explicitly addresses the sign of the wTO values, and is thus especially valuable for the analysis of gene regulatory networks. The package includes the calculation of p-values (raw and adjusted) for each pairwise gene score. Our package also allows the calculation of networks from time series (without replicates). Since networks from independent datasets (biological repeats or related studies) are not the same due to technical and biological noise in the data, we additionally, incorporated a novel method for calculating a consensus network (CN) from two or more networks into our R package. To graphically inspect the resulting networks, the R package contains a visualization tool, which allows for the direct network manipulation and access of node and link information. When testing the package on a standard laptop computer, we can conduct all calculations for systems of more than 20,000 genes in under two hours. We compare our new wTO package to state of art packages and demonstrate the application of the wTO and CN functions using 3 independently derived datasets from healthy human pre-frontal cortex samples. To showcase an example for the time series application we utilized a metagenomics data set.

CONCLUSION

In this work, we developed a software package that allows the computation of wTO networks, CNs and a visualization tool in the R statistical environment. It is publicly available on CRAN repositories under the GPL -2 Open Source License ( https://cran.r-project.org/web/packages/wTO/ ).

The modular network structure of the mutational landscape of Acute Myeloid Leukemia

Acute myeloid leukemia (AML) is associated with the sequential accumulation of acquired genetic alterations. Although at diagnosis cytogenetic alterations are frequent in AML, roughly 50% of patients present an apparently normal karyotype (NK), leading to a highly heterogeneous prognosis. Due to this significant heterogeneity, it has been suggested that different molecular mechanisms may trigger the disease with diverse prognostic implications. We performed whole-exome sequencing (WES) of tumor-normal matched samples of de novo AML-NK patients lacking mutations in NPM1, CEBPA or FLT3-ITD to identify new gene mutations with potential prognostic and therapeutic relevance to patients with AML. Novel candidate-genes, together with others previously described, were targeted resequenced in an independent cohort of 100 de novo AML patients classified in the cytogenetic intermediate-risk (IR) category. A mean of 4.89 mutations per sample were detected in 73 genes, 35 of which were mutated in more than one patient. After a network enrichment analysis, we defined a single in silico model and established a set of seed-genes that may trigger leukemogenesis in patients with normal karyotype. The high heterogeneity of gene mutations observed in AML patients suggested that a specific alteration could not be as essential as the interaction of deregulated pathways.

Typing tumors using pathways selected by somatic evolution

Many recent efforts to analyze cancer genomes involve aggregation of mutations within reference maps of molecular pathways and protein networks. Here, we find these pathway studies are impeded by molecular interactions that are functionally irrelevant to cancer or the patient’s tumor type, as these interactions diminish the contrast of driver pathways relative to individual frequently mutated genes. This problem can be addressed by creating stringent tumor-specific networks of biophysical protein interactions, identified by signatures of epistatic selection during tumor evolution. Using such an evolutionarily selected pathway (ESP) map, we analyze the major cancer genome atlases to derive a hierarchical classification of tumor subtypes linked to characteristic mutated pathways. These pathways are clinically prognostic and predictive, including the TP53-AXIN-ARHGEF17 combination in liver and CYLC2-STK11-STK11IP in lung cancer, which we validate in independent cohorts. This ESP framework substantially improves the definition of cancer pathways and subtypes from tumor genome data.

Informative pathways driving cancer pathogenesis and subtypes can be difficult to identify in the presence of many gene interactions irrelevant to cancer. Here, the authors describe an approach for cancer gene pathway analysis based on key molecular interactions that drive cancer in relevant tissue types, and they assemble a focused map of Evolutionarily Selected Pathways (ESP) with interactions supported by both protein–protein binding and genetic epistasis during somatic tumor evolution.

Analyzing of Molecular Networks for Human Diseases and Drug Discovery

Molecular networks represent the interactions and relations of genes/proteins, and also encode molecular mechanisms of biological processes, development and diseases. Among the molecular networks, protein-protein Interaction Networks (PINs) have become effective platforms for uncovering the molecular mechanisms of diseases and drug discovery. PINs have been constructed for various organisms and utilized to solve many biological problems. In human, most proteins present their complex functions by interactions with other proteins, and the sum of these interactions represents the human protein interactome. Especially in the research on human disease and drugs, as an emerging tool, the PIN provides a platform to systematically explore the molecular complexities of specific diseases and the references for drug design. In this review, we summarized the commonly used approaches to aid disease research and drug discovery with PINs, including the network topological analysis, identification of novel pathways, drug targets and sub-network biomarkers for diseases. With the development of bioinformatic techniques and biological networks, PINs will play an increasingly important role in human disease research and drug discovery.

Pathway-based gene-gene interaction network modelling to predict potential biomarkers of essential hypertension

Essential hypertension (EH) is a major risk factor for cardiovascular disease. Despite considerable efforts to elucidate the pathogenesis of EH, there is an imperious need for novel indicators of EH. This study aimed to develop a method to predict potential biomarkers of EH from the point of view of network. A pathway-based gene-gene interaction (GGI) network model was constructed and analyzed, containing 116 nodes and 1272 connections. The nodes represented EH-related genes, and that connections represented their interactions. The network showed a small-world property and uneven degree distribution, suggesting that a few highly interconnected hubs played a vital role in EH. An inherent hierarchy and assortative mixing pattern were also observed in the network. GNAS, GNB3, PF4 and PPBP showed the highest values of degrees and centrality indices, and were chosen as potential biomarkers of EH. A two-mode network model based on the potential biomarkers demonstrated that hemostasis and GPCR ligand binding pathway were key pathways contributing to EH. Results of this study improve our current understanding of the molecular mechanisms driving EH. The selected genes and pathways have the potential to be used in the diagnosis and treatment of EH. Moreover, the combination of pathway analysis and complex network methodology provides a novel strategy for searching new genetic indicators of complex diseases.

Protein-Protein Interaction Network Analysis of Salivary Proteomic Data in Oral Cancer Cases

Background:

Oral cancer is a frequently encountered neoplasm of the head and neck region, being the eight most common type of human malignancy worldwide. Despite improvement in its control, morbidity and mortality rates have improved little in the past decades. Therefore, prevention and/or early detection are a high priority. Proteomics with network analysis have emerged as a powerful tool to identify important proteins associated with cancer development and progression that can be potential targets for early diagnosis. In the present study, network- based protein- protein interactions (PPI) for oral cancer were identified and then analyzed for use as key proteins/potential biomarkers.

Material and Methods:

Gene expression data in articles which focused on saliva proteomics of oral cancer were collected and 74 candidate genes or proteins were extracted. Related protein networks of differentially expressed proteins were explored and visualized using cytoscape software. Further PPI analysis was performed by Molecular Complex Detection (MCODE) and BiNGO methods.

Results:

Network analysis of genes/proteins related to oral cancer identified kininogen-1, angiotensinogen, annexin A1, IL-8, IgG heavy variable and constant chains, CRP, collagen alpha-1 and fibronectin as 9 hub-bottleneck proteins. In addition, based on clustering with the MCODE tool, vitronectin, collagen alpha-2, IL-8 and integrin alpha-v were established as 5 distinct seed proteins.

Conclusion:

A hub-bottleneck protein panel may offer a potential /candidate biomarker pattern for diagnosis and treatment of oral cancer disease. Further investigation and validation of these proteins are warranted.

Network biology discovers pathogen contact points in host protein-protein interactomes

In all organisms, major biological processes are controlled by complex protein-protein interactions networks (interactomes), yet their structural complexity presents major analytical challenges. Here, we integrate a compendium of over 4300 phenotypes with Arabidopsis interactome (AI-1MAIN). We show that nodes with high connectivity and betweenness are enriched and depleted in conditional and essential phenotypes, respectively. Such nodes are located in the innermost layers of AI-1MAIN and are preferential targets of pathogen effectors. We extend these network-centric analyses to Cell Surface Interactome (CSILRR) and predict its 35 most influential nodes. To determine their biological relevance, we show that these proteins physically interact with pathogen effectors and modulate plant immunity. Overall, our findings contrast with centrality-lethality rule, discover fast information spreading nodes, and highlight the structural properties of pathogen targets in two different interactomes. Finally, this theoretical framework could possibly be applicable to other inter-species interactomes to reveal pathogen contact points.

iCAVE: an open source tool for visualizing biomolecular networks in 3D, stereoscopic 3D and immersive 3D

Visualizations of biomolecular networks assist in systems-level data exploration in many cellular processes. Data generated from high-throughput experiments increasingly inform these networks, yet current tools do not adequately scale with concomitant increase in their size and complexity. We present an open source software platform, interactome-CAVE (iCAVE), for visualizing large and complex biomolecular interaction networks in 3D. Users can explore networks (i) in 3D using a desktop, (ii) in stereoscopic 3D using 3D-vision glasses and a desktop, or (iii) in immersive 3D within a CAVE environment. iCAVE introduces 3D extensions of known 2D network layout, clustering, and edge-bundling algorithms, as well as new 3D network layout algorithms. Furthermore, users can simultaneously query several built-in databases within iCAVE for network generation or visualize their own networks (e.g., disease, drug, protein, metabolite). iCAVE has modular structure that allows rapid development by addition of algorithms, datasets, or features without affecting other parts of the code. Overall, iCAVE is the first freely available open source tool that enables 3D (optionally stereoscopic or immersive) visualizations of complex, dense, or multi-layered biomolecular networks. While primarily designed for researchers utilizing biomolecular networks, iCAVE can assist researchers in any field.