- Research article
- Open Access

# Inferring yeast cell cycle regulators and interactions using transcription factor activities

- Young-Lyeol Yang
^{1}, - Jason Suen
^{1}, - Mark P Brynildsen
^{1}, - Simon J Galbraith
^{1}and - James C Liao
^{1}Email author

**6**:90

https://doi.org/10.1186/1471-2164-6-90

© Yang et al; licensee BioMed Central Ltd. 2005

**Received:**13 January 2005**Accepted:**10 June 2005**Published:**10 June 2005

## Abstract

### Background

Since transcription factors are often regulated at the post-transcriptional level, their activities, rather than expression levels may provide valuable information for investigating functions and their interactions. The recently developed Network Component Analysis (NCA) and its generalized form (gNCA) provide a robust framework for deducing the transcription factor activities (TFAs) from various types of DNA microarray data and transcription factor-gene connectivity. The goal of this work is to demonstrate the utility of TFAs in inferring transcription factor functions and interactions in *Saccharomyces cerevisiae* cell cycle regulation.

### Results

Using gNCA, we determined 74 TFAs from both wild type and *fkh1 fkh2* deletion mutant microarray data encompassing 1529 ORFs. We hypothesized that transcription factors participating in the cell cycle regulation exhibit cyclic activity profiles. This hypothesis was supported by the TFA profiles of known cell cycle factors and was used as a basis to uncover other potential cell cycle factors. By combining the results from both cluster analysis and periodicity analysis, we recovered nearly 90% of the known cell cycle regulators, and identified 5 putative cell cycle-related transcription factors (Dal81, Hap2, Hir2, Mss11, and Rlm1). In addition, by analyzing expression data from transcription factor knockout strains, we determined 3 verified (Ace2, Ndd1, and Swi5) and 4 putative interaction partners (Cha4, Hap2, Fhl1, and Rts2) of the forkhead transcription factors. Sensitivity of TFAs to connectivity errors was determined to provide confidence level of these predictions.

### Conclusion

By subjecting TFA profiles to analyses based upon physiological signatures we were able to identify cell cycle related transcription factors consistent with current literature, transcription factors with potential cell cycle dependent roles, and interactions between transcription factors.

## Keywords

- Transcription Factor Activity
- Forkhead Transcription Factor
- Connectivity Information
- Connectivity Error
- Putative Cell

## Background

Transcription factor activities (TFAs) rather than levels of transcription factor expression mediate transcriptional regulations. Various post-transcriptional and post-translational modifications abolish significant correlations between TFAs and the level of transcription factor expression. Therefore, a strategy to determine the physiological functions and interactions of transcription factors based on the analysis of TFA profiles would be more fundamentally sound than one based on transcript levels *per se*. However, owing to post-translational modifications, TFAs are difficult to measure experimentally and therefore, many analyses infer TFAs by computational analysis of target gene expression levels either from single regulatory factors alone or in combination [1–10]. Among these, Network Component Analysis (NCA) [3] and generalized NCA (gNCA) [2] provide a robust framework for deducing TFAs based on DNA microarray data, promoter connectivity, and genetic regulatory constraints imposed by regulatory knock-out experiments. Using TFA profiles deduced by NCA and gNCA, we are now in a position to assign transcription factor functions and reconstruct functional interactions between transcription factors.

The goal of this work is to demonstrate how we analyze TFA profiles to determine physiologically relevant characteristics of transcription factors. Specifically, by analyzing TFA profiles of the *S. cerevisiae* cell cycle we identify cell cycle regulators and putative interaction partners of two forkhead transcription factors (Fkh1 and Fkh2) that are responsible for expression of a gene cluster within the M/G1 interval and involved in mitotic exit [11].

Once TFA profiles are determined, the transcription factor regulators and interactions are deduced based on the following hypotheses: 1) Transcription factors with similar activity patterns function together; 2) Transcription factors involved in cell cycle regulation exhibit oscillatory activity patterns, 3) Activities of transcription factors which functionally interact with each other will be disturbed if one of the interacting partners is deleted. These hypotheses are in spirit similar to the traditional analysis of gene expression. In contrast to traditional analysis of gene expression, our analyses are constructed and applied directly to TFA profiles.

## Results

### Overall strategy

*fkh1 fkh2*mutant strain under the same experimental conditions were statistically analyzed to deduce transcription factors whose profiles were significantly perturbed, thus identifying putative functional interaction partners of the fork-head transcription factors.

### Application of gNCA to the combined wild-type and mutant data set

We used gNCA to analyze both the wild-type [13] and the *fkh1 fkh2* mutant data [11]. The wild-type data and the TF-knockout data were organized in different columns of the **E** matrix (expression data matrix where *E*_{
ij
}corresponds to the log of gene expression ratio of gene *i* evaluated at experiment *j*, see Methods). The elements in **P** (TFA matrix where *P*_{
ij
}corresponds to the log of relative transcription factor activity of TF *i* at experiment *j*, see Methods) corresponding to the specific knockout strain were kept at zero, in addition to the zero constraints placed in **A** (control strength matrix where *A*_{
ij
}denotes the control strength of transcription factor *j* on gene *i*, see Methods) according to the transcription factor-gene connectivity information.

*L*) in each sub-network be less than that of data points (

*M*) in the data matrix. In the combined cell cycle microarray data set,

*M*= 69, and thus the number of transcription factor in gNCA in each analysis cannot exceed 69. Thus, subsets of the 104 transcription factors included in the genome-wide location analysis [14] were selected to form sub-networks for gNCA. Overlapping random sub-networks were generated and screened using the gNCA criteria. Among them, four that satisfy the gNCA criteria were selected in order to show that using multiple sub-networks makes it possible to determine TFAs more than the number of data points (M = 69). Each sub-network contained 40 transcription factors, but together, a total of 74 transcription factors were analyzed. The transcription factors involved in each of the 4 sub-networks are shown in Figure 2. Only 16 transcription factors among 74 transcription factors were fully overlapped among all 4 sub-networks.

*cdc15*mutation). The dynamics of overlapping TFA profiles between the four different sub-networks are very similar. This indicates that the deduced TFAs were insensitive to the sub-network selection. In addition, combining the

*fkh1 fkh2*mutant data with the wild-type data yields similar TFAs predicted from wild-type data alone [3], suggesting the consistency between the two data sets.

It is recognized that most transcription factors are regulated in post-transcriptional processes. Thus, their activity may deviate significantly from the expression level as noted previously [3]. For 27 transcription factors with both gene expression level and TFA available, the correlation coefficients ranged from - 0.5 to 0.6. The wide range difference in the correlation coefficients confirms that TFAs and expression levels exhibit very different dynamic profiles and therefore TFAs cannot be substituted by expression levels in analysis.

### Sensitivity Analysis

Sensitivity of TFA profiles to errors in connectivity from ChIP-chip assay. Connectivity errors were simulated by randomly inserting and deleting connections in the connectivity graph generated from the ChIP-chip assay. gNCA was performed on the perturbed connectivity, and the resultant TFAs were compared to those from un-perturbed connectivity via the correlation coefficient. This procedure was performed on each of the 4 subnetworks 100 times. The average correlation coefficients for transcription factors in the 4 subnetworks are presented here. A low correlation coefficient (<0.5) suggests sensitivity to network connectivity error.

Subnetwork 1 | Mean Correlation Coefficient | Subnetwork 2 | Mean Correlation Coefficient | Subnetwork 3 | Mean Correlation Coefficient | Subnetwork 4 | Mean Correlation Coefficient |
---|---|---|---|---|---|---|---|

ABF1 | 0.941 | ABF1 | 0.951 | ABF1 | 0.94 | ACE2 | 0.912 |

ACE2 | 0.953 | ACE2 | 0.857 | ACE2 | 0.897 | ARO80 | 0.657 |

ARO80 | 0.645 | ARG81 | 0.286 | CHA4 | 0.362 | CIN5 | 0.601 |

CAD1 | 0.511 | CAD1 | 0.659 | CIN5 | 0.59 | DIG1 | 0.718 |

CBF1 | 0.401 | CHA4 | 0.448 | CRZ1 | 0.313 | FHL1 | 0.992 |

CIN5 | 0.651 | DAL81 | 0.581 | FHL1 | 0.993 | FKH2 | 0.669 |

DIG1 | 0.686 | FHL1 | 0.99 | FKH2 | 0.657 | GAT3 | 0.493 |

FHL1 | 0.992 | FKH2 | 0.635 | GAL4 | 0.492 | GCN4 | 0.601 |

FKH1 | 0.614 | GAT3 | 0.457 | GAT3 | 0.511 | GCR1 | 0.321 |

FKH2 | 0.535 | GCR2 | 0.425 | GCR1 | 0.306 | GLN3 | 0.402 |

GAT3 | 0.446 | GRF10(Pho2) | 0.565 | GCR2 | 0.446 | GTS1 | 0.227 |

HAP4 | 0.706 | HIR1 | 0.732 | HAP4 | 0.769 | HAP2 | 0.375 |

HIR1 | 0.472 | HIR2 | 0.357 | HIR1 | 0.465 | HAP3 | 0.47 |

IXR1 | 0.511 | IME4 | 0.357 | HSF1 | 0.981 | HAP4 | 0.856 |

MBP1 | 0.946 | INO4 | 0.419 | IME4 | 0.442 | HIR1 | 0.32 |

MCM1 | 0.558 | LEU3 | 0.575 | MBP1 | 0.979 | IXR1 | 0.475 |

MET31 | 0.436 | MAC1 | 0.723 | MCM1 | 0.572 | LEU3 | 0.518 |

MET4 | 0.381 | MBP1 | 0.985 | MET31 | 0.435 | MBP1 | 0.951 |

MIG1 | 0.336 | MCM1 | 0.701 | MIG1 | 0.354 | MCM1 | 0.705 |

MSN4 | 0.444 | MET31 | 0.507 | NDD1 | 0.965 | MIG1 | 0.33 |

MTH1 | 0.463 | MOT3 | 0.413 | NRG1 | 0.421 | MSS11 | 0.404 |

NDD1 | 0.958 | NDD1 | 0.971 | PDR1 | 0.366 | NDD1 | 0.962 |

NRG1 | 0.495 | PHD1 | 0.541 | PHD1 | 0.583 | NRG1 | 0.563 |

PDR1 | 0.334 | PHO4 | 0.599 | RAP1 | 0.917 | PDR1 | 0.372 |

PHD1 | 0.569 | RAP1 | 0.856 | RCS1 | 0.48 | PHD1 | 0.64 |

RAP1 | 0.917 | RME1 | 0.407 | REB1 | 0.576 | RAP1 | 0.897 |

REB1 | 0.643 | RTG1 | 0.413 | RME1 | 0.415 | RCS1 | 0.424 |

RFX1 | 0.493 | RTS2 | 0.302 | RPH1 | 0.294 | RLM1 | 0.553 |

RLM1 | 0.519 | SFP1 | 0.564 | SFP1 | 0.582 | RME1 | 0.389 |

RME1 | 0.367 | SKN7 | 0.536 | SKN7 | 0.669 | SKN7 | 0.609 |

SKN7 | 0.545 | STE12 | 0.518 | SMP1 | 0.619 | SMP1 | 0.642 |

SMP1 | 0.614 | STP2 | 0.322 | SWI4 | 0.553 | STB1 | 0.576 |

STE12 | 0.79 | SWI4 | 0.592 | SWI5 | 0.987 | STE12 | 0.738 |

SWI4 | 0.747 | SWI5 | 0.963 | SWI6 | 0.61 | SUM1 | 0.823 |

SWI5 | 0.991 | SWI6 | 0.648 | UGA3 | 0.392 | SWI4 | 0.843 |

SWI6 | 0.606 | UGA3 | 0.383 | YAP5 | 0.512 | SWI5 | 0.991 |

YAP5 | 0.863 | YAP1 | 0.565 | YAP6 | 0.404 | SWI6 | 0.589 |

YAP6 | 0.56 | YAP5 | 0.657 | YFL044C | 0.578 | YAP5 | 0.572 |

YFL044C | 0.583 | YFL044C | 0.546 | ZAP1 | 0.442 | YAP6 | 0.595 |

YJL206C | 0.396 | ZAP1 | 0.527 | ZMS1 | 0.295 | YJL206C | 0.455 |

### Determination of cell cycle-dependent regulators based on TFA dynamics

*et al*. [14] categorized Dal81 and Mss11 as metabolism regulators, Hir2 as a DNA/RNA/Protein biosynthesis regulator, and Rlm1 as environmental response regulators. These results suggested that transcription factors related to many other cellular processes may be involved in or dependent on cell cycle regulation to coordinate cellular processes [14].

### Functional interaction of FKH1 FKH2 with other transcription factors

Comparing the TFAs derived from the wild-type strain and the *fkh1 fkh2* mutant under the same experimental conditions allows the determination of functional interactions between Fkh1 Fkh2 with other transcription factors. If such interaction exists, the activities of the interaction partners will change in the *fkh1 fkh2* mutant. In general the *fkh1 fkh2* mutant showed reduced oscillation amplitude as compared to the wild type (Figure 3). In particular, Ace2, Swi5, and Ndd1 showed significantly reduced oscillation in the *fkh1 fkh2* mutant, suggesting potential interactions between these transcription factors and the two fork-head transcription factors. Indeed, it was found that Fkh2 forms a complex with Mcm1 and Ndd1 [16, 17] which controls G2/M genes and also binds the promoters of Ace2 and Swi5 to activate genes at the following M/G1 phase. The reduced amplitude in TFAs of Ace2 and Swi5 in the *fkh1 fkh2* mutant support that there is a cascade interaction between the two forkhead transcription factors and Ace2 and Swi5 [16].

TFA profiles of the wild-type strain (Figure 3) show that Ace2 and Swi5 were most active around M/G1 phase while Fkh1, Fkh2, Mcm1, and Ndd1 were most active at G2/M, G2, M/G1, and G2/M phase, respectively. These results are well consistent with previous results that Mcm1, Ace2, and Swi5 are M/G1 transcriptional regulators while the complex of Mcm1, Fkh2 and Ndd1 is a G2/M activator. Therefore, this result shows that TFA in a regulator knockout strain can be used to determine functional interactions.

By utilizing the Pearson correlation coefficient and deviation coefficient defined in the Methods section, Figure 8 shows that most transcription factors were aggregated in one cluster. Outside of this cluster are transcription factors which potentially were affected by the *fkh1 fkh2* mutation and thus are functional interaction partners of the forkhead transcription factors. These interaction partners include Ace2, Ndd1, and Swi5, which are known to interact with Fkh2, and Hap2, Rts2, Cha4, and Fhl1, which were not known to interact with the forkhead transcription factors. In particular, Hap2 was also identified in the above analysis as a putative cell cycle-related regulator, supporting their functional interaction with the forkhead transcription factors. Cha4 and Hap2 were classified as metabolism-related factors, Fhl1 DNA/RNA/Protein synthesis, and Rts2 cell cycle and data processing [14]. Their modes of interaction with the forkhead transcription factors remain unknown.

## Discussion

Transcriptional regulators are commonly modified at the post-transcriptional level, and consequently their biological activities do not correlate significantly with expression levels. Previous work infer TFAs from expression levels of genes regulated by single factors or in combination [2, 5–8, 10]. In general, the major difference between previous work and that of NCA is that the former require an explicit quantification of the control strengths (the A matrix in NCA) *a priori*. Bussemaker et al. [7] defined the control strength as the motif copy number in corresponding promoters and found that there was no statistical benefit to model expression with more than single factors and thus deduced single TFAs. Wang et al. [6] considered an expression-weighted motif to find potential target genes of single transcription factors. Gao et al. [18] used ChIP-chip log occupancy ratios as a surrogate for transcription factor binding affinity. In contrast, NCA explicitly models combinatorial regulation of gene expression, and allows both the control strengths and the TFAs to be deduced simultaneously with given network connectivity. In this approach, the lack of connectivity in specific pairs of transcription factor and promoter is used to provide constraints for data decomposition in order to obtain unique solutions when specific criteria are satisfied [3]. gNCA expands these capabilities by allowing incorporation of constraints onto the deduced TFAs, such as transcription factor knockout experiments, which offer a rich source of data and biochemical information. Development of these methodologies significantly expands the capabilities of transcriptional regulation analysis. With gNCA, we analyzed the combined wild-type and *fkh1 fkh2* mutant data set and showed that gNCA can be used to identify TFAs which are consistent with cell physiology, transcription factors with potential cell cycle dependent roles, as well as interactions between transcription factors.

On the basis that transcription factors exhibiting similar activity patterns function together, we identified 11 transcription factors that clustered closely with the known cell cycle regulators. We performed a periodicity test to determine the TFA profiles that exhibit periodic behavior, and by combining the sets of transcription factors collected by these two methods, we identified 5 putative cell cycle-related regulators: Dal81, Hap2, Hir2, Mss11, and Rlm1. These transcription factors may participate in functions driven by cell cycles, or may regulate cell cycle directly or indirectly.

Our comparison between the wild-type TFAs and the mutant TFAs confirmed that the forkhead transcription factors interact with Ace2, Ndd1, and Swi5. This result is consistent with previous reports [16, 17]. Using this approach, we identified 4 additional transcription factors that may functionally interact with Fkh1 Fkh2 directly or indirectly: Hap2, Rts2, Cha4, and Fhl1. Most of these transcription factors are not known to be related to cell cycle, suggesting that cell cycle regulation interacts with other physiological functions.

It is worth noting that our analysis can be sensitive to errors in the connectivity graph. Through a sensitivity analysis we determined that all the known cell cycle regulators and all the known forkhead interaction partners have TFAs that exhibit low sensitivity to the connectivity network when using 0.5 as a correlation coefficient threshold. These results suggest that our analysis is robust to errors in connectivity. With the same sensitivity criterion, 2 (Dal81 and Rlm1) of the 5 putative cell cycle-related regulators and 1 (Fhl1) of the 4 putative forkhead interaction partners were determined to be robust to connectivity errors. The lack of sensitivity to error increases confidence in these predictions.

## Conclusion

Protein function and interaction are often deduced computationally through analysis of genomic sequence, protein sequence, domain architecture, phylogenic profile, or gene expression level. Gene expression analysis represents the only current technique that takes into account dynamic behavior to assign interaction and function. However, it has been shown that for proteins significantly regulated post-transcriptionally (e.g. transcription factors) their gene expression levels do not correlate well with their activity. Here a method for screening the physiological roles of transcription factors and their functional interactions based on their dynamic activity profiles was introduced. Our method first determines TFAs, which were then further analyzed with cluster and periodicity analysis to determine transcription factor function. By combining the results from both cluster analysis and periodicity analysis, we recovered more than 90% of transcription factors that are proposed to be involved in cell cycle regulation. In addition, we discovered 5 putative transcription factors that may be related to cell cycle regulation. Functional interactions between transcription factors were determined by isolating statistically significant perturbations in TFA patterns between mutant Δ*fkh1* Δ*fkh2* and wild-type experiments carried out under the same experimental conditions. This method allowed for the identification of 4 novel and 3 previously verified fork-head transcription factor interaction partners. We recognize that the Chip-chip data may be condition-dependent and noisy. However, together with microarray data, they provide useful information for functional deduction of transcription regulators.

## Methods

### Gene expression data and connectivity information

In this study, we utilized microarray data sets which were taken from wild-type *S. cerevisiae* cultures synchronized by three independent methods, α-factor arrest, elutriation, and arrest of a *cdc15* temperature-sensitive mutant [13] and a mutant that lacks in two forkhead transcription factors [11], Fkh1 and Fkh2 synchronized by α-factor arrest. The two data sets were combined after imputation [21] to fill in missing points, and common genes without missing data points between the two were selected for NCA. The wild type and the mutant data were organized side-by-side in different columns of the data matrix. The total number of microarray experiments in the combined dataset was 69.

The connectivity information between transcription factors and their regulated genes comes from the genome-wide location or ChIP-chip assay [14]. In this study, the *p*-value threshold used was 0.001 for determining connectivity. As with any other method for determining genome wide transcriptional networks, the ChIP-chip assay suffers from miss-connectivity issues. Both under- and over-prediction errors in connectivity can be rationalized to occur in genome wide location analysis. Since NCA deductions are dependent on the connectivity graph generated by ChIP-chip assays, we endeavoured to determine the effects of connectivity errors on our TFA profiles.

### gNCA

NCA and its extension, gNCA, are a recent-developed technique that is capable of deducing TFAs based on gene expression. Compared to NCA, a distinct feature of gNCA is the capability of imposing constraints on both the regulatory network and regulatory signals. This considerable expands the range of network topologies capable of being analyzed with NCA, and thus allows application of the method to a much larger expanse of scenarios. The analysis starts with defining a model that expresses gene expressions as a function of TFA and the corresponding control strengths (CS). The resulting model can be written in a log-log form. In matrix notation, this model is:

**E = A × P + Γ**, **(1)**

where *E*_{
ij
}= log(*g*_{
i
}(*t*_{
j
}) / *g*_{
i
}(*t*_{0})) and *g*_{
i
}(*t*_{
j
}) is the expression of *i*-th gene evaluated at time *t*_{
j
}**, A** is the regulatory network of control strengths where *A*_{
ij
}denotes the control strength of transcription factor *j* on gene *i* (*A*_{
ij
}*= CS*_{
ij
}), *P*_{
ij
}= log(TFA_{
i
}(*t*_{
j
}) / TFA_{
i
}(*t*_{0})) with TFA_{
i
}(*t*_{
j
}) being the *i*-th transcription factor activity evaluated at time *t*_{
j
}, and **Γ** represents external stimulus and stochastic background noise from the DNA microarray. The dimensions of **E**, **A** and **P** are (*N* × *M*), (*N* × *L*) and (*L* × *M*), respectively, where *N* is the number of genes in the network, *M* is the number of data points or experiments conducted, and *L* is the number of transcription factors used in the study. Since most genes are regulated by only a subset of transcription factors, **A** is generally sparse. This sparse connectivity pattern in **A** is defined by *Z*_{A}, which can be obtained from ChIP-chip assay or existing databases. In addition, gNCA can impose constraints on **P**. For example, when a gene that codes for a specific transcription factor is deleted, the TFA profile (represented in the P matrix) is set to zero. The corresponding zero pattern of **P** is denoted as *Z*_{P}. If *Z*_{A} and *Z*_{P} satisfy a given set of conditions, then gNCA can decompose the data **E** into **A** and **P** up to a diagonal scaling factor. This condition is called *essential uniqueness* [3, 22]. Additional technical details may be found elsewhere.

### Selection of multiple sub-networks for gNCA

Among the set of required conditions for an *essentially unique* decomposition, a necessary condition is that the number of transcription factors (*L*) must be less than the number of data points (*M*) in a given network. Therefore, to analyze a large network in yeast, we dissected the regulatory network into multiple sub-networks and recombine them at the end of the analysis. Each sub-network was constructed such that the number of transcription factors was less than the number of time data points.

Satisfaction of the NCA identifiablity criteria [3, 22] for each sub-network requires that each sub-network contain some overlapping TFs to ensure *consistency* among the computed dynamics. This is possible because the nature of the transcriptional regulatory network as determined by genome-wide location analysis data reveals that some transcript factors highly overlap the set of genes on the microarray [14]. Including these transcript factors in each sub-network decomposition step achieves the required consistency. First we select randomly *k* transcript factors from the column of A as an initial sub-network (*S'*). If this sub-network fails a rank check we then remove the *l* TFs that cause the rank deficiency, and replace them with *l* new TFs that were not used to create *S'*. The top *l* TFs are chosen by ranking the TFs by the total number of genes they control in descending order. If this continues to fail greater than a specified number of iterations, then we change *k* (usually to a smaller number) and repeat the process. For large enough *k* (40 TFs) it is likely that the many genes overlap because of the *a priori* network structure. This is by no means the only way to construct sub-networks, nor is it the optimal way. We do not guarantee that this will work in all cases, but for our purposes it is sufficient to generate consistent TFAs.

Among the combined data set with imputed missing data points, a total of 1818 genes with connectivity information of 104 TFs from genome wide location analysis were used to generate random sub-networks. Sub-networks with 40 transcription factors from total 104 TFs were generated and examined for gNCAuniqueness conditions. Four gNCA-compliant sub-networks, with 1110, 847, 1015, and 795 genes respectively, were chosen such that total number of TFAs deduced from 4 random sub-networks could be greater than the number of data points (M = 69). For each sub-network containing Fkh1 or Fkh2, the corresponding TFAs in the **P** matrix corresponding to the deletion mutant were constrained to zero.

### Cluster analysis of TFA profiles

The identification of potential new cell cycle factors was conducted by using hierarchical clustering on the TFA profiles. Specifically, hierarchical clustering based on absolute value of Pearson correlation coefficient, which measures the strength and direction of a linear relationship between two variables, was applied to TFA profiles deduced by gNCA. We hypothesize that transcription factors with similar activities function together. Therefore, transcription factors which are clustered together with 11 known cell cycle factors could be considered as possible cell cycle-dependent transcription factors. Publicly available software packages, Cluster and Treeview [12] were used for this test.

### Periodicity analysis of TFA profiles

We identify periodic functions from stochastic background noise by employing a statistical analysis technique described previously [4]. This technique is based on rejecting the null hypothesis of a purely random process through considering the power spectrum of a time dependent signal. First, the technique computes the power spectrum of a given signal and evaluates the *g*-statistics expressed as the contribution of the power spectrum at a specific frequency to the total intensity of the power spectrum. Large values of *g*-statistics suggest that the underlying process is periodic. Then, the signals are screened for periodic motion using multiple testing under the criterion of false discovery rate (FDR). The R package GeneTS. [4] was used for this study.

### Statistical analysis for interaction determination

The effect of *fkh1* and *fkh2* mutations on TFAs is elucidated by examining the correlation between TFAs from the wild type and the mutant under the same experimental condition. If there is a significant difference in TFAs between the wild type and mutant, the correlation coefficient should be low. However, low correlation coefficients may be caused by low overall activities, since in those situations the TFAs are dominated by noise. Therefore, we need to define a new deviation coefficient of TFAs between two strains to systematically identify transcription factors whose activity profiles are affected by *fkh1 fkh2* mutation. This new deviation coefficient is defined as,

where superscripts *wt* and *mt* represent wild type and *fkh1 fkh2* mutant, respectively, and the overhead bars represent the vector magnitude of TFA. Small values of the deviation coefficient suggest activities dominated by noise. While large values suggest significant activities and significant differences in TFA profiles between mutant and wild-type.

## Declarations

### Acknowledgements

This work is supported by NSF grant BES-0120359, and Center for Cell Mimetic Space Exploration (CMISE), NASA University Research, Engineering and Technology Institute (URETI), under award number #NCC 2-1364.

## Authors’ Affiliations

## References

- Nelson DM, Ye X, Hall C, Santos H, Ma T, Kao GD, Yen TJ, Harper JW, Adams PD: Coupling of DNA synthesis and histone synthesis in S phase independent of cyclin/cdk2 activity. Mol Cell Biol. 2002, 22: 7459-7472. 10.1128/MCB.22.21.7459-7472.2002.PubMedPubMed CentralView ArticleGoogle Scholar
- Tran LM, Brynildsen MP, Kao KC, Suen JK, Liao JC: gNCA: a framework for determining transcription factor activity based on transcriptome: identifiability and numerical implementation. Metab Eng. 2005, 7: 128-141. 10.1016/j.ymben.2004.12.001.PubMedView ArticleGoogle Scholar
- Liao JC, Boscolo R, Yang YL, 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
- Wichert S, Fokianos K, Strimmer K: Identifying periodically expressed transcripts in microarray time series data. Bioinformatics. 2004, 20: 5-20. 10.1093/bioinformatics/btg364.PubMedView ArticleGoogle Scholar
- Wang W, Cherry JM, Nochomovitz Y, Jolly E, Botstein D, Li H: Inference of combinatorial regulation in yeast transcriptional networks: a case study of sporulation. Proc Natl Acad Sci U S A. 2005, 102: 1998-2003. 10.1073/pnas.0405537102.PubMedPubMed CentralView ArticleGoogle Scholar
- Wang W, Cherry JM, Botstein D, Li H: A systematic approach to reconstructing transcription networks in Saccharomycescerevisiae. Proc Natl Acad Sci U S A. 2002, 99: 16893-16898. 10.1073/pnas.252638199.PubMedPubMed CentralView ArticleGoogle Scholar
- Bussemaker HJ, Li H, Siggia ED: Regulatory element detection using correlation with expression. Nat Genet. 2001, 27: 167-171. 10.1038/84792.PubMedView ArticleGoogle Scholar
- Conlon EM, Liu XS, Lieb JD, Liu JS: Integrating regulatory motif discovery and genome-wide expression analysis. Proc Natl Acad Sci U S A. 2003, 100: 3339-3344. 10.1073/pnas.0630591100.PubMedPubMed CentralView ArticleGoogle Scholar
- Bluthgen N, Kielbasa SM, Herzel H: Inferring combinatorial regulation of transcription in silico. Nucleic Acids Res. 2005, 33: 272-279. 10.1093/nar/gki167.PubMedPubMed CentralView ArticleGoogle Scholar
- Keles S, van der Laan M, Eisen MB: Identification of regulatory elements using a feature selection method. Bioinformatics. 2002, 18: 1167-1175. 10.1093/bioinformatics/18.9.1167.PubMedView ArticleGoogle Scholar
- Zhu G, Spellman PT, Volpe T, Brown PO, Botstein D, Davis TN, Futcher B: Two yeast forkhead genes regulate the cell cycle and pseudohyphal growth. Nature. 2000, 406: 90-94. 10.1038/35021046.PubMedView ArticleGoogle Scholar
- Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci U S A. 1998, 95: 14863-14868. 10.1073/pnas.95.25.14863.PubMedPubMed CentralView ArticleGoogle Scholar
- Spellman PT, Sherlock G, Zhang MQ, Iyer VR, Anders K, Eisen MB, Brown PO, Botstein D, Futcher B: Comprehensive identification of cell cycle-regulated genes of the yeast Saccharomyces cerevisiae by microarray hybridization. Mol Biol Cell. 1998, 9: 3273-3297.PubMedPubMed CentralView ArticleGoogle Scholar
- Lee TI, Rinaldi NJ, Robert F, Odom DT, Bar-Joseph Z, Gerber GK, Hannett NM, Harbison CT, Thompson CM, Simon I, Zeitlinger J, Jennings EG, Murray HL, Gordon DB, Ren B, Wyrick JJ, Tagne JB, Volkert TL, Fraenkel E, Gifford DK, Young RA: Transcriptional regulatory networks in Saccharomyces cerevisiae. Science. 2002, 298: 799-804. 10.1126/science.1075090.PubMedView ArticleGoogle Scholar
- Spector MS, Raff A, DeSilva H, Lee K, Osley MA: Hir1p and Hir2p function as transcriptional corepressors to regulate histone gene transcription in the Saccharomyces cerevisiae cell cycle. Mol Cell Biol. 1997, 17: 545-552.PubMedPubMed CentralView ArticleGoogle Scholar
- Futcher B: Transcriptional regulatory networks and the yeast cell cycle. Curr Opin Cell Biol. 2002, 14: 676-683. 10.1016/S0955-0674(02)00391-5.PubMedView ArticleGoogle Scholar
- Simon I, Barnett J, Hannett N, Harbison CT, Rinaldi NJ, Volkert TL, Wyrick JJ, Zeitlinger J, Gifford DK, Jaakkola TS, Young RA: Serial regulation of transcriptional regulators in the yeast cell cycle. Cell. 2001, 106: 697-708. 10.1016/S0092-8674(01)00494-9.PubMedView ArticleGoogle Scholar
- Gao F, Foat BC, Bussemaker HJ: Defining transcriptional networks through integrative modeling of mRNA expression and transcription factor binding data. BMC Bioinformatics. 2004, 5: 31-10.1186/1471-2105-5-31.PubMedPubMed CentralView ArticleGoogle Scholar
- Ubersax JA, Woodbury EL, Quang PN, Paraz M, Blethrow JD, Shah K, Shokat KM, Morgan DO: Targets of the cyclin-dependent kinase Cdk1. Nature. 2003, 425: 859-864. 10.1038/nature02062.PubMedView ArticleGoogle Scholar
- Doree M, Hunt T: From Cdc2 to Cdk1: when did the cell cycle kinase join its cyclin partner?. J Cell Sci. 2002, 115: 2461-2464.PubMedGoogle Scholar
- Herrero J, Al-Shahrour F, Diaz-Uriarte R, Mateos A, Vaquerizas JM, Santoyo J, Dopazo J: GEPAS: A web-based resource for microarray gene expression data analysis. Nucleic Acids Res. 2003, 31: 3461-3467. 10.1093/nar/gkg591.PubMedPubMed CentralView ArticleGoogle Scholar
- Boscolo R, Sabatti C, Liao JC, Roychowdhury VP: Reconstructing hidden regualtory layers by network component analysis: Theory and application. IEEE Trans on Computational Biology and Bioinformatics. 2004Google Scholar

## Copyright

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.