<title>Background</title> <p>DNA sequence motifs that recruit the transcription factors necessary to regulate the expression of a gene, are most commonly found in the flanking DNA regions and provide specificity to the core transcriptional machinery [##REF##12777501##1##]. In plants with annotated whole genome sequence such as Arabidopsis [##REF##11130711##2##], flanking DNA sequences upstream of the coding region can easily be defined. Such sequences are commonly referred to as the promoter and while they can be difficult to delineate in the absence of experimental characterisation, they can be defined as the intergenic sequence upstream of the ATG, and often limited to a defined length eg. 3 kb [##REF##16524982##3##]. In this definition the promoter fragment includes the 5' untranslated region (5'UTR).</p> <p>DNAse I footprinting [##REF##212715##4##] and electrophoretic or gel mobility shift assays [##REF##6269071##5##] have been extensively used to characterise <italic>cis</italic>-regulatory elements. Both methods rely on the direct interaction between DNA fragments that contain the DNA-binding region and the corresponding transcription factor. More recently, ChIP-microarray (also known as ChIP-chip) has been used to immunoprecipitate DNA associated with a TF of interest. The DNA from this complex is then used to probe a genomic DNA microarray [##UREF##0##6##,##REF##14986705##7##]. Studies which have used ChIP to identify TF binding sites include the analysis of the <italic>AGAMOUS </italic>[##REF##15634696##8##], <italic>AGL15</italic>, [##REF##12743119##9##] and the <italic>FLOWERING LOCUS C PROTEIN </italic>(<italic>FLC</italic>) [##REF##16623882##10##] MADS box genes from Arabidopsis, all of which have been shown to bind to a CArG box contained in the promoter of the target gene.</p> <p>Surface plasmon resonance (SPR) is an emerging technology that allows the characterisation of protein DNA interactions <italic>in-vitro </italic>[##REF##17964928##11##]. Importantly, this technique allows an assessment of DNA-protein kinetics, affinity and specificity in real time. A number of plant TF binding sites have been investigated using this technique including, ZPT2-2 from petunia [##REF##11454859##12##] and VRN-1 from Arabidopsis [##REF##12114624##13##].</p> <p>Transcription factor-DNA interactions do not however infer transcriptional activation, for example the <italic>Antirrhinum </italic>MYB305 protein has been shown to bind the CHS promoter in gel-shift analysis but failed to induce transcriptional activation of the gene in yeast-1-hybrid assays [##REF##8837506##14##]. In addition, these experimental approaches rely on the need to purify TF protein beforehand and will only reflect <italic>in-vitro </italic>binding. Often these associations require co-factors or additional transcription factors that facilitate the interaction of a protein to its <italic>cis</italic>-regulatory regions [##REF##15361138##15##,##REF##15708343##16##]. Alternatively, yeast-1-hybrid assays determine protein-DNA interactions through transcriptional activation of several reporter genes: <italic>HIS3</italic>, <italic>URA3 </italic>and <italic>LEU2 in-vivo </italic>[##REF##8266075##17##,##REF##8321284##18##]. While these assays are effective at analysing simple protein-DNA interactions, the absence of any plant-derived factors other than the TF under investigation, can limit the applicability of this technique. The limited number of well characterised TF binding sites highlights the difficulty in adopting these approaches for large-scale characterisation of <italic>cis-</italic>regulatory sites in plants [##REF##16524982##3##].</p> <p>Whilst there are relatively few confirmed <italic>cis</italic>-regulatory sites in relation to the number of known transcription factors, a number of TF classes have consensus binding sites proposed. The bZIP class of transcription factors have been shown to preferentially bind palindromic sites such as the G-box (CACGTG) [##REF##11722549##19##], A box (TACGTA) and C-box (GACGTC) [##REF##8119490##20##,##REF##18287490##21##]. Several plant MADS box genes including; <italic>DEFICIENS </italic>(<italic>DEF</italic>) and <italic>GLOBOSA </italic>(<italic>GLO</italic>) from <italic>Antirrhinum </italic>[##REF##1346760##22##], <italic>APETALA-1 </italic>(<italic>AP1</italic>), <italic>APETALA-3 </italic>(<italic>AP3</italic>), <italic>PISTILLATA </italic>(<italic>PI</italic>) [##REF##8643482##23##] and <italic>AGAMOUS (AG</italic>) [##REF##8106084##24##] in Arabidopsis have been shown to bind variations of a CArG motif, and a consensus CArG sequence has been described as CC(A/T)<sub>6</sub>GG [##REF##16616581##25##]. <italic>LEAFY (LFY) </italic>controls the switch from vegetative to reproductive development in Arabidopsis [##REF##1350515##26##] and interacts with the consensus <italic>LFY </italic>binding site (CCANTG) to activate <italic>AP1 </italic>in the meristem identity pathway and the floral homeotic <italic>AG </italic>gene [##REF##10417387##27##]. The WRKY TF class has been implicated in responses such as pathogen defence, senescence and trichome development [##REF##10785665##28##] and bind to a conserved W box TTTGAC(C/T) motif contained in their respective target promoters [##REF##7969025##29##, ####REF##8541496##30##, ##REF##8896462##31##, ##REF##8972846##32####8972846##32##]. MYB transcription factors regulate a diverse range of pathways including secondary metabolism, signal transduction and defence responses [##REF##11597504##33##]. Two MYB binding site sequence variants described in plants by Romero <italic>et al</italic>. (1998) are the type II (GTT(A/T)GTT(G/A) and IIG G(G/T)T(A/T)GGT(G/A) sites common to a number of genes in the phenylpropanoid pathway [##REF##9628022##34##]. A third conserved sequence (A/C)ACC(A/T)A(A/C)C, has been shown to be bound by the flavonoid regulator MYB305 from <italic>Antirrhinum majus </italic>[##REF##8306956##35##].</p> <title>PAP and anthocyanin biosynthesis</title> <p><italic>PRODUCTION OF ANTHOCYANIN PIGMENT 1 (PAP1</italic>) is an R2–R3 MYB gene from Arabidopsis that is responsible for the co-ordinated up-regulation of genes in the anthocyanin pathway [##REF##11148285##36##]. The anthocyanin biosynthetic pathway has been well characterised at both the biochemical and regulatory level. While over-expression of single enzyme components of the flavonoid pathway does not significantly alter the amount of anthocyanin in plants; over-expression of the <italic>PAP1 </italic>gene activates components of the biosynthetic pathway enabling increases in anthocyanin accumulation [##REF##11148285##36##].</p> <p>Microarray studies of transgenic Arabidopsis over-expressing <italic>PAP1 </italic>identified a list of 38 genes that were selected as significantly changing in expression [##REF##15807784##37##]. This study examined constitutive expression effects in a mature plant so some genes with an altered expression profile may be due to indirect effects of gene over-expression, such as alteration in cell physiology or metabolite partitioning in response to the increase in anthocyanins. Or the result of other transcription factors that activate the expression of a different set of genes.</p> <p>Here we describe a novel method that used transient assays to identify and validate a <italic>cis</italic>-regulatory motif that is necessary for transactivation by PAP1. Candidate genes were selected from microarray analysis of <italic>PAP1 </italic>over-expressing transgenic plants. We identified a sub-set of promoters that were directly transactivated by PAP1 and used this information to identify a sequence motif that was conserved within the promoter regions of these unrelated genes. Deletion and mutation of this candidate <italic>cis</italic>-regulatory element in two promoters led to significant reductions in the level of transactivation by PAP1. Taken together our results demonstrate that validation of microarray data by transactivation assays provides a powerful way of elucidating conserved motifs within co-ordinately regulated genes.</p>
<title>Materials and methods</title> <title>Plant material and growth conditions</title> <p>The 35S-<italic>PAP1 </italic>construct was generated by inserting a genomic clone of the Arabidopsis <italic>PAP1 </italic>(At5g56650) gene into a nos-kanamycin containing vector pGreenII 0029-62-SK as previously described in Hellens <italic>et al. </italic>(2005) [##REF##16359558##38##]. Constructs were electroporated into <italic>Agrobacterium tumefaciens </italic>GV3101 (MP90) then transformed into <italic>Arabidopsis thaliana </italic>col-1 plants using the floral dip method [##REF##10069079##48##]. These plants, and vector-only controls, were grown together in either a greenhouse under short day conditions (8 h light/16 h dark, 21°C) or a growth room (constant light, 25°C). For the transient assays, <italic>Nicotiana benthamiana </italic>plants were grown, and transient leaf assays carried out as described in Hellens <italic>et al. </italic>(2005) [##REF##16359558##38##]. The LUC/REN ratio was used to quantify promoter activity and is a measure of luciferase expression relative to the expression of 35S-<italic>Renilla </italic>also contained on the same reporter plasmid. Background levels of promoter activity were assessed using only the promoter-LUC-35S-REN constructs (no transcription factor) [##REF##16359558##38##].</p> <title>Microarray analysis</title> <p>RNA was extracted from seedling and mature Arabidopsis plants according to Chang <italic>et al. </italic>(1993) [##UREF##1##49##]. RNA was quantified for integrity and concentration using a 2100 BioAnalyzer (Agilent technologies). RNA was labelled with Cy 3 and Cy 5 fluorescent dyes (GE Healthcare) as previously described [##REF##17556515##50##]. All analysis compared 35S-<italic>PAP1 </italic>plants with plants containing vector only. Each condition was repeated twice with a dye swap comparison for each repeated sample (4 arrays).</p> <p>Arabidopsis full genome 27 K oligo microarrays (Operon) were spotted onto epoxy coated slides (MWG) in a 150 mM phosphate buffer, pH 8.5, using a Biorobotics MicroGrid robot and Biorobotics 100 μM pins. Microarrays were hybridised as previously described [##REF##17556515##50##] except the 16-hour hybridisations were carried out at 60°C rather than 45°C. Arrays were scanned using a Genepix 4000 scanner and spots were aligned using Genepix 5 software. All data were processed in R using the Bioconductor limma package [##UREF##2##51##]. Genes were selected as significant using a False Discovery Rate (FDR) of 0.05 [##UREF##3##52##].</p> <title>Promoter cloning and plasmid constructs</title> <p>Promoter sequences were defined according to TIGR 6.0 annotation of the Arabidopsis genome. A 1 kb upstream fragment and the 5'UTR, where present, was amplified by two oligonucleotide primers, one which flanked the ATG start codon and one 1 kb upstream (Additional file ##SUPPL##1##2##). The primers introduced <italic>Xma </italic>I and <italic>Not </italic>I restriction sites into the amplification product respectively, to facilitate directional cloning. Promoter fragments were cloned into a pGem-T easy (Promega Madison, WI) and directionally subcloned into a pGreenII-0800-LUC [##REF##16359558##38##] using the <italic>Xma</italic>I and <italic>Not</italic>I restriction sites and verified by sequencing.</p> <p>Motif deletions and mutations were created by designing divergent PCR primers that flanked or spanned the predicted motifs in At5g17220 and At4g14090 promoters (Additional file ##SUPPL##2##3##). PCR was performed on the corresponding pGem-T easy clone of the promoter fragments using Prime Star polymerase (Takara Shiga, Japan). Blunt-ended PCR products were phosphorylated with 1 mM ATP, 10U T4 Polynucleotide Kinase (New England Biolabs Ipswich MA), and 1× Polynucleotide Kinase Buffer for 1 h at 37°C then re-ligated using the Rapid DNA ligation kit (Roche Mannheim Germany) for 2 h at room temperature to recreate the vector. Modified promoters were sequence verified and directionally subcloned as above.</p> <title>Identification of PAP1 cis-regulatory elements</title> <p>Conserved motifs were identified using the MEME motif search programme [##REF##7584402##42##], with default variables of the following parameters: 1) Any number of repetitions of motif per sequence, 2) motif length min = 6 bp, max = 10 bp, 3) maximum of 5 motifs searched. Only motifs that were represented a least once in each promoter were considered as potential PAP1 <italic>cis</italic>-regulatory elements. The motif search programme COSMO [##REF##17402923##43##] was also used to identify conserved motifs. Default variables were used with motif length min = 6 bp and max = 10 bp.</p>
<title>Results and Discussion</title> <title>Selection of differentially expressed genes resulting from the over-expression of <italic>PAP1</italic></title> <p>A 35S-<italic>PAP1 </italic>over-expression cassette was used to generate a number of stable transgenic Arabidopsis lines with varying levels of anthocyanin accumulation. These plants showed elevated levels of anthocyanin accumulation in all plant parts when grown under standard conditions. The level of expression of the transgene correlated strongly with the levels of anthocyanin. One line homozygous for the transgene and with consistently high levels of anthocyanins was chosen for this study (Fig. ##FIG##0##1A##). Microarray analysis of labelled RNA from this <italic>PAP1 </italic>over-expression line was compared with transgenic lines that contained a vector-only control construct [##REF##16359558##38##]. Four microarrays were hybridised with RNA from seedlings, and four microarrays were hybridised with labelled RNA from mature plants (two biological replicates repeated in a dye swap). Using a false rate discovery threshold of 0.01, 1744 genes were identified that showed significantly different expression levels, in both seedlings and mature plants, between the plants containing the 35S-<italic>PAP1 </italic>gene and those that contained the 35S control construct only (Additional file ##SUPPL##0##1##).</p> <p>We do not believe that all the gene expression changes seen in both our microarray analysis and those previously published were direct targets of the <italic>PAP1 </italic>gene. It is likely that pleiotropic expression changes will arise from effects such as alterations in the cell physiology and downstream regulation by transcription factors. A total of 35 genes on the array list were annotated as regulatory genes, consistent with the hypothesis that many gene expression changes observed are the result of the secondary effects.</p> <title>A subset of up-regulated genes are also transactivated in leaf infiltration assays</title> <p>From the 1744 genes that significantly changed in transcript level due to over-expression of the <italic>PAP1 </italic>gene, 33 were selected for further investigation; based on those genes that showed the greatest increase in expression compared with the control vector. The gene set comprised 17 that increased only in mature plants, 7 that increased only in seedlings and 9 that changed in both tissue types. A previous study using Affymetrix Arabidopsis genome arrays, identified a subset of 39 genes that were up-regulated in response to <italic>PAP1 </italic>over-expression [##REF##15807784##37##]. Of these 39 genes, 7 were also selected in our gene set. The 33 promoter fragments were cloned into a transactivation assay vector, pGreen 0800-LUC [##REF##16359558##38##] (Fig. ##FIG##0##1B## and ##FIG##0##1C##). To simplify the cloning process and to minimise PCR induced errors we used 1 kb of upstream sequence plus the 5' UTR (where annotated) to create the promoter constructs. <italic>Agrobacterium </italic>containing the cloned promoter constructs were infiltrated into the leaves of <italic>N. benthamiana </italic>either with or without <italic>Agrobacterium </italic>containing the 35S-<italic>PAP1 </italic>cassette used to generate the transgenic plants (Fig. ##FIG##1##2A##). From this initial screening, eight promoters showed a statistically significant increase in relative luciferase (LUC) activity when co-infiltrated with the 35S-<italic>PAP1 </italic>cassette, compared with the promoter-only controls (Fig. ##FIG##1##2A##). While most of the promoters gave a low level of relative LUC activity in the absence of the <italic>PAP1 </italic>gene, the promoter for At4g09820 <italic>TRANSPARENT TESTA 8 </italic>(<italic>TT8</italic>) gave a high level of activity in the absence of <italic>PAP1</italic>. Interestingly, TT8 is a bHLH transcription factor that is an important co-factor in the regulation of both anthocyanins and condensed tannins [##REF##11041882##39##]. The result here implies that, in our tobacco assay, the <italic>TT8 </italic>promoter fragment is relatively active transcriptionally.</p> <p>Eight PAP1 transactivated promoters and two non-responsive promoters were re-assayed and the transactivation confirmed (Fig. ##FIG##1##2B##). Of the eight trans-activated promoters identified by our analysis, six were from genes whose encoded proteins have a role in the anthocyanin and proanthocyanin biosynthetic pathways. In addition to these, a gene corresponding to a lipid transfer protein precursor (At5g59310) and a MYB transcription factor (At1g66380) were also identified in our microarray experiment and confirmed in transactivation assays. Lipid transfer proteins are a class of small basic soluble proteins capable of binding fatty acids and acyl CoA esters [##REF##15012303##40##]. As malonyl-CoA is an early precursor of the anthocyanin biosynthetic pathway [##REF##15608338##41##], it is possible that lipid transfer proteins may be acting in a transport role or as co-factors for the conversion of these intermediates to anthocyanins. The transient assay data also indicate that PAP1 was able to transactivate the transcription of a second MYB like gene <italic>MYB114 </italic>(At1g66380). In Arabidopsis <italic>MYB114 </italic>belongs to a tandem repeat with two other MYB genes, <italic>MYB113 </italic>and <italic>MYB90 </italic>(<italic>PAP2</italic>). All these genes show significant sequence similarity to <italic>PAP1 </italic>(<italic>MYB75</italic>), although only the <italic>PAP1 </italic>and <italic>PAP2 </italic>genes have been reported to regulate anthocyanin biosynthesis. This observation supports a potential role for the <italic>PAP1 </italic>gene in a feed-forward regulation of at least one related MYB gene.</p> <p>One promoter, for the gene encoding dihydroflavanol reductase (<italic>DFR</italic>; At5g42800), showed a 122-fold elevation in relative LUC activity in the transactivation assay. Three promoters from the genes encoding glutathione S-transferase (<italic>GST</italic>; At5g17220), leucoanthocyanidin dioxygenase (<italic>LDOX</italic>; At4g22880) and UDP flavonoid 5-O-glycosyltransferase (<italic>UFGT</italic>; At4g14090) showed a 38- to 60-fold elevation. The remaining four promoters derived from the genes encoding chalcone synthase (<italic>CHS</italic>; At5g13930), flavanone 3-hydroxylase (<italic>F3H</italic>; At3g51240), <italic>MYB114 </italic>(At1g66380) and a non-specific lipid transfer protein precursor (<italic>LTP</italic>; At5g59310), showed a PAP1-dependent elevation of 3- to 7-fold (Table ##TAB##0##1##).</p> <title>Motif analysis of transactivated promoter sequence</title> <p>The naive motif search programme MEME [##REF##7584402##42##], was used to search the DNA sequences for sequence motifs that were common to all eight promoters transactivated by the <italic>PAP1 </italic>gene. With the default settings and a maximum of five output motifs, only one 10 bp motif was present in all eight of the transactivated promoters (Table ##TAB##0##1##). This conserved motif was found in both plus and minus orientations and upstream of the 5' UTR where annotation is available. This motif was absent (<italic>P</italic>-value &lt; 0.119) from the 25 promoter fragments that were initially screened and did not alter relative LUC activity in the presence of the <italic>PAP1 </italic>gene. From these predictions, we hypothesise that this conserved motif (C/T)CNCCAC(A/G)(A/T)(G/T) is a PAP1 <italic>cis</italic>-regulatory element (PCE). Searches performed on the same promoter set using the search program COSMO [##REF##17402923##43##] yielded a related motif with a common core and related flanking sequences (C/T)(A/C)NCCACN(G/T)(G/T). When MEME analysis was conducted on the top 10 up-regulated genes from both this study and those previously published, neither the PCE nor any other motif was identified. This demonstrates the benefit of using only direct targets identified in the transient assay.</p> <title>Fold change and PCE frequency</title> <p>The level of transactivation from the co-infiltration of 35S-<italic>PAP1 </italic>varied between promoter-LUC cassettes (Table ##TAB##0##1##). Notably the promoter with the highest transactivation values At5g42800 (<italic>DFR</italic>), which showed a 122-fold increase in luciferase activation when co-infiltrated with 35S-<italic>PAP1</italic>, contains two PCEs within the promoter region used in the assay. Three promoters showed between 38- to 60-fold increases in relative LUC activity (<italic>GST</italic>; At5g17220, <italic>LDOX</italic>; At4g22880 and <italic>UFGT</italic>; At4g14090) and contained only a single PCE. Four other promoters (<italic>CHS</italic>; At5g13930, <italic>LTP</italic>; At5g59310, <italic>MYB114</italic>; At1g66380 and <italic>F3H</italic>; At5g51240) had much smaller 3- to 7-fold increases in relative LUC activity. This lower activation may be explained by the C to T change in the first base of the highly conserved CCAC core of the PCE motif in three of the promoters with the lowest transactivation values. However, this does not explain the data from the At3g51240 promoter, which had a low transactivation value and a fully conserved PCE motif (<italic>P</italic>-value = 5.56e<sup>-08</sup>). As the levels of expression from the transiently infiltrated 35S-<italic>PAP</italic>1 gene were higher than under normal physiological conditions, an alternative explanation for these lower transactivation values is that high levels of <italic>PAP1 </italic>expression may result in non-specific binding to low affinity sites in the promoter. However, as the 35S-<italic>PAP1 </italic>and promoter-LUC fusion were infiltrated in a ratio of 9:1, these transient leaf infiltration assays may more closely resemble the physiological ratio of TF to promoter than the over-expression of a TF in transgenic plants.</p> <title>Validation of PCE by transient leaf infiltration</title> <p>We tested the integrity of the PCE by deleting and mutating the sequence from the At5g17220 (<italic>GST</italic>) and At4g14090 (<italic>UFGT</italic>) promoter fragments. PCE promoter deletions had the 10 bp motif excised from the At5g17220 and At4g14090 promoter fragments. These were labelled At5g17220D and At4g14090D respectively. The sequence from At4g42800 (<italic>DFR</italic>) was not chosen for deletion or mutation analysis due to the presence of two PCEs which complicates base modification by PCR. Promoter mutations were generated by altering four bases of the PCE at positions 1 (C to T), 4 (C to T), 7 (C to T) and 10 (G to A) in the At5g17220 promoter. In promoter At4g14090, six base changes were made at positions 1 (C to T), 4 (C to T), 5 (C to A), 7 (C to T) 8 (A to T) and 10 (G to T). These modified promoters were labelled At5g17220M and At4g14090M respectively. Transactivation assays were used to compare these modified promoters with the unmodified At5g17220 and At4g14090 promoters (Fig. ##FIG##2##3##). At4g17220D had a significantly reduced LUC/REN ratio compared with the wild type sequence in the presence of 35S-<italic>PAP1</italic>. At4g17220M also showed a reduction in LUC/REN ratio in the presence of 35S-<italic>PAP1</italic>, compared with the wild type sequence. Modification of the PCE from At4g14090 showed similar reductions in relative LUC/REN ratio for At4g14090D and At4g140909M compared with the wild type. The six base changes in the PCE of At5g17220M appeared to have more of an effect than the 10 base changes introduced into the PCE region of the At4g14090 promoter. This may be due to different base substitutions at position 10, (At5g17220; G to A change, At4g14090; G to T change). A more detailed analysis of each of the conserved nucleotides within the PCE will be needed to fully interpret the significance of these results.</p> <p>Interestingly, five of the PCE-containing promoters identified in this study also contained a perfect G-box site (CACGTG) adjacent to the PCE sequence. A number of plant promoters regulated by diverse signals contain G-box elements [##REF##8127687##44##]. At least 2 classes of TFs are capable of binding G-boxes: the basic leucine zipper class (bZIP) and the basic helix-loop-helix proteins (bHLH) [##REF##11722549##19##]. Extensive genetic and protein studies have shown a close functional relationship between MYBs regulating anthocyanin accumulation and bHLH proteins [##REF##15361138##15##]. In maize, the activation of anthocyanin biosynthetic genes by ZmPl and ZmC1 requires a bHLH protein encoded by a <italic>R/B </italic>gene [##REF##1577278##45##]. In Arabidopsis the bHLH encoding <italic>TT8 </italic>gene and a MYB TF encoded by the <italic>TT2 </italic>gene, act synergistically to direct the expression of the <italic>DFR </italic>and <italic>BANYULS </italic>(<italic>BAN</italic>) flavonoid pathway genes [##REF##11041882##39##]. PAP1 has been demonstrated to interact with the bHLH proteins encoded by <italic>ENHANCER OF GLABRA3 </italic>(<italic>EGL3</italic>) and <italic>GLABRA3 </italic>(<italic>GL3</italic>) genes, and when co-overexpressed these combinations showed far more severe phenotypes than would be expected for additive regulation alone [##REF##12917293##46##]. The presence of the PCE and the G-box may not be coincidental and may not correspond to the binding site of PAP1, but that of the bHLH gene that we assume to be necessary for transactivation. In these transient assays we presume the appropriate endogenous bHLH protein interacts with the transient expressed <italic>PAP1 </italic>gene product.</p> <title>Occurrence of the PCE in 300 genes from microarray results</title> <p>We calculated probability values for putative PCEs occurring in the 2 kb upstream of the ATG of the top 300 up or down regulated genes selected from the combined mature plant and seedling microarrays, using MAST [##REF##9520501##47##]. This combined analysis did not select one of the eight (At4g14090) that was transactivated in the transient assay due to it having low difference of expression in seedlings. MAST is a search program which scans input sequences for known motifs and calculates match scores for each input sequence. From the MAST output a combined <italic>P</italic>-value can be obtained which measures the strength of the match of the sequence to the input motif. These combined <italic>P</italic>-values were plotted against the log-fold change calculated from the microarray results (Fig. ##FIG##3##4##). It was found that of the 6 genes that showed the biggest up-regulation of expression in the 35S-<italic>PAP1 </italic>plants, all had good (&lt;0.119) <italic>P</italic>-values for the PCE and were transactivated in the transient assay. When a less than log 2 fold up-regulation was observed, the occurrence of a PCE and transactivation became less easy to predict. Of the 7 tested genes that showed transactivation, the combined MAST <italic>P</italic>-value ranged from 2.8e<sup>-3 </sup>to 0.119 suggesting that in this range, transactivation could occur. One of the 25 tested promoters (At5g24770; <italic>VEGETATIVE STORAGE PROTEIN 2</italic>) that was not transactivated by the <italic>PAP1 </italic>gene did have a putative PCE (CCATCACAAG), but only with a <italic>P</italic>-value of 0.199. This suggests that this level of identity to the PCE consensus was insufficient to activate the gene. One additional promoter that was not transactivated (At5g05270, <italic>CHALCONE-FLAVANONE ISOMERASE</italic>) contained a PCE within this <italic>P-</italic>value range but located outside the 1 kb region tested for transactivation. Of the top 300 genes that showed changes in expression in the presence of the <italic>PAP1 </italic>transgene, a further 18 untested genes were significantly up-regulated and had a PCE within the transactivation <italic>P</italic>-value range, suggesting that these may also be regulated by PAP1 (Fig. ##FIG##3##4##). In addition to these 18 genes, there were 13 genes (4%) that had a <italic>P-</italic>value within this range, but were down-regulated. This either implies a repressor function or that the presence of the PCE alone is not sufficient to cause activation of these genes.</p>
<title>Results and Discussion</title> <title>Selection of differentially expressed genes resulting from the over-expression of <italic>PAP1</italic></title> <p>A 35S-<italic>PAP1 </italic>over-expression cassette was used to generate a number of stable transgenic Arabidopsis lines with varying levels of anthocyanin accumulation. These plants showed elevated levels of anthocyanin accumulation in all plant parts when grown under standard conditions. The level of expression of the transgene correlated strongly with the levels of anthocyanin. One line homozygous for the transgene and with consistently high levels of anthocyanins was chosen for this study (Fig. ##FIG##0##1A##). Microarray analysis of labelled RNA from this <italic>PAP1 </italic>over-expression line was compared with transgenic lines that contained a vector-only control construct [##REF##16359558##38##]. Four microarrays were hybridised with RNA from seedlings, and four microarrays were hybridised with labelled RNA from mature plants (two biological replicates repeated in a dye swap). Using a false rate discovery threshold of 0.01, 1744 genes were identified that showed significantly different expression levels, in both seedlings and mature plants, between the plants containing the 35S-<italic>PAP1 </italic>gene and those that contained the 35S control construct only (Additional file ##SUPPL##0##1##).</p> <p>We do not believe that all the gene expression changes seen in both our microarray analysis and those previously published were direct targets of the <italic>PAP1 </italic>gene. It is likely that pleiotropic expression changes will arise from effects such as alterations in the cell physiology and downstream regulation by transcription factors. A total of 35 genes on the array list were annotated as regulatory genes, consistent with the hypothesis that many gene expression changes observed are the result of the secondary effects.</p> <title>A subset of up-regulated genes are also transactivated in leaf infiltration assays</title> <p>From the 1744 genes that significantly changed in transcript level due to over-expression of the <italic>PAP1 </italic>gene, 33 were selected for further investigation; based on those genes that showed the greatest increase in expression compared with the control vector. The gene set comprised 17 that increased only in mature plants, 7 that increased only in seedlings and 9 that changed in both tissue types. A previous study using Affymetrix Arabidopsis genome arrays, identified a subset of 39 genes that were up-regulated in response to <italic>PAP1 </italic>over-expression [##REF##15807784##37##]. Of these 39 genes, 7 were also selected in our gene set. The 33 promoter fragments were cloned into a transactivation assay vector, pGreen 0800-LUC [##REF##16359558##38##] (Fig. ##FIG##0##1B## and ##FIG##0##1C##). To simplify the cloning process and to minimise PCR induced errors we used 1 kb of upstream sequence plus the 5' UTR (where annotated) to create the promoter constructs. <italic>Agrobacterium </italic>containing the cloned promoter constructs were infiltrated into the leaves of <italic>N. benthamiana </italic>either with or without <italic>Agrobacterium </italic>containing the 35S-<italic>PAP1 </italic>cassette used to generate the transgenic plants (Fig. ##FIG##1##2A##). From this initial screening, eight promoters showed a statistically significant increase in relative luciferase (LUC) activity when co-infiltrated with the 35S-<italic>PAP1 </italic>cassette, compared with the promoter-only controls (Fig. ##FIG##1##2A##). While most of the promoters gave a low level of relative LUC activity in the absence of the <italic>PAP1 </italic>gene, the promoter for At4g09820 <italic>TRANSPARENT TESTA 8 </italic>(<italic>TT8</italic>) gave a high level of activity in the absence of <italic>PAP1</italic>. Interestingly, TT8 is a bHLH transcription factor that is an important co-factor in the regulation of both anthocyanins and condensed tannins [##REF##11041882##39##]. The result here implies that, in our tobacco assay, the <italic>TT8 </italic>promoter fragment is relatively active transcriptionally.</p> <p>Eight PAP1 transactivated promoters and two non-responsive promoters were re-assayed and the transactivation confirmed (Fig. ##FIG##1##2B##). Of the eight trans-activated promoters identified by our analysis, six were from genes whose encoded proteins have a role in the anthocyanin and proanthocyanin biosynthetic pathways. In addition to these, a gene corresponding to a lipid transfer protein precursor (At5g59310) and a MYB transcription factor (At1g66380) were also identified in our microarray experiment and confirmed in transactivation assays. Lipid transfer proteins are a class of small basic soluble proteins capable of binding fatty acids and acyl CoA esters [##REF##15012303##40##]. As malonyl-CoA is an early precursor of the anthocyanin biosynthetic pathway [##REF##15608338##41##], it is possible that lipid transfer proteins may be acting in a transport role or as co-factors for the conversion of these intermediates to anthocyanins. The transient assay data also indicate that PAP1 was able to transactivate the transcription of a second MYB like gene <italic>MYB114 </italic>(At1g66380). In Arabidopsis <italic>MYB114 </italic>belongs to a tandem repeat with two other MYB genes, <italic>MYB113 </italic>and <italic>MYB90 </italic>(<italic>PAP2</italic>). All these genes show significant sequence similarity to <italic>PAP1 </italic>(<italic>MYB75</italic>), although only the <italic>PAP1 </italic>and <italic>PAP2 </italic>genes have been reported to regulate anthocyanin biosynthesis. This observation supports a potential role for the <italic>PAP1 </italic>gene in a feed-forward regulation of at least one related MYB gene.</p> <p>One promoter, for the gene encoding dihydroflavanol reductase (<italic>DFR</italic>; At5g42800), showed a 122-fold elevation in relative LUC activity in the transactivation assay. Three promoters from the genes encoding glutathione S-transferase (<italic>GST</italic>; At5g17220), leucoanthocyanidin dioxygenase (<italic>LDOX</italic>; At4g22880) and UDP flavonoid 5-O-glycosyltransferase (<italic>UFGT</italic>; At4g14090) showed a 38- to 60-fold elevation. The remaining four promoters derived from the genes encoding chalcone synthase (<italic>CHS</italic>; At5g13930), flavanone 3-hydroxylase (<italic>F3H</italic>; At3g51240), <italic>MYB114 </italic>(At1g66380) and a non-specific lipid transfer protein precursor (<italic>LTP</italic>; At5g59310), showed a PAP1-dependent elevation of 3- to 7-fold (Table ##TAB##0##1##).</p> <title>Motif analysis of transactivated promoter sequence</title> <p>The naive motif search programme MEME [##REF##7584402##42##], was used to search the DNA sequences for sequence motifs that were common to all eight promoters transactivated by the <italic>PAP1 </italic>gene. With the default settings and a maximum of five output motifs, only one 10 bp motif was present in all eight of the transactivated promoters (Table ##TAB##0##1##). This conserved motif was found in both plus and minus orientations and upstream of the 5' UTR where annotation is available. This motif was absent (<italic>P</italic>-value &lt; 0.119) from the 25 promoter fragments that were initially screened and did not alter relative LUC activity in the presence of the <italic>PAP1 </italic>gene. From these predictions, we hypothesise that this conserved motif (C/T)CNCCAC(A/G)(A/T)(G/T) is a PAP1 <italic>cis</italic>-regulatory element (PCE). Searches performed on the same promoter set using the search program COSMO [##REF##17402923##43##] yielded a related motif with a common core and related flanking sequences (C/T)(A/C)NCCACN(G/T)(G/T). When MEME analysis was conducted on the top 10 up-regulated genes from both this study and those previously published, neither the PCE nor any other motif was identified. This demonstrates the benefit of using only direct targets identified in the transient assay.</p> <title>Fold change and PCE frequency</title> <p>The level of transactivation from the co-infiltration of 35S-<italic>PAP1 </italic>varied between promoter-LUC cassettes (Table ##TAB##0##1##). Notably the promoter with the highest transactivation values At5g42800 (<italic>DFR</italic>), which showed a 122-fold increase in luciferase activation when co-infiltrated with 35S-<italic>PAP1</italic>, contains two PCEs within the promoter region used in the assay. Three promoters showed between 38- to 60-fold increases in relative LUC activity (<italic>GST</italic>; At5g17220, <italic>LDOX</italic>; At4g22880 and <italic>UFGT</italic>; At4g14090) and contained only a single PCE. Four other promoters (<italic>CHS</italic>; At5g13930, <italic>LTP</italic>; At5g59310, <italic>MYB114</italic>; At1g66380 and <italic>F3H</italic>; At5g51240) had much smaller 3- to 7-fold increases in relative LUC activity. This lower activation may be explained by the C to T change in the first base of the highly conserved CCAC core of the PCE motif in three of the promoters with the lowest transactivation values. However, this does not explain the data from the At3g51240 promoter, which had a low transactivation value and a fully conserved PCE motif (<italic>P</italic>-value = 5.56e<sup>-08</sup>). As the levels of expression from the transiently infiltrated 35S-<italic>PAP</italic>1 gene were higher than under normal physiological conditions, an alternative explanation for these lower transactivation values is that high levels of <italic>PAP1 </italic>expression may result in non-specific binding to low affinity sites in the promoter. However, as the 35S-<italic>PAP1 </italic>and promoter-LUC fusion were infiltrated in a ratio of 9:1, these transient leaf infiltration assays may more closely resemble the physiological ratio of TF to promoter than the over-expression of a TF in transgenic plants.</p> <title>Validation of PCE by transient leaf infiltration</title> <p>We tested the integrity of the PCE by deleting and mutating the sequence from the At5g17220 (<italic>GST</italic>) and At4g14090 (<italic>UFGT</italic>) promoter fragments. PCE promoter deletions had the 10 bp motif excised from the At5g17220 and At4g14090 promoter fragments. These were labelled At5g17220D and At4g14090D respectively. The sequence from At4g42800 (<italic>DFR</italic>) was not chosen for deletion or mutation analysis due to the presence of two PCEs which complicates base modification by PCR. Promoter mutations were generated by altering four bases of the PCE at positions 1 (C to T), 4 (C to T), 7 (C to T) and 10 (G to A) in the At5g17220 promoter. In promoter At4g14090, six base changes were made at positions 1 (C to T), 4 (C to T), 5 (C to A), 7 (C to T) 8 (A to T) and 10 (G to T). These modified promoters were labelled At5g17220M and At4g14090M respectively. Transactivation assays were used to compare these modified promoters with the unmodified At5g17220 and At4g14090 promoters (Fig. ##FIG##2##3##). At4g17220D had a significantly reduced LUC/REN ratio compared with the wild type sequence in the presence of 35S-<italic>PAP1</italic>. At4g17220M also showed a reduction in LUC/REN ratio in the presence of 35S-<italic>PAP1</italic>, compared with the wild type sequence. Modification of the PCE from At4g14090 showed similar reductions in relative LUC/REN ratio for At4g14090D and At4g140909M compared with the wild type. The six base changes in the PCE of At5g17220M appeared to have more of an effect than the 10 base changes introduced into the PCE region of the At4g14090 promoter. This may be due to different base substitutions at position 10, (At5g17220; G to A change, At4g14090; G to T change). A more detailed analysis of each of the conserved nucleotides within the PCE will be needed to fully interpret the significance of these results.</p> <p>Interestingly, five of the PCE-containing promoters identified in this study also contained a perfect G-box site (CACGTG) adjacent to the PCE sequence. A number of plant promoters regulated by diverse signals contain G-box elements [##REF##8127687##44##]. At least 2 classes of TFs are capable of binding G-boxes: the basic leucine zipper class (bZIP) and the basic helix-loop-helix proteins (bHLH) [##REF##11722549##19##]. Extensive genetic and protein studies have shown a close functional relationship between MYBs regulating anthocyanin accumulation and bHLH proteins [##REF##15361138##15##]. In maize, the activation of anthocyanin biosynthetic genes by ZmPl and ZmC1 requires a bHLH protein encoded by a <italic>R/B </italic>gene [##REF##1577278##45##]. In Arabidopsis the bHLH encoding <italic>TT8 </italic>gene and a MYB TF encoded by the <italic>TT2 </italic>gene, act synergistically to direct the expression of the <italic>DFR </italic>and <italic>BANYULS </italic>(<italic>BAN</italic>) flavonoid pathway genes [##REF##11041882##39##]. PAP1 has been demonstrated to interact with the bHLH proteins encoded by <italic>ENHANCER OF GLABRA3 </italic>(<italic>EGL3</italic>) and <italic>GLABRA3 </italic>(<italic>GL3</italic>) genes, and when co-overexpressed these combinations showed far more severe phenotypes than would be expected for additive regulation alone [##REF##12917293##46##]. The presence of the PCE and the G-box may not be coincidental and may not correspond to the binding site of PAP1, but that of the bHLH gene that we assume to be necessary for transactivation. In these transient assays we presume the appropriate endogenous bHLH protein interacts with the transient expressed <italic>PAP1 </italic>gene product.</p> <title>Occurrence of the PCE in 300 genes from microarray results</title> <p>We calculated probability values for putative PCEs occurring in the 2 kb upstream of the ATG of the top 300 up or down regulated genes selected from the combined mature plant and seedling microarrays, using MAST [##REF##9520501##47##]. This combined analysis did not select one of the eight (At4g14090) that was transactivated in the transient assay due to it having low difference of expression in seedlings. MAST is a search program which scans input sequences for known motifs and calculates match scores for each input sequence. From the MAST output a combined <italic>P</italic>-value can be obtained which measures the strength of the match of the sequence to the input motif. These combined <italic>P</italic>-values were plotted against the log-fold change calculated from the microarray results (Fig. ##FIG##3##4##). It was found that of the 6 genes that showed the biggest up-regulation of expression in the 35S-<italic>PAP1 </italic>plants, all had good (&lt;0.119) <italic>P</italic>-values for the PCE and were transactivated in the transient assay. When a less than log 2 fold up-regulation was observed, the occurrence of a PCE and transactivation became less easy to predict. Of the 7 tested genes that showed transactivation, the combined MAST <italic>P</italic>-value ranged from 2.8e<sup>-3 </sup>to 0.119 suggesting that in this range, transactivation could occur. One of the 25 tested promoters (At5g24770; <italic>VEGETATIVE STORAGE PROTEIN 2</italic>) that was not transactivated by the <italic>PAP1 </italic>gene did have a putative PCE (CCATCACAAG), but only with a <italic>P</italic>-value of 0.199. This suggests that this level of identity to the PCE consensus was insufficient to activate the gene. One additional promoter that was not transactivated (At5g05270, <italic>CHALCONE-FLAVANONE ISOMERASE</italic>) contained a PCE within this <italic>P-</italic>value range but located outside the 1 kb region tested for transactivation. Of the top 300 genes that showed changes in expression in the presence of the <italic>PAP1 </italic>transgene, a further 18 untested genes were significantly up-regulated and had a PCE within the transactivation <italic>P</italic>-value range, suggesting that these may also be regulated by PAP1 (Fig. ##FIG##3##4##). In addition to these 18 genes, there were 13 genes (4%) that had a <italic>P-</italic>value within this range, but were down-regulated. This either implies a repressor function or that the presence of the PCE alone is not sufficient to cause activation of these genes.</p>
<title>Conclusion</title> <p>The power of bioinformatics and the availability of whole genome sequence has enabled a comprehensive description of transcription factor families in plants. There is much less known about the mechanism that these genes employ to effect co-ordinated regulation. While methods that assay the direct interaction between DNA and proteins have proved effective in the characterisation of some of these genes, this is often limited to those proteins that can be easily purified, form simple complexes or have very high affinity for the target DNA. In addition to the core binding sites that seem to be associated with transcription factor families, there may a degree of subtlety in the <italic>cis</italic>-regulatory elements necessary for transcription factors to facilitate their unique regulatory effects. These gene-specific <italic>cis</italic>-regulatory regions may function through the recruitment of TF combinations or through a DNA motif consensus that is difficult to determine using conventional methods. Here we have used transient infiltration assays to analyse several promoters from unrelated genes that have a co-ordinated up-regulation in response to the over-expression of the MYB transcription factor PAP1. Using computer-based motif searches, we were able to identify a conserved region common to all promoters that were transactivated by the <italic>PAP1 </italic>gene product. While it is not necessarily the PAP1 binding domain, it is a region that is necessary for PAP1 regulation and as such, this method provides an effective tool to complement DNA-protein interaction assays in the effort to elucidate <italic>cis</italic>-regulatory domains of transcription factors. It is also worth noting that this assay uses a heterologous system based on the expression of Arabidopsis genes in tobacco, it is therefore possible that this expression pattern may differ from the native responses in Arabidopsis.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>Transcription factors (TFs) co-ordinately regulate target genes that are dispersed throughout the genome. This co-ordinate regulation is achieved, in part, through the interaction of transcription factors with conserved <italic>cis</italic>-regulatory motifs that are in close proximity to the target genes. While much is known about the families of transcription factors that regulate gene expression in plants, there are few well characterised <italic>cis</italic>-regulatory motifs.</p> <p>In Arabidopsis, over-expression of the MYB transcription factor <italic>PAP1 </italic>(<italic>PRODUCTION OF ANTHOCYANIN PIGMENT 1</italic>) leads to transgenic plants with elevated anthocyanin levels due to the co-ordinated up-regulation of genes in the anthocyanin biosynthetic pathway. In addition to the anthocyanin biosynthetic genes, there are a number of un-associated genes that also change in expression level. This may be a direct or indirect consequence of the over-expression of PAP1.</p> <title>Results</title> <p>Oligo array analysis of PAP1 over-expression Arabidopsis plants identified genes co-ordinately up-regulated in response to the elevated expression of this transcription factor. Transient assays on the promoter regions of 33 of these up-regulated genes identified eight promoter fragments that were transactivated by PAP1. Bioinformatic analysis on these promoters revealed a common <italic>cis</italic>-regulatory motif that we showed is required for PAP1 dependent transactivation.</p> <title>Conclusion</title> <p>Co-ordinated gene regulation by individual transcription factors is a complex collection of both direct and indirect effects. Transient transactivation assays provide a rapid method to identify direct target genes from indirect target genes. Bioinformatic analysis of the promoters of these direct target genes is able to locate motifs that are common to this sub-set of promoters, which is impossible to identify with the larger set of direct and indirect target genes. While this type of analysis does not prove a direct interaction between protein and DNA, it does provide a tool to characterise <italic>cis</italic>-regulatory sequences that are necessary for transcription activation in a complex list of co-ordinately regulated genes.</p>
<title>Competing interests</title> <p>The authors declare that they have no competing interests.</p> <title>Authors' contributions</title> <p>APD for experimental design, promoter cloning, data collection, sequence analysis and manuscript preparation. RJS carried out microarray analysis and contributed to manuscript preparation. KLW for expression analysis of transgenic plants. ACA transgenic plant growth and generation. RPH contributed to experimental design and manuscript preparation.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>We thank William Laing, Andrew Gleave and Cathie Martin for comments on this manuscript and Ariel Liu for technical assistance.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p>Analysis of transgenic Arabidopsis over-expressing <italic>PAP1</italic>. A, Plants at the 21 day stage containing vector only (top) and 35S-<italic>PAP1 </italic>(bottom) B, Comparison of intensity-ratio plot of gene expression in transgenic plants; Blue dot represents genes that change in expression with an adjusted <italic>P </italic>value &lt;0.05, red dot represents 33 genes cloned for promoter analysis and green dot represents genes with significant transactivation when co-infiltrated with the 35S-<italic>PAP1 </italic>construct. The grey dot represents <italic>PAP1 </italic>that was over-expressed. C, Comparison of the 33 genes selected for promoter analysis. Ratios generated between wild type and <italic>PAP1 </italic>over-expression plants are shown. Ratios from seedlings (Seedling), mature plants (Mature) and a combined dataset of mature plants and seedlings (Combined) are given. Ratios with * are not selected.</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p>Transient assay data for promoters up-regulated in <italic>PAP1 </italic>over-expression transgenic plants. A, High-throughput screen of 33 promoters assayed with (black bars) or without (white bars) co-infiltration of <italic>PAP1</italic>. B, Re-analysis and transactivation standard error estimation for eight promoters identified in primary screen and two non-responsive control promoters (At2g37040 and At1g56650). Average transactivation values were calculated from 6 fold experimental replicates. Standard errors were calculated for error bars.</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p>Deletion and mutation of PAP1 <italic>cis</italic>-regulatory elements from At5g17220 and At4g14090 promoters. Transactivation data of At5g17220 and At4g14090 promoters displaying the effect of deleting and mutating the 10 bp PCE from At5g17220 and At4g14090 promoters. Wild-type (WT) promoter, PCE deletion promoter (Del) and PCE mutation promoter (Mut) are shown and transactivation data as black bars (with PAP1), and white bars (without PAP1).</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p>Occurrence of the PAP1 <italic>cis</italic>-regulatory element in the top 300 microarray selected genes. A 3 kb promoter fragment for each of the top 300 genes was analysed for the presence of PCEs using MAST. The average fold change from all microarrays was compared with the <italic>P</italic>-value obtained from the MAST analysis. Genes highlighted in green represent those genes that have a PCE and were transactivated by PAP1. Blue represents those that have a putative PCE but not in the first 1 kb tested. Yellow represents a gene with a PCE in the first 1 kb but was not transactivated. Red represents genes with no PCE and not transactivated. Grey represents untested genes.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>The PAP1 <italic>cis</italic>-regulatory element in eight transactivated promoters</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left"><bold>EST Identifier</bold></td><td align="center"><bold>Gene </bold><break/><bold>description</bold></td><td align="center"><bold>PCE sequence</bold></td><td align="center"><bold>Fold increase </bold><break/><bold>LUC:REN</bold></td></tr></thead><tbody><tr><td align="left">At5g42800</td><td align="center"><italic>DFR</italic></td><td align="center">CCA<underline>CCAC</underline>GTG</td><td align="center">121.6</td></tr><tr><td/><td/><td align="center">TCC<underline>CCAC</underline>GTG</td><td/></tr><tr><td align="left">At5g17220</td><td align="center"><italic>GST</italic></td><td align="center">CCA<underline>CCAC</underline>ATG</td><td align="center">60</td></tr><tr><td align="left">At4g22880</td><td align="center"><italic>LDOX</italic></td><td align="center">TCT<underline>CCAC</underline>GTG</td><td align="center">37.8</td></tr><tr><td align="left">At4g14090</td><td align="center"><italic>UFGT</italic></td><td align="center">CCT<underline>CCAC</underline>AAG</td><td align="center">40.5</td></tr><tr><td align="left">At5g13930</td><td align="center"><italic>CHS</italic></td><td align="center">TCA<underline>TCAC</underline>ATT</td><td align="center">4.1</td></tr><tr><td align="left">At5g59310</td><td align="center"><italic>LTPP</italic></td><td align="center">CCA<underline>TCAC</underline>GTT</td><td align="center">7.4</td></tr><tr><td align="left">At1g66380</td><td align="center"><italic>MYB114</italic></td><td align="center">CCG<underline>TCAC</underline>GTG</td><td align="center">3.1</td></tr><tr><td align="left">At3g51240</td><td align="center"><italic>F3H</italic></td><td align="center">CCG<underline>CCAC</underline>GTG</td><td align="center">6.3</td></tr><tr><td align="center" colspan="2">PCE consensus sequence</td><td align="center">YCN<underline>CCAC</underline>RWK</td><td/></tr></tbody></table></table-wrap>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p>1744 genes up-regulated on microarray of 35S-<italic>PAP1 </italic>Arabidopsis.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S2"><caption><title>Additional file 2</title><p>Oligonucleotide primer pairs used to PCR amplify 1 kb of Arabidopsis promoter sequence.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S3"><caption><title>Additional file 3</title><p>Oligonucleotide primer pairs used to delete or mutate the PCE in the At5g17220 and At4g14090 promoters.</p></caption></supplementary-material>
<table-wrap-foot><p>Conserved CCAC core is underlined. Fold increase in LUC/REN ratio over promoter only control is also shown.</p></table-wrap-foot>
<graphic xlink:href="1746-4811-4-17-1"/> <graphic xlink:href="1746-4811-4-17-2"/> <graphic xlink:href="1746-4811-4-17-3"/> <graphic xlink:href="1746-4811-4-17-4"/>
<media xlink:href="1746-4811-4-17-S1.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="1746-4811-4-17-S2.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="1746-4811-4-17-S3.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Eckardt"], "given-names": ["NA"], "article-title": ["Light regulation of plant development: HY5 Genomic binding sites"], "source": ["Plant Cell"], "year": ["2007"], "volume": ["19"], "fpage": ["727"], "lpage": ["729"]}, {"surname": ["Chang", "Puryear", "Cairney"], "given-names": ["S", "J", "J"], "article-title": ["A simple and efficient method for isolating RNA from pine trees."], "source": ["Plant Mol Biol Rep"], "year": ["1993"], "volume": ["11"], "fpage": ["114"], "lpage": ["117"]}, {"collab": ["Smyth GK", "Speed T"], "article-title": ["Normalization of cDNA microarray data"], "source": ["Methods (Orlando)"], "year": ["2003"], "volume": ["31"], "fpage": ["265"], "lpage": ["273"]}, {"surname": ["Benjamini", "Hochberg"], "given-names": ["Y", "Y"], "article-title": ["Controlling the false discovery rate: a practical and powerful approach to multiple testing"], "source": ["J R Stat Soc [Ser A]"], "year": ["1995"], "volume": ["57"], "fpage": ["289"], "lpage": ["300"]}]
Plant Methods. 2008 Jul 7; 4:17