Initial characterization of the human central proteome
BMC Systems Biology volume 5, Article number: 17 (2011)
On the basis of large proteomics datasets measured from seven human cell lines we consider their intersection as an approximation of the human central proteome, which is the set of proteins ubiquitously expressed in all human cells. Composition and properties of the central proteome are investigated through bioinformatics analyses.
We experimentally identify a central proteome comprising 1,124 proteins that are ubiquitously and abundantly expressed in human cells using state of the art mass spectrometry and protein identification bioinformatics. The main represented functions are proteostasis, primary metabolism and proliferation. We further characterize the central proteome considering gene structures, conservation, interaction networks, pathways, drug targets, and coordination of biological processes. Among other new findings, we show that the central proteome is encoded by exon-rich genes, indicating an increased regulatory flexibility through alternative splicing to adapt to multiple environments, and that the protein interaction network linking the central proteome is very efficient for synchronizing translation with other biological processes. Surprisingly, at least 10% of the central proteome has no or very limited functional annotation.
Our data and analysis provide a new and deeper description of the human central proteome compared to previous results thereby extending and complementing our knowledge of commonly expressed human proteins. All the data are made publicly available to help other researchers who, for instance, need to compare or link focused datasets to a common background.
The understanding of living cells at a systemic level is being recognized more and more as an important component of biology and medicine research [1–9]. Biological pathways and networks of protein interactions are key paradigms to link molecules to biological functions and by so doing bridging the genotype-to-phenotype gap as well as understanding properties of the organization of biological matter [10–13]. In this work we aim at answering three simple but fundamental questions: i) What is the complement of human proteins expressed ubiquitously and abundantly in different cell types? ii) Does this central proteome (C.Prot)  display properties that are distinct from the rest? iii) Can one identify global features of this central proteome?
Gene expression microarrays allow analyzing a large variety of transcriptomes  and several studies using mRNA detection or abundance as a proxy for protein expression or concentration have revealed important properties of gene sets related to tissue specificity [16–18]. Recently, Bossi and Lehner  showed that tissue-specific proteins are less interacting but bind to core cellular components and common proteins. Domains enriched in tissue-specific genes tend to be metazoan-specific and are non-essential . It is also known that widely expressed genes encode protein domains involved in protein degradation, cytoskeleton or RNA-binding .
It is well known that correlation between transcripts and protein abundance is variable  and, as a general rule of thumb, a good correlation is observed in one third of the observed entities only. Subsequent mechanisms of regulation can significantly decouple protein and transcript abundance . For this reason, we believe that it is important to study the central proteome from proteomics data directly. As our data show, mass spectrometry sensitivity has achieved a level that permits such direct approaches. Similar work was conducted by Schirle, et al. , who first coined the term central proteome and used human cell lines as we did, though they limited their analysis to technical aspects related to the proteomics technology. Kislinger, et al.  profiled protein expression in six mouse organs. Another related project is the Human Protein Atlas  that maps protein expression in human tissues through a selected set of antibodies.
The focus of our work is different compared to the aforementioned transcriptomics and proteomics studies. After a brief and classical analysis of the functions of the proteins present in the central proteome, which matches gene microarray results, we reveal important new findings regarding the gene structures of genes coding the central proteome, location on pathways in relation with drug targets, and global properties of the interaction network connecting the central proteome. Furthermore, we show how several characteristics of common proteins vary with protein abundance.
The large amount of data generated for this research constitutes a unique and homogeneous dataset that should interest other investigators. Data are made available as supplementary material and are accessible from the ProteomeCommons.org Tranche public repository.
Cell lines, proteomics and protein identifications
We measured the proteomes of seven cell lines from the three germ layers (HaCat, HepG2, K562, HEK293, Namalwa, U937, HeLa) with 1D SDS-Page followed by LC-MS/MS. The proteomes contained between 2031 and 4154 proteins each (see Table 1). Protein identification was accomplished by a bioinformatics platform combining two database search engines, Mascot  and Phenyx , and an innovative and very stringent validation strategy enforcing a maximum false discovery rate (FDR) of 0.25% on protein groups . In addition, protein groups that were not made of alternative splice variants exclusively (2%) were discarded. Specific peptides allowed us to ascertain the presence of some variants.
Each cell line was analyzed twice in technical replicates (merged results in Table 1) and modest variability in the identified proteins was observed (<4%).
The central proteome
A large number of proteins were identified in each cell line (Table 1). We constructed the central proteome (C.Prot) by selecting proteins found in all the 7 cell lines, i.e. 1124 proteins. HEK293 cell line yielded notably more protein identifications, which was observed by others already . In addition, as mentioned above, technical replicates of all the cell lines were highly reproducible. HEK293 higher number of proteins is hence unlikely to have been caused by experimental bias. In fact, HEK293 cells provide a convenient system for expressing many proteins, notably in affinity purification MS experiments .
We compared the proteins identified in the 7 cell lines and C.Prot with all the human proteins listed in UniProtKB/Swiss-Prot  to determine possible experimental biases in molecular weight, isoelectric point, hydrophobicity, and aliphatic index ; see Additional file 1: Suppl. Figure 1. We observed very modest shifts towards larger (+2.3% on average) and more acidic proteins (-5.5%) in C.Prot. These slight biases are consequences of the analytical technologies used and we do not believe they have any significant impact on what follows.
Overlap between C.Prot and the Human Protein Altlas (HPA) was of less than 40%, depending on HPA detection strength. Namely, HPA contains all our cell lines but Namalwa and we considered the intersection of the 6 shared cell lines in HPA. We found 852 proteins strongly expressed (178 of which are in C.Prot, 16% of C.Prot) and 3314 at least weakly expressed (413 in C.Prot, 37%). We also considered all (union) the proteins expressed in all the 46 HPA cell lines, weakly or better, and we found 4314 proteins (445 in C.Prot, 40%).
Accession codes, IDs, description, and gene names of the 1124 proteins found in the central proteome are provided as Additional file 2.
Main functions represented in the central proteome
We classified the main categories of proteins present in C.Prot by using a slimmed gene ontology (generic GO slim)  and KEGG [32, 33]. GO term statistics were obtained via topGO  with a 5% cutoff on the P-values ("weight.log" method).
C.Prot was mainly enriched for vital processes of the cell, see Table 2 and Suppl. Figure 2. The GO categories broadly ranged from proteostasis, such as translation and protein transport, over metabolic processes to cell cycle and death. The metabolic processes mainly consisted of primary metabolism, which is vital for maintenance and proliferation. Major catabolic KEGG pathways covered by C.Prot were proteasome, citrate cycle, oxidative phosphorylation, glycolysis/gluconeogenesis, phosphate pentose pathway, fatty acid degradation and few amino acid degradation pathways. On the biosynthetic side we found translation, most aminoacyl-tRNA biosynthesis enzymes and pyrimidine metabolism. The GO category "cellular component organization and biogenesis" contained important complexes such as signal recognition particle, coatomer protein complex and the splicosome. Finally, the cell cycle together with DNA metabolism and cell death, which are vital functions of living in a cellular population, were also enriched. The same broad coverage could be observed for the molecular functions and the cellular localizations in C.Prot, except for the plasma membrane and the extracellular space due to proteomics technical limitations and use of cell lines, see Suppl. Figures 3 and 4.
It is worth noting that advanced GO analysis supported by the topGO R package, where detailed GO annotations are reported to ancestor terms if they do not yield significant results only, combined with generic GO slim helped considerably dealing with a large dataset such as C.Prot. Classical enrichment analysis methods overvalued very general categories and/or returned numerous very detailed hits. We found the "weight.ratio" method of topGO too stringent and the "classic" one too "verbose".
A well conserved set of proteins coded by exon rich genes
The presence of a protein in many species is an indicator of high conservation and fundamental functional role. Therefore, we queried all orthologs of the human genome from Ensembl and we measured conservation by counting, for each human protein, the number of species that had one ortholog. We found that, on average, C.Prot proteins had 5.9 orthologs more than human proteins taken from SwissProt (P = 0, Wilcoxon), that is they were present in more species thereby indicating their higher degree of conservation.
To contrast our results we wanted to compare with cell specific proteins. Due to the limited number of cell lines available it was not feasible to define such specific proteins. We hence used tissue specific genes as a proxy . In the Su et al. dataset, testis genes contributed for almost one third and hence we defined a specific transcriptome (Spe.Trans) ignoring testis specific. As a matter of fact, such diversity in testis gene expression would hide gene specificity in many other tissues otherwise. Without testis genes, Spe.Trans contained 282 genes that, on average, had -0.6 orthologs (P = 0.06) compared to the reference Swiss-Prot database.
An important feature of eukaryote proteins is that they have the possibility to exist in different splice variants. The total count of exons indicates combination diversity. Remarkably, in C.Prot, the exon count was significantly shifted towards higher values with an average shift as large as +3.7 (P = 0, Wilcoxon). Conversely, Spe.Trans had -0.7 exons on average (P = 0.5). To prevent a potential bias in exon counts through a slight bias in protein sequence lengths in our data, we compared exon counts in C.Prot with Swiss-Prot applying a bootstrap procedure to draw appropriate sequence length distributions from Swiss-Prot.
A well connected and central central proteome
Using an integrated interaction database, we extracted a human interactome comprising 9,495 proteins and 70,083 interactions. We computed 5 commonly used connectivity measures (betweenness and eigenvector centrality, clustering coefficient, degree, k-core score) to reflect various local and global aspects of the interactome topology at the positions occupied by C.Prot. All 5 measures indicated significant positive biases, i.e. more centrality, higher connectivity, and more frequent participation in protein complexes. Remarkably, a further significant increase for abundant C.Prot proteins was observed. Node degree and eigenvector centrality statistics appear in Figures 1A and 1B as examples.
There was a significant association of C.Prot with drug targets as listed in DrugBank  (176 targets among the 1465 listed, P = 1.6E-23, bootstrapped χ2 see Methods). In the human interactome, we observed that drug targets were more central and more connected nodes as reported previously by others . We did not see a consistent increase of this trend with targets restricted to C.Prot. On the contrary, considering relative positions in pathways (see Methods), drug targets in C.Prot displayed a significantly different profile compared to targets in general. Targets not restricted to C.Prot appeared to be preferentially at the beginning of pathways as opposed to more central positions (P < 3.0E-2). When restricted to C.Prot, this trend was much augmented (P < 1.4E-34). NCI-PID , the pathway database we used, included 759 DrugBank targets (21%), 133 of which were in C.Prot (12% of C.Prot).
A central interactome
It was natural to define the central interactome as the network made of direct interactions between C.Prot proteins; this network is likely to exist in all the human cells. Shortest path distance distributions were computed (Figure 2A). We also found both the central and whole interactome to be scale-free , with the central interactome comprising more highly-connected regions (protein complexes), see Figure 2B. The central interactome contained several essential protein complexes, see Figure 2C. The central interactome can be regarded as platform used by biological processes to exchange information through protein interactions. Therefore, we introduced a notion of flux between biological processes. As explained in Methods, fluxes between biological processes (BPs) mediated by protein-protein interactions (PPIs) can be scored and we compared fluxes over the central interactome with fluxes outside the central interactome and fluxes between C.Prot and non C.Prot proteins. In each case, we generated random interaction networks and GO annotations to determine which fluxes were significantly more intense than what would be expected by chance from the network topology and GO terms frequencies. At the 1% significance level, the random networks simulation yielded 57 significant GO BP fluxes within the central interactome, 135 between C.Prot and proteins of the human interactome not in more than 5 cell lines, and 365 outside the central interactome. In every case, the expected number was 12.75. See Figure 3.
The central proteome dataset
We present results characterizing the human central proteome (C.Prot), i.e. the set of proteins commonly expressed by human cells. Although previous related studies have been conducted successfully on the basis of transcriptomics data, we based our work on proteomics data. Proteomics is likely to yield additional insight because it directly measures the entities of interest, provided it reaches sufficient sensitivity and does not introduce excessive experimental biases. Moreover, our results comprise several important new findings never covered by transcriptomics studies.
We defined C.Prot to be the proteins shared by 7 cell lines whose total cell lysates were analyzed by state of the art proteomics and followed by very stringent protein bioinformatic identifications; it contains 1124 proteins. Clearly, the intersection of a much larger number of cell lines would be too stringent a criterion but we found it appropriate with 7 cell lines. Experimental biases were modest in our data (Suppl. Figure 1).
We could not find comparable human datasets from public repositories [38, 39] to complement our data and cover more cell lines or tissues. The dataset of Schirle, et al.  is much smaller after mapping to Swiss-Prot (5-10 times) and covered by our data essentially (>92%). The modest 8% that is not covered can certainly been explained by MS detection and sample preparation variability. Also, the few false positives contained in both datasets reduce the overlap further.
Data from the Human Protein Atlas  (HPA) are not appropriate for completing ours as they derive from a biased a priori selection of proteins and antibody availability. Antibody-based assays can also be very variable in their sensitivity, thus making expression profiles of different proteins difficult to compare. Nonetheless, HPA is a very valuable resource and comparison with our data showed that strong HPA common detections in 6 out of the 7 cell lines of our study available in HPA have similar sizes (852 proteins), whereas including weak detections increased this number dramatically to 3314. The overlap with C.Prot is modest in both cases: 178 and 413 proteins respectively (17% and 37%). This shows higher sensitivity detection of targeted antibody-based assays compared to broad unbiased MS analysis, which is no surprise. This also shows that HPA is not covering so far an important proportion of C.Prot.
HPA data can reveal an important characteristic of our data: if we consider all the proteins expressed in all the 46 HPA cell lines weakly or stronger, we find 4314 proteins (445 in C.Prot, 40%), which is a very small improvement compared to 413 above although HPA detections rose from 3314 to 4314. This indicates that, within the limitation of current MS detectability, we already identify with 7 cell lines only a significant part of the MS measurable central proteome. To improve coverage of the central proteome it is more important to improve MS sensitivity than to increase the number of cell lines.
We compared C.Prot with the most abundant proteins detected in any of the 7 cell lines that were not in C.Prot (see Methods). GO annotations were very different for the two datasets (Suppl. Figures 2-4). In fact, the non central abundant proteins had no strong functional association, which indicates that the sole abundance is not sufficient to associate with a function and is no longer a dominating factor for inclusion in the central proteome.
From Su et al.  microarray data we assembled a central transcriptome (C.Trans) for comparison taking the intersection of genes expressed in each tissue (commonly used MAS5 > 200), which resulted in a list of 2002 genes. C.Trans only covered 501 proteins of C.Prot. Su et al. chip contained 917 out of the 1124 proteins of C.Prot, meaning that 45% of the proteins we detected would be missed by a transcriptomics approach. Releasing the criterion for central transcriptome inclusion by requiring genes in all but 4 samples (instead of all), we increased its size to 3197 genes and cover 683 proteins, i.e. a missing part representing 26%. Such losses were most likely due to low degradation rate proteins that do not require permanent expression of their coding genes.
Functional analyses revealed prevalence for proteins involved in proteostasis, RNA-binding, primary metabolism, cell cycle and cell death (Table 2, Suppl. Figures 2 and 3). The transcriptomics data set of specific genes (Spe.Trans)  was associated with different biological processes related to signaling and multicellular organization, whereas vital biological processes were much less prevalent (Suppl. Figure 2). Spe.Trans thus seems to play an important role in the establishment of distinct anatomies via cell-cell signaling, whereas C.Prot covers the more "archaic" needs, what is supported by the higher degree of conservation of its members. These observations corroborate results found from transcriptomics data previously . Proteins expressed across many cell types require augmented gene expression flexibility to adapt to local conditions. Our finding that C.Prot proteins have on average almost 4 additional exons compared to all human proteins shows that evolution preferred this economic way of gaining flexibility in common parts instead of augmenting their number through duplication events. Spe.Trans reduced exon count contrasts nicely.
Kislinger, et al.  did not extract a central proteome from their data. Taking the intersection of their 6 organ datasets and mapping the mouse proteins to their human orthologs, we obtained 393 Swiss-Prot entries, which are also biased towards high exon counts (+3).
Implication in central processes of the cell's life is naturally reflected by augmented PPIs and participation in complexes formation as well as central positions in interaction networks. Common proteins are more connected and central in the human interactome than on average as shown by the 5 different global and local network topology measures we computed. This trend is reversed for tissue specific genes and, in addition, we observed a significant dependence with protein abundance: the more abundant a protein the stronger the connectivity/centrality and vice versa, see Figures 1A and 1B. Such observations might be thought to be artifactual since common proteins are likely to be over-represented in interaction databases and abundantly expressed proteins might be more successful in PPI experiments. To exclude this potential confounding effect we mapped C.Prot down to yeast where large-scale unbiased data are available. Using our integrated database we obtained a yeast interactome comprising datasets published with at least 5000 physical interactions only and confirmed that C.Prot increased connectivity/centrality is still valid in yeast (P < 2.0E-7 for all 5 measures).
Because human interactome data have been gathered in many cell types and conditions, we do not believe that higher connectivity necessarily implies that C.Prot entities are permanently connected to more partners. We rather consider it as a strong indication that these proteins have gained, during evolution, the capability to bind to more partners when needed; another way to augment adaptability to various environments. This is corroborated by the top over-representation of the GO molecular function "protein binding", see Suppl. Figure 3.
To investigate robustness of our findings with respect to a less stringent definition of the central proteome, we checked the exon count and network statistics biases either using the intersection of any six cell lines or taking proteins found in 6 out of 7 cell lines. In every case, the results were almost unchanged.
After considering PPIs, it is natural to move to a higher degree of organization of the living matter, i.e. biological pathways. We define an intuitive notion of relative position along a pathway (0 = source, 1 = end point, see Methods and Suppl. Figure 5) and we observe that C.Prot positions just follow the typical distributions of human proteins (Figure 1C). This is no longer true for the most abundant common proteins, which tend to occupy central positions. This indicates that abundant proteins are more likely at non-rate limiting central positions while the key steps are under tighter expressional control . This is also supported by more uniform presence of Spe.Trans genes and low abundant common proteins along pathways (Figure 1C).
Where do drugs hit the central proteome? The size of the overlap between drug targets and C.Prot is significantly large (176 proteins, P = 1.6E-23). Generally, drug targets are mainly designed against catalytic proteins, transporters and receptors, see Suppl. Figures 3 and 4. Drug targets in C.Prot are mainly enzymes (146), which play a pivotal role in primary metabolic processes (130). The biological processes of amino acid metabolism, precursor generation and carbohydrate metabolism are strongly targeted. A preference against nucleotide binding proteins, e.g. 50 ATP binders, and electron carrier activity processes is observed. Pathway position analysis showed that drug targets are clearly shifted towards sources of pathways, see Figure 1D. This trend was strongly reinforced considering C.Prot drug targets. Figure 4 features GO analysis localized along pathway positions and shows that less central positions are targeted in all the biological processes. We believe that the reason for an additional shift towards initial pathway positions is that C.Prot targets contain metabolic proteins mostly. As a matter of fact, drugs targeting metabolic pathways are frequently designed against rate limiting enzymes, e.g. ATP- or NAD-binding proteins, which are found at the beginning of pathways.
Within the human interactome, C.Prot entities tend to be both closer - in terms of shortest path distance - to themselves and to the non C.Prot proteins than on average in the whole human interactome, Figure 2A. This bias is increased within abundant C.Prot entities. It indicates that C.Prot is embedded in the interactome rather uniformly and does not constitute an isolated island, which certainly increases the robustness of the communications between C.Prot and the rest of the proteome through PPIs.
To better understand internal and external communication within and with C.Prot, we have introduced the central interactome and measured how it synchronizes certain biological processes (BPs) preferentially. Obviously, all the (true) PPIs are biologically relevant but our analysis aimed at identifying the main streams of communication. High-scoring exchanges (fluxes) between BPs provide a summary of the communications for which the interaction network is the most efficient. From Figure 3A we see that most exchanges in the central interactome are used to synchronize BPs with translation. BP communication between C.Prot and non C.Prot and outside the central interactome is no longer specialized and translation has a marginal role.
How well do we know these essential components found in C.Prot? Surprisingly, C.Prot contains 22 proteins lacking any information in PubMed (http://www.ncbi.nlm.nih.gov/), searching abstracts with gene symbols and their synonyms; 73 proteins appear in no more than 3 abstracts and 112 in no more than 6, see Suppl. Figure 6 and Suppl. Table 1.
We have determined experimentally an approximation of the central human proteome that is suitable for analyzing its global properties. It is made of rather well conserved proteins which have gained additional expression flexibility through the acquisition of additional exons. These proteins are mainly involved in proteostasis, primary metabolism, cell cycle and death. They tend to be well connected with other proteins via PPIs. Random network simulations show that the central interactome was made very efficient, through evolution, to coordinate translation with other biological processes, or the latter ones via translation indirectly. More abundant proteins tend to be located at biological pathway central positions. Such central positions are generally devoid of drug targets, especially when they are part of the central proteome, which underline the fundamental role of these proteins whose activity should not be altered. Surprisingly, 10% of these common proteins are essentially uncharacterized. Our data can help other researchers to prioritize protein characterization or serve as background when analyzing focused datasets. They are made publicly available through the journal web site (Additional file 2) and a public repository: the complete lists of protein/peptide identifications in each cell line with spectra have been deposited in ProteomeCommons.org Tranche (hash= "JUrzEy1ShYDDUoKVrxHaoMrAu/CGbqv3xqOS/zuErFvlD8MOrVPRu5kOSlcxwK+/EYdA9WoLN5eMeprBzh9rPMIuYksAAAAAAAAPLQ=="); a more compact version of all the identifications found in each cell line in tabular format and without the spectra has also been deposited in Tranche (hash=" DrsqOg2DmUzUlJVLom+O6AHQTyJa2v+Ekhbw8az6OfF/JfF4hv51cyWqCKmaOHZZnOKJTUDl9ziTdCKpEzirhmjt9csAAAAAAAANKw==").
Total lysates of wild-type K562, HEK293, Namalwa, HaCat, HepG2, U937, HeLa cells (50 μg total protein) were reduced, alkylated, and separated by 1D SDS-PAGE. After visualization of the proteins by Coomassie Blue, entire gel lanes were sliced into 50 equal pieces and digested in situ with modified porcine trypsin . The resultant peptide mixture was extracted from the gel slices and desalted with customized reversed-phase stage tips . Each cell line was grown to complete confluence in appropriate media.
Approximately 10% of each tryptically-digested sample was analyzed as technical replicates by data-dependent nanocapillary reversed-phase LC-MSMS. Peptide separation was via customized 50 μm inner diameter columns packed with 3 μm diameter C18 Reprosil beads coupled to a hybrid LTQ-Orbitrap XL mass spectrometer (ThermoFisher Scientific, Waltham, MA). Data-dependent acquisition was performed for 100 min using one MS channel for every four MSMS channels and a dynamic exclusion for selected ions of 60 s.
Protein identification combined Mascot  and Phenyx , both with 4 ppm/0.3 Da parent/fragment mass tolerance, maximum 1 missed cleavage, carbamidomethyl cystein as fixed modification, and methionine oxidation variable, minimum peptide length 6 amino acids. Searches were performed against UniProtKB/Swiss-Prot (vers. 56.1)  human section, including all the isoforms. Results of the two engines were parsed separately and a minimum of 2 distinct peptides above a score threshold was required. Single peptide hits (SPHs) were also accepted but above a much higher score threshold and provided the protein sequence coverage was 2.5% or more, see hereafter. We discarded spectra assigned to different peptides by the two search engines during the merge of the results. The identified proteins passing this selection were grouped according to shared peptides, and groups with no specific peptides were discarded (1%), Suppl. Figure S7. Based on reverse database searches, we imposed a protein group FDR of 0.25%, see Suppl. Figure S8, which resulted in the following Mascot thresholds: 2 peptides with ion score 18 or more, export all additional peptides with score 10 or more, SPHs with ion score 50 at least; and Phenyx thresholds: 2 peptides with z-score 4.5 or more and P-value 0.001 or less, export all additional peptides with z-score above 3.5, SPHs z-score above 6.
We characterized the robustness of protein groups by determining the peptide false positive (FP) rate induced by the score thresholds we obtained constraining the protein groups FDR. We observed that this FP rate is < 0.1%. Since we discarded protein groups not having any specific peptide (Suppl. Figure 7), the probability that a group is nonetheless artifactual, but not detected as such because of a single FP specific peptide identification, is less than 0.1%.
We have identified more proteins with Phenyx than with Mascot (+4% on average) and with a better sequence coverage: 30% of the proteins have higher sequence coverage with Phenyx, whereas 15% have higher sequence coverage with Mascot. Phenyx scoring function used for this work is available as Additional File 3.
Protein abundance estimation
It is well-known that spectral and peptide counts can provide reasonable estimates of protein abundance [43–45]. We decided to employ emPAI , which is a modified spectral count taking into account the instrument mass range, because it has been carefully validated by its authors and is well accepted. We used a mass range of 698 to 2370 Da that covered 99% of the detected peptides.
To separate the common presence from the sole abundance in the analyses below, we built a comparison dataset Top.CL representing the abundant proteins not commonly expressed. For each protein, we estimated its abundance in the 7 cell lines with emPAI and, to obtain a single number, we took the median. Then, considering the proteins that did not appear in more than 5 cell lines, Top.CL was defined as the proteins in the top 25% median emPAI (463 proteins). To distinguish between different levels of expression within C.Prot, we followed a similar procedure and defined Low.C.Prot (low 25%, 278 proteins) and Top.C.Prot (top 25%, 281 proteins).
Protein interaction data and pathways
We have integrated five public databases containing human protein-protein interactions (PPIs) (MINT , IntAct , HPRD , BioGRID , protein complexes found in NCI-PID ) mapping all accession codes onto Swiss-Prot. From this integrated database, restricting to physical measurements of PPIs by tandem affinity purification (TAP), we extracted 70083 interactions between 9495 distinct human proteins and we built a network model of the human interactome; 859 proteins of C.Prot (76%) were found in this network.
Many measures have been proposed to characterize connectivity in networks . We retained 5 well-established measures ranging from centrality of nodes to likely involvement in complexes and modules: Betweenness centrality , which measures the ratio of the number of shortest paths through a given node by the total number of shortest paths, i.e. how likely it is that a sequence of protein interactions - a pathway - goes through this node; Eigenvector centrality  measures centrality in dependence of neighbors own centrality, i.e. the centrality of a node increases when it is connect to nodes that are central themselves; Node degree is the number of neighbors of a node; Clustering coefficient  measures the ratio between the maximum number of edges theoretically possible in the sub-network spanned by a node and its immediate neighbors and the actual number of edges; k-core score measures the maximum number k for a node to be in a k -core . A k-core is a sub-network with all nodes connected to at least k nodes of the sub-network.
To compare distributions of topological measures, such as the clustering coefficient above, for different sets of proteins, e.g. C.Prot versus C.Trans, we applied a goodness-of-fit test. Using classical χ2 test yielded excessively small P-values since a lot of data points were available and very little differences became significant. To circumvent this problem, we re-sampled the data 1000 times, using 500 data points from each of the two sets to compare, and we computed the χ2 statistic (10 bins). We thus obtained 1000 χ2 statistics, took the median and obtained the P-value from the χ2 distribution.
To investigate positions of proteins in biological pathways, we used the NCI-PID database , which comprises BioCarta (http://www.biocarta.com/) and Reactome  pathways in addition to NCI-PID unique pathways. We found 573 of our proteins (55%) in these pathways and we computed the relative position of each. The relative position was defined considering the shortest paths to the closest source and end nodes, see Suppl. Figure S5 for an example. If a protein appeared several times in a pathway, relative positions were averaged. We did not distinguish between isolated proteins and proteins in complexes, and occurrences of proteins in several pathways were averaged as well. Replacing the shortest path to the closest source/end nodes by the average over all the shortest paths to all the source/end nodes, which would be a reasonable alternative measure, we obtained nearly identical results (data not shown).
To measure fluxes between GO  biological processes (BPs) mediated by PPIs, we counted pairs of GO terms. Namely, given two proteins in interaction, P1 annotated with BPs A and B and P2 annotated with C, D, and E, we count 1 for the pairs AC, AD, AE, BC, BD, and BE. The counts are summed over all PPIs. To avoid "re-discovering" the GO structure, we removed all common ancestors of GO terms found at each protein.
To identify GO fluxes that were stronger than what would be expected as a result of the frequencies of GO terms and protein connectivity in the interactome, we generated 100 random interaction networks and GO annotations. We used 90 such random annotated networks to determine the means and standard deviations of random fluxes between BPs: meanrandom(BPi,BPj), sdrandom(BPi,BPj). Normalized scores were then defined as NScore(BPi,BPj) = (flux(BPi,BPj) - meanrandom(BPi,BPj))/sdrandom(BPi,BPj) and brought all terms to a common scale, independent of their frequencies. The last 10 random networks were used to learn the normalized scores null distribution, which is bell-shaped. Random networks were generated with topologies, GO term frequencies, and central proteome nodes matching the original data. We considered random networks, where node degrees and GO term frequencies were preserved individually but decoupled.
tandem mass spectrometry
physical protein interaction
single peptide hit
exponentially modified protein abundance index
top 25% most abundant proteins not in the central proteome
top 25% most abundant proteins in the central proteome
low 25% abundant proteins in the central proteome
tissue-specific transcripts excluding testis
Human Protein Atlas
Hood L, Heath JR, Phelps ME, Lin B: Systems biology and new technologies enable predictive and preventative medicine. Science. 2004, 306 (5696): 640-643. 10.1126/science.1104635
Gavin AC, Aloy P, Grandi P, Krause R, Boesche M, Marzioch M, Rau C, Jensen LJ, Bastuck S, Dumpelfeld B, et al.: Proteome survey reveals modularity of the yeast cell machinery. Nature. 2006, 440 (7084): 631-636. 10.1038/nature04532
Krogan NJ, Cagney G, Yu H, Zhong G, Guo X, Ignatchenko A, Li J, Pu S, Datta N, Tikuisis AP, et al.: Global landscape of protein complexes in the yeast Saccharomyces cerevisiae. Nature. 2006, 440 (7084): 637-643. 10.1038/nature04670
Kashtan N, Alon U: Spontaneous evolution of modularity and network motifs. Proc Natl Acad Sci USA. 2005, 102 (39): 13773-13778. 10.1073/pnas.0503610102
de Chassey B, Navratil V, Tafforeau L, Hiet MS, Aublin-Gex A, Agaugue S, Meiffren G, Pradezynski F, Faria BF, Chantier T, et al.: Hepatitis C virus infection protein network. Mol Syst Biol. 2008, 4: 230- 10.1038/msb.2008.66
Goh KI, Cusick ME, Valle D, Childs B, Vidal M, Barabasi AL: The human disease network. Proc Natl Acad Sci USA. 2007, 104 (21): 8685-8690. 10.1073/pnas.0701361104
Pujana MA, Han JD, Starita LM, Stevens KN, Tewari M, Ahn JS, Rennert G, Moreno V, Kirchhoff T, Gold B, et al.: Network modeling links breast cancer susceptibility and centrosome dysfunction. Nat Genet. 2007, 39 (11): 1338-1349. 10.1038/ng.2007.2
Yildirim MA, Goh KI, Cusick ME, Barabasi AL, Vidal M: Drug-target network. Nat Biotechnol. 2007, 25 (10): 1119-1126. 10.1038/nbt1338
Church GM: From systems biology to synthetic biology. Mol Syst Biol. 2005, 1 (2005): 0032-
Brehme M, Hantschel O, Colinge J, Kaupe I, Planyavsky M, Kocher T, Mechtler K, Bennett KL, Superti-Furga G: Charting the molecular network of the drug target Bcr-Abl. Proc Natl Acad Sci USA. 2009,
Bergholdt R, Storling ZM, Lage K, Karlberg EO, Olason PI, Aalund M, Nerup J, Brunak S, Workman CT, Pociot F: Integrative analysis for finding genes and networks involved in diabetes and other complex diseases. Genome Biol. 2007, 8 (11): R253- 10.1186/gb-2007-8-11-r253
Lage K, Hansen NT, Karlberg EO, Eklund AC, Roque FS, Donahoe PK, Szallasi Z, Jensen TS, Brunak S: A large-scale analysis of tissue-specific pathology and gene expression of human disease genes and complexes. Proc Natl Acad Sci USA. 2008, 105 (52): 20870-20875. 10.1073/pnas.0810772105
Pawson T, Nash P: Assembly of cell regulatory systems through protein interaction domains. Science. 2003, 300 (5618): 445-452. 10.1126/science.1083653
Schirle M, Heurtier MA, Kuster B: Profiling core proteomes of human cell lines by one-dimensional PAGE and liquid chromatography-tandem mass spectrometry. Mol Cell Proteomics. 2003, 2 (12): 1297-1305. 10.1074/mcp.M300087-MCP200
Su AI, Cooke MP, Ching KA, Hakak Y, Walker JR, Wiltshire T, Orth AP, Vega RG, Sapinoso LM, Moqrich A, et al.: Large-scale analysis of the human and mouse transcriptomes. Proc Natl Acad Sci USA. 2002, 99 (7): 4465-4470. 10.1073/pnas.012025199
Muller FJ, Laurent LC, Kostka D, Ulitsky I, Williams R, Lu C, Park IH, Rao MS, Shamir R, Schwartz PH, et al.: Regulatory networks define phenotypic classes of human stem cell lines. Nature. 2008, 455 (7211): 401-405. 10.1038/nature07213
Lattin JE, Schroder K, Su AI, Walker JR, Zhang J, Wiltshire T, Saijo K, Glass CK, Hume DA, Kellie S, et al.: Expression analysis of G Protein-Coupled Receptors in mouse macrophages. Immunome Res. 2008, 4 (1): 5- 10.1186/1745-7580-4-5
Walker JR, Su AI, Self DW, Hogenesch JB, Lapp H, Maier R, Hoyer D, Bilbe G: Applications of a rat multiple tissue gene expression data set. Genome Res. 2004, 14 (4): 742-749. 10.1101/gr.2161804
Bossi A, Lehner B: Tissue specificity and the human protein interaction network. Mol Syst Biol. 2009, 5: 260- 10.1038/msb.2009.17
Lehner B, Fraser AG: Protein domains enriched in mammalian tissue-specific or widely expressed genes. Trends Genet. 2004, 20 (10): 468-472. 10.1016/j.tig.2004.08.002
Gry M, Rimini R, Stromberg S, Asplund A, Ponten F, Uhlen M, Nilsson P: Correlations between RNA and protein expression profiles in 23 human cell lines. BMC Genomics. 2009, 10: 365- 10.1186/1471-2164-10-365
Lu P, Vogel C, Wang R, Yao X, Marcotte EM: Absolute protein expression profiling estimates the relative contributions of transcriptional and translational regulation. Nat Biotechnol. 2007, 25 (1): 117-124. 10.1038/nbt1270
Kislinger T, Cox B, Kannan A, Chung C, Hu P, Ignatchenko A, Scott MS, Gramolini AO, Morris Q, Hallett MT, et al.: Global survey of organ and organelle protein expression in mouse: combined proteomic and transcriptomic profiling. Cell. 2006, 125 (1): 173-186. 10.1016/j.cell.2006.01.044
Ponten F, Gry M, Fagerberg L, Lundberg E, Asplund A, Berglund L, Oksvold P, Bjorling E, Hober S, Kampf C, et al.: A global view of protein expression in human cells, tissues, and organs. Mol Syst Biol. 2009, 5: 337- 10.1038/msb.2009.93
Perkins DN, Pappin DJ, Creasy DM, Cottrell JS: Probability-based protein identification by searching sequence databases using mass spectrometry data. Electrophoresis. 1999, 20: 3551-3567. 10.1002/(SICI)1522-2683(19991201)20:18<3551::AID-ELPS3551>3.0.CO;2-2
Colinge J, Masselot A, Giron M, Dessingy T, Magnin J: OLAV: towards high-throughput tandem mass spectrometry data identification. Proteomics. 2003, 3 (8): 1454-1463. 10.1002/pmic.200300485
Colinge J, Bennett KL: Introduction to computational proteomics. PLoS Comput Biol. 2007, 3 (7): e114- 10.1371/journal.pcbi.0030114
Glatter T, Wepf A, Aebersold R, Gstaiger M: An integrated workflow for charting the human interaction proteome: insights into the PP2A system. Mol Syst Biol. 2009, 5: 237- 10.1038/msb.2008.75
Wu CH, Apweiler R, Bairoch A, Natale DA, Barker WC, Boeckmann B, Ferro S, Gasteiger E, Huang H, Lopez R: The Universal Protein Resource (UniProt): an expanding universe of protein information. Nucleic Acids Res. 2006, D187-191. 34 Database,
Gasteiger E, Hoogland C, Gattiker A, Duvaud S, Wilkins MR, Appel RD, Bairoch A: Protein Identification and Analysis Tools on the ExPASy Server. The Proteomics Protocols Handbook. Edited by: Walker JM. 2005, 571-607. full_text. Humana Press,
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al.: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25 (1): 25-29. 10.1038/75556
Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, Katayama T, Kawashima S, Okuda S, Tokimatsu T: KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008, D480-484. 36 Database,
Huang da W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4 (1): 44-57. 10.1038/nprot.2008.211
Alexa A, Rahnenfuhrer J, Lengauer T: Improved scoring of functional groups from gene expression data by decorrelating GO graph structure. Bioinformatics. 2006, 22 (13): 1600-1607. 10.1093/bioinformatics/btl140
Wishart DS, Knox C, Guo AC, Cheng D, Shrivastava S, Tzur D, Gautam B, Hassanali M: DrugBank: a knowledgebase for drugs, drug actions and drug targets. Nucleic Acids Res. 2008, D901-906. 36 Database,
Schaefer CF, Anthony K, Krupa S, Buchoff J, Day M, Hannay T, Buetow KH: PID: the Pathway Interaction Database. Nucleic Acids Res. 2009, D674-679. 37 Database,
Jeong H, Tombor B, Albert R, Oltvai ZN, Barabasi AL: The large-scale organization of metabolic networks. Nature. 2000, 407 (6804): 651-654. 10.1038/35036627
Desiere F, Deutsch EW, King NL, Nesvizhskii AI, Mallick P, Eng J, Chen S, Eddes J, Loevenich SN, Aebersold R: The PeptideAtlas project. Nucleic Acids Res. 2006, 655-10.1093/nar/gkj040. 34 Database,
Jones P, Cote RG, Cho SY, Klie S, Martens L, Quinn AF, Thorneycroft D, Hermjakob H: PRIDE: new developments and new datasets. Nucleic Acids Res. 2008, D878-883. 36 Database,
Hackl H, Burkard TR, Sturn A, Rubio R, Schleiffer A, Tian S, Quackenbush J, Eisenhaber F, Trajanoski Z: Molecular processes during fat cell development revealed by gene expression profiling and functional annotation. Genome Biol. 2005, 6 (13): R108- 10.1186/gb-2005-6-13-r108
Shevchenko A, Wilm M, Vorm O, Mann M: Mass spectrometric sequencing of proteins silver-stained polyacrylamide gels. Anal Chem. 1996, 68 (5): 850-858. 10.1021/ac950914h
Rappsilber J, Ishihama Y, Mann M: Stop and go extraction tips for matrix-assisted laser desorption/ionization, nanoelectrospray, and LC/MS sample pretreatment in proteomics. Anal Chem. 2003, 75 (3): 663-670. 10.1021/ac026117i
Old WM, Meyer-Arendt K, Aveline-Wolf L, Pierce KG, Mendoza A, Sevinsky JR, Resing KA, Ahn NG: Comparison of label-free methods for quantifying human proteins by shotgun proteomics. Mol Cell Proteomics. 2005, 4 (10): 1487-1502. 10.1074/mcp.M500084-MCP200
Rappsilber J, Ryder U, Lamond AI, Mann M: Large-scale proteomic analysis of the human spliceosome. Genome Res. 2002, 12 (8): 1231-1245. 10.1101/gr.473902
Ishihama Y, Oda Y, Tabata T, Sato T, Nagasu T, Rappsilber J, Mann M: Exponentially modified protein abundance index (emPAI) for estimation of absolute protein amount in proteomics by the number of sequenced peptides per protein. Mol Cell Proteomics. 2005, 4 (9): 1265-1272. 10.1074/mcp.M500061-MCP200
Chatr-aryamontri A, Ceol A, Palazzi LM, Nardelli G, Schneider MV, Castagnoli L, Cesareni G: MINT: the Molecular INTeraction database. Nucleic Acids Res. 2007, D572-574. 35 Database,
Kerrien S, Alam-Faruque Y, Aranda B, Bancarz I, Bridge A, Derow C, Dimmer E, Feuermann M, Friedrichsen A, Huntley R: IntAct--open source resource for molecular interaction data. Nucleic Acids Res. 2007, D561-565. 35 Database,
Mishra GR, Suresh M, Kumaran K, Kannabiran N, Suresh S, Bala P, Shivakumar K, Anuradha N, Reddy R, Raghavan TM: Human protein reference database--2006 update. Nucleic Acids Res. 2006, D411-414. 34 Database,
Breitkreutz BJ, Stark C, Reguly T, Boucher L, Breitkreutz A, Livstone M, Oughtred R, Lackner DH, Bahler J, Wood V: The BioGRID Interaction Database: 2008 update. Nucleic Acids Res. 2008, D637-640. 36 Database,
Costa LdFR, Francisco A, Travieso , Gonzalo , Boas , Villas PR: Characterization of complex networks: A survey of measurements. Advances in Physics. 2007, 56: 167-242. 10.1080/00018730601170527.
Brandes U: A Faster Algorithm for Betweenness Centrality. Journal of Mathematical Sociology. 2001, 25 (2): 163-177. 10.1080/0022250X.2001.9990249.
Bonacich P: Some unique properties of eigenvector centrality. Social Networks. 2007, 29 (4): 555-564. 10.1016/j.socnet.2007.04.002.
Bader GD, Hogue CW: An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinformatics. 2003, 4: 2- 10.1186/1471-2105-4-2
Tong AH, Drees B, Nardelli G, Bader GD, Brannetti B, Castagnoli L, Evangelista M, Ferracuti S, Nelson B, Paoluzi S, et al.: A combined experimental and computational strategy to define protein interaction networks for peptide recognition modules. Science. 2002, 295 (5553): 321-324. 10.1126/science.1064987
Matthews L, Gopinath G, Gillespie M, Caudy M, Croft D, de Bono B, Garapati P, Hemish J, Hermjakob H, Jassal B: Reactome knowledgebase of human biological pathways and processes. Nucleic Acids Res. 2009, D619-622. 37 Database,
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13 (11): 2498-2504. 10.1101/gr.1239303
We thank all our CeMM colleagues for help and discussions during this research project. We were supported by APP II, APP III, and BIN III GEN-AU grants of the Austrian Ministry of Research and Education.
TRB and JC performed the bioinformatics analysis. TRB, JC and GSF designed the project and wrote the manuscript. MP and KB realized the proteomics analysis of the samples. TB contributed to the project design and the interpretation of the results. IK did the cell culture work. FB implemented and adapted several in-house databases used for this project. All authors read and approved the final manuscript.
Thomas R Burkard and Jacques Colinge contributed equally to this work.
Electronic supplementary material
Additional file 1: Supplementary material. Supplementary material contains several figures and tables that further support the results discussed in the paper. (PDF 1010 KB)
Additional file 2: Central proteome list. A table listing all the accession codes, IDs, description, and gene names of the 1124 proteins found in the central proteome. (XLS 272 KB)
Additional file 3: Linear trap Phenyx scoring function. The scoring function used in this study that we developed for ThermoFisher linear traps. It can be added to any Phenyx installation. (XML 29 KB)
About this article
Cite this article
Burkard, T.R., Planyavsky, M., Kaupe, I. et al. Initial characterization of the human central proteome. BMC Syst Biol 5, 17 (2011). https://doi.org/10.1186/1752-0509-5-17