 Research
 Open Access
 Published:
Predicting transcription factor activities from combined analysis of microarray and ChIP data: a partial least squares approach
Theoretical Biology and Medical Modelling volume 2, Article number: 23 (2005)
Abstract
Background
The study of the network between transcription factors and their targets is important for understanding the complex regulatory mechanisms in a cell. Unfortunately, with standard microarray experiments it is not possible to measure the transcription factor activities (TFAs) directly, as their own transcription levels are subject to posttranslational modifications.
Results
Here we propose a statistical approach based on partial least squares (PLS) regression to infer the true TFAs from a combination of mRNA expression and DNAprotein binding measurements. This method is also statistically sound for small samples and allows the detection of functional interactions among the transcription factors via the notion of "meta"transcription factors. In addition, it enables false positives to be identified in ChIP data and activation and suppression activities to be distinguished.
Conclusion
The proposed method performs very well both for simulated data and for real expression and ChIP data from yeast and E. Coli experiments. It overcomes the limitations of previously used approaches to estimating TFAs. The estimated profiles may also serve as input for further studies, such as tests of periodicity or differential regulation. An R package "plsgenomics" implementing the proposed methods is available for download from the CRAN archive.
Background
The transcription of genes is regulated by DNA binding proteins that attach to specific DNA promoter regions. These proteins are known as transcriptional regulators or transcription factors and recruit chromatinmodifying complexes and the transcription apparatus to initiate RNA synthesis [1, 2].
In the last few years, considerable efforts have been made by both experimental and computational biologists to identify transcription factors, their target genes and the sensitivity of the regulation mechanism to changes in environment [3–5]. An important technique for the identification of target genes bound in vivo by known transcription factors is the combination of a modified chromatin immunoprecipitation (ChIP) assay with microarray technology, as proposed by Ren et al. [1]. For instance, in the budding yeast Saccharomyces cerevisiae, ChIP experiments have been utilized to elucidate the binding interactions between 6270 genes and 113 preselected transcription factors [2]. However, as physical binding of transcription factors is a necessary but not a sufficient condition for transcription initiation, ChIP data typically suffer from a large proportion of false positives.
Several attempts have also been made to recover the network structure between transcription factors and their targets using only the gene expression levels of both the transcription factors and the targets, either with [6] or without [7] assuming a subset of putative regulators. Such approaches implicitly assume that the measured gene expression levels of the transcription factors reflect their actual activity. However, owing to various complex posttranslational modifications as well as to interactions among transcription factors themselves, regulator transcription levels are generally inappropriate proxies for transcription factor activities (TFA).
In a few recent papers, integrative analysis of gene expression data and ChIP connectivity data has been suggested as a way of overcoming these difficulties [8]. Most prominently, Liao and coworkers have developed the technique of "network component analysis" (NCA) [9, 10], a dimension reduction approach to inferring the true regulatory activities. In NCA one can also incorporate further a priori qualitative knowledge about genetranscription factor interactions [11]. Unfortunately, a major drawback of the original NCA method is that for identifiable reasons it imposes very strong restrictions on the network topologies allowed, which renders application of classic NCA difficult in many practical cases. Alter and Golub [12] introduced an approach for integrating ChIP and microarray data using pseudoinverse projection. Like NCA, this method is based on algebraic matrix decomposition (in this case singular value decomposition). However, this ignores measurement and biological errors present in both connectivity and gene expression data. Kato et al. [13] proposed yet another integrative approach consisting of several steps combining sequence data, ChIP data and gene expression data. However, here gene expression is used only to check the coherence of expression profiles of genes with common sequence motifs, and not to estimate transcription factor activities. Finally, Gao et al. [14] suggested the "MANetworker" algorithm, which employs multivariate regression to estimate TFAs and backward variable selection to identify the active transcription factors. Unlike the other approaches, it takes full account of stochastic error. However, for classical regression theory to be valid it is necessary not only that the number of gene targets is much greater than both the number of samples and the number of transcription factors, but also that the transcription factors are independent of each other. The latter condition in particular is clearly not generally satisfied with genome data.
Here, we suggest an alternative statistical framework to tackle the problem of network component and regulator analysis. Our approach centers around multivariate partial least squares (PLS) regression, a wellknown analysis tool for highdimensional data with many continuous response variables that has been widely applied, especially to chemometric data [15–17]. Using PLS we are able not only both to integrate and generalize previous NCA approaches, but also to overcome their respective limitations. In particular, PLSbased network component analysis offers a computationally highly efficient and statistically sound way to infer true TFAs for any given connectivity matrix. In addition, it allows statistical assessment of the available connectivity information, and also the discovery of interactions and natural groupings among regulatory genes (corresponding to "meta"transcription factors).
Results
Network model
Suppose gene expression data for n genes and m samples (= arrays, tissue types, time points etc.) are collected in a n × m data matrix . Furthermore, let denote the socalled connectivity matrix with n rows and p columns. Each column in describes the strength of interaction between one of p transcription factors and the n considered gene targets. The entries of can either be binary (0–1) or numeric (e.g. ChIP data), with a zero value indicating no physical binding between a transcription factor and a target.
In order to relate expression to connectivity data we consider the linear model
where A is a n × m constant matrix, is a p × m matrix of regression coefficients and E is a n × m matrix containing error terms. A contains the m different offsets, and may be interpreted as the matrix of the true transcription factor activities (TFAs) of the p transcription factors for each of the m samples.
It is worth noting that in this setting, unlike in most other gene expression analysis studies, the number of genes n is considered as the number of cases rather than the number of variables. In the present case the latter corresponds to the number of transcription factors p (hence, in general, p <n).
NCA and MANetworker algorithms
The above model linking TFAs both with gene expression of the regulated genes and external connectivity information has been the subject of a series of recent studies.
In the classic network component analysis approach [9, 10] the offset matrix A is set to zero and the remainder of Eq. 1 is interpreted as a dimension reduction that projects the output layer with m samples on to a "hidden" layer of p <m transcription factors. In the original NCA algorithm the coefficients are obtained via a novel matrix decomposition that respects the zero pattern constraint given in the connectivity matrix . Unfortunately, this also imposes rather strict identifiability conditions. As a consequence, classic NCA may only be employed with certain classes of "NCA compatible" [9].
In contrast, the "MANetworker" algorithm by Gao et al. [14] employs standard multiple leastsquares regression in conjunction with stepwise variable selection to estimate the true transcription factor activities . This requires that the number of target genes is much larger than both the number of transcription factors and the number of samples. More important, however, is that the stepwise model selection procedure employed is only poorly suited if the regulator genes are themselves interacting with each other. This is a major drawback as it is biologically wellknown that transcription factors often work in conjunction with other regulators, and rarely act independently.
Partial least squares regression
Here we propose to employ the method of partial least squares regression [15] to infer true TFAs and the functional interactions of regulators.
PLS is a wellknown analysis tool for highdimensional data with many continuous response variables that has been widely applied, especially to chemometric data [17]. PLS is particularly suited to the case of nonindependent predictors and for smallsample regression settings [16, 18–20]. It is computationally highly efficient, it does not necessitate variable selection, and it additionally infers meaningful structural components.
For these reasons PLS is now being adopted as a standard tool for multivariate microarray data analysis, particularly in classification problems [21–24]. We believe that PLS also provides an excellent framework for integrative network analysis, as it combines dimension reduction with regression and variable selection, the two key elements from both the NCA and the MANetworker approaches.
In a nutshell, the PLS algorithm consists of the following consecutive steps:

1.
First, the data matrices
and are centered to column mean zero, resulting in matrices X and Y, in order to estimate and to remove the offset A. In addition, it is common practice in PLS analysis (and also recommended here) to scale the input matrices to unit variance.

2.
Second, using the linear dimension reduction T = XR, the p predictors in X are mapped onto c ≤ rank(X) ≤ min(p, n) latent components in T (an n × c matrix). See the section "SIMPLS algorithm" below for the precise procedure employed in this paper. The important key idea in PLS is that the weights R (a p × c matrix) are chosen with the response Y explicitly taken into account, so that the predictive performance is maximal even for small c.

3.
Next, assuming the model Y = TQ' + E, Y is regressed by ordinary least squares against the latent components T (also known as Xscores) to obtain the loadings Q (a m × c matrix), i.e. Q = Y'T(T'T)^{1}.

4.
Subsequently, the PLS estimate of the coefficients B in Y = XB + E is computed from estimates of the weight matrix R and the Yloadings Q via B = RQ'.

5.
Finally, the coefficients
for the original Eq. 1 are computed by rescaling B.
Note that it is step 2 that mostly distinguishes PLS from related bilinear regression approaches such as principal and independent components regression (PCR/ICR) and the pseudoinversebased method of Alter and Golub [12]. In the latter approaches the scores T are computed solely on the basis of the data matrix X without considering the response Y [16].
Other quantities often considered in PLS include, e.g., the Xloadings P that are obtained by regressing X against T, i.e. X = TP' + F and P = X'T(T'T)^{1}.
SIMPLS algorithm
PLS aims to find latent variables T that simultaneously explain both the predictors X and the response Y. The original ideas motivating the PLS decomposition were entirely heuristic. As a result, a broad variety of different, but in terms of predictive power equivalent, PLS algorithms have emerged – for an overview see e.g. Martens [17].
For the present application to infer true TFAs, we suggest using the SIMPLS ("Statistically Inspired Modification of PLS") algorithm, which has the following appealing properties [18–20]:

it produces orthogonal, i.e. empirically uncorrelated, latent components;

it allows for a multivariate response; and

it optimizes a simple statistical criterion.
A further added advantage of SIMPLS is that it is also one of the most computationally efficient PLS algorithms.
We note that other PLS variants described in the literature have predictive power comparable to SIMPLS. However, these either provide orthogonal loadings rather than orthogonal latent components T (Martens' PLS), or they do not elegantly extend from 1dimensional to mdimensional responses Y in terms of their optimized objective function (NIPALS).
In SIMPLS, the latent components t_{1}, t_{2},..., t_{ c }of the columns in T are inferred by sequentially estimating the column vectors r_{1},..., r_{ c }of R according to the following criterion [20]:

1.
r_{1} is the unit vector (with r_{1} = 1) maximizing the length Y' Xr_{1} of the m × 1 covariance vector cov(Y, t_{1}).

2.
For all j = 2,...,c, r_{ j }are the unit vectors (with r_{ j } = 1) maximizing the length Y' Xr_{ j } of the vector cov(Y, t_{ j }) subject to the orthogonality constraint for all i = 1,...,j  1.
In the actual SIMPLS procedure, the weights R and the derived quantities T and Q are obtained by a GramSchmidttype algorithm [18].
On a practical note, we would like to mention that in many implementations of SIMPLS (e.g. in the "pls.pcr" R package by Ron Wehrens, University of Nijmegen), conventions different from the above are used. In particular, the Xscores T* returned will often be orthonormal (rather than orthogonal) and consequently the weights R* will not have unit norm as in our case. For conversion, define M = diag(,...,,) and set T = T*M^{1}, R = R*M^{1}, Q = Q*M, and P = P*M. This provides orthogonal scores and unitnorm weights as assumed in our description of SIMPLS.
The resulting estimates of the matrices B, T, and R are now straightforward to interpret in terms of transcriptional regulation. B (and ) give the inferred activities of the p transcription factors in each of the m experiments. The inferred latent components T describe "meta"transcription factors that combine related groups of transcription factors. R reflects the involvement of each of the p regulators in the c metafactors.
Determining the number of PLS components
A remaining aspect of PLS regression analysis is the optimal choice of the number c of latent components. If the maximal value c_{max} = rank(X) is chosen, then PLS becomes equivalent to principal components regression (PCR) with the same number of components, and if additionally n >p both PLS and PCR turn into ordinary leastsquares multiple regression.
Hence, with PLS it is desirable to choose as small a value of c as possible without sacrificing too much predictive power. One straightforward statistical procedure to estimate this minimum value c_{min} is the method of crossvalidation, which proceeds as follows (cf. also refs. [25] and [26]):

1.
Split the set of n genes randomly into 2 sets: a learning set containing 2/3 of the genes and a test set containing the remaining genes.

2.
Use the learning set to determine the matrix of regression coefficients B for different values c = 1, 2,...,c_{max}.

3.
Predict the gene expression of the n/3 genes from the test set using B with the different values of c.

4.
Repeat steps 1–3 K = 100 times and compute the mean squared prediction error for each c.
Subsequently, the value of c yielding the smallest mean squared prediction error is selected.
Alternatively, the optimal number of components may also be determined by considering the value of the criterion Z_{ i }= Y't_{ i } for a given latent component t_{ i }. If Z_{ i }falls below an a priori specified threshold then c_{min} = i is reached.
Discussion
Data sets
Next, we illustrate the versatility of the proposed PLS approach to network component analysis by analyzing several real biomolecular data sets.
First, in order to validate the linear regression approach (Eq. 1) we reanalyzed hemoglobin data from Liao et al. [9]. Second, we analyzed two different S. Cerevisiae gene expression data sets in conjunction with a regulatortarget connectivity matrix from the largescale ChIP experiment of Lee et al. [2]. The yeast expression data investigated comprise a time series experiment from Spellman et al. [27] and a compilation of yeast stress response experiments from Gasch et al. [6, 28]. Finally, we analyzed expression and connectivity data for an E. Coli regulatory network containing 100 genes and 16 transcription factors from Kao et al. [10]. The general characteristics of these four data sets are summarized in Table 1.
The data investigated were preprocessed as follows. The yeast ChIP data set [2] contains proteinDNA interaction data for 6270 genes and 113 transcription factors. It includes missing values that correspond to noninteracting genetranscription factor pairs. Although ChIP data are essentially continuous, it is common practice to dichotomize them according to the pvalues into discrete levels of interaction (0 or 1). In this study, we used data obtained at a pvalue threshold of 0.001, as suggested by Lee et al. [2]. However, note that in contrast to the NCA method, dichotomization of the ChIP data is optional in our approach.
The Spellman et al. [27] microarray data originally contained the gene expression of 4289 genes at 24 time points during the cellcycle. From these genes, a subset of 3638 are also contained in the Lee et al. [2] ChIP data set. Our analysis is based on these 3638 genes. Similarly, the Gasch expression data set [6, 28] contains the expression of 2292 genes for 173 arrays corresponding to different stress conditions (e.g. heat shock, amino acid starvation, nitrogen depletion). Of these 2292 genes, a subset of 1993 overlap with the genes considered in the ChIP data.
The connectivity matrix for the E. coli data was compiled mainly by Kao et al. [10] from the RegulonDB [11] database. In addition, they incorporated a few corrections using literature data. The temporal E. coli expression data for 100 genes across 25 time points was also introduced in Kao et al. [10] and is publicly available at http://www.seas.ucla.edu/~liaoj/.
Validation of the regression approach
The hemoglobin data used in Liao et al. [9] for validation of the classic NCA approach have the advantage that the true coefficients of the network model in Eq. 1 are known, and therefore can be directly compared with the inferred values.
Reanalyzing these data, it is straightforward to show (see Figure 1) that the true regression coefficients can be recovered exactly by multivariate regression (of which PLS is a special case). According to Liao et al. [9], this is also true for classic NCA but not for PCA and ICA interpretations of Eq. 1. This discrepancy can be explained by the fact the neither PCA nor ICA explicitly takes account of the response Y, whereas NCA and PLS do.
PLS components and Yloadings
Subsequently, we determined the minimum number of PLS components for the yeast and E. coli data sets using crossvalidation. The results are plotted in Figure 2 (top) after normalization (the mean crossvalidation error with one PLS component is set to one). As can be seen from Figure 2, the minimal mean crossvalidation error is obtained with 5 PLS components for the Spellman data, 8 PLS components for the Gasch data and 2 PLS components for the E. coli data. For comparison, the (normalized) objective criterion Y't_{ i } of the SIMPLS algorithm is also represented on Figure 2 (bottom) for different numbers of PLS components. These results are in good agreement with the crossvalidation error: it increases when PLS components with a low objective criterion are added.
The Yloadings contained in the m × c matrix Q give the projection of the c "meta"transcription factors for each of the m experiments. As can be seen from Figure 3 for the Spellman data, both the first and the third metafactors explain the periodic part of the expression data, but with different phases. The second metafactor corresponds to small oscillations with very short period, whereas the fourth and fifth metafactors reflect longtime trends (slow and stepwise increasing, respectively). Using Fisher's gtest as proposed in Wichert et al. [29], we detected statistically relevant periodicity for the four first metafactors. In Figure 3, the Yloadings are also represented for the E. coli data. Whereas the projection of the first metafactor is approximately constant over time, the projection of the second metafactor increases strongly and (almost) uniformly. Thus, in both data sets, the PLS algorithm allows us to extract metafactors from the data corresponding to distinct latent trends.
For the Gasch data, the m experiments do not correspond to different time points but to 13 different stress conditions (see Gasch et al. [28] for further details, and Table 2 for the list of conditions). In this case the Yloadings may be interestingly analyzed using Wilcoxon's rank sum test. For each condition k and each metafactor j, we tested the H_{0} hypothesis that the median of the projection of the jth metafactor is the same in condition k as in all the other conditions ({1,..., k  1, k + 1,..., 13}). In this situation, Wilcoxon's rank sum test is preferable to the wellknown twosample ttest, because some of the conditions include only a very small number of experiments. The results obtained with a pvalue threshold of 0.05 are displayed in Table 2. The entries 1 and 0 correspond to significant and nonsignificant (FDR adjusted) pvalues, respectively. As can be seen from Table 2, each PLS component carries a particular pattern of associated significant conditions, indicating that the metafactors capture a distinct direction of the data.
Inferred transcription factor activities
One of the main objectives of our PLSbased approach is to estimate the true transcription factor activities (TFAs). Although all the TFAs can be estimated in the same way for the three data sets, we display only the evolution over time of a few interesting TFAs for the two time series data sets (i.e. the Spellman and the E. coli data).
The TFAs (top) and expression profiles (bottom) of 4 wellknown cellcycle regulators are depicted in Figure 4 for the Spellman data. The TFAs of MCM1, SWI4, SWI5 and ACE2 show highly periodic patterns, which is consistent with common biological knowledge. In contrast, the expression profiles of MCM1 and SWI4 are not periodic (this can be confirmed by Fisher's gtest [29]). On the other hand, the expression profiles of SWI5 and ACE2 are periodic, though not with the same phase as the inferred TFAs. This may indicate either an inhibiting or a phaseshift effect of the transcription factors on the regulated genes.
The remainder of the TFAs and the regulated genes were also tested for periodicity using the gtest [29]. After FDR adjustment of the pvalues, we found that 62 of the 113 transcription factors (= 55%) in the Spellman/Lee data have significantly periodic TFAs at the level 0.05. In contrast, only 804 of the 4289 genes (= 19%) exhibit significantly periodic expression profiles.
For the E. coli data the time profiles of the estimated TFAs of the 16 transcription factors are represented in Figure 5. The TFAs of ArcA, GatR, Lrp, PhoB, PurR, RpoS decrease over time, those of CRP, CysB, FadR, IcIR, NarL, RpoE, TrpR and TyrR remain approximately constant and those of FruR and LeuO increase strongly. This is consistent with previous results obtained by NCA [10]. We point out, however, that unlike NCA our approach may be applied to any arbitrary network topology, whereas the present E. coli network was chosen specifically to meet the NCA compatibility criteria [9].
As can be seen already from the few examples depicted in Figure 4, the TFAs do not always correlate with the respective expression profiles. We tested this for all the transcription factors of which the expression profiles were also included in the data sets. For the Gasch data, we found that only 63 from the 90 available transcription factors exhibit expression profiles that are correlated with TFAs (at the level 0.05 with FDR pvalue adjustment). For the Spellman time series data, none of the 78 available TFAexpression profile pairs are correlated. These results clearly indicate that methods investigating transcriptional regulation with expression data as their sole basis are likely to miss potentially important regulation activities.
Generegulator coupling factors
Another topic of interest is the identification of false positives in ChIP data. Following Gao et al. [14] we investigate this problem using Pearson's correlation test. For each supposed genetranscription factor pair (according to the dichotomized ChIP data) we test if the inferred TFA is significantly correlated with the expression profile of the regulated gene. For the Gasch data, we find that 73% of the 1495 genetranscription factor pairs are correct (i.e. the TFA is significantly correlated with the expression profile at the 0.05 level with FDR pvalue adjustment). The concordance with the ChIP connectivity information is much worse for the Spellman data, where only 32% of the 2535 genetranscription factor pairs are significantly correlated.
We should like to add as a note of caution that the lack of correlation between TFA and target gene needs to be viewed as specific to the microarray study investigated. Other expression experiments may activate different pathways and thus produce different patterns of correlation in conjunction with the ChIP connectivity information.
Conclusion
Network component analysis combines microarray data with ChIP data with the aim of enhancing the estimation of regulator activities and of connectivity strengths. In this paper we have presented an approach to NCA based on partial least squares, a computationally efficient statistical regression tool.
Our PLS framework allows several drawbacks, inherent both in the classic NCA methods based on matrix decomposition and in the MANetworker algorithm, to be overcome. Its simplicity (no iterative step, no variable selection, no stochastic search) and its flexibility (no distributional assumptions, no topological constraints, no conditions on the dimensions) compared to competing approaches make it particularly attractive as an integrative method for analyzing complex regulatory networks. Moreover, the PLS algorithm not only extracts information on generegulator and on TFAexpression profile pairs but also identifies coherent metafactors reflecting the main directions of variation of the data, taking account both of the expression (
) and the connectivity information ().
Our analysis of biological data shows the versatility of our PLS approach and at the same time dramatically confirms the need for a combined expressionChIP analysis for inferring regulation. Particularly striking are the sometimes drastic differences between the measured transcription levels and the PLSinferred transcription activities. According to Segal et al. [6], some transcription factors may also not be active in all conditions. Note that this assumption is also automatically taken into account by our approach.
NCA in general, and the present PLSbased variant in particular, may be criticized for relying on a simple linear model  see Buchler et al. [30] and Setty et al. [31] for counterexamples. Therefore, more elaborate regression approaches such as generalized linear models (GLMs) or generalized additive models (GAMs) may be required to further enhance our current understanding of how best to model the complex structures governing genetic networks.
[Note added in proof: See Yang et al. [32] for a related study in the sister journal BMC Genomics.]
Appendix: Computer program
All algorithms have been implemented in the R language [33]. A corresponding R package "plsgenomics" developed by the authors is available for download from the CRAN archive http://cran.rproject.org.
References
 1.
Ren B, Robert F, Wyrick JJ, Aparicio O, Jennings EG, Simon I, Zeitlinger J, Schreiber J, Hannett N, Kanin E, Volkert TL, Wilson CJ, Bell SP, Young RA: Genomewide location and function of DNA binding proteins. Science. 2000, 290: 23062309. 10.1126/science.290.5500.2306.
 2.
Lee TI, Rinaldi NJ, Robert F, Odom DT, BarJoseph 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: 799804. 10.1126/science.1075090.
 3.
Iyer VR, Horak CE, Scafe CS, Botstein D, Snyder M, Brown PO: Genomic binding sites of the yeast cellcycle transcription factors SBF and MBF. Nature. 2001, 409: 533538. 10.1038/35054095.
 4.
van Steensel B, Delrow J, Bussemaker HJ: Genomewide analysis of Drosophila GAGA factor target genes reveals contextdependent DNAbinding. Proc Natl Acad Sci USA. 2003, 100: 25802585. 10.1073/pnas.0438000100.
 5.
Harbison CT, Gordon DB, Lee TI, Rinaldi NJ, Macisaac KD, Danford TW, Hannett NM, Tagne JB, Reynolds DB, Yoo J, Jennings EG, Zeitlinger J, Pokholok DK, Kellis M, Rolfe PA, Takusagawa KT, Lander ES, Gifford DK, Fraenkel E, Young RA: Transcriptional regulatory code of a eukaryotic genome. Nature. 2004, 431: 99104. 10.1038/nature02800.
 6.
Segal E, Shapira M, Regev A, Pe'er D, Botstein D, Koller D, Friedman N: Module networks: identifing regulatory modules and their conditionspecific regulators from gene expression data. Nature Genetics. 2003, 34: 166176.
 7.
Xiong M, Li J, Fang X: Identification of genetic networks. Genetics. 2004, 166: 10371052. 10.1534/genetics.166.2.1037.
 8.
Li Z, Chan C: Extracting novel information from gene expression data. Trends Biotechnol. 2004, 22: 381383. 10.1016/j.tibtech.2004.06.007.
 9.
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 USA. 2003, 100: 1552215527. 10.1073/pnas.2136632100.
 10.
Kao KC, Yang YL, Boscolo R, Sabatti C, Roychowdhury V, Liao JC: Transcriptomebased determination of multiple transcription regulator activities in Escherichia coli by using network component analysis. Proc Natl Acad Sci USA. 2004, 101: 641646. 10.1073/pnas.0305287101.
 11.
Salgado H, SantosZavaleta A, GamaCastro S, MillanZarate D, DiazPeredo E, SanchezSolano F, PerezRueda E, BonavidesMartinez C, ColladoVides J: RegulonDB (version 3.2): transcriptional regulation and operon organization in Escherichia coli K12. Nucleic Acids Res. 2001, 29: 7274. 10.1093/nar/29.1.72.
 12.
Alter O, Golub GH: Integrative analysis of genomescale data by using pseudoinverse projection predicts novel correlation between DNA replication and RNA transcription. Proc Natl Acad Sci USA. 2004, 101: 1657716582. 10.1073/pnas.0406767101.
 13.
Kato M, Hata N, Banerjee N, Futcher B, Zhang MQ: Identifying combinatorial regulation of transcription factors and binding motifs. Genome Biology. 2004, 5: R5610.1186/gb200458r56.
 14.
Gao F, Foat BC, Bussemaker HJ: Defining transcriptional networks through integrative modeling of mRNA expression and trasncription factor binding data. BMC Bioinformatics. 2004, 5: 3110.1186/14712105531.
 15.
Wold S, Martens H, Wold H: The multivariate calibration method in chemistry solved by the PLS method. Proc. Conf. Matrix Pencils, Lecture Notes in Mathematics. Edited by: Ruhe A, Kagstrom B. 1983, Heidelberg: Springer Verlag, 286293.
 16.
Frank IE, Friedman JH: A statistical view of some chemometrics regression tools. Technometrics. 1993, 35: 109135.
 17.
Martens H: Reliable and relevant modelling of real world data: a personal account of the development of PLS regression. Chemom Intell Lab Syst. 2001, 58: 8595. 10.1016/S01697439(01)001538.
 18.
de Jong S: SIMPLS: An alternative approach to partial least squares regression. Chemom Intell Lab Syst. 1993, 18: 251253. 10.1016/01697439(93)85002X.
 19.
Braak CJFT, de Jong S: The objective function of partial least squares regression. J Chemometrics. 1998, 12: 4154. 10.1002/(SICI)1099128X(199801/02)12:1<41::AIDCEM500>3.0.CO;2F.
 20.
de Jong S, Wise BM, Ricker NL: Canonical partial least squares and continuum power regression. J Chemometrics. 2001, 15: 85100. 10.1002/1099128X(200102)15:2<85::AIDCEM601>3.0.CO;29.
 21.
Datta S: Exploring relationships in gene expressions: a partial least squares approach. Gene Expression. 2001, 9: 249255.
 22.
Nguyen D, Rocke DM: Tumor classification by partial least squares using microarray gene expression data. Bioinformatics. 2002, 18: 3950. 10.1093/bioinformatics/18.1.39.
 23.
Nguyen D, Rocke DM: Partial least squares proportional hazard regression for application to DNA microarray survival data. Bioinformatics. 2002, 18: 16251632. 10.1093/bioinformatics/18.12.1625.
 24.
Boulesteix AL: PLS dimension reduction for classification with microarray data. SAGMB. 2004, 3: 33
 25.
Dudoit S, Fridlyand J, Speed TP: Comparison of discrimination methods for the classification of tumors using gene expression data. J Amer Statist Assoc. 2002, 97: 7787. 10.1198/016214502753479248.
 26.
BragaNeto U, Dougherty ER: Is crossvalidation valid for smallsample microarray classification?. Bioinformatics. 2004, 20: 374380. 10.1093/bioinformatics/btg419.
 27.
Spellman PT, Sherlock G, Zhang MQ, Iyer VR, Anders K, Eisen MB, Brown PO, Botstein D, Futcher B: Comprehensive identification of cell cycleregulated genes of the yeast Saccharomyces cerevisiae by microarray hybridization. Mol Biol Cell. 1998, 9: 32733297.
 28.
Gasch AP, Spellman PT, Kao CM, CarmelHarel O, Eisen MB, Storz G, Botstein D, Brown PO: Genomic expression programs in the response of yeast cells to environmental changes. Mol Biol Cell. 2000, 11: 42414257.
 29.
Wichert S, Fokianos K, Strimmer K: Indentifying periodically expressed transcripts in microarray time series data. Bioinformatics. 2004, 20: 520. 10.1093/bioinformatics/btg364.
 30.
Buchler NE, Gerland U, Hwa T: On schemes of combinatorial transcription logic. Proc Natl Acad Sci USA. 2003, 100: 51365141. 10.1073/pnas.0930314100.
 31.
Setty Y, Mayo AE, Surette MG, Alon U: Detailed map of a cisregulatory input function. Proc Natl Acad Sci USA. 2003, 100: 77027707. 10.1073/pnas.1230759100.
 32.
Yang YL, Suen J, Brynildsen M, Galbraith S, Liao JC: Inferring yeast cell cycle regulators and interactions using transcription factor activities. BMC Genomics. 2005, 6: 9010.1186/14712164690.
 33.
R Development Core Team: R: A language and environment for statistical computing. 2004, R Foundation for Statistical Computing, Vienna, Austria, [ISBN 3900051070].http://www.Rproject.org
Acknowledgements
We thank Eran Segal and James Liao for kindly providing the Saccharomyces cerevisiae data and the hemoglobin data, respectively. We also thank the anonymous referees for useful comments. This work was supported by the Deutsche Forschungsgemeinschaft (DFG) through an EmmyNoether research grant to K.S. and the Sonderforschungsbereich 386.
Author information
Affiliations
Corresponding author
Additional information
Authors' contributions
A.L.B. performed all the data analysis and simulations. Both authors jointly developed the methodology, wrote the manuscript, and approved of the final version.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
About this article
Cite this article
Boulesteix, AL., Strimmer, K. Predicting transcription factor activities from combined analysis of microarray and ChIP data: a partial least squares approach. Theor Biol Med Model 2, 23 (2005). https://doi.org/10.1186/17424682223
Received:
Accepted:
Published:
Keywords
 Partial Little Square
 Principal Component Regression
 Connectivity Matrix
 ChIP Data
 Partial Little Square Component