- Research article
- Open Access
A systems biology model of the regulatory network in Populusleaves reveals interacting regulators and conserved regulation
© Street et al; licensee BioMed Central Ltd. 2011
- Received: 8 October 2010
- Accepted: 13 January 2011
- Published: 13 January 2011
Green plant leaves have always fascinated biologists as hosts for photosynthesis and providers of basic energy to many food webs. Today, comprehensive databases of gene expression data enable us to apply increasingly more advanced computational methods for reverse-engineering the regulatory network of leaves, and to begin to understand the gene interactions underlying complex emergent properties related to stress-response and development. These new systems biology methods are now also being applied to organisms such as Populus, a woody perennial tree, in order to understand the specific characteristics of these species.
We present a systems biology model of the regulatory network of Populus leaves. The network is reverse-engineered from promoter information and expression profiles of leaf-specific genes measured over a large set of conditions related to stress and developmental. The network model incorporates interactions between regulators, such as synergistic and competitive relationships, by evaluating increasingly more complex regulatory mechanisms, and is therefore able to identify new regulators of leaf development not found by traditional genomics methods based on pair-wise expression similarity. The approach is shown to explain available gene function information and to provide robust prediction of expression levels in new data. We also use the predictive capability of the model to identify condition-specific regulation as well as conserved regulation between Populus and Arabidopsis.
We outline a computationally inferred model of the regulatory network of Populus leaves, and show how treating genes as interacting, rather than individual, entities identifies new regulators compared to traditional genomics analysis. Although systems biology models should be used with care considering the complexity of regulatory programs and the limitations of current genomics data, methods describing interactions can provide hypotheses about the underlying cause of emergent properties and are needed if we are to identify target genes other than those constituting the "low hanging fruit" of genomic analysis.
- Drought Stress
- Transcriptional Module
- Nucleosome Assembly
- System Biology Model
- Biotic Infection
Biologists have long been fascinated by the green plant leaf and have tried to understand how leaves are born, live and die. In the last decades, several new approaches to study the structure and function of leaves have emerged: Molecular biology and molecular genetics have, for example, enabled identification of genes that regulate the primary function of the leaf - photosynthesis - and leaf development has been understood in much greater detail; high through-put transcriptomics has identified additional factors influencing leaf function, but traditional transcriptome analyses typically reduces the problem of finding key regulators to detecting differentially expressed genes or computing pair-wise similarity between targets and putative regulators (e.g. hierarchical clustering or co-expression networks). In contrast, systems biology analysis of transcriptional programs treats genes as interacting rather than isolated entities. Thus these methods can begin to understand how so-called emergent properties such as complex phenotypes arise from interacting genes. Whether this can be seen as taking a holistic rather than a reductionistic approach to science has generated quite some debate [1, 2], but systems biology methods account for synergistic and competitive effects between regulators that individually could have low similarity to the target. Methods for reverseengineering the transcriptional network from collections of gene expression data have been pioneered on single-cell organisms, but have increasingly been applied to higher order organisms  including plants [4, 5] where applications of systems biology methods are now emerging. Most systems biology studies have - not surprisingly - utilized using "THE model plant" Arabidopsis thaliana, where large transcriptomics programs have generated adequate quantities of high-quality data to enable systems analysis . For example, Carerra et al.  modeled the transcriptional network of Arabidopsis and identified plant-specific properties such as high connectivity between genes involved in response and adaptation to changing environments. However, not all aspects of plant biology can be studied in Arabidopsis, which in many respects is a rather atypical plant. Indeed, it was not selected as a model system due to its physiological and ecological qualities, but rather for its suitability for genetic and genomic studies. Therefore, it is important to perform parallel studies in plants with other characteristics, as well as developing the methods to allow data from the Arabidopsis system to inform studies in other organisms.
One rapidly emerging plant model system is Populus ; it's interesting biology (a woody perennial) and the access to a sequenced genome  represent an attractive combination. Correspondingly, more advanced data analyses approaches are now being applied in Populus. Populus provides an attractive model system for studies of leaf biology. For example, Sjödin et al.  exploited the fact that mature aspen (Populus tremula) in boreal regions have the rather unique property that all leaves emerge simultaneously from overwintering buds. This provides a synchronized system, resulting in a full temporal separation of the leaf developmental stages and subsequent acclimation that could be exploited using transcriptomics. Access to a centralized repository of much of the Populus cDNA microarray data  and databases for the analysis of gene expression - and other - data  substantially facilitates the ability to perform systems biology studies. For example, Grönlund et al.  induced a co-expression network revealing modular architecture explaining gene function and tissue-specific expression; Street et al.  identified co-expression networks across a large collection of leaf transcriptomics data and found that some network hubs have existing functional evidence in Arabidopsis; Quesada et al.  performed a comparative analysis of the transcriptomes of Populus and Arabidopsis, and found evidence of extensive remodeling of the transcriptional network, although some essential functions showed little divergence. A few studies have also integrated promoter information to study regulatory control in Populus. Shi et al.  identified combinations of xylem-specific motifs in Populus promoters. Another study inferred transcriptional networks in xylem, leaves, and roots, and showed that genes with conserved regulation across tissues are primarily cis-regulated, while genes with tissue-specific regulation are often trans-regulated . All these studies are essentially co-expression networks that visualize expression similarity between pairs of genes, but do not infer complex interactions.
Network inference methods using expression data can be divided into those that aim to model the general influence that genes have on the expression of other genes (gene networks) [17, 18] and methods that aim to model the physical interaction between transcription factors and the regulated genes (gene regulatory networks) . Both approaches employ common network inference methods (see e.g. [20–22]), but those that infer gene regulatory networks also typically integrate motif finding and detection of transcriptional modules [23, 24]. Approaches that describe how the regulatory genome orchestrates dynamic gene expression has developed from Pilpel et al. , who showed that yeast genes sharing pairs of binding sites in their promoters were significantly more likely to be co-expressed than genes sharing only single binding sites, to various machine learning methods that identify modules of co-expressed genes with common motif patterns in their promoters (so-called cis-transcriptional modules)[26–34].
Here we apply a network inference method combining promoter information and expression data to describe the transcriptional network in Populus leaves. Our aims were (1) to detect regulatory hubs in leaves, (2) to describe conservation of transcriptional regulation within Populus and between Populus and Arabidopsis, and (3) to understand the regulatory complexity in leaves by comparing systems biology and traditional bioinformatics as methods for detecting target genes for further analysis. This study goes beyond previous meta-analyses of Populus transcriptome data by taking into account synergistic and competitive interactions between regulators, and by systematically integrating the regulatory genome and the transcriptome to infer networks. We show that our network is robust, explains available gene function information and generalizes to new expression data in both Populus and Arabidopsis. We identify the main regulators of primary processes in leaves, and show how some of these have regulatory partners orchestrating expression either in a synergistic or competitive manner. Such interactions are not considered by pair-wise similarity methods, and thus several of the regulators predicted here would not have been identified by traditional approaches.
Discovered transcriptional modules reflect important processes in leaves
One of the goals of this study was to investigate regulatory complexity. Interesting, very few of the discovered modules are associated with only one sequence motif (Figure 3B). Typically two or three motifs were required to find a significant correspondence between motifs and co-expression, indicating a complex relationship between observed expression and the regulatory genome. To evaluate the biological significance of the discovered modules, and their suggested regulatory control, we used functional annotations from Gene Ontology and KEGG. In general, 71% of the modules had some evidence of biological relevance in terms of over-represented Gene Ontology annotations (23 modules) and KEGG annotations (16 modules). Many of these were related to photosynthesis and ribosomal activity, and thus of relevance to leaf development (Figure 3C). Since all genes in this study were leaf-specific with a corresponding over-representation of leaf-specific annotations , one could argue that any division of these genes into modules would produce relevant annotations. However, in our statistical tests we used only the leaf-specific genes, not the whole genome, as background to avoid that typical leaf-functions show up as significant just because of the bias in the dataset. Hence, the large fraction of significant modules indicates that our division into modules based on common motifs and co-expression is indeed relevant. This was also confirmed by randomization experiments, which invariably resulted in modules with considerably lower significance than reported here.
Regulatory network indicates complex regulations
Predicted regulators of the Populus leaf transcriptional program.
subunit of chloroplast RNA polymerase, response to red and blue light
zinc ion binding
embryonic development ending in seed dormancy abscisic acid biosynthetic process, response to water deprivation, heat and osmotic
stress, xanthophylls biosynthetic process, sugar mediated signaling pathway, response to red light
histone H2A protein, nucleosome assembly
regulates cell growth, nuclear division and stem cell maintenance
basic helix-loop-helix (bHLH) family
histone H2A protein, nucleosome assembly
basic helix-loop-helix (bHLH) family regulation of flower development, meristem
structural organization, abaxial cell fate specification
response to auxin stimulus, lateral root morphogenesis
epidermal cell fate specification, seed coat development
ethylene mediated signaling pathway cinnamic acid biosynthetic process,
response to wounding, salt stress and abscisic and salicylic acid stimulus, negative regulation of metabolic process cell death, response to stress, ethylene
mediated signaling pathway, response to cytokinin stimulus, ethylene stimulus and other organism
involved in trichome and root hair patterning
Our method of increasingly evaluating more complex regulatory mechanism allowed us to quantify the complexity of the regulation in Populus leaves. The distribution of modules over the number of transcription factors in the predicted regulatory mechanism (Figure 5C) roughly follows that of the number of motifs (Figure 3B). Thus, the predictive power of the regulatory mechanisms of most modules benefit significantly from including more than one transcription factor. Both steps in our method predict expression of genes, however, while the module discovery approach finds sequence motifs predictive of gene expression clusters, the network inference approach finds transcription factors predictive of the gene expression in each module. Both approaches are guided by the principle of Occam's razor, that is, that the simplest model explaining the data is the best, and both approaches, as we have seen, result in the same distribution for the number of regulators per module.
The network is fully connected except for a small sub-network of the three nucleosome assembly modules discussed earlier. One of these modules is shown in Figure 2A, and is predicted to be regulated by 268609 (HTA7, closest homolog AT5G27670.1). This factor is a histone protein with a known role in nucleosome assembly (Table 1). The other two modules are predicted to be regulated by 268609 in concert with 232345 (HTA10, closest homolog AT1G51060.1), also a histone protein with a known role in nucleosome assembly. The protein 232345 is itself a member of the example module from Figure 2A. The fact that we did not allow auto-regulations in our inference method might thus be the reason why this module only has one regulator (i.e. 268609). The two modules associated with both factors are the two modules with the strongest competitive regulatory mechanisms in the network (Figure 6). Both these regulators have a significant individual influence on the expression of the modules, but they also have a highly significant negative cross-term indicating the competitive regulation. Intriguingly, these are the only two modules in the network with a significant co-expression during biotic infection, although they are also co-expressed in a number of other experiments.
Regulatory network predicts expression in unseen experiments
Bootstrap analysis is often used in computational studies to evaluate the statistical significance of models such as phylogenetic trees . A bootstrap dataset has the same number of genes and conditions as the original data, but with some conditions occurring several time and some conditions not occurring at all (i.e. drawn with replacement). On average, 36.8% of the conditions will not occur in the bootstrap dataset and we refer to this as the hold-out set. Our network was validated statistically by first inferring a number of networks from different bootstrap dataset, and then (a) assessing the agreement between these bootstrap networks and the original network (stability) and (b) using the regression models from the bootstrap networks to predict expression values in the hold-out sets (predictive power).
Our Populus network models show a remarkable ability to generalize to unseen conditions, although similar predictive capability has been demonstrated also for other organisms [4, 38]. Since we use the expression of a set of transcription factors to predict one expression profile per module, the correlation between observed and predicted expression is limited by the degree of expression similarity of genes within modules. Still, all co-expressed genes in modules had a significant correlation between observed and predicted expression when using the bootstrap networks to predict the expression in the hold-out sets (Figure 7B). In fact, 90% of genes, and all the modules, obtained a correlation above 0.5 (the original threshold for including genes in modules). We also held out entire experiments (e.g. budset, biotic infection, etc.) and used the resulting networks to predict the expression values in the missing experiment (Figure 7C). Since few modules have a significant expression similarity within modules in stress responses (Figure 3A), we are naturally unable to predict the expression in these experiments. However, the regulation of the developmental programs, in particular leaf primordia and budset, can be predicted from the other experiments (Figure 7C). This is also true for drought stress, indicating that regulation of drought response corresponds to the regulation of development in that there is a conserved relationship between regulating transcription factors and regulated gene modules. A notable exception is the nucleosome assembly modules from Figure 2A with a role in water deprivation response. This role is confirmed by the fact that the expression profile of this module cannot be predicted without the drought stress dataset (correlation -0.24 versus 0.56 in the bootstrap analysis).
Several regulatory mechanisms are conserved between Populus and Arabidopsis
Other studies have also investigated the conservation of gene expression across Populus and Arabidopsis. Quesada et al.  reported evidence of extensive evolution of gene expression regulation. Street et al.  identified hub-genes in leaf development, and quantified the fraction of conserved genes to about 60%. Our results seem to imply similar conclusions, although the present study directly identified conserved relationships between transcription factors and gene modules. An interesting question not addressed here is to what degree evolution of gene expression can be explained by divergence in the regulatory regions (promoter sequences) of the two species.
Systems biology predicts new leaf regulators
Our approach describes interactions between regulators by inferring sets of transcription factors that regulate modules in concert. This systems biology approach differs from traditional analysis such as hierarchical clustering or co-expression networks that only consider pair-wise similarity between the regulator and the regulated genes. To compare these two approaches, we also constructed a co-expression network where each module is regulated by the single transcription factor with the most similar expression to that module (Additional file 6). Table 1 lists all transcription factors in our systems biology-based network, and compares these to the regulators in this reductionistic co-expression network. While the co-expression network identifies 8 transcription factors as regulators of Populus leaf transcription, our network includes 20 of the 35 transcription factors in the data. From Figure 6 it was apparent that most collaborative regulations in our network have a master regulator, and this is the regulator typically identified by the co-expression network. Thus, most of the new regulators in our network are due to the fact that collaboration between transcription factors explains more of the expression in modules than single factors. Thus, although the regulations in our network are considerably stronger in terms of prediction power, they rarely exclude the transcription factor found in the co-expression network. However, two modules were predicted by the co-expression network to be regulated by transcription factor 639804 even though this most-similar factor was excluded as a regulator in our network. Somewhat surprisingly, the proposed regulatory mechanisms in these two cases are a weighted sum of two and three transcription factors without a statistically significant synergistic or competitive interaction (i.e. non-significant cross-terms).
One of the aims of systems biology is to model the complex interactions in living cells, describing emerging properties not apparent from studying genes, proteins or metabolites individually. Still, most computational approaches just take pair-wise similarity, not interactions between genes, into account when inferring network from expression data. The reason for this is at least two-fold. First, exploring combinatorics is computationally expensive. For example, there are over 2000 transcription factors in Populus giving rise to over 2 million pairs, 1.3 billion triplets, etc. Second, more complex models (e.g. cross-terms in regression models) imply many more parameters that have to be estimated from data (i.e. the β's in regression models). Since we need more observations than model-parameters to avoid over-fitting the models, the number of required observations grows quadraticly with the number of regulators when considering pairs. This curse of dimensionality represents a huge obstacle to studying interactions in biological systems. Here, we deal with these problems in several different ways. First, we restrict our study to leaf-specific genes rendering far fewer combinations than an unfiltered whole-genome study. Second, rather than considering all regulators at once, we devised a method that starts with single regulators, and then moves to pairs and higher-order combinations. This provides adequate observations to estimate parameters for each model (365 observations versus only four parameters in the case of two regulators), but because we test so many different models it comes with the risk of finding combinations that obtain high predictive power by change (i.e. over-fitting). We deal with this problem by only increasing model complexity if a statistically significant boost in predictive power is observed on unseen data (cross validation). In the statistical test we used the highly-conservative Bonferroni correction where the initial significance threshold (0.05) was divided by the number of transcription factor combinations tested.
For the leaf-specific genes studied here, the systems biology-based network mostly discovered co-regulators to the transcription factors also identified in the co-expression network. That means that although 11 of 38 modules had regulatory mechanisms with a significant interaction term (cross-term, Figure 6), these regulators also had significant individual contributions of which the strongest is detected by pair-wise similarity. A situation where the cross-term is significant, while the individual contributions are non-significant, is not observed in this data. An example of such a regulation is the logical XOR, that is, the regulated module is up-regulated only if one of the regulators is up-regulated (but not both). Whether such regulations exist in Populus leaves cannot be settled from this study considering the limited set of genes included. Interestingly, the interaction term was non-significant in both cases where the best individual regulator was not part of our regulatory mechanism, meaning that the single best regulator was outperformed by a linear combination of other regulators. Such examples demonstrate how systems biology approaches have a better power to dissect regulatory complexity of biological systems than traditional approaches [1, 41, 42]. They also show that systems biology is able to better model the 'real world' as QTL analysis of quantitative traits typically identifies numerous genetic loci, suggesting the involvement of numerous genes.
A particularly appealing feature of regression-based networks is their ability to predict expression of genes based on the expression of transcription factors. We have used this to quantify the stability and predictive power of the network, but also to study module-conservation between experiments in Populus and in Arabidopsis. Several interesting predictions were found when studying modules that are co-expressed and correctly predicted using bootstrap networks, but that lose their predictability in particular experiment when these are entirely removed before network inference. We have already mentioned the nucleosome assembly modules that are predicted to be regulated by histone H2A proteins. The drought response profiles in these modules cannot be predicted by networks not trained on drought stress data. Another module (characterized by motifs AS~TATA-box, AT~TATA-box, BN~TATA-box, PC~Box_4 and ZM~TATA-box) was affected by the removal of the budset data and is predicted to be regulated by factor 725612, a known cell death regulator. The module characterized by motif OS~TGGCA looses predictability without the dynamics of leaf growth-dataset, and is predicted to be controlled by a cell growth regulator (protein id 639804). Genes in this module are also over-represented for carbon fixation. Prediction is a central theme in this study, and we strongly believe that predictive models have a lot to offer experimental biology as hypotheses generators.
The complete and correct regulatory network of an organism cannot be reverse-engineered from a limited collection of gene expression data. However, we believe that such models represent a powerful starting point for further analysis as both hypothesis generation and descriptive tools. The hubs in our network (Table 1) thus represent attractive candidates for Chip-seq analysis, functional knock-down studies and regulon engineering. The network we present here only reflects the best regulators of each module. However, behind each module in this network there is a ranked list of regulatory mechanisms (Additional file 4), and as we have seen through bootstrap analysis, the ranking of these lists is not written in stone. In the future one might hope that additional, and higher quality, data (e.g. RNA-Seq) will enable creation of more robust network models that more accurately reflect the underlying biological truth. Obviously, even a perfect network inference method cannot be better than the data it is modeled on (junk in, junk out). Another route to more reliable networks lies in combining computational inference with experimental testing in an iterative modeling approach. Several studies have shown how systematic perturbation of critical pathway components can be used to refine network representations [43, 44]. In plants, the lignin systems-project is taking this approach to model the lignin biosynthesis pathway http://www.ligninsystems.org. Other sources of information may also be integrated into the network, but were not considered here, including epigenetic signatures such as nucleosome positioning and methylation patterns , predicted binding site strength and transcription factor binding site preference , and miRNA regulation .
We have outlined a systems biology model of the regulatory network of Populus leaves. The approach goes beyond previous analyses of Populus transcriptome data by systematically considering interactions between transcription factors, leading us to predict new regulators of leaf development not found by traditional genomics methods. These regulators orchestrate the transcriptional program in a synergistic or competitive manner, and thus constitute non-obvious targets for further analysis. The model is robust when applied to predict expression levels in new data, and reveals conserved and diverged regulation both in different conditions within Populus and between Populus and Arabidopsis.
Street et al.  identified 562 leaf-specific Populus genes that were profiled using Populus cDNA microarrays in 465 different experimental conditions (data available in UPSC-BASE  and in Additional file 1). These experiments included budset (74 conditions), biotic infection (21), weather dependent gene expression (33), CBF over express/freezing tolerance (17), seasonal leaf growth (37), elevated CO2 (12), PsbS antisense (17), leaf primordia (32), dynamics of leaf growth (21), P. nigra rust infection (24), herbivory/jasmonic acid (36), drought stress (57) and various other conditions (84).
Sequence motifs and promoters
We created a database of 312 non-redundant plant-related transcription factor binding sites from PlantCare , Transfac  and JASPAR plantae . From the initial set of 470 motifs, we iteratively identified the two most similar motifs and removed the longest until no pair had a MotifComparison  distance bellow 0.3. 2000 bp Populus promoters were taken from the PopGenIE online resource http://popgenie.org/ and MotifScanner  was used to scan these promoters for occurrences of the motifs. MotifScanner was run with a second order background model created from all Populus promoters and an a prior probability of finding one instance of the motif equal to 0.2. 307 motifs had hits to at least five genes in the leaf expression dataset.
Transcriptional module discovery
We have previously developed a method for discovering transcriptional modules that uses rule-based machine learning to find combinations of motifs that are predictive of co-expression [29–31]. Here we used this approach to find modules within the leaf experiments. For each gene, we identified all co-expressed genes at different levels of expression similarity and applied the rule learning method to find motif combinations explaining this co-expression pattern. Two genes were deemed co-expressed if their expression profiles had a Spearman correlation coefficient higher than a threshold (calculated based only on the experiments where both genes had measured expression). This threshold was varied from 0.50 to 0.95 in steps of 0.05. Only motif combinations with at least five genes over the co-expression threshold, and no more than 50 genes below the threshold, were considered. P-values for the overlap between genes with the motif combination and co-expressed genes were computed using the hyper-geometric distribution, and only FDR-significant rules (controlled at 0.05) were retained.
Gene function annotations were taken from KEGG  and Gene Ontology (GO) . Since GO do not provide annotations for Populus genes, we took annotations from the five closest proteins in the GO database with BLAST E-value less than 1E-6 or, if BLAST gave no hits, PSI-BLAST E-value less than 1E-6. Using the hyper-geometric distribution, we computed p-values for all annotations (at all levels in GO) with assignments to at least two co-expressed genes in a module, and retained all FDR-significant annotations (controlled at 0.05). We also performed randomization experiments by randomly shuffling promoters among the genes to create 1000 randomized data sets, and then performing module discovery and annotation analysis of each of these.
We used a least square regression model to infer regulators of each transcriptional module. Here, the expression of a module m i was modeled as the weighted sum of the expression of a set of transcription factors m i = β0 + ∑j ∈ Rβ j t j + ∑j, k ∈ R, j < kβ jk t j t k , where t j is the transcription factor with index j and R is the set of transcription factor indices. The best regulators of each module were found by estimating the performance of different sets of possible regulators R. Performance was quantified as the correlation between observed (i.e. measured by cDNA microarray) and predicted expression during cross validation (five iterations of 5-fold cross validation). The order of R was iteratively increased from single transcription factors (order 1), to pairs of transcription factors (order 2), etc. The best set of regulators of order n was selected as the final regulatory mechanism of the module if no set of regulators of order n+1 could predict expression of the module significantly better. Significance was determined by using the Bonferroni corrected p-value (i.e. multiplied by the number of transcription factor combinations tested) calculated using a t-test for the difference between two non-independent Pearson correlations . The expression profile of a module was defined as the concatenation of the expression profiles of each co-expressed gene in the module. The regulatory networks was constructed by using transcription factors and modules as nodes, and drawing an edge between a transcription factor and a module if the transcription factor was part of the best regulatory mechanism for that module.
We drew 100 bootstrap datasets from the original 465 conditions in the leaf dataset (i.e. 100 samples of 465 conditions drawn with replacement) and inferred networks from each of these datasets. The regression model of each module was then used to predict the expression in non-sampled conditions for the co-expressed genes in that module. For each gene, predicted expression values from each condition were averaged across the bootstrap samples, and correlation between observed and predicted expression was calculated. The resulting correlation for a gene was thus only calculated for conditions that were not part of at least one bootstrap sample. We also investigated the stability of the regulations by calculating the fraction of bootstrapped networks that contained each edge in the original network.
Arabidopsis data was taken from the AtGenExpress resource: development (237 conditions) , abiotic stress (298) , biotic stress (108) and light (48). We mapped our Populus proteins to the closest proteins in Arabidopsis as detected by BLAST . We then used regression models trained on the Populus expression data to predict expression in Arabidopsis.
Thanks to Patrik Rydén for advice on some statistical analysis. The High Performance Computing Center North (HPC2N) was utilized for computer intensive calculations.
This work was supported by funds from The Swedish Research Council (VR) and The Swedish Governmental Agency for Innovation Systems (VINNOVA) through the UPSC Berzelii Centre for Forest Biotechnology, and from the Kempe foundation.
- Fischbach MA, Krogan NJ: The next frontier of systems biology: higher-order and interspecies interactions. Genome Biol. 11 (5): 208Google Scholar
- Gatherer D: So what do we really mean when we say that systems biology is holistic. BMC Systems Biology. 2010, 4: 22-10.1186/1752-0509-4-22.PubMedPubMed CentralView ArticleGoogle Scholar
- Carro MS, Lim WK, Alvarez MJ, Bollo RJ, Zhao X, Snyder EY, Sulman EP, Anne SL, Doetsch F, Colman H, et al: The transcriptional network for mesenchymal transformation of brain tumours. Nature. 2009, 463 (7279): 318-325. 10.1038/nature08712.PubMedPubMed CentralView ArticleGoogle Scholar
- Carrera J, Rodrigo G, Jaramillo A, Elena SF: Reverse-engineering the Arabidopsis thaliana transcriptional network under changing environmental conditions. Genome Biol. 2009, 10 (9): R96-10.1186/gb-2009-10-9-r96.PubMedPubMed CentralView ArticleGoogle Scholar
- Needham CJ, Manfield IW, Bulpitt AJ, Gilmartin PM, Westhead DR: From gene expression to gene regulatory networks in Arabidopsis thaliana. BMC Syst Biol. 2009, 3: 85-10.1186/1752-0509-3-85.PubMedPubMed CentralView ArticleGoogle Scholar
- Krishnan A, Pereira A: Integrative approaches for mining transcriptional regulatory programs in Arabidopsis. Brief Funct Genomic Proteomic. 2008, 7 (4): 264-274. 10.1093/bfgp/eln035.PubMedView ArticleGoogle Scholar
- Jansson S, Douglas CJ: Populus: a model system for plant biology. Annu Rev Plant Biol. 2007, 58: 435-458. 10.1146/annurev.arplant.58.032806.103956.PubMedView ArticleGoogle Scholar
- Tuskan GA, Difazio S, Jansson S, Bohlmann J, Grigoriev I, Hellsten U, Putnam N, Ralph S, Rombauts S, Salamov A, et al: The genome of black cottonwood, Populus trichocarpa (Torr. & Gray). Science. 2006, 313 (5793): 1596-1604. 10.1126/science.1128691.PubMedView ArticleGoogle Scholar
- Sjodin A, Wissel K, Bylesjo M, Trygg J, Jansson S: Global expression profiling in leaves of free-growing aspen. BMC Plant Biol. 2008, 8: 61-10.1186/1471-2229-8-61.PubMedPubMed CentralView ArticleGoogle Scholar
- Sjodin A, Bylesjo M, Skogstrom O, Eriksson D, Nilsson P, Ryden P, Jansson S, Karlsson J: UPSC-BASE--Populus transcriptomics online. Plant J. 2006, 48 (5): 806-817. 10.1111/j.1365-313X.2006.02920.x.PubMedView ArticleGoogle Scholar
- Sjödin A, Street NR, Sandberg G, Gustafsson P, Jansson S: The Populus Genome Integrative Explorer (PopGenIE): a new resource for exploring the Populus genome. New Phytol. 2009, 182 (4): 1013-1025.PubMedView ArticleGoogle Scholar
- Grönlund A, Bhalerao RP, Karlsson J: Modular gene expression in Poplar: a multilayer network approach. New Phytol. 2009, 181 (2): 315-322.PubMedView ArticleGoogle Scholar
- Street NR, Sjodin A, Bylesjo M, Gustafsson P, Trygg J, Jansson S: A cross-species transcriptomics approach to identify genes involved in leaf development. BMC Genomics. 2008, 9: 589-10.1186/1471-2164-9-589.PubMedPubMed CentralView ArticleGoogle Scholar
- Quesada T, Li Z, Dervinis C, Li Y, Bocock PN, Tuskan GA, Casella G, Davis JM, Kirst M: Comparative analysis of the transcriptomes of Populus trichocarpa and Arabidopsis thaliana suggests extensive evolution of gene expression regulation in angiosperms. New Phytol. 2008, 180 (2): 408-420. 10.1111/j.1469-8137.2008.02586.x.PubMedView ArticleGoogle Scholar
- Shi R, Sun YH, Li Q, Heber S, Sederoff R, Chiang VL: Towards a systems approach for lignin biosynthesis in Populus trichocarpa: transcript abundance and specificity of the monolignol biosynthetic genes. Plant Cell Physiol. 51 (1): 144-163. 10.1093/pcp/pcp175.Google Scholar
- Drost DR, Benedict CI, Berg A, Novaes E, Novaes CR, Yu Q, Dervinis C, Maia JM, Yap J, Miles B, et al: Diversification in the genetic architecture of gene expression and transcriptional networks in organ differentiation of Populus. Proc Natl Acad Sci USA. 107 (18): 8492-8497. 10.1073/pnas.0914709107.Google Scholar
- Bansal M, Belcastro V, Ambesi-Impiombato A, di Bernardo D: How to infer gene networks from expression profiles. Mol Syst Biol. 2007, 3: 78-PubMedPubMed CentralView ArticleGoogle Scholar
- Marbach D, Prill RJ, Schaffter T, Mattiussi C, Floreano D, Stolovitzky G: Revealing strengths and weaknesses of methods for gene network inference. Proc Natl Acad Sci USA. 107 (14): 6286-6291. 10.1073/pnas.0913357107.Google Scholar
- Styczynski MP, Stephanopoulos G: Overview of computational methods for the inference of gene regulatory networks. Computers & Chemical Engineering. 2005, 29 (3): 519-534.View ArticleGoogle Scholar
- Huynh-Thu VA, Irrthum A, Wehenkel L, Geurts P: Inferring Regulatory Networks from Expression Data Using Tree-Based Methods. PLoS ONE. 2010, 5 (9): e12776-10.1371/journal.pone.0012776.PubMedPubMed CentralView ArticleGoogle Scholar
- Swain MT, Mandel JJ, Dubitzky W: Comparative study of three commonly used continuous deterministic methods for modeling gene regulation networks. BMC Bioinformatics. 11: 459-10.1186/1471-2105-11-459.Google Scholar
- Wilczynski B, Dojer N: BNFinder: exact and efficient method for learning Bayesian networks. Bioinformatics. 2009, 25 (2): 286-287. 10.1093/bioinformatics/btn505.PubMedPubMed CentralView ArticleGoogle Scholar
- Segal E, Widom J: From DNA sequence to transcriptional behaviour: a quantitative approach. Nat Rev Genet. 2009, 10 (7): 443-456. 10.1038/nrg2591.PubMedPubMed CentralView ArticleGoogle Scholar
- Nero D, Katari MS, Kelfer J, Tranchina D, Coruzzi GM: In silico evaluation of predicted regulatory interactions in Arabidopsis thaliana. BMC Bioinformatics. 2009, 10: 435-10.1186/1471-2105-10-435.PubMedPubMed CentralView ArticleGoogle Scholar
- Pilpel Y, Sudarsanam P, Church GM: Identifying regulatory networks by combinatorial analysis of promoter elements. Nat Genet. 2001, 29 (2): 153-159. 10.1038/ng724.PubMedView ArticleGoogle Scholar
- Beer MA, Tavazoie S: Predicting gene expression from sequence. Cell. 2004, 117 (2): 185-198. 10.1016/S0092-8674(04)00304-6.PubMedView ArticleGoogle Scholar
- Segal E, Shapira M, Regev A, Pe'er D, Botstein D, Koller D, Friedman N: Module networks: identifying regulatory modules and their condition-specific regulators from gene expression data. Nat Genet. 2003, 34 (2): 166-176. 10.1038/ng1165.PubMedView ArticleGoogle Scholar
- Segal E, Yelensky R, Koller D: Genome-wide discovery of transcriptional modules from DNA sequence and gene expression. Bioinformatics. 2003, 19 (Suppl 1): I273-I282. 10.1093/bioinformatics/btg1038.PubMedView ArticleGoogle Scholar
- Andersson CR, Hvidsten TR, Isaksson A, Gustafsson MG, Komorowski J: Revealing cell cycle control by combining model-based detection of periodic expression with novel cis-regulatory descriptors. BMC Syst Biol. 2007, 1: 45-10.1186/1752-0509-1-45.PubMedPubMed CentralView ArticleGoogle Scholar
- Hvidsten TR, Wilczynski B, Kryshtafovych A, Tiuryn J, Komorowski J, Fidelis K: Discovering regulatory binding-site modules using rule-based learning. Genome Res. 2005, 15 (6): 856-866. 10.1101/gr.3760605.PubMedPubMed CentralView ArticleGoogle Scholar
- Wilczynski B, Hvidsten TR, Kryshtafovych A, Tiuryn J, Komorowski J, Fidelis K: Using local gene expression similarities to discover regulatory binding site modules. BMC Bioinformatics. 2006, 7: 505-10.1186/1471-2105-7-505.PubMedPubMed CentralView ArticleGoogle Scholar
- Rouault H, Mazouni K, Couturier L, Hakim V, Schweisguth F: Genome-wide identification of cis-regulatory motifs and modules underlying gene coregulation using statistics and phylogeny. Proc Natl Acad Sci USA. 107 (33): 14615-14620. 10.1073/pnas.1002876107.Google Scholar
- Ihmels J, Friedlander G, Bergmann S, Sarig O, Ziv Y, Barkai N: Revealing modular organization in the yeast transcriptional network. Nat Genet. 2002, 31 (4): 370-377.PubMedGoogle Scholar
- Reiss DJ, Baliga NS, Bonneau R: Integrated biclustering of heterogeneous genome-wide datasets for the inference of global regulatory networks. BMC Bioinformatics. 2006, 7: 280-10.1186/1471-2105-7-280.PubMedPubMed CentralView ArticleGoogle Scholar
- Liu ZQ, Gao J, Dong AW, Shen WH: A truncated Arabidopsis NUCLEOSOME ASSEMBLY PROTEIN 1, AtNAP1;3T, alters plant growth responses to abscisic acid and salt in the Atnap1;3-2 mutant. Mol Plant. 2009, 2 (4): 688-699. 10.1093/mp/ssp026.PubMedView ArticleGoogle Scholar
- Das D, Pellegrini M, Gray JW: A primer on regression methods for decoding cis-regulatory logic. PLoS Comput Biol. 2009, 5 (1): e1000269-10.1371/journal.pcbi.1000269.PubMedPubMed CentralView ArticleGoogle Scholar
- Efron B: Bootstrap Methods: Another Look at the Jackknife. The Annals of Statistics. 1979, 7 (1): 1-26. 10.1214/aos/1176344552.View ArticleGoogle Scholar
- Carrera J, Rodrigo G, Jaramillo A: Model-based redesign of global transcription regulation. Nucleic Acids Res. 2009, 37 (5): e38-10.1093/nar/gkp022.PubMedPubMed CentralView ArticleGoogle Scholar
- Tirosh I, Barkai N, Verstrepen KJ: Promoter architecture and the evolvability of gene expression. J Biol. 2009, 8 (11): 95-10.1186/jbiol204.PubMedPubMed CentralView ArticleGoogle Scholar
- Schmid M, Davison TS, Henz SR, Pape UJ, Demar M, Vingron M, Scholkopf B, Weigel D, Lohmann JU: A gene expression map of Arabidopsis thaliana development. Nat Genet. 2005, 37 (5): 501-506. 10.1038/ng1543.PubMedView ArticleGoogle Scholar
- Bhardwaj N, Carson MB, Abyzov A, Yan KK, Lu H, Gerstein MB: Analysis of combinatorial regulation: scaling of partnerships between regulators with the number of governed targets. PLoS Comput Biol. 6 (5): e1000755-10.1371/journal.pcbi.1000755.Google Scholar
- Opper M, Sanguinetti G: Learning combinatorial transcriptional dynamics from gene expression data. Bioinformatics. 26 (13): 1623-1629. 10.1093/bioinformatics/btq244.Google Scholar
- Ideker T, Thorsson V, Ranish JA, Christmas R, Buhler J, Eng JK, Bumgarner R, Goodlett DR, Aebersold R, Hood L: Integrated Genomic and Proteomic Analyses of a Systematically Perturbed Metabolic Network. Science. 2001, 292 (5518): 929-934. 10.1126/science.292.5518.929.PubMedView ArticleGoogle Scholar
- Lorenz DR, Cantor CR, Collins JJ: A network biology approach to aging in yeast. Proc Natl Acad Sci USA. 2009, 106 (4): 1145-1150. 10.1073/pnas.0812551106.PubMedPubMed CentralView ArticleGoogle Scholar
- Feng S, Cokus SJ, Zhang X, Chen PY, Bostick M, Goll MG, Hetzel J, Jain J, Strauss SH, Halpern ME, et al: Conservation and divergence of methylation patterning in plants and animals. Proc Natl Acad Sci USA. 107 (19): 8689-8694. 10.1073/pnas.1002720107.Google Scholar
- Segal E, Raveh-Sadka T, Schroeder M, Unnerstall U, Gaul U: Predicting expression patterns from regulatory sequence in Drosophila segmentation. Nature. 2008, 451 (7178): 535-540. 10.1038/nature06496.PubMedView ArticleGoogle Scholar
- Elemento O, Slonim N, Tavazoie S: A universal framework for regulatory element discovery across all genomes and data types. Mol Cell. 2007, 28 (2): 337-350. 10.1016/j.molcel.2007.09.027.PubMedPubMed CentralView ArticleGoogle Scholar
- Lescot M, Dehais P, Thijs G, Marchal K, Moreau Y, Van de Peer Y, Rouze P, Rombauts S: PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 2002, 30 (1): 325-327. 10.1093/nar/30.1.325.PubMedPubMed CentralView ArticleGoogle Scholar
- Wingender E, Dietze P, Karas H, Knuppel R: TRANSFAC: a database on transcription factors and their DNA binding sites. Nucleic Acids Res. 1996, 24 (1): 238-241. 10.1093/nar/24.1.238.PubMedPubMed CentralView ArticleGoogle Scholar
- Bryne JC, Valen E, Tang MH, Marstrand T, Winther O, da Piedade I, Krogh A, Lenhard B, Sandelin A: JASPAR, the open access database of transcription factor-binding profiles: new content and tools in the 2008 update. Nucleic Acids Res. 2008, D102-106. 36 DatabaseGoogle Scholar
- Thijs G, Moreau Y, De Smet F, Mathys J, Lescot M, Rombauts S, Rouze P, De Moor B, Marchal K: INCLUSive: integrated clustering, upstream sequence retrieval and motif sampling. Bioinformatics. 2002, 18 (2): 331-332. 10.1093/bioinformatics/18.2.331.PubMedView ArticleGoogle Scholar
- Kanehisa M, Goto S: KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28 (1): 27-30. 10.1093/nar/28.1.27.PubMedPubMed CentralView ArticleGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25 (1): 25-29. 10.1038/75556.PubMedPubMed CentralView ArticleGoogle Scholar
- Howell DC: Statistical Methods for Psychology. Wadsworth CENGAGE Learning; 1997.Google Scholar
- Kilian J, Whitehead D, Horak J, Wanke D, Weinl S, Batistic O, D'Angelo C, Bornberg-Bauer E, Kudla J, Harter K: The AtGenExpress global stress expression data set: protocols, evaluation and model data analysis of UV-B light, drought and cold stress responses. Plant J. 2007, 50 (2): 347-363. 10.1111/j.1365-313X.2007.03052.x.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.