- Open Access
SoyDB: a knowledge database of soybean transcription factors
BMC Plant Biology volume 10, Article number: 14 (2010)
Transcription factors play the crucial rule of regulating gene expression and influence almost all biological processes. Systematically identifying and annotating transcription factors can greatly aid further understanding their functions and mechanisms. In this article, we present SoyDB, a user friendly database containing comprehensive knowledge of soybean transcription factors.
The soybean genome was recently sequenced by the Department of Energy-Joint Genome Institute (DOE-JGI) and is publicly available. Mining of this sequence identified 5,671 soybean genes as putative transcription factors. These genes were comprehensively annotated as an aid to the soybean research community. We developed SoyDB - a knowledge database for all the transcription factors in the soybean genome. The database contains protein sequences, predicted tertiary structures, putative DNA binding sites, domains, homologous templates in the Protein Data Bank (PDB), protein family classifications, multiple sequence alignments, consensus protein sequence motifs, web logo of each family, and web links to the soybean transcription factor database PlantTFDB, known EST sequences, and other general protein databases including Swiss-Prot, Gene Ontology, KEGG, EMBL, TAIR, InterPro, SMART, PROSITE, NCBI, and Pfam. The database can be accessed via an interactive and convenient web server, which supports full-text search, PSI-BLAST sequence search, database browsing by protein family, and automatic classification of a new protein sequence into one of 64 annotated transcription factor families by hidden Markov models.
A comprehensive soybean transcription factor database was constructed and made publicly accessible at http://casp.rnet.missouri.edu/soydb/.
Soybean is a great source of protein, as it contains significant amounts of all the essential amino acids, including some that cannot be synthesized by the human body . Soybean has been used as a food and a drug component in China for thousands of years  and over the past 60 years has become a leading crop in many nations around the world . Because of its high value in the agricultural and food industry, soybean has received greater and greater research attention, both to improve soybean agronomic performances and as a model for basic biological studies. In early 2008, the Department of Energy-Joint Genome Institute (DOE-JGI) finished sequencing the soybean genome using a whole-genome shotgun approach , which makes soybean the most complex plant so far ever sequenced . The homology-based gene prediction and annotation produced putative protein sequences [4, 5], which makes it feasible to identify and annotate soybean transcription factors.
Transcription factors (TF) are proteins that bind to DNA sequences (i.e., promoters) and regulate gene expression by one or more DNA binding domains. Virtually all biological processes are directly regulated or influenced by transcription factors . For example, the transcription process in eukaryotes would not occur in the absence of a specific class of transcription factors named "general transcription factors" [7, 8]. Studies have shown that transcription factors are closely involved in the process of cell development, such as cellular division, migration, and differentiation . Transcription factors of Arabidopsis thaliana have been well studied since its genome has been fully sequenced as a model specie [6, 10–14]. This makes it possible to identify and study transcription factors of other newly sequenced species, such as soybean, by homology searching and comparative analysis.
Several databases for soybean genome analysis have been built and made publicly available, such as SoyGD , SoyBase , and SoyXpress . These databases contain a variety of information, such as soybean genome sequences, bacterial artificial chromosome (BAC), expressed sequence tags (EST), and some useful tools including genome browsers, BLAST searching, and pathway searching. However, these databases only contain general annotations for the soybean genome, instead of knowledge specifically targeting the transcription factors. For example, none of them systematically organizes transcription factors into families or clearly points out the DNA binding domains. PlantTFDB  and DBD  are two existent transcription factor databases, which contain knowledge about transcription factors from multiple species. For each soybean transcription factor, PlantTFDB contains information including protein sequence, Gene Ontology annotation , putative binding domains found by InterProScan , and cross-links to external databases, including EMBL , UniProt , RefSeq , and TRANSFAC . DBD contains the amino acid sequence of each transcription factor and external links to Ensembl , Pfam , and SUPERFAMILY . Compared to PlantTFDB, DBD has less external database links, but DBD claims to contain the transcription factors of 927 completely sequenced genomes whereas PlantTFDB covers 22 species. The knowledge in these two databases is very useful; however, they were built based on a relatively older version of soybean sequence data and their annotations are still incomplete. The most important component they lack is the three dimensional structure for each transcription factor, because the visualization of the transcription factor, especially binding sites, can help further understanding the mechanism and functions of transcription factors, which is indispensible to structural genomics [29, 30]. Furthermore, with the complete genome sequences of more and more species available, a computer system is needed that can automatically generate a knowledge database and publish it with a user-friendly interface.
To fill the gap, we developed SoyDB - a comprehensive and integrated database for soybean transcription factors. This database not only contains most of the content and features already existed in PlantTFDB and DBD, but also extends them by containing more comprehensive knowledge and links to more versatile external datasets. The annotations in SoyDB include predicted tertiary structures, protein domains, multiple sequence alignments, DNA binding sites, and web logos and consensus sequences for each family. The SoyDB database also contains links to the homologous EST sequences, and the same or homologous proteins in other databases including PlantTFDB, PDB , Swiss-Prot , TAIR , RefSeq , SMART , Pfam , KEGG , SPRINTS , EMBL , InterPro , PROSITE , and Gene Ontology .
Moreover, our system can automatically execute bioinformatics tools and generate annotations, link to other well-known protein databases, construct MySQL databases, and generate PHP scripts to build its website. This fully automated approach can be used to create a protein annotation database and website for any sequenced organism in the future.
Construction and Content
SoyDB contains the annotations of 5,671 putative transcription factors. These proteins were classified into 64 families (for details see the section "Transcription Factor Family Prediction Using SAM Hidden Markov Models"). Figure 1 illustrates the architecture of the SoyDB website. Users can access the main components from the home page: full-text search, PSI-BLAST sequence search, family classification by hidden Markov model, family browsing, protein browsing, family information, protein information, and FTP site.
The soybean genome sequences and gene model predictions used in this study were acquired from the publicly available database Phytozome . These sequences were generated by the preliminary GenomeScan , FgenesH , and PASA  gene annotations based on the Gm1.01 version of soybean assembly data .
Transcription Factor Identification
We used the standalone versions of InterProScan  to search all the soybean protein sequences against 11 databases integrated in InerPro . These databases and their corresponding scanning methods include: PROSITE (pfscan) , PRINTS (FingerPRINTScan) , Pfam (HMMPfam) , ProDom (ProDomBlast3i) , SMART (HMMSmart) , TIGRFAMs (HMMTigr) , PIR SuperFamily (HMMPIR) , SUPERFAMILY (superfamily) , Gene3D (gene3d) , PANTHER (HMMPanther) , and HAMAP (pfscan) . InterProScan systematically searches each of these databases using their corresponding scanning methods to find domains. The proteins predicted to contain TF related domain(s) were considered as putative transcription factors. Using the Plant Transcription Factor Database (PlnTFDB)  and the classification of Medicago truncatula TF genes  as references, we manually curated the list of putative transcription factors and eliminated any mistakenly identified sequences. In this way, we identified 5,671 putative TF sequences.
Transcription Factor Family Prediction Using SAM Hidden Markov Models
The transcription factors of Arabidopsis thaliana have been well studied and classified into 64 families . This provides a model for us to classify soybean transcription factors. We used MUSCLE  to generate a multiple sequence alignment for each Arabidopsis thaliana TF family. The multiple sequence alignment was then input into SAM 3.5  to build a hidden Markov model (HMM) for each family. Every soybean TF sequence was aligned with each of the 64 HMMs, which outputs an e-value. This e-value can be considered as a fitness score between a TF sequence and a hidden Markov model: lower e-value indicates better fitness. Finally, a transcription factor was classified into the family whose HMM yields the lowest e-value.
Annotations Using Bioinformatics Tools
The standalone versions of several bioinformatics tools were locally installed and executed to generate annotations for soybean transcription factors. An accurate protein structure prediction tool MULTICOM  was used to predict the tertiary structure of each transcription factor when homologous template proteins could be found in PDB. According to the official evaluations of the 8th community-wide Critical Assessment of Techniques for Protein Structure Prediction (CASP8) , MULTICOM was able to predict high-accuracy tertiary structures with an average GDT-TS score  0.87 if suitable templates can be found. GDT-TS score ranges from 0 to 1 measuring the similarities between the predicted and experimental structures, whereas 1 indicates completely the same and 0 completely different. Figure 2 illustrates the predicted tertiary structure of a transcription factor in SoyDB with ID GM00002, and the electrostatic polarization of the predicted structure. The blue area in the electrostatic polarization shows residues positively charged, which is found to be highly identical to the green area in Figure 2(b), which is the putative DNA-binding sites identified by a pair-wise alignment between GM00002 and its template protein 1WID (Figure 2(d)). Since it has been studied and found that the DNA-binding area is positively charged if analyzed by electrostatic potentials , the highly identical area in Figure 2(a) and 2(b) strongly confirms that the predicted structure has the electrostatic properties of a transcription factor. This further confirms the qualities of MULTICOM predictions and the correctness of the predicted binding sites derived from the homology alignment. In SoyDB, a predicted tertiary structure is visualized by Jmol . In order to clearly visualize the tertiary structure of the DNA-binding region, only the segments containing homologous DNA binding domains are visualized by Jmol. Users can view a TF structure from various perspectives in a three-dimensional way and perform many operations including selecting and highlighting interested regions, changing view styles and colors, and measuring atom distances and angles by right clicking on the Jmol console and selecting corresponding menus. Detailed instructions about Jmol menus can be found at Jmol website . During the structure prediction process, MULTICOM generates the sequence alignments between the transcription factor and its homologous templates using PSI-BLAST. These sequence alignments can be used to predict the binding sites of a transcription factor based on the experimentally determined binding sites of its template as shown in Figure 2.
A predicted structure was parsed into domains by Protein Domain Parser (PDP) . Since some transcription factors did not have homologous templates found in PDB, DOMAC , an accurate ab initio domain prediction tool, was also used to predict domains for each transcription factor.
The protein sequences in the same family were aligned by MUSCLE  and visualized by WebLogo . A consensus sequence was derived from the multiple sequence alignment by selecting the most frequently appeared amino acid at each position. The multiple sequence alignments were also used to identify the conserved signatures (likely the DNA binding domains) for each family.
All of the bioinformatics tools incorporated to construct SoyDB can be used to automatically annotate other species in the future.
Links to External Databases and Datasets
In order to annotate the functions of soybean transcription factors, each TF protein sequence was searched against the soybean TF database PlantTFDB, NCBI known EST sequences, and other general protein databases by PSI-BLAST or TBLASTN. The external protein databases include Swiss-Prot , TAIR , RefSeq , SMART , Pfam , KEGG , SPRINTS , EMBL , InterPro , PROSITE , and Gene Ontology . Web links to these databases were created when the same transcription factor or its homologous proteins were found in them; and for each database or EST dataset only the PSI-BLAST or PBLASTN hit with the smallest e-values was listed in SoyDB. To search the known EST sequences, PSI-BLAST was first used to build a position-specific score matrix for each transcription factor. TBLASTN was then used to search each protein sequence against three known EST datasets: EST human, EST mouse, and EST others. GenBank  web page of each EST hit was linked to SoyDB website. The gene expression of a subset of TF genes (about 1,000 TF genes) was recently published . Transcription profile of all soybean TF genes in various conditions is under investigation.
These external links greatly expand the annotation scope of SoyDB providing related knowledge from various perspectives. SoyDB provides a systematic view of a transcription factor -- from the features of the protein itself, to the biological pathway it locates in. The links to the external databases and datasets can be updated by a re-run of PSIBLAST and TBLASTN. Currently, these links are scheduled to be updated once every six months. This time interval can be changed if necessary.
Database and Website Implementation
The programs used to automatically annotate proteins were written in PERL. The relational database was built on MySQL with database schemas automatically generated by programs written in PERL. The website was implemented in PHP. The database and web site were automatically constructed by computer programs with little human intervention.
Utility and Discussion
This component contains the complete annotations for each transcription factor, including protein ID, protein name and description, tools used for TF identification, family ID, family name and description, amino acid sequence, homology domain prediction, ab initio domain prediction, PDB homologous templates, and predicted tertiary structure. This component can be reached by clicking the sequence ID, such as GM00001, or the Phytozome protein name, such as Glyma01g11670.1, at the "Protein Browsing" webpage (for details see the following "Protein Browsing" section). Figure 3 illustrates the protein information page. The sequence ID and family ID, such as GM00001 and GMF0001, are internal indices used by the SoyDB, and the sequence name is the standard soybean TF name used by the soybean genome database Phytozome . We noticed the trend of unifying annotation formats within the soybean community. Therefore, the commonly used TF ID format, such as PTGm00009.1, is also compatible in SoyDB. Details are described in the "Full-Text Search" section below.
This component contains the complete annotations for each TF family, including family ID, family name and description, number of sequences within the family, consensus sequence, consensus signatures (likely the DNA binding regions), web logo of the signature profile, and multiple sequence alignment of the protein sequences within the family. Figure 4 demonstrates a family information web page. This component can be reached from the "Family Browsing" web page.
The transcription factors within a family are listed in the order of sequence IDs. The list contains the thumbnail of tertiary structure, protein ID and name, family ID, and family name of each transcription factor (Figure 5). Users can further view the complete annotations by clicking its sequence ID or the Phytozome protein name. This component can be reached by clicking the number of sequences in the "Family Browsing" or the "Family Information" web page.
A user can browse SoyDB from TF family perspective. The 64 TF families are listed in the order of family IDs. The family ID, family name, and the number of transcription factors within each family are listed. By clicking the family ID or name, users can view the complete annotations for a family, or further browse the sequences within a family by clicking the number of sequences. This component can be reached by clicking "Browse Database" in both the top and bottom menu bars from any SoyDB web pages. Additional file 1 (Figure S1) illustrates the web page showing a TF family list.
Full Text Search
This component allows users to search the entire SoyDB database by a query text, such as protein name or family name. Given input keywords, SoyDB searches all the fields in the database and returns matched transcription factors with links to their annotations. Users can also search SoyDB by the TF IDs used in PlantTFDB. The search component will return the homologous soybean TFs found in SoyDB.
PSI-BLAST Sequence Search
This component allows users to search a query sequence against every TF sequence stored in SoyDB. Users can submit a query sequence and adjust PSI-BLAST parameters from a web page as shown in Additional file 2 (Figure S2). After a PSI-BLAST search is performed, the significant hits, with links to their annotation web pages, are ranked based on the e-values generated by PSI-BLAST. Additional file 3 (Figure S3) illustrates a PSI-BLAST result web page.
Family Classification by Hidden Markov Model
This component classifies a query protein sequence into one of the 64 TF families. Additional file 4 (Figure S4) illustrates the web page for family classification. A submitted query sequence is aligned with each of the 64 hidden Markov models built based on the 64 Arabidopsis thaliana TF families. The query sequence is classified into a family whose hidden Markov model outputs the lowest e-value (correspondingly the highest alignment score or fitness score). More details about family classification can be found at the "Transcription Factor Family Prediction Using SAM Hidden Markov Models" section under "Construction and Content".
All of the information in SoyDB is available for users to download from an FTP site. For example, users can download all of the TF protein sequences in the FASTA format and the multiple sequence alignments for each family in plain text. This makes it possible for other websites to link to SoyDB by performing PSI-BLAST searches on SoyDB sequences, similarly as SoyDB links with other external databases.
Comparisons and Overlapping between SoyDB and PlantTFDB
In total, SoyDB has 5,671 transcription factors - 4,306 of them (75.9%) have hits found in PlantTFDB identified by PSI-BLAST with an e-value threshold 10-3. PlantTFDB has 1,891 soybean transcription factors (based on the FASTA file downloaded from PlantTFDB FTP site), and 1,805 of them (95.5%) have hits found from SoyDB based on a PSI-BLAST search with an e-value threshold 10-3.
Comparisons of Soybean Transcription Factor Family Distributions with Other Plants
The collection and analyses in SoyDB allows us to perform comparisons of TF family distributions across the plant kingdom. The large number of soybean TF genes (5,671) described in this study is likely due to the two soybean whole genome duplication events; one estimated to have occurred 40-50 million years ago (mya) and the most recent one approximately 10-15 million years ago [64, 65]. By comparing the total number of genes in different organisms, it was found that the increase of plant gene number is related to multicellularity and ploidy. For example, compared to the unicellular eukaryote Chlamydomonas reinhardtii where 15,143 genes were predicted , larger numbers of protein-encoding genes were reported in multicellular plant organisms, e.g., Physcomitrella patens (35,938; ), Arabidopsis thaliana (32,944; TAIR ), and the tetraploid Glycine max (66,153, Phytozome). We hypothesize that TF gene number also follows the same trend as land plants, which have a larger number of TF genes compared to algae. To perform comparisons of plant TF genes and their distributions across TF gene families, we mined the last updated DBD database  for 11 plant species (C. reinhardtii, P. patens, Oryza sativa, Zea mays, Sorghum bicolor, Lotus japonicum, Medicago truncatula, A. thaliana, Vinis vinifera, Ricinus communis, and Populus trichocarpa). These species were then compared with the soybean TF genes stored in our SoyDB database.
Our analysis showed that the unicellular C. reinhardtii has the lowest number of TF genes compared to multicellular land plants (the exceptions are L. japonicus and M. truncatula where only a partial genome sequence is available). This trend also reflects the differences of total gene number in the organisms shown in Figure 6. For example, it is interesting to note that homeobox, MYB, NAC, and WRKY TF genes in C. reinhardtii lack or have very low representations compared to the 11 other plant models (Table 1). Previous studies defined a role for homeobox  and WRKY genes  in plant development. Therefore, the occurrence of these genes only in multicellular plants may reflect their special roles in development. In addition, a close relationship between TF gene number and total gene number  is observed when comparing the TF gene numbers of G. max and A. thaliana with their total gene numbers (i.e., G. max encodes 66,153 protein-coding genes including 5,683 TF genes; A. thaliana encodes 32,944 protein-coding genes and 1,738 TF genes). Thus, the family distribution of soybean TF genes is similar to other land plant species, except for P. patens (e.g., AP2 represents 7% of total TF genes in soybean vs. 8-12% for other land plants; bZIP: 3% vs. 3-7%; bHLH: 7% vs. 8-11%; homeobox: 6% vs. 4-7%; MYB: 14% vs. 7-14%; NAC: 4% vs. 4-9%; WRKY: 3% vs. 4-7%; ZF-C2H2: 7% vs. 5-9%) (Figure 6 and Table 1).
Collectively, these results suggest that soybean TF genes were not lost following soybean genome duplication, and may have evolved for specialized functions in plant development or response to the environment.
Future Development Plan
In the future, we plan to link to more soybean database, such as SoyBase, and add a human expert discussion section for each transcription factor where biologists can register, log in, and make comments on any annotation items. Also, we plan to link the protein name, such as Glyma01g11670.1, listed in each protein information page to its entry in Phytozome. By doing this, SoyDB can be linked with other soybean genome annotations. Furthermore, we may identify the binding regions on the soybean DNA sequences, which can further help biologists target the regulated regions on soybean genome.
SoyDB is a comprehensive database for soybean transcription factors. It integrates bioinformatics tools and various external databases to provide rich annotations, which can be browsed and retrieved through convenient web interfaces. The automated process generates annotations and creates database and website, and can be used to annotate other sequenced species.
Availability and Requirements
SoyDB is freely available at http://casp.rnet.missouri.edu/soydb/ for academic use. Based on our test, SoyDB is fully functional with three web browsers: Mozilla Firefox, Internet Explorer, and Safari, and four operating systems: Windows XP, Windows Vista, Linux (Red Hat), and Mac OS. The only system requirement for SoyDB is that JAVA runtime environment (JRE) needs to be installed and set fully functional in order to make Jmol work.
Henkel J: Soy: health claims for soy protein, questions about other components. FDA consumer. 2000, 34:
Han B-Z, Rombouts FM, Nout MJR: A Chinese fermented soybean food. International Journal of Food Microbiology. 2001, 65 (1-2): 1-10. 10.1016/S0168-1605(00)00523-7.
Carpenter J, Gianessi L: Agricultural biotechnology: updated benefit estimates. 2001, Washington, USA National Centre for Food and Agricultural Policy (NCFAP)
Schmutz J, Cannon S, Schlueter J, Ma J, Hyten D, Song Q, Mitros T, Nelson W, May G, Gill N, et al: Genome sequence of the paleopolyploid soybean (Glycine max (L.) Merr.). Nature.
Jakoby M, Weisshaar B, Droge-Laser W, Vicente-Carbajosa J, Tiedemann J, Kroj T, Parcy F: bZIP transcription factors in Arabidopsis. Trends in Plant Science. 2002, 7 (3): 106-111. 10.1016/S1360-1385(01)02223-3.
Reese J: Basal transcription factors. Current opinion in genetics & development. 2003, 13 (2): 114-118. 10.1016/S0959-437X(03)00013-3.
Weinzierl R: Mechanisms of gene expression: structure, function and evolution of the basal transcriptional machinery. 1999, London, UK Imperial College Press
Corrinne C: Transcription factors and mammalian development. Current Topics in Developmental Biology. Edited by: Pedersen RA. 1992, Academic Press, Inc, 27: 351-full_text.
Liu Q, Kasuga M, Sakuma Y, Abe H, Miura S, Yamaguchi-Shinozaki K, Shinozaki K: Two transcription factors, DREB1 and DREB2, with an EREBP/AP2 DNA binding domain separate two cellular signal transduction pathways in drought-and low-temperature-responsive gene expression, respectively, in Arabidopsis. The Plant Cell Online. 1998, 10 (8): 1391-1406. 10.1105/tpc.10.8.1391.
Riechmann J, Heard J, Martin G, Reuber L, Z C, Keddie J, Adam L, Pineda O, Ratcliffe O, Samaha R: Arabidopsis transcription factors: genome-wide comparative analysis among eukaryotes. Science. 2000, 290 (5499): 2105-2110. 10.1126/science.290.5499.2105.
Sakuma Y, Maruyama K, Osakabe Y, Qin F, Seki M, Shinozaki K, Yamaguchi-Shinozaki K: Functional analysis of an Arabidopsis transcription factor, DREB2A, involved in drought-responsive gene expression. The Plant Cell Online. 2006, 18 (5): 1292-1309. 10.1105/tpc.105.035881.
Johnson C, Kolevski B, Smyth D: TRANSPARENT TESTA GLABRA2, a trichome and seed coat development gene of Arabidopsis, encodes a WRKY transcription factor. The Plant Cell. 2002, 14 (6): 1359-1375. 10.1105/tpc.001404.
Ulker B, Somssich I: WRKY transcription factors: from DNA binding towards biological function. Current Opinion in Plant Biology. 2004, 7 (5): 491-498. 10.1016/j.pbi.2004.07.012.
Shultz J, Kurunam D, Shopinski K, Iqbal M, Kazi S, Zobrist K, Bashir R, Yaegashi S, Lavu N, Afzai A, et al: The Soybean Genome Database (SoyGD): a browser for display of duplicated, polyploid, regions and sequence tagged sites on the integrated physical and genetic maps of Glycine max. Nucleic Acids Research. 2006, D758-D765. 10.1093/nar/gkj050. 34 Database
SoyBase and the soybean breeder's toolbox. [http://soybase.org/]
Cheng K, Stromvik M: SoyXpress: a database for exploring the soybean transcriptome. BMC genomics. 2008, 9 (1): 368-10.1186/1471-2164-9-368.
Guo A, Chen X, Gao G, Zhang H, Zhu Q, Liu X, Zhong Y, Gu X, He K, Luo J: PlantTFDB: a comprehensive plant transcription factor database. Nucleic Acids Research. 2008, D966-D969. 36 Database
Wilson D, Charoensawan V, Kummerfeld S, Teichmann S: DBD taxonomically broad transcription factor predictions: new content and functionality. Nucleic Acids Research. 2007, D88-D92. 10.1093/nar/gkm964. 36 Database
Ashburner M, Ball C, Blake J, Botstein D, Butler H, Cherry J, Davis A, Dolinski K, Dwight S, Eppig J: Gene ontology: tool for the unification of biology. Nature Genetics. 2000, 25 (1): 25-29. 10.1038/75556.
Zdobnov E, Apweiler R: InterProScan-an integration platform for the signature-recognition methods in InterPro. Bioinformatics. 2001, 17 (9): 847-848. 10.1093/bioinformatics/17.9.847.
Stoesser G, Tuli M, Lopez R, Sterk P: The EMBL nucleotide sequence database. Nucleic Acids Research. 1999, 27 (1): 18-24. 10.1093/nar/27.1.18.
Apweiler R, Bairoch A, Wu C, Barker W, Boeckmann B, Ferro S, Gasteiger E, Huang H, Lopez R, Magrane M, et al: UniProt: the universal protein knowledgebase. Nucleic Acids Research. 2004, D115-D119. 10.1093/nar/gkh131. 32 Database
Pruitt K, Tatusova T, Maglott D: NCBI Reference Sequence (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Research. 2006, D1-D5. 00 Database
Wingender E, Chen X, Hehl R, Karas H, Liebich I, Matys V, Meinhardt T, Prub M, Reuter I, Schacherer F: TRANSFAC: an integrated system for gene expression regulation. Nucleic Acids Research. 2000, 28 (1): 316-319. 10.1093/nar/28.1.316.
Hubbard T, Barker D, Birney E, Cameron G, Chen Y, Clark L, Cox T, Cuff J, Curwen V, Down T: The Ensembl genome database project. Nucleic Acids Research. 2002, 30 (1): 38-41. 10.1093/nar/30.1.38.
Bateman A, Coin L, Durbin R, Finn R, Hollich V, Griffiths-Jones S, Khanna A, Marshall M, Moxon S, Sonnhammer E: The Pfam protein families database. Nucleic Acids Research. 2004, 32 (1): 276-280. 10.1093/nar/gkh121.
Madera M, Vogel C, Kummerfeld S, Chothia C, Gough J: The SUPERFAMILY database in 2004: additions and improvements. Nucleic Acids Research. 2004, D235-D239. 10.1093/nar/gkh117. 32 Database
Burley S, Almo S, Bonanno J, Capel M, Chance M, Gaasterland T, Lin D, Sali A, Studier F, Swaminathan S: Structural genomics: beyond the human genome project. Nature Genetics. 1999, 23: 151-158. 10.1038/13783.
Burley S: An overview of structural genomics. Nature Structural & Molecular Biology. 2000, 7: 932-934. 10.1038/80697.
Berman H, Westbrook J, Feng Z, Gilliland G, Bhat T, Weissig H, Shindyalov I, Bourne P: The protein data bank. Nucleic Acids Research. 2000, 28 (1): 235-242. 10.1093/nar/28.1.235.
Boeckmann B, Bairoch A, Apweiler R, Blatter M, Estreicher A, Gasteiger E, Martin M, Michoud K, O'Donovan C, Phan I: The SWISS-PROT protein knowledgebase and its supplement TrEMBL in 2003. Nucleic Acids Research. 2003, 31 (1): 365-370. 10.1093/nar/gkg095.
Rhee S, Beavis W, Berardini T, Chen G, Dixon D, Doyle A, Garcia-Hernandez M, Huala E, Lander G, Montoya M, et al: The Arabidopsis Information Resource (TAIR): a model organism database providing a centralized, curated gateway to Arabidopsis biology, research materials and community. Nucleic Acids Research. 2003, 224-228. 10.1093/nar/gkg076. 31
Letunic I, Copley R, Pils B, Pinkert S, Schultz J, Bork P: SMART 5: domains in the context of genomes and networks. Nucleic Acids Research. 2006, D257-D260. 10.1093/nar/gkj079. 34 Database
Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, Katayama T, Kawashima S, Okuda S, Tokimatsu T: KEGG for linking genomes to life and the environment. Nucleic Acids Research. 2008, D480-D484. 36 Database
Attwood T, Blythe M, Flower D, Gaulton A, Mabey J, Maudling N, McGregor L, Mitchell A, Moulton G, Paine K: PRINTS and PRINTS-S shed light on protein ancestry. Nucleic Acids Research. 2002, 30 (1): 239-241. 10.1093/nar/30.1.239.
Angiuoli S, Gussman A, Klimke W, Cochrane G, Field D, Garrity G, Kodira C, Kyrpides N, Madupu R, Markowitz V: Toward an online repository of standard operating procedures (SOPs) for (Meta) genomic annotation. OMICS: A Journal of Integrative Biology. 2008, 12 (2): 137-141. 10.1089/omi.2008.0017.
Mulder N, Apweiler R, Attwood T, Bairoch A, Bateman A, Binns D, Biswas M, Bradley P, Bork P, Bucher P: InterPro: An integrated documentation resource for protein families, domains and functional sites. Briefings in Bioinformatics. 2002, 3 (3): 225-235. 10.1093/bib/3.3.225.
Hulo N, Bairoch A, Bulliard V, Cerutti L, De Castro E, Langendijk-Genevaux P, Pagni M, Sigrist C: The PROSITE database. Nucleic Acids Research. 2006, D227-D230. 10.1093/nar/gkj063. 34 Database
Yeh R, Lim L, Burge C: Computational inference of homologous gene structures in the human genome. Genome research. 2001, 11 (5): 803-816. 10.1101/gr.175701.
Haas B, Delcher A, Mount S, Wortman J, Smith R, Hannick L, Maiti R, Ronning C, Rusch D, Town C: Improving the Arabidopsis genome annotation using maximal transcript alignment assemblies. Nucleic Acids Research. 2003, 31 (19): 5654-5666. 10.1093/nar/gkg770.
Attwood T, Croning M, Flower D, Lewis A, Mabey J, Scordis P, Selley J, Wright W: PRINTS-S: the database formerly known as PRINTS. Nucleic Acids Research. 2000, 28 (1): 225-227. 10.1093/nar/28.1.225.
Corpet F, Gouzy J, Kahn D: Recent improvements of the ProDom database of protein domain families. Nucleic Acids Research. 1999, 27 (1): 263-267. 10.1093/nar/27.1.263.
Haft D, Loftus B, Richardson D, Yang F, Eisen J, Paulsen I, White O: TIGRFAMs: a protein family resource for the functional identification of proteins. Nucleic Acids Research. 2001, 29 (1): 41-43. 10.1093/nar/29.1.41.
Wu C, Huang H, Yeh L, Barker W: Protein family classification and functional annotation. Computational Biology and Chemistry. 2003, 27 (1): 37-47. 10.1016/S1476-9271(02)00098-1.
Gough J, Karplus K, Hughey R, Chothia C: Assignment of homology to genome sequences using a library of hidden Markov models that represent all proteins of known structure. Journal of Molecular Biology. 2001, 313 (4): 903-919. 10.1006/jmbi.2001.5080.
Buchan D, Shepherd A, Lee D, Pearl F, Rison S, Thornton J, Orengo C: Gene3D: structural assignment for whole genes and genomes using the CATH domain structure database. Genome Research. 2002, 12 (3): 503-514. 10.1101/gr.213802.
Mi H, Lazareva-Ulitsky B, Loo R, Kejariwal A, Vandergriff J, Rabkin S, Guo N, Muruganujan A, Doremieux O, Campbell M: The PANTHER database of protein families, subfamilies, functions and pathways. Nucleic Acids Research. 2005, D284-D288. 33 Database
Lima T, Auchincloss A, Coudert E, Keller G, Michoud K, Rivoire C, Bulliard V, de Castro E, Lachaize C, Baratin D: HAMAP: a database of completely sequenced microbial proteome sets and manually curated microbial protein families in UniProtKB/Swiss-Prot. Nucleic Acids Research. 2009, D471-D478. 10.1093/nar/gkn661. 37 Database
Riano-Pachon D, Ruzicic S, Dreyer I, Mueller-Roeber B: PlnTFDB: an integrative plant transcription factor database. BMC Bioinformatics. 2007, 8 (1): 42-10.1186/1471-2105-8-42.
Kakar K, Wandrey M, Czechowski T, Gaertner T, Scheible W, Stitt M, Torres-Jerez I, Xiao Y, Redman J, Wu H: A community resource for high-throughput quantitative RT-PCR analysis of transcription factor gene expression in Medicago truncatula. Plant Methods. 2008, 4 (1): 18-10.1186/1746-4811-4-18.
Edgar R: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Research. 2004, 32 (5): 1792-1797. 10.1093/nar/gkh340.
Cheng J: A multi-template combination algorithm for protein comparative modeling. BMC Structural Biology. 2008, 8 (1): 18-10.1186/1472-6807-8-18.
CASP8 Group Performance. [http://www.predictioncenter.org/casp8/results.cgi]
Zemla A: LGA: a method for finding 3D similarities in protein structures. Nucleic Acids Research. 2003, 31 (13): 3370-3374. 10.1093/nar/gkg571.
Jones S, Shanahan H, Berman H, Thornton J: Using electrostatic potentials to predict DNA-binding sites on DNA-binding proteins. Nucleic Acids Research. 2003, 31 (24): 7189-7198. 10.1093/nar/gkg922.
Jmol: an open-source Java viewer for chemical structures in 3D. [http://jmol.sourceforge.net/]
Alexandrov N, Shindyalov I: PDP: protein domain parser. 2003, Oxford Univ Press, 19: 429-430.
Cheng J: DOMAC: an accurate, hybrid protein domain prediction server. Nucleic Acids Research. 2007, W354-W356. 10.1093/nar/gkm390. 35 Web Server
Crooks G, Hon G, Chandonia J, Brenner S: WebLogo: a sequence logo generator. Genome Research. 2004, 14 (6): 1188-1190. 10.1101/gr.849004.
Benson D, Boguski M, Lipman D, Ostell J, Ouellette B, Rapp B, Wheeler D: GenBank. Nucleic Acids Research. 1999, 27 (1): 12-17. 10.1093/nar/27.1.12.
Libault M, Joshi T, Takahashi K, Hurley-Sommer A, Puricelli K, Blake S, Xu D, Nguyen H, Stacey G: Large-scale analysis of putative soybean regulatory gene expression identifies a Myb gene involved in soybean nodule development. Plant Physiology. 2009, 151: 1207-1220. 10.1104/pp.109.144030.
Schlueter J, Lin J, Schlueter S, Vasylenko-Sanders I, Deshpande S, Yi J, O'Bleness M, Roe B, Nelson R, Scheffler B: Gene duplication and paleopolyploidy in soybean and the implications for whole genome sequencing. BMC genomics. 2007, 8 (1): 330-10.1186/1471-2164-8-330.
Schlueter J, Dixon P, Granger C, Grant D, Clark L, Doyle J, Shoemaker R: Mining EST databases to resolve evolutionary events in major crop species. Genome. 2004, 47 (5): 868-876. 10.1139/g04-047.
Merchant S, Prochnik S, Vallon O, Harris E, Karpowicz S, Witman G, Terry A, Salamov A, Fritz-Laylin L, Marechal-Drouard L: The chlamydomonas genome reveals the evolution of key animal and plant functions. Science. 2007, 318 (5848): 245-250. 10.1126/science.1143609.
Rensing S, Lang D, Zimmer A, Terry A, Salamov A, Shapiro H, Nishiyama T, Perroud P, Lindquist E, Kamisugi Y: The physcomitrella genome reveals evolutionary insights into the conquest of land by plants. Science. 2008, 319 (5859): 64-69. 10.1126/science.1150646.
Scofield S, Murray J: KNOX gene function in plant stem cell niches. Plant Molecular Biology. 2006, 60 (6): 929-946. 10.1007/s11103-005-4478-y.
Libault M, Joshi T, Benedito V, Xu D, Udvardi M, Stacey G: Legume transcription factor genes; what makes legumes so special?. Plant Physiology. 2009
Yan C, Terribilini M, Wu F, Jernigan R, Dobbs D, Honavar V: Predicting DNA-binding sites of proteins from amino acid sequence. BMC Bioinformatics. 2006, 7 (1): 262-10.1186/1471-2105-7-262.
Wang L, Brown S: BindN: a web-based tool for efficient prediction of DNA and RNA binding sites in amino acid sequences. Nucleic Acids Research. 2006, W243-W248. 10.1093/nar/gkl298. 34 Web Server
Hwang S, Gou Z, Kuznetsov I: DP-Bind: a web server for sequence-based prediction of DNA-binding residues in DNA-binding proteins. Bioinformatics. 2007, 23 (5): 634-636. 10.1093/bioinformatics/btl672.
Yamasaki K, Kigawa T, Inoue M, Tateno M, Yamasaki T, Yabuki T, Aoki M, Seki E, Matsuda T, Tomo Y: Solution structure of the B3 DNA binding domain of the Arabidopsis cold-responsive transcription factor RAV1. The Plant Cell Online. 2004, 16 (12): 3448-3459. 10.1105/tpc.104.026112.
Kaplan W, Littlejohn T: Swiss-PDB viewer (deep view). Briefings in Bioinformatics. 2001, 2 (2): 195-197. 10.1093/bib/2.2.195.
The PyMOL molecular graphics system. [http://pymol.sourceforge.net/]
This work was partially funded by a University of Missouri (MU) start-up grant, a MU research board grant, and a MU research council grant to JC. The worked conducted by ML and GS was funded by a grant from the National Science Foundation (Plant Genome Program, #DBI-0421620). Other funding includes grants from the Missouri Soybean Merchandising Council to HN, DX, and GS.
JC conceived the idea of automating the process of protein annotation and database construction, and coordinated the project. JC and ZW designed the database and formulated the content of the annotations. JC developed or executed the tools that generate structures, templates, domains, DNA binding sites, and multiple sequence alignments. ZW implemented the computer system that automatically gathers data, constructs database, builds web site, and links to external databases and data sources. ZW, TJ, ML and JC performed the soybean TF family classification. ZW, JC, ML, and TJ wrote the first draft of the manuscript. ML analyzed and compared the distributions of transcription factors in multiple species. TJ and ML identified the list of transcription factors. GS, DX, HN, and BV contributed some annotation ideas. JC, GS, DX, and HN provided financial support for the project. All edited and approved the manuscript.
Electronic supplementary material
Additional file 1: Figure S1 The SoyDB web page showing a list of transcription factor families. TF families are shown with their family ID, family name, and number of sequences within the family. Click on the family ID can further view the detailed information about the family as shown in Figure 4, and click on the number of sequences can open the webpage showing all the transcription factors within the family, as shown in Figure 5. (PDF 2 MB)
Additional file 2: Figure S2 The PSI-BLAST search web page. Users can paste or type in a query amino acid sequence and specify PSI-BLAST parameters on the web page. Click on the "Run" button will execute PSI-BLAST. (PDF 1001 KB)
Additional file 3: Figure S3 The result web page of PSI-BLAST search. PSI-BLAST result page shows the hit TF sequence ID, and the PSI-BLAST score and E-value. The hits are listed in a decreasing order of the PSI-BLAST score. Click on the sequence ID can open the web page showing detailed TF information, as shown in Figure 3. (PDF 2 MB)
Additional file 4: Figure S4 The HMM family classification web page. Users can paste or type in a query amino acid sequence. Click on the "Predict" button will execute family classification by HMM. (PDF 1 MB)
Authors’ original submitted files for images
About this article
Cite this article
Wang, Z., Libault, M., Joshi, T. et al. SoyDB: a knowledge database of soybean transcription factors. BMC Plant Biol 10, 14 (2010). https://doi.org/10.1186/1471-2229-10-14
- Hide Markov Model
- Protein Data Bank
- Transcription Factor Gene
- Transcription Factor Family
- Soybean Genome