Skip to main content
  • Research article
  • Open access
  • Published:

Identification of interacting transcription factors regulating tissue gene expression in human

Abstract

Background

Tissue gene expression is generally regulated by multiple transcription factors (TFs). A major first step toward understanding how tissues achieve their specificity is to identify, at the genome scale, interacting TFs regulating gene expression in different tissues. Despite previous discoveries, the mechanisms that control tissue gene expression are not fully understood.

Results

We have integrated a function conservation approach, which is based on evolutionary conservation of biological function, and genes with highest expression level in human tissues to predict TF pairs controlling tissue gene expression. To this end, we have identified 2549 TF pairs associated with a certain tissue. To find interacting TFs controlling tissue gene expression in a broad spatial and temporal manner, we looked for TF pairs common to the same type of tissues and identified 379 such TF pairs, based on which TF-TF interaction networks were further built. We also found that tissue-specific TFs may play an important role in recruiting non-tissue-specific TFs to the TF-TF interaction network, offering the potential for coordinating and controlling tissue gene expression across a variety of conditions.

Conclusion

The findings from this study indicate that tissue gene expression is regulated by large sets of interacting TFs either on the same promoter of a gene or through TF-TF interaction networks.

Background

Transcriptional regulation in eukaryotic organisms is a fundamental process to determine a gene's spatial and temporal expression. One of the main events involved in this process is the binding of TFs to short DNA motifs, called transcription factor binding sites (TFBSs), on the promoter regions of genes, activating or repressing the transcription machinery. In mammalian tissues most TFs do not act alone, but work through combinatorial regulation [1, 2], in which two or more TFs work synergistically to control individual gene expression. This combinatorial regulation is able to increase the specificity and flexibility of genes in controlling tissue development and differentiation. Therefore, one of the major first steps toward understanding how tissues achieve their specificity is to identify interacting TFs regulating gene expression in different tissues.

Early attempts to identify interacting TFs controlling tissue gene expression came from the use of experimental approaches such as gel retardation assay [3], site-directed mutagenesis [4], chromatin immunoprecipitation [5, 6], and genomic microarrays [5, 6] in tissues such as liver [3, 5–8], pancreas [6], immune systems [9, 10], muscle [11–13], and neural stem cells [14]. In these studies, interactions between TFs were discovered on a limited scale. To overcome this limitation, some researchers built models to predict tissue-specific cis-regulatory modules in liver [15, 16] and muscle [17] tissues. Taking advantage of the unprecedented amount of sequence and gene expression information from the most recent technical and experimental advances, a few researchers have developed computational approaches to predict tissue-specific TFs and cis-regulatory modules based on recognizable sequence features from either highly expressed genes [18] or genes expressed only in a particular tissue [19–21] derived from genome-wide gene-expression profiling. Some of these researchers have defined tissue-specific enhancers by combining gene-expression profiling, genome comparison, and TFBS analyses [18] or have predicted TF synergy using the relative position and co-occurrence of TFBSs in the promoters of genes expressed only in a particular tissue [19]. Others have looked for tissue-specific cis-regulatory modules by enrichment analysis for motifs discovered de novo in tissue-specific promoters relative to other promoters from the same species [21]. Despite all these efforts, the mechanism that determines tissue development and differentiation is still not fully understood, as the regulation of tissue gene expression involves complex combinatorial interactions between TFs.

In this study, rather than using sequence features of promoters from genes that are expressed only in a particular tissue [19–21], we used our function conservation approach [22] to predict interacting TFs from the most highly expressed genes in each of 79 human tissues [23]. Our approach predicts interacting TFs by integrating the function conservation of interacting TFs from both their binding sites and target genes between closely related species, which are based on the following two assumptions. The first is based on the strong possibility that functional TFBS pairs have more distance constraint than random co-occurrence of TFBSs. The second relies on the biological assumption that while a TF pair plays the same role in regulating gene expression between closely related species, the occurrence of its binding sites is expected to be more highly enriched in promoter sequences of orthologous genes than in promoter sequences of non-orthologous genes. Other than function conservation, the use of highly expressed genes in a tissue allows one to avoid the elimination of common genes contributing to tissue development and differentiation between tissues, especially for closely related tissues (e.g. skeletal muscle and heart), when compared to the use of tissue-specific genes [19–21] which are expressed only in a particular tissue. To our knowledge, this is the first use of a function conservation approach and highly expressed genes in tissues for interacting TF prediction. Therefore, the findings provide novel insight into how tissue gene expression is controlled.

The application of the function conservation approach to the most highly expressed genes has led to the prediction of hundreds of interacting TFs from each of the 79 human tissues. Based on these predictions, TF pairs associated with a certain tissue were identified. The validity of these discovered TF pairs has been evaluated by both known interacting and liver-specific TFs. We further extended our study to find interacting TFs controlling gene expression in a broad spatial and temporal manner by looking for TF pairs common to the same type of tissues, from which TF-TF interaction networks were further built. As a first step to elucidating cis-regulatory modules involved in tissue gene regulation, we also performed analysis to identify interactions of 3 TFs.

Results

Overall analysis procedures

The overall analysis procedures are shown in Additional File 1 and Figure 1. We employed our previously developed function conservation approach [22] to first search for TF pairs using the top 300 expressed genes (referred to tissue-expressed genes) [18] in each of the 79 human tissues from the GNF Atlas2 gene expression database (gnfAtlas2) [23] and their corresponding mouse orthologous genes (Additional File 1). We also utilized promoter sequences from 1018 human housekeeping genes [24] and their mouse orthologous genes to predict interacting TFs playing ubiquitous roles in different tissues (see Methods).

Figure 1
figure 1

Flowchart of analysis procedures. (a) Identification of tissue and tissue-type TF pairs. (b) Prediction of interactions of three TFs. (c) Reconstruction of tissue-type TF-TF interaction networks. HK: housekeeping genes; hs: human; mm: mouse.

We next filtered out the TF pairs in a particular tissue common to those from housekeeping genes (Figure 1a). The remaining TF pairs (referred to tissue TF pairs) in each tissue were more tissue-specific. The rationale for this filtering is that in each tissue some of the interacting TFs play general roles, since all tissues possess common mechanisms to control the fundamental biological processes. To find interacting TFs controlling tissue gene expression in a broad spatial and temporal manner, we extended the analysis to identify tissue TF pairs common to the same type of tissues (referred to tissue-type TF pairs) as well as interactions of 3 TFs. For the former, we looked for common tissue TF pairs in at least 50% tissues of the same type (Figure 1a). We also built TF-TF interaction networks by joining 2 or more tissue-type TF pairs with one shared TF between TF pairs in the same tissue. TF-TF interaction networks with the same topology in at least 2 tissues from the same tissue type were defined as "tissue-type TF-TF interaction networks" (Figure 1c). Finally, a two-step analysis of TFBS conservation and enrichment of overlapping TF target orthologous genes was performed to predict interactions of 3 TFs (Figure 1b).

Identification of tissue TF pairs

Using the function conservation approach [22] and tissue-expressed genes in each of the 79 human tissues, we were able to identify a few hundred TF pairs for each tissue, for which BM-CD71+early erythroid has the largest number of 383 TF pairs, and the ovary tissue has the smallest number of at 230. We also identified 647 TF pairs from housekeeping genes. Filtering TF pairs of housekeeping genes from those in each tissue has greatly reduced the number of TF pairs in each tissue, ranging from 39% TF pairs for lymph node to 59% TF pairs for BM-CD105+endothelial, indicating that a large portion of the TF pairs performs ubiquitous roles across different tissues. The resulting tissue TF pairs range from 111 to 176 for different tissues. We also searched for TF pairs specific to one tissue (referred to tissue-unique TF pairs) and obtained from 2 to 20 such TF pairs for different tissues. The number of tissue TF pairs and tissue-unique TF pairs for each of the 79 human tissues are summarized in Table 1. The top 5 tissue TF pairs that have the most significant correlations between enriched TFBS pairs and enriched overlapping orthologous genes are also listed for each tissue.

Table 1 Summary of the identified tissue and tissue-unique TF pairs as well as top 5 tissue TF pairs in the 79 human tissues.

Overall, we identified 2549 tissue and 803 tissue-unique TF pairs for the 79 human tissues. These results indicate that tissue gene expression is regulated by large sets of interacting TFs. Furthermore, the relative small number of tissue-unique TF pairs out of all tissue TF pairs suggests that identical tissue TF pairs in different tissues may play different functional roles, which prompted us to investigate their biological function. For this purpose, we used Gene2go http://www.ncbi.nlm.nih.gov/ to annotate human genes whose promoters contained the target TFBS pairs, as TFs control cellular biological processes via transcriptional regulation of groups of genes with similar functions. Significant (q-value < 0.1) biological processes for tissue TF pairs were obtained by comparing the number of TF target genes involved in a particular biological process to the number of genes for the same biological process in the whole human genome (Fisher's exact test; p = 2.4 × 10-4 to 8 × 10-28). All tissue and tissue-unique TF pairs as well as their potential biological functions are listed in Additional File 2.

Evaluation by known interacting TFs

Although the function conservation approach has been proven to be a successful means for predicting interacting TFs [22], we sought to assess the validity of the identified tissue TF pairs. We first used TRANSCompel® 10.4 [25] to determine if known interacting TFs were statistically enriched in the predicted tissue TF pairs. The TRANSCompel® database contains 180 experimentally proven composite elements of two or more binding sites which were previously identified by individual wet lab studies from others. Of the180 composite elements, 105 were mapped to the 23,005 (214*215/2) possible combinations of 2 TFs from the 214 non-redundant position weight matrices (PWMs). We first investigated the statistical significance for the occurrence of known interacting TFs in both predicted TF pairs (before filtering) and tissue TF pairs in each of the 79 tissues. Figure 2 shows that known interacting TFs display enrichment more in the tissue TF pairs then in predicted TF pairs for both the number of tissues (37 vs. 9) and the degree of enrichment (Binomial test: p = 3.2 × 10-2 to 6.8 × 10-6 vs. p = 4 × 10-2 to 3.4 × 10-4). We also computed the occurrence of known interacting TFs in all predicted TFs pairs (before filtering) and all tissue TF pairs from the 79 tissues. We found that 40 (38.1%) of the 105 known interacting TFs were in both predicted TF pairs (Binomial test; p = 6.4 × 10-9) and tissue TF pairs (Binomial test; p = 5.4 × 10-11).

Figure 2
figure 2

Statistical significance for the occurrence of known interacting TFs. Validation results (-log10(p-values)) for predicted TF pairs (orange) and tissue TF pairs (green) in the 79 human tissues. Dash line indicates significant p-value cutoff which is < 0.05 above the dash line. Known interacting TFs display enrichment more in the tissue TF pairs than in predicted TF pairs for both the number of tissues (37 vs. 9) and the degree of enrichment (Binomial test: p = 3.2 × 10-2 to 6.8 × 10-6 vs. p = 4 × 10-2 to 3.4 × 10-4).

To further verify our prediction, we next compared the tissue TF pairs to known tissue-specific TFs from liver, for which the cis-regulatory systems for both individual TF binding and synergistic actions have been thoroughly studied [3, 4, 25, 26]. These studies found 40 liver-specific single TFs and 27 liver-specific interacting TFs. We first computed the tissue TF pairs whose two TFs were all liver-specific based on the 820 (40*41/2) possible combinations of 2 TFs from the 40 liver-specific single TFs. Out of 162 tissue TF pairs from liver tissue, we were able to obtain 30 (Binomial test; p = 2.3 × 10-14) where both TFs were liver-specific. For the 27 liver-specific known interacting TFs, we found 8 (30%) in both the predicted TF pairs (Binomial test; p = 3.6 × 10-9) and tissue TF pairs (Binomial test; p = 2.9 × 10-11) from liver tissue. These include HNF4ALPHA:HNF4ALPHA, NF1:COUP_DR1, CEBPGAMMA:HNF4ALPHA, CEBPA:HNF3B, HNF3:HNF4ALPHA, HNF3:PPARA, CEPBA:GATA4, and HNF1:OCT1. All of them are key elements in liver specific transcriptional regulation. GO enrichment analyses indicated that genes whose promoters contained the predicted liver-specific TFBS pairs were mainly involved in liver specific functions [27, 28], including oxidation reduction, acute-phase response, gluconegnesis, and lipoprotein & lipid metabolic processes (Additional File 2). Further analysis of the binding sites on the promoter sequence of individual genes indicated that we were able to reliably identify interacting TFs similar to those previously reported. One of the examples was the APOA1 gene, which was well-characterized to be synergistically bound by HNF3 and HNF4 [8]. Our prediction was able to successfully identify the HNF3:HNF4 binding sites on its promoter. A closer examination shows that our predicted HNF3 and HNF4 binding sites for the APOA1 gene are exactly those experimentally proven, liver tissue-specific HNF3 and HNF4 binding site combinations [8], which are highly conserved between human and mouse in regards to both nucleotide sequence and spacing between each binding site (Figure 3).

Figure 3
figure 3

Conservation of identified HNF3 and HNF4ALPHA binding sites in human and mouse APOA1 genes. Both schematic and sequence alignments for the predicted HNF3 and HNF4ALPHA binding sites between human and mouse promoter sequences are depicted. In the sequence alignment the core motifs are shown in upper case letters and the distances between adjacent binding sites are shown in brackets. Also shown are the locations of each binding site in relation to the transcriptional starting site.

It is important to note that the 79 human tissues represent only part of the temporal and spatial conditions from which the 105 known interacting TFs were discovered, and therefore it is unlikely to have all known interacting TFs in our predicted list. Nevertheless, our results indicate that the use of function conservation approach and tissue-expressed genes was able to reliably identify to a great extent known interacting TFs, thus presenting very strong evidence for the validity of the identified tissue TF pairs. These results also indicate that filtering the TF pairs of housekeeping genes from those in each tissue is an important step to eliminate TFs playing a ubiquitous role, thereby the resulting TF pairs are more tissue-specific.

Identification of tissue-type TF pairs

One of the goals of this study is to find interacting TFs controlling gene expression in a broad spatial and temporal manner such as interacting TFs common to the same type of tissues. This can be achieved by searching tissue TF pairs common across all tissues of the same type such as the 7 muscle tissues. However, the use of all tissues may reduce the power for tissue-type TF pair identification, since the contents of tissue TF pairs and even the function of a common tissue TF pair may be different between tissues of the same type. Therefore, we sought to first classify tissues into smaller but more closely related groups based on tissue TF pairs, from which representative tissues for the same tissue type could be obtained. Accordingly, we used hierarchical clustering to group tissues, as no a priori knowledge was available for the number of groups for each tissue type. The results are shown in Figure 4, where tissues of the same type are generally grouped together such as testis, liver, pancreas, and brain. There are, however, exceptions for other tissue types which are grouped into either distinct groups or into groups with other types of tissues such as muscle and immune systems. While the muscle tissues are classified into two distinct groups, of which one contains skeletal muscle, heart, and tongue and the other contains smooth muscle and cardiac myocytes, tissues for immune systems are classified into a few groups, one of which displays tighten link with a few cancer tissues.

Figure 4
figure 4

Hierarchical clustering over tissue TF pairs from the 79 human tissues. The distance matrix was built using the "binary" method, and hierarchical clustering was performed using the "complete" agglomeration method. Arrows and numbers indicate the selected tissue groups for further analysis.

We extended our analysis to investigate conservation for tissue TF pairs between tissues of the two muscle groups. We computed overlap for both tissue TF pairs and their biological functions between tissues using hypergeometric distribution. We found little or no overlap for both tissue TF pairs and their functions among tissues between these two groups, which was especially true for the function of tissue TF pairs (data not shown). On the other hand, both tissue TF pairs and their functions showed significant overlap between tissues within the same group (Figure 5). These results not only demonstrate the validity of our tissue classification but also indicated that tissues from the same type (here 5 muscle tissues) may have great difference in both the contents of tissue TF pairs and TF functional roles.

Figure 5
figure 5

Overlapping tissue TF pairs and their functions between muscle tissues. (a) Venn diagram displaying the significant overlap for both tissue TF pairs and their functions in the group of skeletal muscle. (b) Significant overlap for both tissue TF pairs and their functions in the group of smooth muscle. Each circle indicates the number of tissue TF pairs. The number of overlapping tissue TF pairs and TF functions between two tissues is indicated in bold (# function/# TF). Also shown are their corresponding p-values from hypergeometric tests.

Based on the clustering results, we selected 11 tissue-type groups, each having 2 to 16 tissues, for tissue-type TF pair discovery. A cutoff threshold of tissue TF pairs common in at least 50% tissues from the same group was set up for searching tissue-type TF pairs. In addition to the TF level, we also searched for tissue-type TF pairs based on their function using the same criteria of > 50% tissues in the same group. To this end, we were able to identify tissue-type pairs for all tissue groups as listed in Table 2. Whereas the number of tissue-type TF pairs ranges from 17 for immune/cancer group to 74 for testis, those at the functional level have relatively smaller numbers, ranging from 3 for thyroid to 40 for testis. All (379) tissue-type TF pairs as well as their corresponding functions for the 11 tissue-type groups are listed in Additional File 3.

Table 2 Number of tissue-type TF pairs in the selected 11 tissue groups.

Reconstruction of tissue-type TF-TF interaction networks

In an effort to reveal TF relationships in controlling tissue gene expression, we performed analysis to reconstruct TF-TF interaction networks. Using tissue-type TF pairs, we first looked for those with one shared TF between each other in the same tissue, from which TF-TF interaction networks were built by joining 2 or more TF pairs (Figure 1c). TF-TF interaction networks with the same topology in at least 2 tissues from the same tissue-type group were then selected as tissue-type TF-TF interaction networks, which are multi-input network motifs consisting of at least 3 TFs that bind to a set of gene promoters. A total of 84 tissue-type TF-TF interaction networks were identified for the 11 tissue-type groups, ranging from 1 for immune/cancer to 22 for testis (Additional File 4). Sixty two of these tissue-type TF-TF interaction networks have a linear relationship between TFs with 1 to 4 internal TFs (i.e. TF connecting to 2 other TFs), indicating that the majority of the TF-TF regulatory networks have simple TF relationships for controlling tissue gene expression. Figure 6a shows a multi-input network motif from liver tissues, in which FOXJ2, HNF1, and TTF1 regulate 6 genes in a combinatorial manner by either 2 or 3 TFs. The remaining 22 tissue-type TF-TF interaction networks display more complex interacting structures with some of the internal TFs connecting to 3 or more TFs.

Figure 6
figure 6

Tissue-type TF-TF interaction networks for liver tissues. (a) Multi-input network motif for TFs of FOXJ2, HNF1, and TTF1 and their target genes which are regulated by either two or three TFs. (b) Topology of seven tissue-type TF-TF interaction networks for liver tissues. Each node represents a TF and each edge (link) represents a significant synergy between two TFs from tissue TF pairs. Previously known liver-specific TFs are labeled with an asterisk.

Unlike the tissue TF pairs, we did not find any common tissue-type TF-TF interaction networks between different tissue types. In light of this, we performed a search to see if any single TFs played central roles in controlling tissue gene expression across different tissues, and looked for internal TFs in multiple tissue-type TF-TF interaction networks. To this end, we found that AP2, PPARA, PAX4, FAC1, ZIC3, and SPZ1 served as internal TFs in 8, 8, 8, 6, 5, and 4 tissue-type TF-TF interaction networks, respectively, suggesting their role as central hubs in tissue-type TF-TF interaction networks. Whereas FAC1 acts as the internal TF in 6 tissue-type TF-TF interaction networks from immune systems and cancer, SPZ1 mainly serves as the internal TF in tissue-type TF-TF interaction networks from testis, and the rest in 5 to 6 tissue-type TF-TF interaction networks from different tissue types. These results indicate that FAC1, when serving as the internal TF, is restricted to the two related tissue types, and that SPZ1, a bHLH-Zip protein, has an important role in testis [29, 30]. The rest have more diversified roles for coordinating network TFs in controlling tissue gene expression.

It is interesting to note that no single TFs serve as the central hub for tissue-type TF-TF interaction networks from liver tissue. However, we observed that 6 of 7 tissue-type TF-TF interaction networks had at least one known liver-specific TF serving as the internal TF as shown in Figure 6b. To investigate if this distribution pattern of liver-specific TFs in the TF-TF interaction networks had any biological meaning, we randomly sampled TFs from the 214 PWMs to build TF-TF networks, each having the same size and order as the real TF-TF interaction networks. The simulated TF-TF networks were then compared to tissue-type TF-TF interaction networks to estimate the statistical significance for the distribution of liver-specific TFs. The results indicate that known liver-specific TFs were significantly enriched as internal TFs for these 7 tissue-type TF-TF interaction networks (bootstrap analysis; p < 10-20). By contrast, the total number of liver-specific TFs in these 7 tissue-type TF-TF interaction networks was not enriched (bootstrap analysis; p = 0.11). These results suggest that liver-specific TFs, other than initiating liver-specific transcriptional event, may play an important role in recruiting non-liver-specific TFs to the tissue-type TF-TF interaction network, thus offering the potential for coordinating and controlling gene expression across a variety of conditions.

Prediction of multiple interacting TFs

As a first step to elucidate cis-regulatory modules involved in tissue gene regulation, we extended our analysis to the interactions of 3 TFs (named as multiple interacting TFs). Using tissue TF pairs from each of the 79 tissues, we performed a two-step analysis of TFBS conservation and enrichment of overlapping orthologous genes between human and mouse (see Methods). Although it is likely that multiple interacting TFs may be under estimation by the use of tissue TF pairs instead of all predicted TF pairs, the predicted multiple interacting TFs are tissue-specific. Therefore, these predictions most likely represent cis-regulatory modules involved in tissue gene regulation. To this end, we identified 1735 unique interactions of 3 TFs for the 79 human tissues, ranging from 9 multiple interacting TFs for testis interstitial to 72 multiple interacting TFs for caudate nucleus (Additional File 5).

The validity of these predicted multiple interacting TFs was assessed by using liver-specific single TFs from TRANSCFAC11.4 [25], as few known cis-regulatory modules were available. We performed analysis to see if known liver-specific TFs were statistically enriched in 30 predicted multiple interacting TFs from liver tissue. We found 4 of them (bootstrap analysis; p < 10-3) whose 3 TFs were all liver-specific, 18 (bootstrap analysis; p < 10-8) with at least 2 liver-specific TFs, and 28 (bootstrap analysis; p < 10-5) with at least 1 liver-specific TF. These results provide evidence for the enrichment of liver-specific TFs in the predicted multiple interacting TFs, which in turn demonstrated the validity of the prediction.

We next searched for all predicted multiple interacting TFs and their potential functions that are common between tissues. The results indicated that, although common multiple interacting TFs existed between most tissues, the highest overlap was within brain tissues and between brain and gland tissues. By contrast, there was little overlap for the functions of multiple interacting TFs, except within brain and cancer and between these 2 tissue types (Additional File 6). The latter is especially interesting to us, as cancer cells have a global effect on immune systems, which in turn control and shape developing cancer [31]. Six multiple interacting TFs were found to have common functions between immune systems and cancer tissues, including CEBPGAMMA:NKX25:PLZF, CEBPGAMMA:PAX4:PLZF, CP2:NFY:PAX4, FOXJ2:PAX4:POU3F2, CEBPGAMMA:PAX4:PLZF, and FOXJ2:HNF3:PAX4. These results revealed not only the common mechanisms for transcriptional regulation but also the common functional role of multiple interacting TFs between cancer and immune systems, including cell cycle, cell division, DNA replication, mitosis, phosphoinositide-mediated signaling, and immune response. These findings therefore provide new insight into the molecular interplay between cancer and immune systems.

Discussion

Tissue gene expression is generally regulated by multiple transcription factors. A major first step toward understanding how tissues achieve their specificity is to identify interacting TFs regulating gene expression in different tissues. Previous computational approaches to predict interacting TFs were mainly based on recognizable sequence features of tissue-specific [19–21] genes derived from genome-wide gene-expression profiling. Despite these studies, the mechanisms controlling tissue gene expression are still not fully understood.

In this study, we utilized our previously developed function conservation approach, which, based on this and a prior study [22], was shown to successfully predict interacting TFs from tissue-expressed genes. Based on the predictions, tissue TF pairs were identified. The advantage of our approach lies in the fact that it does not depend solely on sequence features of genes but rather function conservation of interacting TFs from both their binding sites and putative target genes between closely related species. Other than function conservation, the use of tissue-expressed genes would allow one to avoid the elimination of common genes contributing to tissue development and differentiation between tissues, especially for these closely related tissues (e.g. skeletal muscle and heart) when compared to the use of tissue-specific genes [19–21] which are expressed in a particular tissue. Therefore, the utilization of our function conservation approach and tissue-expressed genes provides an alternative way for tissue interacting TF discovery.

One of the findings of our study indicates that tissue gene expression is controlled by large sets of tissue TF pairs, which is in agreement with previously reported findings from an approach using sequence features of tissue-specific genes by Yu et al. [19]. We were curious to know the differences of interacting TFs identified by the two different approaches, and selected the liver tissue for comparison. For the 8 known liver-specific interacting TFs that were successfully predicted by our approach in the 162 liver tissue TF pairs, we found that HNF3:HNF4ALPHA was in the liver-specific TF pairs predicted by Yu et al. However, we did not find the other 7 known liver-specific interacting TFs predicted in our 162 tissue TF pairs from Yu et al. On the other hand, 6 of the 27 known liver-specific interacting TFs were correctly predicted by Yu et al but were not in our tissue TF pairs from liver tissue. A closer examination shows that liver-tissue TF pairs from our prediction are enriched with CEBP, HNF3, and HNF4, and that liver-specific TF pairs from Yu et al are enriched with HNF1 and HNF4. All these TFs are known liver-specific TFs such as HNF3 [32], which initiates the liver transcriptional event, and HNF1 [33], which interacts with other important TFs to establish transcriptional hierarchy in liver tissues. These results demonstrate that different methods were able to identify interacting TFs from different angles. Therefore, the findings from our study provide new insight into the mechanism controlling tissue gene expression.

Filtering TF pairs of housekeeping genes from those of tissue-expressed genes is an important step to eliminate TF pairs which play general but not tissue-specific roles in individual tissues. The filtering process reduced the number of predicted TF pairs from 3024 to 2549 (15.7%) for all 79 tissues. This reduction for TF pairs was, however, significantly larger when individual tissues were concerned (39% to 59%), indicating that a large number of overlapping TF pairs had ubiquitous roles among different tissues. The remaining interacting TFs in each tissue were more tissue-specific, which was best evidenced by the result that the predicted TF pairs from liver tissue contained the same number of known liver-specific interacting TFs before and after the filtering. The relative small number of tissue-unique TF pairs out of all tissue TF pairs and the findings from conservation analysis for the functions of tissue TF pairs between tissues of two muscle groups from this study also indicate that tissue TF pair with identical 2-TF combination might play different functional roles in different tissues.

Our findings show that tissue gene expression is controlled by a variety of interacting TFs either on the promoter of a gene or through TF-TF interaction networks. These identified TF interactions may constitute a large part of interacting TFs in each tissue but is not a complete list. To fully understand the mechanisms controlling tissue gene expression requires additional study, which has been best evidenced from the comparison of interacting TFs in liver tissue between Yu et al. [19] and ours. Other than the prediction methods, the target gene selection can contribute greatly to tissue TF identification. Our prediction picked up 8 of the 27 known liver-specific interacting TFs in liver tissues. A couple factors might be responsible for not identifying the other known liver-specific interacting TFs. First, these known liver-specific TF interactions were discovered from broad spatial and temporal conditions. The selected liver genes in this study however represented only one of many conditions under which liver-specific TFs play their roles. This was exemplified by known liver-specific interacting TFs in tissue TF pairs from liver and fetal liver tissues from our prediction. Whereas tissue TF pairs from liver tissue contained 8 known live-specific interacting TFs, fetal liver contained 3 known live-specific interacting TFs with 2 common to those in liver, demonstrating the impact of temporal conditions on tissue TF discoveries. Second, it is unlikely for the top 300 tissue-expressed genes from a single condition to all have information for tissue interacting TF prediction. The choice of the top 300 tissue-expressed genes was based on the report of Pennacchio et al. [18] who have successfully used them to predict tissue-specific enhancers. Increasing the size of genes however would increase the chance of bringing noise to the prediction. Therefore, other than different computational approaches, selecting a proper list of tissue-expressed genes would have a great impact on the prediction of tissue TF pairs.

One of the goals of this study was to find interacting TFs controlling tissue gene expression in a broad spatial and temporal manner. We performed analysis to identify tissue-type TF pairs for 11 selected tissue-type groups. While, as described above, each specific tissue may reflect only a small portion of all spatial and temporal conditions where tissue TF pairs play their regulation roles, tissue-type TF interactions provide a general view of their roles in multiple conditions. The analysis process has also led to other findings that the same type of tissues may have significant differences in both the contents of tissue TF pairs and the TF functional roles, which has been demonstrated by the conservation analysis of tissue TF pairs and their functions from muscle tissues. Tissue-type TF-TF interaction networks have provided not only lines of information on how tissue transcriptional programs are constructed but also new findings of potential roles for tissue-specific TFs in TF-TF interaction networks from liver tissue.

Conclusions

In this study, we successfully employed our previously developed function conservation approach [22], to predict functional TF pairs from tissue-expressed genes in 79 human tissues. Based on the predictions, tissue TF pairs were identified. Our analyses led to the discovery of 2549 unique tissue TF pairs for the 79 human tissues. The validity of the discovered tissue TF pairs has been demonstrated by both known interacting and liver-specific TFs. We also extended our study to find interacting TFs controlling gene expression in a broad temporal and spatial manner and identified 379 tissue-type TF pairs from 11 tissue-type groups, from which tissue-type TF-TF interaction networks have been built. The results also indicated that tissue-specific TFs may play an important role in recruiting non-tissue-specific TFs to the TF-TF interaction network, offering the potential for coordinating and controlling tissue gene expression across a variety of conditions. In summary, our findings have shown that tissue gene expression is regulated by large sets of interacting TFs either on the same promoter of a gene or through TF-TF interaction networks.

Methods

Promoter sequences for housekeeping and tissue-expressed genes

The GNF Atlas2 gene expression database (gnfAtlas2) [34], which contains gene expression data from 79 human tissues, was used for the selection of genes. Based on the report of Pennacchio et al[18] we selected in each tissue the top 300 expressed genes (referred to tissue-expressed genes), which have been used and proven to successfully predict tissue-specific enhancers. Housekeeping genes are the 1018 genes defined by Farre et al[24]. Redundant genes in each group were first removed. Although regulatory elements can exist anywhere in the genome, they are more concentrated around the transcriptional start sites [35]. To reduce false predictions we focused on the proximal promoters which have been proven to successfully predict tissue-specific regulatory elements [21, 36]. It is however worthy to note that the use of 1 kb promoter sequences has limitation for the prediction of tissue TF pairs when compared to the experimental approaches such as ChIP-chip experiment, in which TF pairs can be detected anywhere in the genome. Considering no benchmark promoter sequence dataset is currently available for computational prediction of functional TF pairs, the use of 1 kb promoter sequences and our computational approach nevertheless provide an alternative way for tissue interacting TF discovery. Promoter sequences within 1 kb upstream of transcriptional starting sites for both human and corresponding mouse orthologous genes were extracted from the UCSC genome browser (hg18 March 2006 assembly, mm9 July 2007 assembly). Orthologous genes with promoter sequences from both human and mouse were selected for further analysis. This procedure resulted in 208 to 278 orthologous promoter sequences for tissue-expressed genes and 986 orthologous promoter sequences for housekeeping genes.

Prediction of TF pairs and tissue TF pairs

The procedures for predicting TF pair are basically the same as previously described [22] (Additional File 1). Briefly, background sequences were created by shuffling the DNA sequences within each promoter by either mixing completely or keeping dinucleotides together. These background sequences are preferable to using intergenic sequences which usually are AT-rich or exonic sequences whose nucleotide distributions tend to be biased, when compared to the test promoter sequences. The resulting shuffled sequences from human and mouse, together with the original promoter sequences, were employed for TFBS detection using the Match® program [37], for which the profile parameter was set to "minimize the sum of false positives and negatives", and the 214 non-redundant vertebrate PWMs from the professional TRANSFAC11.4 database [25]. To detect enriched TF pairs out of 23,005 (214*215/2) possible combinations of 2 TFs, distance constraints were first applied for the selection of co-occurring TFBSs with a defined maximum distance between 2 TFBSs. A total of 10 distances were defined, ranging from the smallest 20 bp to the largest 200 bp with a 20 bp increment. The assumption behind the distance constraint is that functional TFBS pairs are more distance-restricted than random co-occurrence of TFBSs [19, 38]. This is true not only in human, for which we found that functional TF pairs were enriched within 200 bp distance ranges [22], but also in Drosophila, in which short-range linkages (< 50 bp) between TFs was overrepresented but mid-range distances (100-500 bp) between TFs was depleted [39]. Enrichment of TFBS pairs for each distance constraint was achieved by computing the ratio of counts for a particular TFBS pair in real promoter sequences vs. the counts of the same TFBS pair in background sequences. To reduce noise while keeping as many as TFBS pairs for the integration of function conservation analysis described below, TFBS pairs with ratio > 1 in more than 5 distance constraints were selected.

A two-step analysis procedure was employed to compute the enrichment of overlapping orthologous genes for a particular TFBS pair. First, a cutoff threshold of at least 10% overlapping orthologous genes between mouse and human was set up for selecting genes whose promoters contained the TFBS pair. The enrichment of overlapping orthologous genes was then estimated by computing the ratio of overlapping orthologous genes from real promoter sequences against those from shuffled sequences. This analysis was performed for each distance constraint. The integration of function conservation for each TF pair was achieved by estimating the correlation (Pearson correlation coefficients) between the 10 enriched TFBS pairs and 10 corresponding enriched overlapping orthologous genes from the same distance constraint. Permutation tests were employed to estimate the statistical significance of correlation by randomly matching the 10 TFBS pair ratios with the 10 overlapping orthologous gene ratios. For multiple test correction, a cutoff threshold of q-value < 0.05 was applied. TF pairs are those passing the cutoff and common between human and mouse.

We next filtered TF pairs of housekeeping genes from those in each tissue (Figure 1a). This was done by removing TF pairs in a particular tissue common to those from housekeeping genes. The remaining TF pairs in each tissue were more tissue-specific, and therefore, were defined as tissue TF pairs. Similar results were obtained from using background sequences of either completely mixed nucleotides or keeping dinucleotides together or completely mixing nucleotides. The results from completely mixed nucleotides were used.

Clustering analysis

To group tissues based on their tissue TF pairs, a 2549 (tissue TF pairs) × 79 (tissues) matrix with binary numbers was first built for all tissue TF pairs in the 79 human tissues. The presence of a tissue TF pair in the matrix was labeled with 1 and the absence was labeled with 0. A distance matrix was then built using the "binary" method, and hierarchical clustering was subsequently performed using the "complete" agglomeration method. All analysis was performed using the R statistical package [40].

Predicting multiple interacting TFs

A two-step analysis of TFBS conservation and enrichment of overlapping orthologous genes was performed to predict interactions of 3 TFs. For TFBS conservation, the identified tissue TFBS pairs were first used to construct all possible 3-TFBS combinations by searching paired tissue TFBS pairs with one shared TFBS between each other on exactly the same location of a gene's promoter (Figure 1b) in a particular tissue. Orthologous gene pairs containing conserved 3-TFBS combination between human and mouse were then selected. Conserved 3-TFBS combinations are those whose 3 TFBSs have the same order and orientation on the promoter sequences between human and mouse orthologous genes. For enrichment of overlapping orthologous genes in a tissue, however, multiple interacting TFs from different orthologous gene pairs were considered to be the same as long as they contained the same 3 TFs. Enriched multiple interacting TFs are those with 3-TFBS combinations occurring on at least 10 orthologous gene promoters and with their target orthologous genes displaying significant overlap between human and mouse (p = 3 × 10-2 to < 10-36 and q < 0.05).

Statistical methods for enrichment analyses

Two main statistical methods were employed for estimating the significance of enrichment in this study. For validating predicted TF pairs and tissue TF pairs by known interacting TFs, the binomial distribution probability, as shown below, was used to determine if known interacting TFs were present more often in the predicted TF pairs or tissue TF pairs than in a randomly selected group from a given list of TFs.

For example, in the case of estimating the statistical significance of known liver-specific interacting TFs in our predicted tissue TF pairs from liver tissue, the n is the number of known liver-specific interacting TFs in the predicted tissue TF pairs from this study; N the number of tissue TF pairs from liver tissue; and p f the background probability of liver-specific TF pairs in all possible combinations of 2 TFs from 214 PWMs.

The statistical significance was computed using the hypergeometric distribution to estimate (1) the enrichment of overlapping tissue TF pairs and their overlapping functions between muscle tissues, and (2) overlapping orthologous genes in predicting multiple interacting TFs.

In the case of overlapping orthologous genes in predicting multiple interacting TFs, for example, c is the number of orthologous gene pairs containing conserved 3-TFBS combination between human and mouse; N the number of tissue-expressed genes for a particular tissue; S1 and S2 are the numbers of tissue-expressed genes with 3-TFBS combinations corresponding to those in c for human and mouse, respectively. The resulting p-value is the probability of observing c or more orthologous gene pairs containing conserved 3-TFBS combination from two sets of size S1 and S2 drawn from a set of N tissue-expressed genes.

Abbreviations

TF:

transcription factor

TFBS:

transcription factor binding site

PWM:

position weight matrices.

References

  1. Arnone MI, Davidson EH: The hardwiring of development: organization and function of genomic regulatory systems. Development (Cambridge, England). 1997, 124 (10): 1851-1864.

    CAS  Google Scholar 

  2. Odom DT, Dowell RD, Jacobsen ES, Nekludova L, Rolfe PA, Danford TW, Gifford DK, Fraenkel E, Bell GI, Young RA: Core transcriptional regulatory circuitry in human hepatocytes. Molecular systems biology. 2006, 2: 2006 0017-10.1038/msb4100059.

    Article  PubMed Central  PubMed  Google Scholar 

  3. Metzger S, Halaas JL, Breslow JL, Sladek FM: Orphan receptor HNF-4 and bZip protein C/EBP alpha bind to overlapping regions of the apolipoprotein B gene promoter and synergistically activate transcription. The Journal of biological chemistry. 1993, 268 (22): 16831-16838.

    CAS  PubMed  Google Scholar 

  4. Costa RH, Grayson DR: Site-directed mutagenesis of hepatocyte nuclear factor (HNF) binding sites in the mouse transthyretin (TTR) promoter reveal synergistic interactions with its enhancer region. Nucleic Acids Res. 1991, 19 (15): 4139-4145. 10.1093/nar/19.15.4139.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  5. Rada-Iglesias A, Wallerman O, Koch C, Ameur A, Enroth S, Clelland G, Wester K, Wilcox S, Dovey OM, Ellis PD: Binding sites for metabolic disease related transcription factors inferred at base pair resolution by chromatin immunoprecipitation and genomic microarrays. Human molecular genetics. 2005, 14 (22): 3435-3447. 10.1093/hmg/ddi378.

    Article  CAS  PubMed  Google Scholar 

  6. Odom DT, Zizlsperger N, Gordon DB, Bell GW, Rinaldi NJ, Murray HL, Volkert TL, Schreiber J, Rolfe PA, Gifford DK: Control of pancreas and liver gene expression by HNF transcription factors. Science. 2004, 303 (5662): 1378-1381. 10.1126/science.1089769.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  7. Rubins NE, Friedman JR, Le PP, Zhang L, Brestelli J, Kaestner KH: Transcriptional networks in the liver: hepatocyte nuclear factor 6 function is largely independent of Foxa2. Molecular and cellular biology. 2005, 25 (16): 7069-7077. 10.1128/MCB.25.16.7069-7077.2005.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  8. Harnish DC, Malik S, Karathanasis SK: Activation of apolipoprotein AI gene transcription by the liver-enriched factor HNF-3. The Journal of biological chemistry. 1994, 269 (45): 28220-28226.

    CAS  PubMed  Google Scholar 

  9. Hernandez-Munain C, Krangel MS: Regulation of the T-cell receptor delta enhancer by functional cooperation between c-Myb and core-binding factors. Molecular and cellular biology. 1994, 14 (1): 473-483.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  10. John S, Reeves RB, Lin JX, Child R, Leiden JM, Thompson CB, Leonard WJ: Regulation of cell-type-specific interleukin-2 receptor alpha-chain gene expression: potential role of physical interactions between Elf-1, HMG-I(Y), and NF-kappa B family proteins. Molecular and cellular biology. 1995, 15 (3): 1786-1796.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  11. Sartorelli V, Webster KA, Kedes L: Muscle-specific expression of the cardiac alpha-actin gene requires MyoD1, CArG-box binding factor, and Sp1. Genes Dev. 1990, 4 (10): 1811-1822. 10.1101/gad.4.10.1811.

    Article  CAS  PubMed  Google Scholar 

  12. Naidu PS, Ludolph DC, To RQ, Hinterberger TJ, Konieczny SF: Myogenin and MEF2 function synergistically to activate the MRF4 promoter during myogenesis. Molecular and cellular biology. 1995, 15 (5): 2707-2718.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  13. Moore ML, Wang GL, Belaguli NS, Schwartz RJ, McMillin JB: GATA-4 and serum response factor regulate transcription of the muscle-specific carnitine palmitoyltransferase I beta in rat heart. The Journal of biological chemistry. 2001, 276 (2): 1026-1033. 10.1074/jbc.M009352200.

    Article  CAS  PubMed  Google Scholar 

  14. Bailey PJ, Klos JM, Andersson E, Karlen M, Kallstrom M, Ponjavic J, Muhr J, Lenhard B, Sandelin A, Ericson J: A global genomic transcriptional code associated with CNS-expressed genes. Experimental cell research. 2006, 312 (16): 3108-3119. 10.1016/j.yexcr.2006.06.017.

    Article  CAS  PubMed  Google Scholar 

  15. Krivan W, Wasserman WW: A predictive model for regulatory sequences directing liver-specific transcription. Genome Res. 2001, 11 (9): 1559-1566. 10.1101/gr.180601.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  16. Chen X, Blanchette M: Prediction of tissue-specific cis-regulatory modules using Bayesian networks and regression trees. BMC bioinformatics. 2007, 8 (Suppl 10): S2-10.1186/1471-2105-8-S10-S2.

    Article  Google Scholar 

  17. Wasserman WW, Fickett JW: Identification of regulatory regions which confer muscle-specific gene expression. Journal of molecular biology. 1998, 278 (1): 167-181. 10.1006/jmbi.1998.1700.

    Article  CAS  PubMed  Google Scholar 

  18. Pennacchio LA, Loots GG, Nobrega MA, Ovcharenko I: Predicting tissue-specific enhancers in the human genome. Genome Res. 2007, 17 (2): 201-211. 10.1101/gr.5972507.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  19. Yu X, Lin J, Zack DJ, Qian J: Computational analysis of tissue-specific combinatorial gene regulation: predicting interaction between transcription factors in human tissues. Nucleic Acids Res. 2006, 34 (17): 4925-4936. 10.1093/nar/gkl595.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  20. Smith AD, Sumazin P, Xuan Z, Zhang MQ: DNA motifs in human and mouse proximal promoters predict tissue-specific expression. Proc Natl Acad Sci USA. 2006, 103 (16): 6275-6280. 10.1073/pnas.0508169103.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  21. Smith AD, Sumazin P, Zhang MQ: Tissue-specific regulatory elements in mammalian promoters. Molecular systems biology. 2007, 3: 73-

    PubMed Central  PubMed  Google Scholar 

  22. Hu Z, Hu B, Collins JF: Prediction of synergistic transcription factors by function conservation. Genome Biol. 2007, 8 (12): R257-10.1186/gb-2007-8-12-r257.

    Article  PubMed Central  PubMed  Google Scholar 

  23. Su AI, Wiltshire T, Batalov S, Lapp H, Ching KA, Block D, Zhang J, Soden R, Hayakawa M, Kreiman G: A gene atlas of the mouse and human protein-encoding transcriptomes. Proc Natl Acad Sci USA. 2004, 101 (16): 6062-6067. 10.1073/pnas.0400782101.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  24. Farre D, Bellora N, Mularoni L, Messeguer X, Alba MM: Housekeeping genes tend to show reduced upstream sequence conservation. Genome Biol. 2007, 8 (7): R140-10.1186/gb-2007-8-7-r140.

    Article  PubMed Central  PubMed  Google Scholar 

  25. Matys V, Kel-Margoulis OV, Fricke E, Liebich I, Land S, Barre-Dirrie A, Reuter I, Chekmenev D, Krull M, Hornischer K: TRANSFAC and its module TRANSCompel: transcriptional gene regulation in eukaryotes. Nucleic Acids Res. 2006, D108-110. 10.1093/nar/gkj143. 34 Database

  26. Tuteja G, Jensen ST, White P, Kaestner KH: Cis-regulatory modules in the mammalian liver: composition depends on strength of Foxa2 consensus site. Nucleic Acids Res. 2008, 36 (12): 4149-4157. 10.1093/nar/gkn366.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  27. Parker GA, Picut CA: Liver immunobiology. Toxicol Pathol. 2005, 33 (1): 52-62. 10.1080/01926230590522365.

    Article  CAS  PubMed  Google Scholar 

  28. Nguyen P, Leray V, Diez M, Serisier S, Le Bloc'h J, Siliart B, Dumon H: Liver lipid metabolism. J Anim Physiol Anim Nutr (Berl). 2008, 92 (3): 272-283. 10.1111/j.1439-0396.2007.00752.x.

    Article  CAS  Google Scholar 

  29. Hsu SH, Shyu HW, Hsieh-Li HM, Li H: Spz1, a novel bHLH-Zip protein is specifically expressed in testis. Mechanisms of development. 2001, 100 (2): 177-187. 10.1016/S0925-4773(00)00513-X.

    Article  CAS  PubMed  Google Scholar 

  30. Hrabchak C, Varmuza S: Identification of the spermatogenic zip protein Spz1 as a putative protein phosphatase-1 (PP1) regulatory protein that specifically binds the PP1cgamma2 splice variant in mouse testis. The Journal of biological chemistry. 2004, 279 (35): 37079-37086. 10.1074/jbc.M403710200.

    Article  CAS  PubMed  Google Scholar 

  31. Finn OJ: Cancer immunology. The New England journal of medicine. 2008, 358 (25): 2704-2715. 10.1056/NEJMra072739.

    Article  CAS  PubMed  Google Scholar 

  32. Friedman JR, Kaestner KH: The Foxa family of transcription factors in development and metabolism. Cell Mol Life Sci. 2006, 63 (19-20): 2317-2328. 10.1007/s00018-006-6095-6.

    Article  CAS  PubMed  Google Scholar 

  33. Zhou DX, Yen TS: The ubiquitous transcription factor Oct-1 and the liver-specific factor HNF-1 are both required to activate transcription of a hepatitis B virus promoter. Molecular and cellular biology. 1991, 11 (3): 1353-1359.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  34. Su AI, Cooke MP, Ching KA, Hakak Y, Walker JR, Wiltshire T, Orth AP, Vega RG, Sapinoso LM, Moqrich A: Large-scale analysis of the human and mouse transcriptomes. Proc Natl Acad Sci USA. 2002, 99 (7): 4465-4470. 10.1073/pnas.012025199.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  35. Cooper SJ, Trinklein ND, Anton ED, Nguyen L, Myers RM: Comprehensive analysis of transcriptional promoter structure and function in 1% of the human genome. Genome Res. 2006, 16 (1): 1-10. 10.1101/gr.4222606.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  36. Martinez MJ, Smith AD, Li B, Zhang MQ, Harrod KS: Computational prediction of novel components of lung transcriptional networks. Bioinformatics. 2007, 23 (1): 21-29. 10.1093/bioinformatics/btl531.

    Article  CAS  PubMed  Google Scholar 

  37. Kel AE, Gossling E, Reuter I, Cheremushkin E, Kel-Margoulis OV, Wingender E: MATCH: A tool for searching transcription factor binding sites in DNA sequences. Nucleic Acids Res. 2003, 31 (13): 3576-3579. 10.1093/nar/gkg585.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  38. Yu X, Lin J, Masuda T, Esumi N, Zack DJ, Qian J: Genome-wide prediction and characterization of interactions between transcription factors in Saccharomyces cerevisiae. Nucleic Acids Res. 2006, 34 (3): 917-927. 10.1093/nar/gkj487.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  39. Papatsenko D, Goltsev Y, Levine M: Organization of developmental enhancers in the Drosophila embryo. Nucleic Acids Res. 2009, 37 (17): 5665-5677. 10.1093/nar/gkp619.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  40. The R Project for Statistical Computing. [http://www.r-project.org/]

  41. Sterneck E, Muller C, Katz S, Leutz A: Autocrine growth induced by kinase type oncogenes in myeloid cells requires AP-1 and NF-M, a myeloid specific C/EBP-like factor. The EMBO journal. 1992, 11 (1): 115-126.

    CAS  PubMed Central  PubMed  Google Scholar 

  42. Ribeiro A, Pastier D, Kardassis D, Chambaz J, Cardot P: Cooperative binding of upstream stimulatory factor and hepatic nuclear factor 4 drives the transcription of the human apolipoprotein A-II gene. The Journal of biological chemistry. 1999, 274 (3): 1216-1225. 10.1074/jbc.274.3.1216.

    Article  CAS  PubMed  Google Scholar 

  43. Tsukada J, Misago M, Serino Y, Ogawa R, Murakami S, Nakanishi M, Tonai S, Kominato Y, Morimoto I, Auron PE: Human T-cell leukemia virus type I Tax transactivates the promoter of human prointerleukin-1beta gene through association with two transcription factors nuclear factor-interleukin-6 and Spi-1. Blood. 1997, 90 (8): 3142-3153.

    CAS  PubMed  Google Scholar 

  44. Janson L, Pettersson U: Cooperative interactions between transcription factors Sp1 and OTF-1. Proc Natl Acad Sci USA. 1990, 87 (12): 4732-4736. 10.1073/pnas.87.12.4732.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  45. Du W, Thanos D, Maniatis T: Mechanisms of transcriptional synergism between distinct virus-inducible enhancer elements. Cell. 1993, 74 (5): 887-898. 10.1016/0092-8674(93)90468-6.

    Article  CAS  PubMed  Google Scholar 

  46. Leger H, Sock E, Renner K, Grummt F, Wegner M: Functional interaction between the POU domain protein Tst-1/Oct-6 and the high-mobility-group protein HMG-I/Y. Molecular and cellular biology. 1995, 15 (7): 3738-3747.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  47. Yoo JY, Wang W, Desiderio S, Nathans D: Synergistic activity of STAT3 and c-Jun at a specific array of DNA elements in the alpha 2-macroglobulin promoter. The Journal of biological chemistry. 2001, 276 (28): 26421-26429. 10.1074/jbc.M009935200.

    Article  CAS  PubMed  Google Scholar 

  48. Zhang X, Wrzeszczynska MH, Horvath CM, Darnell JE: Interacting regions in Stat3 and c-Jun that participate in cooperative transcriptional activation. Molecular and cellular biology. 1999, 19 (10): 7138-7146.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  49. Dwivedi PP, Omdahl JL, Kola I, Hume DA, May BK: Regulation of rat cytochrome P450C24 (CYP24) gene expression. Evidence for functional cooperation of Ras-activated Ets transcription factors with the vitamin D receptor in 1,25-dihydroxyvitamin D(3)-mediated induction. The Journal of biological chemistry. 2000, 275 (1): 47-55. 10.1074/jbc.275.1.47.

    Article  CAS  PubMed  Google Scholar 

  50. Verrijzer CP, van Oosterhout JA, Vliet van der PC: The Oct-1 POU domain mediates interactions between Oct-1 and other POU proteins. Molecular and cellular biology. 1992, 12 (2): 542-551.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  51. Maschek U, Pulm W, Hammerling GJ: Altered regulation of MHC class I genes in different tumor cell lines is reflected by distinct sets of DNase I hypersensitive sites. The EMBO journal. 1989, 8 (8): 2297-2304.

    CAS  PubMed Central  PubMed  Google Scholar 

  52. Tian G, Erman B, Ishii H, Gangopadhyay SS, Sen R: Transcriptional activation by ETS and leucine zipper-containing basic helix-loop-helix proteins. Molecular and cellular biology. 1999, 19 (4): 2946-2957.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  53. Dang W, Sun XH, Sen R: ETS-mediated cooperation between basic helix-loop-helix motifs of the immunoglobulin mu heavy-chain gene enhancer. Molecular and cellular biology. 1998, 18 (3): 1477-1488.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  54. Magne S, Caron S, Charon M, Rouyez MC, Dusanter-Fourt I: STAT5 and Oct-1 form a stable complex that modulates cyclin D1 expression. Molecular and cellular biology. 2003, 23 (24): 8934-8945. 10.1128/MCB.23.24.8934-8945.2003.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  55. Gegonne A, Bosselut R, Bailly RA, Ghysdael J: Synergistic activation of the HTLV1 LTR Ets-responsive region by transcription factors Ets1 and Sp1. The EMBO journal. 1993, 12 (3): 1169-1178.

    CAS  PubMed Central  PubMed  Google Scholar 

  56. Meyer WK, Reichenbach P, Schindler U, Soldaini E, Nabholz M: Interaction of STAT5 dimers on two low affinity binding sites mediates interleukin 2 (IL-2) stimulation of IL-2 receptor alpha gene transcription. The Journal of biological chemistry. 1997, 272 (50): 31821-31828. 10.1074/jbc.272.50.31821.

    Article  CAS  PubMed  Google Scholar 

  57. Yang TT, Chow CW: Transcription cooperation by NFAT.C/EBP composite enhancer complex. The Journal of biological chemistry. 2003, 278 (18): 15874-15885. 10.1074/jbc.M211560200.

    Article  CAS  PubMed  Google Scholar 

  58. Rivera RR, Stuiver MH, Steenbergen R, Murre C: Ets proteins: new factors that regulate immunoglobulin heavy-chain gene expression. Molecular and cellular biology. 1993, 13 (11): 7163-7169.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  59. Shapiro LH: Myb and Ets proteins cooperate to transactivate an early myeloid gene. The Journal of biological chemistry. 1995, 270 (15): 8763-8771.

    Article  CAS  PubMed  Google Scholar 

  60. Bert AG, Burrows J, Hawwari A, Vadas MA, Cockerill PN: Reconstitution of T cell-specific transcription directed by composite NFAT/Oct elements. J Immunol. 2000, 165 (10): 5646-5655.

    Article  CAS  PubMed  Google Scholar 

  61. Duncliffe KN, Bert AG, Vadas MA, Cockerill PN: A T cell-specific enhancer in the interleukin-3 locus is activated cooperatively by Oct and NFAT elements within a DNase I-hypersensitive site. Immunity. 1997, 6 (2): 175-185. 10.1016/S1074-7613(00)80424-0.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The authors would like to thank Dr. Thomas Furlani for critical reading of the manuscript. This work was partially supported by NIH grant AI28304.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Zihua Hu.

Additional information

Authors' contributions

ZH initiated and designed the study, conceived the analysis procedure, carried out data analysis, and wrote the manuscript. SG helped write Perl scripts to process data. Both authors read and approved the final manuscript.

Electronic supplementary material

Additional file 1: Flowchart of analysis procedure for TF pair prediction. (PDF 15 KB)

12864_2009_2643_MOESM2_ESM.XLS

Additional file 2: Lists tissue and tissue-unique TF pairs and their potential functions for the 79 human tissues. (XLS 1 MB)

Additional file 3: Lists tissue-type TF pairs for the 11 selected tissue groups. (XLS 67 KB)

Additional file 4: Shows the 84 tissue-type TF-TF interaction networks from the 11 tissue-type groups. (PDF 97 KB)

12864_2009_2643_MOESM5_ESM.xls

Additional file 5: Lists multiple interacting TFs (3 TFs) and their potential functions for the 79 human tissues. (XLS 284 KB)

12864_2009_2643_MOESM6_ESM.pdf

Additional file 6: Overlap matrix for multiple interacting TFs. The overlap of multiple interacting TFs between the 79 human tissues is depicted in the upper right panel and overlap of function for multiple interacting TFs in the lower left panel. The degree of overlap is indicated by color with red showing the greatest overlap and yellow showing less overlap. (PDF 31 KB)

Authors’ original submitted files for images

Rights and permissions

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.

Reprints and permissions

About this article

Cite this article

Hu, Z., Gallo, S.M. Identification of interacting transcription factors regulating tissue gene expression in human. BMC Genomics 11, 49 (2010). https://doi.org/10.1186/1471-2164-11-49

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-11-49

Keywords