Differential network biology

Network-based inference framework for identifying cancer genes from gene expression data

Great efforts have been devoted to alleviate uncertainty of detected cancer genes as accurate identification of oncogenes is of tremendous significance and helps unravel the biological behavior of tumors. In this paper, we present a differential network-based framework to detect biologically meaningful cancer-related genes. Firstly, a gene regulatory network construction algorithm is proposed, in which a boosting regression based on likelihood score and informative prior is employed for improving accuracy of identification. Secondly, with the algorithm, two gene regulatory networks are constructed from case and control samples independently. Thirdly, by subtracting the two networks, a differential-network model is obtained and then used to rank differentially expressed hub genes for identification of cancer biomarkers. Compared with two existing gene-based methods (t-test and lasso), the method has a significant improvement in accuracy both on synthetic datasets and two real breast cancer datasets. Furthermore, identified six genes (TSPYL5, CD55, CCNE2, DCK, BBC3, and MUC1) susceptible to breast cancer were verified through the literature mining, GO analysis, and pathway functional enrichment analysis. Among these oncogenes, TSPYL5 and CCNE2 have been already known as prognostic biomarkers in breast cancer, CD55 has been suspected of playing an important role in breast cancer prognosis from literature evidence, and other three genes are newly discovered breast cancer biomarkers. More generally, the differential-network schema can be extended to other complex diseases for detection of disease associated-genes.

The effect of network biology on drug toxicology

INTRODUCTION

The high failure rate of drug candidates due to toxicity, during clinical trials, is a critical issue in drug discovery. Network biology has become a promising approach, in this regard, using the increasingly large amount of biological and chemical data available and combining it with bioinformatics. With this approach, the assessment of chemical safety can be done across multiple scales of complexity from molecular to cellular and system levels in human health. Network biology can be used at several levels of complexity.

AREAS COVERED

This review describes the strengths and limitations of network biology. The authors specifically assess this approach across different biological scales when it is applied to toxicity.

EXPERT OPINION

There has been much progress made with the amount of data that is generated by various omics technologies. With this large amount of useful data, network biology has the opportunity to contribute to a better understanding of a drug's safety profile. The authors believe that considering a drug action and protein's function in a global physiological environment may benefit our understanding of the impact some chemicals have on human health and toxicity. The next step for network biology will be to better integrate differential and quantitative data.

A network integration approach to predict conserved regulators related to pathogenicity of influenza and SARS-CoV respiratory viruses

Respiratory infections stemming from influenza viruses and the Severe Acute Respiratory Syndrome corona virus (SARS-CoV) represent a serious public health threat as emerging pandemics. Despite efforts to identify the critical interactions of these viruses with host machinery, the key regulatory events that lead to disease pathology remain poorly targeted with therapeutics. Here we implement an integrated network interrogation approach, in which proteome and transcriptome datasets from infection of both viruses in human lung epithelial cells are utilized to predict regulatory genes involved in the host response. We take advantage of a novel “crowd-based” approach to identify and combine ranking metrics that isolate genes/proteins likely related to the pathogenicity of SARS-CoV and influenza virus. Subsequently, a multivariate regression model is used to compare predicted lung epithelial regulatory influences with data derived from other respiratory virus infection models. We predicted a small set of regulatory factors with conserved behavior for consideration as important components of viral pathogenesis that might also serve as therapeutic targets for intervention. Our results demonstrate the utility of integrating diverse ‘omic datasets to predict and prioritize regulatory features conserved across multiple pathogen infection models.

Disturbed protein-protein interaction networks in metastatic melanoma are associated with worse prognosis and increased functional mutation burden

For disseminated melanoma, new prognostic biomarkers and therapeutic targets are urgently needed. The organization of protein-protein interaction networks was assessed via the transcriptomes of four independent studies of metastatic melanoma and related to clinical outcome and MAP-kinase pathway mutations (BRAF/NRAS). We also examined patient outcome-related differences in a predicted network of microRNAs and their targets. The 32 hub genes with the most reproducible survival-related disturbances in co-expression with their protein partner genes included oncogenes and tumor suppressors, previously known correlates of prognosis, and other proteins not previously associated with melanoma outcome. Notably, this network-based gene set could classify patients according to clinical outcomes with 67-80% accuracy among cohorts. Reproducibly disturbed networks were also more likely to have a higher functional mutation burden than would be expected by chance. The disturbed regions of networks are therefore markers of clinically relevant, selectable tumor evolution in melanoma which may carry driver mutations.

Mapping molecular association networks of nervous system diseases via large-scale analysis of published research

Network medicine has been applied successfully to elicit the structure of large-scale molecular interaction networks. Its main proponents have claimed that this approach to integrative medical investigation should make it possible to identify functional modules of interacting molecular biological units as well as interactions themselves. This paper takes a significant step in this direction. Based on a large-scale analysis of the nervous system molecular medicine literature, this study analyzes and visualizes the complex structure of associations between diseases on the one hand and all types of molecular substances on the other. From this analysis it then identifies functional co-association groups consisting of several types of molecular substances, each consisting of substances that exhibit a pattern of frequent co-association with similar diseases. These groups in turn exhibit interlinking in a complex pattern, suggesting that such complex interactions between functional molecular modules may play a role in disease etiology. We find that the patterns exhibited by the networks of disease - molecular substance associations studied here correspond well to a number of recently published research results, and that the groups of molecular substances identified by statistical analysis of these networks do appear to be interesting groups of molecular substances that are interconnected in identifiable and interpretable ways. Our results not only demonstrate that networks are a convenient framework to analyze and visualize large-scale, complex relationships among molecular networks and diseases, but may also provide a conceptual basis for bridging gaps in experimental and theoretical knowledge.

Cancer systems biology in the genome sequencing era: part 1, dissecting and modeling of tumor clones and their networks

Recent tumor genome sequencing confirmed that one tumor often consists of multiple cell subpopulations (clones) which bear different, but related, genetic profiles such as mutation and copy number variation profiles. Thus far, one tumor has been viewed as a whole entity in cancer functional studies. With the advances of genome sequencing and computational analysis, we are able to quantify and computationally dissect clones from tumors, and then conduct clone-based analysis. Emerging technologies such as single-cell genome sequencing and RNA-Seq could profile tumor clones. Thus, we should reconsider how to conduct cancer systems biology studies in the genome sequencing era. We will outline new directions for conducting cancer systems biology by considering that genome sequencing technology can be used for dissecting, quantifying and genetically characterizing clones from tumors. Topics discussed in Part 1 of this review include computationally quantifying of tumor subpopulations; clone-based network modeling, cancer hallmark-based networks and their high-order rewiring principles and the principles of cell survival networks of fast-growing clones.

A systematic approach for the genetic dissection of protein complexes in living cells

Cells contain many important protein complexes involved in performing and regulating structural, metabolic, and signaling functions. One major challenge in cell biology is to elucidate the organization and mechanisms of robustness of these complexes in vivo. We developed a systematic approach to study structural dependencies within complexes in living cells by deleting subunits and measuring pairwise interactions among other components. We used our methodology to perturb two conserved eukaryotic complexes: the retromer and the nuclear pore complex. Our results identify subunits that are critical for the assembly of these complexes, reveal their structural architecture, and uncover mechanisms by which protein interactions are modulated. Our results also show that paralogous proteins play a key role in the robustness of protein complexes and shape their assembly landscape. Our approach paves the way for studying the response of protein interactomes to mutations and enhances our understanding of genotype-phenotype maps.

Differential network analysis for the identification of condition-specific pathway activity and regulation

MOTIVATION

Identification of differential expressed genes has led to countless new discoveries. However, differentially expressed genes are only a proxy for finding dysregulated pathways. The problem is to identify how the network of regulatory and physical interactions rewires in different conditions or in disease.

RESULTS

We developed a procedure named DINA (DIfferential Network Analysis), which is able to identify set of genes, whose co-regulation is condition-specific, starting from a collection of condition-specific gene expression profiles. DINA is also able to predict which transcription factors (TFs) may be responsible for the pathway condition-specific co-regulation. We derived 30 tissue-specific gene networks in human and identified several metabolic pathways as the most differentially regulated across the tissues. We correctly identified TFs such as Nuclear Receptors as their main regulators and demonstrated that a gene with unknown function (YEATS2) acts as a negative regulator of hepatocyte metabolism. Finally, we showed that DINA can be used to make hypotheses on dysregulated pathways during disease progression. By analyzing gene expression profiles across primary and transformed hepatocytes, DINA identified hepatocarcinoma-specific metabolic and transcriptional pathway dysregulation.

AVAILABILITY

We implemented an on-line web-tool http://dina.tigem.it enabling the user to apply DINA to identify tissue-specific pathways or gene signatures.

CONTACT

dibernardo@tigem.it

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Impact of natural genetic variation on gene expression dynamics

DNA sequence variation causes changes in gene expression, which in turn has profound effects on cellular states. These variations affect tissue development and may ultimately lead to pathological phenotypes. A genetic locus containing a sequence variation that affects gene expression is called an “expression quantitative trait locus” (eQTL). Whereas the impact of cellular context on expression levels in general is well established, a lot less is known about the cell-state specificity of eQTL. Previous studies differed with respect to how “dynamic eQTL” were defined. Here, we propose a unified framework distinguishing static, conditional and dynamic eQTL and suggest strategies for mapping these eQTL classes. Further, we introduce a new approach to simultaneously infer eQTL from different cell types. By using murine mRNA expression data from four stages of hematopoiesis and 14 related cellular traits, we demonstrate that static, conditional and dynamic eQTL, although derived from the same expression data, represent functionally distinct types of eQTL. While static eQTL affect generic cellular processes, non-static eQTL are more often involved in hematopoiesis and immune response. Our analysis revealed substantial effects of individual genetic variation on cell type-specific expression regulation. Among a total number of 3,941 eQTL we detected 2,729 static eQTL, 1,187 eQTL were conditionally active in one or several cell types, and 70 eQTL affected expression changes during cell type transitions. We also found evidence for feedback control mechanisms reverting the effect of an eQTL specifically in certain cell types. Loci correlated with hematological traits were enriched for conditional eQTL, thus, demonstrating the importance of conditional eQTL for understanding molecular mechanisms underlying physiological trait variation. The classification proposed here has the potential to streamline and unify future analysis of conditional and dynamic eQTL as well as many other kinds of QTL data.

Complex physiological traits are affected through subtle changes of molecular traits like gene expression in the relevant tissues, which in turn are caused by genetic variation. A genetic locus containing a sequence variation affecting gene expression is called an expression quantitative trait locus (eQTL). Understanding the tissue and cell type specificity of eQTL effects is essential for revealing the molecular mechanisms underlying disease phenotypes. However, so far the cell-state dependence of eQTL is poorly understood. In order to systematically assess the importance of cell state-specific eQTL, we propose to distinguish static, conditional and dynamic eQTL and suggest strategies for mapping these eQTL classes. We applied our framework to mouse gene expression data from four hematopoietic stages and related cellular traits. The different eQTL classes, although derived from the same expression data, represent functionally distinct types of eQTL. Importantly, conditional eQTL are well correlated with relevant hematological traits. These findings emphasize the condition specificity of many regulatory relationships, even if the conditions under study are related. This calls for due caution when transferring conclusions about regulatory mechanisms across cell types or tissues. The proposed classification will also help to unravel dynamic behaviors in many other kinds of QTL data.

Systematic triple-mutant analysis uncovers functional connectivity between pathways involved in chromosome regulation

Genetic interactions reveal the functional relationships between pairs of genes. In this study, we describe a method for the systematic generation and quantitation of triple mutants, termed triple-mutant analysis (TMA). We have used this approach to interrogate partially redundant pairs of genes in S. cerevisiae, including ASF1 and CAC1, two histone chaperones. After subjecting asf1Δ cac1Δ to TMA, we found that the Swi/Snf Rdh54 protein compensates for the absence of Asf1 and Cac1. Rdh54 more strongly associates with the chromatin apparatus and the pericentromeric region in the double mutant. Moreover, Asf1 is responsible for the synthetic lethality observed in cac1Δ strains lacking the HIRA-like proteins. A similar TMA was carried out after deleting both CLB5 and CLB6, cyclins that regulate DNA replication, revealing a strong functional connection to chromosome segregation. This approach can reveal functional redundancies that cannot be uncovered through traditional double-mutant analyses.

Structure and dynamics of molecular networks: a novel paradigm of drug discovery: a comprehensive review

Despite considerable progress in genome- and proteome-based high-throughput screening methods and in rational drug design, the increase in approved drugs in the past decade did not match the increase of drug development costs. Network description and analysis not only give a systems-level understanding of drug action and disease complexity, but can also help to improve the efficiency of drug design. We give a comprehensive assessment of the analytical tools of network topology and dynamics. The state-of-the-art use of chemical similarity, protein structure, protein-protein interaction, signaling, genetic interaction and metabolic networks in the discovery of drug targets is summarized. We propose that network targeting follows two basic strategies. The "central hit strategy" selectively targets central nodes/edges of the flexible networks of infectious agents or cancer cells to kill them. The "network influence strategy" works against other diseases, where an efficient reconfiguration of rigid networks needs to be achieved by targeting the neighbors of central nodes/edges. It is shown how network techniques can help in the identification of single-target, edgetic, multi-target and allo-network drug target candidates. We review the recent boom in network methods helping hit identification, lead selection optimizing drug efficacy, as well as minimizing side-effects and drug toxicity. Successful network-based drug development strategies are shown through the examples of infections, cancer, metabolic diseases, neurodegenerative diseases and aging. Summarizing >1200 references we suggest an optimized protocol of network-aided drug development, and provide a list of systems-level hallmarks of drug quality. Finally, we highlight network-related drug development trends helping to achieve these hallmarks by a cohesive, global approach.

Identifiability and inference of pathway motifs by epistasis analysis

The accuracy of genetic network inference is limited by the assumptions used to determine if one hypothetical model is better than another in explaining experimental observations. Most previous work on epistasis analysis-in which one attempts to infer pathway relationships by determining equivalences among traits following mutations-has been based on Boolean or linear models. Here, we delineate the ultimate limits of epistasis-based inference by systematically surveying all two-gene network motifs and use symbolic algebra with arbitrary regulation functions to examine trait equivalences. Our analysis divides the motifs into equivalence classes, where different genetic perturbations result in indistinguishable experimental outcomes. We demonstrate that this partitioning can reveal important information about network architecture, and show, using simulated data, that it greatly improves the accuracy of genetic network inference methods. Because of the minimal assumptions involved, equivalence partitioning has broad applicability for gene network inference.

Multidimensional control of cell structural robustness

Ample adaptive and functional opportunities of a living cell are determined by the complexity of its structural organisation. However, such complexity gives rise to a problem of maintenance of the coherence of inner processes in macroscopic interims and in macroscopic volumes which is necessary to support the structural robustness of a cell. The solution to this problem lies in multidimensional control of the adaptive and functional changes of a cell as well as its self-renewing processes in the context of environmental conditions. Six mechanisms (principles) form the basis of this multidimensional control: regulatory circuits with feedback loops, redundant inner diversity within a cell, multilevel distributed network organisation of a cell, molecular selection within a cell, continuous informational flows and functioning with a reserve of power. In the review we provide detailed analysis of these mechanisms, discuss their specific functions and the role of the superposition of these mechanisms in the maintenance of cell structural robustness in a wide range of environmental conditions.

Teaching the bioinformatics of signaling networks: an integrated approach to facilitate multi-disciplinary learning

The number of bioinformatics tools and resources that support molecular and cell biology approaches is continuously expanding. Moreover, systems and network biology analyses are accompanied more and more by integrated bioinformatics methods. Traditional information-centered university teaching methods often fail, as (1) it is impossible to cover all existing approaches in the frame of a single course, and (2) a large segment of the current bioinformation can become obsolete in a few years. Signaling network offers an excellent example for teaching bioinformatics resources and tools, as it is both focused and complex at the same time. Here, we present an outline of a university bioinformatics course with four sample practices to demonstrate how signaling network studies can integrate biochemistry, genetics, cell biology and network sciences. We show that several bioinformatics resources and tools, as well as important concepts and current trends, can also be integrated to signaling network studies. The research-type hands-on experiences we show enable the students to improve key competences such as teamworking, creative and critical thinking and problem solving. Our classroom course curriculum can be re-formulated as an e-learning material or applied as a part of a specific training course. The multi-disciplinary approach and the mosaic setup of the course have the additional benefit to support the advanced teaching of talented students.

Genome-scale modeling of human metabolism - a systems biology approach

Altered metabolism is linked to the appearance of various human diseases and a better understanding of disease-associated metabolic changes may lead to the identification of novel prognostic biomarkers and the development of new therapies. Genome-scale metabolic models (GEMs) have been employed for studying human metabolism in a systematic manner, as well as for understanding complex human diseases. In the past decade, such metabolic models - one of the fundamental aspects of systems biology - have started contributing to the understanding of the mechanistic relationship between genotype and phenotype. In this review, we focus on the construction of the Human Metabolic Reaction database, the generation of healthy cell type- and cancer-specific GEMs using different procedures, and the potential applications of these developments in the study of human metabolism and in the identification of metabolic changes associated with various disorders. We further examine how in silico genome-scale reconstructions can be employed to simulate metabolic flux distributions and how high-throughput omics data can be analyzed in a context-dependent fashion. Insights yielded from this mechanistic modeling approach can be used for identifying new therapeutic agents and drug targets as well as for the discovery of novel biomarkers. Finally, recent advancements in genome-scale modeling and the future challenge of developing a model of whole-body metabolism are presented. The emergent contribution of GEMs to personalized and translational medicine is also discussed.

Clinical cancer genome and precision medicine

Revolutionary sequencing technologies have changed biomedical research and life science exponentially. Revealing the whole landscape of causal somatic and inherited mutations underlying individual patient's cancer sample by whole-genome sequencing (WGS) and whole-exome sequencing (WES) can lead to not only a new mutations-based taxonomy of solid tumors (Stratton, Science 331:1553-1558, 2011). But also shapes a roadmap for precision medicine (Roychowdhury et al., Sci Transl Med 3:111ra121, 2011; Roukos, Expert Rev Mol Diagn 12:215-218, 2012; Mirnezami et al., N Engl J Med 366:489-491, 2012). This inevitable approach for personalized diagnostics in concert with free-falling genome sequencing costs raises now the question of applying next-generation sequencing (NGS) technology in the clinic. In the pragmatic clinical world and in contrast to innovative research, is NGS-based clinical evidence sufficient for decision-making on tailoring the best available treatment to the individual cancer patient?

An integrated proteomics and bioinformatics analyses of hepatitis B virus X interacting proteins and identification of a novel interactor apoA-I

HBx is well-known to be a multifunctional protein encoded by HBV and its biological functions are mainly dependent on pleiotropic protein-protein interactions (PPIs); however, the global mapping of HBx-interactome has not been established so far. Thus, in this study, we have identified 127 HBx-interacting proteins by a profound GST pull-down assay coupled with mass spectrometry, and constructed an HBx-interactome network and core apoA-I pathways with a series of bioinformatics approaches. One of the identified HBx-binding partners is apolipoprotein A-I (apoA-I), which has a specific role in lipid and cholesterol metabolism. The HBx-apoA-I protein interaction was confirmed by both GST pull-down and co-immunoprecipitation. The ectopic overexpression of apoA-I can lead to a significant inhibition on HBV secretion concomitant with the reduction of cellular cholesterol level. In addition, HBV can modulate the function of apoA-I through HBx which might interact with the 44-189 residues of apoA-I and result in dysfunction of apoA-I such as decreased self-association ability, increased carbonyl level and impaired lipid-binding ability. Our results demonstrate an integrated physical association of HBx and host proteins, especially a novel interactor apoA-I that may influence the HBV secretion, which would shed new light on exploring the complicated mechanisms of HBV manipulation on host cellular functions.

BIOLOGICAL SIGNIFICANCE

HBx is well-known to be a multifunctional protein encoded by HBV and its biological functions are mainly dependent on pleiotropic protein-protein interactions. Although a series of HBx-interacting proteins have been identified, a global characterization of HBx interactome has not been reported. In this study, we have identified a total of 127 HBx-interacting proteins by a profound GST pull-down assay coupled with mass spectrometry, and constructed an HBx-interactome network with a series of bioinformatics approaches. Our results demonstrate an integrated physical association of HBx and host proteins which may help us explore the complicated mechanisms of HBV manipulation on host cellular functions. In addition, we validated one of the identified HBx-binding partners, apolipoprotein A-I (apoA-I), which played a significant inhibitory effect on HBV secretion, indicating a crucial role of the HBx-apoA-I axis in HBV life cycle.

PodNet, a protein-protein interaction network of the podocyte

Interactions between proteins crucially determine cellular structure and function. Differential analysis of the interactome may help elucidate molecular mechanisms during disease development; however, this analysis necessitates mapping of expression data on protein-protein interaction networks. These networks do not exist for the podocyte; therefore, we built PodNet, a literature-based mouse podocyte network in Cytoscape format. Using database protein-protein interactions, we expanded PodNet to XPodNet with enhanced connectivity. In order to test the performance of XPodNet in differential interactome analysis, we examined podocyte developmental differentiation and the effect of cell culture. Transcriptomes of podocytes in 10 different states were mapped on XPodNet and analyzed with the Cytoscape plugin ExprEssence, based on the law of mass action. Interactions between slit diaphragm proteins are most significantly upregulated during podocyte development and most significantly downregulated in culture. On the other hand, our analysis revealed that interactions lost during podocyte differentiation are not regained in culture, suggesting a loss rather than a reversal of differentiation for podocytes in culture. Thus, we have developed PodNet as a valuable tool for differential interactome analysis in podocytes, and we have identified established and unexplored regulated interactions in developing and cultured podocytes.

Systems biology and "omics" tools: a cooperation for next-generation mycorrhizal studies

Omics tools constitute a powerful means of describing the complexity of plants and soil-borne microorganisms. Next generation sequencing technologies, coupled with emerging systems biology approaches, seem promising to represent a new strategy in the study of plant-microbe interactions. Arbuscular mycorrhizal fungi (AMF) are ubiquitous symbionts of plant roots, that provide their host with many benefits. However, as obligate biotrophs, AMF show a genetic, cellular and physiological complexity that makes the study of their biology as well as their effective agronomical exploitation rather difficult. Here, we speculate that the increasing availability of omics data on mycorrhiza and of computational tools that allow systems biology approaches represents a step forward in the understanding of arbuscular mycorrhizal symbiosis. Furthermore, the application of this study-perspective to agriculturally relevant model plants, such as tomato and rice, will lead to a better in-field exploitation of this beneficial symbiosis in the frame of low-input agriculture.

Systems immunogenetics of vaccines

Vaccines are the most cost effective public health measure for preventing viral infection and limiting epidemic spread within susceptible populations. However, the efficacy of current protective vaccines is highly variable, particularly in aging populations. In addition, there have been a number of challenges in the development of new vaccines due to a lack of detailed understanding of the immune correlates of protection. To identify the mechanisms underlying the variability of the immune response to vaccines, system-level tools need to be developed that will further our understanding of virus-host interactions and correlates of vaccine efficacy. This will provide critical information for rational vaccine design and allow the development of an analog to the "precision medicine" framework (already acknowledged as a powerful approach in medicine and therapeutics) to be applied to vaccinology.

From systems to structure: bridging networks and mechanism

There is a wide gap between the generation of large-scale biological data sets and more-detailed, structural and mechanistic studies. However, recent studies that explicitly combine data from systems and structural biological approaches are having a profound effect on our ability to predict how mutations and small molecules affect atomic-level mechanisms, disrupt systems-level networks, and ultimately lead to changes in organismal fitness. In fact, we argue that a shared framework for analysis of nonadditive genetic and thermodynamic responses to perturbations will accelerate the integration of reductionist and global approaches. A stronger bridge between these two areas will allow for a deeper and more-complete understanding of complex biological phenomenon and ultimately provide needed breakthroughs in biomedical research.

Networks in proteomics analysis of cancer

Proteomics provides direct biological information on proteins but is still a limited platform. Borrowing from genomics, its cancer-specific applications can be broadly categorized as (1) pure diagnostics, (2) biomarkers, (3) identification of root causes and (4) identification of cancer-specific network rewirings. Biological networks capture complex relationships between proteins and provide an appropriate means of contextualization. While playing significantly larger roles, especially in 1 and 3, progress in proteomics-specific network-based methods is lagging as compared to genomics. Rapid hardware advances and improvements in proteomic identification and quantification have given rise to much better quality data alongside advent of new network-based analysis methods. However, a tighter integration between analytics and hardware is still essential for network analysis to play more significant roles in proteomics analysis.

Canonical correlation analysis for RNA-seq co-expression networks

Digital transcriptome analysis by next-generation sequencing discovers substantial mRNA variants. Variation in gene expression underlies many biological processes and holds a key to unravelling mechanism of common diseases. However, the current methods for construction of co-expression networks using overall gene expression are originally designed for microarray expression data, and they overlook a large number of variations in gene expressions. To use information on exon, genomic positional level and allele-specific expressions, we develop novel component-based methods, single and bivariate canonical correlation analysis, for construction of co-expression networks with RNA-seq data. To evaluate the performance of our methods for co-expression network inference with RNA-seq data, they are applied to lung squamous cell cancer expression data from TCGA database and our bipolar disorder and schizophrenia RNA-seq study. The preliminary results demonstrate that the co-expression networks constructed by canonical correlation analysis and RNA-seq data provide rich genetic and molecular information to gain insight into biological processes and disease mechanism. Our new methods substantially outperform the current statistical methods for co-expression network construction with microarray expression data or RNA-seq data based on overall gene expression levels.

Current challenges and future potential of tomato breeding using omics approaches

As tomatoes are one of the most important vegetables in the world, improvements in the quality and yield of tomato are strongly required. For this purpose, omics approaches such as metabolomics and transcriptomics are used not only for basic research to understand relationships between important traits and metabolism but also for the development of next generation breeding strategies of tomato plants, because an increase in the knowledge improves the taste and quality, stress resistance and/or potentially health-beneficial metabolites and is connected to improvements in the biochemical composition of tomatoes. Such omics data can be applied to network analyses to potentially reveal unknown cellular regulatory networks in tomato plants. The high-quality tomato genome that was sequenced in 2012 will likely accelerate the application of omics strategies, including next generation sequencing for tomato breeding. In this review, we highlight the current studies of omics network analyses of tomatoes and other plant species, in particular, a gene coexpression network. Key applications of omics approaches are also presented as case examples to improve economically important traits for tomato breeding.

Over-expression of heat shock factor 1 phenocopies the effect of chronic inhibition of TOR by rapamycin and is sufficient to ameliorate Alzheimer's-like deficits in mice modeling the disease

Rapamycin, an inhibitor of target-of-rapamycin, extends lifespan in mice, possibly by delaying aging. We recently showed that rapamycin halts the progression of Alzheimer's (AD)-like deficits, reduces amyloid-beta (Aβ) and induces autophagy in the human amyloid precursor protein (PDAPP) mouse model. To delineate the mechanisms by which chronic rapamycin delays AD we determined proteomic signatures in brains of control- and rapamycin-treated PDAPP mice. Proteins with reported chaperone-like activity were overrepresented among proteins up-regulated in rapamycin-fed PDAPP mice and the master regulator of the heat-shock response, heat-shock factor 1, was activated. This was accompanied by the up-regulation of classical chaperones/heat shock proteins (HSPs) in brains of rapamycin-fed PDAPP mice. The abundance of most HSP mRNAs except for alpha B-crystallin, however, was unchanged, and the cap-dependent translation inhibitor 4E-BP was active, suggesting that increased expression of HSPs and proteins with chaperone activity may result from preferential translation of pre-existing mRNAs as a consequence of inhibition of cap-dependent translation. The effects of rapamycin on the reduction of Aβ, up-regulation of chaperones, and amelioration of AD-like cognitive deficits were recapitulated by transgenic over-expression of heat-shock factor 1 in PDAPP mice. These results suggest that, in addition to inducing autophagy, rapamycin preserves proteostasis by increasing chaperones. We propose that the failure of proteostasis associated with aging may be a key event enabling AD, and that chronic inhibition of target-of-rapamycin may delay AD by maintaining proteostasis in brain. Read the Editorial Highlight for this article on doi: 10.1111/jnc.12098.

Protein complex-based analysis framework for high-throughput data sets

Analysis of high-throughput data increasingly relies on pathway annotation and functional information derived from Gene Ontology. This approach has limitations, in particular for the analysis of network dynamics over time or under different experimental conditions, in which modules within a network rather than complete pathways might respond and change. We report an analysis framework based on protein complexes, which are at the core of network reorganization. We generated a protein complex resource for human, Drosophila, and yeast from the literature and databases of protein-protein interaction networks, with each species having thousands of complexes. We developed COMPLEAT (http://www.flyrnai.org/compleat), a tool for data mining and visualization for complex-based analysis of high-throughput data sets, as well as analysis and integration of heterogeneous proteomics and gene expression data sets. With COMPLEAT, we identified dynamically regulated protein complexes among genome-wide RNA interference data sets that used the abundance of phosphorylated extracellular signal-regulated kinase in cells stimulated with either insulin or epidermal growth factor as the output. The analysis predicted that the Brahma complex participated in the insulin response.

Concise review: new paradigms for Down syndrome research using induced pluripotent stem cells: tackling complex human genetic disease

Down syndrome (DS) is a complex developmental disorder with diverse pathologies that affect multiple tissues and organ systems. Clear mechanistic description of how trisomy of chromosome 21 gives rise to most DS pathologies is currently lacking and is limited to a few examples of dosage-sensitive trisomic genes with large phenotypic effects. The recent advent of cellular reprogramming technology offers a promising way forward, by allowing derivation of patient-derived human cell types in vitro. We present general strategies that integrate genomics technologies and induced pluripotent stem cells to identify molecular networks driving different aspects of DS pathogenesis and describe experimental approaches to validate the causal requirement of candidate network defects for particular cellular phenotypes. This overall approach should be applicable to many poorly understood complex human genetic diseases, whose pathogenic mechanisms might involve the combined effects of many genes.

Interpreting the regulatory genome: the genomics of transcription factor function in Drosophila melanogaster

Researchers have now had access to the fully sequenced Drosophila melanogaster genome for over a decade, and the sequenced genomes of 11 additional Drosophila species have been available for almost 5 years, with more species' genomes becoming available every year [Adams MD, Celniker SE, Holt RA, et al. The genome sequence of Drosophila melanogaster. Science 2000;287:2185-95; Clark AG, Eisen MB, Smith DR, et al. Evolution of genes and genomes on the Drosophila phylogeny. Nature 2007;450:203-18]. Although the best studied of the D. melanogaster transcription factors (TFs) were cloned before sequencing of the genome, the availability of sequence data promised to transform our understanding of TFs and gene regulatory networks. Sequenced genomes have allowed researchers to generate tools for high-throughput characterization of gene expression levels, genome-wide TF localization and analyses of evolutionary constraints on DNA elements across multiple species. With an estimated 700 DNA-binding proteins in the Drosophila genome, it will be many years before each potential sequence-specific TF is studied in detail, yet the last decade of functional genomics research has already impacted our view of gene regulatory networks and TF DNA recognition.

Analysis of protein-protein interactions in cross-talk pathways reveals CRKL protein as a novel prognostic marker in hepatocellular carcinoma

Deciphering the network of signaling pathways in cancer via protein-protein interactions (PPIs) at the cellular level is a promising approach but remains incomplete. We used an in situ proximity ligation assay to identify and quantify 67 endogenous PPIs among 21 interlinked pathways in two hepatocellular carcinoma (HCC) cells, Huh7 (minimally migratory cells) and Mahlavu (highly migratory cells). We then applied a differential network biology analysis and determined that the novel interaction, CRKL-FLT1, has a high centrality ranking, and the expression of this interaction is strongly correlated with the migratory ability of HCC and other cancer cell lines. Knockdown of CRKL and FLT1 in HCC cells leads to a decrease in cell migration via ERK signaling and the epithelial-mesenchymal transition process. Our immunohistochemical analysis shows high expression levels of the CRKL and CRKL-FLT1 pair that strongly correlate with reduced disease-free and overall survival in HCC patient samples, and a multivariate analysis further established CRKL and the CRKL-FLT1 as novel prognosis markers. This study demonstrated that functional exploration of a disease network with interlinked pathways via PPIs can be used to discover novel biomarkers.

Using an adaptive gene network model for self-organizing multicellular behavior

Using the transient interleukin (IL)-2 secretion of effector T helper (T(eff)) cells as an example, we show that self-organizing multicellular behavior can be modeled and predicted by an adaptive gene network model. Incorporating an adaptation algorithm we established previously, we construct a network model that has the parameter values iteratively updated to cope with environmental change governed by diffusion and cell-cell interactions. In contrast to non-adaptive models, we find that the proposed adaptive model for individual T(eff) cells can generate transient IL-2 secretory behavior that is observed experimentally at the population level. The proposed adaptive modeling approach can be a useful tool in the study of self-organizing behavior observed in other contexts in biology, including microbial pathogenesis, antibiotic resistance, embryonic development, tumor formation, etc.

Studying the differential co-expression of microRNAs reveals significant role of white matter in early Alzheimer's progression

MicroRNAs (miRNAs) are a class of short non-coding RNAs, which show tissue-specific regulatory activity on genes. Expression profiling of miRNAs is an important step for understanding the pathology of Alzheimer's disease (AD), a neurodegenerative disorder originating in the brain. Recent studies highlight that miRNAs enriched in gray matter (GM) and white matter (WM) of AD brains show differential expression. However, no in-depth study has yet been conducted on analysing the differential co-expression of pairs of miRNAs over GM and WM. Two genes (or miRNAs) are said to be co-expressed if their expression profiles change similarly over a number of samples. A pair of co-expressed genes under a condition type (or phenotype) may not remain co-expressed, or get contra-expressed, under another condition. Such pairs of genes are referred to as differentially co-expressed. Such an investigation in the early stage of AD is reported in this article. A network of differentially co-expressed miRNAs in GM and WM is first built. Analysis of the differential co-expression property reveals that such a network can not have any cycle. We use the notion of switching to distinguish two distinct types of differential co-expression patterns - a pair of miRNAs that are highly co-expressed in GM but does not remain so in WM, and vice versa. Based on this, we find the substructures, referred to as differentially co-expressed switching tree (DCST), that throughout have similar pattern of switching. The miR-423-5p emerges as a hub of the network. We extract subtrees of these DCSTs that have similar switching pattern throughout. These substructures are found to be both statistically and biologically significant. A large number of miRNAs obtained from the DCSTs are found to have association with AD, most of which are enriched in WM. This computational study therefore indicates a significant role of WM in early AD progression, a hitherto less acknowledged fact.

Ripples in the pond--using a systems approach to decipher the cellular functions of membrane microdomains

Membrane microdomains such as lipid rafts and caveolae regulate a myriad of cellular functions including cell signalling, protein trafficking, cell viability, and cell movement. They have been implicated in diseases such as cancer, diabetes and Alzheimer's disease, highlighting the essential role they play in cell processes. Despite much research and debate on the size, composition and dynamics of membrane microdomains, the molecular mechanism(s) of their action remain poorly understood. Most studies have dealt solely with the content and properties of the membrane microdomain as an entity in itself. However, recent work shows that membrane microdomain disruption has wide ranging effects on other subcellular compartments, and the cell as a whole. Hence we propose that a systems approach incorporating many cellular attributes such as subcellular localisation is required in order to understand the global impact of microdomains on cell function. Although analysis of sub-proteome changes already provides additional insight, we further propose biological network analysis of functional proteomics data to capture effects at the systems level. In this review, we highlight the use of protein-protein interactions networks and mixed networks to portray and visualize the relationships between proteins within and between subcellular fractions. Such a systems analysis will be required to improve our understanding of the full cellular function of membrane microdomains.

Differential RNAi screening provides insights into the rewiring of signalling networks during oxidative stress

Reactive Oxygen Species (ROS) are a natural by-product of cellular growth and proliferation, and are required for fundamental processes such as protein-folding and signal transduction. However, ROS accumulation, and the onset of oxidative stress, can negatively impact cellular and genomic integrity. Signalling networks have evolved to respond to oxidative stress by engaging diverse enzymatic and non-enzymatic antioxidant mechanisms to restore redox homeostasis. The architecture of oxidative stress response networks during periods of normal growth, and how increased ROS levels dynamically reconfigure these networks are largely unknown. In order to gain insight into the structure of signalling networks that promote redox homeostasis we first performed genome-scale RNAi screens to identify novel suppressors of superoxide accumulation. We then infer relationships between redox regulators by hierarchical clustering of phenotypic signatures describing how gene inhibition affects superoxide levels, cellular viability, and morphology across different genetic backgrounds. Genes that cluster together are likely to act in the same signalling pathway/complex and thus make "functional interactions". Moreover we also calculate differential phenotypic signatures describing the difference in cellular phenotypes following RNAi between untreated cells and cells submitted to oxidative stress. Using both phenotypic signatures and differential signatures we construct a network model of functional interactions that occur between components of the redox homeostasis network, and how such interactions become rewired in the presence of oxidative stress. This network model predicts a functional interaction between the transcription factor Jun and the IRE1 kinase, which we validate in an orthogonal assay. We thus demonstrate the ability of systems-biology approaches to identify novel signalling events.

Visualization and analysis of biological networks

The study of the interactome-the totality of the protein-protein interactions taking place in a cell-has experienced an enormous growth in the last few years. Biological networks representation and analysis has become an everyday tool for many biologists and bioinformatics, as these interaction graphs allow us to map and characterize signaling pathways and predict the function of unknown proteins. However, given the size and complexity of interactome datasets, extracting meaningful information from interaction networks can be a daunting task. Many different tools and approaches can be used to build, represent, and analyze biological networks. In this chapter, we will use a practical example to guide novice users through this process. We will be making use of the popular open source tool Cytoscape and of other resources such as : the PSICQUIC client to access several protein interaction repositories and the BiNGO plugin to perform GO enrichment analysis of the resulting network.

From a biological hypothesis to the construction of a mathematical model

Mathematical models serve to explain complex biological phenomena and provide predictions that can be tested experimentally. They can provide plausible scenarios of a complex biological behavior when intuition is not sufficient anymore. The process from a biological hypothesis to a mathematical model might be challenging for biologists that are not familiar with mathematical modeling. In this chapter we discuss a possible workflow that describes the steps to be taken starting from a biological hypothesis on a biochemical cellular mechanism to the construction of a mathematical model using the appropriate formalism. An important part of this workflow is formalization of biological knowledge, which can be facilitated by existing tools and standards developed by the systems biology community. This chapter aims at introducing modeling to experts in molecular biology that would like to convert their hypotheses into mathematical models.

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.