Differential network biology

Genome scale transcriptional response diversity among ten ecotypes of Arabidopsis thaliana during heat stress

In the scenario of global warming and climate change, heat stress is a serious threat to crop production worldwide. Being sessile, plants cannot escape from heat. Plants have developed various adaptive mechanisms to survive heat stress. Several studies have focused on diversity of heat tolerance levels in divergent Arabidopsis thaliana (A. thaliana) ecotypes, but comprehensive genome scale understanding of heat stress response in plants is still lacking. Here we report the genome scale transcript responses to heat stress of 10 A. thaliana ecotypes (Col, Ler, C24, Cvi, Kas1, An1, Sha, Kyo2, Eri, and Kond) originated from different geographical locations. During the experiment, A. thaliana plants were subjected to heat stress (38°C) and transcript responses were monitored using Arabidopsis NimbleGen ATH6 microarrays. The responses of A. thaliana ecotypes exhibited considerable variation in the transcript abundance levels. In total, 3644 transcripts were significantly heat regulated (p < 0.01) in the 10 ecotypes, including 244 transcription factors and 203 transposable elements. By employing a systems genetics approach- Network Component Analysis (NCA), we have constructed an in silico transcript regulatory network model for 35 heat responsive transcription factors during cellular responses to heat stress in A. thaliana. The computed activities of the 35 transcription factors showed ecotype specific responses to the heat treatment.

Large Scale Chemical Cross-linking Mass Spectrometry Perspectives

The spectacular heterogeneity of a complex protein mixture from biological samples becomes even more difficult to tackle when one’s attention is shifted towards different protein complex topologies, transient interactions, or localization of PPIs. Meticulous protein-by-protein affinity pull-downs and yeast-two-hybrid screens are the two approaches currently used to decipher proteome-wide interaction networks. Another method is to employ chemical cross-linking, which gives not only identities of interactors, but could also provide information on the sites of interactions and interaction interfaces. Despite significant advances in mass spectrometry instrumentation over the last decade, mapping Protein-Protein Interactions (PPIs) using chemical cross-linking remains time consuming and requires substantial expertise, even in the simplest of systems. While robust methodologies and software exist for the analysis of binary PPIs and also for the single protein structure refinement using cross-linking-derived constraints, undertaking a proteome-wide cross-linking study is highly complex. Difficulties include i) identifying cross-linkers of the right length and selectivity that could capture interactions of interest; ii) enrichment of the cross-linked species; iii) identification and validation of the cross-linked peptides and cross-linked sites.

In this review we examine existing literature aimed at the large-scale protein cross-linking and discuss possible paths for improvement. We also discuss short-length cross-linkers of broad specificity such as formaldehyde and diazirine-based photo-cross-linkers. These cross-linkers could potentially capture many types of interactions, without strict requirement for a particular amino-acid to be present at a given protein-protein interface. How these shortlength, broad specificity cross-linkers be applied to proteome-wide studies? We will suggest specific advances in methodology, instrumentation and software that are needed to make such a leap.

Metabolic flux-based modularity using shortest retroactive distances

Background

Graph-based modularity analysis has emerged as an important tool to study the functional organization of biological networks. However, few methods are available to study state-dependent changes in network modularity using biological activity data. We develop a weighting scheme, based on metabolic flux data, to adjust the interaction distances in a reaction-centric graph model of a metabolic network. The weighting scheme was combined with a hierarchical module assignment algorithm featuring the preservation of metabolic cycles to examine the effects of cellular differentiation and enzyme inhibitions on the functional organization of adipocyte metabolism.

Results

Our analysis found that the differences between various metabolic states primarily involved the assignment of two specific reactions in fatty acid synthesis and glycerogenesis. Our analysis also identified cyclical interactions between reactions that are robust with respect to metabolic state, suggesting possible co-regulation. Comparisons based on cyclical interaction distances between reaction pairs suggest that the modular organization of adipocyte metabolism is stable with respect to the inhibition of an enzyme, whereas a major physiological change such as cellular differentiation leads to a more substantial reorganization.

Conclusion

Taken together, our results support the notion that network modularity is influenced by both the connectivity of the network’s components as well as the relative engagements of the connections.

DiffCorr: an R package to analyze and visualize differential correlations in biological networks

Large-scale "omics" data, such as microarrays, can be used to infer underlying cellular regulatory networks in organisms, enabling us to better understand the molecular basis of disease and important traits. Correlation approaches, such as a hierarchical cluster analysis, have been widely used to analyze omics data. In addition to the changes in the mean levels of molecules in the omics data, it is important to know about the changes in the correlation relationship among molecules between 2 experimental conditions. The development of a tool to identify differential correlation patterns in omics data in an efficient and unbiased manner is therefore desirable. We developed the DiffCorr package, a simple method for identifying pattern changes between 2 experimental conditions in correlation networks, which builds on a commonly used association measure, such as Pearson's correlation coefficient. DiffCorr calculates correlation matrices for each dataset, identifies the first principal component-based "eigen-molecules" in the correlation networks, and tests differential correlation between the 2 groups based on Fisher's z-test. We illustrated its utility by demonstrating biologically relevant, differentially correlated molecules in transcriptome coexpression and metabolite-to-metabolite correlation networks. DiffCorr can explore differential correlations between 2 conditions in the context of post-genomics data types, namely transcriptomics and metabolomics. DiffCorr is simple to use in calculating differential correlations and is suitable for the first step towards inferring causal relationships and detecting biomarker candidates. The package can be downloaded from the following website: http://diffcorr.sourceforge.net/.

Network analysis: tackling complex data to study plant metabolism

Incomplete knowledge of biochemical pathways makes the holistic description of plant metabolism a non-trivial undertaking. Sensitive analytical platforms, which are capable of accurately quantifying the levels of the various molecular entities of the cell, can assist in tackling this task. However, the ever-increasing amount of high-throughput data, often from multiple technologies, requires significant computational efforts for integrative analysis. Here we introduce the application of network analysis to study plant metabolism and describe the construction and analysis of correlation-based networks from (time-resolved) metabolomics data. By investigating the interactions between metabolites, network analysis can help to interpret complex datasets through the identification of key network components. The relationship between structural and biological roles of network components can be evaluated and employed to aid metabolic engineering.

Next-generation proteomics: towards an integrative view of proteome dynamics

Next-generation sequencing allows the analysis of genomes, including those representing disease states. However, the causes of most disorders are multifactorial, and systems-level approaches, including the analysis of proteomes, are required for a more comprehensive understanding. The proteome is extremely multifaceted owing to splicing and protein modifications, and this is further amplified by the interconnectivity of proteins into complexes and signalling networks that are highly divergent in time and space. Proteome analysis heavily relies on mass spectrometry (MS). MS-based proteomics is starting to mature and to deliver through a combination of developments in instrumentation, sample preparation and computational analysis. Here we describe this emerging next generation of proteomics and highlight recent applications.

Cancer omics: from regulatory networks to clinical outcomes

Current limitation in cancer genomic studies is a lack of the integration of various omics data generated through next generation sequencing technologies, as well as a lack of the sounding and comprehensive epigenomic and genomic information about a particular cancer cell type. In this review, we will discuss main aspects of current genomics research with its application in cancer topics. We will first overview the next-generation sequencing technologies, then outline the major computational approaches, particularly focusing on ChIP-based omics data, and list several remaining open questions facing computational biologists, further present regulatory network analysis inferred from the ChIP-based omics data; finally implicate the clinical outcomes from the network and pathway analysis.

Differential network entropy reveals cancer system hallmarks

The cellular phenotype is described by a complex network of molecular interactions. Elucidating network properties that distinguish disease from the healthy cellular state is therefore of critical importance for gaining systems-level insights into disease mechanisms and ultimately for developing improved therapies. By integrating gene expression data with a protein interaction network we here demonstrate that cancer cells are characterised by an increase in network entropy. In addition, we formally demonstrate that gene expression differences between normal and cancer tissue are anticorrelated with local network entropy changes, thus providing a systemic link between gene expression changes at the nodes and their local correlation patterns. In particular, we find that genes which drive cell-proliferation in cancer cells and which often encode oncogenes are associated with reductions in network entropy. These findings may have potential implications for identifying novel drug targets.

qPCA: a scalable assay to measure the perturbation of protein-protein interactions in living cells

One of the most important challenges in systems biology is to understand how cells respond to genetic and environmental perturbations. Here we show that the yeast DHFR-PCA, coupled with high-resolution growth profiling (DHFR-qPCA), is a straightforward assay to study the modulation of protein-protein interactions (PPIs) in vivo as a response to genetic, metabolic and drug perturbations. Using the canonical Protein Kinase A (PKA) pathway as a test system, we show that changes in PKA activity can be measured in living cells as a modulation of the interaction between its regulatory (Bcy1) and catalytic (Tpk1 and Tpk2) subunits in response to changes in carbon metabolism, caffeine and methyl methanesulfonate (MMS) treatments and to modifications in the dosage of its enzymatic regulators, the phosphodiesterases. Our results show that the DHFR-qPCA is easily implementable and amenable to high-throughput. The DHFR-qPCA will pave the way to the study of the effects of drug, genetic and environmental perturbations on in vivo PPI networks, thus allowing the exploration of new spaces of the eukaryotic interactome.

A multidimensional matrix for systems biology research and its application to interaction networks

A multidimensional matrix containing 76 parameters from 21 transcriptomics, proteomics, interactomics, phenotypic and sequence-based data sets, in which each data set covered most of the Saccharomyces cerevisiae proteome, was compiled for systems biology research. The maximal information coefficient (MIC) was used to measure correlations between every pair of parameters. Out of 2850 possible comparisons, 340 pairs of variables (12%) showed statistically significant MIC scores. There were 321 relationships that were expected; these included relationships within physicochemical parameters of proteins, between abundance levels of genes/proteins and expression noise, and between different types of intracellular networks. We found 19 potentially novel relationships between different types of "-omics" data. The strongest of these involved genetic interaction networks, which were correlated with pleiotropy and cell-to-cell variability in protein expression. Protein disorder also showed a number of significant relationships with protein abundance, signaling and regulatory networks. Significant cross-talk was seen between the signaling and kinase interaction networks. Investigation of this revealed densely connected kinase clusters and significant signaling between them, along with signaling centers that act as integrators or broadcasters of intracellular information. These centers may allow for redundancy and a means of dampening noise in networks under a variety of genetic or environmental perturbations.

Structural and dynamical analysis of biological networks

Biological networks are currently being studied with approaches derived from the mathematical and physical sciences. Their structural analysis enables to highlight nodes with special properties that have sometimes been correlated with the biological importance of a gene or a protein. However, biological networks are dynamic both on the evolutionary time-scale, and on the much shorter time-scale of physiological processes. There is therefore no unique network for a given cellular process, but potentially many realizations, each with different properties as a consequence of regulatory mechanisms. Such realizations provide snapshots of a same network in different conditions, enabling the study of condition-dependent structural properties. True dynamical analysis can be obtained through detailed mathematical modeling techniques that are not easily scalable to full network models.

Perspectives on systems biology applications in diabetic kidney disease

Diabetic kidney disease (DKD) is a microvascular complication of type 1 and 2 diabetes with a devastating impact on individuals with the disease, their families, and society as a whole. DKD is the single most frequent cause of incident chronic kidney disease cases and accounts for over 40% of the population with end-stage renal disease. Contributing factors for the high prevalence are the increase in obesity and subsequent diabetes combined with an improved long-term survival with diabetes. Environment and genetic variations contribute to DKD susceptibility and progressive loss of kidney function. How the molecular mechanisms of genetic and environmental exposures interact during DKD initiation and progression is the focus of ongoing research efforts. The development of standardized, unbiased high-throughput profiling technologies of human DKD samples opens new avenues in capturing the multiple layers of DKD pathobiology. These techniques routinely interrogate analytes on a genome-wide scale generating comprehensive DKD-associated fingerprints. Linking the molecular fingerprints to deep clinical phenotypes may ultimately elucidate the intricate molecular interplay in a disease stage and subtype-specific manner. This insight will form the basis for accurate prognosis and facilitate targeted therapeutic interventions. In this review, we present ongoing efforts from large-scale data integration translating "-omics" research efforts into improved and individualized health care in DKD.

Integration of biological networks and pathways with genetic association studies

Millions of genetic variants have been assessed for their effects on the trait of interest in genome-wide association studies (GWAS). The complex traits are affected by a set of inter-related genes. However, the typical GWAS only examine the association of a single genetic variant at a time. The individual effects of a complex trait are usually small, and the simple sum of these individual effects may not reflect the holistic effect of the genetic system. High-throughput methods enable genomic studies to produce a large amount of data to expand the knowledge base of the biological systems. Biological networks and pathways are built to represent the functional or physical connectivity among genes. Integrated with GWAS data, the network- and pathway-based methods complement the approach of single genetic variant analysis, and may improve the power to identify trait-associated genes. Taking advantage of the biological knowledge, these approaches are valuable to interpret the functional role of the genetic variants, and to further understand the molecular mechanism influencing the traits. The network- and pathway-based methods have demonstrated their utilities, and will be increasingly important to address a number of challenges facing the mainstream GWAS.

Multiple genetic interaction experiments provide complementary information useful for gene function prediction

Genetic interactions help map biological processes and their functional relationships. A genetic interaction is defined as a deviation from the expected phenotype when combining multiple genetic mutations. In Saccharomyces cerevisiae, most genetic interactions are measured under a single phenotype - growth rate in standard laboratory conditions. Recently genetic interactions have been collected under different phenotypic readouts and experimental conditions. How different are these networks and what can we learn from their differences? We conducted a systematic analysis of quantitative genetic interaction networks in yeast performed under different experimental conditions. We find that networks obtained using different phenotypic readouts, in different conditions and from different laboratories overlap less than expected and provide significant unique information. To exploit this information, we develop a novel method to combine individual genetic interaction data sets and show that the resulting network improves gene function prediction performance, demonstrating that individual networks provide complementary information. Our results support the notion that using diverse phenotypic readouts and experimental conditions will substantially increase the amount of gene function information produced by genetic interaction screens.

IntScore: a web tool for confidence scoring of biological interactions

Knowledge of all molecular interactions that potentially take place in the cell is a key for a detailed understanding of cellular processes. Currently available interaction data, such as protein–protein interaction maps, are known to contain false positives that inevitably diminish the accuracy of network-based inferences. Interaction confidence scoring is thus a crucial intermediate step after obtaining interaction data and before using it in an interaction network-based inference approach. It enables to weight individual interactions according to the likelihood that they actually take place in the cell, and can be used to filter out false positives. We describe a web tool called IntScore which calculates confidence scores for user-specified sets of interactions. IntScore provides six network topology- and annotation-based confidence scoring methods. It also enables the integration of scores calculated by the different methods into an aggregate score using machine learning approaches. IntScore is user-friendly and extensively documented. It is freely available at http://intscore.molgen.mpg.de.

Network analysis has diverse roles in drug discovery

Computational biologists use network analysis to uncover relationships between various data types of interest for drug discovery. For example, signalling and metabolic pathways are commonly used to understand disease states and drug mechanisms. However, several other flavours of network analysis techniques are also applicable in a drug discovery context. Recent advances include networks that encompass relationships between diseases, molecular mechanisms and gene targets. Even social networks that mirror interactions within the scientific community are helping to foster collaborations and novel research. We review how these different types of network analysis approaches facilitate drug discovery and their associated challenges.

Analyzing protein-protein interaction networks

The advent of the "omics" era in biology research has brought new challenges and requires the development of novel strategies to answer previously intractable questions. Molecular interaction networks provide a framework to visualize cellular processes, but their complexity often makes their interpretation an overwhelming task. The inherently artificial nature of interaction detection methods and the incompleteness of currently available interaction maps call for a careful and well-informed utilization of this valuable data. In this tutorial, we aim to give an overview of the key aspects that any researcher needs to consider when working with molecular interaction data sets and we outline an example for interactome analysis. Using the molecular interaction database IntAct, the software platform Cytoscape, and its plugins BiNGO and clusterMaker, and taking as a starting point a list of proteins identified in a mass spectrometry-based proteomics experiment, we show how to build, visualize, and analyze a protein-protein interaction network.