Genomic Resources of Magnaporthe oryzae (GROMO): A comprehensive and integrated database on rice blast fungus
© Thakur et al; licensee BioMed Central Ltd. 2009
Received: 05 November 2008
Accepted: 15 July 2009
Published: 15 July 2009
Magnaporthe oryzae, rice blast fungus, is the most devastating pathogen of rice. It has emerged as a model phytopathogen for the study of host-pathogen interactions. A large body of data has been generated on different aspects of biology of this fungus and on host-pathogen interactions. However, most of the data is scattered and is not available as a single resource for researchers in this field.
Genomic Resources of Magnaporthe oyzae (GROMO), is a specialized, and comprehensive database for rice blast fungus, integrating information from several resources. GROMO contains information on genomic sequence, mutants available, gene expression, localization of proteins obtained from a variety of repositories, as primary data. In addition, prediction of domains, pathways, protein-protein interactions, sumolyation sites and biochemical properties that were obtained after computational analysis of protein sequences have also been included as derived data. This database has an intuitive user interface that shall prompt the user to explore various possible information resources available on a given gene or a protein, from a single source.
Currently, information on M. oryzae is available from different resources like BROAD MIT Magnaporthe database, Agrobacterium tumefaciens-mediated transformation (ATMT) M. oryzae database, Magnaporthe grisea – Oryza sativa (MGOS) and Massive Parallel Signature Sequencing (MPSS) databases. In the GROMO project, an effort has been made to integrate information from all these databases, derive some new data based on the available information analyzed by relevant programs and make more insightful predictions to better understand the biology of M. oryzae. The database is currently available at: http://gromo.msubiotech.ac.in/
M. oryzae, previously called M. grisea, is the most devastating fungal pathogen of rice, accounting for more than 10 million tons of yield loss every year . This fungus is a haploid filamentous ascomycete (class Pyrenomycetes). Members of M. oryzae species complex are also reported to cause disease on many other economically important crops such as barley, wheat, and millet . The rice blast fungus invades rice plants in a manner typical of many foliar pathogens by producing specialized infection structures called 'appressoria'. In nature, rice blast fungus attacks all above-ground parts of rice plants, and seedlings can be killed during epidemics. Under laboratory conditions, root infection of wheat and rice seedlings by M. oryzae has also been reported [3, 4].
Rice blast pathosystem has emerged as a model to study plant-pathogen interactions because genome sequence information on both the host  as well as the pathogen  is available and they are both amenable to genetic manipulation. M. oryzae was the first plant fungal pathogen to be sequenced. Analysis of M. oryzae genome is providing valuable insights on fungal pathogenesis .
M. oryzae genome is predicted to contain 12,841  and 11,074 (as per BROAD MIT database; version 6) genes , but the function of more than 70% of genes is still unknown. A few avirulence/pathogenicity genes have been characterized recently . One of the challenges in post-genomic era is to identify all transcribed regions and experimentally assign gene functions to M. oryzae genome. Although computational programs have played an important role in genome annotation, experimental evidence is needed to validate predicted functions. Availability of complete genome sequence of M. oryzae provides an opportunity to design insightful experiments so as to understand molecular mechanisms of pathogenesis. Expression of various genes has been analyzed in mycelium and appressorium using MPSS, robust-long serial analysis of gene expression (RL-SAGE) and oligoarray methods .
We have developed a comprehensive database, GROMO, integrating available information on genome, transcriptome, and proteome, thus, providing a useful resource for studies on functional genomics in the rice blast fungus. Data on this fungus were obtained from various resources like National Center for Biotechnology Information (NCBI) , BROAD MIT , MGOS , Kyoto Encyclopedia of Genes and Genomes (KEGG) , DNA Data Bank of Japan (DDBJ) , M. oryzae ATMT database , etc. Programs written using Perl language were used to analyze the data obtained from various resources and extract necessary information. The programs were also used to establish a connection with various databases for accessing the data using web services provided by the databases. All extracted information from the analyzed data was compiled in a back-end database created using MySQL for construction of GROMO. Some of the information incorporated in GROMO was also generated by correlating the extracted information from different resources with one another. Thus, a comprehensive resource on M. oryzae was built on protein domains, pathways, protein localization, presence of sumoylation sites, expression data, protein interacting partners, availability of mutants, etc. This database not only provides researchers an opportunity to extract detailed biological information on a specific gene or protein from a single resource but also prompts the researcher to explore new territories in fungal genomics.
Construction and content
Sequence information on M. oryzae proteins as available on BROAD, MIT database  was downloaded. Sequences from this database have unique locus tags which were used during the analysis for distinguishing sequences from each other. A large number of sequences obtained were described as either hypothetical or predicted. Naming convention for the sequences has been described at BROAD, MIT website.
Information on various M. oryzae mutants was obtained from MGOS database and ATMT database for M. oryzae which contains information on 21,070 hygromycin-resistant M. oryzae (strain KJ201) mutants generated through a large-scale insertional mutagenesis using ATMT . Information on 321 and 299 locus tags and their phenotypes was available for mutants from ATMT and MGOS respectively.
MGOS is a web-based database that contains expression profile data on rice-blast fungus interactions. It provides information on gene expression for both rice as well as M. oryzae during compatible interaction. Expression profiles obtained using techniques like Expressed Sequence Tags (ESTs) and Serial Analysis of Gene Expression (SAGE) on M. oryzae were downloaded and incorporated in our MySQL database file as a part of the back-end database. MPSS expression profile data on M. oryzae genes was obtained from MPSS database . Total number of genes for which EST contigs were found was 2,629. MPSS tags of 17-bp and 20-bp signatures were found in 3,911 genes and 5,095 genes were found to have SAGE tags.
Results of domain analysis performed using SMART and Pfam
No. of Genes
Putative pathways were predicted for M. oryzae protein sequences using KEGG Pathway database. KEGG PATHWAY is a collection of manually drawn pathway maps representing information on the molecular interaction and metabolic networks. In KEGG, 74,729 pathways are generated from 361 reference pathways. KEGG Application Programming Interface (API) was used for retrieving information from the database for analysis of putative pathways in M. oryzae. M. oryzae proteins were compared to Swiss-Prot database [23, 24] using BlastP  with an E-value cut-off of 10-3. Top 30 hits from blast result of each protein sequence were selected and checked in KEGG for prediction of biological pathways in which they were involved. Each query protein sequence from M. oryzae was assigned probable pathways on the basis of pathways obtained from KEGG database for their homologous protein sequences in other species (See Additional file 1: Schematic representation of pathway prediction in M. oryzae ). Approximately 219 different pathways were predicted for 2,477 proteins in M. oryzae.
Probable protein-protein interactions in M. oryzae were predicted using interactome information for Saccharomyces cerevisiae from Center of Cancer Systems Biology (CCSB) Yeast Interactome Database at Dana-Farber Cancer Institute and Harvard Medical School . Interaction data from the database was obtained as downloadable tab-delimited files. Protein interactions in M. oryzae were predicted on the basis of ortholog interaction data in S. cerevisiae. KEGG database was scanned to find orthologs of M. oryzae in S. cerevisiae and vice versa with an E-value cut off of 10-3. A total of 2,083 M. oryzae orthologs in S. cerevisiae that resulted from this study were subjected to scan against interactome data from CCSB Yeast Interactome database for prediction on protein-protein interactions (See Additional file 2: Schematic representation showing prediction of protein-protein interactions in M. oryzae).
Putative sumoylation sites in M. oryzae proteins were predicted using SUMOsp 2.0 software for sumoylation site prediction by Cuckoo workgroup . The non-redundant training data in software contained 279 sumoylation sites from 166 distinct proteins. SUMOsp 2.0 predicted 6,000 sumoylation sites for 4,494 protein sequences in M. oryzae at a high cut-off value.
Biochemical properties of the protein sequences were calculated using Pepstats program  from European Molecular Biology Open Software Suite (EMBOSS) package . EMBOSS is a free Open Source software analysis package specially developed for the needs of the molecular biology user community. Pepstats was programmatically linked and used to predict biochemical properties of M. oryzae proteins. Pepstats calculated molecular weight, isoelectric point, charge, size of protein, extinction coefficient and average residue weight for all the proteins.
The GROMO database is currently available at: http://gromo.msubiotech.ac.in/
Architecture and Design of GROMO
The architecture and design of GROMO (See Additional file 3: Overview of GROMO architecture and Design) consists of four tiers. T1: User Interface developed using HTML, T2: Apache web server  and programs for sequence analysis and information retrieval, T3: MySQL database  storing analysis data and T4: Perl CGI scripts  for retrieving and displaying analysis results for selected gene(s). All scripts and program in T2 and T4 are accessible from the analysis information page generated for each gene.
Tier 1: User Interface
User interface provides the user access to GROMO using various input queries and provide links to additional information pages which guide the user during browsing of GROMO. The query inputs from user interface are sent to programs in layer T2 via post method.
Tier 2: Programs for sequence analysis and information retrieval
Apache web server receives query request from user interface and sends it to Perl CGI script in T2 for retrieving locus tag and gene description information from MySQL database. Program in T2 also use BLAST program obtained from NCBI ftp site for sequence based analysis and parse BLAST result to represent necessary information on browser.
Tier 3: Database Schema
The Relational Database Management System MySQL  was used to store data integrated in GROMO. DBD::mysql (Database Driver) and DBI (Database Interface) module are used in Perl CGI scripts for accessing data from MySQL database. Database schema was divided into different tables depending on the type of data incorporated. The description of different tables of schema is given in the additional file 4: Description of information contained in different tables of relational database.
Tier 4: Scripts retrieving and displaying analysis result for selected gene
Perl CGI scripts in this layer retrieve analyzed information from MySQL database for selected locus tag from result page and display it on browser. "Export" option for downloading analyzed information in the form of text file is provided on summary section for selected gene.
Web Interface Access
Data stored in GROMO can be accessed through web interfaces which have been generated using HTML. It gives entry point to explore the data stored in GROMO by following two different kinds of links provided: (i) Search using Keyword (ii) Link to BLAST search.
Search using Keyword
BLAST-based search allows user to browse GROMO using sequence in FASTA format. With a query sequence, BLAST generates a table containing best hits in M. oryzae which are organized according to percent identity between the sequences.
Representation of analysis results
Other web interfaces
Other web interfaces includes "Help" section which provides information about the inputs accepted and output generated. "Link" section provides links to external web resources and groups working on M. oryzae. "Feedback" form provides user with an option to post their comments and queries about the database. "GROMO" section gives overall diagrammatic representation of complete database and its contents. User can submit any valid new information available about genes and proteins of M. oryzae using the format given on home page of the database.
The queries provided by GROMO are focused on retrieving details available from various databases along with queried information for particular a gene or protein in M. oryzae. Currently, information about this fungus is dispersed in different resources at different locations. Among the existing databases, (i) BROAD MIT provides sequence information on M. oryzae genes and proteins with some information on protein domains; (ii) M. oryzae ATMT database provides information on some mutants and (iii) Information on gene expression during host-pathogen interactions is available from MGOS and MPSS databases.
In the GROMO project, an effort has been made to bring together information on M. oryzae from various resources and develop a comprehensive database which consists of experimental and computational data. GROMO database provides researcher information not only on gene and protein sequences but also on possible domains present in a protein, predicted pathways, probable interacting partners, sub-cellular localization, protein sumoylation sites and even biochemical properties of any protein. In addition to a common blast search, GROMO provides the user the possibility of keyword search using the options like locus tag, domain name, pathway, localization, SignalP prediction, SAGE tag, MPSS tag, EST ID and mutants. Moreover, some of the experimental data obtained from external resources are represented in more interpretable form, thus, providing researchers with a better understanding about the fungus in order to help design critical experiments to gain deep insights into biology of fungi in general and rice blast in particular. In order to incorporate newer findings, the database will be updated every 6 months. The directly submitted data will be cross-checked against all existing databases and will be uploaded to GROMO with all the possible predicted features for a specific protein.
GROMO is a composite and unique resource database on rice blast fungus. In this project, an effort has been made to integrate information from various public databases, derive some new data based on available information resources and make more intuitive predictions to better understand the biology of this fungus. Many scientists are working actively on functional characterization of genes of the blast fungus, expression profiling during host-pathogen interactions and also in the field of proteomics. Updating the current information on blast functional genomics using newer findings coupled with the use of bioinformatics tools will be an area of future interest. Magnaporthe research community is encouraged to submit its research findings periodically so that this unique resource database can be kept up to date.
Availability and requirements
Project Name: GROMO – Genomic Resources of Magnaporthe oryzae
Project homepage: The database is currently available at http://gromo.msubiotech.ac.in/
Operating system(s): Platform independent (web server)
Programming language(s): HTML, Perl, CGI
License: Free for academics, Authorization is needed for commercial use
(Please contact the corresponding author for more details)
We are thankful to Mr. Jay Kumar Shah for designing the background of the web server and his help in analyzing the pathway-related data. Thanks to Mr. A. Varghese and Mr. N. D. Bhatt for their help with the computing facilities. This work was supported by the grants from Department of Biotechnology (DBT), Ministry of Science and Technology, Govt. of India.
- Xu J-R, Zhao X, Dean RA: From genes to genomes; a new paradigm for studying fungal pathogenesis in Magnaporthe oryzae. Adv Genet. 2007, 57: 175-218.View ArticlePubMedGoogle Scholar
- Valent B, Chumley FG: Molecular genetic analysis of the. rice blast fungus Magnaporthe grisea. Annu Rev Phytopathol. 1991, 29: 443-467. 10.1146/annurev.py.29.090191.002303.View ArticlePubMedGoogle Scholar
- Dufresne M, Osbourn AE: Definition of tissue-specific and general requirements for plant infection in a phytopathogenic fungus. Mol Plant Microbe Interact. 2001, 14 (3): 300-307. 10.1094/MPMI.2001.14.3.300.View ArticlePubMedGoogle Scholar
- Sesma A, Osbourn AE: The rice leaf blast pathogen undergoes developmental processes typical of root-infecting fungi. Nature. 2004, 431 (7008): 582-586. 10.1038/nature02880.View ArticlePubMedGoogle Scholar
- International Rice Genome Sequencing Project: The map-based sequence of the rice genome. Nature. 2005, 436 (7052): 793-800. 10.1038/nature03895.View ArticleGoogle Scholar
- Dean RA, Talbot NJ, Ebbole DJ, Farman ML, Mitchell TK, Orbach MJ, Thon M, Kulkarni R, Xu J-R, Pan H, et al: The genome sequence of the rice blast fungus Magnaporthe grisea. Nature. 2005, 434 (7036): 980-986. 10.1038/nature03449.View ArticlePubMedGoogle Scholar
- Park J, Park B, Jung K, Jang S, Yu K, Choi J, Kong S, Park J, Kim S, Kim H: CFGP: a web-based comparative fungal genomics platform. Nucl Acid Res. 2008, D562-D571. 36 DatabaseGoogle Scholar
- BROAD MIT Magnaporthe database. [http://www.broad.mit.edu/annotation/genome/magnaporthe_grisea/Home.html]
- Talbot NJ: On the trail of a cereal killer: Exploring the biology of Magnaporthe grisea. Annu Rev Microbiol. 2003, 57: 177-202. 10.1146/annurev.micro.57.030502.090957.View ArticlePubMedGoogle Scholar
- Gowda M, Venu RC, Raghupathy MB, Nobuta K, Li H, Wing R, Stahlberg E, Couglan S, Haudenschild CD, Dean R, et al: Deep and comparative analysis of the mycelium and appressorium transcriptomes of Magnaporthe grisea using MPSS, RL-SAGE and oligoarray methods. BMC Genomics. 2006, 7: 310-10.1186/1471-2164-7-310.PubMed CentralView ArticlePubMedGoogle Scholar
- National Center for Biotechnology Information. [http://www.ncbi.nlm.nih.gov]
- Soderlund C, Haller K, Pampanwar V, Ebbole D, Farman M, Orbach MJ, Wang GL, Wing R, Xu JR, Brown D, et al: MGOS: A resource for studying Magnaporthe grisea and Oryza sativa interactions. Mol Plant Microbe Interact. 2006, 19 (10): 1055-1061. 10.1094/MPMI-19-1055.View ArticlePubMedGoogle Scholar
- Ogata H, Goto S, Sato K, Fujibuchi W, Bono H, Kanehisa M: KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucl Acid Res. 1999, 27 (1): 29-34. 10.1093/nar/27.1.29.View ArticleGoogle Scholar
- DNA Data Bank of Japan. [http://www.ddbj.nig.ac.jp/]
- Jeon J, Park S-Y, Chi M-H, Choi J, Park J, Rho H-S, Kim S, Goh J, Yoo S, Choi J: Genome-wide functional analysis of pathogenicity genes in the rice blast fungus. Nat Genet. 39 (4): 561-565. 10.1038/ng2002.Google Scholar
- Hedeler C, Wong HM, Cornell MJ, Alam I, Soanes DM, Rattray M, Hubbard SJ, Talbot NJ, Oliver SG, Paton NW: e-Fungi: a data resource for comparative analysis of fungal genomes. BMC Genomics. 426-8 DatabaseGoogle Scholar
- Nakai K, Horton P: PSORT: a program for detecting sorting signals in proteins and predicting their subcellular localization. Trend Biochem Sci. 1999, 24 (1): 34-35. 10.1016/S0968-0004(98)01336-X.View ArticlePubMedGoogle Scholar
- Horton P, Park KJ, Obayashi T, Fujita N, Harada H, Adams-Collier CJ, Nakai K: WoLF PSORT: Protein Localization Predictor. Nucleic Acids Res. 2007, 35 (Web Server): W585-W587. 10.1093/nar/gkm259.PubMed CentralView ArticlePubMedGoogle Scholar
- Nielsen H, Engelbrecht J, Brunak S, Heijne G: Identification of prokaryotic and eukaryotic signal peptides and prediction of their cleavage sites. Protein Eng. 1997, 10 (1): 1-6. 10.1093/protein/10.1.1.View ArticlePubMedGoogle Scholar
- Schultz J, Copley RR, Doerks T, Ponting CP, Bork P: A Web-based tool for the study of genetically mobile domains. Nucleic Acids Res. 2000, 28 (1): 231-234. 10.1093/nar/28.1.231.PubMed CentralView ArticlePubMedGoogle Scholar
- Finn RD, Mistry J, Schuster-Böckler B, Griffiths-Jones S, Hollich V, Lassmann T, Moxon S, Marshall M, Khanna A, Durbin R: Pfam: clans, web tools and services. Nucl Acid Res. 2006, D247-D251. 10.1093/nar/gkj149. 34 DatabaseGoogle Scholar
- Servant F, Bru C, Carrere S, Courcelle E, Gouzy J, Peyruc D, Kahn D: ProDom: Automated clustering of homologous domains. Brief Bioinform. 2002, 3 (3): 246-251. 10.1093/bib/3.3.246.View ArticlePubMedGoogle Scholar
- Bairoch A, Apweiler R: The SWISS-PROT protein sequence data bank and its supplement TrEMBL. Nucleic Acids Res. 1997, 25 (1): 31-36. 10.1093/nar/25.1.31.PubMed CentralView ArticlePubMedGoogle Scholar
- Bairoch A, Apweiler R: The SWISS-PROT protein sequence data bank and its supplement TrEMBL in 1998. Nucleic Acids Res. 1998, 26 (1): 38-42. 10.1093/nar/26.1.38.PubMed CentralView ArticlePubMedGoogle Scholar
- 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. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.PubMed CentralView ArticlePubMedGoogle Scholar
- Yu H, Braun P, Yildirim MA, Lemmens I, Venkatesan K, Sahalie J, Hirozane-Kishikawa T, Gebreab F, Li N, et al: High-Quality Binary Protein Interaction Map of the Yeast Interactome Network. Science. 2008, 322 (5898): 104-110. 10.1126/science.1158684.PubMed CentralView ArticlePubMedGoogle Scholar
- Xue Y, Zhou F, Fu C, Xu Y, Yao X: SUMOsp: a web server for sumoylation site prediction. Nucl Acid Res. 2006, W254-W257. 10.1093/nar/gkl207. 34 Web ServerGoogle Scholar
- Rice P, Longden I, Bleasby A: EMBOSS: The European Molecular Biology Open Software Suite. Trends Genet. 2000, 16 (6): 276-277. 10.1016/S0168-9525(00)02024-2.View ArticlePubMedGoogle Scholar
- The European Molecular Biology Open Software Suite. [http://emboss.sourceforge.net/]
- Apache web server. [http://httpd.apache.org/]
- MySQL Database Sever. [http://www.mysql.com/]
- Perl CGI scripts. [http://www.activestate.com/Products/activeperl/index.mhtml]
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.