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

Integrating bioinformatic resources to predict transcription factors interacting with cis-sequences conserved in co-regulated genes

Abstract

Background

Using motif detection programs it is fairly straightforward to identify conserved cis-sequences in promoters of co-regulated genes. In contrast, the identification of the transcription factors (TFs) interacting with these cis-sequences is much more elaborate. To facilitate this, we explore the possibility of using several bioinformatic and experimental approaches for TF identification. This starts with the selection of co-regulated gene sets and leads first to the prediction and then to the experimental validation of TFs interacting with cis-sequences conserved in the promoters of these co-regulated genes.

Results

Using the PathoPlant database, 32 up-regulated gene groups were identified with microarray data for drought-responsive gene expression from Arabidopsis thaliana. Application of the binding site estimation suite of tools (BEST) discovered 179 conserved sequence motifs within the corresponding promoters. Using the STAMP web-server, 49 sequence motifs were classified into 7 motif families for which similarities with known cis-regulatory sequences were identified. All motifs were subjected to a footprintDB analysis to predict interacting DNA binding domains from plant TF families. Predictions were confirmed by using a yeast-one-hybrid approach to select interacting TFs belonging to the predicted TF families. TF-DNA interactions were further experimentally validated in yeast and with a Physcomitrella patens transient expression system, leading to the discovery of several novel TF-DNA interactions.

Conclusions

The present work demonstrates the successful integration of several bioinformatic resources with experimental approaches to predict and validate TFs interacting with conserved sequence motifs in co-regulated genes.

Background

In recent years large numbers of novel cis-regulatory sequences have been described that are conserved in stress-response genes, but the role of these elements and their binding transcription factors (TFs) remains unknown [1–3]. The identification of the binding TFs is therefore a major challenge for bioinformaticians and experimentalists. Using database-assisted analysis of cis-sequences it is possible to generate hypotheses on the nature of the binding TF [4], but the experimental validation of these predictions is often missing. The current work investigates how efficient a database-assisted approach leads to the prediction of the correct TF or TF-families that bind to conserved cis-sequences in co-regulated genes.

The identification of cis-regulatory sequences has been facilitated using bioinformatic and web-queryable resources [4–9]. One approach is the detection of cis-elements annotated in existing resources by performing a database-assisted analysis [10–12]. Another approach is the de novo discovery of conserved sequence patterns in sets of co-regulated genes without knowing if these sequences have been associated with a function before. This involves detection of over-represented sequences in the upstream region of co-regulated genes by using pattern mining programs such as MEME, AlignACE, CONSENSUS, Co-Bind, BioProspector, MITRA, or integrative frameworks such as BEST [13–19].

To discover conserved sequence motifs in promoters of up-regulated A. thaliana genes, BEST was applied in the work presented here. Co-regulated genes were identified using drought responsive microarray expression data annotated to the PathoPlant database [20, 21]. The identified motifs were classified with STAMP [22] and compared to known cis-regulatory sequences annotated to the AthaMap, PLACE and AGRIS databases [23–28]. Furthermore, the newly developed footprintDB repository, built on top of 3D-footprint, was employed to predict interacting TFs [29, 30]. To confirm bioinformatic predictions, 15 cis-sequences were used for the isolation of interacting TFs in a yeast one-hybrid screening. The specificity of each TF-DNA interaction was further validated in yeast. Furthermore, two TFs were used for the generation of synthetic factors that activate reporter gene expression under the control of synthetic promoters harbouring corresponding cis-sequences in a P. patens expression system.

Results

Identification of seven motif families in drought-responsive A. thalianagenes

The goal of this work was the identification of transcription factors binding to cis-regulatory sequences conserved in promoters of drought-responsive genes. Thus, the first part aimed to identify sequence motifs harbouring such cis-sequences. Figure 1 gives an overview of the workflow. In a first step, microarray experiments annotated to the PathoPlant database were employed to identify gene sets two- to tenfold up-regulated by drought in roots and shoots. A total of 32 queries were performed. Additional file 1 lists the parameters for these queries and also the number of induced genes obtained in each of the 32 queries. The 32 queries yielded 32 up-regulated gene groups containing up to 34 co-regulated genes. In a second step, conserved sequence motifs within the upstream region of the 32 co-regulated gene groups were identified with the software package BEST [19]. For motif detection, 1,000 bp upstream of all genes in a co-induced gene group were screened with this software. This analysis resulted in 179 sequence motifs, 22 from roots, 46 from shoots, and 111 from combinations of timepoints or tissues. In a third step, all sequences from the 179 motifs were subjected to an in silico expression analysis [31]. With this tool, motifs specific for drought responsive genes were identified. A total of 49 motifs containing cis-sequences that pass this analysis remained, 5 from roots, 14 from shoots and 30 from combinations of timepoints or tissues. These sequences and the alignments generating the 49 sequence motifs are shown in Additional file 2. Single sequences that pass the in silico expression analysis criteria are shown in bold in this file.

Figure 1
figure 1

Workflow to identify sequence motifs conserved in drought responsive genes.

These 49 sequence motifs were further analysed using STAMP (Step 4, Figure 1). STAMP determines the relationship between the motifs and also calculates similarities to known cis-regulatory sequences [22]. To illustrate the similarities among the 49 sequence motifs, a tree was generated and illustrated by the program MEGA [32], which is shown in Figure 2. Based on the relationship between the motifs, 7 motif groups were designated. STAMP was also employed for the identification of motif similarities by comparing the 49 motifs with known cis-elements from plant databases AthaMap, AGRIS and PLACE [23–25]. Additional file 3 summarises these results. Several motifs from group I, that harbours 29 motifs, show significant similarities (low E-value) to a G-box motif and to abscisic acid-responsive elements (ABRE) known as bZIP binding sites. For example with the PLACE database eight motifs show similarities with the GBOXLERBCS sequence [33] and nine motifs show similarities to four different ABREs from PLACE (ABREATRD22, ABREAZMRAB28, ABRETAEM, ABREMOTIFIIIOSRAB16B). These ABREs are associated with abiotic stress response [34–37]. Twenty motifs cluster within the remaining groups (II to VII) and display less significant similarities (higher E-value) to a wide array of cis-elements, with many of them so far not being associated with abiotic stress response (Additional file 3).

Figure 2
figure 2

Relationship tree for 49 sequence motifs identified bioinformatically in co-regulated gene groups.

Bioinformatic prediction of TFs interacting with conserved cis-sequences

The collection of sequence motifs detected with BEST were further analysed with footprintDB, which currently contains over 5,000 unique TFs and their experimentally-associated DNA motifs [30]. The flowchart in Figure 3 illustrates this approach. This strategy retrieved A. thaliana candidates for all motif groups but group V. These predictions further confirm that 29 group I motifs display obvious similarities to G-box, ABRE-binding motifs recognised by TFs of the ABF family (with bZIP domains), and in general reproduce the STAMP analysis explained earlier. The remaining groups do not have a clear resemblance to any documented DNA binding specificity, and thus proteins from a variety of families were retrieved. In total, 865 unique A. thaliana protein sequences were retrieved for the 49 drought-related motifs (Additional file 4, sheet [A.th.homologues]).

Figure 3
figure 3

Workflow to predict binding TFs for 49 cis -elements/sequence motifs.

Isolation of TFs interacting with conserved cis-sequences

Sequences from the identified DNA motifs, which are conserved in the promoter of genes simultaneously induced by drought were cloned (six repeats) into the yeast pHISi vector with the aim to carry out yeast one-hybrid (Y1H) screens, in order to identify TFs that are able to activate the transcription from these putative cis-regulatory sequences. Out of these 49 sequence motifs, 15 individual cis-sequences shown in Table 1 and representing all seven motifs groups, were used for Y1H. Additional file 5 shows the result of the Y1H screen. A total of 49 different TFs belonging to 12 different families were identified. Amongst the identified factors, 11 were already described as playing a role in various stress responses. Five were found to be associated with abiotic stresses (e.g. drought, salt), three with biotic stresses (e.g. pathogens), and three with both types of stress (Additional file 5). Moreover, this analysis led in most of the cases to the identification of specific associations between a motif group and a TF family. Table 2 and Figure 4A summarise these findings. Table 3 shows all selected TFs that were able to activate the screened sequence but unable to activate the mutated sequence in yeast (Figure 4B, C, D, and E, left panel). Sequence 2 from motif group I and sequence 10 from motif group V preferably select bHLH TFs (Table 2). Instead, both sequences from motif group II preferably select NAC factors. Interestingly, sequences from motif group III and VI, which are closely related (Figure 2), preferably select MYB TFs. For the other sequences, no preferred TF family members were selected (Table 2).

Table 1 Cis -sequences selected for further experimental and bioinformatics analysis
Table 2 Yeast One-Hybrid Screening with cis -sequences
Figure 4
figure 4

Identification and functional characterization of TFs interacting with conserved cis -sequences. (A) Simplified similarity tree (based on STAMP analysis) displaying the main associations observed between the different groups of conserved cis-sequences and the different class of TFs identified in yeast one-hybrid (Y1H) experiments. (B-E) Transcriptional activity and interaction of selected TFs with sequences 3 (B), 7 (C), 5 (D) and 9 (E) in Y1H experiments (left panels). Six repeats of the tested sequence, or their mutated versions (bold characters indicate the mutated nucleotides), fused to the HIS3 auxotrophic gene were tested. Upper part, growth on control media deprived of W amino acid, allowing the selection of yeast that express the studied TFs (i.e. NAC, NAC42, PIF7, MYB61 and MYB31). Lower part, growth on selective media deprived of W and H amino acids. Physcomitrella patens protoplasts transient expression assays (right panels). Green fluorescent protein (GFP) intensity measured in P. patens protoplasts cotransfected with sequences 3, 7, 5 and 9 fused to the 35S cauliflower mosaic virus minimal promoter and the GFP reporter gene with either, the NAC (sequences 3 and 9) or the MYB31 (sequences 7 and 5) DNA binding domain fused to the VP16 activation domain. Error bars ± SE. t-test significance: *, P < 0.05 and ***, P < 0.001.

Table 3 Summary of the interactions identified and verified in yeast one-hybrid experiments

Within motif group I, two different sequences were assayed. For the first one, AACGTGGG, no specific class of TF was found to preferentially activate transcription in yeast. In this assay, five different TFs, belonging to five different families (i.e. bHLH, GATA, MADS, 3R-MYB and NAC), were identified (Table 3). Conversely, seven bHLH TFs (i.e. bHLH31, 48, 64, 72, 74, 79 and 115) were able to activate the transcription from the second sequence, GCACGTGGAG, revealing some transcriptional specificity (Table 3). Interestingly, both elements contain the ACGT core sequence, which once mutated (AACCTCGG and GCAGCTGGAG) prevented the identified TFs to activate transcription in yeast, indicating its preponderant role in these interactions. The Y1H analysis of group V motifs, which belong to the cluster showing the highest level of sequence similarity with group I, led to the identification of four bHLH TFs. In fact, these bHLH factors were found to activate transcription in yeast from the ATGTGATGC sequence, but not from the GCATCACCC sequence, for which TFs from two other families (i.e. G2-like and Zn finger) were identified. The careful analysis of these group V motif sequences and their corresponding mutated versions ATGCAATGC and GCATAGCCC) exposed that when the ATGTGATGC sequence is concatenated, an E-box (CANNTG) is created, whereas this is not the case with the GCATCACCC sequence. The fact that sequences from group I and V display strong similarities with either the well described G ([GC] ACGT [GC]) or E boxes probably explains why bHLH TFs were found to be associated with these two motif groups. Group II is closely related, at the sequence level, to group I and V. However, unlike groups I and V, only two NAC TFs were found to interact with the two sequences tested (i.e. GCTGCCGGAGA and GCCACGTCAGC), with one (At3g12910) being identical for both sequences (Table 3).

From group III, four sequences were analysed, namely CACCTAAC, ACCACAACC (GGTTGTGGT), ACCAAACAT and TCGGACCAA. From these AC-rich sequences only R2R3-MYB TFs were identified as able to activate transcription in yeast (Table 3), with the exception of TCGGACCAA, for which no interacting TF was identified (Table 3).

The analysis of the CTCTCTCAC element from group IV concluded with the identification of three factors: one NAC and two bZIP (subgroup I) TFs ([38, 39]; Figure 4E, Table 3).

Finally, three different sequences belonging to group VII were also analysed, namely TCTCCTCCAC, CTTTCCCC and CCTCCTTCT, from which no clear pattern of binding emerged. The screening of the TCTCCTCCAC sequence led to the identification of three different classes of TFs (bHLH, bZIP and NAC). Three additional factors belonging to three further classes (C3HC4, MYB and WRKY) where identified with the CTTTCCCC element. Finally, no binding TFs for the CCTCCTTCT element were retrieved from Regia in yeast (Table 3).

Validation of transcription factor predictions with yeast-one-hybrid experimental data

In order to benchmark the accuracy of bioinformatically predicted candidate TFs, they were compared to those experimentally identified by Y1H screening. Key search parameters were tuned as explained in Methods and in Figure 3, and a trade-off between precision and sensitivity could be achieved. In summary, we were able to lower the BLASTP E-value threshold to 1E-10 without a significant sensitivity loss, but fixed the STAMP cutoff to 1E-3 to avoid a major sensitivity drop while ensuring reliable short cis-element alignments [30]. We also found that filtering out candidate TFs with poorly conserved interface residues, those directly contacting DNA nitrogen bases, was helpful in order to discard false predictions. A cutoff of 50% interface similarity was observed to be adequate, increasing specificity to 0.921, as shown in Additional file 4 [PredictionSummary]. Tested cis-elements were classified as predictable if the footprintDB repository contained at least one annotated TF with a significantly similar motif (STAMP E-value < 1E-3), which was homologous to the experimental Y1H-inferred factors. Thirteen of the 15 cis-elements tested in the Y1H experiment were accordingly called ‘predictable', but only 5 of them (2, 5, 6, 7, 12) yielded true positives among footprintDB results. The rest had significantly divergent DNA motifs annotated in the database, with E-values over the cutoff, and therefore could not be identified (see Figure 5). It is important to note that in all predictable cases the correct binding TFs were identified by homology to annotated plant regulatory proteins, including A. thaliana, Petunia hybrida, Oryza sativa or Phaseolus vulgaris. On the other hand, group II (sequence 4), III (sequence 8) and VII (sequence 15) elements could not be successfully matched as the IDEF2 and GCN4-like motifs (see Discussion) were not part of footprintDB when the predictions were made. Overall, after conducting the benchmark experiment using as input the 15 cis-sequences in Table 1, we conclude that we correctly predicted 56 out of 81 (69%) of the TFs isolated in the Y1H assays, finding on average 1 validated TF every 24 predictions. If only predictable sites are considered, these numbers improve to 56 out of 58 (97%) and 1 confirmed TF every 13 predictions, respectively. Finally, when the complete analysis was repeated with input DNA motifs (as produced by BEST) instead of individual cis-sequences, the fraction of successfully recovered TFs dropped to 39 out of 81 (48%), yielding 1 confirmed TF every 35 predictions. These results and the complete benchmark reports can be found in Additional file 4 [PredictionSummary].

Figure 5
figure 5

Comparison between bioinformatically predicted and experimentally isolated TFs. Cis-elements were classified as ‘predictable’ (left) or as unpredictable (right) when the footprintDB repository contained at least one annotated TF homologous to the experimental Y1H-inferred binding proteins. Left vertical axis: columns report the number of TFs assigned to each tested cis-element by either footprintDB (blue) or the Y1H assay (purple), showing in green those assigned by both methods. Right vertical axis: line shows STAMP E-values for all 15 cis-elements aligned to the first hit validated by Y1H assays. Note that both vertical axes are in logarithmic scale.

Functional analysis of TFs interacting with conserved cis-sequences

In this study, some factors were found to be able to activate transcription in yeast from various cis-sequences conserved in the promoters of genes that are simultaneously induced by drought. Amongst these, two were selected for more detailed analysis, NAC protein At3g12910 and MYB31 At1g74650. These selected TFs were able to induce transcription in Y1H experiments from three (group II, GCTGCCGGAGA and GCCACGTCAGC; group IV, CTCTCTCAC; and group VII, TCTCCTCCAC) and two (group III, CACCTAAC and ACCACAACC; and group VI, CCAACTAA) different motif groups, respectively (Table 3).

The ability of these two factors to interact with their DNA target in vivo was analysed in transient expression assays, using Physcomitrella patens protoplasts. In this experiment, the DNA binding domain (DBD) of the NAC protein and MYB31 was fused to the VP16 activation domain, resulting in the formation of two synthetic TFs. These chimeric proteins were then assayed against different conserved cis-sequences, namely CTCTCTCAC and GCTGCCGGAGA (TCTCCGGCAGC), and CACCTAAC and ACCACAACC, for the NAC and MYB31 DBD, respectively (Figure 4).

This analysis confirms the results obtained in Y1H experiments. Interestingly, the transcriptional activity supported by the NAC DBD was a lot stronger with the GCTGCCGGAGA sequence than with CTCTCTCAC, strengthening the idea that the GCTGCCGGAGA element is most probably a new NAC target sequence (Figure 4B, E). Conversely, similar activity was observed when MYB31 DBD was assayed against either, CACCTAAC or ACCACAACC, suggesting that MYB31 could activate the transcription from a wide variety of AC-rich cis-sequences (Figure 4C, D).

The cis-sequences and transcription factors identified in the present study may play a role in drought-responsive gene expression in A. thaliana. In the context of this assumption, several cis-sequences were used to generate synthetic promoters, which were tested in transgenic A. thaliana plants. However, no drought-responsive reporter gene expression has been observed in these transgenic lines so far (data not shown). It will be interesting to see if overexpression or knock-down mutations of the identified transcription factors will have any effect on drought-responsiveness of the plants.

Discussion

The application of bioinformatic resources has led to the identification of large numbers of novel cis-regulatory sequences conserved in stress response genes [1–3]. Often the role of these elements and their binding transcription factors has remained unknown. The primary goal of the present work was the identification of TFs interacting with cis-sequences conserved in co-regulated genes. A total of 49 sequence motifs classified into 7 motif families were identified from microarray data on drought responsive gene expression. These motifs were bioinformatically predicted to interact with TFs belonging to specific TF families. These predictions were then tested experimentally using yeast one-hybrid screenings leading to the identification of novel TF-DNA interactions. Initially, the two bioinformatic approaches employed, preferably predicted TFs interacting with putative binding sites for motif group I that harbours 29 of a total of 49 motifs. Using the STAMP web server as well as footprintDB, these show significant similarities (low E-value) to G-box and ABRE elements. However, no consistency in TF predictions was obtained with these database-assisted analyses using the other motifs. This first set of results exposed the natural limitations of this approach, mainly associated with the data content of the underlying scanned databases. For example, predictions using the STAMP webserver depend on previously annotated data derived from databases such as TRANSFAC, JASPAR, PLACE, AGRIS, and AthaMap [23, 28, 40–42]. The content of these repositories is primarily derived from functionally known regulatory sequences, which are expected to be only a tiny fraction of all cis-elements. Similarly, footprintDB (http://floresta.eead.csic.es/footprintdb) contains TFs and their experimentally-associated DNA motifs from these and other databases. A main difference to the STAMP webserver is that footprintDB motifs are associated to their cognate binding TFs, and their DNA-binding interfaces are annotated using 3D-footprint structural data [29]. However, STAMP as well as footprintDB detect and predict TF-DNA interactions based on known data. Therefore, TF-DNA interactions not readily predicted with STAMP and footprintDB may represent novel or to date unknown TF-DNA interactions. This was substantiated in the present study by using sequences enriched in promoters of drought-responsive genes, but not readily associated with interacting TFs for experimental analysis using Y1H-screenings. Towards these ends two different sequences from motif group I were assayed. For the first one, AACGTGGG, no specific class of TF was found to preferentially induce transcription in yeast. Conversely, seven bHLH TFs were able to activate the transcription from the second sequence, GCACGTGGAG. Both elements contain the ACGT core sequence, which once mutated (AACCTCGG and GCAGCTGGAG) prevented transcription to be activated. These observations suggest that sequences AACGTGGG and GCACGTGGAG are bona fide binding sites for bHLH31, 48, 64, 72, 74, 79 and 115. Motif group V showed the highest levels of sequence similarity with group I and drove the identification of four bHLH TFs. In fact, the bHLH TFs were found to prompt transcription in yeast from cis-element ATGTGATGC, but not from GCATCACCC. When sequence ATGTGATGC is concatenated, an E-box (CANNTG) is created (i.e. ATGTGATGCATGTGATGC), whereas this is not the case with GCATCACCC. The fact that sequences from group I and V contain or display strong similarities with either the well described G-box ([GC] ACGT [GC]) or E-box probably explains why bHLH TFs were found to be associated with these two motif groups. In fact, numerous studies have shown that bHLH TFs regulate the activity of their target gene promoters through these types of cis-regulatory sequences [43, 44]. However, no bZIP TFs were identified as interacting with the motif that contains the ACGT core, suggesting that the interactions with the bHLH TFs are quite specific and that the full-length sequence of the motifs may have a role in this specificity [38]. Group II is closely related, at the sequence level, to group I and V. However, unlike group I and V, only two NAC TFs were found to interact with the two sequences tested (GCTGCCGGAGA and GCCACGTCAGC), with one (At3g12910) being identical for both sequences (Table 3). A wide variety of NAC binding sites have been identified so far in various plant species, revealing the large diversity of the sequences that can be recognised by this group of TFs [45–47]. The GCCACGTCAGC sequence displays similarities with the core sequence CA [AC] G [TC] [TCA] [TCA] that is recognised by the rice (Oryza sativa) and barley (Hordeum vulgare L.) IDEF2 NAC proteins, suggesting that the identified NAC could activate the transcription in yeast from it [47]. This assertion was supported by the lack of transcriptional activity associated with the mutagenesis of this element (GCCATATCAGC). Interestingly, the GCTGCCGGAGA sequence appears to be quite divergent when compared to the other described NAC binding sites even though both NAC TFs were able to activate the transcription in yeast from it, suggesting that the GCTGCCGGAGA DNA sequence could be a new cis-regulatory element that could be targeted by NAC proteins (Figure 4B).

From group III, the AC-rich sequences CACCTAAC, ACCACAACC (GGTTGTGGT), ACCAAACAT and TCGGACCAA identified only R2R3-MYB TFs to be able to activate transcription in yeast (Table 3). This result is in agreement with previous studies that have demonstrated that numerous R2R3-MYB proteins from different plant species are able to bind to and induce transcription from AC-rich sequences [48–50]. One of the studied sequences (CACCTAAC) contains the well-described AC-I cis-regulatory element (ACCTAAC) to which the A. thaliana MYB61 TF (R2R3-MYB subgroup 13) was shown to bind to in vitro and activate transcription from it in yeast [50]. This result was confirmed in our screen (Figure 4C). Similarly, we found that two other R2R3-MYB subgroup 13 members (composed of four TF genes), namely MYB50 and MYB86, were also able to activate transcription in yeast from the same element [51]. The overlap between the R2R3-MYB proteins that interact with the CACCTAAC (16 interacting R2R3-MYB) and the ACCACAACC (14 interacting R2R3-MYB) sequences was quite high, as 12 out of the 18 interacting R2R3-MYBs (67%) were common to both elements (Table 3, Figure 4C, D). Interestingly, only four R2R3-MYBs prompted transcription in yeast from the ACCAAACAT element, from which three belong to the same R2R3-MYB subgroup, namely subgroup 20 (composed of six members) [51]. Group VI, such as Goup III, is composed of AC-rich motifs. From this motif group only one putative cis-regulatory sequence was assayed, namely CCAACTAA, from which only R2R3-MYB TFs (17 genes) were also identified (Table 3). Eleven (64,7%) and 12 (70,6%) of the identified R2R3-MYBs were identical with the R2R3-MYBs found to interact with the CACCTAAC and ACCACAACC sequences, respectively. When concatenated, this motif partially contains the R2R3-MYB targeted AC-II element (ACCAACC) from which the identified TFs are likely to induce transcription in yeast, based on the absence of activation with the mutated version (GTTAAC) of the element [48].

Surprisingly, as no ACGT core sequence (A-box, C-box or G-box) is present in this motif, the analysis of the CTCTCTCAC element from group IV leads to the identification of two bZIP (subgroup I) factors ([38, 39]; Figure 4E, Table 3). Nevertheless, a detailed analysis of this sequence showed that a bZIP-like DNA target (GGTGAGAGAG) similar to two GCN4-like (GGTGAG and TGTGTGACA) motifs found in the promoter of the wheat storage protein genes was present, suggesting that this sequence could be mainly recognised by this class of bZIP factors [52].

While these Y1H results on their own already provide valuable biological insight, they also served in this work as a reference dataset. Indeed they were also used to benchmark the predictive value of the footprintDB approach for recognizing putative binding TFs with different parameter settings. Two STAMP E-value cutoffs were tested. While a stringent 1E-5 threshold yielded the most specific predictions, it was at the cost of a significantly reduced sensitivity (52%). On the contrary, a relaxed 1E-3 cut-off increased the sensitivity (78%) by compromising specificity. However, this loss of specificity could be corrected by lowering the BLASTP E-value cut-off to 1E-10. In addition, we observed a further gain in specificity by applying an interface similarity filter for candidate TFs, requiring them to have at least 50% similar interface residues to be considered. Altogether, these settings yield a 69% sensitivity and a 92% specificity. Another important observation regarding the performance was that individual input cis-elements worked much better than sequence motifs, which presumably increase the chance of obtaining significant but irrelevant alignments of DNA sequences, which are the cornerstone of this approach.

In addition to the Y1H screenings, the in vivo interaction of two TFs with their DNA target was analysed in transient expression assays using Physcomitrella patens protoplasts. In this experiment, two synthetic TFs were generated by fusing the DNA binding domain (DBD) of the NAC protein (At3g12910) and MYB31 (At1g74650) to the VP16 activation domain. These chimeric proteins were then assayed against the conserved cis-sequences, CTCTCTCAC and GCTGCCGGAGA, and CACCTAAC and ACCACAACC, for the NAC and MYB31 DBD, respectively (Figure 4). Interestingly, the transcriptional activity supported by the NAC DBD was a lot stronger with the GCTGCCGGAGA sequence than with CTCTCTCAC one, strengthening the idea that the GCTGCCGGAGA element is most probably a new NAC target sequence (Figure 4B, E). To date the NAC protein encoded by the Arabidopsis gene At3g12910 had not been investigated experimentally. Comparing the proposed novel NAC binding site GCTGCCGGAGA with known Arabidopsis NAC binding sites reveals a large variability of binding sites recognised by NAC proteins. For example, the NAM (At1g52880) and NAP (At1g69490) proteins bind to the cis-sequences AAGGGATGA and CACGTAAGT, respectively [53, 54]. The high variability of NAC binding sites is also illustrated by the A-rich binding site of ATAF2 [55]. NAC transcription factors are involved in a wide array of abiotic stress responses [56]. Interestingly, the NAC transcription factors ANAC019 (At1g52890), ANAC055 (At3g15500), and ANAC072 (At4g27410) bind to the sequence ANNNNNTCNNNNNNNACACGCATGT, a drought responsive cis-sequence [45]. Although this sequence is part of a drought responsive promoter (ERD1), overexpression of ANAC019, ANAC055, and ANAC072 in transgenic plants did not up-regulate ERD1 expression indicating that other interacting factors may be necessary for the induction of the ERD1 gene [45]. In this context it is interesting to note that none of the cis-sequences investigated here, although enriched in drought responsive genes, individually confer drought responsive gene expression (data not shown). Among these sequences was also the novel NAC binding site GCTGCCGGAGA. Similarly it may be suggested that the investigated cis-sequences are each part of a combinatorial element that requires a second binding site for an interacting transcription factor in close vicinity for its functionality [57].

When the MYB31 DBD was assayed against CACCTAAC and ACCACAACC in Physcomitrella patens protoplasts using a synthetic TF, similar activity was observed with both cis-sequences, suggesting that MYB31 could induce transcription from a wide variety of AC-rich cis-sequences (Figure 4C, D). MYB31 belongs to subgroup 1 of R2R3MYB TFs [51]. The gene is known to be up-regulated by chitin but no further functional analysis has been carried out to date [58]. There are 126 R2R3-MYB TFs in the genome of Arabidopsis[51] and random binding site selection experiments with MYB15, 61, 77, and 84 revealed that these factors always bind to cis-sequences with one or two conserved ACC and/or AAC core sequences [59, 60]. Consequently, the adjacent nucleotides can vary significantly and Arabidopsis MYB factors have a relatively degenerate binding site recognition. Consistent with known MYB binding sites, the synthetic TF containing the MYB31 DBD activates gene expression from two synthetic promoters harbouring the ACC and AAC motif (CACCTAAC and ACCACAACC). Although the MYB31 binding site CACCTAAC does not confer drought responsive reporter gene expression in transgenic Arabidopsis (data not shown), this motif has been found in Arabidopsis in the promoter of three MYB61 target genes that are involved in vasculature development, one of the main components of the transpiration stream, which indirectly may participate in plant adaptation to drought stress [50]. Such observations reinforce the postulate that a combination of regulatory elements is required to integrate the environmental signals on a specific gene promoter [57].

Conclusions

The work presented demonstrates the successful integration of several bioinformatic resources to predict and validate TFs interacting with conserved sequence motifs in co-regulated genes. Predictions were confirmed by using a yeast-one-hybrid approach to identify interacting TFs belonging to the predicted TF families. TF-DNA interactions were further experimentally validated in yeast and with a Physcomitrella patens transient expression system, leading to the discovery of several novel TF-DNA interactions. Our work establishes a novel approach to identify TFs interacting with conserved cis-sequences. This approach may facilitate the experimental identification of TFs because a candidate TF-family can be predicted bioinformatically using the footprintDB database. Although the limitation of this approach is the content of the databases used, footprintDB is particularly valuable because it contains data for TF-DNA interactions from many different species. Therefore, this database may help to identify a DNA binding protein domain for any submitted sequence. Afterwards one may check plant TF databases for the predicted DNA binding domain to identify novel sequence-specific TF-DNA interactions.

Methods

Bioinformatic analysis to identify conserved cis-sequences in co-regulated genes

In silico identification and bioinformatic analysis of sequence motifs overrepresented in drought stress-induced A. thaliana genes was based on microarray expression data from the AtGenExpress global stress expression data set [61]. To identify genes up-regulated by drought, Affymetrix ATH1 microarray raw expression data were downloaded from NASCArrays (NASCARRAYS-141) [62], normalised using MAS5 algorithm [63], and scaled to a TGT of 100. An untreated control data set (NASCARRAYS-137) was also downloaded and identically processed. Array elements were assigned to genes according to an assignment table based on TAIR release 8 annotations [64]. Expression data were imported into the PathoPlant database [20]. The annotation procedure of cDNA microarray data and Affymetrix ATH1 data has been described earlier [3, 21]. All data as well as links to the microarray source of the expression set can be found on the PathoPlant site at http://www.pathoplant.de/documentation.php.

Expression data was used to identify genes up-regulated upon drought stress. Genes showing at least a 2-fold induction compared to untreated control were defined as up-regulated. Using an SQL server query tool, 32 drought stress combinations at different time points and in different tissues were queried. PathoPlant’s ‘Microarray expression’ online tool displays a similar functionality to the query tool described above and can be used to determine sets of genes co-regulated upon drought.

For promoter analyses of drought-induced gene sets, sequences 1,000 nt upstream of the transcription start site, if known as for the majority of genes, otherwise 1,000 nt upstream of the ATG start codon were extracted using TAIR release 8 sequence annotation and converted into FASTA format. To identify overrepresented motif sequences within these promoters, the BEST software package [19], locally installed on a Linux SuSE9.2 system was used [3]. The package combines 4 different motif-finding programs (MEME, AlignACE, CONSENSUS, BioProspector) and an optimization step. BEST was run with default parameters and predefined motif lengths of 5 to 10 nucleotides. The application of these parameters had previously shown to yield optimal results with promoter sequences from A. thaliana[3]. Overrepresented sequence motifs identified by BEST were further used if detected by at least 2 out of the 4 motif finding programs.

Such a BEST analysis with 33 gene sets yields a high number of enriched sequence motifs. In order to pre-validate these bioinformatically, a so-called in silico expression analysis was performed with all motif sequences (http://www.pathoplant.de/expression_analysis.php). This analysis is based on PathoPlant’s microarray expression database and correlates motif sequence occurrences with stress-specific gene expression data resulting in an evaluation of the sequence to identify the stress this sequence is most likely to be responsive for [31]. Strict correlation criteria were applied by only considering motif sequences that displayed at least two stresses associated with drought among the first 3 highest ranked stresses.

Identified sequence motifs in which at least one sequence was found to be enriched in promoters of drought stress responsive genes were submitted to the STAMP web server applying the recommended default parameters [22, 40]. STAMP classified all motifs based on matrix alignment to a similarity tree given in Newick tree format that was displayed using MEGA [32]. Groups containing similar motifs were generally defined by clustering single motifs on branches with lengths <0.05. In one exceptional case a motif (shoots3hr_1) was grouped into motif group II based on its conserved core sequence GCC that is characteristic for AP2/EREBP binding proteins. This core sequence is also conserved in the other two motifs of motif group II. STAMP was also employed for the identification of motif similarities by comparison with known cis-elements from plant databases AthaMap, AGRIS and PLACE [23–25].

Bioinformatic prediction of candidate transcription factors for selected cis-elements

Sequence motifs obtained for drought stress related genes were used as queries for footprintDB (http://floresta.eead.csic.es/footprintdb). This repository currently contains over 5,000 unique TFs and their DNA binding preferences annotated as position-specific scoring matrices (PSSMs) that capture the occurrence of nucleotides in different positions of the DNA binding site [65]. When building footprintDB for this analysis, data were extracted from the literature and other repositories such as TRANSFAC, JASPAR, 3D-footprint, UniPROBE, AthaMap, DBTBS, RegulonDB [12, 29, 40, 41, 66–68], and several papers reporting human, murine, and Drosophila motifs [69, 70].

The footprintDB engine searches for similar motifs using STAMP [22]. The main difference to the STAMP webserver is that footprintDB motifs are associated to their cognate binding TFs, and their DNA-binding interfaces are annotated using 3D-footprint structural data [29]. Interface residues are defined for being located within 4.5 Angstroms with respect to at least one DNA nitrogen base in homologous structural complexes deposited at the Protein Data Bank. By taking advantage of these data, footprintDB is able to predict A. thaliana homologous TFs that are candidates to bind the input cis-sequences, calculating interface similarity with a custom scoring matrix motifs [30]. The search parameters were: 1) STAMP E-value: 1E-3, 2) BLASTP E-value: 1E-10, 3) Arabidopsis protein source: TAIR9, containing 33,410 protein sequences [64], 4) Interface similarity cut-off: 50%.

Yeast one-hybrid (Y1H) experiments

All the primers used in this section were purchased from Sigma-Aldrich (Lyon, France) and are listed in Additional file 6.

The putative cis-regulatory sequences identified in this study were synthesised as hexamers, cloned (EcoRI and XbaI) into the pHISi vector (Clontech, Saint-Germain-en-Laye, France) and stably transformed into yeast (Saccharomyces cerevisiae, EGY48α-type mating strain) at the URA3 locus. Fifteen different cis-sequences, from seven motif subgroups, that were not inducing self-activation in yeast and that were preferentially not already described in the literature as associated with drought stress response, were selected for experiments. For this purpose, a yeast (YM4271, a-type mating strain) normalised cDNA library of A. thaliana TFs cloned into pDEST22™ (Invitrogen, Saint Aubin, France) was used (REGIA library [71]).

Following the mating of both yeast strains, diploid colonies growing on a medium lacking the histidine amino acid were considered as positive clones expressing the candidate TFs interacting with the studied sequences. Then, a verification step was carried out, which consisted in testing the candidate TF’s activity from both, their target sequence, and their corresponding mutated version. When yeast growth was observed from the unmodified sequence and not from the mutated version, the corresponding TF was considered as positive interactor. In order to increase the stringency of the screens, various concentrations of 3-aminotriazol (3-AT) were used in all the Y1H experiments (from 15 to 60 mM). The details of the yeast transformation, mating and selection are reported elsewhere [72].

Finally, for each of the confirmed interactions, the expressed cDNA was amplified from the diploid colonies by PCR and the resulting amplicons were subsequently sequenced in order to confirm the identity of the identified TFs. In order to facilitate the PCR reaction, the yeast cell wall was hydrolysed by a lyticase (Sigma-Aldrich, Lyon, France) treatment. Briefly, yeast cells were suspended in 15 μl lyticase solution (2 mg/ml in 0.1 M sodium phosphate buffer, pH 7.5), incubated 30 min at 37°C and then 10 min at 95°C (heat inactivation), and finally diluted by adding 85 μl of sterile water. Five μl lysate was then used per 50 μl PCR reaction.

Bioinformatic validation of transcription factor predictions with yeast-one-hybrid experimental data

Y1H experimental results allowed us to assess the performance of footprintDB transcription factor predictions and to tune the search parameters to increase its effectiveness. The set of 15 cis-sequences used for Y1H screenings was analysed in order to predict putative A. thaliana binding TFs. Different strategies and parameter combinations were tested and their impact on the results was evaluated, since we observed that the systematic comparison of short DNA sequences, such as cis-elements, was troublesome [73]. Two A. thaliana protein sequence sources were tested: the list of 995 TFs cloned in the Y1H library (from REGIA) and the TAIR9 longest transcript proteome [64]. The resulting optimum strategy was employed in order to predict putative binding TFs. Parameter setting evaluation and the complete predictions are reported in Additional file 4. This file contains several sheets: 1) the data supporting Figure 5; 2) a table with performance measurements of footprintDB with different parameter settings [PredictionSummary]; 3) the complete set of predictions for input motifs [Prediction (motifs)]; 4) the complete set of predictions for input individual DNA sites [Prediction (sites)]; 5) A file with all A. thaliana homologous TFs and their annotations [A.th.homologues]; and 6) a summary of significantly similar motifs found within footprintDB for each input DNA motif [MotifSimilarity].

Physcomitrella patenstransfection assays

All the primers used in this section were purchased from Sigma-Aldrich (Lyon, France) and are listed in Additional file 6.

Four sequences belonging to three different motif groups were assayed against DNA binding domains from two TFs (NAC, At3g12910 and MYB31, At1g74650) fused to the VP16 activation domain in P. patens protoplast transient expression assays. These sequences were synthesised as hexamers fused to the 35S cauliflower mosaic virus minimal promoter and recombined into the pBS TPp-B vector [74]. The NAC and MYB31 VP16-fusions were carried out by PCR using the high-fidelity Phusion DNA polymerase with the HF buffer (Thermo Fisher Scientific, Villebon sur Yvette, France) prior to recombination into the pBS TPp-A vector [74]. Gateway® recombination, protoplast transformation and quantitative analysis of reporter gene expression are described in detail elsewhere [74].

References

  1. Harb A, Krishnan A, Ambavaram MM, Pereira A: Molecular and physiological analysis of drought stress in Arabidopsis reveals early responses leading to acclimation in plant growth. Plant Physiol. 2010, 154: 1254-1271. 10.1104/pp.110.161752.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  2. Zou C, Sun K, Mackaluso JD, Seddon AE, Jin R, Thomashow MF, Shiu SH: Cis-regulatory code of stress-responsive transcription in Arabidopsis thaliana. Proc Natl Acad Sci U S A. 2011, 108: 14992-14997. 10.1073/pnas.1103202108.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  3. Koschmann J, Machens F, Becker M, Niemeyer J, Schulze J, Bülow L, Stahl DJ, Hehl R: Integration of bioinformatics and synthetic promoters leads to the discovery of novel elicitor-responsive cis-regulatory sequences in Arabidopsis. Plant Physiol. 2012, 160: 178-191. 10.1104/pp.112.198259.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  4. Hehl R, Wingender E: Database-assisted promoter analysis. Trends Plant Sci. 2001, 6: 251-255. 10.1016/S1360-1385(01)01954-9.

    Article  CAS  PubMed  Google Scholar 

  5. Hehl R, Steffens NO, Wingender E: Isolation and analysis of gene regulatory sequences. Handbook of Plant Biotechnology. Edited by: Christou P, Klee H. 2004, Weinheim Germany: John Wiley and Sons Ltd, 81-102.

    Google Scholar 

  6. Hehl R, Bülow L: Internet resources for gene expression analysis in Arabidopsis thaliana. Curr Genomics. 2008, 9: 375-380. 10.2174/138920208785699535.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  7. Brady SM, Provart NJ: Web-queryable large-scale data sets for hypothesis generation in plant biology. Plant Cell. 2009, 21: 1034-1051. 10.1105/tpc.109.066050.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  8. Priest HD, Filichkin SA, Mockler TC: Cis-Regulatory elements in plant cell signaling. Curr Opin Plant Biol. 2009, 12: 643-649. 10.1016/j.pbi.2009.07.016.

    Article  CAS  PubMed  Google Scholar 

  9. Usadel B, Obayashi T, Mutwil M, Giorgi FM, Bassel GW, Tanimoto M, Chow A, Steinhauser D, Persson S, Provart NJ: Co-expression tools for plant biology: opportunities for hypothesis generation and caveats. Plant Cell Environ. 2009, 32: 1633-1651. 10.1111/j.1365-3040.2009.02040.x.

    Article  CAS  PubMed  Google Scholar 

  10. Galuschka C, Schindler M, Bülow L, Hehl R: AthaMap web-tools for the analysis and identification of co-regulated genes. Nucleic Acids Res. 2007, 35: D857-D862. 10.1093/nar/gkl1006.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  11. Chang WC, Lee TY, Huang HD, Huang HY, Pan RL: PlantPAN: Plant promoter analysis navigator, for identifying combinatorial cis-regulatory elements with distance constraint in plant gene groups. BMC Genomics. 2008, 9: 561-10.1186/1471-2164-9-561.

    Article  PubMed Central  PubMed  Google Scholar 

  12. Bülow L, Engelmann S, Schindler M, Hehl R: AthaMap, integrating transcriptional and post-transcriptional data. Nucleic Acids Res. 2009, 37: D983-D986. 10.1093/nar/gkn709.

    Article  PubMed Central  PubMed  Google Scholar 

  13. Bailey TL, Elkan C: The value of prior knowledge in discovering motifs with MEME. Proc Int Conf Intell Syst Mol Biol. 1995, 3: 21-29.

    CAS  PubMed  Google Scholar 

  14. Roth FP, Hughes JD, Estep PW, Church GM: Finding DNA regulatory motifs within unaligned noncoding sequences clustered by whole-genome mRNA quantitation. Nat Biotechnol. 1998, 16: 939-945. 10.1038/nbt1098-939.

    Article  CAS  PubMed  Google Scholar 

  15. Hertz GZ, Stormo GD: Identifying DNA and protein patterns with statistically significant alignments of multiple sequences. Bioinformatics. 1999, 15: 563-577. 10.1093/bioinformatics/15.7.563.

    Article  CAS  PubMed  Google Scholar 

  16. GuhaThakurta D, Stormo GD: Identifying target sites for cooperatively binding factors. Bioinformatics. 2001, 17: 608-621. 10.1093/bioinformatics/17.7.608.

    Article  CAS  PubMed  Google Scholar 

  17. Liu X, Brutlag DL, Liu JS: BioProspector: discovering conserved DNA motifs in upstream regulatory regions of co-expressed genes. Pac Symp Biocomput. 2001, 127-138.

    Google Scholar 

  18. Eskin E, Pevzner PA: Finding composite regulatory patterns in DNA sequences. Bioinformatics. 2002, 18 (Suppl 1): S354-S363. 10.1093/bioinformatics/18.suppl_1.S354.

    Article  PubMed  Google Scholar 

  19. Che D, Jensen S, Cai L, Liu JS: BEST: binding-site estimation suite of tools. Bioinformatics. 2005, 21: 2909-2911. 10.1093/bioinformatics/bti425.

    Article  CAS  PubMed  Google Scholar 

  20. Bülow L, Schindler M, Choi C, Hehl R: PathoPlant: A database on plant-pathogen interactions. In Silico Biol. 2004, 4: 529-536.

    PubMed  Google Scholar 

  21. Bülow L, Schindler M, Hehl R: PathoPlant: a platform for microarray expression data to analyze co-regulated genes involved in plant defense responses. Nucleic Acids Res. 2007, 35: D841-D845. 10.1093/nar/gkl835.

    Article  PubMed Central  PubMed  Google Scholar 

  22. Mahony S, Benos PV: STAMP: a web tool for exploring DNA-binding motif similarities. Nucleic Acids Res. 2007, 35: W253-W258. 10.1093/nar/gkm272.

    Article  PubMed Central  PubMed  Google Scholar 

  23. Higo K, Ugawa Y, Iwamoto M, Korenaga T: Plant cis-acting regulatory DNA elements (PLACE) database: 1999. Nucleic Acids Res. 1999, 27: 297-300. 10.1093/nar/27.1.297.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  24. Davuluri RV, Sun H, Palaniswamy SK, Matthews N, Molina C, Kurtz M, Grotewold E: AGRIS: Arabidopsis Gene Regulatory Information Server, an information resource of Arabidopsis cis-regulatory elements and transcription factors. BMC Bioinformatics. 2003, 4: 25-10.1186/1471-2105-4-25.

    Article  PubMed Central  PubMed  Google Scholar 

  25. Steffens NO, Galuschka C, Schindler M, Bülow L, Hehl R: AthaMap: an online resource for in silico transcription factor binding sites in the Arabidopsis thaliana genome. Nucleic Acids Res. 2004, 32: D368-D372. 10.1093/nar/gkh017.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  26. Bülow L, Steffens NO, Galuschka C, Schindler M, Hehl R: AthaMap: from in silico data to real transcription factor binding sites. In Silico Biol. 2006, 6: 243-252. 10.1007/3-540-28185-1_10.

    Article  PubMed  Google Scholar 

  27. Palaniswamy SK, James S, Sun H, Lamb RS, Davuluri RV, Grotewold E: AGRIS and AtRegNet. a platform to link cis-regulatory elements and transcription factors into regulatory networks. Plant Physiol. 2006, 140: 818-829. 10.1104/pp.105.072280.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  28. Yilmaz A, Mejia-Guerra MK, Kurz K, Liang X, Welch L, Grotewold E: AGRIS: the Arabidopsis Gene Regulatory Information Server, an update. Nucleic Acids Res. 2011, 39: D1118-D1122. 10.1093/nar/gkq1120.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  29. Contreras-Moreira B: 3D-footprint: a database for the structural analysis of protein-DNA complexes. Nucleic Acids Res. 2010, 38: D91-D97. 10.1093/nar/gkp781.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  30. Sebastian A, Contreras-Moreira B: footprintDB: a database of transcription factors with annotated cis elements and binding interfaces. Bioinformatics. 2014, 30: 258-265. 10.1093/bioinformatics/btt663.

    Article  CAS  PubMed  Google Scholar 

  31. Bolívar JC, Machens F, Brill Y, Romanov A, Bülow L, Hehl R: 'In silico expression analysis’, a novel PathoPlant web-tool to identify abiotic and biotic stress conditions associated with specific cis-regulatory sequences. Database (Oxford). 2014, 2014: bau030-10.1093/database/bau030.

    Article  Google Scholar 

  32. Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.

    Article  CAS  PubMed  Google Scholar 

  33. Giuliano G, Pichersky E, Malik VS, Timko MP, Scolnik PA, Cashmore AR: An evolutionarily conserved protein binding sequence upstream of a plant light-regulated gene. Proc Natl Acad Sci U S A. 1988, 85: 7089-7093. 10.1073/pnas.85.19.7089.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  34. Guiltinan MJ, Marcotte WR, Quatrano RS: A plant leucine zipper protein that recognizes an abscisic acid response element. Science. 1990, 250: 267-271. 10.1126/science.2145628.

    Article  CAS  PubMed  Google Scholar 

  35. Iwasaki T, Yamaguchi-Shinozaki K, Shinozaki K: Identification of a cis-regulatory region of a gene in Arabidopsis thaliana whose induction by dehydration is mediated by abscisic acid and requires protein synthesis. Mol Gen Genet. 1995, 247: 391-398. 10.1007/BF00293139.

    Article  CAS  PubMed  Google Scholar 

  36. Ono A, Izawa T, Chua NH, Shimamoto K: The rab16B promoter of rice contains two distinct abscisic acid-responsive elements. Plant Physiol. 1996, 112: 483-491. 10.1104/pp.112.2.483.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  37. Busk PK, Pages M: Protein binding to the abscisic acid-responsive element is independent of VIVIPAROUS1 in vivo. Plant Cell. 1997, 9: 2261-2270. 10.1105/tpc.9.12.2261.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  38. Jakoby M, Weisshaar B, Dröge-Laser W, Vicente-Carbajosa J, Tiedemann J, Kroj T, Parcy F: bZIP transcription factors in Arabidopsis. Trends Plant Sci. 2002, 7: 106-111. 10.1016/S1360-1385(01)02223-3.

    Article  CAS  PubMed  Google Scholar 

  39. Correa LG, Riano-Pachon DM, Schrago CG, dos Santos RV, Mueller-Roeber B, Vincentz M: The role of bZIP transcription factors in green plant evolution: adaptive features emerging from four founder genes. PLoS One. 2008, 3: e2944-10.1371/journal.pone.0002944.

    Article  PubMed Central  PubMed  Google Scholar 

  40. Matys V, Fricke E, Geffers R, Gossling E, Haubrock M, Hehl R, Hornischer K, Karas D, Kel AE, Kel-Margoulis OV, Kloos DU, Land S, Lewicki-Potapov B, Michael H, Münch R, Reuter I, Rotert S, Saxel H, Scheer M, Thiele S, Wingender E: TRANSFAC: transcriptional regulation, from patterns to profiles. Nucleic Acids Res. 2003, 31: 374-378. 10.1093/nar/gkg108.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  41. Sandelin A, Alkema W, Engstrom P, Wasserman WW, Lenhard B: JASPAR: an open-access database for eukaryotic transcription factor binding profiles. Nucleic Acids Res. 2004, 32: D91-D94. 10.1093/nar/gkh012.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  42. Bülow L, Bolívar JC, Ruhe J, Brill Y, Hehl R: MicroRNA Targets’A new AthaMap web-tool for genome-wide identification of miRNA targets in Arabidopsis thaliana. Bio Data Min. 2012, 5: 7-

    Article  Google Scholar 

  43. Xu W, Grain D, Le Gourrierec J, Harscoet E, Berger A, Jauvion V, Scagnelli A, Berger N, Bidzinski P, Kelemen Z, Salsac F, Baudry A, Routaboul JM, Lepiniec L, Dubos C: Regulation of flavonoid biosynthesis involves an unexpected complex transcriptional regulation of TT8 expression, in Arabidopsis. New Phytol. 2013, 198: 59-70. 10.1111/nph.12142.

    Article  CAS  PubMed  Google Scholar 

  44. Zhang Y, Mayba O, Pfeiffer A, Shi H, Tepperman JM, Speed TP, Quail PH: A Quartet of PIF bHLH Factors provides a transcriptionally centered signaling hub that regulates seedling morphogenesis through differential expression-patterning of shared target genes in Arabidopsis. PLoS Genet. 2013, 9: e1003244-10.1371/journal.pgen.1003244.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  45. Tran LS, Nakashima K, Sakuma Y, Simpson SD, Fujita Y, Maruyama K, Fujita M, Seki M, Shinozaki K, Yamaguchi-Shinozaki K: Isolation and functional analysis of Arabidopsis stress-inducible NAC transcription factors that bind to a drought-responsive cis-element in the early responsive to dehydration stress 1 promoter. Plant Cell. 2004, 16: 2481-2498. 10.1105/tpc.104.022699.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  46. Olsen AN, Ernst HA, Lo Leggio L, Skriver K: DNA-binding specificity and molecular functions of NAC transcription factors. Plant Sci. 2005, 169: 785-797. 10.1016/j.plantsci.2005.05.035.

    Article  CAS  Google Scholar 

  47. Ogo Y, Kobayashi T, Nakanishi Itai R, Nakanishi H, Kakei Y, Takahashi M, Toki S, Mori S, Nishizawa NK: A novel NAC transcription factor, IDEF2, that recognizes the iron deficiency-responsive element 2 regulates the genes involved in iron homeostasis in plants. J Biol Chem. 2008, 283: 13407-13417. 10.1074/jbc.M708732200.

    Article  CAS  PubMed  Google Scholar 

  48. Patzlaff A, Newman LJ, Dubos C, Whetten RW, Smith C, McInnis S, Bevan MW, Sederoff RR, Campbell MM: Characterisation of Pt MYB1, an R2R3-MYB from pine xylem. Plant Mol Biol. 2003, 53: 597-608.

    Article  CAS  PubMed  Google Scholar 

  49. Prouse MB, Campbell MM: The interaction between MYB proteins and their target DNA binding sites. Biochim Biophys Acta. 1819, 2012: 67-77.

    Google Scholar 

  50. Romano JM, Dubos C, Prouse MB, Wilkins O, Hong H, Poole M, Kang KY, Li E, Douglas CJ, Western TL, Mansfield SD, Campbell MM: AtMYB61, an R2R3-MYB transcription factor, functions as a pleiotropic regulator via a small gene network. New Phytol. 2012, 195: 774-786. 10.1111/j.1469-8137.2012.04201.x.

    Article  CAS  PubMed  Google Scholar 

  51. Dubos C, Stracke R, Grotewold E, Weisshaar B, Martin C, Lepiniec L: MYB transcription factors in Arabidopsis. Trends Plant Sci. 2010, 15: 573-581. 10.1016/j.tplants.2010.06.005.

    Article  CAS  PubMed  Google Scholar 

  52. Albani D, Hammond-Kosack MC, Smith C, Conlan S, Colot V, Holdsworth M, Bevan MW: The wheat transcriptional activator SPA: a seed-specific bZIP protein that recognizes the GCN4-like motif in the bifactorial endosperm box of prolamin genes. Plant Cell. 1997, 9: 171-184. 10.1105/tpc.9.2.171.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  53. Duval M, Hsieh TF, Kim SY, Thomas TL: Molecular characterization of AtNAM: a member of the Arabidopsis NAC domain superfamily. Plant Mol Biol. 2002, 50: 237-248. 10.1023/A:1016028530943.

    Article  CAS  PubMed  Google Scholar 

  54. Zhang K, Gan SS: An abscisic acid-AtNAP transcription factor-SAG113 protein phosphatase 2C regulatory chain for controlling dehydration in senescing Arabidopsis leaves. Plant Physiol. 2012, 158: 961-969. 10.1104/pp.111.190876.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  55. Wang X, Culver JN: DNA binding specificity of ATAF2, a NAC domain transcription factor targeted for degradation by Tobacco mosaic virus. BMC Plant Biol. 2012, 12: 157-10.1186/1471-2229-12-157.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  56. Nakashima K, Takasaki H, Mizoi J, Shinozaki K, Yamaguchi-Shinozaki K: NAC transcription factors in plant abiotic stress responses. Biochim Biophys Acta. 2012, 1819: 97-103. 10.1016/j.bbagrm.2011.10.005.

    Article  CAS  PubMed  Google Scholar 

  57. Singh KB: Transcriptional regulation in plants: the importance of combinatorial control. Plant Physiol. 1998, 118: 1111-1120. 10.1104/pp.118.4.1111.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  58. Wan J, Zhang XC, Neece D, Ramonell KM, Clough S, Kim SY, Stacey MG, Stacey G: A LysM receptor-like kinase plays a critical role in chitin signaling and fungal resistance in Arabidopsis. Plant Cell. 2008, 20: 471-481. 10.1105/tpc.107.056754.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  59. Romero I, Fuertes A, Benito MJ, Malpica JM, Leyva A, Paz-Ares J: More than 80R2R3-MYB regulatory genes in the genome of Arabidopsis thaliana. Plant J. 1998, 14: 273-284. 10.1046/j.1365-313X.1998.00113.x.

    Article  CAS  PubMed  Google Scholar 

  60. Prouse MB, Campbell MM: Interactions between the R2R3-MYB Transcription Factor, AtMYB61, and Target DNA Binding Sites. PLoS One. 2013, 8: e65132-10.1371/journal.pone.0065132.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  61. Kilian J, Whitehead D, Horak J, Wanke D, Weinl S, Batistic O, D’Angelo C, Bornberg-Bauer E, Kudla J, Harter K: The AtGenExpress global stress expression data set: protocols, evaluation and model data analysis of UV-B light, drought and cold stress responses. Plant J. 2007, 50: 347-363. 10.1111/j.1365-313X.2007.03052.x.

    Article  CAS  PubMed  Google Scholar 

  62. Craigon DJ, James N, Okyere J, Higgins J, Jotham J, May S: NASCArrays: a repository for microarray data generated by NASC's transcriptomics service. Nucleic Acids Res. 2004, 32: D575-D577. 10.1093/nar/gkh133.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  63. Hubbell E, Liu WM, Mei R: Robust estimators for expression analysis. Bioinformatics. 2002, 18: 1585-1592. 10.1093/bioinformatics/18.12.1585.

    Article  CAS  PubMed  Google Scholar 

  64. Lamesch P, Berardini TZ, Li D, Swarbreck D, Wilks C, Sasidharan R, Muller R, Dreher K, Alexander DL, Garcia-Hernandez M, Karthikeyan AS, Lee CH, Nelson WD, Ploetz L, Singh S, Wensel A, Huala E: The Arabidopsis Information Resource (TAIR): improved gene annotation and new tools. Nucleic Acids Res. 2012, 40: D1202-D1210. 10.1093/nar/gkr1090.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  65. Stormo GD: DNA binding sites: representation and discovery. Bioinformatics. 2000, 16: 16-23. 10.1093/bioinformatics/16.1.16.

    Article  CAS  PubMed  Google Scholar 

  66. Robasky K, Bulyk ML: UniPROBE, update 2011: expanded content and search tools in the online database of protein-binding microarray data on protein-DNA interactions. Nucleic Acids Res. 2011, 39: D124-D128. 10.1093/nar/gkq992.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  67. Sierro N, Makita Y, de Hoon M, Nakai K: DBTBS: a database of transcriptional regulation in Bacillus subtilis containing upstream intergenic conservation information. Nucleic Acids Res. 2008, 36: D93-D96. 10.1093/nar/gkn421.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  68. Salgado H, Peralta-Gil M, Gama-Castro S, Santos-Zavaleta A, Muniz-Rascado L, Garcia-Sotelo JS, Weiss V, Solano-Lira H, Martinez-Flores I, Medina-Rivera A, Salgado-Osorio G, Alquicira-Hernandez S, Alquicira-Hernandez K, Lopez-Fuentes A, Porron-Sotelo L, Huerta AM, Bonavides-Martinez C, Balderas-Martinez YI, Pannier L, Olvera M, Labastida A, Jimenez-Jacinto V, Vega-Alvarado L, Del Moral-Chavez V, Hernandez-Alvarez A, Morett E, Collado-Vides J: RegulonDB v8.0: omics data sets, evolutionary conservation, regulatory phrases, cross-validated gold standards and more. Nucleic Acids Res. 2013, 41: D203-D213. 10.1093/nar/gks1201.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  69. Down TA, Bergman CM, Su J, Hubbard TJ: Large-scale discovery of promoter motifs in Drosophila melanogaster. PLoS Comput Biol. 2007, 3: e7-10.1371/journal.pcbi.0030007.

    Article  PubMed Central  PubMed  Google Scholar 

  70. Jolma A, Yan J, Whitington T, Toivonen J, Nitta KR, Rastas P, Morgunova E, Enge M, Taipale M, Wei G, Palin K, Vaquerizas JM, Vincentelli R, Luscombe NM, Hughes TR, Lemaire P, Ukkonen E, Kivioja T, Taipale J: DNA-binding specificities of human transcription factors. Cell. 2013, 152: 327-339. 10.1016/j.cell.2012.12.009.

    Article  CAS  PubMed  Google Scholar 

  71. Paz-Ares J, Regia Consortium: REGIA, an EU project on functional genomics of transcription factors from Arabidopsis Thaliana. Comp Funct Genomics. 2002, 3: 102-108. 10.1002/cfg.146.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  72. Castrillo G, Turck F, Leveugle M, Lecharny A, Carbonero P, Coupland G, Paz-Ares J, Onate-Sanchez L: Speeding cis-trans regulation discovery by phylogenomic analyses coupled with screenings of an arrayed library of Arabidopsis transcription factors. PLoS One. 2011, 6: e21524-10.1371/journal.pone.0021524.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  73. Sebastian A, Contreras-Moreira B: The twilight zone of cis element alignments. Nucleic Acids Res. 2013, 41: 1438-1449. 10.1093/nar/gks1301.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  74. Thevenin J, Dubos C, Xu W, Le Gourrierec J, Kelemen Z, Charlot F, Nogue F, Lepiniec L, Dubreucq B: A new system for fast and quantitative analysis of heterologous gene expression in plants. New Phytol. 2012, 193: 504-512. 10.1111/j.1469-8137.2011.03936.x.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Part of the results have been achieved within the framework of the Transnational (Germany, France, Spain) Cooperation STREG within the PLANT-KBBE Initiative, with funding from Ministerio de Ciencia e Innovación, Agence Nationale de la Recherche (ANR) and the German Federal Ministry for Education (BMBF). The work of Z.K., G.H., A.S., F.S., C.B., and L.B. was supported by the STREG project. W.X.’s work was supported by the China Scholarship Council (CSC). The authors would like to acknowledge the support of Jasmin Huebner and Nina Schmidt in the generation and analysis of transgenic A. thaliana. We also acknowledge F. Turck (MPI, Köln) and V. Colot (ENS, Paris) for providing the REGIA library.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Reinhard Hehl.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

CD, ZK, GH, WX, DG, FS, and CB performed the experiments. AS, BCM, and LB performed the bioinformatic analysis. CD, ZK, GH, AS, BCM, LB, and RH wrote the paper. CD, LL, BW, BCM, LB, and RH designed the work. All authors read and approved the final manuscript.

Christian Dubos, Zsolt Kelemen contributed equally to this work.

Electronic supplementary material

12864_2013_5999_MOESM1_ESM.xls

Additional file 1: Parameters for PathoPlant database queries and number of induced genes obtained in each of the 32 queries.(XLS 22 KB)

Additional file 2: Sequences and their alignments generating the 49 sequence motifs.(XLS 140 KB)

12864_2013_5999_MOESM3_ESM.xls

Additional file 3: Similarities of motifs to known cis -regulatory sequences in the AthaMap, AGRIS, and PLACE databases.(XLS 36 KB)

Additional file 4: Results of the footprintDB analyses.(XLS 11 MB)

Additional file 5: Results of the Y1H screen.(XLS 66 KB)

Additional file 6: Primers used for the Y1H screen and the Physcomitrella patens transfection assays.(XLS 39 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( https://creativecommons.org/publicdomain/zero/1.0/ ) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Dubos, C., Kelemen, Z., Sebastian, A. et al. Integrating bioinformatic resources to predict transcription factors interacting with cis-sequences conserved in co-regulated genes. BMC Genomics 15, 317 (2014). https://doi.org/10.1186/1471-2164-15-317

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-15-317

Keywords