Volume 17 Supplement 5
Inferring active regulatory networks from gene expression data using a combination of prior knowledge and enrichment analysis
© Chouvardas et al. 2016
Published: 6 June 2016
Under both physiological and pathological conditions gene expression programs are shaped through the interplay of regulatory proteins and their gene targets, interactions between which form intricate gene regulatory networks (GRN). While the assessment of genome-wide expression for the complete set of genes at a given condition has become rather straight-forward and is performed routinely, we are still far from being able to infer the topology of gene regulation simply by analyzing its “descendant” expression profile. In this work we are trying to overcome the existing limitations for the inference and study of such regulatory networks. We are combining our approach with state-of-the-art gene set enrichment analyses in order to create a tool, called Regulatory Network Enrichment Analysis (RNEA) that will prioritize regulatory and functional characteristics of a genome-wide expression experiment.
RNEA combines prior knowledge, originating from manual literature curation and small-scale experimental data, to construct a reference network of interactions and then uses enrichment analysis coupled with a two-level hierarchical parsing of the network, to infer the most relevant subnetwork for a given experimental setting. It is implemented as an R package, currently supporting human and mouse datasets and was herein tested on one test case for each of the two organisms. In both cases, RNEA’s gene set enrichment analysis was comparable to state-of-the-art methodologies. Moreover, through its distinguishing feature of regulatory subnetwork reconstruction, RNEA was able to define the key transcriptional regulators for the studied systems as supported from the literature.
RNEA constitutes a novel computational approach to obtain regulatory interactions directly from a genome-wide expression profile. Its simple implementation, with minimal requirements from the user is coupled with easy-to-parse enrichment lists and a subnetwork file that may be readily visualized to reveal the most important components of the regulatory hierarchy. The combination of prior information and novel concept of a hierarchical reconstruction of regulatory interactions makes RNEA a very useful tool for a first-level interpretation of gene expression profiles.
KeywordsTranscriptional regulation Gene regulatory network Gene set enrichment analysis
The advent of high-throughput genomics that started with DNA microarrays and is now rapidly shifting to next-generation-sequencing, has been producing a vast amount of information regarding a variety of cellular functions. In the context of gene expression measurements, genome-wide profiling approaches through RNASeq have made possible the monitoring of gene expression at unprecedented resolution, allowing not only the detection of genes present in the cell in only a few mRNA copies, but also revealing the transcriptional complexity reflected in the use of alternative transcript isoforms [1–3]. In this sense, the output of all genome-wide expression profiling approaches, summarized in lists of differentially expressed genes, may be seen as an accurate reflection of the intricate regulatory dynamics that reshape the expression programs of a cell even, under the most subtle perturbations. Such differentially expressed (DE) gene lists are often quite extended, including a great number of genes, for which there is little if any knowledge related to the system under study. In this regard, considerable effort has been directed towards methods for the efficient analysis and interpretation of whole transcriptome read-outs [4–6]. Such analyses focus mostly on the testing of DE genes for enrichment in various functional groupings, such as Gene Ontology (GO) terms , or molecular pathways such as those compiled by the Kyoto Encyclopedia of Genes and Genomes (KEGG) . Biologists have thus to choose from a variety of existing tools for data analysis and interpretation.
Over the years, the accumulation of genome-wide data has increased the possible gene groupings and categorizations, for which enrichment analyses may be conducted. These may now include protein families, molecular signatures defined under certain pathological conditions, chromosomal territories or co-expressed gene clusters obtained through meta-analyses of publicly available datasets [9, 10]. Among the various available gene categorizations, those referring to gene regulation are of particular interest, not only because of their wealth, encompassing predicted and experimental transcription factor binding sites and miRNA targets, but mainly because of their potential in inferring the gene regulatory program responsible for the observed expression profile.
Gene regulation takes place in various stages among which, transcription is the one most readily analyzed and easy to quantify. Given a certain stimulus or under specific conditions, the relative abundance of a great number of mRNA species may vary due to both orchestrated changes resulting from the activation of a particular gene expression program and random noise. The main goal of a functional analysis at the regulatory level will thus be to distinguish between the two and, moreover, to propose a hierarchy for the gene regulators involved in the system under examination. The concept of hierarchical regulatory interactions between genes is not new. Master regulators are important drivers of gene expression  and defining them is of primary interest at both experimental  and theoretical levels [13, 14].
In spite of the increase of available information, the problem of determining the hierarchy of transcriptional regulators involved in given conditions remains an open question. The definition of regulatory networks of interaction is a complex difficult task that may only be achieved through the integration of multiple datasets from various sources (TF binding, miRNA expression, gene expression etc.). Currently, the reconstruction of such global networks has been limited in the context of large genomic consortia (e.g. the ENCODE Project Consortium) [15, 16], or small unicellular organisms [17, 18], but even in these cases the resulting networks are extensive and difficult to interpret. At the same time, gene expression profiles are rapidly accumulating, exploring a vast amount of possible regulatory patterns and pressing for more efficient analysis. In this sense, it becomes plausible to seek ways to predict the regulatory network using only gene expression data, in attempts to treat whole genome expression profiles as a detailed reflection of the underlying regulatory program.
Accumulating genome-wide data, coupled with detailed studies has led to the creation of large compendia of well-defined regulatory interactions for a number of model species, compendia that are now being compiled in specialized databases. The use, however, of the reported resources requires filtering of noisy or trivial information. HTRIdb  contains a large number of interactions (>50000), the largest part of which are inferred from ChIP experiments that are known to be extremely noisy. ORegAnno  contains a more moderate number of regulatory interactions but in many cases these are reported as based on “unknown evidence” or refer to unknown genes. Smaller databases such as TRED , or TFactS , on the other hand, are built through a more thorough process that involves manual curation of literature and public datasets. Finally, there are databases employing intermediate approaches such as TRRUST , which makes use of a text-mining algorithm coupled with manual curation of the results to populate a database of ~8000 interactions.
We use the proposed methodology to infer regulatory networks for two test cases, one for each organism for which we have compiled information (human, mouse). We show RNEA able to reconstruct networks that are in agreement with existing knowledge of the systems under study, while at the same time they provide lists of additional candidate genes being involved in key processes. Our implementation shows RNEA to be a very useful resource for a first-level analysis of gene expression datasets in order to gain insight in the system of study and to obtain leads for genes and proteins of primary importance.
The Regulatory Network Enrichment Analysis tool (RNEA) is based on a collection of regulatory interactions compiled from manually curated databases. RNEA uses prior knowledge, coupled with standard statistical methods for the inference of active regulators, miRNAs and functional categories. Most importantly, RNEA extracts the presumably active regulatory subnetwork from a global gene regulatory network (GRN) based on the calculated transcriptional regulator - gene target enrichments, showing how significantly enriched regulators interact with their target genes and between each other. RNEA receives a gene expression profile, in the form of a complete list of differential expression values as input. This should include gene name identifiers coupled with differential expression values and significance p-values. Based on this list of differential expression, fold-change values and their accompanying p-values, it calculates enrichments for particular gene groupings. The output of RNEA is dual; on one hand it produces lists of over-represented gene categories in the form of current state-of-the-art approaches, but most importantly, it provides the user with a regulatory subnetwork file where the relevant gene interactions are registered. The lists can help identify important regulators and functions, while the created regulatory subnetwork, provides a view of the transcriptional regulation at a system level and may enhance the interpretation of a genome-wide expression experiment.
RNEA has been developed in R, aiming to be a cross-platform and easy-to-use tool. It is compatible with widely used differential expression analysis software such as Cufflinks , EdgeR  and DESeq  and may therefore be easily incorporated in already existing pipelines. Its results are displayed as HTML files with sortable tables, which include the corrected p-values for the functional and the regulatory groups respectively for greater ease of use. The regulatory subnetwork is extracted in a tab separated file format in order to be compatible with typical network visualization software, such as Cytoscape . The source code, alongside the reference networks for human and mouse and detailed documentation may be found at: https://sites.google.com/a/fleming.gr/rnea/.
In order to create a highly-confident dataset of transcriptional regulator–gene target interactions we searched for databases, which mainly contain experimentally validated or manually curated interactions. With these criteria, TRED , TFactS , Oreganno  and TRRUST  were chosen and lists of human and mouse regulator–gene-target interactions were created separately. TRED is a database designed as a resource for gene regulation. It has gathered data for many elements of regulation, such as promoters, many of which are annotated with computational tools produced by the same group. Their analysis is done genome-wide for human, mouse and rat. In this way, TRED combines promoter annotation with experimental results to assign target genes to transcriptional regulators, assignments, further refined through manual curation of the results and validation from the literature. TRED was organized in modular gene regulatory networks (GRNs) that were created and uploaded in the form of network figures (https://cb.utdallas.edu/TRED/GRN/grn.htm). These GRNs include visual representation of 34 transcription factor (TF) families. All images for human and mouse were downloaded and, in a thorough and time-consuming procedure, each interaction was recorded in a tab-separated text file.
Most of the interactions TFactS includes, overlap the dataset compiled by TRED. Most of the non-overlapping interactions are based on manual curation of literature articles, missing from the TRED reference database. We only kept species-specific interactions because in spite of the extended conservation between human and mouse at gene level, one cannot rule out significant differences existing at the level of protein-protein and protein-DNA interaction hierarchies between the two species. This was observed in the case of TRED, where in many cases the GRNs in human and mouse significantly differed for the same TF family.
Interactions between miRNA and protein coding genes were also compiled in a way that focused on experimentally validated datasets. Data were retrieved from Tarbase , currently the biggest repository containing experimentally validated miRNA-gene interactions. Tarbase also includes various types of information for each interaction, and computational predictions of microT , a computational approach for miRNA target prediction that is based on a combination of experimental data and sequence conservation. The extraction of interactions from Tarbase was performed for human and mouse in a way that excluded all interactions which were not experimentally validated, which resulted in two lists containing 1573 and 407 interactions respectively.
RNEA also provides standard enrichment tables for Gene Ontology and Pathway annotations, the latest versions of which were downloaded from the corresponding web resources. Data for GO were obtained from the Gene Ontology Consortium (http://geneontology.org/) and Biological Pathway annotations were retrieved from the Kyoto Encyclopedia of Genes and Genomes (KEGG) (http://www.genome.jp/kegg/).
RNEA performs analysis at two levels. After performing a typical enrichment analysis at the levels of GO, KEGG Pathways, miRNA and protein transcriptional regulators (TF), it couples the latter with a search in the corresponding reference regulatory network, in order to extract a subset of interactions and to reconstruct the most informative regulatory subnetwork.
At a first level, RNEA employs a standard over-representation analysis in order to calculate the enrichment of deregulated genes among certain categories. In this way, it requires defined fold-change and p-value thresholds, which may be provided by the user. As this is a typical example of a statistical experiment of drawing an object from a finite population which has two distinct states without replacement, the hypergeometric test is used to calculate the significance of the enrichment. RNEA differs from most tools in the sense that it automatically performs three different types of enrichments, aiming to find regulatory (TF, miRNA) or functional (KEGG pathways, GO terms) components whose members are over-represented among a) overall DE genes, b) strictly over-expressed or c) strictly under-expressed genes. It thus also processes information separately for genes whose expression is increased or decreased. Reporting is performed through ranked lists of enrichment, after conducting a suitable FDR correction for multiple testing. The output is organized in a set of prioritized lists of regulators and functional groups.
A profile is created for each transcriptional regulator containing all of its target genes. These first-level regulators are called “parent” regulators.
If any of the target genes is also a regulator, it is linked to both its “parent” regulator and its “daughter” targets, thus creating a series of second-level interactions.
This two level profile is then used in order to extract nodes and interactions from the reference network according to three simple rules:
o A regulator is included in the subnetwork if it is differentially expressed.
o First-level targets of the regulator are included in the subnetwork if they are differentially expressed.
o First-level targets of the regulator are included in the subnetwork even if they are NOT differentially expressed as long as their “parent” regulator and a second-level “daughter” target are both differentially expressed
The size of the profile and the number of deregulated genes in the profile are used for the statistical calculation of the enrichment as described above.
Our final goal is to capture the hierarchical structure of gene expression regulation and impart more depth in the regulatory network extracted. We chose to constrain our analysis to two hierarchical levels for two reasons: One is that a significant proportion of interactions (~70 % for human) are 1st-level interactions that lead directly from a regulator to a leaf end node in the network. It is thus reasonable to expect that a two-level approach will encompass the overwhelming majority of interactions. The second reason is that, by attempting to track the small proportion of higher-level (>2) interactions, we would inevitably incorporate a great number of cyclic-interactions in the parsing of the network. This would bring about a significant slowing down of the whole process. Restricting our analysis to two levels we thus achieve the most efficient ratio of retrieved information over processing time. RNEA progressively builds a flatfile, which contains a subset of the initial reference “super-network”. This may be directly visualized through open-source network visualization and analysis platforms. Images presented in this paper were produced with Cytoscape . One of the advantages of RNEA compared to other similar approaches is that it may incorporate protein and miRNA regulators in the same network. This is achieved via the incorporation of information from both transcriptional regulator and miRNA interactions.
RNEA is implemented as an R package which may be downloaded from https://sites.google.com/a/fleming.gr/rnea/home. It only requires one additional R package called “SortableHTMLTables” which is used in order to report the results in HTML format. RNEA currently accepts HUGO and RefSeq gene names as gene identifiers. It is also advised that users report differential expression as log2(fold-change) values in accordance with standard software. Nevertheless, the user may analyze data sets with non-standard differential expression values as long as he suitably adjusts the corresponding parameters. Both fold-change and p-value thresholds are set by the user, since quite often criteria need to be relaxed or made stricter in order to result in a reasonable number of differentially expressed genes that is sufficient for a statistical functional analysis.
A “global” network that includes TF, miRNA and functional categories connected with their respective gene members that are enriched in differentially expressed genes. The goal of this approach is to extract central regulatory or functional components with the use of Network properties. The combined regulatory/functional network is extracted with the use of the network = “global” argument.
Only functional or only regulatory networks can be extracted with the use of the respective setting “functional”, or “regulatory” of the network argument. Given that RNEA’s main aim is producing concise and summarized regulatory networks, “regulatory” is the default argument.
Other parameters include the species from which expression values have been obtained (“Human”, “Mouse”) and the output type (“html” or “csv”).
More details may be found here: https://sites.google.com/a/fleming.gr/rnea/manual.
Inferring important regulators from genome-wide expression experiments is a complex problem. There are only a couple of available tools which may help in the definition of master regulators mostly by finding over-represented binding motifs in deregulated genes. TFactS is the only tool using a similar enrichment approach to ours, but it has significant limitations in terms of the number of studied regulators. TRRUST  represents a recent attempt to provide a golden standard, against which regulatory networks may be tested. Nevertheless, the variability of expression programs is immense and the underlying complexity of gene regulation suggests that very different networks may be produced with even mild changes in cellular conditions. In this sense, already available networks can only serve as providing the “reference” interactions, among which each condition may choose combinatorially. RNEA’s aim is providing a framework for revealing such combinations of known interactions. Its distinguishing features are that a) it focuses on well-documented regulatory interactions b) it aims at capturing the hierarchical structure of the network by the two-level scanning of regulator profiles and c) that it incorporates miRNA and protein regulators in a common regulatory network (see Methods). From the application’s point of view, RNEA is rather straightforward and in principle can be run with only a small list of prerequisites.
RNEA is able to infer, in a single run, both regulatory and functional enrichments from raw differential expression data. To date there are only a couple of similar methodologies whose scope however differs from RNEA. SPIA  implements a perturbation analysis to infer the significance of a given pathway based on the differential expression of its genes, while taking into account the topology of the pathway network. In this sense, it aims at a better and more accurate assessment of pathway deregulation based on the reported interactions and its final output is a list of deregulated pathways. PARADIGM  employs a similar approach through the additional incorporation of multiple omics and genetic data. Finally GGEA  is similar to our approach in terms that it combines gene expression and regulatory interactions aiming at an initial refinement of differentially expressed genes list, which it then uses for a gene set enrichment analysis. Compared to the above, our method’s distinguishing characteristic is that it projects DE gene lists on a reference map of regulatory interactions to infer a subnetwork of relevance to the particular gene expression profile. Providing the subnetwork reconstruction as primary output is thus RNEA’s particular feature, but the lack of similar approaches makes its cross-validation quite difficult. In order to assess its predictive power we have here applied it in cases of well-defined systems, for which there is documented knowledge of the underlying regulatory program. This is an approach also undertaken by the aforementioned methods [27, 29]. In this context, we chose to validate our methodology in two publicly available datasets originating from the two species for which RNEA provides information, a genome wide expression profile of cancerous versus normal human colon tissue samples and a gene profile of mouse RAW264.7 macrophage cells, stimulated with LPS.
In the following section, we briefly present the results of the analysis of the 2 test cases conducted with RNEA. As our method focuses primarily on transcriptional regulation we were more interested in assessing the robustness of the regulatory networks. In this sense, we attempted to validate the primary nodes of the deduced networks from the existing literature.
Human test case
Colorectal cancer is the third most common type of cancer and the second most common cause of cancer death in the human population . Despite a considerable amount of evidence delineating biological pathways related to the disease, the characterization of regulatory networks in cancer remains an open problem due to the great number of disease subtypes and the overall variability of the phenotypes.
In an attempt to predict a regulatory subnetwork for colorectal cancer, we obtained genome-wide expression data from 104 patients and 46 healthy individuals, that were normalized with robust multi-array average (RMA) and presented in a form suitable for RNEA analysis (log2(FC) of gene expression alongside the corresponding p-values) from  (GEO Accession: GSE21510). Patients with distant metastases were selected in order to assess the metastatic potential conferred to the disease by candidate gene markers. Using standard fold-change and p-value thresholds we obtained 148 DE genes. The functional analysis conducted with RNEA resulted in a number of deregulated GO terms and KEGG pathways that may be accessed in Additional file 1. Below we focus on the regulatory analysis.
At the periphery of our network, apart from a set of rather expected regulators whose relation to cancer is well known (FOS, JUN, EGR1 and MYC) we find a set of secondary metabolism related genes including HSD17B2, a gene that is involved in lipid biosynthesis that has been shown to have a prognostic role in colorectal cancer . At the same time all of the reported miRNA in the network have been reported to have a role in colorectal cancers. In particular, hsa-miR-26b-5p has been shown to exhibit a tumor suppressive role , hsa-miR-590-3p has been also known to have increased expression in colorectal cancers  while hsa-miR-374a-5p has been reported to be beneficial for the prognosis . An important aspect of RNEA may be seen here. RNEA by default performs a dual enrichment analysis separately for over- and under-expressed gene targets. As expected, all miRNA-genes with significant enrichments in this setting are enriched primarily towards under-expressed targets. We consider this feature of great importance since a number of existing methods merely report enrichment on the basis of differential expression, regardless of the direction of regulation (activation or suppression). In certain cases however (miRNA regulation being a very clear one) this direction should be taken into account.
Mouse test case
We next analyzed a test case for the mouse genome. We chose a well-described process of an external stimulus that is expected to show an inflammatory response, which we indeed observe through our functional analyses. We examined an experiment that aimed to analyze the differences in gene expression of the inflammatory response in RAW264.7 murine macrophage cells under stimulation by LPS  (GEO Accession: GSE63889). Using relatively stringent criteria on differential expression fold-change and p-value (absolute log2FC > =1, p-value < =0.05) we ended up with 121 differentially expressed genes, which were shifted towards over-expression with a ratio close to 2:1. This is to be expected given that LPS stimulation of immune cells is known to bring about an acute generalized response through the activation of a number of pathways. Functional analysis at the level of GO and KEGG pathways showed these pathways to match the expectations. In this way, inflammatory pathways including the TNF-NFkB signaling axis, cytokine and TLR signaling and a range of various infection-responsive pathways were strongly enriched in upregulated genes (see Additional file 2 for the complete lists).
Reconstructing the hierarchy of gene regulation is an open question of great importance in modern molecular biology, one that has not been effectively answered by the advent of next-generation sequencing technologies and genome-wide interrogation of gene expression. The main reasons for this, besides the inherent complexity of this dynamic process have been a) the scarcity of high-quality data (or -in contrast- the considerable amount of noisy information) connecting transcription factors and their gene targets and b) the lack of user-friendly methods summarizing the output of genome-wide expression experiments in a way that will allow for a fast and meaningful, first-level inspection of the data. The main goal of the presented work lies on these two axes: On one hand it aims to produce a compendium of transcription factor–target genes interactions that would be as reliable as possible, on the other, to implement these relationships, alongside other meaningful, functional information in a simple approach that would allow biologists to perform a quick overview of a gene expression experiment, prioritizing their results and putting the spotlight on highly significant regulatory interactions.
A major limitation of genome-scale experiments has, since very early, been a lack of summarization in their analysis. The production of enormous lists of genes, enriched processes and functions often results in complicating, rather than assisting in, their interpretation. Key aspects in the process of extracting knowledge from a large-scale experiment are related to a) efficient summarization, i.e. compiling the useful and significant information from a functional point of view and b) prioritization, the ranking, that is, of the relevant pieces of information in a way that will help focus on the most important facets of the results. RNEA addresses the problem of extensive lists by producing regulatory networks that combine prioritization and summarization of the observed enriched relationships. In this way, it reveals the genes that are important in the process of transcriptional regulation in the particular experimental setting. This may be done either through the assessment of the statistical enrichment of TF targets’ profiles or by evaluating network characteristics of the extracted regulatory subnetwork. Our work’s distinguishing feature is the output of an “active regulatory subnetwork”, which constitutes the most probable network of transcription factors being active in the studied condition, based on the gene expression values of their target genes. This also represents the major advantage of our pipeline, whose main goals are to detect and report regulatory information in the shape of regulatory networks that are of modest sizes, therefore allowing the user to easily interpret results and plan follow-up experiments.
When compared to other similar approaches (such as TFactS or Enrichr) our method’s main differences are related with the way it infers the regulatory hierarchy. RNEA does this by employing a two-level hierarchical approach, by adding at each TF profile the targets of its targets, if available. In this way TFs which are “higher” in the known hierarchy will have more targets in their profile and if an enrichment is found, it will add significance to the TF lying higher in the hierarchy. In addition, by checking for enrichment in both up-regulated and down-regulated genes in two different calculations, RNEA enables the investigation of possible dual roles for given TF, while at the same time to safeguard against contradictory annotation evidence that often assign inconsistent roles for a TF based on the literature.
Currently RNEA only supports human and mouse datasets. Given their relevance from the biomedical perspective these two organisms (human for obvious reasons and mouse due to the fact of being the most widely-used mammalian model organism) represent more than 90 % of the public repositories of gene interaction data. Their share of mammalian genome-wide expression profiles in gene expression databases is probably even higher. It was therefore reasonable to aim at the construction of reference networks for these two at a first level. Nevertheless, the incorporation of information on other organisms within the RNEA framework is expected to be quite straight-forward, given that a sufficient number of experimentally verified regulatory interactions are reported.
As a final comment, we should point out that RNEA deduces gene regulatory interactions directly from their noisy, highly complex end-product which is the relative abundance of mRNA molecules. In this regard, it is expected that the predicted networks will also be partly noisy and contain a reasonable amount of false positives. We, nonetheless believe, that being able to visualize a relevant network of interactions in a single step from your differential expression experiment makes up for a positive trade-off.
RNEA is a framework for functional analysis of gene expression experiments, with a primary focus on gene regulatory relationships. It is easy to apply on standard gene expression read-outs, readily producing ranked lists of various functional groupings. Its key idea, though, is the derivation of a network of regulatory interactions. By creating regulatory subnetworks, RNEA enables a better overview of the regulatory process through direct visualization. RNEA benefits from (and also depends on) the accuracy of the prior knowledge used and the originality of network reconstruction. Most of the existing functional analysis tools mainly rely on computational predictions (through PWM) for the calculation of transcription factor target enrichments, while very few also employ experimental data from ChIP that are, however, still limited.
The main advantage of RNEA is the originality of the network approach. To our knowledge this is one of the few functional analysis tools that aims at the extraction of a regulatory subnetwork. Most of the existing approaches in this regard have been implemented in a sort of ‘personalized’ way, aiming at the interpretation of specific experiments, instead of proposing a generalized approach. Further validation of RNEA can mostly take place with extended use by the community, while refinements in the original reference networks are bound to increase its potential.
Availability and requirements
Project name: Regulatory Network Enrichment Analysis (RNEA)
Project home page: https://sites.google.com/a/fleming.gr/rnea/home
Operating system(s): Multiple platforms
Programming language: R
Other requirements: R package “SortableHTMLTables” (http://cran.r-project.org/web/packages/SortableHTMLTables/index.html)
License: GNU GPLv3
Any restrictions to use by non-academics: None
Τhis work has been supported by IMI-funded project BeTheCure (BTCure, 115142–2), GSRT project INNATE FIBLROBLAST (ERC06, co-financed by the ESF and NSRF 2007–2013) and Advanced ERC grant MCs-inTEST (340217) to GK. The authors would like to thank the organizers of the SMODIA 2014 (Statistical Methods for Omics Data Integration and Analysis) Workshop for giving them the opportunity to present this work to a broader audience. The authors would like to thank three anonymous reviewers for their fruitful comments and recommendations for additional resources to be used in the formation of the reference networks.
Publication charges for this article have been funded by IMI project BeTheCure (BTCure, 115142-2).This article has been published as part of BMC Bioinformatics Volume 17 Supplement 5, 2016: Selected articles from Statistical Methods for Omics Data Integration and Analysis 2014. The full contents of the supplement are available online at http://bmcbioinformatics.biomedcentral.com/articles/supplements/volume-17-supplement-5.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Clark MB, Amaral PP, Schlesinger FJ, Dinger ME, Taft RJ, Rinn JL, et al. The reality of pervasive transcription. PLoS Biol. 2011;9(7):e1000625. discussion e1001102.View ArticlePubMedPubMed CentralGoogle Scholar
- Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28(5):511–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Rapaport F, Khanin R, Liang Y, Pirun M, Krek A, Zumbo P, et al. Comprehensive evaluation of differential gene expression analysis methods for RNA-seq data. Genome Biol. 2013;14(9):R95.View ArticlePubMedPubMed CentralGoogle Scholar
- Al-Shahrour F, Minguez P, Vaquerizas JM, Conde L, Dopazo J. BABELOMICS: a suite of web tools for functional annotation and analysis of groups of genes in high-throughput experiments. Nucleic Acids Res. 2005;33(Web Server issue):W460–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Reimand J, Kull M, Peterson H, Hansen J, Vilo J. G:Profiler-a web-based toolset for functional profiling of gene lists from large-scale experiments. Nucleic Acids Res. 2007;35(Web Server issue):W193–200.View ArticlePubMedPubMed CentralGoogle Scholar
- Huang DW, Sherman BT, Lempicki R. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37(1):1–13.View ArticleGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000;25(1):25–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Ogata H, Goto S, Sato K, Fujibuchi W, Bono H, Kanehisa M. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 1999;27:29–34.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen EY, Tan CM, Kou Y, Duan Q, Wang Z, Meirelles GV, et al. Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinformatics. 2013;14:128.View ArticlePubMedPubMed CentralGoogle Scholar
- Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL. Gene set enrichment analysis : A knowledge-based approach for interpreting genome-wide. Proc Natl Acad Sci U S A. 2005;102(43):15545–50.View ArticlePubMedPubMed CentralGoogle Scholar
- Ohno S. Major sex-determining genes. Am J Hum Genet. 1980;32(3):464.Google Scholar
- Milo R, Shen-Orr S, Itzkovitz S, Kashtan N, Chklovskii D, Alon U. Network motifs: simple building blocks of complex networks. Science. 2002;298(5594):824–7.View ArticlePubMedGoogle Scholar
- Kashtan N, Alon U. Spontaneous evolution of modularity and network motifs. Proc Natl Acad Sci U S A. 2005;102(39):13773–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Bovolenta LA, Acencio ML, Lemke N. HTRIdb: an open-access database for experimentally verified human transcriptional regulation interactions. BMC Genomics. 2012;13:405.View ArticlePubMedPubMed CentralGoogle Scholar
- Essaghir A, Toffalini F, Knoops L, Kallin A, van Helden J, Demoulin J-B. Transcription factor regulation can be accurately predicted from the presence of target gene signatures in microarray gene expression data. Nucleic Acids Res. 2010;38:e120.View ArticlePubMedPubMed CentralGoogle Scholar
- Han H, Shim H, Shin D, Shim JE, Ko Y, Shin J, et al. TRRUST: a reference database of human transcriptional regulatory interactions. Sci Rep. 2015;5:11432.View ArticlePubMedPubMed CentralGoogle Scholar
- Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.View ArticleGoogle Scholar
- Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.View ArticlePubMedGoogle Scholar
- Griffith OL, Montgomery SB, Bernier B, Chu B, Kasaian K, Aerts S, et al. ORegAnno: an open-access community-driven resource for regulatory annotation. Nucleic Acids Res. 2007;36(Database issue):D107–13.View ArticlePubMedPubMed CentralGoogle Scholar
- Jiang C, Xuan Z, Zhao F, Zhang MQ. TRED: a transcriptional regulatory element database, new entries and other development. Nucleic Acids Res. 2007;35(Database issue):D137–40.View ArticlePubMedPubMed CentralGoogle Scholar
- Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11(10):R106.View ArticlePubMedPubMed CentralGoogle Scholar
- Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.View ArticlePubMedPubMed CentralGoogle Scholar
- Barabási A-L, Albert R. Emergence of scaling in random networks. Science. 1999;286(5439):509–12.View ArticlePubMedGoogle Scholar
- Barabási A-L, Oltvai ZN. Network biology: understanding the cell’s functional organization. Nat Rev Genet. 2004;5(2):101–13.View ArticlePubMedGoogle Scholar
- Vlachos IS, Paraskevopoulou MD, Karagkouni D, Georgakilas G, Vergoulis T, Kanellos I, et al. DIANA-TarBase v7.0: indexing more than half a million experimentally supported miRNA:mRNA interactions. Nucleic Acids Res. 2014;43(Database issue):D153–9.PubMedPubMed CentralGoogle Scholar
- Maragkakis M, Reczko M, Simossis VA, Alexiou P, Papadopoulos GL, Dalamagas T, et al. DIANA-microT web server: elucidating microRNA functions through target prediction. Nucleic Acids Res. 2009;37(Web Server issue):W273–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Tarca AL, Draghici S, Khatri P, Hassan SS, Mittal P, Kim JS, et al. A novel signaling pathway impact analysis. Bioinformatics. 2009;25(1):75–82.View ArticlePubMedGoogle Scholar
- Vaske CJ, Benz SC, Sanborn JZ, Earl D, Szeto C, Zhu J, et al. Inference of patient-specific pathway activities from multi-dimensional cancer genomics data using PARADIGM. Bioinformatics. 2010;26(12):i237–45.View ArticlePubMedPubMed CentralGoogle Scholar
- Geistlinger L, Csaba G, Küffner R, Mulder N, Zimmer R. From sets to graphs: towards a realistic enrichment analysis of transcriptomic systems. Bioinformatics. 2011;27(13):i366–73.View ArticlePubMedPubMed CentralGoogle Scholar
- Ricchi P, Zarrilli R, Di Palma A, Acquaviva AM. Nonsteroidal anti-inflammatory drugs in colorectal cancer: from prevention to therapy. Br J Cancer. 2003;88(6):803–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Tsukamoto S, Ishikawa T, Iida S, Ishiguro M, Mogushi K, Mizushima H, et al. Clinical significance of osteoprotegerin expression in human colorectal cancer. Clin Cancer Res. 2011;17(8):2444–50.View ArticlePubMedGoogle Scholar
- Gordziel C, Bratsch J, Moriggl R, Knösel T, Friedrich K. Both STAT1 and STAT3 are favourable prognostic determinants in colorectal carcinoma. Br J Cancer. 2013;109(1):138–46.View ArticlePubMedPubMed CentralGoogle Scholar
- Iacopetta B. TP53 mutation in colorectal cancer. Hum Mutat. 2003;21(3):271–6.View ArticlePubMedGoogle Scholar
- Hu R, Zuo Y, Zuo L, Liu C, Zhang S, Wu Q, et al. KLF4 Expression Correlates with the Degree of Differentiation in Colorectal Cancer. Gut Liver. 2011;5(2):154–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Lee H-Y, Ahn JB, Rha SY, Chung HC, Park KH, Kim TS, et al. High KLF4 level in normal tissue predicts poor survival in colorectal cancer patients. World J Surg Oncol. 2014;12:232.View ArticlePubMedPubMed CentralGoogle Scholar
- Lee Y-E, He H-L, Shiue Y-L, Lee S-W, Lin L-C, Wu T-F, et al. The prognostic impact of lipid biosynthesis-associated markers, HSD17B2 and HMGCS2, in rectal cancer treated with neoadjuvant concurrent chemoradiotherapy. Tumour Biol. 2015;36(10):7675–83.View ArticlePubMedGoogle Scholar
- Zeitels LR, Acharya A, Shi G, Chivukula D, Chivukula RR, Anandam JL, et al. Tumor suppression by miR-26 overrides potential oncogenic activity in intestinal tumorigenesis. Genes Dev. 2014;28(23):2585–90.View ArticlePubMedPubMed CentralGoogle Scholar
- Cummins JM, He Y, Leary RJ, Pagliarini R, Diaz LA, Sjoblom T, et al. The colorectal microRNAome. Proc Natl Acad Sci U S A. 2006;103(10):3687–92.View ArticlePubMedPubMed CentralGoogle Scholar
- Slattery ML, Herrick JS, Mullany LE, Valeri N, Stevens J, Caan BJ, et al. An evaluation and replication of miRNAs with disease stage and colorectal cancer-specific mortality. Int J Cancer. 2015;137(2):428–38.View ArticlePubMedGoogle Scholar
- Sabari BR, Tang Z, Huang H, Yong-Gonzalez V, Molina H, Kong HE, et al. Intracellular Crotonyl-CoA Stimulates Transcription through p300-Catalyzed Histone Crotonylation. Mol Cell. 2015;58(2):203–15.View ArticlePubMedPubMed CentralGoogle Scholar
- Shi L, Kishore R, McMullen MR, Nagy LE. Lipopolysaccharide stimulation of ERK1/2 increases TNF-alpha production via Egr-1. Am J Physiol Cell Physiol. 2002;282(6):C1205–11.View ArticlePubMedGoogle Scholar
- Gomard T, Michaud H-A, Tempé D, Thiolon K, Pelegrin M, Piechaczyk M. An NF-kappaB-dependent role for JunB in the induction of proinflammatory cytokines in LPS-activated bone marrow-derived dendritic cells. PLoS ONE. 2010;5(3):e9585.View ArticlePubMedPubMed CentralGoogle Scholar
- Eissmann P, Evans JH, Mehrabi M, Rose EL, Nedvetzki S, Davis DM. Multiple mechanisms downstream of TLR-4 stimulation allow expression of NKG2D ligands to facilitate macrophage/NK cell crosstalk. J Immunol. 2010;184(12):6901–9.View ArticlePubMedGoogle Scholar
- Bazzoni F, Rossato M, Fabbri M, Gaudiosi D, Mirolo M, Mori L, et al. Induction and regulatory function of miR-9 in human monocytes and neutrophils exposed to proinflammatory signals. Proc Natl Acad Sci U S A. 2009;106:5282–7.View ArticlePubMedPubMed CentralGoogle Scholar