- Open Access
Integration with the human genome of peptide sequences obtained by high-throughput mass spectrometry
- Frank Desiere†1, 2,
- Eric W Deutsch†2,
- Alexey I Nesvizhskii†2,
- Parag Mallick†2,
- Nichole L King2,
- Jimmy K Eng2,
- Alan Aderem2,
- Rose Boyle2,
- Erich Brunner2, 3,
- Samuel Donohoe2,
- Nelson Fausto4,
- Ernst Hafen3,
- Lee Hood2,
- Michael G Katze5,
- Kathleen A Kennedy2,
- Floyd Kregenow2,
- Hookeun Lee2,
- Biaoyang Lin2,
- Dan Martin2,
- Jeffrey A Ranish2,
- David J Rawlings6,
- Lawrence E Samelson7,
- Yuzuru Shiio2,
- Julian D Watts2,
- Bernd Wollscheid2,
- Michael E Wright2,
- Wei Yan2,
- Lihong Yang8,
- Eugene C Yi2,
- Hui Zhang2 and
- Ruedi Aebersold2, 9Email author
© Desiere et al.; licensee BioMed Central Ltd. 2004
Received: 1 September 2004
Accepted: 17 November 2004
Published: 10 December 2004
A crucial aim upon the completion of the human genome is the verification and functional annotation of all predicted genes and their protein products. Here we describe the mapping of peptides derived from accurate interpretations of protein tandem mass spectrometry (MS) data to eukaryotic genomes and the generation of an expandable resource for integration of data from many diverse proteomics experiments. Furthermore, we demonstrate that peptide identifications obtained from high-throughput proteomics can be integrated on a large scale with the human genome. This resource could serve as an expandable repository for MS-derived proteome information.
The recent definition of the complete nucleotide sequence of the human genome [1, 2] has motivated the full annotation of the sequence. The true promise of the human genome project, to become the foundation for medical and biological research benefiting human health and quality of life , can only be realized if the coding sequences are conclusively identified, intron/exon structures are accurately described and the potential protein products from each gene in different tissues and cellular states are determined. Current methods for gene-prediction provide useful information but are still limited . It is not presently possible to predict all features of the genome from its sequence alone. Therefore, the value of the human genome sequence can be enhanced through the collection of different types of experimental data and its integration and validation in a genomic context .
Current use of expressed sequence tags (EST) and full coding DNA (cDNA) sequences is extremely helpful in achieving complete genome annotation [6–9]. However, these data are not sufficient to unequivocally predict which proteins (and with what covalent structure) are expressed in a given tissue. The complete characterization of all proteins across disease states, tissues and stages of development can now be addressed through experimental protein identifications generated by proteomic methods. Experiments carried out over the past years have illustrated that peptides resulting from proteolytic digests of complex protein mixtures can be identified in a high-throughput mode using a combination of liquid chromatography (LC) and tandem mass spectrometry (MS/MS) (LC-MS/MS) [10–15]. Peptides are thus useful as the currency of MS/MS-based protein identification . By combining a large number of experiments sampling different cell and tissue types, the observed peptides can be mapped onto the genome covering a significant part of its chromosomes.
Results and discussion
To begin annotating the human genome with protein-level information, we have built PeptideAtlas. The generally applicable procedure to annotate eukaryotic genomes with peptide sequences can be applied when datasets are acquired using different experimental protocols. In each case, sample proteins were first proteolytically cleaved into peptides using the enzyme trypsin. The resulting peptide mixture was then subjected to chromatographic separation by strong cation exchange and reverse-phase capillary chromatography. In addition, those experiments using the ICAT (isotope-coded affinity tag) reagent for quantification included an avidin affinity-purification step to select peptides containing biotinylated, stable-isotope-tagged cysteines . The resulting peptide pools were then analyzed by electrospray ionization (ESI)-MS/MS. The database search program SEQUEST  was used to assign the resulting MS/MS spectra to a peptide sequence. The confidence of these peptide assignments was evaluated using PeptideProphet . All of the experimental data products, including PeptideProphet probability scores, are loaded into SBEAMS - Proteomics, a proteomics analysis database built as a module under the Systems Biology Experiment Analysis Management System (SBEAMS) framework. All of the identifications above a certain probability threshold within a specific set of experiments are extracted from the main database tables into another set of tables containing the attributes of each distinct peptide.
Summary of PeptideAtlas results
Ensembl gene predictions
13525 from Release 3.1 FlyBase
Ensembl gene transcripts
Number of experiments
PeptideProphet probability threshold
PeptideAtlas mapped peptides
PeptideAtlas mapped proteins
PeptideAtlas mapped genes
Percentage of the genome
The current build of PeptideAtlas contains peptide sequences identified in 52 proteomic experiments in which proteins were extracted from a particular cell or tissue type, digested with trypsin and analyzed with a mass spectrometer. The 52 proteomic experiments comprised 14 published as well as 38 unpublished human datasets from various cell types such as T cells, B cells, lymphocytes, lymphoblasts, hepatocytes, intestinal cells, hepatoma cells and others. The 14 published datasets contain 47% of the distinct peptides in PeptideAtlas. A full listing of all the experiments and samples currently in PeptideAtlas can be found at the project website . The raw data for all published datasets is also provided in a repository there.
Applying our pipeline described in Materials and methods, 25,754 of the 26,840 distinct peptides in PeptideAtlas were mapped to 9,747 (28.6 %) of the 34,091 human Ensembl proteins (version 22.34d.1, 2004-06-02). These proteins represent unique proteins or splice forms from 6,423 genes (27%) of human genes in Ensembl.
Some peptides have indistinguishable, perfect protein sequence matches to multiple proteins. These proteins are typically paralogs (protein families), protein isoforms or repeated protein domains in the human genome. We identified 3,718 proteins unambiguously by one or more 'discrete peptides' - peptides that map uniquely to a single protein - in the current build of PeptideAtlas. Those peptides are marked in the genome browser as 'discrete peptide'. 'Degenerate peptides' that map to several protein isoforms are also used to identify proteins. It would thus be more accurate to state that a product of a certain gene, rather than a certain protein, has been identified [20–22]. Moreover, the experimental data from those degenerate peptides generally do not allow differentiation between the sequence alternatives that exist in Ensembl. In fact, not all splicing variants that are in Swiss-Prot are also present in Ensembl and, therefore, it is impossible to ascertain the number of unambiguous identifications at the moment. This limitation underscores the requirement for mapping large-scale proteomic data to the human genome, such as presented in this report to aid in the generation of unambiguous sequence databases.
A significant number of distinct peptides (1,086), assigned by SEQUEST/ProteinProphet from over 5,000 MS/MS spectra, could not be mapped to Ensembl database version 22.34d.1. These peptides were identified by SEQUEST searches against the IPI database  or ABCC non-redundant protein database (NCI) . These peptides are of special interest as they often document interesting biological phenomena such as single-nucleotide polymorphisms (SNPs) and novel splice variants, demonstrating the need for annotating the human genome sequence with high-quality experimental data obtained from expressed proteins. The existence of these sequences also illustrates the flux in the genome annotation and sequence databases. For example, in Ensembl version 18.34.1, only 92% of genes from the previous build were transferred across to the new build. The missing 8% were predominantly inappropriate protein-coding genes coming from large-scale cDNA projects, which have a number of artifactual errors, or from chimeric cDNA clones from cancer cell lines. Experimentally observed, unmapped peptides are an ideal source of information for refining genome assembly and gene prediction.
The absence of Ensembl matches does raise the question of whether these peptides are false positives or whether real proteins are missing in the Ensembl database. When these peptides were investigated in more detail it was found that nearly 100 were identified 10 or more times in several different experiments, and that many had protein sequence matches for Swiss-Prot entries. They are therefore likely to be true peptide attributions. For example, peptide PAp00000363 (AGKPVICATQMLESMIK) was identified 626 times at different charge states and with different mass modifications in 22 distinct experiments and mapped to KPY1_HUMAN, a pyruvate kinase M1 isozyme. Interestingly, the protein appears to have a likely SNP, which mutates the valine present in the Ensembl genome sequence to the isoleucine observed in PAp00000363.
The development of PeptideAtlas and a method for mapping observed peptides to the genome allows us to determine the distribution of multiple peptide hits to specific proteins and the distribution of peptide sequences that are present in multiple proteins. Also, in some cases splice junctions and gene boundaries could be confirmed. Our method allows us also to identify peptides corresponding to abundant proteins such as actin, elongation factor and glyceraldehyde-3-phosphate dehydrogenase, which are commonly identified in high-throughput LC-MS/MS experiments. These proteins are products of housekeeping genes, which are expressed most of the time in almost every tissue , or are structural proteins which are also known to be abundant in cells.
The need for public proteomics data repositories is recognized  and we intend PeptideAtlas to become a growing database and public resource. We have structured the system in a way that allows scientists to submit their own MS data for incorporation into PeptideAtlas, thus increasing the number of experiments and identified peptides. Naturally, to be useful for the project, inclusion of third-party data is dependent upon data compatibility and consistent data quality. Consequently, only data with accurate statistical measures of confidence computed by, for example, PeptideProphet, or another published and tested statistical algorithm, will be included. Datasets for which such statistical analyses have been performed can be submitted for incorporation following the procedure detailed at the PeptideAtlas website. Alternatively, data contributors can submit raw MS/MS data directly. This information should preferably be formatted into mzXML  or mzData (HUPO Proteomics Standards Initiative) which are open file formats for the representation of MS data. Other traditionally used data formats are accepted as well.
This data will then be searched by the PeptideAtlas curators using SEQUEST to correlate MS/MS spectra of peptides with amino-acid sequences using protein databases such as IPI, and the results will be further analyzed with PeptideProphet. An effort to add support for additional search engines is underway. This procedure will ensure the highest degree of consistency for the data in PeptideAtlas. In the future, the pipeline in general and the data submission process in particular, can be further improved and make compliant with the community accepted statistical data-validation standards and data file formats when such standards emerge . Please see the submission section on the PeptideAtlas web-site for the most up-to-date submission methods and curator contact information. With an increasing number of included peptides, the utility of the resource will improve, as increasing numbers of genes, exons, transcripts and variant transcripts in many tissues and developmental stages will be verified on the protein level.
All MS/MS spectra are stored in the SBEAMS - Proteomics database, from which PeptideAtlas is derived. While at present it is not possible easily to access the MS/MS spectra starting from the public PeptideAtlas interface, this possibility could be added in the future. All spectra for published experiments are available in the mzXML files in the repository. Access to raw spectra can be beneficial for many applications not related to the main purpose of PeptideAtlas. Furthermore, because peptide modifications (for example, phosphorylation) are stored, this information could be displayed as well.
It is well understood and discussed in the literature  that all large-scale datasets obtained using high-throughput methods inherently contain a certain fraction of false-positive data. Thus, estimation of false-positive error rates is a very important but often challenging task. One significant advantage of the high-throughput pipeline implemented in this work is that computed peptide probabilities (here produced by PeptideProphet) allow estimation of the upper bound (most conservative estimate) of the false-positive identification error rates for any dataset submitted to PeptideAtlas. As the main purpose of PeptideAtlas is to map peptide identifications to the genome, the most relevant estimate of the false-positive error rates is the one at the level of distinct peptide assignments that have a defined mapping to Ensembl.
Comparison of different probability thresholds that were applied to the MS results
Total number of passing spectra
Distinct peptides with protein sequence matches
Number of mapped proteins
Number of simple reduced proteins
False-positive estimate MS/MS spectra
False-positive estimate with protein sequence matches
To assess the effect of using a particular probability threshold on the number of peptides in the atlas, we ran the PeptideAtlas pipeline using probability thresholds P ≥ 0.7, 0.9, 0.95 and 0.99. Decreasing the probability threshold increases the number of peptides, both correctly and incorrectly identified, and the corresponding proteins (Table 2). The most stringent threshold of P ≥ 0.99 produced 21,030 peptides with protein sequence matches (4,845 protein identifications), almost 8,400 fewer than the lowest threshold of 0.7 (2,252 fewer protein identifications). The P ≥ 0.9 threshold yielded 25,754 peptides with protein sequence matches at an estimated false-positive rate of less than 7%, and we selected this as an acceptable level for the default PeptideAtlas. The number of false-positive identifications could be reduced by selecting a higher threshold; however, a significant number of correct peptides and proteins would then also be eliminated. The additional peptides resulting from the low-probability threshold were valuable for adding additional peptide evidence in combination with higher-probability peptides corresponding to the same protein (peptides corresponding to proteins to which other peptides correspond are more likely to be correct than their probability value indicates ). We provide at our website the option for users to browse or download versions of the Atlas generated with the other P thresholds, which might be useful for some applications.
To validate our approach for general use in eukaryote genomes, we have extended our methods to peptides obtained from Drosophila melanogaster LC-MS/MS experiments. We collected data obtained from cytoplasmic, nuclear and membrane fractions derived from a Drosophila S2 Schneider cell line. The resulting 4,406 different peptides with P > 0.9 were compared to the 18,289 proteins (Ensembl fly database version 18.3a.1, 2003-07-01) using the same pipeline as described for human. From the fly, 3,107 proteins could be validated, representing 1,876 (14%) of the fly's genes. These results show that our method could easily be adapted to other organisms, thus opening up the way for comparative proteome-level evaluations of eukaryotic organisms.
We have annotated the human genome with protein evidence for nearly 10,000 proteins. Although this number only represents a fraction of the genome and still contains some erroneous identifications, it is a first step towards the final goal: to fully annotate eukaryotic genomes via validation of expressed proteins. PeptideAtlas provides a method and a framework to accommodate proteome information generated by high-throughput proteomics technologies and is able to efficiently disseminate experimental data in the public domain. Its significance continues to grow as more data are submitted.
Moreover, PeptideAtlas also allows one to address the important question of how big the human proteome is. Due to the technical limitations of current proteomics technologies, it is not possible yet to determine the complete proteome in one experiment. However, if the data from diverse experiments, using different cellular compartments and enrichment methods were combined, the determination of the complete proteome could eventually be achieved. PeptideAtlas offers the framework to answer this question accurately and to determine the size of the complete human proteome using pooled experimental data. Furthermore, PeptideAtlas provides a resource for the development of new avenues of research. The dataset will provide a rich source of data for computational scientists to develop and test new algorithms for proteomic analysis, gene discovery and splice-variant prediction.
The methods described here, combined with the ever-increasing power of proteomics and bioinformatics technologies, will facilitate the determination or characterization of protein-coding genes, their features, and their processing and expression in relationship to the sequence of the human genome, thus contributing significantly to our understanding of genome structure.
Materials and methods
The assembly of experimentally derived distinct peptides is mapped to the human genome in the following way. First, we use BLAST  to match the peptides to the Ensembl human protein database. The Ensembl database project  provides a bioinformatics framework to organize biology around the sequences of large genomes and, furthermore, extensive resources and visualization options as well as remote access to the underlying relational databases . The human genome sequence (release 22.34d.1, 2004-06-02) contains 23,758 genes and 34,091 gene transcripts. Second, complete matches, spanning each peptide's complete length, were used to determine human chromosomal coordinates. The method for retrieving chromosomal coordinates within the human genome accounts for splice junctions; in cases where a peptide maps onto a splice junction, it is projected to both parts of the chromosome, generating multiple sets of coordinates. Third, the results are loaded into a relational database. This database schema (available at the project website ) is able to accommodate data for different PeptideAtlas builds, for different organisms or different reference protein sequence sets as starting material and is thus extremely versatile. Fourth, visualization of the results was achieved using the Distributed Annotation System (DAS) (Figure 2) in conjunction with the Ensembl database. DAS allows sequence annotations to be decentralized among multiple third-party annotators and integrated on an as-needed basis by the Ensembl genome browser .
LC-MS/MS analysis was performed on LCQ, Ion-trap (Thermo Finnigan LCQ) and Q-Tof (Micromass Waters) instruments.
To estimate the false-positive error rate on the level of distinct peptide identifications, we first note that there is an almost 10-fold difference between the number of peptide assignments to MS/MS spectra and the number of resulting distinct peptide identifications. This can be explained by the fact that many peptides were sequenced multiple times, with some of the most abundant peptides sequenced more than 1,000 times (for example, peptides PAp00004784, PAp00003568, PAp00026910). While many correct peptide assignments to MS/MS spectra represent the same peptide sequence, the majority of incorrect peptide assignments are expected to be single identifications. As a result, the false-positive error rate on the level of distinct peptides is higher than that on the level of peptide assignments to MS/MS spectra.
Second, it should also be taken into account that a considerable fraction of all distinct peptides did not match any Ensembl entry. This is due to the fact that MS/MS spectra were searched against larger databases, such as human IPI, which contained a number of protein sequences not present in Ensembl. The fraction of all distinct peptide identifications that did not map to any Ensembl entry can be estimated using information provided in Table 2. Among peptides with probability of being correct of 0.99 or greater, only 2.6% of all distinct peptides did not map to any Ensembl sequence. The fraction of unmapped distinct peptides increases to 8.3% among peptides in the 0.95-0.99 probability range, 12.9% in the 0.9-0.95 range and 18.2% in the 0.7-0.9 range, reflecting the increase in the number of incorrect peptide identifications among peptides with lower probabilities. Thus, one can estimate that at least 18.2% of all incorrectly identified peptides did not map to any entry in Ensembl.
The false-positive error rate among distinct peptides that mapped to Ensembl (peptides with protein sequence match) can then be estimated to be not higher than the maximum possible number of distinct incorrect peptide identifications that have protein sequence matches (computed by multiplying the total number of peptide assignments to MS/MS spectra by the corresponding false-positive error rate and applying an 18.2% correction to account for peptides with no mapping to Ensembl) divided by the total number of peptides with protein sequence matches. The corresponding estimates are 16%, 6%, 3% and 0.8% in the case of minimum-probability thresholds 0.7, 9, 0.95 and 0.99, respectively (Table 2). It should be noted that these are conservative (upper bound) estimates and the actual error rate may be significantly smaller.
Population of the database
The PeptideAtlas pipeline begins with the download of the Ensembl human protein database from . Release 22.34d.1 (2004-06-02) was used here. PeptideAtlas peptides were then searched against the human proteins using BLAST with the following parameters adapted for searching small peptides : -E 1 -W 2 -M PAM30 -G 9 -e 10 -K 50 -b 50 -F F. The BLAST results were then filtered for identical matches and mapped into chromosomal coordinates using Bio::EnsEMBL and Bioperl  Perl modules. The results are uploaded into the PeptideAtlas database and then the Ensembl genome browser. The PeptideAtlas database can handle different PeptideAtlas builds, different organisms and different versions of underlying genome data for maximum flexibility.
We thank Mike Carlson, B. Brett Finlay, Philip R. Hardwidge, Stephen D. Hauschka, Charis L. Himeda, and Priska von Haller for making the raw data from their published results available for this study. This project has been funded in part with federal funds from the National Heart, Lung, and Blood Institute, National Institutes of Health, under contract No. N01-HV-28179, and from the National Institute on Drug Abuse under contract P30DA015625. Some unpublished data contributing to this work was supported in part by a grant from the NIH to J.W. (RO1-AI-51344-01).
- Venter JC, Adams MD, Myers EW, Li PW, Mural RJ, Sutton GG, Smith HO, Yandell M, Evans CA, Holt RA, et al: The sequence of the human genome. Science. 2001, 291: 1304-1351. 10.1126/science.1058040.View ArticleGoogle Scholar
- Lander ES, Linton LM, Birren B, Nusbaum C, Zody MC, Baldwin J, Devon K, Dewar K, Doyle M, FitzHugh W, et al: Initial sequencing and analysis of the human genome. Nature. 2001, 409: 860-921. 10.1038/35057062.View ArticleGoogle Scholar
- Collins FS, Green ED, Guttmacher AE, Guyer MS: A vision for the future of genomics research. Nature. 2003, 422: 835-847. 10.1038/nature01626.View ArticleGoogle Scholar
- Pennisi E: Bioinformatics. Gene counters struggle to get the right answer. Science. 2003, 301: 1040-1041. 10.1126/science.301.5636.1040.View ArticleGoogle Scholar
- Birney E, Clamp M, Hubbard T: Databases and tools for browsing genomes. Annu Rev Genomics Hum Genet. 2002, 3: 293-310. 10.1146/annurev.genom.3.030502.101529.View ArticleGoogle Scholar
- Imanishi T, Itoh T, Suzuki Y, O'Donovan C, Fukuchi S, Koyanagi KO, Barrero RA, Tamura T, Yamaguchi-Kabata Y, Tanino M, et al: Integrative annotation of 21,037 human genes validated by full-length cDNA clones. PLoS Biol. 2004, 2: 856-875. 10.1371/journal.pbio.0020162.View ArticleGoogle Scholar
- Furuno M, Kasukawa T, Saito R, Adachi J, Suzuki H, Baldarelli R, Hayashizaki Y, Okazaki Y: CDS annotation in full-length cDNA sequence. Genome Res. 2003, 13: 1478-1487. 10.1101/gr.1060303.View ArticleGoogle Scholar
- Numata K, Kanai A, Saito R, Kondo S, Adachi J, Wilming LG, Hume DA, Hayashizaki Y, Tomita M, RIKEN GER Group, GSL Members: Identification of putative noncoding RNAs among the RIKEN mouse full-length cDNA collection. Genome Res. 2003, 13: 1301-1306. 10.1101/gr.1011603.View ArticleGoogle Scholar
- de Souza SJ, Camargo AA, Briones MR, Costa FF, Nagai MA, Verjovski-Almeida S, Zago MA, Andrade LE, Carrer H, El-Dorry HF, et al: Identification of human chromosome 22 transcribed sequences with ORF expressed sequence tags. Proc Natl Acad Sci USA. 2000, 97: 12690-12693. 10.1073/pnas.97.23.12690.View ArticleGoogle Scholar
- Washburn MP, Wolters D, Yates JR: Large-scale analysis of the yeast proteome by multidimensional protein identification technology. Nat Biotechnol. 2001, 19: 242-247. 10.1038/85686.View ArticleGoogle Scholar
- Florens L, Washburn MP, Raine JD, Anthony RM, Grainger M, Haynes JD, Moch JK, Muster N, Sacci JB, Tabb DL, et al: A proteomic view of the Plasmodium falciparum life cycle. Nature. 2002, 419: 520-526. 10.1038/nature01107.View ArticleGoogle Scholar
- Lasonder E, Ishihama Y, Andersen JS, Vermunt AM, Pain A, Sauerwein RW, Eling WM, Hall N, Waters AP, Stunnenberg HG, Mann M: Analysis of the Plasmodium falciparum proteome by high-accuracy mass spectrometry. Nature. 2002, 419: 537-542. 10.1038/nature01111.View ArticleGoogle Scholar
- Kuster B, Mortensen P, Andersen JS, Mann M: Mass spectrometry allows direct identification of proteins in large genomes. Proteomics. 2001, 1: 641-650. 10.1002/1615-9861(200104)1:5<641::AID-PROT641>3.3.CO;2-I.View ArticleGoogle Scholar
- Choudhary JS, Blackstock WP, Creasy DM, Cottrell JS: Interrogating the human genome using uninterpreted mass spectrometry data. Proteomics. 2001, 1: 651-667. 10.1002/1615-9861(200104)1:5<651::AID-PROT651>3.3.CO;2-E.View ArticleGoogle Scholar
- Han DK, Eng J, Zhou H, Aebersold R: Quantitative profiling of differentiation-induced microsomal proteins using isotope-coded affinity tags and mass spectrometry. Nat Biotechnol. 2001, 19: 946-951. 10.1038/nbt1001-946.View ArticleGoogle Scholar
- Aebersold R, Mann M: Mass spectrometry-based proteomics. Nature. 2003, 422: 198-207. 10.1038/nature01511.View ArticleGoogle Scholar
- Eng J, McCormack AL, Yates JR: An approach to correlate tandem mass spectral data of peptides with amino acid sequences in a protein database. J Am Soc Mass Spectrom. 1994, 5: 976-989. 10.1016/1044-0305(94)80016-2.View ArticleGoogle Scholar
- Keller A, Nesvizhskii AI, Kolker E, Aebersold R: Empirical statistical model to estimate the accuracy of peptide identifications made by MS/MS and database search. Anal Chem. 2002, 74: 5383-5392. 10.1021/ac025747h.View ArticleGoogle Scholar
- PeptideAtlas home. [http://www.peptideatlas.org]
- Rappsilber J, Mann M: What does it mean to identify a protein in proteomics?. Trends Biochem Sci. 2002, 27: 74-78. 10.1016/S0968-0004(01)02021-7.View ArticleGoogle Scholar
- Nesvizhskii AI, Aebersold R: Analysis, statistical validation and dissemination of large-scale proteomics datasets generated by tandem MS. Drug Discov Today. 2004, 9: 173-181. 10.1016/S1359-6446(03)02978-7.View ArticleGoogle Scholar
- Nesvizhskii AI, Keller A, Kolker E, Aebersold R: A statistical model for identifying proteins by tandem mass spectrometry. Anal Chem. 2003, 75: 4646-4658. 10.1021/ac0341261.View ArticleGoogle Scholar
- Kersey PJ, Duarte J, Williams A, Karavidopoulou Y, Birney E, Apweiler R: The International Protein Index: an integrated database for proteomics experiments. Proteomics. 2004, 4: 1985-1988. 10.1002/pmic.200300721.View ArticleGoogle Scholar
- NRP (Non-Redundant Protein) Database: National Cancer Institute Advanced Biomedical Computing Center, 2004. [AU: please give a fuller ftp address for this database], [ftp://ftp.ncifcrf.gov/pub/nonredun]
- Eisenberg E, Levanon EY: Human housekeeping genes are compact. Trends Genet. 2003, 19: 362-365. 10.1016/S0168-9525(03)00140-9.View ArticleGoogle Scholar
- Aebersold R: Constellations in a cellular universe. Nature. 2003, 422: 115-116. 10.1038/422115a.View ArticleGoogle Scholar
- Machiels BM, Zorenc AH, Endert JM, Kuijpers HJ, van Eys GJ, Ramaekers FC, Broers JL: An alternative splicing product of the lamin A/C gene lacks exon 10. J Biol Chem. 1996, 271: 9249-9253. 10.1074/jbc.271.16.9249.View ArticleGoogle Scholar
- Prince JT, Carlson MW, Wang R, Lu P, Marcotte EM: The need for a public proteomics repository. Nat Biotechnol. 2004, 22: 471-472. 10.1038/nbt0404-471.View ArticleGoogle Scholar
- Pedrioli PG, Eng JK, Hubley R, Vogelzang M, Deutsch EW, Raught B, Pratt B, Nilsson E, Angeletti RH, Apweiler R, et al: A common open representation of mass spectrometry data and its application to proteomics research. Nat Biotechnol. 2004, 22: 1459-1466. 10.1038/nbt1031.View ArticleGoogle Scholar
- Orchard S, Hermjakob H, Julian RK, Runte K, Sherman D, Wojcik J, Zhu W, Apweiler R: Common interchange standards for proteomics data: Public availability of tools and schema. Proteomics. 2004, 4: 490-491. 10.1002/pmic.200300694.View ArticleGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410. 10.1006/jmbi.1990.9999.View ArticleGoogle Scholar
- Birney E, Andrews D, Bevan P, Caccamo M, Cameron G, Chen Y, Clarke L, Coates G, Cox T, Cuff J, et al: Ensembl 2004. Nucleic Acids Res. 2004, 32 (Database issue): D468-D470. 10.1093/nar/gkh038.View ArticleGoogle Scholar
- Clamp M, Andrews D, Barker D, Bevan P, Cameron G, Chen Y, Clark L, Cox T, Cuff J, Curwen V, et al: Ensembl 2002: accommodating comparative genomics. Nucleic Acids Res. 2003, 31: 38-42. 10.1093/nar/gkg083.View ArticleGoogle Scholar
- Dowell RD, Jokerst RM, Day A, Eddy SR, Stein L: The distributed annotation system. BMC Bioinformatics. 2001, 2: 7-10.1186/1471-2105-2-7.View ArticleGoogle Scholar
- Ensembl. [http://www.ensembl.org]
- Stajich JE, Block D, Boulez K, Brenner SE, Chervitz SA, Dagdigian C, Fuellen G, Gilbert JG, Korf I, Lapp H, et al: The Bioperl toolkit: Perl modules for the life sciences. Genome Res. 2002, 12: 1611-1618. 10.1101/gr.361602.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.