Network component analysis provides quantitative insights on an Arabidopsis transcription factor-gene regulatory network
© Misra and Sriram; licensee BioMed Central Ltd. 2013
Received: 24 August 2013
Accepted: 5 November 2013
Published: 14 November 2013
Gene regulatory networks (GRNs) are models of molecule-gene interactions instrumental in the coordination of gene expression. Transcription factor (TF)-GRNs are an important subset of GRNs that characterize gene expression as the effect of TFs acting on their target genes. Although such networks can qualitatively summarize TF-gene interactions, it is highly desirable to quantitatively determine the strengths of the interactions in a TF-GRN as well as the magnitudes of TF activities. To our knowledge, such analysis is rare in plant biology. A computational methodology developed for this purpose is network component analysis (NCA), which has been used for studying large-scale microbial TF-GRNs to obtain nontrivial, mechanistic insights. In this work, we employed NCA to quantitatively analyze a plant TF-GRN important in floral development using available regulatory information from AGRIS, by processing previously reported gene expression data from four shoot apical meristem cell types.
The NCA model satisfactorily accounted for gene expression measurements in a TF-GRN of seven TFs (LFY, AG, SEPALLATA3 [SEP3], AP2, AGL15, HY5 and AP3/PI) and 55 genes. NCA found strong interactions between certain TF-gene pairs including LFY → MYB17, AG → CRC, AP2 → RD20, AGL15 → RAV2 and HY5 → HLH1, and the direction of the interaction (activation or repression) for some AGL15 targets for which this information was not previously available. The activity trends of four TFs - LFY, AG, HY5 and AP3/PI as deduced by NCA correlated well with the changes in expression levels of the genes encoding these TFs across all four cell types; such a correlation was not observed for SEP3, AP2 and AGL15.
For the first time, we have reported the use of NCA to quantitatively analyze a plant TF-GRN important in floral development for obtaining nontrivial information about connectivity strengths between TFs and their target genes as well as TF activity. However, since NCA relies on documented connectivity information about the underlying TF-GRN, it is currently limited in its application to larger plant networks because of the lack of documented connectivities. In the future, the identification of interactions between plant TFs and their target genes on a genome scale would allow the use of NCA to provide quantitative regulatory information about plant TF-GRNs, leading to improved insights on cellular regulatory programs.
Gene expression is a complex process regulated by the interactions of proteins and other molecules with genes. This regulation occurs at multiple levels, giving rise to gene regulatory networks (GRNs) that define the regulatory programs for the expression of specific genes in response to specific cues . One of the biggest challenges of systems biology is deciphering the organization of GRNs [2, 3]. This task is further complicated by feedback- and feedforward-type interactions of a multitude of genes and their protein products upon themselves and others. GRNs are usually modeled as graphs with nodes representing system components (e.g. molecules) and edges indicating interactions between components [1, 4, 5]. Various methodologies have been developed for the analysis of GRNs including directed graphs, Boolean networks, Bayesian networks and differential equations [2, 6–11]. An important subset of GRNs models gene expression as a result of the action of transcription factors (TFs) upon their target genes. In these models, directed edges from TFs to their target genes represent transcriptional regulation, and constitute a hierarchical network governing gene expression [2, 12]. The reconstruction of TF-GRNs involves the identification of genes that encode the TFs and the identification of the target genes of the TFs.
There is a considerable amount of information available on TF-gene interactions in microbes which is housed in databases. For example, RegulonDB and DBTBS are extensively curated databases containing information on transcriptional regulation in the bacteria Escherichia coli and Bacillus subtilis respectively [13, 14]. The RegPrecise database contains similar information for many other prokaryotes , as does the YEASTRACT database for Saccharomyces cerevisiae. The availability of such resources permits accurate reconstruction of TF-GRNs, and consequent network analyses to obtain insights on regulatory capabilities of the organism of interest. For plants, such information is comparatively sparse, with most regulatory studies directed at inferring GRNs in isolated organs such as roots or leaves, or processes such as development or abiotic stress response [9, 17, 18]. Large-scale TF-gene interaction data are only available for Arabidopsis thaliana and housed in the Arabidopsis Gene Regulatory Information Server (AGRIS) .
Although the establishment of TF-GRN connectivity (i.e. which TF regulates which gene) is very useful, the information contained in such connectivity maps is binary and not quantitative. Understanding quantitative changes in gene expression would provide deeper insights into gene regulation and perhaps even enable predictive modeling of cellular regulatory programs. This would, however, require significant mathematical processing of high-throughput gene expression datasets . Under a given condition, gene expression would depend on the strength of the interaction between a TF and its target gene as well as the activity of the TF at that condition. Therefore, given the connectivity of a TF-GRN and gene expression values under a set of conditions, the next set of questions that need to be answered are: (i) Is it possible to obtain connectivity strengths (CS) of TF-gene interactions for the network and (ii) Can we quantify how TF activity varies across conditions? Estimating the CS between a TF and its target gene may be possible computationally by determining the decrease in free energy for binding between the TF and the DNA region of the target gene it binds to [21, 22]. A higher free energy change would indicate stronger binding and a lower free energy change weaker binding [21, 23]. However, thermodynamic calculations for determining changes in free energy are nontrivial and would require knowledge of binding thermodynamics of many TFs and their target genes. The CS between a TF and a gene can also be determined experimentally by using binding assays for determining parameters such as the dissociation constant or changes in free energy and enthalpy . Although parameters derived from such TF-gene binding assays are available in some databases, it would be a laborious exercise to obtain these values for every TF-gene pair . For estimating changes in TF activity, experimental assays may be employed based on the binding of the active form of the TF with a target reporter molecule. However, such assays are only available for a limited number of TFs and would have to be conducted for each condition. Additionally, the experimental approaches for determining TF-gene CS and TF activities suffer from the drawback of being in vitro studies. Consequently, the values determined may not represent the in vivo interactions of the TFs and genes wherein multiple TFs can act on a single gene. It may appear that changes in the expression levels of the genes corresponding to the TFs could be used as surrogates for TF activities. However, a shortcoming of this approach is that TF activity could be considerably affected by post-transcriptional and post-translational modifications such as phosphorylation and acetylation, and can therefore, differ substantially from the expression levels of corresponding genes.
To deduce such quantitative information about TF-GRNs, researchers have developed methodologies like network component analysis (NCA) and regulatory element detection using correlation with expression (REDUCE) [26–29]. NCA, in particular, models gene expression to be the result of the connectivity strength between TF-gene pairs and TF activity . The strength of the TF-gene interaction indicates the extent of the control of a TF over the transcription of a target gene, whereas the TF activity quantifies how active the TF is in regulating its target genes either via activation or repression. NCA uses connectivity information about the underlying network and gene expression data to obtain non trivial information about TF activity and TF-gene connectivity strength. Because the TF activity provides a measure for the TF in its final state, it includes information about the post-transcriptional and post-translational modifications. Compared to experimental approaches for obtaining similar information, NCA allows the deduction of such important regulatory information by a much simpler approach involving the measurement of gene expression for the set of genes in a network. The other input for NCA, the connectivity between TFs and genes, is available for many organisms in databases. Consequently, NCA provides an additional layer of regulatory information without the use of sophisticated experimental measurements .
Here, [G]m×n is a matrix representing an experimental gene expression dataset consisting of the expression of m genes across n conditions; [log G]m×n is its log-transformed version. Similarly, [TFA]p×n is a matrix of the activities of p TFs across the n conditions; [log TFA]p×n is its log-transformed version. These two matrices are linked by [CS]m×p, which consists of the CS of p TFs on m genes.
The log-linear relationship used in NCA allows the benefits of linearization during the decomposition while capturing non-linear network behavior to a limited extent. Besides, since high-throughput gene expression data are usually expressed relative to a control condition, the log-linear relationship is convenient while working with relative gene expression data . The NCA decomposition is unique up to a scaling factor, when the [CS] and [TFA] matrices satisfy a set of criteria termed “NCA-compliance” criteria . The originally reported NCA algorithm  required the presence of as many gene expression data points as regulators for the decomposition. However, a more recent modification of that algorithm  permits the analysis of limited microarray datasets, thus widening the applicability of NCA. A detailed analysis of the original NCA algorithm and the modified algorithm are provided in the respective publications [26, 30].
NCA has been previously applied for the analysis of microbial and mammalian transcriptional networks. Liao et al.  first used NCA to study cell cycle regulation in S. cerevisiae, and specifically to quantify the activities of different TFs during various stages of the cell cycle, thus gaining insight on the regulatory roles of specific TFs at each stage. Kao et al.  investigated the effect of a glucose-to-acetate carbon source transition on the activity of TFs in E. coli. They observed specific trends in the changes in activities of several TFs (CRP, FadR, IclR, and Cra) important during this transition. In a further extension of this study, they investigated the growth lag that resulted by the deletion of the ppsA gene in E. coli during this carbon source transition . By using NCA, they deduced the activities of TFs that were affected by the deletion and proposed a mechanism for explaining the growth lag. A set of twin studies investigating the effect of the reactive nitrogen species, nitric oxide and S-nitrosoglutathione, on E. coli identified important TFs involved in response to the respective treatments [31, 32]. The first study identified 13 important TFs of which ten have not been previously documented to be involved in response to nitric oxide . The subsequent study with S-nitrosoglutathione identified four novel TFs (CysB, SF, FlhDC, and TTA) involved in response to the treatment . The use of NCA in combination with transcriptome data allowed the construction of models depicting the response process for both studies. Brynildsen et al. investigated the isobutanol response network in E. coli and identified the ArcA-ArcB system to be a major regulator of the response via a loss of quinone function . They also compared differences in TF activities in response to isobutanol with those seen for butanol and ethanol, and identified 6 TFs with differing activities for butanol, and 19 TFs with differing activities for ethanol compared to isobutanol. In another study , Buescher et al. performed genome wide TF-gene analysis of B. subtilis during a change in carbon substrate from glucose to malate and vice versa, and determined CS for 2900 TF-gene interactions. They deduced TF activities for 154 TFs out of which 127 TFs were found to change their activities significantly. Interestingly, many of these changes in TF activity were not seen at the mRNA level thus implicating the role of posttranslational modifications for the changes in TF activities. In mammalian systems, Sriram et al. studied the effect of overexpressing the glycerol kinase gene in rat hepatoma cells using a network of 62 genes and 9 TFs . They found an increase in the TF activity for 7 of the TFs (ChREBP, Sp1, HNF1α, HNF4α, PPARα, LXRα, and glucocorticoid receptor [GR]) and a decrease in activity for the remaining 2 TFs (SREBP1a and CEBPβ). The increased activity of GR was hypothesized to be a result of the moonlighting nature of the glycerol kinase enzyme . Sriram et al. experimentally verified the NCA-deduced change in TF activity of GR in the glycerol kinase-overexpressing cell line, thus demonstrating the power of NCA for deducing TF activities from gene expression data in a mammalian network. In a recent study , Tran et al. studied the TFs directly downstream of PTEN (phosphatase and tensin homologue deleted on chromosome 10), which is an important tumor suppressor gene. They identified 20 TFs whose activities were altered significantly by the expression of PTEN even when the mRNA levels of the corresponding genes did not alter significantly. They found many of the identified TFs varied in murine and human cancer models, and provided a signature for identifying the status of PTEN in cancers caused by PTEN loss.
In this article, we report the application of NCA on a plant TF-GRN using available regulatory information from AGRIS. Starting with a set of TFs known to be important in floral development, we mined AGRIS to establish a network consisting of confirmed TF-gene connectivities in this developmental event. We used previously published gene expression data  for four types of cells isolated from the shoot apical meristem, which is known to initiate the growth of floral organs. By using the connectivity information and gene expression datasets, we used NCA to deduce activities for the NCA-compliant TFs, and numerical values of CS between the TFs and their target genes. To the best of our knowledge, this is the first study to apply NCA to dissect a plant TF-GRN.
NCA deduces the strengths of TF-gene interactions
Gene expression levels simulated by NCA agree well with the originally measured gene expression levels
TF activities deduced for LFY, AG, HY5 and AP3/PI agree well with expression levels of genes encoding these TFs
Normalized plots of TF activities and gene expression values showed a good fit for LFY, AG, HY5 and AP3
TF-GRNs, which model interactions between TFs and their target genes, are an important class of cellular networks that define regulatory programs leading to gene expression [2, 12]. TF-GRNs provide Boolean information about the regulation of genes by TFs, with meticulously compiled data available in databases like RegulonDB, YEASTRACT and AGRIS [13, 16, 19]. To deduce further quantitative information about the connectivities between TFs and their target genes, methodologies such as NCA and REDUCE have been developed [26, 29]. Given the underlying network connectivity information, NCA can provide information on the connectivity strength between a TF and its target gene as well as the TF activity by using gene expression data [26, 30, 40]. Through such nontrivial, quantitative information, NCA can provide important parameters about a TF-GRN. In this study, we sought to apply the NCA approach to analyze a network comprising TFs important for floral development and their targets using underlying connectivity information available in the AGRIS database.
Floral development is one of the best characterized processes in plants with multiple studies providing much information at the molecular genetic level [41–43]. The most widely used model for explaining the initial development of the organs of a flower is the ABC model and its variants . The model predicts floral development to result from the concerted action of multiple TF-encoding genes. For this study, we constructed a plant TF-GRN consisting of ten TFs, known to be involved in floral development, (LFY, AG, SEPALLATA3 (SEP3), AP2, AGL15, HY5, AP3/PI, FD, WUS and BLR) and 57 target genes with verified interactions obtained from AGRIS. LFY is known to be a master TF that regulates important events in the transition from vegetative to reproductive growth, and has another important role in the activation of floral homeotic genes [44–46]. Some of its downstream targets are known to be TFs that are important in flower morphogenesis. The other TFs included in our original network are important factors in floral development: AG, SEP3 and AGL15 are MADS domain TFs; AP2 belongs to the AP2/EREBP (ethylene responsive element binding protein) class of TFs; HY5 and FD are basic leucine zipper TFs that regulate flower development; AP3/PI is a member of the NAC TF family that is expressed in floral primordia and WUS and BLR are homeobox TFs . We were unable to include some of the other TFs (AP1, FT and AGL20) important in the process due to a lack of sufficient confirmed targets for them in AGRIS for NCA compliance. We used gene expression data from a study by Yadav et al.  that analyzed the expression patterns across four different types of cells (named CLV3n, CLV3p, FILp and WUSp) isolated from shoot apical meristems of A. thaliana. The study isolated protoplasts of the cells by using fluorescent markers unique to them, and revealed a strong expression of the LFY gene across all cell types.
During preparation for NCA, three of the TFs (FD, WUS and BLR) and their corresponding gene connections had to be removed as they were not NCA-compliant. The final NCA-compliant network consisted of the remaining 7 TFs and 55 genes. For the NCA, we assumed same connectivity strengths for TF with their target genes across all cell lines, which is a reasonable assumption. NCA provided CS for all TF-gene pairs. However, after NCA decomposition, the CS needed to be checked for their signs (a positive sign signifies activation and a negative sign signifies repression). This is done by comparing the CS with the initial connectivity matrix, and especially the connectivity directions of well-established TF-gene pairs. We found that the TF activities and CS for the AG, HY5, SEP3 and AP2 TFs needed to be corrected for their signs. The TF-gene pairs showing strong CS represent strong binding between a TF and its target. However, many TF-gene pairs showed very low CS, so that their documented regulatory connection would be worth re-examining . Interestingly, AGRIS did not list the direction of interaction between AGL15 and four of the genes regulated by it (AGL22, AGL25, EDF4 and RAV2). NCA deduced AGL15 to be a strong repressor of AGL22, strong activator of RAV2, moderate activator for AGL25 and very weak repressor for EDF4. Thus, given verified information about the sign of a TF-gene interaction, NCA can deduce whether the TF is an activator or repressor of other target genes based on gene expression data. We should point out though that the strength of NCA is the deduction of quantitative information about a TF-GRN based on verified information about the underlying connections and gene expression data for the network. AGL22, also known as Short Vegetative Phase (SVP) encodes a TF that can repress flowering time in addition to other genes AGL15, AGL18 and FLM[48–50]. Based on our NCA, we determined that AGL22 is repressed much more strongly by AGL15 compared to SEP3. Interestingly, though, the gene expression of AGL22 increased several-fold compared to the control across all four cell types. This might be explained by the observation that even though the TF activity of SEP3 increases relative to the control, the TF activity of AGL15 is reduced compared to the control by a similar extent. As AGL15 controls the repression of AGL22 more strongly compared to SEP3, the gene expression of AGL22 compared to the control increases. Two other genes, HLH1 and RD20, are regulated by the same TFs, HY5 (activation) and AP2 (repression). NCA determined HLH1 to have similar connectivity strengths to both HY5 and AP2 but of opposite signs while HLH1 gene expression was found to be slightly higher compared to the control strain. This could be because of the slightly higher TF activity of HY5 compared to AP2 as deduced by NCA. RD20, on the other hand, was found to be mildly repressed across the four cell types compared to the control. This could be because it is more strongly repressed by AP2 compared to activation by HY5.
Of the different TFs included in our study, LFY plays the role of master regulator during floral development. Out of the direct targets of LFY included in our network, MYB17 or late meristem identity 2 is very important in meristem identity transition . MYB17 was found to be very strongly activated by LFY. This, combined with high TF activity of LFY would explain the high expression levels seen for the MYB17 gene from mRNA analysis. We were unable to include AP1, which is another important TF in the meristem identity pathway that is known to interact in a positive feedback network with LFY and MYB17. We can, however, deduce that the AP1 TF would have higher activity across the four cell types compared to the control based on strong activities of LFY and MYB17. In fact, the reproductive phase in Arabidopsis involves the transition of the SAM to an inflorescence meristem and then to a floral meristem . The floral meristem identity proteins in Arabidopsis  include the TFs that were found to be upregulated from our analysis (LFY and SEP3) which seems to indicate that the cells were isolated from a floral and not a vegetative meristem.
We compared the TF activities obtained by NCA with the expression values for their corresponding genes. TF activities can in general be expected to be proportional to the expression levels of the corresponding genes. However, TFs that need to undergo extensive post-translational modification to be active can be exceptions to this expected trend. Our analysis showed that the profiles of TF activities obtained from NCA compared well with the expression levels of the genes coding for these TFs in the case of the majority of TFs (LFY, AG, HY5, AP3/PI and SEP3 (in two out of four cell types). However AP2 and AGL15 are exceptions. The discrepancy for AP2 and AGL15 could quite possibly be because of the large error in the measurement of the microarray replicates leading to problems with the NCA. A repeat of the gene expression analysis with better control on the replicates may provide a better answer to this. If a discrepancy is still observed, this would indicate a change in TFs due to post-transcriptional and post-translational modifications. NCA thus allows the generation of newer hypotheses relating to the conversion of a gene product to an active TF based on how well the gene expression results agree with the deduced activities of their corresponding TFs. As a further step, we compared normalized values for both, using maximum or minimum values for TF activity or gene expression across the four cell types to allow better comparison between them. We found a very good correlation for LFY; decent matches for AG, SEP3, HY5 and AP3/PI; and poor matches for AP2 and AGL15 from this analysis.
The application of NCA to microbial and mammalian systems has provided interesting insights into gene regulation by TFs. As previously described, the applications of NCA to microbial systems include the following: (i) investigation of TF changes during cell cycle regulation in S. cerevisiae (ii) analysis of changes in TF activities in E. coli during the change from a glycolytic carbon source (glucose) to a gluconeogenic carbon source (acetate)  (iii) studying the effects of reactive nitrogen species on a TF network in E. coli[31, 32] (iv) identification of TFs important in the isobutanol response network in E. coli and (v) determining TF-gene interactions in B. subtilis during a carbon source transition from glucose to malate and vice-versa , Applications of NCA to mammalian systems are more recent (i) studying the effects of overexpression of the glycerol kinase gene in rat hepatoma cells  and (ii) identifying TFs with altered activity in response to PTEN expression .
These studies of TF-GRNs have revealed the strengths of NCA in providing insights about the regulatory aspects of a system given the basic structural information about the underlying network. In the case of plants, there is lesser information available about TF-gene interactions. The AtRegNet database from AGRIS, which is the most comprehensive resource for such information, contains 768 confirmed TF-gene interactions for 46 TFs in A. thaliana, which is estimated to contain more than 1700 TFs . In our NCA of a network derived from AGRIS, the original network consisting of 10 TFs and 57 genes reduced to 7 TFs and 55 genes for NCA compliance. This is because of the absence of sufficient regulatory information about the three TFs that had to be removed. NCA requires that any TF in a network regulate at least two genes. The availability of more information about TF-gene interactions would overcome this issue of NCA non-compliant TFs.
NCA uses gene expression data and underlying network connectivity during its analysis; consequently, the quantitative measures provided by NCA are dependent on the accuracy of the underlying network. For example, many of the genes considered in this study have unconfirmed interactions with other TFs. If any of these interactions were confirmed, the current NCA could be rerun to account for the effect of additional TFs on expression of the target genes. Thus, having correct prior connectivity information about a network would increase the accuracy of NCA substantially. Such information on TF-gene interactions is obtained mainly through ChIP-CHIP or ChIP-SEQ experiments that allow the detection of binding patterns of TFs with DNA sequences. In fact, a lot of the confirmed interactions between TFs and genes listed on AGRIS are derived from such papers investigating binding targets for particular TFs .
Another limitation of NCA is its inability to model feedback and feedforward regulations between TFs. TF-GRNs are cascades of TFs regulating genes where the product of many genes are TFs that regulate downstream genes. However, for NCA, if a TF is included as a regulator in a network, the gene corresponding to it cannot be included in the network. As a result, NCA cannot determine how strongly other TFs influence the expression of the corresponding gene. In our original network, AG was included as a TF and also present as a gene regulated by LFY, AG, SEP3, AP2, WUS and BLR. We had to remove the AG gene during the NCA because of the presence of AG as a regulatory TF. This limits the application of NCA to non TF target genes in many instances.
Additionally, the NCA decomposition suffers from some variability in estimating CS and TF activity from gene expression data. This is because the NCA decomposition is unique to a scaling factor which can be different for each TF and vary during different data decomposition of the same set of gene expression values and initial connectivity matrix. NCA uses a two-step least squares approach to minimize the difference between experimental and NCA reconstructed gene expression data. As a result, based on the scaling factor chosen, the same gene expression data and initial connectivity matrix could give slightly differing TF activities and CS. In addition, the decomposition process might introduce some variability in estimating TF activities and CS. For the NCA decomposition of the floral TF-GRN used in this study, we found differences in TF activities and CS during repeat runs (Additional file 3). For this network, the LFY TF shows very little variability across the different runs while the other TFs have greater degree of variability. Thus, while the TF activity and CS obtained from NCA decomposition provide quantitative measures for the underlying network, they should be treated not as absolute but relative parameters.
Another drawback that all approaches for modeling gene expression of eukaryotic organisms suffer from, is the inability to include all the factors that regulate gene expression . Most of the current modeling approaches depict gene expression to result from the effect of some of these factors alone, which is not the case . For example, microRNAs play a very important role in gene regulation at the post-transcriptional level similar to the TF regulation at the transcriptional level [54–56]. In humans, microRNAs have been found to use two modes for gene regulation – the first mode is rapid and modulated by homoclusters; the second is delayed and mediated by heteroclusters of microRNAs. Of the two, heteroclusters have been found to indirectly influence gene regulation in tandem with TFs . In addition to microRNAs, other factors including chromatin structure and nucleosome sliding would affect gene expression especially in eukaryotes . Consequently, an accurate model for depicting gene regulation in eukaryotes would have to include all these interactions to capture the true picture of genetic regulation.
Despite these limitations, NCA can provide very interesting hypotheses and insights about regulatory signals in a TF-GRN. Previous applications have shown its utility in understanding microbial systems whose regulatory networks are well characterized, and mammalian sytems to some extent. Plants and eukaryotes operate more complex regulatory mechanisms. Additionally, complicated post-translational modifications can alter the activity of a TF compared to its mRNA transcript level. Consequently, the application of NCA to plant systems would provide interesting insights about these. Hence, there is a need for applying significant efforts in obtaining information about interactions between TFs and genes in plants for constructing TF-GRNs. Such information coupled with NCA would allow the determination of underlying properties of the system and establish paradigms for predicting cellular behavior.
In this work, we applied constructed a plant TF-GRN important in flower development using regulatory information from the AGRIS database. The initial network consisting of 10 TFs and 57 genes was found to be NCA-compliant for 7 TFs and 55 genes. We applied NCA to the reduced network to obtain CS between TF-gene pairs and TF activities. The CS showed strong connectivity between certain TF-gene pairs including LFY → MYB17, LFY → TLP8, AP2 → HLH1, AP2 → RD20, AGL15 → AGL22, AGL15 → RAV2, HY5 → HLH1 and HY5 → RD20, among others. For some of the co-regulated genes, we were able to determine the extent of transcriptional control of different TFs on a target gene using the CS. Additionally, we were able to determine TF activities for all TFs. Good agreement was seen for the changes in TF activities for multiple TFs and their corresponding gene expression levels. However, for some of the TFs (AP2, SEP3 and AGL15), the change in TF activities did not match with changes in gene expression levels. There could be multiple reasons for this discrepancy including post translation modifications which significantly alter the activity of a TF; noisy data or the small size of the network among others.
Our study is the first application of NCA to a plant TF-GRN and demonstrates the power of NCA for determining nontrivial information about a network based solely on gene expression data and underlying network connectivity. NCA has been widely used to decipher interesting insights about microbial TF-GRNs. However, since NCA relies on underlying network connectivity, incomplete information about the network hinders the accuracy of NCA. Plant TF-GRNs are poorly documented with sparse data about specific sets of TFs and processes. As more information about TF-GRNs is uncovered in plants, similar analysis using NCA would provide profound insights regarding the role of TFs in various cellular processes.
TF-gene network reconstruction
We obtained TF-gene connectivity information from AGRIS (http://arabidopsis.med.ohio-state.edu) . For the GRN analysis, we selected 10 TFs known to be important in floral development and listed in AGRIS. We selected 57 genes that were documented in AGRIS to be the targets of these TFs (Additional file 1, Sheet: AGRIS TF-gene verification). We constructed an initial connectivity matrix to map the TF-gene interactions documented in AGRIS (Additional file 1, Sheet: Initial connectivity matrix). Entries in this matrix were 1 (indicating a documented activation interaction), –1 (indicating a documented repression interaction) or 0 (indicating no documented interaction). Documented TF-gene interactions for which the type of interaction (activation or repression) were not known were assigned an entry of 1 (highlighted cells).
Gene expression data
We used the Botany Array Resource (http://www.bar.utoronto.ca)  for obtaining gene expression data pertinent to the TFs and genes in our network during floral development. This database provided gene expression data from the study by Yadav et al. that provided expression levels of the genes of interest across four SAM cell types. The original and log transformed gene expression data are summarized in Additional file 1 (Sheet: Original microarray data, and Sheet: Log transformed microarray data, respectively).
We used the NCA toolbox (http://www.seas.ucla.edu/~liaoj/downloads.html) [26, 30] in conjunction with the initial TF-gene connectivity matrix (Additional file 1, Sheet: Initial connectivity matrix) for decomposing the gene expression data. We independently analyzed the gene expression dataset corresponding to each biological replicate of each cell line. On completion, NCA provided TF activities for each replicate of each cell line (Additional file 2, Sheet: TFA and mRNA) as well as TF-gene CS common to all cell lines (Additional file 2, Sheet: Connectivity strengths).
The authors wish to thank Sam Prager (Department of Chemical and Biomolecular Engineering, University of Maryland) for his assistance with the data analysis. This work was funded by the U.S. National Science Foundation (award number IOS-0922650).
- van Someren E, Wessels L, Backer E, Reinders M: Genetic network modeling. Pharmacogenomics. 2002, 3: 507-525. 10.1517/146224184.108.40.2067.PubMedView ArticleGoogle Scholar
- Karlebach G, Shamir R: Modelling and analysis of gene regulatory networks. Nat Rev Mol Cell Biol. 2008, 9: 770-780. 10.1038/nrm2503.PubMedView ArticleGoogle Scholar
- Markowetz F, Spang R: Inferring cellular networks – a review. BMC Bioinformatics. 2007, 8 (Suppl 6): S5-10.1186/1471-2105-8-S6-S5.PubMedPubMed CentralView ArticleGoogle Scholar
- Moreno-Risueno MA, Busch W, Benfey PN: Omics meet networks – using systems approaches to infer regulatory networks in plants. Curr Opin Plant Biol. 2010, 13: 126-131. 10.1016/j.pbi.2009.11.005.PubMedPubMed CentralView ArticleGoogle Scholar
- Schlitt T, Brazma A: Current approaches to gene regulatory network modelling. BMC Bioinformatics. 2007, 8 (Suppl 6): S9-10.1186/1471-2105-8-S6-S9.PubMedPubMed CentralView ArticleGoogle Scholar
- De Jong H: Modeling and simulation of genetic regulatory systems: a literature review. J Comput Biol. 2002, 9: 67-103. 10.1089/10665270252833208.PubMedView ArticleGoogle Scholar
- Bezerianos A, Maraziotis IA: Computational models reconstruct gene regulatory networks. Mol Biosyst. 2008, 4: 993-1000. 10.1039/b800446n.PubMedView ArticleGoogle Scholar
- Li Z, Shaw SM, Yedwabnick MJ, Chan C: Using a state-space model with hidden variables to infer transcription factor activities. Bioinformatics. 2006, 22: 747-754. 10.1093/bioinformatics/btk034.PubMedView ArticleGoogle Scholar
- Long TA, Brady SM, Benfey PN: Systems approaches to identifying gene regulatory networks in plants. Annu Rev Cell Dev Biol. 2008, 24: 81-103. 10.1146/annurev.cellbio.24.110707.175408.PubMedPubMed CentralView ArticleGoogle Scholar
- Mussel C, Hopfensitz M, Kestler HA: BoolNet–an R package for generation, reconstruction and analysis of Boolean networks. Bioinformatics. 2010, 26: 1378-1380. 10.1093/bioinformatics/btq124.PubMedView ArticleGoogle Scholar
- Baldan P, Cocco N, Marin A, Simeoni M: Petri nets for modelling metabolic pathways: a survey. Nat Comput. 9: 955-989.
- Babu MM, Luscombe NM, Aravind L, Gerstein M, Teichmann SA: Structure and evolution of transcriptional regulatory networks. Curr Opin Struct Biol. 2004, 14: 283-291. 10.1016/j.sbi.2004.05.004.PubMedView ArticleGoogle Scholar
- Gama-Castro S, Salgado H, Peralta-Gil M, Santos-Zavaleta A, Muniz-Rascado L, Solano-Lira H, Jimenez-Jacinto V, Weiss V, Garcia-Sotelo JS, Lopez-Fuentes A, Porron-Sotelo L, Alquicira-Hernandez S, Medina-Rivera A, Martinez-Flores I, Alquicira-Hernandez K, Martinez-Adame R, Bonavides-Martinez C, Miranda-Rios J, Huerta AM, Mendoza-Vargas A, Collado-Torres L, Taboada B, Vega-Alvarado L, Olvera M, Olvera L, Grande R, Morett E, Collado-Vides J: RegulonDB version 7.0: transcriptional regulation of Escherichia coli K-12 integrated within genetic sensory response units (Gensor Units). Nucleic Acids Res. 2011, 39 (Database issue): 98-105.View ArticleGoogle Scholar
- Makita Y, Nakao M, Ogasawara N, Nakai K: DBTBS: database of transcriptional regulation in Bacillus subtilis and its contribution to comparative genomics. Nucleic Acids Res. 2004, 32 (Database issue): D75-D77.PubMedPubMed CentralView ArticleGoogle Scholar
- Novichkov PS, Laikova ON, Novichkova ES, Gelfand MS, Arkin AP, Dubchak I, Rodionov DA: RegPrecise: a database of curated genomic inferences of transcriptional regulatory interactions in prokaryotes. Nucleic Acids Res. 2010, 38 (suppl 1): D111-D118.PubMedPubMed CentralView ArticleGoogle Scholar
- Abdulrehman D, Monteiro PT, Teixeira MC, Mira NP, Lourenço AB, dos Santos SC, Cabrito TR, Francisco AP, Madeira SC, Aires RS, Oliveira AL, Sá-Correia I, Freitas AT: YEASTRACT: providing a programmatic access to curated transcriptional regulatory associations in Saccharomyces cerevisiae through a web services interface. Nucleic Acids Res. 2011, 39 (suppl 1): D136-D140.PubMedPubMed CentralView ArticleGoogle Scholar
- Middleton AM, Farcot E, Owen MR, Vernoux T: Modeling regulatory networks to understand plant development: small is beautiful. Plant Cell Online. 2012, 24: 3876-3891. 10.1105/tpc.112.101840.View ArticleGoogle Scholar
- Alvarez-Buylla ER, Benitez M, Davila EB, Chaos A, Espinosa-Soto C, Padilla-Longoria P: Gene regulatory network models for plant development. Curr Opin Plant Biol. 2007, 10: 83-91. 10.1016/j.pbi.2006.11.008.PubMedView ArticleGoogle Scholar
- Palaniswamy SK, James S, Sun H, Lamb RS, Davuluri RV, Grotewold E: AGRIS and AtRegNet. A platform to link cis-regulatory elements and transcription factors into regulatory networks. Plant Physiol. 2006, 140: 818-829. 10.1104/pp.105.072280.PubMedPubMed CentralView ArticleGoogle Scholar
- Bintu L, Buchler NE, Garcia HG, Gerland U, Hwa T, Kondev J, Phillips R: Transcriptional regulation by the numbers: models. Curr Opin Genet Dev. 2005, 15: 116-124. 10.1016/j.gde.2005.02.007.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhao Y, Granas D, Stormo GD: Inferring binding energies from selected binding sites. PLoS Comput Biol. 2009, 5: e1000590-10.1371/journal.pcbi.1000590.PubMedPubMed CentralView ArticleGoogle Scholar
- Turner D, Kim R, Guo J: TFinDit: transcription factor-DNA interaction data depository. BMC Bioinformatics. 2012, 13: 220-10.1186/1471-2105-13-220.PubMedPubMed CentralView ArticleGoogle Scholar
- He X, Samee MAH, Blatti C, Sinha S: Thermodynamics-based models of transcriptional regulation by enhancers: the roles of synergistic activation. Cooperative binding and short-range repression. PLoS Comput Biol. 2010, 6: e1000935-10.1371/journal.pcbi.1000935.PubMedPubMed CentralView ArticleGoogle Scholar
- Geertz M, Maerkl SJ: Experimental strategies for studying transcription factor-DNA binding specificities. Briefings Funct Genomics. 2010, 9: 362-373. 10.1093/bfgp/elq023.View ArticleGoogle Scholar
- Prabakaran P, An J, Gromiha MM, Selvaraj S, Uedaira H, Kono H, Sarai A: Thermodynamic database for protein–nucleic acid interactions (ProNIT). Bioinformatics. 2001, 17: 1027-1034. 10.1093/bioinformatics/17.11.1027.PubMedView ArticleGoogle Scholar
- Liao JC, Boscolo R, Yang Y-L, Tran LM, Sabatti C, Roychowdhury VP: Network component analysis: reconstruction of regulatory signals in biological systems. Proc Natl Acad Sci U S A. 2003, 100: 15522-15527. 10.1073/pnas.2136632100.PubMedPubMed CentralView ArticleGoogle Scholar
- Kao KC, Yang Y-L, Boscolo R, Sabatti C, Roychowdhury V, Liao JC: Transcriptome-based determination of multiple transcription regulator activities in Escherichia coli by using network component analysis. Proc Natl Acad Sci. 2004, 101: 641-646. 10.1073/pnas.0305287101.PubMedPubMed CentralView ArticleGoogle Scholar
- Kao KC, Tran LM, Liao JC: A global regulatory role of gluconeogenic genes in escherichia coli revealed by transcriptome network analysis. J Biol Chem. 2005, 280: 36079-36087. 10.1074/jbc.M508202200.PubMedView ArticleGoogle Scholar
- Roven C, Bussemaker HJ: REDUCE: an online tool for inferring cis-regulatory elements and transcriptional module activities from microarray data. Nucleic Acids Res. 2003, 31: 3487-3490. 10.1093/nar/gkg630.PubMedPubMed CentralView ArticleGoogle Scholar
- Galbraith SJ, Tran LM, Liao JC: Transcriptome network component analysis with limited microarray data. Bioinformatics. 2006, 22: 1886-1894. 10.1093/bioinformatics/btl279.PubMedView ArticleGoogle Scholar
- Hyduke DR, Jarboe LR, Tran LM, Chou KJY, Liao JC: Integrated network analysis identifies nitric oxide response networks and dihydroxyacid dehydratase as a crucial target in Escherichia coli. Proc Natl Acad Sci. 2007, 104: 8484-8489. 10.1073/pnas.0610888104.PubMedPubMed CentralView ArticleGoogle Scholar
- Jarboe LR, Hyduke DR, Tran LM, Chou KJY, Liao JC: Determination of the Escherichia coli S-Nitrosoglutathione response network using integrated biochemical and systems analysis. J Biol Chem. 2008, 283: 5148-5157.PubMedView ArticleGoogle Scholar
- Brynildsen MP, Liao JC: An integrated network approach identifies the isobutanol response network of Escherichia coli. Mol Syst Biol. 2009, 5: 277-PubMedPubMed CentralView ArticleGoogle Scholar
- Buescher JM, Liebermeister W, Jules M, Uhr M, Muntel J, Botella E, Hessling B, Kleijn RJ, Chat LL, Lecointe F, Mäder U, Nicolas P, Piersma S, Rügheimer F, Becher D, Bessieres P, Bidnenko E, Denham EL, Dervyn E, Devine KM, Doherty G, Drulhe S, Felicori L, Fogg MJ, Goelzer A, Hansen A, Harwood CR, Hecker M, Hubner S, Hultschig C, et al: Global network reorganization during dynamic adaptations of Bacillus subtilis metabolism. Science. 2012, 335: 1099-1103. 10.1126/science.1206871.PubMedView ArticleGoogle Scholar
- Sriram G, Parr LS, Rahib L, Liao JC, Dipple KM: Moonlighting function of glycerol kinase causes systems-level changes in rat hepatoma cells. Metab Eng. 2010, 12: 332-340. 10.1016/j.ymben.2010.04.001.PubMedPubMed CentralView ArticleGoogle Scholar
- Sriram G, Martinez JA, McCabe ERB, Liao JC, Dipple KM: Single-gene disorders: what role could moonlighting enzymes play?. Am J Hum Genet. 2005, 76: 911-924. 10.1086/430799.PubMedPubMed CentralView ArticleGoogle Scholar
- Tran LM, Chang C-J, Plaisier S, Wu S, Dang J, Mischel PS, Liao JC, Graeber TG, Wu H: Determining PTEN functional status by network component deduced transcription factor activities. PLoS One. 2012, 7: e31053-10.1371/journal.pone.0031053.PubMedPubMed CentralView ArticleGoogle Scholar
- Yadav RK, Girke T, Pasala S, Xie M, Reddy GV: Gene expression map of the Arabidopsis shoot apical meristem stem cell niche. Proc Natl Acad Sci. 2009, 106: 4941-4946. 10.1073/pnas.0900843106.PubMedPubMed CentralView ArticleGoogle Scholar
- Toufighi K, Brady SM, Austin R, Ly E, Provart NJ: The botany array resource: e-northerns, expression angling, and promoter analyses. Plant J. 2005, 43: 153-163. 10.1111/j.1365-313X.2005.02437.x.PubMedView ArticleGoogle Scholar
- Chang C, Ding Z, Hung YS, Fung PCW: Fast network component analysis (FastNCA) for gene regulatory network reconstruction from microarray data. Bioinformatics. 2008, 24: 1349-1358. 10.1093/bioinformatics/btn131.PubMedView ArticleGoogle Scholar
- Weigel D, Meyerowitz EM: The ABCs of floral homeotic genes. Cell. 1994, 78: 203-209. 10.1016/0092-8674(94)90291-7.PubMedView ArticleGoogle Scholar
- Causier B, Schwarz-Sommer Z, Davies B: Floral organ identity: 20 years of ABCs. Semin Cell Dev Biol. 2010, 21: 73-79. 10.1016/j.semcdb.2009.10.005.PubMedView ArticleGoogle Scholar
- Theißen G: Development of floral organ identity: stories from the MADS house. Curr Opin Plant Biol. 2001, 4: 75-85. 10.1016/S1369-5266(00)00139-4.PubMedView ArticleGoogle Scholar
- Siriwardana NS, Lamb RS: The poetry of reproduction: the role of LEAFY in Arabidopsis thaliana flower formation. Int J Dev Biol. 2012, 56: 207-221. 10.1387/ijdb.113450ns.PubMedView ArticleGoogle Scholar
- William DA, Su Y, Smith MR, Lu M, Baldwin DA, Wagner D: Genomic identification of direct target genes of LEAFY. Proc Natl Acad Sci U S A. 2004, 101: 1775-1780. 10.1073/pnas.0307842100.PubMedPubMed CentralView ArticleGoogle Scholar
- Weigel D, Alvarez J, Smyth DR, Yanofsky MF, Meyerowitz EM: LEAFY controls floral meristem identity in Arabidopsis. Cell. 1992, 69: 843-859. 10.1016/0092-8674(92)90295-N.PubMedView ArticleGoogle Scholar
- Swarbreck D, Wilks C, Lamesch P, Berardini TZ, Garcia-Hernandez M, Foerster H, Li D, Meyer T, Muller R, Ploetz L, Radenbaugh A, Singh S, Swing V, Tissier C, Zhang P, Huala E: The Arabidopsis Information Resource (TAIR): gene structure and function annotation. Nucleic Acids Res. 2008, 36 (suppl 1): D1009-D1014.PubMedPubMed CentralGoogle Scholar
- Adamczyk BJ, Lehti-Shiu MD, Fernandez DE: The MADS domain factors AGL15 and AGL18 act redundantly as repressors of the floral transition in Arabidopsis. Plant J Cell Mol Biol. 2007, 50: 1007-1019. 10.1111/j.1365-313X.2007.03105.x.View ArticleGoogle Scholar
- Gregis V, Andrés F, Sessa A, Guerra RF, Simonini S, Mateos JL, Torti S, Zambelli F, Prazzoli GM, Bjerkan KN, Grini PE, Pavesi G, Colombo L, Coupland G, Kater MM: Identification of pathways directly regulated by Short Vegetative Phase during vegetative and reproductive development in Arabidopsis. Genome Biol. 2013, 14: R56-10.1186/gb-2013-14-6-r56.PubMedPubMed CentralView ArticleGoogle Scholar
- Gregis V, Sessa A, Colombo L, Kater MM: AGAMOUS-LIKE24 and SHORT VEGETATIVE PHASE determine floral meristem identity in Arabidopsis. Plant J. 2008, 56: 891-902. 10.1111/j.1365-313X.2008.03648.x.PubMedView ArticleGoogle Scholar
- Pastore JJ, Limpuangthip A, Yamaguchi N, Wu M-F, Sang Y, Han S-K, Malaspina L, Chavdaroff N, Yamaguchi A, Wagner D: LATE MERISTEM IDENTITY2 acts together with LEAFY to activate APETALA1. Development. 2011, 138: 3189-3198. 10.1242/dev.063073.PubMedPubMed CentralView ArticleGoogle Scholar
- Riechmann JL, Ratcliffe OJ: A genomic perspective on plant transcription factors. Curr Opin Plant Biol. 2000, 3: 423-434. 10.1016/S1369-5266(00)00107-2.PubMedView ArticleGoogle Scholar
- Wilczynski B, Furlong EEM: Challenges for modeling global gene regulatory networks during development: Insights from Drosophila. Dev Biol. 2010, 340: 161-169. 10.1016/j.ydbio.2009.10.032.PubMedView ArticleGoogle Scholar
- Wang J, Haubrock M, Cao K-M, Hua X, Zhang C-Y, Wingender E, Li J: Regulatory coordination of clustered microRNAs based on microRNA-transcription factor regulatory network. BMC Syst Biol. 2011, 5: 199-10.1186/1752-0509-5-199.PubMedPubMed CentralView ArticleGoogle Scholar
- Lin C-C, Chen Y-J, Chen C-Y, Oyang Y-J, Juan H-F, Huang H-C: Crosstalk between transcription factors and microRNAs in human protein interaction network. BMC Syst Biol. 2012, 6: 18-10.1186/1752-0509-6-18.PubMedPubMed CentralView ArticleGoogle Scholar
- Croft L, Szklarczyk D, Jensen LJ, Gorodkin J: Multiple independent analyses reveal only transcription factors as an enriched functional class associated with microRNAs. BMC Syst Biol. 2012, 6: 90-10.1186/1752-0509-6-90.PubMedPubMed CentralView 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.