Skip to main content

Development of a novel data mining tool to find cis-elements in rice gene promoter regions

Abstract

Background

Information on more than 35 000 full-length Oryza sativa cDNAs, together with associated microarray gene expression data collected under various treatment conditions, has made it feasible to identify motifs that are conserved in gene promoters and may act as cis-regulatory elements with key roles under the various conditions.

Results

We have developed a novel tool that searches for cis-element candidates in the upstream, downstream, or coding regions of differentially regulated genes. The tool first lists cis-element candidates by motif searching based on the supposition that if there are cis-elements playing important roles in the regulation of a given set of genes, they will be statistically overrepresented and will be conserved. Then it evaluates the likelihood scores of the listed candidate motifs by association rule analysis. This strategy depends on the idea that motifs overrepresented in the promoter region could play specific roles in the regulation of expression of these genes. The tool is designed so that any biological researchers can use it easily at the publicly accessible Internet site http://hpc.irri.cgiar.org/tool/nias/ces. We evaluated the accuracy and utility of the tool by using a dataset of auxin-inducible genes that have well-studied cis-elements. The test showed the effectiveness of the tool in identifying significant relationships between cis-element candidates and related sets of genes.

Conclusion

The tool lists possible cis-element motifs corresponding to genes of interest, and it will contribute to the deeper understanding of gene regulatory mechanisms in plants.

Background

With the completion of rice genome sequencing by the International Rice Genome Sequencing Project [1], the Beijing Genomics Institute (BGI) [2], and Syngenta [3], many rice functional genomic resources have become available, including whole genome sequences from ssp. japonica 'Nipponbare' and ssp. indica line 93-11; a set of rice full-length cDNA clones and their complete and partial end sequences [4, 5], microarray gene expression systems based on full-length cDNA sequences, ESTs (Expressed Sequence Tag), MPSS (Massively Parallel Signature Sequencing), SAGE (Serial Analysis of Gene Expression), and predicted genes in the genome sequences; and many kinds of insertion mutants with Tos17, Ac-Ds, and T-DNAs [6]. As analytical technology progresses, the database continues to be upgraded and serves as a useful resource for studying mechanisms that regulate gene expression.

Cis-elements in the promoter regions of genes and trans-acting transcription factors are major biological features to be characterized if we are to achieve an understanding of the systems that regulate gene expression. Identification of candidate cis-elements corresponding to genes is now practicable through the use of available sequence and genome mapping information, combined with information about the responses of genes to specific experimental conditions; such responses have been elucidated by using gene expression profiles now publicly available.

Exhaustive sequence analysis by using available public databases can identify cis-element candidate motifs for further examination, but such approaches are not quite efficient. One confounding factor is that public databases are independently constructed and not generally optimized to facilitate integration of information from many sources with local experimental data. A more perplexing issue for experimental researchers who are not very familiar with bioinformatics techniques is the challenge of finding unknown but biologically notable relationships among genes, cis-elements, and experimental conditions from the huge number of possible combinations generated by large experimental data sets.

To resolve some of these issues, we developed a novel data mining tool to identify cis-elements in the rice genome. It performs the complex bioinformatics analysis mentioned above, then lists cis-element candidates for genes. The genes can be grouped by similarity of expression profiles and other criteria for assessment by researchers, then the tool annotates them with related public database information.

Similar tools have been developed previously. Helden released RSAT, which includes a program that can detect over-represented motifs in upstream regions of co-regulated genes [7]. Holt et al. established CoReg, which links the hierarchical clustering of co-expressed gene sets with frequency tables of promoter elements [8]. Zhao et al. established TRED, which integrates a database and a system for predicting cis- and trans-elements in mammals [9]. Galuschka et al. developed AthaMAP, which includes a program for comparative analysis of cis-elements in sets of co-transcribed genes of Arabidopsis thaliana [10].

Our tool is distinguished by several points: (i) It focuses on the rice genome, being based on full-length cDNAs, and is designed to pick up cis-element candidates associated with genes that users designate. (ii) It evaluates the likelihood score of cis-element candidates by comparing frequency counts in the user-selected gene set and a reference gene set. (iii) It can evaluate previously known cis-element sequences as well as user-specified sequences prepared by other analysis tools, and it can examine several cis-elements together.

The tool carries out both ab initio motif searches of promoter sequences and searches against known plant cis-elements, then performs a likelihood analysis of identified cis-elements on the basis of their presence in a significant proportion of the promoters of a given set of genes. This evaluation is achieved by an association rule analysis.

Here, we present technical details of the tool and demonstrate the practical assessment of its utility with a biologically relevant sample data set.

Implementation

The tool, called Rice Cis-Element Searcher (RiCES), consists of a cis-element searching pipeline, controlled via a Web-based user interface. Fig. 1 summarizes the procedure. The pipeline first reads a list of gene identifiers from the user, which it uses to retrieve the promoter sequences corresponding to the listed genes. Then a preliminary list of cis-element candidates is built by aligning information from the built-in list of plausible motifs, or by ab initio motif searching of the sequence data. Association rule analysis is carried out and reported to support the candidacy of the resulting cis-element list.

Figure 1
figure 1

Features of RiCES.

Gene list

RiCES assumes that a user has already identified genes of interest from experimental analysis (e.g. clusters of coordinately regulated genes). The list of identifiers is input into a Web-based data entry form. RiCES recognizes GenBank accession numbers, identifiers of transcription units (TUs) as defined in the TIGR pseudomolecular assemblies [11], and several other major gene identification systems. Using the list, it retrieves the set of associated upstream, downstream, or coding region sequences flanking the specified genes from available genomic sequence data.

Preliminary cis-element candidate list

The second step of the analysis is the compilation of a list of motifs as candidate cis-elements. At present RiCES supports two methods to achieve this.

The first method depends on ab initio motif searching based on the supposition that if there are cis-elements playing important roles in the regulation of a given set of genes, they will be statistically overrepresented in the associated promoter sequences as conserved motifs that can be identified by using a suitable motif search program. There are several programs implementing several algorithms. We have chosen to use MEME, which is a publicly available motif discovery program [12] supporting an expectation maximization algorithm. In our analysis algorithm, MEME is invoked to identify motifs 6 to 8 bp long that look highly conserved among promoter sequences of the selected genes. Users can modify some of the search parameters of the MEME program via the Web form.

The second method relies on the hypothesis that common, known cis-elements play important roles under the experimental conditions that gave rise to the list of genes specified by the user. Therefore, RiCES searches for matches to a pre-compiled list of known cis-elements.

Several databases of plant cis-elements are publicly available. PLACE [13] is one of the most popular databases of known cis-elements in plant genomes. AtcisDB, a part of AGRIS [14], includes information on cis-elements involved in gene regulation in Arabidopsis thaliana.

Although these databases are extremely useful resources, it is not straightforward to cross-link information from them directly to the researcher's own data. Current databases are not exhaustive enough to distinguish 'core' motifs, which decide the function of cis-elements, from co-existing sequences in neighboring regions. As a result, many cis-element sequence data in these databases include superficial core motifs for which no evidence of functionality has been obtained. The use of such data prohibits effective informatic analysis.

We compiled a novel database of known cis-elements and incorporated it into RiCES [See Additional file 1]. The cis-elements are collected from reports of experiments such as gel shift assays and footprint analyses, categorized by transcription factor, and documented with respect to known activity in the plant genome. Some cis-elements known only in organisms other than plants are also listed, in consideration of their possible, albeit unknown, roles in plants. The database includes four types of cis-elements: (1) G-box and E-box, which bind to common sequences such as bHLH or bZIP in many organisms; (2) A-box, T-box, and GGTTTAG repeats, which bind to common sequences in many organisms, such as homeodomain and Myb; (3) CArG boxes and GCC-box, which bind to plant MADS, zinc finger, and AP2/EREBP elements; and (4) other cis-elements, binding only in animals, such as HSF, PcG, and HMG.

Association rule analysis

The third step of the analysis is the likelihood evaluation of the cis-element candidates by association rule analysis, which is a data mining method designed to discover significant relationships between pairs of characteristics observed in data sets. Candidates showing the highest likelihood (specificity) are retained in the final cis-element candidate list.

Association rule analysis has been applied to mechanisms that regulate gene expression [e.g. [15, 16]]. We used it to find relationships between identified cis-elements and gene expression profiles. The strategy depends on the idea that motifs overrepresented in the promoter region of the genes of interest could play specific roles in regulation of the expression of those genes.

Implied cause-and-effect relationships documented as 'rules' are evaluated by using several well-known indices of likelihood, including support, confidence, and lift [15]. On the basis of sample data sets, the lift index appeared to best discriminate significant relationships between experimental conditions and cis-element candidates.

In a rule described as

the presence of motif X in a gene implies that the gene is a member of group Y,

lift is the ratio of the posterior probability (the probability that the gene is in group Y if it possess motif X) to the prior probability (the probability of X possession, irrespective of the membership of Y). When lift > 1.0, the coexistence of X and Y is not a random occurrence, but suggests some causal relationship between them. If lift < 1.0, it is not considered probabilistically significant. Consequently, we set the default threshold of lift to 1.0, and the cis-element candidates are included in the final candidate list only if their lift value is higher than this threshold.

RiCES also evaluates pairwise combinations of motifs in the preliminary candidate list (upper right-hand box in Fig. 1), in consideration of possible protein-protein interactions of multiple transcription elements binding cis-elements, as illustrated by experimental evidence [17, 18].

Output

The final cis-element candidate list is presented as an association table with the identifier of the submitted genes (TU identifiers based on TIGR gene model annotation are used in the current version) annotated with any available corresponding information from RiceCyc [19] and Gene Ontology [20]. RiCES also provides information on candidate motifs, including the positions of the element in the promoter regions of corresponding TUs, the sequence, and related information from AtcisDB [14]. The position of the cis-element candidates is also presented in both text and graphics.

Validation

To test whether or not the output of RiCES was meaningful, we validated it with a list of auxin-inducible genes with known characteristics, compiled from RiceTFDB 2.0 [21]. First, Aux/IAA genes stored in RiceTFDB were applied as queries in a BLASTN search [22] of GenBank, returning a list containing 28 rice TUs [See Additional file 2]. These genes were fed into the pipeline. When the MEME program was called, the length of target motifs was set to 6, 7, or 8 bases, the number of occurrences of each motif was set to 7, 14, or 21, and the search algorithm was set to 'zoops' to check zero or one occurrence per sequence. The outputs of each option setting were merged but not otherwise filtered.

Results and Discussion

Many Aux/IAA genes are auxin-inducible [23] and contain the TGTCTC element [24]. This element is commonly found in the upstream region of auxin-responsive genes. Thus, the detection of all instances of the motif by the pipeline could serve as a validation of the pipeline algorithm. The auxin-responsive element (AuxRE) containing the TGTCTC motif in some cases requires another proximal AuxRE for biological activity [17, 25]. In other contexts, AuxRE functions only when it occurs with its palindromic components separated by 7 or 8 nucleotides [26].

In our validation test, MEME listed 7514 motifs in total from 1000 bp of the upstream sequences [See Additional file 3], of which 4128 showed a high lift value (>1.0) [See Additional file 4]. A search of AtcisDB for these motifs returned 4 showing a partial match to the record of 'PRHA binding sites' (Table 1), which is derived from the report of Plesch et al. [27], describing auxin-induced expression of the Arabidopsis prha homeobox gene. Another 4 motifs contained the TGTCTC element. The result was consistent with previous work, as TGTCTC was listed as a candidate in the single motif search of Aux/IAA genes.

Table 1 Cis-element candidate motifs from Aux/IAA genes and suggested to be auxin-induction related according to ATCIS.

Table 2 shows the result of the validation test with a pre-compiled cis-element list generated by the test gene list. The analysis returned 22 cis-element candidates with lift > 1.0 [See Additional file 5 and 6]. Some of these candidates were suggested by previous studies to have some kind of relationship to auxin response. For example, RAV1 was found in the promoter region of ABP, which encodes an auxin-binding protein [28]. Expression of LEAFY (LFY) is affected by the auxin gradient in Arabidopsis [29]. ETT is another auxin response factor [30], and LFY and ETT expression are closely correlated [18, 31].

Table 2 Cis-element candidates selected from the pre-compiled list, likely corresponding to Aux/IAA genes.

The position of a cis-element is important information to consider in relation to the function of the cis-element. For biological activity to occur, the distance of some cis-elements from the coding region or other collaborating elements is constrained. To this end, RiCES highlights the distribution of cis-element candidates. It provides tables of identified cis-element motifs and graphical motif maps to help researchers grasp positional relationships among the candidate elements.

The positions of the listed elements, some of which include TGTCTC, varied among upstream regions of genes (Fig. 2), and it was hard to detect any skewed distribution of motifs. Goda et al. [32] studied the distribution of TGTCTC motifs in the genome of A. thaliana, and pointed out that 25% of investigated genes had TGTCTC motifs in the upstream region within 1000 bp of the start codon, and 14% within 500 bps. Our results do not seem in conflict of theirs.

Figure 2
figure 2

Distribution of the 15 Aux/IAA-related cis -element candidates. The presence of the motifs of candidates with high lift values (see 4th column in Table 1) was searched in the 1000-bp upstream region of genes, and frequency was counted in segmented regions at an interval of 10 bp. The X-axis represents the position in the upstream region, and the bars designate frequency of motifs (counted after distribution of multiple regions was merged).

TGTCTC motifs are scattered over wide regions of many plant species (Table 3). It is possible that the variety of the roles of genes reflects the variety of mechanisms regulating gene expression and positions of cis-elements, even if the genes in question can be classified as 'auxin-responsive genes' in a larger sense.

Table 3 Representative plant genes possessing TGTCTC element in corresponding upstream region.

A major research concern is how to pick up cis-element candidates worthy of further experimentation. Computational and manual selection of cis-element candidates should play complementary roles to resolve this issue. It should be emphasized that cis-element candidates listed by RiCES are rated according to the likelihood provided by association rule analysis. On the other hand, researchers can check the significance of candidates in detail by using related information derived from several databases. The supported databases include AGRIS, Gene Ontology, and RiceCyc, as well as the map information described above.

Fig. 3 is an example of the output for the TGTCTC motif. The outputs are not only easily accessible in a Web browser, but are also usable in further statistical or bioinformatics analysis, as they are also provided in XML format (Fig. 3A), which is a tagged plain-text format compatible with various computer programs.

Figure 3
figure 3

Snapshots of representative outputs of RiCES. A: List of cis-element candidate motifs including related information. B: Mapping image of cis-element candidate motifs.

In some cases, the results of the analysis from the pre-compiled list of elements will be easily comparable with prior knowledge. In other cases involving solely ab initio evidence from MEME, the results of motif searches should be interpreted carefully, because the result will change considerably in accordance with the options selected. An appropriate set of motif search options should be determined each time, by trial and error. However, as described above, a motif search can find cis-element candidates of which the sequences do not exactly match those of known cis-elements.

Although RiCES is focused on the role of cis-elements in Oryza sativa ssp. japonica, the methodology can be applied easily to studies of other plant species, or of other genome sequence motifs involving gene expression regulation, such as motifs in coding regions of genes or downstream of the gene sequence. Such work can be made possible by replacing the reference data set containing whole genes of rice with other data sets.

Conclusion

We presented here a newly developed tool to search for cis-element candidates in a list of genes. A case study showed the applicability of the tool. The tool is easy to use and publicly available. We expect that its use will deepen understanding of the mechanisms that regulate gene expression in plants.

Availability and requirements

RiCES is accessible at http://hpc.irri.cgiar.org/tool/nias/ces by any JavaScript-capable browsers.

Project Name: Generation Challenge Programme Subprogramme 4

Project Home Page: http://www.generationcp.org/subprogramme4.php

Operating system(s): Platform independent

Other requirements: None

Programming language: Perl

License: Freely available for use

Any restrictions to use by non-academics: None

References

  1. IRGSP Project: The map-based sequence of the rice genome. Nature. 2005, 436: 793-800. 10.1038/nature03895.

    Article  Google Scholar 

  2. Yu J, Hu S, Wang J, Wong GK, Li S, Liu B, Deng Y, Dai L, Zhou Y, Zhang X, Cao M, Liu J, Sun J, Tang J, Chen Y, Huang X, Lin W, Ye C, Tong W, Cong L, Geng J, Han Y, Li L, Li W, Hu G, Huang X, Li W, Li J, Liu Z, Li L, Liu J, Qi Q, Liu J, Li L, Li T, Wang X, Lu H, Wu T, Zhu M, Ni P, Han H, Dong W, Ren X, Feng X, Cui P, Li X, Wang H, Xu X, Zhai W, Xu Z, Zhang J, He S, Zhang J, Xu J, Zhang K, Zheng X, Dong J, Zeng W, Tao L, Ye J, Tan J, Ren X, Chen X, He J, Liu D, Tian W, Tian C, Xia H, Bao Q, Li G, Gao H, Cao T, Wang J, Zhao W, Li P, Chen W, Wang X, Zhang Y, Hu J, Wang J, Liu S, Yang J, Zhang G, Xiong Y, Li Z, Mao L, Zhou C, Zhu Z, Chen R, Hao B, Zheng W, Chen S, Guo W, Li G, Liu S, Tao M, Wang J, Zhu L, Yuan L, Yang H: A draft sequence of the rice genome (Oryza sativa L. ssp. indica). Science. 2002, 296: 79-92. 10.1126/science.1068037.

    Article  PubMed  Google Scholar 

  3. Goff SA, Ricke D, Lan T, Presting G, Wang R, Dunn M, Glazebrook J, Sessions A, Oeller P, Varma H, Hadley D, Hutchison D, Martin C, Katagiri F, Lange BM, Moughamer T, Xia Y, Budworth P, Zhong J, Miguel T, Paszkowski U, Zhang S, Colbert M, Sun W, Chen L, Cooper B, Park S, Wood TC, Mao L, Quail P, Wing R, Dean R, Yu Y, Zharkikh A, Shen R, Sahasrabudhe S, Thomas A, Cannings R, Gutin A, Pruss D, Reid J, Tavtigian S, Mitchell J, Eldredge G, Scholl T, Miller RM, Bhatnagar S, Adey N, Rubano T, Tusneem N, Robinson R, Feldhaus J, Macalma T, Oliphant A, Briggs S: A draft sequence of the rice genome (Oryza sativa L. ssp. japonica). Science. 2002, 296: 92-100. 10.1126/science.1068275.

    Article  PubMed  Google Scholar 

  4. Rice Full-Length cDNA Consortium: Collection, mapping, and annotation of over 28,000 cDNA clones from japonica rice. Science. 2003, 301: 376-379. 10.1126/science.1081288.

    Article  Google Scholar 

  5. Satoh K, Doi K, Nagata T, Kishimoto N, Suzuki K, Otomo Y, Kawai J, Nakamura M, Hirozane-Kishikawa T, Kanagawa S, Arakawa T, Takahashi-Iida J, Murata M, Ninomiya N, Sasaki D, Fukuda S, Tagami M, Yamagata H, Kurita K, Kamiya K, Yamamoto M, Kikuta A, Bito T, Fujitsuka N, Ito K, Kanamori H, Choi I, Nagamura Y, Matsumoto T, Murakami K, Matsubara K, Carninci P, Hayashizaki Y, Kikuchi S: Gene organization in rice revealed by full-length cDNA mapping and gene expression analysis through microarray. PLoS One. 2007, 2: e1235-10.1371/journal.pone.0001235.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Hirochika H, Guiderdoni E, An G, Hsing Y, Eun MY, Han C, Upadhyaya N, Ramachandran S, Zhang Q, Pereira A, Sundaresan V, Leung H: Rice mutant resources for gene discovery. Plant Mol Biol. 2004, 54: 325-334. 10.1023/B:PLAN.0000036368.74758.66.

    Article  PubMed  Google Scholar 

  7. van Helden J: Regulatory sequence analysis tools. Nucleic Acids Res. 2003, 31: 3593-3596. 10.1093/nar/gkg567.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Holt KE, Millar AH, Whelan J: ModuleFinder and CoReg: alternative tools for linking gene expression modules with promoter sequences motifs to uncover gene regulation mechanisms in plants. Plant Methods. 2006, 2: 8-10.1186/1746-4811-2-8.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Zhao F, Xuan Z, Liu L, Zhang MQ: TRED: a Transcriptional Regulatory Element Database and a platform for in silico gene regulation studies. Nucleic Acids Res. 2005, 33: D103-D107. 10.1093/nar/gki004.

    Article  PubMed  PubMed Central  Google Scholar 

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

    Article  PubMed  PubMed Central  Google Scholar 

  11. Ouyang S, Zhu W, Hamilton J, Lin H, Campbell M, Childs K, Thibaud-Nissen F, Malek RL, Lee Y, Zheng L, Orvis J, Haas B, Wortman J, Buell CR: The TIGR Rice Genome Annotation Resource: improvements and new features. Nucleic Acids Res. 2007, 35: D883-D887. 10.1093/nar/gkl976.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Bailey TL, Elkan C: Fitting a mixture model by expectation maximization to discover motifs in biopolymers. Proc Int Conf Intell Syst Mol Biol. 1994, 2: 28-36.

    PubMed  Google Scholar 

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

    Article  PubMed  PubMed Central  Google Scholar 

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

    Article  PubMed  PubMed Central  Google Scholar 

  15. Carmona-Saez P, Chagoyen M, Rodriguez A, Trelles O, Carazo JM, Pascual-Montano A: Integrated analysis of gene expression by association rules discovery. BMC Bioinformatics. 2006, 7: 54-10.1186/1471-2105-7-54.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Conklin D, Jonassen I, Aasland R, Taylor WR: Association of nucleotide patterns with gene function classes: application to human 3' untranslated sequences. Bioinformatics. 2002, 18: 182-189. 10.1093/bioinformatics/18.1.182.

    Article  PubMed  Google Scholar 

  17. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Ulmasov T, Hagen G, Guilfoyle TJ: Dimerization and DNA binding of auxin response factors. Plant J. 1999, 19: 309-319. 10.1046/j.1365-313X.1999.00538.x.

    Article  PubMed  Google Scholar 

  19. Gramene pathway tools (RiceCyc). [http://www.gramene.org/pathway/]

  20. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene Ontology: tool for the unification of biology. Nat Genet. 2000, 25: 25-29. 10.1038/75556.

    Article  PubMed  PubMed Central  Google Scholar 

  21. RiceTFDB. [http://ricetfdb.bio.uni-potsdam.de/]

  22. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic Local Alignment Search Tool. J Mol Biol. 1990, 215: 403-410.

    Article  PubMed  Google Scholar 

  23. Reed JW: Roles and activities of Aux/IAA proteins in Arabidopsis. Trends Plant Sci. 2001, 6: 420-425. 10.1016/S1360-1385(01)02042-8.

    Article  PubMed  Google Scholar 

  24. Tiwari SB, Wang XJ, Hagen G, Guilfoyle TJ: AUX/IAA proteins are active repressors, and their stability and activity are modulated by auxin. Plant Cell. 2001, 13: 2809-2822. 10.1105/tpc.13.12.2809.

    Article  PubMed  PubMed Central  Google Scholar 

  25. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  26. 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.

    Article  PubMed  Google Scholar 

  27. Plesch G, Stoermann K, Torres JT, Walden R, Somssich IE: Developmental and auxin-induced expression of the Arabidopsis prha homeobox gene. Plant J. 1997, 12: 635-647. 10.1046/j.1365-313X.1997.00635.x.

    Article  PubMed  Google Scholar 

  28. Kagaya Y, Ohmiya K, Hattori T: RAV1, a novel DNA-binding protein, binds to bipartite recognition sequence through two distinct DNA-binding domains uniquely found in higher plants. Nucleic Acids Res. 1999, 27: 470-478. 10.1093/nar/27.2.470.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Ezhova TA, Soldatova OP, Kalinina AIu, Medvedev SS: Interaction of ABRUPTUS/PINOID and LEAFY genes during floral morphogenesis in Arabidopsis thaliana (L.) Heynh. Genetika. 2000, 36: 1682-1687.

    PubMed  Google Scholar 

  30. Sessions A, Nemhauser JL, McColl A, Roe JL, Feldmann KA, Zambryski PC: ETTIN patterns the Arabidopsis floral meristem and reproductive organs. Development. 1997, 124: 4481-4491.

    PubMed  Google Scholar 

  31. Remington DL, Vision TJ, Guilfoyle TJ, Reed JW: Contrasting modes of diversification in the Aux/IAA and ARF gene families. Plant Physiol. 2004, 135: 1738-1752. 10.1104/pp.104.039669.

    Article  PubMed  PubMed Central  Google Scholar 

  32. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Nag R, Maity MK, Dasgupta M: Dual DNA binding property of ABA insensitive 3 like factors targeted to promoters responsive to ABA and auxin. Plant Mol Biol. 2005, 59: 821-838. 10.1007/s11103-005-1387-z.

    Article  PubMed  Google Scholar 

  34. Yang G, Nakamura H, Ichikawa H, Kitano H, Komatsu S: OsBLE3, a brassinolide-enhanced gene, is involved in the growth of rice. Phytochemistry. 2006, 67: 1442-1454. 10.1016/j.phytochem.2006.05.026.

    Article  PubMed  Google Scholar 

  35. Hsu CY, Jenkins J, Saha S, Ma DP: Transcriptional regulation of the lipid transfer protein gene LTP3 in cotton fibers by a novel MYB protein. Plant Sci. 2005, 168: 167-181. 10.1016/j.plantsci.2004.07.033.

    Article  Google Scholar 

  36. Bai F, Watson JC, Walling J, Weeden N, Santner AA, DeMason DA: Molecular characterization and expression of PsPK2, a PINOID-like gene from pea (Pisum sativum). Plant Sci. 2005, 168: 1281-1291. 10.1016/j.plantsci.2005.01.005.

    Article  Google Scholar 

  37. Szopa J, Lukaszewicz M, Aksamit A, Korobczak A, Kwiatkowska D: Structural organisation, expression, and promoter analysis of a 16R isoform of 14-3-3 protein gene from potato. Plant Physiol Biochem. 2003, 41: 417-423. 10.1016/S0981-9428(03)00048-2.

    Article  Google Scholar 

  38. Navarro-Avino JP, Bennett AB: Role of a Ca2+ -ATPase induced by ABA and IAA in the generation of specific Ca2+ signals. Biochem Biophys Res Commun. 2005, 329: 406-415. 10.1016/j.bbrc.2005.01.142.

    Article  PubMed  Google Scholar 

  39. Ishiki Y, Oda A, Yaegashi Y, Orihara Y, Arai T, Hirabayashi T, Nakagawa H, Sato T: Cloning of an auxin-responsive 1-aminocyclopropane-1-carboxylate synthase gene (CMe-ACS2) from melon and the expression of ACS genes in etiolated melon seedlings and melon fruits. Plant Sci. 2000, 159: 173-181. 10.1016/S0168-9452(00)00298-3.

    Article  PubMed  Google Scholar 

  40. Ge L, Chen H, Jiang JF, Zhao Y, Xu ML, Xu YY, Tan KH, Xu ZH, Chong K: Overexpression of OsRAA1 causes pleiotropic phenotypes in transgenic rice plants, including altered leaf, flower, and root development and root response to gravity. Plant Physiol. 2004, 135: 1502-1513. 10.1104/pp.104.041996.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Borisov AY, Madsen LH, Tsyganov VE, Umehara Y, Voroshilova VA, Batagov AO, Sandal N, Mortensen A, Schauser L, Ellis N, Tikhonovich IA, Stougaard J: The Sym35 gene required for root nodule development in pea is an ortholog of Nin from Lotus japonicus. Plant Physiol. 2003, 131: 1009-1017. 10.1104/pp.102.016071.

    Article  PubMed  PubMed Central  Google Scholar 

  42. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Carrasco JL, Ancillo G, Mayda E, Vera P: A novel transcription factor involved in plant defense endowed with protein phosphatase activity. EMBO J. 2003, 22: 3376-3384. 10.1093/emboj/cdg323.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Esmon CA, Tinsley AG, Ljung K, Sandberg G, Hearne LB, Liscum E: A gradient of auxin and auxin-dependent transcription precedes tropic growth responses. Proc Natl Acad Sci USA. 2006, 103: 236-241. 10.1073/pnas.0507127103.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Okumoto S, Schmidt R, Tegeder M, Fischer WN, Rentsch D, Frommer WB, Koch W: High affinity amino acid transporters specifically expressed in xylem parenchyma and developing seeds of Arabidopsis. J Biol Chem. 2002, 277: 45338-45346. 10.1074/jbc.M207730200.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

This work was supported by a grant from the Generation Challenge Programme SP4 2005–32 project.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Shoshi Kikuchi.

Additional information

Authors' contributions

KD designed the algorithm, did all the programming, and performed the feasibility test of the tool. AH helped to prepare test data sets and the literature search. TN supplied the inner database of known cis-elements to which the tool refers. KSa and KSu prepared the reference data. RM, MJM, and RB made many technical suggestions on the implementation and set up the host computer. RB also corrected the English of this manuscript. SK conceived the study and participated in its design and coordination. All authors read and approved the manuscript.

Electronic supplementary material

Additional file 1: Known plant cis-elements listed for analysis by RiCES. See text for further details. (XLS 104 KB)

12870_2007_228_MOESM2_ESM.csv

Additional file 2: Transcription units (TUs) used in the feasibility test. Auxin-inducible genes were picked up from RiceTFDB 2.0 (1st column). Corresponding full-length cDNAs were designated by BLASTN (2nd column) and translated to TUs defined in Pseudomolecule ver. 4 (3rd column). (CSV 2 KB)

12870_2007_228_MOESM3_ESM.txt

Additional file 3: Preliminary list of cis-element candidates listed by MEME analysis for TUs shown in Supplementary Table S2. See text for further details. (TXT 61 KB)

12870_2007_228_MOESM4_ESM.csv

Additional file 4: Result of association rule analysis of cis-element candidates listed by MEME. 1st column: examined sequence. 2nd column: number of TUs possessing the designated motif within 28 TUs of the target gene list. 3rd column: number of TU possessing the designated motif within 22 943 TUs stored in KOME database. 4th column: lift value. (CSV 100 KB)

12870_2007_228_MOESM5_ESM.csv

Additional file 5: Result of sequence search for motifs shown in Supplementary Table S1 in 22 943 TUs stored in KOME database. 1st column: examined TUs. 2nd column: motifs found in upstream region of TU. Other columns: position of motifs within the upstream region of each TU. (CSV 467 KB)

12870_2007_228_MOESM6_ESM.csv

Additional file 6: Result of association rule analysis after sequence search shown in Supplementary Table S6. 1st column: examined sequence. 2nd column: number of TUs possessing the designated motif within 28 TUs of the target gene list. 3rd column: number of TU possessing the designated motif within 22 943 TUs stored in KOME database. 4th column: lift value. (CSV 1007 bytes)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Doi, K., Hosaka, A., Nagata, T. et al. Development of a novel data mining tool to find cis-elements in rice gene promoter regions. BMC Plant Biol 8, 20 (2008). https://doi.org/10.1186/1471-2229-8-20

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2229-8-20

Keywords