EPIC-DB: a proteomics database for studying Apicomplexan organisms
© Madrid-Aliste et al; licensee BioMed Central Ltd. 2009
Received: 05 September 2008
Accepted: 21 January 2009
Published: 21 January 2009
High throughput proteomics experiments are useful for analyzing the protein expression of an organism, identifying the correct gene structure of a genome, or locating possible post-translational modifications within proteins. High throughput methods necessitate publicly accessible and easily queried databases for efficiently and logically storing, displaying, and analyzing the large volume of data.
EPICDB is a publicly accessible, queryable, relational database that organizes and displays experimental, high throughput proteomics data for Toxoplasma gondii and Cryptosporidium parvum. Along with detailed information on mass spectrometry experiments, the database also provides antibody experimental results and analysis of functional annotations, comparative genomics, and aligned expressed sequence tag (EST) and genomic open reading frame (ORF) sequences. The database contains all available alternative gene datasets for each organism, which comprises a complete theoretical proteome for the respective organism, and all data is referenced to these sequences. The database is structured around clusters of protein sequences, which allows for the evaluation of redundancy, protein prediction discrepancies, and possible splice variants. The database can be expanded to include genomes of other organisms for which proteome-wide experimental data are available.
EPICDB is a comprehensive database of genome-wide T. gondii and C. parvum proteomics data and incorporates many features that allow for the analysis of the entire proteomes and/or annotation of specific protein sequences. EPICDB is complementary to other -genomics- databases of these organisms by offering complete mass spectrometry analysis on a comprehensive set of all available protein sequences.
High throughput proteomics experiments (mass spectrometry) provide unique insight into the protein expression profile of an organism and can be useful in identifying the correct gene structure of a genome, particularly when splice variants or alternative gene predictions exist, or to elucidate the roles of post-translational modifications within proteins . Indeed, large scale proteomics approaches have been used in recent years to analyze genomes of various organisms such as S. cerevisiae , M. mobile , C. parvum , T. gondii , and S. luteogriseus , as well as various subproteomes [7–9]. Recent advances in high throughput methods as well as the increased emphasis on and amount of high throughput data necessitate publicly accessible and easily queried databases for efficiently storing, displaying, and analyzing the large volume of data. We have established a proteomics database, EPICDB (Experimental ProteomICs Database) to accommodate experimental data from the large-scale proteomics exploration of two Apicomplexan organisms; Toxoplasma gondii and Cryptosporidium parvum.
T. gondii is an obligate intracellular protozoan, belonging to the phylum Apicomplexa, and is an important pathogen in both immune competent and immune compromised humans. The parasite causes chronic infection in adults and is present in an estimated 22.5% of people older than 12 in the United States  and up to 90% of the population in other regions of the world . T. gondii clinical disease is most typical in immune compromised individuals and congenitally infected children and is a common opportunistic pathogen associated with AIDS. Also, T. gondii is an important model system for the phylum Apicomplexa , which includes, among others, Plasmodium (malaria) species. Unlike many other Apicomplexa, which are experimentally intractable, T. gondii is easily cultured in vitro, has well established experimental protocols for genetic manipulation, and has a well characterized mouse model . Cryptosporidium parvum is an Apicomplexan parasite that infects the epithelial cells of the microvillus border of the gastrointestinal tract where it resides in a unique vacuole bellow the host cell membrane but outside the host cell cytoplasm . This organism forms resistant oocysts which are transmitted effectively by food or water. Ingestion of these oocysts results in infection of the gastrointestinal track with the development of diarrhea, which can be profuse especially in immune compromised hosts. Unlike T. gondii, in vitro culture is difficult and there are no established experimental protocols for genetic manipulation. Due to waterborne outbreaks associated with T. gondii  and C. parvum , these organisms are classified by the National Institute of Allergy and Infectious Diseases as Category B priority pathogens.
EPICDB contains high throughput mass spectrometry data that is cross referenced to all available computationally and experimentally derived protein sequences for T. gondii and C. parvum. The database can be queried to find proteins that are experimentally verified by proteomics data and, in doing, can aid in the discovery of sequences that have not been fully annotated but are supported by experimental data. Further, since all the proteomics data are comprehensively searched against all available protein sequences, EPICDB is a good resource for assessing discrepancies among sequences predicted by different algorithms and for examining splice variants of the same gene. Additionally, the format is suitable to accommodate additional experimental data or other organisms and can be used as a model system for other genomes for which high throughput proteomics data will become available.
The primary focus of EPICDB is to provide a comprehensive proteome-based description of the selected pathogens. This distinguishes it from other, genome-based databases, which exist for various apicomplexan and other organisms, such as, for example, ToxoDB http://www.toxodb.org, CryptoDB http://www.cryptodb.org, or TrichDB http://www.trichdb.org.
Construction and content
EPICDB is a relational database using a MySQL database management system and consisting of 40 tables. Within the database, query-response time is optimized by reading a summary table containing the results of a complex "JOIN" query.
EPICDB is a collection of experimental and computational data that are referenced to all available protein sequences for T. gondii and C. parvum, which represent the theoretical proteomes of the respective organisms. The available T. gondii protein sequences were compiled from five datasets; TigrScan , TwinScan , GlimmerHMM , ToxoDB.org Release4, and experimental sequences from the NCBI protein databases. TigrScan, TwinScan, and GlimmerHMM are computational gene finder algorithms that were employed by The Institute for Genomic Research (now the J. Craig Venter Institute) on the ME49 strain of T. gondii. The Release4 dataset is an annotation of the ME49 strain with sequences predicted by the GLEAN  algorithm. The dataset of experimental sequences was obtained from the NCBI Entrez Protein Database, which was filtered, by the organism name, for Toxoplasma gondii (As of July 2008, predicted genes/proteins from T. gondii genome analysis have not been deposited to NCBI or GenBank databases). The five datasets were combined to provide a comprehensive set of all available protein sequences for the Toxoplasma gondii proteome and comprised 30,197 sequences. While T. gondii has an intron rich genome and, consequently, presents substantial difficulty for gene prediction programs to properly identify gene coding regions, in C. parvum there are almost no introns and gene prediction is rather straightforward . The available C. parvum sequences were obtained from CryptoDB http://www.cryptodb.org and from the NCBI Entrez Protein Database, which was filtered, by the organism name, for Cryptosporidium Parvum. The combined C. parvum dataset comprises 8,316 protein sequences. In the case of both organisms the sequence sets are redundant due to computational algorithms predicting the same (or essentially same) sequences or the same sequences appearing more than once as experimentally derived sequences in the NCBI database.
Because of the redundancy among the alternate gene models, the database is organized in clusters of the protein sequences that share at least 90% sequence identity for the overlapping parts . Clustering allows the user to evaluate redundancy or prediction differences among different sequences. However, since the clustering scheme allows sequences to be grouped if they share local sequence similarity and not only sequences that are exact matches for the entire length, in many cases, clustering also offers the possibility of exploring alternative splicing events in the genome, when proteins are predicted from the same genomic location and differ by the inclusion, exclusion, elongation, or truncation of some introns or exons. The T. gondii dataset contains 30,197 alternative protein sequences, or possible splice variants, that are clustered into 14,983 groups of sequences, or possible protein coding genomic regions. The C. parvum dataset contains 8,316 sequences that are clustered into 3,852 groups.
EPICDB currently contains the results of 183 high throughput, tandem mass spectrometry experiments that were performed on T. gondii plasma membrane, cytoskeletal, cyst wall, and cytosolic protein preparations. These experiments produced 33,045 proteolytic peptides that, when searched against a database of the T. gondii protein sequences, experimentally verified 8,372 sequences (that can be grouped into 3,233 clusters). Further, EPICDB contains data from 88 high throughput, tandem mass spectrometry experiments that were performed on C. parvum. These experiments produced 3,049 proteolytic peptides that experimentally verified 1,241 C. parvum protein sequences (481 clusters). The data contained within EPICDB is a comprehensive summary of the mass spectrometry experiments and the MASCOT http://www.matrixscience.com searches of the data against the theoretical proteomes. For each protein, a list of the experiments that identified a proteolytic peptide that maps to that protein is included, along with the type of cell fraction analyzed, the type of mass spectrometry experiment conducted, and the number of peptide hits that were obtained. The mass spectrometry data and the results of the MASCOT search are provided and publicly available in five formats; CSV (comma separated value) formatted mass spectrometry files, RAW unformatted mass spectrometry data files, MGF (MASCOT Generic Format) and DTA files, containing peak lists of mass and intensity values, and MASCOT search summary files. Data for peptide scores, expectation values, location in the protein, and protein coverages are included for each proteolytic peptide that is mapped to the protein.
As an example of the addition of other types of experimental data, EPICDB also currently stores the results of antibody experiments that were performed on 52 T. gondii and 34 C. parvum protein sequences, including information on the peptide sequences used for immunization in the production of the antibodies, the amount of sera produced for each peptide mixture, and the results of Enzyme-Linked ImmunoSorbent Assays (ELISA), Immuno-Fluorescence Assays (IFA), and Immunoblot experiments, with images if available.
Additional mass spectrometry and antibody experiments are being continuously performed by the Albert Einstein Biodefense Proteomics Research Center and the results are incrementally added to EPICDB.
EPICDB contains functional annotations and comparative genomics for the theoretical proteomes of T. gondii and C. parvum. Transmembrane segments and signal peptides were predicted using the Phobius program . All known PFAM domains  were also identified. For each sequence in EPICDB, orthologous proteins were identified from sequences in the human genome, in other Apicomplexan genomes, and for the complete NCBI non-redundant protein database. Aligned genomic data (cDNA (ESTs) and ORFs) are also presented for each T. gondii protein sequence within the database. EST sequences were obtained from the NCBI EST "others" database and filtered for the organism name. ORF sequence data was obtained from ToxoDB.org version Release4. All experimental details and analysis of how mass spectrometry peptide hits validate alternative gene predictions can be found in a separate research article (J. Dybas et al. Computational Analysis and Experimental Validation of Gene Predictions in Toxoplasma gondii PLOSone. 2008;3(12):e3899. Epub 2008).
Utility and discussion
Querying the database
The database can be queried in a variety of ways. The Experimental Characterization query option selects the subset of protein sequences that are experimentally validated by mass spectrometry and/or antibody data. For the mass spectrometry experiments it is possible to select a subset of proteins based on the type of experimental cell fraction, such as membrane, cytosolic, or deglycosylated cell wall. The Annotations query is used to search for those proteins that are functionally characterized in terms of transmembrane domains, signal peptides and/or Pfam domains. Within the Pfam domains, it is possible to search by a specific domain name, such as "ACC_central", or by a description such as "Acetyl-CoA". The Comparative Genomics query locates predicted proteins that share a detectable sequence similarity (based on a BLAST search ) with other Apicomplexan genomes, with the human genome, or with any known proteins (referencing against the complete NCBI non-redundant protein database). Annotations from similar proteins are inherited and stored for each EPICDB protein sequence and the database search can be refined by a search for specific keywords, such as "transporter". Any or all of the aforementioned parameters can be selected to be combined in a single query.
In addition to seeing the query results for the entire theoretical genome, the Gene Predictions and Experimental Datasets query allows for the selection of proteins based on the type of gene prediction method or by only experimentally derived proteins. Within this query option one can search the database for a specific gene name, such as "TgTigrScan_1043".
Rather than querying the database based on experimental characteristics, it is possible to perform a BLAST search for a user-input protein sequence, or sequence fragment, against the entire EPICDB. The sensitivity of the BLAST search can be adjusted by defining the corresponding e-value cutoff.
Parsing the query results
Viewing experimental and computational data
Most options in the query page contain a link to a help page with a brief description of the specific parameter or option and, where applicable, a link to an associated web page. In addition, most navigating options, category headings, and table titles, on the query page, results page, and data pages contain pin-point helps providing a brief explanation of that entity. The upper right corner of the front (query) page has a link to three Flash formatted video tutorials that describe the options for querying the database and for understanding the query results page and the various data pages.
All available protein sequences for T. gondii and C. parvum that emerged from various gene prediction approaches or were experimentally derived are compiled into a hypothetical proteomes for each of the respective organisms in order to be able to search the proteomics data with the most comprehensive and unbiased dataset possible. Therefore, while some redundancy is obviously introduced, it also allows for the experimental verification and, potentially, the added scrutiny and study of many more sequences than what would be possible if only one set of protein predictions were adequately searched with the proteomics data.
EPICDB's functionality provides the options for a systems-level analysis as well as specific protein-level analysis. The proteome can be searched for any sequences that have experimental validation. Subsets of sequences can be found for those containing desired characteristics such as functional annotations or comparative genomics. Thus, the proteome can be studied by isolating specific sets of sequences that are of unique interest to the user. Within the obtained sequence sets, individual sequences can be analyzed. The individual sequence can be compared to similar sequences, via the clustering analysis, to determine any possible splice variability of the protein or instances in which there are disparities in the sequence predictions. Proteomics data can be analyzed to assess the validity of the predictions and/or splice variants by studying the distribution of assigned peptides on the amino acid sequence. In some cases a mass spectrometry peptide can indicate the possibility of an additional exon or the presence of incorrect splicing. Finally, EPICDB allows for the user to enter an amino acid sequence of interest and find any similar sequences in the database. This functionality provides the opportunity to examine any experimental evidence for a specific gene that the user is attempting to study and characterize.
At the present time the types of data contained within EPICDB are mass spectrometry and antibody experiments, EST and ORF sequences, functional annotations, and comparative genomics. However, the database can be easily expanded. We are currently planning to integrate genome-wide microarray expression data and ChIP-on-chip data into the database. Also, with little effort, other genomes can be added to the system.
While EPICDB provides a variety of automated annotations and all validation information, one of the major challenges is to properly annotate and analyze predicted genes and proteins on an individual basis. This type of work is best approached by a community effort where various groups have specific areas of expertise, perhaps regarding a specific protein or groups of proteins. Therefore we plan to establish a wiki-like option so experts can add information on annotations or functions of proteins.
High throughput mass spectrometry data has a variety of important applications to genomics and proteomics. Vast amount of genomic sequence data with comparatively small amounts of experimentally derived protein sequence data necessitate the need for computational gene finders and high throughput mass spectrometry data can inform and support these predictions. High throughput mass spectrometry is an emerging efficient method toward a comprehensive, proteome-wide analysis of an organism, which is an important step in identifying protein interactions, studying protein expression levels, elucidating alternative splice sites, or predicting potential chemotherapeutic targets. EPICDB is a comprehensive database for organizing, querying, and displaying proteome-wide proteomics data and incorporates many features that allow for the analysis of the entire proteomes and/or annotation of specific protein sequences for characterizing the proteomes of T. gondii and C. parvum.
Availability and requirements
EPICDB can be accessed by the website http://www.fiserlab.org/epicdb. The database and all data is publicly available. There are, on average, 6 unique users per day (statistics do not count access from Albert Einstein College of Medicine computers) and 36.4 Mb of data transferred.
Financial support was provided by NIH-NIAID HHSN266200400054C. AF was supported by the HuBi (Hungarian Bioinformatics) project, in the framework of the European Community's "Structuring the European Research Area" programme.
- Fauquenoy S, Morelle W, Hovasse A, Bednarczyk A, Slomianny C, Schaeffer C, Van Dorsselaer A, Tomavo S: Proteomics and glycomics analyses of N-glycosylated structures involved in Toxoplasma gondii–host cell interactions. Mol Cell Proteomics. 2008, 7 (5): 891-910.View ArticlePubMedGoogle Scholar
- Washburn MP, Wolters D, Yates JR: Large-scale analysis of the yeast proteome by multidimensional protein identification technology. Nat Biotechnol. 2001, 19 (3): 242-247.View ArticlePubMedGoogle Scholar
- Jaffe JD, Berg HC, Church GM: Proteogenomic mapping as a complementary method to perform genome annotation. Proteomics. 2004, 4 (1): 59-77.View ArticlePubMedGoogle Scholar
- Sanderson SJ, Xia D, Prieto H, Yates J, Heiges M, Kissinger JC, Bromley E, Lal K, Sinden RE, Tomley F: Determining the protein repertoire of Cryptosporidium parvum sporozoites. Proteomics. 2008, 8 (7): 1398-1414.PubMed CentralView ArticlePubMedGoogle Scholar
- Xia D, Sanderson SJ, Jones AR, Prieto JH, Yates JR, Bromley E, Tomley FM, Lal K, Sinden RE, Brunk BP: The proteome of Toxoplasma gondii: integration with the genome provides novel insights into gene expression and annotation. Genome Biol. 2008, 9 (7): R116-PubMed CentralView ArticlePubMedGoogle Scholar
- Wang R, Prince JT, Marcotte EM: Mass spectrometry of the M. smegmatis proteome: protein expression levels correlate with function, operons, and codon bias. Genome Res. 2005, 15 (8): 1118-1126.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhou XW, Kafsack BF, Cole RN, Beckett P, Shen RF, Carruthers VB: The opportunistic pathogen Toxoplasma gondii deploys a diverse legion of invasion and survival proteins. J Biol Chem. 2005, 280 (40): 34233-34244.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhou XW, Blackman MJ, Howell SA, Carruthers VB: Proteomic analysis of cleavage events reveals a dynamic two-step mechanism for proteolysis of a key parasite adhesive complex. Mol Cell Proteomics. 2004, 3 (6): 565-576.View ArticlePubMedGoogle Scholar
- Bradley PJ, Ward C, Cheng SJ, Alexander DL, Coller S, Coombs GH, Dunn JD, Ferguson DJ, Sanderson SJ, Wastling JM: Proteomic analysis of rhoptry organelles reveals many novel constituents for host-parasite interactions in Toxoplasma gondii. J Biol Chem. 2005, 280 (40): 34245-34258.View ArticlePubMedGoogle Scholar
- Jones JL, Kruszon-Moran D, Wilson M, McQuillan G, Navin T, McAuley JB: Toxoplasma gondii infection in the United States: seroprevalence and risk factors. Am J Epidemiol. 2001, 154 (4): 357-365.View ArticlePubMedGoogle Scholar
- Sousa OE, Saenz RE, Frenkel JK: Toxoplasmosis in Panama: a 10-year study. Am J Trop Med Hyg. 1988, 38 (2): 315-322.PubMedGoogle Scholar
- Roos DS, Crawford MJ, Donald RG, Fohl LM, Hager KM, Kissinger JC, Reynolds MG, Striepen B, Sullivan WJ: Transport and trafficking: Toxoplasma as a model for Plasmodium. Novartis Found Symp. 1999, 226: 176-195. discussion 195-178PubMedGoogle Scholar
- Kim K, Weiss LM: Toxoplasma gondii: the model apicomplexan. Int J Parasitol. 2004, 34 (3): 423-432.PubMed CentralView ArticlePubMedGoogle Scholar
- Fayer R, Xiao L, eds: Cryptosporidium and Cryptosporidosis. 2008, Boca Raton FL: CRC Press, 2
- Bowie WR, King AS, Werker DH, Isaac-Renton JL, Bell A, Eng SB, Marion SA: Outbreak of toxoplasmosis associated with municipal drinking water. The BC Toxoplasma Investigation Team. Lancet. 1997, 350 (9072): 173-177.View ArticlePubMedGoogle Scholar
- Mac Kenzie WR, Hoxie NJ, Proctor ME, Gradus MS, Blair KA, Peterson DE, Kazmierczak JJ, Addiss DG, Fox KR, Rose JB: A massive outbreak in Milwaukee of cryptosporidium infection transmitted through the public water supply. N Engl J Med. 1994, 331 (3): 161-167.View ArticlePubMedGoogle Scholar
- Majoros WH, Pertea M, Salzberg SL: TigrScan and GlimmerHMM: two open source ab initio eukaryotic gene-finders. Bioinformatics. 2004, 20 (16): 2878-2879.View ArticlePubMedGoogle Scholar
- Korf I, Flicek P, Duan D, Brent MR: Integrating genomic homology into gene structure prediction. Bioinformatics. 2001, 17 (Suppl 1): S140-148.View ArticlePubMedGoogle Scholar
- Elsik CG, Mackey AJ, Reese JT, Milshina NV, Roos DS, Weinstock GM: Creating a honey bee consensus gene set. Genome Biol. 2007, 8 (1): R13-PubMed CentralView ArticlePubMedGoogle Scholar
- Abrahamsen MS, Templeton TJ, Enomoto S, Abrahante JE, Zhu G, Lancto CA, Deng M, Liu C, Widmer G, Tzipori S: Complete genome sequence of the apicomplexan, Cryptosporidium parvum. Science. 2004, 304 (5669): 441-445.View ArticlePubMedGoogle Scholar
- Li W, Godzik A: Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006, 22 (13): 1658-1659.View ArticlePubMedGoogle Scholar
- Kall L, Krogh A, Sonnhammer EL: A combined transmembrane topology and signal peptide prediction method. J Mol Biol. 2004, 338 (5): 1027-1036.View ArticlePubMedGoogle Scholar
- Finn RD, Mistry J, Schuster-Bockler B, Griffiths-Jones S, Hollich V, Lassmann T, Moxon S, Marshall M, Khanna A, Durbin R: Pfam: clans, web tools and services. Nucleic Acids Res. 2006, D247-251. 34 Database
- 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 (17): 3389-3402.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.