- Open Access
Computational analysis of auxin responsive elements in the Arabidopsis thaliana L. genome
BMC Genomics volume 15, Article number: S4 (2014)
Auxin responsive elements (AuxRE) were found in upstream regions of target genes for ARFs (Auxin response factors). While Chip-seq data for most of ARFs are still unavailable, prediction of potential AuxRE is restricted by consensus models that detect too many false positive sites. Using sequence analysis of experimentally proven AuxREs, we revealed both an extended nucleotide context pattern for AuxRE itself and three distinct types of its coupling motifs (Y-patch, AuxRE-like, and ABRE-like), which together with AuxRE may form the composite elements. Computational analysis of the genome-wide distribution of the predicted AuxREs and their impact on auxin responsive gene expression allowed us to conclude that: (1) AuxREs are enriched around the transcription start site with the maximum density in 5'UTR; (2) AuxREs mediate auxin responsive up-regulation, not down-regulation. (3) Directly oriented single AuxREs and reverse multiple AuxREs are mostly associated with auxin responsiveness. In the composite AuxRE elements associated with auxin response, ABRE-like and Y-patch are 5'-flanking or overlapping AuxRE, whereas AuxRE-like motif is 3'-flanking. The specificity in location and orientation of the coupling elements suggests them as potential binding sites for ARFs partners.
The hormone auxin is a major regulator of plant growth and development. The influence of auxin on gene transcription is primarily mediated by its binding to TIR1/AFB receptors [1, 2]. Auxin changes the conformation of receptors and thereby promotes their interaction with Aux/IAA proteins. These proteins form heterodimers with Auxin Response Transcription Factors (ARFs) making them inactive. TIR1/AFB receptors trigger ubiquitination of Aux/IAA proteins followed by their proteasome degradation. ARFs free from their co-repressors activate or repress transcription of their target genes. ARFs bind in target promoters to the specific sites called AuxREs (Auxin Response Elements) with the TGTCNN (most frequently TGTCTC) consensus core sequence [3–5]. Recently, it has been shown, that the TGTCGG motif is more effective in binding ARF1 and ARF5 .
About half of the Arabidopsis thaliana genes have at least one TGTCTC in any orientation within the first 1000 nt of their promoter regions . Thus, TGTCTC consensus is really not a reliable method for AuxRE identification. Distribution and auxin-responsiveness of TGTCGG motifs in the plant genome have not been studied to date. While single TGTCTC hexamer does not confer auxin inducibility , this is provided by multimerized , or composite AuxREs . Two copies of the TGTCTC element oriented as a palindrome or as a direct repeat (even if both sequences are inverse) are sufficient to provide auxin response [8, 9]. ARFs bind as dimers on palindromic AuxREs where TGTCTC or TGTCGG are spaced by 5 to 9 nucleotides [3, 6].
In the composite AuxREs, TGTCNC adjoins or overlaps with coupling (constitutive) elements [4, 5]. For example, the coupling element CACGCAAT alone confers constitutive expression to a minimal promoter and the reporter gene shows no auxin response . In functional composite AuxREs TGTCNC were found both in direct and reverse orientations [5, 9, 10].
Bioinformatic analysis of cis-elements in Arabidopsis and rice disclosed bZIP- and MYB-related binding sites as potential AuxRE coupling elements . In this study, promoters (-1000 nt from TSS) of the auxin inducible genes in comparison to a randomized full genomic promoter dataset were enriched in AuxRE-, bZIP- and MYB-related elements and some of their composite modules in both genomes. The G-box (CACGTG), a type of bZIP-related motifs, binds BZR1 and PIF4 transcription factors [12, 13], which recently were shown to heterodimerize with ARF6 . ARF6 Chip-Seq analysis  revealed that composite AuxRE/G-box and AuxRE/HUD elements are highly enriched in ARF6 binding regions. The coupling motifs tend to be located close to the core AuxRE, mostly within 20 base pairs. There is also evidence of the functionality of AuxRE/MYB elements. First, the existence of ARF/MYB heterodimers was shown in vitro and in vivo for MYB77 and ARF7 transcription factors . Second, auxin induction of seven well known auxin responsive genes having multiple putative MYB binding motifs in their promoters was greatly attenuated or abolished in both myb77 and arf7 knockout mutants. Third, the proximity of the MYB and AuxRE binding sites in one of them (IAA19) suggested that ARF7/MYB77 heterodimers may bind the composite elements in the promoter of certain auxin-responsive genes.
Also, there is evidence of other transcription factors forming heterodimers with ARFs. The BIGPETALp (BPEp) basic helix-loop-helix (bHLH) transcription factor and ARF8 interact through their C-terminal domains . Double mutant analysis revealed synergic action of both transcription factors in limiting the first cell divisions during petal development and cell growth at later stages. The mutant defects were associated with changes in expression of the early auxin-responsive genes. Both transcription factors ARF5 and BES1, BZR1 homolog crucial for brassinosteroid response, bind to the SAUR15 promoter region containing composite AuxRE/HUD element .
Here we searched for the footprints of ARF dimerization with other transcription factors in the nucleotide context surrounding AuxRE. As a first approach, using information on experimentally proven AuxREs, we recognized potential AuxREs in the Arabidopsis thaliana genome. Together with the other AuxRE-related motifs [5, 18] we analyzed their association with auxin-responsive expression. This allowed us to reveal some specificity in orientation and location of potential AuxREs in the auxin responsive promoters. Second, we performed a context analysis of the flanks in experimentally proven AuxREs and found three distinct types of potential coupling motifs (Y-patch, AuxRE-like, and ABRE-like). The analysis of a number of microarray datasets assured us that the composite elements with a specific orientation of AuxRE and the coupling motifs and the certain range of spacer length between them were associated with auxin responsiveness.
Recognition of potential Auxin Responsive Elements
To apply methods for AuxRE recognition, we collected from published papers (a) the nucleotide sequences of experimentally proven AuxREs (Training set; 25 sequences 106-nt in length with centrally located TGTCNN hexamer, Additional file 1) and (b) extended promoters [-2000;-1] of auxin-regulated genes (Positive set; 44 sequences). We applied site recognition tools oPWM and SiteGA  to the Training set data. While the oPWM approach found some additional context features around TGTCNN core (Figure 1A-B, Additional file 2, Additional file 3), SiteGA revealed dependencies between locally positioned dinucleotides along the 24 nt AuxRE region (Figure 1C, Additional file 4). Next, we established the thresholds for genome-wide AuxRE recognition as those allowed to predict potential AuxREs in about 30% of promoters from the Positive set by both oPWM and SiteGA. Finally, potential AuxREs were predicted in the whole genome of Arabidopsis thaliana by: (a) oPWM/SiteGA combination (AuxREP&S); (b) TGTSTSBC consensus ; (c) TGTCTC consensus , and (d) TGTCGG consensus .
To find potential coupling elements for AuxREs, we performed the motif search in the Training set using the MotiGA software (Additional file 5). MotiGA identified a number of motifs flanking TGTCNN core. Among the most overrepresented ones three main classes of similar motifs were distinguished (Additional file 3). We investigated further (see Section 4 of the results) one representative of each class: (a) pyrimidine-rich Y-patch, (b) ABRE-like motif, and (c) AuxRE-like motif (Figure 1D-F).
Genome-wide statistics for potential AuxREs
We explored the abundances of four variants of potential AuxREs (TGTCTC, TGTCGG, TGTSTSBC and AuxREP&S) in different genomic regions of Arabidopsis thaliana: intergenic and transcript (5'UTR, exons, introns and 3'UTR) (Figure 2A). Unexpectedly, all the AuxRE variants were found to be the most abundant in 5'UTR and exons. The density of AuxREP&S were additionally enriched in introns, while TGTCGG was depleted there. TGTSTSBC, TGTCGG and AuxREP&S have the above-average distribution density in proximal [-300; +1] promoters, AuxREP&S was the most enriched among them.
The distribution densities of all the four AuxRE variants in the extended [-2000; +1] promoters were similar to the average in the genome (Figure 2A). The numbers of genes having these potential sites in the extended upstream regions were quite high: over 65% of genes contained TGTCTC, over 35% TGTSTSBC, 27.5% TGTCGG and about 10% AuxREP&S (Figure 2B; Table 1). However, the distribution of all the four AuxRE variants along upstream regions is nonuniform, having a maximum around the Transcription Start Site (TSS) (Figure 2C). It has been shown that TGTSTSBC are enriched in the proximal [-250; +1] promoter of the auxin responsive genes . Our results demonstrated that potential AuxREs are by more than 1.5 times dense around TSS than on average (except TGTCGG, which had lower density, but in a wider region). The regions of high distribution density spread both to the proximal promoter and 5'UTR. These results have an implication: 5'UTR should be considered in analysis of the auxin responsive regions.
AuxREP&Spredicts auxin-regulated genes better than other models
To test which AuxRE model predicts better auxin responsive elements, we performed meta-analysis of publicly available microarray experiments with auxin treatments (16 microarrays, Table 2). First, we created a list of auxin-regulated genes which significantly changed their expression (by more than 1.5-fold, p < 0.05) in over three microarrays experiments. The threshold for the number of microarrays was set by the binomial trial estimate (see Methods). The resulting list contained 1965 up-regulated and 1916 down-regulated genes. Second, the fractions of the significantly up- or down-regulated genes with an AuxRE variant in their promoter and 5'UTR regions were compared with that for all the genes tested in the experiment (Table 1, Additional file 6). The statistical significance of the difference between the fractions was estimated by the t-test for proportions (See Methods). The analysis showed that the targets with AuxREs predicted by simple models (TGTCTC and TGTSTSBC and AuxREPWM) were not enriched in the auxin-regulated genes, whereas those predicted by more complex models (AuxRESGA and AuxREP&S) were significantly enriched (p < 0.05). Notably, the fraction of the auxin up-regulated genes among those with predicted AuxREP&S (10.7%) was higher than the respective fractions for AuxRESGA (10.3%) and AuxREPWM (9.8%), which substantiates that usage of two models based on different principles improves prediction (Table 1). No statistical significance was found for the down-regulated by auxin genes.
The applicability of different models for AuxRE recognition was also tested using the data on phylogenetic footprinting (See Methods). For the analysis we used the data from Vista tool  on the whole genome alignments between four plant species. We tested how many potential AuxREs were conserved, namely located in an alignment with a high score and had no substitutions in the TGTC core (See Methods). About one third of AuxRE variants were found to be conserved for each of the AuxRE variant. Then, we analyzed how many of the auxin-regulated genes with a potential AuxRE obtained a conserved one (Table 1, Additional file 6). The portions were again the highest for AuxREP&S (82.7% of the up-regulated genes and 68.7% of the down-regulated). The difference in comparison with an average conservation in the regulatory regions (See Methods) was significant only for the up-regulated genes with conserved TGTSTSBC (76.7%, p < 0.05), AuxREPWM (78%, p < 0.05) and AuxREP&S (82.7%, p < 0.005).
The analysis provided proof that, indeed, AuxREP&S predicts auxin responsive elements better than the other existing models.
The preferences in location and orientation of potential AuxREs associated with auxin response
A more detailed meta-analysis of microarray data was performed to analyze the associations between the presence of AuxREP&S in a gene promoter (Additional file 7) and its auxin responsiveness. Independently of each of the 16 microarrays (Table 2), we tested the statistical significance between the fractions of the up- or down-regulated genes having and not having AuxREP&S (See Methods). The association between the presence of a potential AuxRE and auxin-responsive gene expression was considered as significant, if it was revealed in more than three microarrays, the threshold was set by the binomial trial estimate (See Methods). A significantly higher portion of auxin upregulated genes was found in five and seven microarrays for the genes with AuxREP&S predicted in the extended and proximal upstream regions, respectively (Table 3; Additional file 8). No significant association was found between the presence of a single AuxREP&S and auxin down-regulation. Interestingly, the significant associations were found mainly in the microarrays, where an early auxin response (2 h treatment and less) was investigated (Additional file 8). No association with auxin dosage was found.
Notably stronger association with auxin upregulation was found for directly oriented single AuxREP&S than for reverse ones (Table 3). However, among the genes with multiple AuxREP&S, significant association was found for reverse AuxREP&S, rather than direct AuxREP&S.
The same analysis for TGTCTC, TGTCGG and TGTSTSBC confirmed the results described above for AuxREP&S (Additional file 8) and allowed us to conclude that functional AuxREs: (1) tend to be located in proximal promoters; (2) mediate transcriptional activation in response to auxin; (3) single AuxREs preferably have direct orientation, while multiple AuxREs are reversely oriented.
Functional annotation of composite auxin responsive elements
We found three types of the coupling motifs (Figure 1D-G) on the flanks of experimentally proven AuxREs. Pyrimidine-rich TCTCCTYT motif (Figure 1D) in its sequence and distribution resembles Y-patch  or TCT-motifs , both related to transcription initiation. We designate further this coupling motif as Y-patch. HDSDYSKS motif (Figure 1F) was found as significantly matching with the EmBP-1 binding site (Additional file 9) by the TomTom tool . It was shown that the transcription factor Embp-1 interacts with the ABA-responsive element (ABRE) in the 5'-regulatory region of a wheat gene , which is why we further refer to HDSDYSKS as ABRE-like. The third type of AuxRE-flanking motif DTTBTCTH meets AuxRE itself (Figure 1G). Since B=T/G/C and H=T,G,C, a third part of the AuxRE/AuxRE-like composite elements recognized genome-wide have two TGTCNN core motifs. As AuxRE repeats were identified in promoters of the first analyzed auxin-responsive genes [5, 25], this was the expected result.
We recognized three types of composite AuxREs in A. thaliana extended upstream regions and analyzed (a) their distribution relative to TSS; (b) the functional annotation of the genes having different composite elements in their upstream regions in the AgriGO database ; (c) their impact on auxin responsive gene expression.
Composite AuxRE/Y-patch elements are asymmetrically distributed around TSS, having a maximum close to TSS with a shift to 5'UTR (Figure 3A). This distribution is a characteristic feature of the composite elements where both AuxREP&S and Y-patch were identified in the direct orientation (AuxRE+/Y-patch+). Genes having AuxRE+/Y-patch+ where functionally annotated were significantly (p < 0.002, Hochberg FDR test) related to auxin response, regulation of metabolic and cellular processes, regulation of transcription (Figure 3B). Moreover, in half of tested microarrays AuxRE+/Y-patch+ composite elements were significantly associated with auxin up-regulation (Figure 4A).
Composite element AuxRE/ABRE-like as AuxRE/Y-patch is also enriched around TSS, and the elements with the AuxRE in the direct orientation (AuxRE+/ABRE-like+ and AuxRE+/ABRE-like-) impact on this asymmetry (Additional file 10). Functional annotation and microarray data analysis revealed enrichments for AuxRE+/ABRE-like elements (Additional file 11; Additional file 9). In microarray data meta-analysis, AuxRE+/ABRE-like-elements were the most attributed to auxin response among AuxRE+/ABRE-like elements (Figure 4A; Additional file 9). Functional annotation of the genes having AuxRE+/ABRE-like composite elements also revealed their significant enrichment (p < 0.002, Hochberg FDR test) in the genes involved in regulation of various biological processes, such as development, regulation of transcription and response to external stimulus (Additional file 11).
AuxRE/AuxRE-like with direct orientation of both coupling motifs are enriched around TSS as the other two composite elements are (Additional file 10). In the analysis of microarray data, AuxRE+/AuxRE-like+ composite elements showed the best results, revealing significant associations with auxin upregulated genes in ten out of sixteen microarray experiments (Figure 4A). The functional annotation of the genes with predicted AuxRE+/AuxRE-like+ elements in their upstream regions also revealed the significant relation (p < 0.002, Hochberg FDR test) with growth and developmental processes (Additional file 12). Based on this analysis, AuxRE+/AuxRE-like+ were also found associated with shoot morphogenesis.
An important role of a spacer in the composite elements
The above analysis based on microarray data did not take into account the length of the spacer between potential AuxRE and the coupling motif. To do this, we considered six locations for the coupling motifs relative to the center of potential AuxRE either in 5' or 3' directions: overlapped (5 nt between the centers or less) and with short (from 6 to 25 nt) or long (from 26 to 50 nt) spacers (Figure 4). The subsets of genes carrying the composite element with the specific orientation and location of the coupling motif relative to AuxREP&S were analyzed separately. The composite element variant was considered as influential, if it was found significantly associated with auxin response in at least four microarrays (see Methods). In Figures 4B-C only the influential composite elements are shown. Three types of composite AuxREP&S had different preferences in the coupling motif positioning, but they all had directly oriented core AuxREP&S (Figure 4C). AuxRE-like motif tends to be located also directly on the 3'flank with no preferences to the length of the spacer. Y-patch needs to be located closely. We found two different types of AuxRE/Y-patch elements: Y-patch overlapping AuxRE core sequence at the 3' side and Y-patch located 5' to the core with a short spacer. ABRE-like motifs are located on the 5' flank of AuxRE with a long spacer.
The auxin signaling pathway has multiple players with diverse interactions. In Arabidopsis, the ARF family of transcription factors consists of 23 members, forming heterodimers with 29 Aux/IAA co-repressors . ARFs also interact with other transcriptional co-regulators like TOPLESS , SEUSS  and with other transcription factors like MYB77  and BPEp . The diversity in the interaction map for ARF transcription factors may be in some way footprinted in the regulatory regions of auxin-responsive genes, for example, as composite AuxRE elements. Indeed, ARF homodimers need two closely located binding sites with a spacer of a certain length .
Here we performed sequence analysis of the experimentally proven AuxREs with their flanks (Additional file 1) to reveal: (1) a nucleotide context of AuxRE itself; (2) coupling motifs, which together with AuxRE may form composite elements.
In search of new AuxRE variants
There were three AuxRE variants revealed to date - TGTCTC , TGTSTSBC  and TGTCGG . Distribution of TGTCGG in regulatory regions and its relation to auxin response has not been studied to date. The TGTCTC motif is present in upstream regions of approximately 65% of the A. thaliana genes (Figure 2B), which makes it over predictive for AuxRE. TGTSTSBC is less sensitive in recognizing AuxRE, but more specific than the canonical TGTCTC . Despite the gain in specificity, the number of potential AuxRE identified by TGTSTSBC consensus was still too high (Figure 2B). However, not all early auxin responsive genes have any of these two AuxRE variants in their proximal promoters. If we still suppose that auxin response of these genes is mediated by ARFs, we should expect either the existence of other AuxRE variants or another location for these AuxREs.
To search for new AuxRE variants, we collected sequences of all AuxREs experimentally confirmed to date (Additional file 1). Application of two recognition models, oPWM and SiteGA, that are based on different approaches , allowed us to reveal context features of AuxRE (Figure 1; Additional file 2, Additional file 3, Additional file 4) and to deduce the new variant for AuxRE (AuxREP&S). The genes having AuxREP&S in their extended promoters (Additional file 7) only partly overlap with those having TGTCTC and/or TGTSTSBC in the same region (data not shown). Microarray data analysis showed that the presence of AuxREP&S in the extended promoters correlates with auxin up-regulation better than those for TGTSTSBC, TGTCGG or TGTCTC (Additional file 8). All the four investigated AuxRE types showed good correlations with auxin up-regulation when they located in the proximal promoters.
We also tested the relation between conservation of an AuxRE variant in a gene promoter and its auxin responsiveness. As a result, we established that among the auxin up-regulated genes with potential AuxREP&S the highest gene fraction was obtained for the conserved ones (Table 1, Additional file 6). These genes maintain the core TGTC sequence in a conserved regulatory region among four plant species: A. thaliana, A. lyrata, T. halophile, C. papaya, and C. clementine. Among the auxin up-regulated genes of the other AuxRE variants there were lower fractions of conserved sites.
AuxREs mediate auxin responsive up-regulation
Extensive analysis of association between auxin response and the presences of each AuxRE type in the regulatory regions of various lengths allowed us to conclude that AuxREs mainly regulate auxin responsive activation of transcription. The influence of AuxREs on gene down-regulation is expected to be minor, as among thousands of tested AuxRE types (different variants of AuxREs, single/multiple AuxREs, various composite elements) none was found significantly associated with down-regulation (Additional file 9). This finding well agrees with the previously published analysis of auxin responsive transcriptome , where it was shown that the primary auxin response was almost entirely restricted to up-regulated genes. Also, the prevalence of AuxRE-dependent transcriptional activation in primary auxin response came from the analysis of Aux/IAA-ARF interactome, which revealed that sensitive to auxin Aux/IAA proteins mainly form dimers with ARF activators .
Functional AuxREs appeared to be enriched around TSS
We showed here that potential AuxREs are not uniformly distributed in the Arabidopsis genome with the highest density in 5'UTRs (Figure 2A). Aligning the AuxREs densities along the promoter and 5'UTR regions, we found that all the three AuxRE composite elements are enriched around TSS (Figure 2C). It has been demonstrated earlier  that AuxRE are enriched in proximal promoters of the auxin regulated genes. We confirmed this statement and extended it to the gene regions located downstream TSS. Two following independent results of our computational analysis support the idea that functional AuxREs appeared to be enriched around TSS.
First, microarray data analysis provides an insight that there is a strong association between the presence of potential AuxREs in the proximal promoter and 5'UTR and auxin responsiveness. Analyzing auxin responsive expression of the genes having predicted AuxREP&S, TGTCGG, TGTCTC or TGTSTSBC in [-300; 5'UTR] regions, we found among them a significantly higher portion of auxin up-regulated genes (comparing with the whole microarray gene set) in more than seven of the tested sixteen microarrays (Additional file 8). For the extended [-1500 bp; 5'UTR] regions, the significant differences were revealed in fewer microarrays: five for AuxREP&S (Table 1, Additional file 6), TGTCGG and TGTCTC and three for TGTSTSBC.
Second, all the three types of composite AuxRE elements were enriched around TSS (Figure 3). The composite elements accumulated around TSS, namely AuxRE+/Y-patch+, AuxRE+/ABRE-like and AuxRE+/AuxRE-like+, showed significant association with auxin response in microarray data analysis (Figure 4A). Functional annotation of the genes having the AuxRE+/Y-patch+ elements also showed significant enrichment of the genes related to response to auxin stimulus (Figure 3B). Functional annotation of genes carried out for other enriched around TSS composite elements (AuxRE+/AuxRE-like+, AuxRE+/ABRE-like) also showed significant results for developmental and regulatory processes (Additional file 11, Additional file 12).
Enrichment of AuxREs in 5'UTR has not been revealed earlier. However, in some auxin responsive genes potential AuxRE were described in 5'UTR. Composite AuxRE elements (AuxRE in the close vicinity to MYB binding sites) were found in 5'UTR of the IAA19 gene . ChIP analysis revealed potential AuxRE in 5'UTR of the TMO5 gene . There are two explanations of this bias for potential AuxREs towards downstream regulatory regions. The first is poor annotation of TSS in the A. thaliana genome or the existence of multiple and alternative TSS sites. The second relies on the mechanism of primary auxin response, which probably needs binding of ARF transcription factors in the close proximity to RNA polymerase complex, without specificity in location relative to TSS. In animals, there are some examples of transcription factors, whose binding sites distribution spreads below TSS, for example p53  and GAGA . Hence, it may be assumed that the distribution of some transcription factor binding sites around TSS may be a fundamental feature of the eukaryotic transcription initiation machinery.
Orientation of AuxRE influences auxin response
As there was no difference in auxin responsiveness of DR5 reporter versions with direct or reverse repeats of CCTTTTGTCTC motif , it has been proposed that there is no specificity in orientation of all AuxRE elements . And, indeed, reversely oriented AuxREs were experimentally proven in several auxin responsive genes, for example AtLBD29  and OsCRL1 .
However, by genome-wide analysis of associations between the presence of potential AuxREP&S and auxin responsive gene expression, we found significant difference for direct and reverse AuxREP&S. Single AuxREP&S significantly associated with up-regulation were directly oriented, while multiple AuxREP&S were reversed (Table 1, Additional file 6). All the composite elements, which were found associated with auxin response, had direct AuxREP&S core (Figure 4, Additional file 8).
Composite Auxin Responsive Elements
A composite element contains two functionally linked binding sites for distinct transcription factors and in this way provides the basis for integrating the inputs of two signaling pathways and enabling specificity of expression . In the simplest case, the composite element corresponds to a pair of individual binding motifs located at a particular distance from each other and involved in formation of specific DNA-protein-protein-DNA complexes . In the COMPEL database, several hundred composite elements are described with evidences that their regulatory functions differ from those provided by their individual components . As ARF transcription factors are known to form homo- and hetero-dimers, by sequence analysis of experimentally proven AuxREs we expected to find several AuxRE coupling motifs. Unlike the previous study of coupled regulatory motifs in promoters of auxin responsive genes [11, 41], we searched only for the composite elements with closely located coupling motifs (within 50 nt). As a result, we found and analyzed three distinct types of composite elements: AuxRE/Y-patch, AuxRE/ABRE-like and AuxRE/AuxRE-like (Figure 1D). Using microarray data analysis and functional annotation, we showed that among the genes having the composite elements in their upstream regions there were significantly higher portions of auxin regulated genes in comparison with those having a single potential AuxRE (Additional file 8).
AuxRE/ABRE-like composite element
The first discovered composite AuxRE named D1-4 consists of 11 nucleotides CCTCGTGTCTC, where 3'located TGTCTC overlaps by 3 nucleotides with 5'located reverse oriented ABA response element, ABRE . ABRE elements contain G-box (ACGA), which is bound by bZIP transcription factor EmBP-1 . EmBP-1 homologs specifically bind auxin-responsive D1-4 . In this work we encountered again on AuxRE/ABRE-like composite elements. The degenerative ABRE-like motif was discovered on the flanks of experimentally proven AuxREs in the regulatory regions of GH3, BRX, TMO7, IPT5 and DRN A. thaliana genes (Additional file 2). Functional annotation and expression level analysis for the genes with AuxRE/ABRE-like composite elements in the promoter and 5'UTR regions allowed us to specify the structure of functional AuxRE/ABRE-like. The composite element having the major impact on auxin responsiveness consists of direct AuxREP&S with 5' flanking ABRE-like with a long spacer (Figure 4B-C). In three microarrays, we also found a significant association of auxin response with the presence of the reversed AuxREs overlapped from the 5' side with reversed ABRE-like motif (Additional file 8). Despite, the result was insignificant, it appeared of interest, because the latter AuxRE/ABRE-like type resembles D1-4 composite element .
AuxRE/Y-patch composite element
The second AuxRE coupling motif recognized on the flanks of experimentally proven AuxREs was CT-rich motif similar to Y-patch (Additional file 3). Y-patch sequence was found significantly enriched around TSS in different plant genomes [21, 42]. Whereas TATA box is located within −45 to −18 in 25% of A.thaliana promoters, Y-Patch is located around TSS (from −50 to +50) in 49% of genes , which suggested their important role in transcription initiation . Here we found significant association between the presence of AuxRE/Y-patch composite element in gene upstream regions and its auxin responsiveness (Additional file 8).
The composite elements with the direct orientations of the coupling AuxRE and Y-patch were found significant by functional annotation and microarray data analysis (Figure 3A-B; 4A; Additional file 8). Both AuxRE and Y-patch were enriched around TSS, so the composite AuxRE/Y-patch might be just a consequence of their co-localization. However, the analysis of the relative positions and orientations of the AuxRE/Y-patch coupling elements suggests that they are not just located in the same region, but also need a specific structure, namely directly oriented Y-patch overlapped with the direct AuxRE on its 3' flank or located 3' with a short spacer (Figure 4). All the other combinatorial variants were found to be irrelevant to the auxin response.
Y-patch functionality has not been revealed yet , our data provide an idea that a protein, which binds to Y-patch, may interact with an ARF transcription factor and impact thus on auxin response.
AuxRE/AuxRE-like composite element
In auxin responsive gene promoters one usually may find multiple copies of potential AuxREs. The multimerized AuxREs were experimentally shown as providing auxin responsive gene expression . For example, the promoter of the widely used auxin sensor DR5 consists of 7 AuxREs . Recently, it has been shown that AuxRE-AuxRE reverted repeats are bound by ARF homodimers with a spacer-sensitive specificity . According to , 42% of auxin responsive promoters contain at least one pair of TGTC sites within 50 bp. Analysis of auxin responsive regions with multiple AuxREs suggested that they cooperate in a plant promoter providing for synergic auxin inductiveness . Here we found the AuxRE/AuxRE-like composite element consisting of AuxRE and its degenerative partner. The presence of this element in plant promoters significantly increase the feasibility for these genes to be involved in the plant developmental processes and to be auxin regulated, both was shown by functional annotation study (Additional file 12) and microarray data analysis (Additional file 8). We found that among AuxRE/AuxRE-like of the all possible structures, mainly those having direct AuxRE-like motif on the 3' flank of the direct AuxRE impact on auxin response (Additional file 12).
The AuxREs repeats can be the binding sites for ARF homodimers and heterodimers between different ARFs. They also may serve as the so called shadow enhancers. The shadow enhancers consisting of transcription factor binding site tandems, usually degenerative, were described as important cis-regulatory elements in a number of developmental genes in animals [44–46]. In the case of AuxRE/AuxRE-like composite elements, multiplication of weak sites into tandem clusters could make binding with ARF factors highly cooperative and strong.
The auxin signaling mechanisms in plants are very diverse, engaging many transcription factors and their co-regulators. Deciphering auxin responsive elements code may reveal some basic features in ARF-Aux/IAA machinery and also disclose additional players in auxin signaling. Here we presented a comprehensive computational study on distribution and annotation of simple and composite auxin responsive elements in A. thaliana. We revealed some associations in location, orientation and composite structure of potential AuxREs with auxin response.
Data sets for the bioinformatic analysis
The following datasets were used for AuxRE context analysis and recognition:
Training set. 25 experimentally proven AuxRE sites (Additional file 1) were collected from published data [5, 11, 17, 33, 37, 47–58]. Among them 19 were from Arabidopsis thaliana, 2 from Pisum sativum L., 2 from Glycine max L. and 1 from Oryza sativa L. and Withania somnifera L.. The sequences of the training set were aligned relative to TGTCNN with 50 nt flanks on both sides. The dataset was used to (1) deduce AuxREP&S; (2) search motifs on the flanks of the core hexamer.
Positive set contained promoter sequences of 44 auxin responsive genes, taken from .
Whole genome dataset and genomic annotation data for Arabidopsis thaliana were taken from the Plant Ensembl database (MySQL server, http://www.ensembl.org/info/data/mysql.html)
Whole-genome multispecies alignment for A. thaliana, A. lyrata, T. halophile, C. papaya and C. clementine was downloaded from http://pipeline.lbl.gov/downloads.shtml of Vista Tools for comparative genomics .
The analysis of experimentally proven AuxREs set was done with two approaches: Optimized PWM (oPWM) and SiteGA . The combination of both methods has been shown to be effective for large-scale genome transcription factors binding sites recognition . Optimized PWM (oPWM) utilizes positional weight matrix for regulatory elements analysis and discovery. The procedure of matrix optimization, which implies the search of the best matrix location and length, defined the optimum length for AuxRE sequence as 25 nt with the core hexamer located in the center (Figure 1A-B, Additional file 1).
The second approach, SiteGA  used the same dataset for training. SiteGA is based on the genetic algorithm (GA) involving a discriminant function of locally positioned dinucleotides (LPD). The chosen set of LPD defines the SiteGA recognition function and allows computing the distribution of significant dependencies along the site sequence (Figure 1C) as follows. Initially we selected a significance level (e.g. p < 0.05). Next, we tested all the possible pairs of LPDs and compiled only those pairs that conformed significant positive or negative correlations, i.e. significant dependencies.
Potential AuxRE sites were predicted separately by both oPWM and SiteGA methods at specified stringencies. The thresholds 0.78/0.936 for recognition methods oPWM/SiteGA were chosen using the positive set, such as upstream regions of 32% of genes contain at least one potential AuxRE. If potential AuxRE passed these thresholds, we defined it as an AuxREP&S variant.
The density of the site distribution was computed as the number of potential AuxREs per total number of tested positions in a sequence dataset. The abundance of potential AuxREs was estimated as the ratio of the AuxREs density in the analyzed region to that for the whole genome.
AuxRE coupling motifs search
Motif search in the neighborhood of the experimentally proven AuxREs was done with the MotiGA algorithm, an implementation of genetic algorithm, based on computation of the p-value for a definite weight matrix score  and the KDIC measure for frequency matrix . The description of MotiGA is given in Additional File 5. As described above, 106 nt long sequences from the training set were used for the coupling elements search. To avoid an influence of AuxRE consensus, the hexamer core TGTCNN was masked. For each overrepresented motif, the MotiGA method computed matrix of the frequencies and weights. Each sequence from the training set was assigned by the p-values and the PWM score for the coupling motif. Finally, the threshold of PWM was chosen by a critical p-value for each motif separately (Additional file 3).
Microarray data analysis
Sixteen microarray datasets on auxin treatment influences from published experiments (Table 2) were taken from the GEO database. The data were normalized using the RMA method and converted to log2 scale. The gene was considered as significantly up/down regulated by auxin, if the fold change between treatment and control was higher/lower than 1.5 and significant (p < 0.05) by the t-test (except GDS1408 which contained only one replica). Additionally, we required that at least 100 genes in each microarray were significantly up- or down-regulated. Associations between the presence of single or composite AuxRE in a gene promoter and its expression level were analyzed as follows. The portion of auxin up (down)-regulated genes observed in the experiment was compared with the respective portion observed for the subset of genes having a certain type of predicted AuxRE in a specified region. Namely, in this analysis we considered AuxREs recognized in the extended [-1500; 5'UTR] and proximal [-300; 5'UTR] regulatory regions. We considered an association as significant, if it was found at least in four microarray experiments (see below).
To confirm statistical significance of the differences calculated for the gene subsets, we applied the following t-test for proportions. The proportions were computed as the ratio of (a) the number of auxin regulated genes (significantly up- or down-regulated genes, genes from the positive set) to (b) the total number of genes in the set. The first proportion p 1 was computed for the set of genes, which had at least one potential AuxRE. The second proportion p 2 was calculated for the whole genome. In microarray analysis, p 2 was calculated for the genes which had a probe on the microchip. The angular (arcsine square-root) transformation z(p) was used to compute the t-test for two proportions, .
The threshold for the number of microarray experiments whereas we found statistically significant association was set by the binomial trial estimate . Namely, if we expect «success» in each of the sixteen tests with the probability p < 0.05, then we expect at least four «successes» in sixteen attempts with the probability P(16,4) < 0.01.
Only continuous fragments (more than 100 nt in length) of alignment from Vista tool  conserved for four plant species were taken in analysis. A potential AuxRE site was considered as conserved if it was located to the existing alignment and its core sequence TGTC was conserved among the four plant species. To estimate an average conservation of the regulatory regions we mapped the extended promoter regions [-1500; 5'UTR] to the existing alignments. Among all 21098 genes analyzed in microarray data, we defined only 15060 with any conserved nt in the regulatory region.
Functional annotation was performed with the use of the AgriGO GO analysis toolkit . For statistical analysis, we used the Fisher test method, for multiple hypothesis testing, we applied Hochberg FDR test  with a marginal p-value of 0.05. As input data we used a list of AGI codes of genes that contain certain elements in their regulatory regions. As the successful result, we took only those terms of functional annotation, which had undergone the mentioned above multiple testing.
Auxin Responsive Element
AuxRE recognized in the present work by both oPWM and SiteGA methods
Auxin Response Factor
Transcription Start Site
Pierre-Jerome E, Moss BL, Nemhauser JL: Tuning the auxin transcriptional response. J Exp Bot. 2013, 64: 2557-2563. 10.1093/jxb/ert100.
Sauer M, Robert S, Kleine-Vehn J: Auxin: simply complicated. J Exp Bot. 2013, 64: 2565-2577. 10.1093/jxb/ert139.
Ulmasov T, Hagen G, Guilfoyle TJ: ARF1, a Transcription Factor That Binds to Auxin Response Elements. Science. 1997, 276: 1865-1868. 10.1126/science.276.5320.1865.
Guilfoyle TJ, Ulmasov T, Hagen G: The ARF family of transcription factors and their role in plant hormone-responsive transcription. Cell Mol Life Sci. 1998, 54: 619-627. 10.1007/s000180050190.
Ulmasov T, Liu ZB, Hagen G, Guilfoyle TJ: Composite structure of auxin response elements. Plant Cell. 1995, 7: 1611-1623. 10.1105/tpc.7.10.1611.
Boer DR, Freire-Rios A, van den Berg WAM, Saaki T, Manfield IW, Kepinski S, López-Vidrieo I, Franco-Zorrilla JM, de Vries SC, Solano R, Weijers D, Coll M: Structural Basis for DNA Binding Specificity by the Auxin-Dependent ARF Transcription Factors. Cell. 2014, 156: 577-89. 10.1016/j.cell.2013.12.027.
Goda H, Sawa S, Asami T, Fujioka S, Shimada Y, Yoshida S: Comprehensive Comparison of Auxin-Regulated and Brassinosteroid-Regulated Genes in Arabidopsis. Plant Physiol. 2004, 134: 1555-1573. 10.1104/pp.103.034736.
Ulmasov T, Murfett J, Hagen G, Guilfoyle TJ: Aux/IAA Proteins Repress Expression of Reporter Genes Containing Natural and Highly Active Synthetic Auxin Response Elements. Plant Cell. 1997, 9: 1963-1971. 10.1105/tpc.9.11.1963.
Hagen G, Guilfoyle TJ: Auxin-responsive gene expression: genes, promoters and regulatory factors. Plant Mol Biol. 2002, 49: 373-385. 10.1023/A:1015207114117.
Liu ZB, Hagen G, Guilfoyle TJ: A G-Box-Binding Protein from Soybean Binds to the E1 Auxin-Response Element in the Soybean GH3 Promoter and Contains a Proline-Rich Repression Domain. Plant Physiol. 1997, 115: 397-407. 10.1104/pp.115.2.397.
Berendzen KW, Weiste C, Wanke D, Kilian J, Harter K, Dröge-Laser W: Bioinformatic cis-element analyses performed in Arabidopsis and rice disclose bZIP- and MYB-related binding sites as potential AuxRE-coupling elements in auxin-mediated transcription. BMC Plant Biol. 2012, 12: 125-10.1186/1471-2229-12-125.
Sun Y, Fan XY, Cao DM, Tang W, He K, Zhu JY, He JX, Bai MY, Zhu S, Oh E, Patil S, Kim TW, Ji H, Wong WH, Rhee SY, Wang ZY: Integration of brassinosteroid signal transduction with the transcription network for plant growth regulation in Arabidopsis. Dev Cell. 2010, 19: 765-77. 10.1016/j.devcel.2010.10.010.
Oh E, Zhu JY, Wang ZY: Interaction between BZR1 and PIF4 integrates brassinosteroid and environmental responses. Nat Cell Biol. 2012, 14: 802-9. 10.1038/ncb2545.
Oh E, Zhu J-Y, Bai M-Y, Arenhart RA, Sun Y, Wang Z-Y: Cell elongation is regulated through a central circuit of interacting transcription factors in the Arabidopsis hypocotyl. Elife. 2014, e03031-
Shin R, Burch AY, Huppert K a, Tiwari SB, Murphy AS, Guilfoyle TJ, Schachtman DP: The Arabidopsis transcription factor MYB77 modulates auxin signal transduction. Plant Cell. 2007, 19: 2440-2453. 10.1105/tpc.107.050963.
Varaud E, Brioudes F, Szécsi J, Leroux J, Brown S, Perrot-Rechenmann C, Bendahmane M: AUXIN RESPONSE FACTOR8 regulates Arabidopsis petal growth by interacting with the bHLH transcription factor BIGPETALp. Plant Cell. 2011, 23: 973-983. 10.1105/tpc.110.081653.
Walcher CL, Nemhauser JL: Bipartite promoter element required for auxin response. Plant Physiol. 2012, 158: 273-282. 10.1104/pp.111.187559.
Keilwagen J, Grau J, Paponov I, Posch S, Strickert M, Grosse I: De-novo discovery of differentially abundant transcription factor binding sites including their positional preference. PLoS Comput Biol. 2011, 7: e1001070-10.1371/journal.pcbi.1001070.
Levitsky VG, Ignatieva E V, Ananko EA, Turnaev II, Merkulova TI, Kolchanov NA, Hodgman TC: Effective transcription factor binding site prediction using a combination of optimization, a genetic algorithm and discriminant analysis to capture distant interactions. BMC Bioinformatics. 2007, 8: 481-10.1186/1471-2105-8-481.
Ratnere I, Dubchak I: Obtaining comparative genomic data with the VISTA family of computational tools. Curr Protoc Bioinformatics. 2009, 10: 10.6-
Yamamoto YY, Ichida H, Abe T, Suzuki Y, Sugano S, Obokata J: Differentiation of core promoter architecture between plants and mammals revealed by LDSS analysis. Nucleic Acids Res. 2007, 35: 6219-6226. 10.1093/nar/gkm685.
Parry TJ, Theisen JWM, Hsu J-Y, Wang Y-L, Corcoran DL, Eustice M, Ohler U, Kadonaga JT: The TCT motif, a key component of an RNA polymerase II transcription system for the translational machinery. Genes Dev. 2010, 24: 2013-2018. 10.1101/gad.1951110.
Bailey TL, Boden M, Buske FA, Frith M, Grant CE, Clementi L, Ren J, Li WW, Noble WS: MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res. 2009, W202-8. 37
Guiltinan MJ, Marcotte WR, Quatrano RS: A Plant Leucine Zipper Protein That Recognizes Abscisic Acid Response Element. Science. 1990, 250: 267-272. 10.1126/science.2145628.
Li Y, Liu ZB, Shi X, Hagen G, Guilfoyle TJ: An auxin-inducible element in soybean SAUR promoters. Plant Physiol. 1994, 106: 37-43. 10.1104/pp.106.1.37.
Zhou S, Zhou X, Ling Y, Zhang Z, Zhen D: agriGO: a GO analysis toolkit for the agricultural community. Nucleic Acids Res. 2010, 38: 64-70. 10.1093/nar/gkq310.
Guilfoyle TJ, Hagen G: Auxin response factors. Curr Opin Plant Biol. 2007, 10: 453-60. 10.1016/j.pbi.2007.08.014.
Szemenyei H, Hannon M, Long JA: TOPLESS mediates auxin-dependent transcriptional repression during Arabidopsis embryogenesis. Science. 2008, 319: 1384-1386. 10.1126/science.1151461.
Pfluger J, Zambryski P: The role of SEUSS in auxin response and floral organ patterning. Development. 2004, 131: 4697-4707. 10.1242/dev.01306.
Ulmasov T, Hagen G, Guilfoyle TJ: Activation and repression of transcription by auxin-response factors. Proc Natl Acad Sci. 1999, 96: 5844-5849. 10.1073/pnas.96.10.5844.
Paponov IA, Paponov M, Teale W, Menges M, Chakrabortee S, Murray JAH, Palme K: Comprehensive transcriptome analysis of auxin responses in Arabidopsis. Mol Plant. 2008, 1: 321-337. 10.1093/mp/ssm021.
Vernoux T, Brunoud G, Farcot E, Morin V, Van den Daele H, Legrand J, Oliva M, Das P, Larrieu A, Wells D, Guédon Y, Armitage L, Picard F, Guyomarc'h S, Cellier C, Parry G, Koumproglou R, Doonan JH, Estelle M, Godin C, Kepinski S, Bennett M, De Veylder L, Traas J: The auxin signalling network translates dynamic input into robust patterning at the shoot apex. Mol Syst Biol. 2011, 7: 508-
Schlereth A, Möller B, Liu W, Kientz M, Flipse J, Rademacher EH, Schmid M, Jürgens G, Weijers D, Mo B, Ju G: MONOPTEROS controls embryonic root initiation by regulating a mobile transcription factor. Nature. 2010, 464: 913-916. 10.1038/nature08836.
Riley T, Sontag E, Chen P, Levine A: Transcriptional control of human p53-regulated genes. Nat Rev Mol Cell Biol. 2008, 9: 402-412. 10.1038/nrm2395.
Omelina ES, Baricheva EM, Oshchepkov DY, Merkulova TI: Analysis and recognition of the GAGA transcription factor binding sites in Drosophila genes. Comput Biol Chem. 2011, 35: 363-370. 10.1016/j.compbiolchem.2011.10.008.
Okushima Y, Mitina I, Quach HL, Theologis A: AUXIN RESPONSE FACTOR 2 (ARF2): a pleiotropic developmental regulator. Plant J. 2005, 43: 29-46. 10.1111/j.1365-313X.2005.02426.x.
Inukai Y, Sakamoto T, Ueguchi-tanaka M, Shibata Y, Gomi K, Umemura I, Hasegawa Y, Ashikari M, Kitano H, Matsuoka M: Crown rootless1, Which Is Essential for Crown Root Formation in Rice, Is a Target of an AUXIN RESPONSE FACTOR in Auxin Signaling. Plant Cell. 2005, 17: 1387-1396. 10.1105/tpc.105.030981.
Kel O V., Romaschenko AG, Kel AE, Wingender E, Kolchanov NA: A compilation of composite regulatory elements affecting gene transcription in vertebrates. Nucleic Acids Res. 1995, 23: 4097-4103. 10.1093/nar/23.20.4097.
Makeev VJ, AP L, AG N, DA. P: Distance preferences in the arrangement of binding motifs and hierarchical levels in organization of transcription regulatory information. Nucleic Acids Res. 2003, 31: 6016-6026. 10.1093/nar/gkg799.
Heinemeyer T, Wingender E, Reuter I, Hermjakob H, Kel AE, Kel O V, Ignatieva E V, Ananko EA, Podkolodnaya OA, Kolpakov FA, Podkolodny NL, Kolchanov NA: Databases on transcriptional regulation: TRANSFAC, TRRD and COMPEL. Nucleic Acids Res. 1998, 26: 362-367. 10.1093/nar/26.1.362.
Nemhauser JL, Mockler TC, Chory J: Interdependency of brassinosteroid and auxin signaling in Arabidopsis. PLoS Biol. 2004, 2: E258-10.1371/journal.pbio.0020258.
Yamamoto YY, Yoshitsugu T, Sakurai T, Seki M, Shinozaki K, Obokata J: Heterogeneity of Arabidopsis core promoters revealed by high-density TSS analysis. Plant J. 2009, 60: 350-362. 10.1111/j.1365-313X.2009.03958.x.
Mironova VV, Omelyanchuk NA, Savina MS, Ponomarenko PM, Ponomarenko MP, Likhoshvai VA, Kolchanov NA: How multiple auxin responsive elements may interact in plant promoters: a reverse problem solution. J Bioinform Comput Biol. 2013, 11: 1340011-10.1142/S0219720013400118.
Papatsenko DA, Makeev VJ, Lifanov AP, Regnier M, Nazina AG, Desplan C: Extraction of Functional Binding Sites from Unique Regulatory Regions: The Drosophila Early Developmental Enhancers. Genome Res. 2002, 12: 470-481. 10.1101/gr.212502. Article published online before print in February 2002.
Lifanov AP, Makeev VJ, Nazina AG, Papatsenko : Homotypic Regulatory Clusters in Drosophila. Genome Res. 2003, 13: 579-588. 10.1101/gr.668403.
Gotea V, Visel A, Westlund JM, Nobrega MA, Pennacchio LA, Ovcharenko I: Homotypic clusters of transcription factor binding sites are a key component of human promoters and enhancers. Genome Res. 2010, 20: 565-77. 10.1101/gr.104471.109.
Ballas N, Wong LM, Ke M, Theologis A: Two auxin-responsive domains interact positively to induce expression of the early indoleacetic acid-inducible gene PS-IAA4/5. Proc Natl Acad Sci. 1995, 92: 3483-3487. 10.1073/pnas.92.8.3483.
Okushima Y, Overvoorde PJ, Arima K, Alonso JM, Chan A, Chang C, Ecker JR, Hughes B, Lui A, Nguyen D, Onodera C, Quach H, Smith A: Functional Genomic Analysis of the AUXIN RESPONSE FACTOR Gene Family Members in Arabidopsis thaliana: Unique and Overlapping Functions of ARF7 and ARF19. Plant Cell. 2005, 17: 444-463. 10.1105/tpc.104.028316.
Cole M, Chandler J, Weijers D, Jacobs B, Comelli P, Werr W: DORNROSCHEN is a direct target of the auxin response factor MONOPTEROS in the Arabidopsis embryo. Development. 2009, 136: 1643-1651. 10.1242/dev.032177.
Ismail IO: Function and Regulation of Xylem Cysteine Protease 1 and Xylem Cysteine Protease 2 in Arabidopsis. Dissertation. 2004, Virginia Polytechnic Institute
Donner TJ, Sherr I, Scarpella E: Regulation of preprocambial cell state acquisition by auxin signaling in Arabidopsis leaves. Development. 2009, 136: 3235-3246. 10.1242/dev.037028.
Scacchi E, Salinas P, Gujas B, Santuari L, Krogan N, Ragni L, Berleth T, Hardtke CS: Spatio-temporal sequence of cross-regulatory events in root meristem growth. Proc Natl Acad Sci. 2010, 107: 22734-22739. 10.1073/pnas.1014716108.
Zhao Z, Andersen SU, Ljung K, Dolezal K, Miotk A, Schultheiss SJ, Lohmann JU: Hormonal control of the shoot stem-cell niche. Nature. 2010, 465: 1089-1092. 10.1038/nature09126.
Hirota A, Kato T, Fukaki H, Aida M, Tasaka M: The auxin-regulated AP2/EREBP gene PUCHI is required for morphogenesis in the early lateral root primordium of Arabidopsis. Plant Cell. 2007, 19: 2156-2168. 10.1105/tpc.107.050674.
Zhu C, Perry SE: Control of expression and autoregulation of AGL15, a member of the MADS-box family. Plant J. 2005, 41: 583-594.
Sakamoto T, Fujioka S: Auxins increase expression of the brassinosteroid receptor and brassinosteroid-responsive genes in Arabidopsis. Plant Signal Behav. 2013, 8: e23509-10.4161/psb.23509.
Khan MR, Hu J, Ali GM: Reciprocal loss of CArG-boxes and auxin response elements drives expression divergence of MPF2-Like MADS-box genes controlling calyx inflation. PLoS One. 2012, 7: e42781-10.1371/journal.pone.0042781.
Cheng ZJ, Wang L, Sun W, Zhang Y, Zhou C, Su YH, Li W, Sun TT, Zhao XY, Li XG, Cheng Y, Zhao Y, Xie Q, Zhang XS: Pattern of auxin and cytokinin responses for shoot meristem induction results from the regulation of cytokinin biosynthesis by AUXIN RESPONSE FACTOR3. Plant Physiol. 2013, 161: 240-251. 10.1104/pp.112.203166.
Touzet H, Varré J-S: Efficient and accurate P-value computation for Position Weight Matrices. Algorithms Mol Biol. 2007, 2: 15-10.1186/1748-7188-2-15.
Kulakovskiy IV, Boeva VA, Favorov AV, Makeev VJ: Deep and wide digging for binding motifs in ChIP-Seq data. Bioinformatics. 2010, 26: 2622-2623. 10.1093/bioinformatics/btq488.
Levitsky VG: DNA motif search by genetic algorithm. Proceedings of the eighth International Conference on Bioinformatics of genome regulation and Structure\Systems Biology. 2012, Novosibirsk
Sokal RR, James RF: Biometry: the principles and practice of statistics in biological research. SERBIULA (sistema Libr 20). 2013
Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society. 1995, 1: 289-300. http://www.jstor.org/discover/10.2307/2346101?uid=368418201&uid=3738936&uid=368323351&uid=2&uid=3&uid=67&uid=62&sid=21104612809181,
Overvoorde PJ, Okushima Y, Alonso M, Chan A, Chang C, Ecker JR, Hughes B, Liu A, Onodera C, Quach H, Smith A: Functional Genomic Analysis of the AUXIN/INDOLE-3-ACETIC ACID Gene Family Members in Arabidopsis thaliana. Plant Cell. 2005, 17: 3282-3300. 10.1105/tpc.105.036723.
Redman JC, Haas BJ, Tanimoto G, Town CD: Development and evaluation of an Arabidopsis whole genome Affymetrix probe array. Plant J. 2004, 38: 545-61. 10.1111/j.1365-313X.2004.02061.x.
Stepanova AN, Yun J, Likhacheva A V, Alonso JM: Multilevel interactions between ethylene and auxin in Arabidopsis roots. Plant Cell. 2007, 19: 2169-85. 10.1105/tpc.107.052068.
Nagpal P, Ellis CM, Weber H, Ploense SE, Barkawi LS, Guilfoyle TJ, Hagen G, Alonso JM, Cohen JD, Farmer EE, Ecker JR, Reed JW: Auxin response factors ARF6 and ARF8 promote jasmonic acid production and flower maturation. Development. 2005, 132: 4107-18. 10.1242/dev.01955.
Vanneste S, De Rybel B, Beemster GTS, Ljung K, De Smet I, Van Isterdael G, Naudts M, Iida R, Gruissem W, Tasaka M, Inzé D, Fukaki H, Beeckman T: Cell cycle progression in the pericycle is not sufficient for SOLITARY ROOT/IAA14-mediated lateral root initiation in Arabidopsis thaliana. Plant Cell. 2005, 17: 3035-50. 10.1105/tpc.105.035493.
Armstrong JI, Yuan S, Dale JM, Tanner VN, Theologis A: Identification of inhibitors of auxin transcriptional activation by means of chemical genetics in Arabidopsis. Proc Natl Acad Sci USA. 2004, 101: 14978-83. 10.1073/pnas.0404312101.
Bargmann BOR, Vanneste S, Krouk G, Nawy T, Efroni I, Shani E, Choe G, Friml J, Bergmann DC, Estelle M, Birnbaum KD: A map of cell type-specific auxin responses. Mol Syst Biol. 2013, 9: 688-
We thank to Professor Charlie Hodgman, University of Nottingham for his encouragement in initiation of this project, Olga Chernikh for critical reading of the manuscript, Mikhail Ponomarenko and Vadim Efimov for valuable assistance in statistical analysis. The work was supported by the Russian Foundation for Basic Research [12-04-33112-mol-a-ved] to VVM, DSW and VGL; the Russian Science Foundation [14-14-00734] to NAO; Dynasty Foundation grant for young biologists to VVM.
Publication of this article has been funded by the Russian Science Foundation.
This article has been published as part of BMC Genomics Volume 15 Supplement 12, 2014: Selected articles from the IX International Conference on the Bioinformatics of Genome Regulation and Structure\Systems Biology (BGRS\SB-2014): Genomics. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcgenomics/supplements/15/S12.
The authors declare that they have no competing interests.
Conceived and designed the study: VVM, NAO and VGL. Performed the calculations: DSW and VGL. Analyzed the data: VVM and NAO. Wrote the paper: VVM, NAO and VGL. Coordinated the project: VVM. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 2: The frequencies of the dinucleotides that were used to construct oPWM model for AuxRE prediction. x axis on the figure denotes the position relative to the centrally located AuxRE core sequence, indicated by the capital letters in the frame. Y axis marks the dinucleotide frequencies, i.e. counts of the specific dinucleotides in the certain position of the sequence alignment. (TIF 8 MB)
Additional file 3: Description of Y-patch, AuxRE-like and ABRE-like coupling motifs, revealed by de novo motif search on the flanks of experimentally proven AuxREs. For each motif we compiled the frequency matrix with the estimated threshold; #seq denotes the number of the sequence in the Additional file 1; position, strand, PWM score, and sequences that passed the threshold. Logos constructed for these three sequence sets are shown on Figure 1D-F. (DOCX 14 KB)
Additional file 4: The significant correlations between the frequencies of the locally positioned dinucleotides (LPDs) for AuxREs from the Training set (Additional file 1). Correlations were deduced from the SiteGA model  for AuxRE recognition. Each horizontal strip depicts one correlation between two LPDs. A - positive correlations; B - negative correlations. The analyzed region consisted of 25 nt located [-12;+13] relative to the centrally positioned AuxRE core hexamer. (TIF 3 MB)
Additional file 6: The association of different AuxRE variants with auxin up-regulation. A. The portion of genes significantly up-regulated by auxin in more than three microarray experiments (Table 2) among the genes with a predicted AuxRE variant in [-1500; 5'UTR] regions. The basal level (dashed line) presents the average portion of auxin regulated genes in the whole genome. Namely, 1965 genes were significantly (>1,5 fold, p < 0,05) up-regulated in more than three microarrays (Table 1) among 21098 genes, which were detected by the microarray platform. B. The portion of up-regulated by auxin genes with conservative AuxRE among the up-regulated genes having an AuxRE. The plot shows that AuxREP&S related to auxin up-regulation are more conservative than other AuxRE variants. The basal conservation level denotes the portion of up-regulated genes which have any conservative nucleotide in the alignment of [-1500; 5'UTR] regions. Namely, 1442 out of 1965 genes. Statistics was calculated by t-test for proportions, * - p < 0.05; ** - p < 0.01. (TIF 799 KB)
Additional file 8: The summary table on microarray data analysis for the genes having in their regulatory regions predicted AuxRE variants. The portion of auxin up/down-regulated genes observed in the microarray was compared with the respective portion observed for the subset of genes having the certain variant of predicted AuxRE in a specified region. Significant increase in the portion by t-test for proportions in at least three microarrays was considered as significant association (see Methods for details). The specified regions are: [-1500; +5'UTR], [-300; +5'UTR]. The AuxRE variants are: (1) single/multiple; (2) predicted by SiteGA, oPWM, or by both methods (P&S), by TGTCTC or TGTSTSBC consensuses; (3) located in direct (+) or reverse (-) strand; (4) simple or composite; (5) composite AuxREs with a certain location of Y-Patch, ABRE-like or AuxRE-like relative to predicted AuxRE. (XLS 248 KB)
Additional file 9: The analysis of ABRE-like motif frequency matrix by the motif comparison tool TOMTOM . The search revealed the best match to MA0128.1 matrix for EmBP-1 (p < 0.005) is shown on the left, two logos on the right correspond to MA0128.1 and ABRE-like motif (HDSDYKS, see Additional file 3 for the frequency matrix). (TIF 36 KB)
Additional file 11: Functional annotation of the genes with predicted composite element AuxRE+/ ABRE-like in their regulatory [-1500; 5'UTR] region. Functional annotation was performed using the singular enrichment analysis from AgriGO analysis tool  under the Hochberg FDR multitest adjustment method. Each node in the graph represents a definite GO term with the value of FDR in the brackets. The yellow nodes represent significantly enriched terms, p < 0.002 (see annotation of the Figure 3). (TIF 1 MB)
Additional file 12: Functional annotation of the genes with predicted composite element AuxRE+/AuxRE-like+ in their regulatory [-1500; 5'UTR] region. Functional annotation was performed using the singular enrichment analysis from AgriGO analysis tool  under the Hochberg FDR multitest adjustment method. Each node in the graph represents a definite GO term with the value of FDR in the brackets. The yellow nodes represent significantly enriched terms, p < 0.002 (see annotation of the Figure 3). (TIF 240 KB)
About this article
Cite this article
Mironova, V.V., Omelyanchuk, N.A., Wiebe, D.S. et al. Computational analysis of auxin responsive elements in the Arabidopsis thaliana L. genome. BMC Genomics 15, S4 (2014) doi:10.1186/1471-2164-15-S12-S4
- Plant hormone
- primary response
- transcription regulation
- composite element