DATFAP: A Database of Primers and Homology Alignments for Transcription Factors from 13 Plant Species
© Fredslund; licensee BioMed Central Ltd. 2008
Received: 25 October 2007
Accepted: 25 March 2008
Published: 25 March 2008
Many well-known transcription factor databases do not provide PCR primers for the sequences. Often, transcription factors from the same family have a very high sequence similarity, and so specific primer sets that only amplify their targets and none of the other family members may be hard to design manually. Also, it may often be useful to have one general primer set targeting the slightly different homologs of some transcription factor from several species.
DATFAP is a free, web-based, very user-friendly browsing tool based on a new database of more than 55,000 EST sequences from 13 plant species, classified as transcription factors. Further, the database offers primers designed for RealTime-PCR as well as homology alignments and phylogenies for the sequences. The provided PCR primers are designed so that they have a perfect sequence alignment to their target only and not to any other sequences in the database from the same species. Via a direct link to a helper tool, the user may also design a general primer set targeting all sequences in any hand-picked group of homologs. A sophisticated search facility enables the user to find exactly the relevant sequences which subsequently may be easily downloaded. All homologies among the more than one billion possible pairwise sequence comparisons of DATFAP have been charted in advance. Thus, the user may quickly display the alignment of any sequence and all its database homologs.
Because of the comprehensive homology analysis, it is very easy to find related transcription factors from different species, i.e. to navigate the network of inter-related transcription factors from the different species, and to find specific or general primers for them. DATFAP is found at the project homepage.
Transcription factors are proteins acting to regulate expression of genes by binding to their regulatory regions. Activated in signal transduction pathways, they are essential in the production of the right cellular response to some stimulus. Hence, specific transcription factors are the topic of study in hundreds of research projects, including projects dealing with plant breeding. Such studies can be very different in nature, ranging from the precise, functional characterization of certain proteins  over identification and classification of families of related proteins  to genetic modification .
When examining cDNA samples for the presence or expression levels of transcripts, one very often uses the Polymerase Chain Reaction, PCR. As many transcription factor families exist with minor variations among members, detecting one particular member requires uniquely specific PCR primers capable of binding to the target only, and avoiding the (often highly similar) relatives. Under other circumstances it might be relevant to have one general set of primers targeting a certain transcription factor such that the same primer set would amplify the correct target in several different species; a sort of orthogonal primer design problem.
In the search of new or unannotated transcription factors, sequence similarity to known transcription factors is a strong indicator. A large number of families have been identified, among plants in particular in rice and the model plant Arabidopsis thaliana, and homologs in other species to proteins from these families have been found.
This paper presents DATFAP (Database of Transcription Factors with Alignments and Primers), a resource where transcription factors from 12 plants and the green algae Chlamydomonas reinhardtii are collected and made available together with specific primers, homology alignments and sequence phylogenies.
Construction and Content
The EST sequences in DATFAP originate primarily from the Institute of Genomic Research (consensus sequences of EST alignments ), but also from the National Center for Biotechnology Information , the Database of Arabidopsis Transcription Factors , PlnTFDB , and PlantTFDB . A sequence is classified as a transcription factor in DATFAP if one or more of the following conditions are met:
• [keyword] The phrase "transcription factor" or one of the known transcription factor family names defined for rice or Arabidopsis is found in its header.
• [external database] It is listed in one of the above-mentioned external databases of transcription factors.
• [gene ontology] It has a Gene Ontology  annotation containing the phrase "transcription factor" or one of the known transcription factor family names defined for rice or Arabidopsis.
• [homology] It is very similar to a DATFAP sequence from another species classified by the above criteria as a transcription factor. Similarity is determined using a standalone version of the program BLAST  with e-value 10-7.
Besides its header, each sequence is associated with a classification record which states the reason(s) for its classification as a transcription factor in DATFAP.
The transcription factor family keywords used are: abi3, abi3/vp1, abi3vp1, alfin, alfin-like, ap2, ap2-erebp, ap2/erebp, arf, arid, arr-b, as2, at-hook, atrkd, aux, aux-iaa, bbr, bbr-bpc, bbr/bpc, bes1, bhlh, bpc, bzip, bzr, c2c2-co, c2c2-co-like, c2c2-dof, c2c2-gata, c2c2-yabby, c2h2, c3h, c3h-type 1, c3h-type 2, camta, cbf5, ccaat-dr1, ccaat-hap2, ccaat-hap3, ccaat-hap5, cpp, dbp, dp, e2f, e2f-dp, e2f/dp, eil, erebp, erf, fha, g2, g2-like, garp-arr-b, garp-g2-like, gata, gebp, gif, gras, grf, hb, hmg, hmg-box, homeobox, hrt, hrt-like, hsf, iaa, jumonji, lfy, lim, lim-domain, lug, mads, mbf1, myb, myb-related, nac, nin-like, nlp, nozzle, nsp1, nsp2, nzz, orphan, orphans, pbf-2-like, pcg, phd, phd-finger, platz, pseudo arr-b, pti4, pti5, pti6, rav, rem, rwp-rk, s1fa-like, sap, sbp, sigma70-like, sir2, srs, sw13, swi4/swi6, taz, tcp, tga3, tlp, trihelix, tub, ult, vip3, voz, voz-9, vp1, whirly, wrky, zf-hd, zim [6–8, 11, 12].
Utility and Discussion
Following the above-mentioned criteria, 55,361 out of more than 750,000 downloaded sequences were classified as transcription factors and curated in DATFAP. The access point for DATFAP is the front page where the user gives directives to the search engine via input to four text boxes. First of all, in the leftmost selection box, one or more species must be chosen among
• Medicago truncatula
• Zea mays (maize)
• Oryza sativa (rice)
• Lycopersicon esculentum (tomato)
• Chlamydomonas reinhardtii
• Solanum tuberosum (potato)
• Populus trichocarpa (poplar)
• Glycine max (soybean)
• Hordeum vulgare (barley)
• Lotus japonicus
• Arabidopsis thaliana
• Gossypium hirsutum (cotton)
• Triticum aestivum (wheat)
In the text box entitled "Search keywords", the user may type any keyword(s) of particular interest, e.g. homeobox. It is also possible simply to type an ID such as aw694210. Only sequences whose headers or classification records contain at least one of the search keywords will be selected. The search is case-insensitive.
Complex search queries can be composed using OR, AND, NOT, *, single quotes and parentheses; e.g.:
'zinc finger' and ('purine-binding' or CCCH) and not (partial or putative)
Some further details about the search facility:
• Parentheses are used to build sub-expressions
• Single quotes are used around search words or phrases containing spaces
• The wildcard character * is used to search for parts of words. E.g., trans* matches any word beginning with the substring "trans"; *nin* matches any string containing the substring "nin".
• AND, NOT and OR have their obvious meanings. E.g., not trans* matches anything not containing a word beginning with "trans".
Another example of a search phrase might be
*trans* and not transcription
This matches sequences containing the substring "trans" but not the word "transcription". The and is important here; if left out, the phrase is interpreted as *trans* or not transcription – and that matches all sequences! The Help pages show several more examples of search phrases (available from the front page).
In the Classification criteria selection box, the search may be restricted such that only transcription factors classified by certain criteria are searched.
The classification record of a sequence states its reason for being in DATFAP; in other words, it states by which criteria this sequence has been classified as a transcription factor (see Materials and Methods).
In case primers are provided for a sequence, the Primers table column contains a link to a page giving details about them, including GC content, melting temperature, length and sequence of each primer, as well as the expected amplicon length and a score for the primer set (since primers are designed for ESTs and not full gene sequences, they may cross intron/exon boundaries. However, in practice this is probably not the reason why a primer pair might fail to yield a specific product, as observed in ). The score of a primer set depends on whether the individual primers have 1- or 2-mismatch alignments to any non-target sequences. Such "close match" alignments may be considered problematic since they represent potential mis-primings. A score of 3000 indicates that the primer pair has no 1- or 2-mismatch alignments to non-target sequences. Further details on the score can be found in .
All DATFAP primers are automatically designed using a standalone version of the tool primique  which designs primers for groups of sequences such that each primer pair uniquely and specifically amplifies its target sequence and no other sequence from the group. Thus, no primer pair in DATFAP has a perfect match among non-target DATFAP sequences from the target species.
Selected parameter settings for round one in the primer design process.
18 – 22 nucleotides
80 – 120 nucleotides
Primer melting temperature
58.0 – 60.0 °
Max. primer melting temperature diff.
Enforce G/C terminal
40.0 – 60.0%
Max. base repeat
2 (e.g., GGG in a primer is disallowed)
Force specificity for both primers
Yes (i.e., none of the two primers has a perfect alignment to a non-target)
Parameter settings relaxed for round two in the primer design process.
Product length, round 2
80 – 160 nucleotides
Primer melting temperature, round 2
50.0 – 60.0 °
Max. base repeat, round 2
3 (e.g., GGG in a primer is allowed)
Parameter settings relaxed for round three in the primer design process.
Primer length, round 3
16 – 28 nucleotides
Product length, round 3
80 – 200 nucleotides
Max. primer melting temperature diff., round 3
Max. base repeat, round 3
4 (e.g., GGGG in a primer is allowed, although not in the 3'-end)
Thus, first all primers found in round 1 are stored. Then for round 2, primique is run with relaxed criteria on the sequences for which no valid primer sets were found in round 1. And finally, sequences still without primers are run through round 3. Some sequences are still without primers after round 3, most often because no primers can be found that align perfectly only to their target, or because the allowed product length is quite limited. The fraction of sequences for which primers could be found ranges from for 72.5% for Arabidopsis to 94.5% for poplar.
From the main results page (Figure 1), the user also has the opportunity to Create new primers for all sequences using primique. Clicking this link uploads all the result sequences to primique and lets the user decide all primer design parameter settings freely. The primers thus found are guaranteed to be specific among the uploaded sequences only, not among the entire set of same-species sequences from DATFAP as are the primers provided a priori with DATFAP.
The homologies between all sequences in DATFAP have been charted. This means that for any sequence resulting from a search, all its homologs can be displayed. For DATFAP purposes, two sequences A and B are defined to be homologous if and only if:
• blastn with e-value 10-7 finds significant sequence similarity between A and B
• the total length of all regions in A that are similar to regions in B must be at least 60% of the total length of A, or vice versa.
Thus, in order for A and B to be defined as homologs, one must "cover" at least 60% of the other with high similarity. This measure naturally induces a graph structure over all sequences in DATFAP, called the TF graph, where the sequences are the nodes of the graph, and an edge between two nodes represents a homology relationship between them. Starting at any node in the graph, one can therefore follow its edges to reach its homologs – if it has any. DATFAP provides an easy means of navigating the graph this way, described below.
Phylogeny images are created by PhyFi . By clicking the image, the user can open PhyFi in a new window, resize or otherwise manipulate the image and save it in a number of formats.
After the user has performed some search, frame 1 of Figure 5 shows the resulting sequences from Medicago (orange balls), cotton (purple balls), tomato (green balls) and rice (light green balls). The user may now choose to view the specific primers provided for, say, one of the rice sequences (frame 2). Next, going back to the main results page (frame 1), the user may click to view the homologs of the same rice sequence.
This leads to a page (frame 3) showing the rice sequence (light green ball) and its four DATFAP homologs from Chlamydomonas (ice blue ball), Arabidopsis (red ball) and maize (yellow balls). In fact, what is shown is a part of one of the connected components of the huge TF graph, namely the rice node and the four nodes it has edges to. Frame 3 shows other edges and nodes in the graph in shaded gray.
Now the user may choose to create a search from the group of five homologs currently showing. This leads to a new main results page (frame 4) with the rice, Chlamydomonas, Arabidopsis and maize sequences. Again, the user can choose to view the homologs of some selected sequence, this time a maize sequence (upper yellow ball), and frame 5 then shows this sequence and its five homologs from Arabidopsis (red ball), wheat (lilac ball), tomato (green balls) and rice (previously selected light green ball). Comparing to frame 3, this is simply a different part (subtree) of the TF graph with the yellow maize node as root.
Finally, the user can select some of the current sequences (rice, wheat and tomato sequences; yellow, lilac and green balls) and design a general set of primers targeting all of these sequences (frame 6). This is helpful if the user is studying homologs of the same transcription factor in several species: Then one such general primer set would suffice to amplify the gene in all species.
Since sequences are collected from several sources, duplicates do exist and have been removed from DATFAP. Two sequences are considered duplicates either if they are identical, or if they are at least 99.9% similar (according to blastn with e-value 10-7) while one is completely covered by the other (minus at most 10 nucleotides). Priorities are given to the classification criteria: if sequences A and B are considered duplicates while A originates from an external transcription factor database and B is solely in DATFAP because of homology to another DATFAP sequence, A is kept and B is removed since higher "transcription factor credibility" is given to external database origin than to homology.
The primer pairs are designed according to stringent criteria regarding melting temperature, GC content etc., and first and foremost according to the principle that a primer pair must have a perfect alignment to its target sequence and to no other sequences from the same species in the database. Primer design and PCR is an experimental science, and no guarantees can be given that primers designed by a computer program will work, regardless of the program. Good algorithms show a very high success rate, though; e.g., in the extensive studies by Czechowski and colleagues  and Caldana and colleagues , 96% and 97% of all reactions yielded specific products, respectively. The primers of DATFAP are designed in a way which is similar to what was done in these studies. This author has designed several free, computerized, web-based primer design tools in the past according to exactly the same principles as those used for the DATFAP primers (e.g. [19, 15]). One, PriFi, consistently has 25–50 users each day from more than 30 countries all over the world, most of whom are returning customers. In this author's view, the satisfaction of these users warrants the approach.
The projects described in [13, 19] are Real-Time PCR platforms for Arabidopsis and Rice transcription factors, respectively, offering primers for the purposes of high-throughput expression profiling. Moreover, extensive experimental work utilizing the primers resulted in important insight and new knowledge regarding, e.g., organ-specific expression of certain genes. Such were not the primary goals of the DATFAP project. Through its suite of 13 species and easy means of searching and relating sequences within or between species, and through its predefined specific primers complemented by its option for creating cross-species primers, DATFAP is a general-purpose tool, hopefully paving the way for many different kinds of research centered around transcription factors.
DATFAP is a user-friendly, web-based, graphical resource freely available to anyone with an Internet browser. Equipped with a sophisticated search facility and specific primers for (almost) all sequences, DATFAP constitutes a valuable tool to researchers in all areas of plant molecular biology working with transcription factors. No other multi-species transcription factor database offers such easy inter- and intra-species navigation in the network of related transcription factors.
The DATFAP web site includes an elaborate Help page which illustrations, details and examples of the search facility.
Availability and Requirements
Project name: "DATFAP: A Database of Primers and Homology Alignments for Transcription Factors from 13 Plant Species"
Project home page: http://cgi-www.daimi.au.dk/cgi-chili/datfap/frontdoor.py
Operating system(s): Platform independent (web-based)
Programming language: Python
Licence: Free to all non-commercial users
Any restrictions to use by non-academics: Licence needed.
I would like to thank Lene H. Madsen, Jens Stougaard, Niels Sandal and Mette Lange for their advice and support with this work. Mostly I wish to thank everyone at BiRC for five great years.
- Miller G, Mittler R: Could heat shock transcription factors function as hydrogen peroxide sensors in plants?. Ann Bot (Lond). 2006, 98 (2): 279-88. 10.1093/aob/mcl107.View ArticleGoogle Scholar
- Toledo-Ortiz G, Huq E, Quail PH: The Arabidopsis basic/helix-loop-helix transcription factor family. Plant Cell. 2003, 15 (8): 1749-70. 10.1105/tpc.013839.PubMedPubMed CentralView ArticleGoogle Scholar
- Memelink J, Kijne JW, van der Heijden R, Verpoorte R: Genetic modification of plant secondary metabolite pathways using transcriptional regulators. Adv Biochem Eng Biotechnol. 2001, 72: 103-25.PubMedGoogle Scholar
- Lee Y, Tsai J, Sunkara S, Karamycheva S, Pertea G, Sultana R, Antonescu V, Chan A, Cheung F, Quackenbush J: The TIGR Gene Indices: clustering and assembling EST and known genes and integration with eukaryotic genomes. Nucleic Acids Research. 2005, D71-D74. 33 DatabaseGoogle Scholar
- National Center for Biotechnology Information. [http://www.ncbi.nlm.nih.gov/]
- Guo A, He K, Liu D, Bai S, Gu X, Wei L, Luo J: DATF: a Database of Arabidopsis Transcription Factors. Bioinformatics. 2005, 21: 2568-2569. 10.1093/bioinformatics/bti334.PubMedView ArticleGoogle Scholar
- Riano-Pachon DM, Ruzicic S, Dreyer I, Mueller-Roeber B: PlnTFDB: An integrative plant transcription factor database. BMC Bioinformatics. 2007, 8: 42-10.1186/1471-2105-8-42.PubMedPubMed CentralView ArticleGoogle Scholar
- Guo AY, Chen X, Gao G, Zhang H, Zhu QH, Liu XC, Zhong YF, Gu X, He K, Luo J: PlantTFDB: a comprehensive plant transcription factor database. Nucl Acids Res. 2008, D966-D969. 36 DatabaseGoogle Scholar
- Gene Ontology. [http://geneontology.org/]
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucl Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.PubMedPubMed CentralView ArticleGoogle Scholar
- Iida K, Seki M, Sakurai T, Satou M, Akiyama K, Toyoda T, Konagaya A, Shinozaki K: RARTF: Database and Tools for Complete Sets of Arabidopsis Transcription Factors. DNA Research. 2005, 12: 247-256. 10.1093/dnares/dsi011.PubMedView ArticleGoogle Scholar
- Riechmann JL, Heard J, Martin G, Reuber L, Jiang C, Keddie J, Adam L, Pineda O, Ratcliffe OJ, Samaha RR, Creelman R, Pilgrim M, Broun P, Zhang JZ, Ghandehari D, Sherman BK, Yu G: Arabidopsis transcription factors: genome-wide comparative analysis among eukaryotes. Science. 2000, 290: 2105-2110. 10.1126/science.290.5499.2105.PubMedView ArticleGoogle Scholar
- Czechowski T, Bari RP, Stitt M, Scheible WR, Udvardi M: Real-time RT-PCR profiling of over 1400 Arabidopsis transcription factors: unprecedented sensitivity reveals novel root- and shoot-specific genes. Plant J. 2004, 38 (2): 366-79. 10.1111/j.1365-313X.2004.02051.x.PubMedView ArticleGoogle Scholar
- Fredslund J, Lange M: Primique: automatic design of specific PCR primers for each sequence in a family. BMC Bioinformatics. 2007, 8: 369-10.1186/1471-2105-8-369.PubMedPubMed CentralView ArticleGoogle Scholar
- Corpet F: Multiple sequence alignment with hierarchical clustering. Nucl Acids Res. 1988, 16 (22): 10881-10890. 10.1093/nar/16.22.10881.PubMedPubMed CentralView ArticleGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignments through sequence weighting, position specific gap penalties and weight matrix choice. Nucl Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.PubMedPubMed CentralView ArticleGoogle Scholar
- Fredslund J: PHYFI: fast and easy online creation and manipulation of phylogeny color figures. BMC Bioinformatics. 2006, 7: 315-10.1186/1471-2105-7-315.PubMedPubMed CentralView ArticleGoogle Scholar
- Fredslund J, Schauser L, Madsen LH, Sandal N, Stougaard J: PriFi – Using a Multiple Alignment of Related Sequences to Find Primers for Amplification of Homologs. Nucl Acids Res. 2005, W516-W520. 10.1093/nar/gki425. doi:10.1093/nar/gki425, 33 Web ServerGoogle Scholar
- Caldana C, Scheible W-R, Mueller-Roeber B, Ruzicic S: A quantitative RT-PCR platform for high-throughput expression profiling of 2500 rice transcription factors. Plant Methods. 2007, 3: 7-10.1186/1746-4811-3-7.PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.