LncRNA and mRNA integration network reconstruction reveals novel key regulators in esophageal squamous-cell carcinoma

Many experimental and computational studies have identified key protein coding genes in initiation and progression of esophageal squamous cell carcinoma (ESCC). However, the number of researches that tried to reveal the role of long non-coding RNAs (lncRNAs) in ESCC has been limited. LncRNAs are one of the important regulators of cancers which are transcribed dominantly in the genome and in various conditions. The main goal of this study was to use a systems biology approach to predict novel lncRNAs as well as protein coding genes associated with ESCC and assess their prognostic values. By using microarray expression data for mRNAs and lncRNAs from a large number of ESCC patients, we utilized “Weighted Gene Co-expression Network Analysis” (WGCNA) method to make a big coding-non-coding gene co-expression network, and discovered important functional modules. Gene set enrichment and pathway analysis revealed major biological processes and pathways involved in these modules. After selecting some protein coding genes involved in biological processes and pathways related to cancer, we used “LncTar”, a computational tool to predict potential interactions between these genes and lncRNAs. By combining interaction results with Pearson correlations, we introduced some novel lncRNAs with putative key regulatory roles in the network. Survival analysis with Kaplan-Meier estimator and Log-rank test statistic confirmed that most of the introduced genes are associated with poor prognosis in ESCC. Overall, our study reveals novel protein coding genes and lncRNAs associated with ESCC, along with their predicted interactions. Based on the promising results of survival analysis, these genes can be used as good estimators of patients' survival, or even can be analyzed further as new potential signatures or targets for the therapy of ESCC disease.     

Deduction of intracellular sub-systems from a topological description of the network

Non-linear behaviour of biochemical networks, such as intracellular gene, protein or metabolic networks, is commonly represented using graphs of the underlying topology. Nodes represent abundance of molecules and edges interactions between pairs of molecules. These graphs are linear and thus based on an implicit linearization of the kinetic reactions in one or several dynamic modes of the total system. It is common to use data from different sources -- experiments conducted under different conditions or even on different species -- meaning that the graph will be a superposition of linearizations made in many different modes. The mixing of different modes makes it hard to identify functional modules, that is sub-systems that carry out a specific biological function, since the graph will contain many interactions that do not naturally occur at the same time. The ability to establish a boundary between the sub-system and its environment is critical in the definition of a module, contrary to a motif in which only internal interactions count. Identification of functional modules should therefore be done on graphs depicting the mode in which their function is carried out, i.e. graphs that only contain edges representing interactions active in the specific mode. In general, when an interaction between two molecules is established, one should always state the mode of the system in which it is active.     

Edge-based scoring and searching method for identifying condition-responsive protein-protein interaction sub-network

MOTIVATION: Current high-throughput protein-protein interaction (PPI) data do not provide information about the condition(s) under which the interactions occur. Thus, the identification of condition-responsive PPI sub-networks is of great importance for investigating how a living cell adapts to changing environments. RESULTS: In this article, we propose a novel edge-based scoring and searching approach to extract a PPI sub-network responsive to conditions related to some investigated gene expression profiles. Using this approach, what we constructed is a sub-network connected by the selected edges (interactions), instead of only a set of vertices (proteins) as in previous works. Furthermore, we suggest a systematic approach to evaluate the biological relevance of the identified responsive sub-network by its ability of capturing condition-relevant functional modules. We apply the proposed method to analyze a human prostate cancer dataset and a yeast cell cycle dataset. The results demonstrate that the edge-based method is able to efficiently capture relevant protein interaction behaviors under the investigated conditions. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.     

Multilevel regulation of protein-protein interactions in biological circuitry

Protein-protein interactions are central to biology and, in this 'post-genomic era', prediction of these interactions has become the goal of many computational efforts. Close inspection of even relatively simple biological regulatory circuitry reveals multiple levels of control of the contributing protein interactions. The fundamental probability that an interaction will occur under a given set of conditions is difficult to predict because the relationship between structure and energy is not known. Layered on this basic difficulty are allosteric control mechanisms involving post-translational modification or small ligand binding. In addition, many biological processes involve multiple protein-protein interactions, some of which may be cooperative or even competitive. Finally, although the emphasis in predicting protein interactions is based on equilibrium thermodynamic principles, kinetics can be a major controlling feature in these systems. This complexity reinforces the necessity of performing detailed quantitative studies of the component interactions of complex biological regulatory systems. Results of such studies will help us to bridge the gap between our knowledge of structure and our understanding of functional biology.     

High‐order combination effects and biological robustness

Biological systems are robust, in that they can maintain stable phenotypes under varying conditions or attacks. Biological systems are also complex, being organized into many functional modules that communicate through interlocking pathways and feedback mechanisms. In these systems, robustness and complexity are linked because both qualities arise from the same underlying mechanisms. When perturbed by multiple attacks, such complex systems become fragile in both theoretical and experimental studies, and this fragility depends on the number of agents applied. We explore how this relationship can be used to study the functional robustness of a biological system using systematic high‐order combination experiments. This presents a promising approach toward many biomedical and bioengineering challenges. For example, high‐order experiments could determine the point of fragility for pathogenic bacteria and might help identify optimal treatments against multi‐drug resistance. Such studies would also reinforce the growing appreciation that biological systems are best manipulated not by targeting a single protein, but by modulating the set of many nodes that can selectively control a system's functional state.     

Module-Based Association Analysis for Omics Data with Network Structure

Module-based analysis (MBA) aims to evaluate the effect of a group of biological elements sharing common features, such as SNPs in the same gene or metabolites in the same pathways, and has become an attractive alternative to traditional single bio-element approaches. Because bio-elements regulate and interact with each other as part of network, incorporating network structure information can more precisely model the biological effects, enhance the ability to detect true associations, and facilitate our understanding of the underlying biological mechanisms. How-ever, most MBA methods ignore the network structure information, which depicts the interaction and regulation relationship among basic functional units in biology system. We construct the con-nectivity kernel and the topology kernel to capture the relationship among bio-elements in a mod-ule, and use a kernel machine framework to evaluate the joint effect of bio-elements. Our proposed kernel machine approach directly incorporates network structure so to enhance the study effi-ciency; it can assess interactions among modules, account covariates, and is computational effi-cient. Through simulation studies and real data application, we demonstrate that the proposed network-based methods can have markedly better power than the approaches ignoring network information under a range of scenarios.     

Identification of Risk Pathways and Functional Modules for Coronary Artery Disease Based on Genome-wide SNP Data

Coronary artery disease(CAD) is a complex human disease, involving multiple genes and their nonlinear interactions, which often act in a modular fashion. Genome-wide single nucleotide polymorphism(SNP) profiling provides an effective technique to unravel these underlying genetic interplays or their functional involvements for CAD. This study aimed to identify the susceptible pathways and modules for CAD based on SNP omics. First, the Wellcome Trust Case Control Consortium(WTCCC) SNP datasets of CAD and control samples were used to assess the jointeffect of multiple genetic variants at the pathway level, using logistic kernel machine regression model. Then, an expanded genetic network was constructed by integrating statistical gene–gene interactions involved in these susceptible pathways with their protein–protein interaction(PPI)knowledge. Finally, risk functional modules were identified by decomposition of the network. Of 276 KEGG pathways analyzed, 6 pathways were found to have a significant effect on CAD. Other than glycerolipid metabolism, glycosaminoglycan biosynthesis, and cardiac muscle contraction pathways, three pathways related to other diseases were also revealed, including Alzheimer's disease, non-alcoholic fatty liver disease, and Huntington's disease. A genetic epistatic network of 95 genes was further constructed using the abovementioned integrative approach. Of 10 functional modules derived from the network, 6 have been annotated to phospholipase C activity and cell adhesion molecule binding, which also have known functional involvement in Alzheimer's disease.These findings indicate an overlap of the underlying molecular mechanisms between CAD and Alzheimer's disease, thus providing new insights into the molecular basis for CAD and its molecular relationships with other diseases.     

Identifying responsive modules by mathematical programming: an application to budding yeast cell cycle

High-throughput biological data offer an unprecedented opportunity to fully characterize biological processes. However, how to extract meaningful biological information from these datasets is a significant challenge. Recently, pathway-based analysis has gained much progress in identifying biomarkers for some phenotypes. Nevertheless, these so-called pathway-based methods are mainly individual-gene-based or molecule-complex-based analyses. In this paper, we developed a novel module-based method to reveal causal or dependent relations between network modules and biological phenotypes by integrating both gene expression data and protein-protein interaction network. Specifically, we first formulated the identification problem of the responsive modules underlying biological phenotypes as a mathematical programming model by exploiting phenotype difference, which can also be viewed as a multi-classification problem. Then, we applied it to study cell-cycle process of budding yeast from microarray data based on our biological experiments, and identified important phenotype- and transition-based responsive modules for different stages of cell-cycle process. The resulting responsive modules provide new insight into the regulation mechanisms of cell-cycle process from a network viewpoint. Moreover, the identification of transition modules provides a new way to study dynamical processes at a functional module level. In particular, we found that the dysfunction of a well-known module and two new modules may directly result in cell cycle arresting at S phase. In addition to our biological experiments, the identified responsive modules were also validated by two independent datasets on budding yeast cell cycle.     

Building functional modules from molecular interactions

The main reaction pathways in the living cell are carried out by functional modules--namely, macromolecular machines with compact structure or ensembles that change their composition and/or organization during function. Modules define themselves by spatial sequestration, chemical specificity and a characteristic time domain within which their function proceeds. On receiving a specific input, modules go through functional cycles, with phases of increasing and decreasing complexity of molecular interactions. Here, we discuss how such modules are formed and the experimental and theoretical approaches that can be used to investigate them, using examples from polynucleotide-protein interactions, vesicle transport and signal transduction to illustrate the underlying principles. Further progress in this field, where systems biology and biochemistry meet, will depend on iterative validation of the experimental and theoretical approaches.     

Systematic identification of common functional modules related to heart failure with different etiologies

The development of heart failure (HF) is a complex process that can be initiated by multiple etiologies. Identifying common functional modules associated with HF is a challenging task. Here, we developed a systems method to identify these common functional modules by integrating multiple expression profiles, protein interactions from four species, gene function annotations, and text information. We identified 1439 consistently differentially expressed genes (CDEGs) across HF with different etiologies by applying three meta-analysis methods to multiple HF-related expression profiles. Using a weighted human interaction network constructed by combining interaction data from multiple species, we extracted 60 candidate CDEG modules. We further evaluated the functional relevance of each module by using expression, interaction network, functional annotations, and text information together. Finally, five functional modules with significant biological relevance were identified. We found that almost half of the genes in these modules are hubs in the weighted network, and that these modules can accurately classify HF patients from healthy subjects. We also identified many significantly enriched biological processes that contribute to the pathophysiology of HF, including two new ones, RNA splicing and vesicle-mediated protein transport. In summary, we proposed a novel framework to analyze common functional modules related to HF with different etiologies. Our findings provide important insights into the complex mechanism of HF. Further biological experimentations should be required to validate these novel biological processes.     

Protein complexes are central in the yeast genetic landscape

If perturbing two genes together has a stronger or weaker effect than expected, they are said to genetically interact. Genetic interactions are important because they help map gene function, and functionally related genes have similar genetic interaction patterns. Mapping quantitative (positive and negative) genetic interactions on a global scale has recently become possible. This data clearly shows groups of genes connected by predominantly positive or negative interactions, termed monochromatic groups. These groups often correspond to functional modules, like biological processes or complexes, or connections between modules. However it is not yet known how these patterns globally relate to known functional modules. Here we systematically study the monochromatic nature of known biological processes using the largest quantitative genetic interaction data set available, which includes fitness measurements for ～5.4 million gene pairs in the yeast Saccharomyces cerevisiae. We find that only 10% of biological processes, as defined by Gene Ontology annotations, and less than 1% of inter-process connections are monochromatic. Further, we show that protein complexes are responsible for a surprisingly large fraction of these patterns. This suggests that complexes play a central role in shaping the monochromatic landscape of biological processes. Altogether this work shows that both positive and negative monochromatic patterns are found in known biological processes and in their connections and that protein complexes play an important role in these patterns. The monochromatic processes, complexes and connections we find chart a hierarchical and modular map of sensitive and redundant biological systems in the yeast cell that will be useful for gene function prediction and comparison across phenotypes and organisms. Furthermore the analysis methods we develop are applicable to other species for which genetic interactions will progressively become more available.     

Using pathway modules as targets for assay development in xenobiotic screening

Toxicology and pharmaceutical research is increasingly making use of high throughout-screening (HTS) methods to assess the effects of chemicals on molecular pathways, cells and tissues. Whole-genome microarray analysis provides broad information on the response of biological systems to chemical exposure, but is not practical to use when thousands of chemicals need to be evaluated at multiple doses and time points, as well as across different tissues, species and life-stages. A useful alternative approach is to identify a focused set of genes that can give a coarse picture of systems-level responses and that can be scaled to the evaluation of thousands of chemicals and diverse biological contexts. We demonstrate a computational approach to select in vitro expression assay targets that are informative and broadly distributed in biological pathway space, using the concept of pathway modularity. Canonical pathways are decomposed into subnetworks (modules) of functionally-related genes based on rules such as co-regulated expression, protein-protein interactions, and coordinated physiological activity. Pathway modules are constructed using these rules but are then restricted by the bounds of canonical pathways. We demonstrate this approach using a subset of genes associated with tumor development and cancer progression. Target genes were identified for assay development, and then validated by using independent, published microarray data. The result is a targeted set of genes that are sensitive predictors of whether a chemical will perturb each pathway module. These selected genes could then form the basis for a battery to test for pathway-chemical interactions under many biological contexts using throughput expression-based assays.     

Novel insights into DAPK autophagic signalling using peptide aptamer combinatorial protein-interaction screens

DAPK represents a relatively unique enzyme in the protein kinase superfamily whose major biological functions are linked to both autophagy and signal-mediated apoptosis. However, genetic studies have not yet uncovered how DAPK integrates into the core autophagy-related (Atg) machinery since DAPK is not present in a genetically tractable eukaryotic cell such as yeast. Furthermore, there have been no definitive DAPK binding proteins identified in metazoan systems that play a direct role in cooperating with DAPK in autophagy. We have utilized a growing concept in systems biology that invokes linear peptide-motifs as a fundamental mechanism driving protein-protein interactions and as a key switch underlying the dynamics of a signal transduction pathway. By using peptide combinatorial libraries as an assay that reflects the diversity of the linear peptide motif repertoire in the mammalian proteome, we identified microtubule-associated protein 1B (MAP1B) as a novel DAPK interacting protein that stimulates DAPK-dependent membrane blebbing and autophagy. MAP1B has previously been shown to form a functional interaction with the autophagosomal protein Atg8 (LC3). Together these studies define a genetic interaction between DAPK-MAP1B in the regulation of autophagy that may have particular relevance to cellular signalling pathways that regulate cell survival or cell death under distinct environmental stresses.     

Evaluation of Physical and Functional Protein-Protein Interaction Prediction Methods for Detecting Biological Pathways

Background

Cellular activities are governed by the physical and the functional interactions among several proteins involved in various biological pathways. With the availability of sequenced genomes and high-throughput experimental data one can identify genome-wide protein-protein interactions using various computational techniques. Comparative assessments of these techniques in predicting protein interactions have been frequently reported in the literature but not their ability to elucidate a particular biological pathway.

Methods

Towards the goal of understanding the prediction capabilities of interactions among the specific biological pathway proteins, we report the analyses of 14 biological pathways of Escherichia coli catalogued in KEGG database using five protein-protein functional linkage prediction methods. These methods are phylogenetic profiling, gene neighborhood, co-presence of orthologous genes in the same gene clusters, a mirrortree variant, and expression similarity.

Conclusions

Our results reveal that the prediction of metabolic pathway protein interactions continues to be a challenging task for all methods which possibly reflect flexible/independent evolutionary histories of these proteins. These methods have predicted functional associations of proteins involved in amino acids, nucleotide, glycans and vitamins & co-factors pathways slightly better than the random performance on carbohydrate, lipid and energy metabolism. We also make similar observations for interactions involved among the environmental information processing proteins. On the contrary, genetic information processing or specialized processes such as motility related protein-protein linkages that occur in the subset of organisms are predicted with comparable accuracy. Metabolic pathways are best predicted by using neighborhood of orthologous genes whereas phyletic pattern is good enough to reconstruct central dogma pathway protein interactions. We have also shown that the effective use of a particular prediction method depends on the pathway under investigation. In case one is not focused on specific pathway, gene expression similarity method is the best option.     

Protein interaction networks in plants

Protein–protein interactions are fundamental to virtually every aspect of cellular functions. With the development of high-throughput technologies of both the yeast two-hybrid system and tandem mass spectrometry, genome-wide protein-linkage mapping has become a major objective in post-genomic research. While at least partial “interactome” networks of several model organisms are already available, in the plant field, progress in this respect is slow. However, even with comprehensive protein interaction data still missing, substantial recent advance in the graph-theoretical functional interpretation of complex network architectures might pave the way for novel approaches in plant research. This article reviews current progress and discussions in network biology. Emphasis is put on the question of what can be learned about protein functions and cellular processes by studying the topology of complex protein interaction networks and the evolutionary mechanisms underlying their development. Particularly the intermediate and local levels of network organization—the modules, motifs and cliques—are increasingly recognized as the operational units of biological functions. As demonstrated by some recent results from systematic analyses of plant protein families, protein interaction networks promise to be a valuable tool for a molecular understanding of functional specificities and for identifying novel regulatory components and pathways.     

Computational Prediction of Protein–Protein Interaction Networks: Algo-rithms and Resources

Protein interactions play an important role in the discovery of protein functions and pathways in biological processes. This is especially true in case of the diseases caused by the loss of specific protein-protein interactions in the organism. The accuracy of experimental results in finding protein-protein interactions, however, is rather dubious and high throughput experimental results have shown both high false positive beside false negative information for protein interaction. Computational methods have attracted tremendous attention among biologists because of the ability to predict protein-protein interactions and validate the obtained experimental results. In this study, we have reviewed several computational methods for protein-protein interaction prediction as well as describing major databases, which store both predicted and detected protein-protein interactions, and the tools used for analyzing protein interaction networks and improving protein-protein interaction reliability.     

Probabilistic model of the human protein-protein interaction network

A catalog of all human protein-protein interactions would provide scientists with a framework to study protein deregulation in complex diseases such as cancer. Here we demonstrate that a probabilistic analysis integrating model organism interactome data, protein domain data, genome-wide gene expression data and functional annotation data predicts nearly 40,000 protein-protein interactions in humans-a result comparable to those obtained with experimental and computational approaches in model organisms. We validated the accuracy of the predictive model on an independent test set of known interactions and also experimentally confirmed two predicted interactions relevant to human cancer, implicating uncharacterized proteins into definitive pathways. We also applied the human interactome network to cancer genomics data and identified several interaction subnetworks activated in cancer. This integrative analysis provides a comprehensive framework for exploring the human protein interaction network.