Gene Ontology annotation of the rice blast fungus, Magnaporthe oryzae
© Meng et al; licensee BioMed Central Ltd. 2009
Published: 19 February 2009
Magnaporthe oryzae, the causal agent of blast disease of rice, is the most destructive disease of rice worldwide. The genome of this fungal pathogen has been sequenced and an automated annotation has recently been updated to Version 6 http://www.broad.mit.edu/annotation/genome/magnaporthe_grisea/MultiDownloads.html. However, a comprehensive manual curation remains to be performed. Gene Ontology (GO) annotation is a valuable means of assigning functional information using standardized vocabulary. We report an overview of the GO annotation for Version 5 of M. oryzae genome assembly.
A similarity-based (i.e., computational) GO annotation with manual review was conducted, which was then integrated with a literature-based GO annotation with computational assistance. For similarity-based GO annotation a stringent reciprocal best hits method was used to identify similarity between predicted proteins of M. oryzae and GO proteins from multiple organisms with published associations to GO terms. Significant alignment pairs were manually reviewed. Functional assignments were further cross-validated with manually reviewed data, conserved domains, or data determined by wet lab experiments. Additionally, biological appropriateness of the functional assignments was manually checked.
In total, 6,286 proteins received GO term assignment via the homology-based annotation, including 2,870 hypothetical proteins. Literature-based experimental evidence, such as microarray, MPSS, T-DNA insertion mutation, or gene knockout mutation, resulted in 2,810 proteins being annotated with GO terms. Of these, 1,673 proteins were annotated with new terms developed for Plant-Associated Microbe Gene Ontology (PAMGO). In addition, 67 experiment-determined secreted proteins were annotated with PAMGO terms. Integration of the two data sets resulted in 7,412 proteins (57%) being annotated with 1,957 distinct and specific GO terms. Unannotated proteins were assigned to the 3 root terms. The Version 5 GO annotation is publically queryable via the GO site http://amigo.geneontology.org/cgi-bin/amigo/go.cgi. Additionally, the genome of M. oryzae is constantly being refined and updated as new information is incorporated. For the latest GO annotation of Version 6 genome, please visit our website http://scotland.fgl.ncsu.edu/smeng/GoAnnotationMagnaporthegrisea.html. The preliminary GO annotation of Version 6 genome is placed at a local MySql database that is publically queryable via a user-friendly interface Adhoc Query System.
Our analysis provides comprehensive and robust GO annotations of the M. oryzae genome assemblies that will be solid foundations for further functional interrogation of M. oryzae.
Magnaporthe oryzae, the rice blast fungus, infects rice and other agriculturally important cereals, such as wheat, rye and barley. The pathogen is found throughout the world and each year is estimated to destroy enough rice to feed more than 60 million people . A comprehensive understanding of the genetic makeup of the rice blast fungus is crucial in efforts to understand the biology and develop effective disease management strategies of this destructive pathogen.
The rice blast fungus has been the focus of intense investigation and a number of genomic resources have been generated. These include a genome sequence , genome-wide expression from microarray  and massive parallel signature sequencing (MPSS) , as well as large bank of T-DNA insertion mutants [5, 6]. In addition, numerous genes have been functionally characterized by targeted knockout [7–18]. While these resources are of tremendous utility, much of the genome remains unexplored. Till now, only automated annotations of the predicted genes have been available. In order to maximize the utility of the genome sequence, we have developed manual curations of the predicted genes.
Providing functionality through careful and comprehensive application of a standardized vocabulary, such as the Gene Ontology (GO) requires manual curation. The GO has evolved into a reliable and rapid means of assigning functional information [19–22]. There are two types of GO annotations. One is referred to as similarity-based GO annotation, and the other is literature-based GO annotation. Similarity-based GO annotation applies computational approaches to match characterized gene products and their associated GO terms to gene products under study. Orthology-based GO annotation is a more stringent application of similarity-based GO annotation. Literature-based GO annotation involves reviewing published work and then manually assigning GO terms to characterized gene products. Currently, similarity-based GO annotation predominates since it is rapid and relatively inexpensive [21, 23]. On the other hand, although literature-based annotation is time consuming, it is considered more reliable and provides a mechanism to assign previously unassigned GO terms or newly developed GO terms to proteins. Here, we present an overview of the strategy and results obtained from the integration of both approaches to assign GO terms to Versions 5 of M. oryzae genome.
M. oryzae genome sequence
This paper summarizes manual annotation of Version 5 of the M. oryzae genome sequence. At the time of submission of this manuscript, Version 6 of the genome assembly of M. oryzae was released by the Broad Institute. Version 6 will be annotated using the same methodology described here. A preliminary GO annotation of the Version 6 genome sequence, based on the Version 5 annotation, has been placed at our local MySQL database at http://scotland.fgl.ncsu.edu/smeng/GoAnnotationMagnaporthegrisea.html.
Sequence similarity-based GO annotation
Predicted proteins of Version 5 of the M. oryzae genome sequence were downloaded from the Broad Institute at http://www.broad.mit.edu/annotation/genome/magnaporthe_grisea/MultiDownloads.html. GO-annotated proteins were downloaded from the Gene Ontology (GO) database at http://www.Geneontology.org/GO.downloads.database.shtml. These GO-annotated proteins were from about 50 organisms with published association with GO terms. Only three of the 50 organisms are fungi. They are Candida albicans, Saccharomyces cerevisiae, and Schizosaccharomyces pombe. Other organisms are from bacteria, plants, or animals etc. Proteins of these non-fungal organisms were retained to increase the number of proteins with validated functions available for matching to M. oryzae.
Possible ortholog pairs between GO proteins and predicted proteins from M. oryzae genome sequence Version 5 were estimated by searching for reciprocal best hits using BLASTP (e-value < 10-3) .
Significant alignment pairs with 80% or better coverage of both query and subject sequences, 10-20 or less BLASTP E-value, and 40% or higher of amino acid identity (pid) were manually reviewed.
The functions of significantly matched GO proteins were manually cross- validated using data from wet lab experiments, and the NCBI Conserved Domain Database (CDD) .
If the functions suggested from different sources were consistent with each other, and with available M. oryzae data, the functions of the experimentally characterized, significantly matched GO proteins, were transferred to the M. oryzae proteins in our study, and given the evidence code ISS (Inferred from Sequence Similarity) [26, 27].
The information was recorded into a gene association file following the format standard at http://www.geneontology.org/GO.format.annotation.shtml.
Literature-based GO annotation
Literature at public databases such as PubMed [a database of biomedical literature citations and abstracts at the National Center for Biotechnology Information (NCBI)] were searched using key words, including alternative species names for the organism such as Magnaporthe grisea and Pyricularia oryzae.
Relevant published papers were read and genes or gene products and their functions were identified.
Where necessary, gene IDs and sequences at public databases, such as the NCBI protein database were identified.
Based on the functions identified in the paper(s), appropriate GO terms were found using AmiGO, the GO-supported tool for searching and browsing the Gene Ontology database.
Evidence codes were assigned following the guide at http://www.geneontology.org/GO.evidence.shtml.
Data were recorded into the gene association file manually or using custom PERL scripts for large gene sets with the same biological process.
Integration of the results from the two types of GO annotations
Similarity-based annotations were replaced with literature-based annotations, where redundant, using custom PERL scripts.
Custom PERL scripts were used to annotate each protein with GO terms from the three ontologies using the following protocol. Any protein not annotated with a GO term following similarity-based and literature-based GO annotations was annotated with the three root GO terms, GO:0005575 (Cellular Component), GO:0003674 (Molecular Function), and GO:0008150 (Biological Process). Additionally, if any protein was lacking annotation from any of the three GO categories, Cellular Component, Molecular Function, or Biological Process, the protein was annotated with the root GO terms of the missing GO categories.
Errors in the gene association file were checked using the script, filter-gene-association.pl, which was downloaded from the GO database at ftp://ftp.geneontology.org/pub/go/software/utilities/filter-gene-association.pl.
The gene association file for Version 5 of the M. oryzae genome sequence was uploaded to the GO database at http://www.geneontology.org/GO.current.annotations.shtml. Many protocols and scripts were created for generating and parsing the data. For example, a protocol and five scripts were developed to replace redundant similarity-based annotation with literature-based annotation. Furthermore, a protocol and eight scripts were developed to provide each gene with a GO term from the three ontologies. In addition, a PERL script to record many genes into the gene association file was developed. This script, with slight modification, easily recorded different types of data, such as microarray expression, MPSS, or T-DNA insertion mutation, etc., into the gene association file. These protocols and scripts are available upon request from the corresponding or the first author.
Computational GO annotation
A total of 67 secreted proteins of M. oryzae was experimentally demonstrated to be secreted through cloning into an overexpression vector and expressed in M. oryzae transformants (Ebbole and Dean, unpublished data). These 67 secreted proteins were annotated with a biological process term GO:0009306 ("protein secretion") and a cellular component term GO:0005576 (extracellular region). An evidence code IDA was assigned to annotations of these 67 proteins since function was determined through direct assay.
A total of 128 curated cytochrome P450's of M. oryzae were validated by comparison and analysis of gene location and structure, clustering of genes, and phylogenetic reconstruction . Different subsets of these proteins were annotated with different GO terms. For example, 75 of the 128 P450 proteins were annotated with the molecular function term GO:0005506 ("iron ion binding"), and 40 P450 proteins with the molecular function term GO:0016491 ("oxidoreductase activity"). An evidence code IGC was assigned to annotations of these P450 proteins since annotations were based on genomic context.
A total of 428 putative transcription factors of M. oryzae were validated by integrated computational analysis of whole genome microarray expression data, and matches to InterPro, pfam, and COG . Again, different subsets of the 428 proteins were annotated with different GO terms. For example, 36 proteins were annotated with GO:0005975 ("carbohydrate metabolic process"), and 12 proteins were annotated with GO:0006520 ("amino acid metabolic process"). An evidence code RCA was assigned to annotations of the 428 transcription factors since the annotations were based on reviewed computational analysis.
A total of 2,548 conserved domains from NCBI CDD were used as evidence for cross-checking putative functions, but no GO annotation was made based solely on identification of these domains.
In addition, the evidence code ISS was assigned to annotations of 216 M. oryzae proteins for the following reasons: 1) These proteins have significant similarity to experimentally-characterized homologs over the majority (at least 80%) of the full length sequences. 2) The pairwise alignments of good matches between the characterized proteins and the proteins of M. oryzae were manually reviewed. 3) Functional domains were conserved between the M. oryzae proteins and their homologs. 4) The GO assignments from the characterized match proteins to the M. oryzae proteins were manually determined to be biologically relevant.
The remaining 1,343 proteins with a reciprocal BLASTP best match of e-value > 10-20 and pid < 40% were assigned GO terms from their characterized matches, but the evidence codes were identified as IEA (Inferred from Electronic Annotation).
In sum, GO terms were assigned to 6,286 proteins of M. oryzae. Among the 6,286 proteins, 2,732 hypothetical proteins, 125 predicted proteins, and 14 unknown proteins were assigned functions.
Literature-based GO annotation
More than 400 research articles were read, and 71 genes with gene knockout mutations and with accession numbers and sequences deposited in public databases such as NCBI were manually annotated using GO terms, including newly developed Plant-Associated Microbe Gene Ontology (PAMGO) terms. Gene products were annotated with GO terms relevant to their biological functions. For example, 6 genes were annotated with GO:0000187 ("activation of MAPK activity"), 5 genes with GO:0075053 ("formation of symbiont penetration peg for entry into host"), 14 genes with GO:0044409 ("entry into host"), 8 genes with GO:0044412 ("growth or development of symbiont within host"), and 43 genes with GO:0009405 ("pathogenesis"). The evidence code IMP (inferred from Mutant Phenotype) was assigned to these annotations since gene-knockout mutants were generated in order to determine functions of these genes.
A total of 210 genes were annotated on the basis of published microarray studies . Again, gene products were annotated with GO terms, including PAMGO terms, relevant to their biological functions. For example, 67 genes were annotated with GO:0044271 ("nitrogen compound biosynthetic process"), 27 genes with GO:0075005 ("spore germination on or near host"), 26 genes with GO:0075035 ("maturation of appressorium on or near host"), and 114 genes with GO:0075016 ("appressorium formation on or near host"). The evidence code IEP (Inferred from expression Pattern) was assigned to these annotations on the basis that the genes were up-regulated by at least 10-fold in association with the particular biological process. A further 2,433 genes were annotated on the basis of published Massively Parallel Signature Sequencing ( MPSS) studies , including 1,041 genes annotated with GO:0043581 ("mycelium development"), and 1,392 genes annotated with GO:0075016 ("appressorium formation on or near host"). The evidence code IEP was also assigned to these annotations since the genes were up-regulated only during a certain biological process, such as mycelium formation, and the fold change was equal to or greater than 10.
On the basis of whole genome T-DNA insertion mutation data , 120 genes were annotated with relevant GO terms and PAMGO terms. For instance, 43 genes were annotated with GO:0030437 ("ascospore formation"), 14 genes with GO:0009847 ("spore germination"), 64 genes with GO:0075016 ("appressorium formation on or near host"), and 106 genes with GO:0009405 ("pathogenesis"). An evidence code IMP (inferred from mutant phenotype) was assigned to these annotations.
In total, 2,810 proteins were annotated based on experimental data from published peer-reviewed literature. Of these, 1,673 proteins were annotated with terms created by the PAMGO consortium to describe interactions between symbionts and their hosts.
Integration of results from the two types of GO annotations
Integration of the similarity-based and literature-based annotation resulted in 7,412 proteins being annotated with specific GO terms, covering more than 57% of the inferred proteome. The remaining 5,464 predicted proteins, not having high similarity to GO-annotated proteins, were annotated with three general GO terms. GO:0005575 (Cellular Component), GO:0003674 (Molecular Function), and GO:0008150 (Biological Process). Therefore, our GO annotation provides an annotation of the entire 12,832 proteins predicted in M. oryzae, and each protein being annotated with GO terms from the three GO categories.
The GO annotation of Version 5 of the genome sequence of Magnaporthe oryzae is available at the GO Consortium database http://www.geneontology.org/GO.current.annotations.shtml.
Here, we present a detailed protocol for integrating the results of similarity-based annotation with a literature-based annotation of the predicted proteome of Version 5 of the genome sequence of the rice blast fungus M. oryzae. Through careful manual inspection of these annotations, we are able to provide a reliable and robust GO annotation for more than half of the predicted gene products. Of 6,286 proteins receiving computational annotations, only 1,343 did not exceed our stringent match criteria upon manual review and so were assigned the evidence code IEA. It should be noted that annotations with the IEA evidence code are retained in the GO database for only one year, and then the GO Consortium will remove them from a gene association file. To be retained, IEA annotations must be manually reviewed in order to be assigned an upgraded evidence code such as ISS (Inferred from Sequence or Structural Similarity). Currently, there is no recognized standard to assign the ISS code. We recommend the following criteria for assigning the ISS code:
The functions of the proteins from which the annotation will be transferred must be experimentally characterized.
The similarity between the characterized proteins and the proteins under study must be significant. For example, we used ≥ 80% coverage of both query and subject sequences, ≤ 10-20 E-value, and ≥ 40% percentage of identity (pid) as cutoff criteria in our similarity-based GO annotation. Ideally, orthology should be established by phylogenetic analysis.
The pairwise alignment between the characterized proteins and the proteins under study should be manually reviewed and cross-validated with characterized or reviewed data of other resources such as functional domains, active sites, and sequence patterns etc.
Biological appropriateness of all assigned GO terms should be manually reviewed.
All authors read and approved the final manuscript. We thank Michelle Gwinn Giglio, Brett Tyler, and Candace Collmer for their comments and suggestions in annotating the genome of the rice blast fungus Magnaporthe grisea with GO terms, and Brett Tyler for editing of the manuscript. This work is a part of the PAMGO project, which is supported by National Research Initiative of the USDA Cooperative State Research, Education and Extension Service (grant number 2005-35600-16370), and the National Science Foundation (grant number EF-0523736).
This article has been published as part of BMC Microbiology Volume 9 Supplement 1, 2009: The PAMGO Consortium: Unifying Themes In Microbe-Host Associations Identified Through The Gene Ontology. The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2180/9?issue=S1.
- Brüssow H: The quest for food. 2007, Springer, New YorkGoogle Scholar
- Dean RA, Talbot NJ, Ebbole DJ, Farman ML, Mitchell TK, Orbach MJ, Thon M, Kulkarni R, Xu J-R, Pan H, Read ND, Lee Y-H, Carbone I, Brown D, Oh YY, Donofrio N, Jeong JS, Soanes DM, Djonovic S, Kolomiets E, Rehmeyer C, Li W, Harding M, Kim S, Lebrun M-H, Bohnert H, Coughlan S, Butler J, Calvo S, Ma L-J, Nicol R, Purcell S, Nusbaum C, Galagan JE, Birren BW: The genome sequence of the rice blast fungus Magnaporthe grisea. Nature. 2005, 434: 980-986. 10.1038/nature03449.View ArticlePubMedGoogle Scholar
- Oh YY, Donofrio N, Pan H, Coughlan S, Brown DE, Meng S, Mitchell T, Dean RA: Transcriptome analysis reveals new insight into appressorium formation and function in the rice blast fungus Magnaporthe oryzae. Genome Biol. 2008, 9 (5): R85-10.1186/gb-2008-9-5-r85.PubMed CentralView ArticlePubMedGoogle Scholar
- Gowda M, Venu RC, Raghupathy , Mohan B, Nobuta K, Li H, Wing R, Stahlberg E, Couglan S, Haudenschild , Christian D, Dean R, Nahm B-H, Meyers BC, Wang G-L: 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
- Jeon J, Park SY, Chi MH, Choi J, Park J, Rho HS, Kim S, Goh J, Yoo S, Choi J, Park JY, Yi M, Yang S, Kwon MJ, Han SS, Kim BR, Khang CH, Park B, Lim SE, Jung K, Kong S, Karunakaran M, Oh HS, Kim H, Kim S, Park J, Kang S, Choi WB, Kang S, Lee YH: Genome-wide functional analysis of pathogenicity genes in the rice blast fungus. Nat Genet. 2007, 39 (4): 561-565. 10.1038/ng2002.View ArticlePubMedGoogle Scholar
- Choi J, Park J, Jeon J, Chi MH, Goh J, Yoo SY, Park J, Jung K, Kim H, Park SY, Rho HS, Kim S, Kim BR, Han SS, Kang S, Lee YH: Genome-wide analysis of T-DNA integration into the chromosomes of Magnaporthe oryzae. Mol Microbiol. 2007, 66 (2): 371-382. 10.1111/j.1365-2958.2007.05918.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu S, Dean RA: G protein a subunit genes control growth, development, and pathogenicity of Magnaporthe grisea. Mol Plant-Micro Interact. 1997, 10 (9): 1075-1086. 10.1094/MPMI.19184.108.40.2065.View ArticleGoogle Scholar
- Choi W, Dean RA: The adenylate cyclase gene MACI of Magnaporthe grisea controls appressorium formation and other aspects of growth and development. Plant Cell. 1997, 9: 1973-1983. 10.1105/tpc.9.11.1973.PubMed CentralView ArticlePubMedGoogle Scholar
- Kulkarni RD, Dean RA: Identification of proteins that interact with two regulators of appressorium development, adenylate cyclase and cAMP-dependent protein kinase A, in the rice blast fungus Magnaporthe grisea. Mol Genet Genomics. 2004, 270: 497-508. 10.1007/s00438-003-0935-y.View ArticlePubMedGoogle Scholar
- Wilson RA, Jenkinson JM, Gibson RP, Littlechild JA, Wang Z-Y, Talbot NJ: Tps1 regulates the pentose phosphate pathway, nitrogen metabolism and fungal virulence. EMBO J. 2007, 26: 3673-3685. 10.1038/sj.emboj.7601795.PubMed CentralView ArticlePubMedGoogle Scholar
- Gilbert MJ, Thornton CR, Wakley GE, Talbot NJ: A P-type ATPase required for rice blast disease and induction of host resistance. Nature. 2006, 440: 535-539. 10.1038/nature04567.View ArticlePubMedGoogle Scholar
- Veneault-Fourrey C, Barooah M, Egan M, Wakley G, Talbot NJ: Autophagic fungal cell death is necessary for infection by the rice blast fungus. Science. 2006, 312: 580-583. 10.1126/science.1124550.View ArticlePubMedGoogle Scholar
- Li L, Ding S-L, Sharon A, Orbach M, Xu J-R: Mir1 is highly upregulated and localized to nuclei during infectious hyphal growth in the rice blast fungus. Mol Plant-Micro Interact. 2007, 20 (4): 448-458. 10.1094/MPMI-20-4-0448.View ArticleGoogle Scholar
- Zhao X, Xu J-R: A highly conserved MAPK-docking site in Mst7 is essential for Pmk1 activation in Magnaporthe grisea. Mol Microbiol. 2007, 63 (3): 881-894. 10.1111/j.1365-2958.2006.05548.x.PubMedGoogle Scholar
- Orbach MJ, Farrall L, Sweigard JA, Chumley FG, Valent B: A telomeric avirulence gene determines efficacy for the rice blast resistance gene Pi-ta. Plant Cell. 2000, 12 (11): 2019-32. 10.1105/tpc.12.11.2019.PubMed CentralView ArticlePubMedGoogle Scholar
- DeZwaan TM, Carroll AM, Valent B, Sweigard JA: Magnaporthe grisea pth11p is a novel plasma membrane protein that mediates appressorium differentiation in response to inductive substrate cues. Plant Cell. 1999, 11 (10): 2013-2030. 10.1105/tpc.11.10.2013.PubMed CentralView ArticlePubMedGoogle Scholar
- Vergne E, Ballini E, Marques S, Mammar BS, Droc G, Gaillard S, Bourot S, DeRose R, Tharreau D, Nottéghem J-L, Lebrun M-H, Morel J-B: Early and specific gene expression triggered by rice resistance gene Pi33 in response to infection by ACE1 avirulent blast fungus. New phytol. 2007, 174: 159-171. 10.1111/j.1469-8137.2007.01971.x.View ArticlePubMedGoogle Scholar
- Gupta A, Chattoo BB: A novel gene MGA1 is required for appressorium formation in Magnaporthe grisea. Fungal Genet Biol. 2007, 44: 1157-1169. 10.1016/j.fgb.2007.02.014.View ArticlePubMedGoogle Scholar
- The Gene Ontology Consortium: Gene Ontology: tool for the unification of biology. Nat Genet. 2000, 25: 25-29. 10.1038/75556.PubMed CentralView ArticleGoogle Scholar
- Gene Ontology Consortium: The Gene Ontology (GO) project in 2006. Nucleic Acids Res. 2006, 34: D322-D326. 10.1093/nar/gkj021.PubMed CentralView ArticleGoogle Scholar
- The Gene Ontology Consortium: The Gene Ontology project in 2008. Nucleic Acids Res. 2008, 36: D440-D444. 10.1093/nar/gkm883.PubMed CentralView ArticleGoogle Scholar
- Hong EL, Balakrishnan R, Dong Q, Christie KR, Park J, Binkley G, Costanzo MC, Dwight SS, Engel SR, Fisk DG, Hirschman JE, Hitz BC, Krieger CJ, Livstone MS, Miyasato SR, Nash RS, Oughtred R, Skrzypek MS, Weng S, Wong ED, Zhu KK, Dolinski K, Botstein D, Cherry JM: Gene ontology annotations at SGD: new data sources and annotation methods. Nucleic Acids Res. 2008, 36: D577-D581. 10.1093/nar/gkm909.PubMed CentralView ArticlePubMedGoogle Scholar
- Thomas PD, Mi H, Lewis S: Ontology annotation: mapping genomic regions to biological function. Curr Opin Chem Biol. 2007, 11: 4-11. 10.1016/j.cbpa.2006.11.039.View ArticlePubMedGoogle Scholar
- Moreno-Hagelsieb G, Latimer K: Choosing BLAST options for better detection of orthologs as reciprocal best hits. Bioinformatics. 2007, 319-324.Google Scholar
- Marchler-Bauer A, Anderson JB, Derbyshire MK, DeWeese-Scott C, Gonzales NR, Gwadz M, Hao L, He S, Hurwitz DI, Jackson JD, Ke Z, Krylov D, Lanczycki C, Liebert CA, Liu C, Lu F, Marchler GH, Mullokandov M, Song JS, Thanki N, Yamashita RA, Yin JJ, Zhang D, Bryant SH: CDD: a conserved domain database for interactive domain family analysis. Nucleic Acids Res. 2007, D237-40. 10.1093/nar/gkl951. 35 DatabaseGoogle Scholar
- Tatusov RL, Koonin EV, Lipman DJ: A Genomic perspective on protein families. Science. 1997, 278: 631-637. 10.1126/science.278.5338.631.View ArticlePubMedGoogle Scholar
- Hartl DL, Jones EW: Genetics: analysis of genes and genomes. 2004, Jones & Bartlett Publishers, SixthGoogle Scholar
- Deng J, Carbone I, Dean RA: The evolutionary history of cytochrome P450 genes in four filamentous Ascomycetes. BMC Evol Biol. 2007, 7: 30-10.1186/1471-2148-7-30.PubMed CentralView ArticlePubMedGoogle 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.