Skip to main content

What makes species unique? The contribution of proteins with obscure features

Abstract

Background

Proteins with obscure features (POFs), which lack currently defined motifs or domains, represent between 18% and 38% of a typical eukaryotic proteome. To evaluate the contribution of this class of proteins to the diversity of eukaryotes, we performed a comparative analysis of the predicted proteomes derived from 10 different sequenced genomes, including budding and fission yeast, worm, fly, mosquito, Arabidopsis, rice, mouse, rat, and human.

Results

Only 1,650 protein groups were found to be conserved among these proteomes (BLAST E-value threshold of 10-6). Of these, only three were designated as POFs. Surprisingly, we found that, on average, 60% of the POFs identified in these 10 proteomes (44,236 in total) were species specific. In contrast, only 7.5% of the proteins with defined features (PDFs) were species specific (17,554 in total). As a group, POFs appear similar to PDFs in their relative contribution to biological functions, as indicated by their expression, participation in protein-protein interactions and association with mutant phenotypes. However, POF have more predicted disordered structure than PDFs, implying that they may exhibit preferential involvement in species-specific regulatory and signaling networks.

Conclusion

Because the majority of eukaryotic POFs are not well conserved, and by definition do not have defined domains or motifs upon which to formulate a functional working hypothesis, understanding their biochemical and biological functions will require species-specific investigations.

Background

Comparative analysis of eukaryotic genomes provides an unprecedented opportunity to investigate what makes a given species unique. The genetic mechanisms that generate species-specific differences include positively selected mutations that confer a fitness advantage, random fixation of selectively neutral mutations, and acquisition of new genes [1, 2]. Divergence among species, therefore, includes variation in gene sequences, in particular those of regulatory genes, features of non-coding sequences and repetitive DNA, and gene number and repertoire [3, 4]. To date, comparative genomics in eukaryotes has focused largely on genes that encode proteins with experimentally defined domains or motifs (proteins with defined features (PDFs)) [5, 6]. Because the analysis of PDFs revealed a high degree of similarity among different species, it has been accepted widely that the uniqueness of a particular species was driven by changes in regulatory genes or elements [1–6], as opposed to the divergence of established coding sequences or the creation of new genes. This has led to a widespread perspective that just a few model organisms can provide the experimental foundation to assign functions to nearly every eukaryotic gene.

Noticeably lacking from the comparative analysis of eukaryotic genomes to date, however, is an analysis of the origins and functions of genes encoding proteins that currently lack defined motifs or domains (proteins with obscure features (POFs)) [7, 8]. Expression profiling studies in different organisms suggested that POFs play an important role in many different biological processes. Nevertheless, their biological roles and origins remain poorly understood and elucidating their functions is currently a major goal of biological research in almost all organisms studied [7, 8]. In this paper we examine the possibility that genes encoding POFs, which account for approximately one-quarter of all eukaryotic genes, play a role in determining differences among species. By analogy to the expectation that PDFs are often conserved among species [4–6], one might expect that POFs would show a parallel pattern of phylogenetic conservation. To test this assumption, we performed a comparative analysis of 10 different eukaryotic proteomes, including budding and fission yeast, worm, fruit fly, mosquito, Arabidopsis, rice, mouse, rat, and human. Surprisingly, in contrast to PDFs, we found that POFs include a much larger percentage of proteins that are highly divergent. Our results underscore the importance of delineating the origins and functions of POFs as an underlying cause of species specificity.

Results

Approximately one-quarter of eukaryotic proteins are POFs

Proteins were analyzed from ten different model proteomes and classified as POFs if they lacked an established domain or motif including domains of unknown function. The ten model proteomes were derived from gene models based on the genome sequences of Saccharomyces cerevisiae (Sc) and Schizosaccharomyces pombe (Sp), Arabidopsis thaliana (At), Oryza sativa (Os), Drosophila melanogaster (Dm), Anopheles gambiae (Ag), Caenorhabditis elegans (Ce), Mus musculus (Mm), Rattus norvegicus (Rn), and Homo sapiens (Hs). As shown in Figure 1, between 18% and 38% of all proteins (average 26%) predicted from each genome were classified as POFs.

Figure 1
figure 1

Representation of POFs in 10 different eukaryotic genomes. POFs represent 18% to 38% of the proteins in 10 different eukaryotic proteomes (S. cerevisiae, Sc; S. pombe, Sp; A. thaliana, At; O. sativa, Os; D. melanogaster, Dm; A. gambiae, Ag; M. musculus, Mm; R. norvegicus, Rn; H. sapiens, Hs; C. elegans, Ce). Proteomes were obtained and analyzed as described in Materials and methods.

POFs are more divergent than PDFs

To evaluate the diversity among PDFs and POFs in different proteomes, we compared their sequence relatedness to each other using BLAST (Figure 2; Figures 1S to 3S in Additional data file 1). The percentage of related proteins was plotted as a function of similarity cutoff thresholds that ranged from non-stringent (BLAST E-values greater than 10-6) to stringent (from 10-9 to 10-80 or less). This method of plotting similarity differences permits the visualization of reproducible differences between PDFs and POFs across a wide range of cutoff thresholds. Unless noted otherwise, a BLAST similarity of greater than 10-6 was used as the cutoff threshold for classifying a sequence as related. Using this similarity threshold, a total of 1,650 protein groups were found to be conserved among all 10 proteomes (Tables 1S and 2S in Additional data file 1). Surprisingly, only 3 of those (<0.2%) were POFs (as represented in S. cerevisiae by gi|6319274|, gi|6320573| and gi|6324048|).

Figure 2
figure 2

POFs are more divergent than PDFs in different eukaryotic proteomes.(a) Relative similarity among PDFs and POFs in S. cerevisiae (Sc) and S. pombe (Sp). (b) Relative similarity among PDFs and POFs in M. musculus (Mm) and R. norvegicus (Rn). (c) Relative similarity among total or essential PDFs and POFs in S. cerevisiae (Sc) and C. elegans (Ce). BLAST comparisons were performed as described in Materials and methods.

Among the 10 proteomes, POFs always showed significantly more divergence than PDFs, as illustrated for S. cerevisiae and S. pombe, or M. musculus and R. norvegicus in Figure 2a, b, respectively (and documented for all proteomes in Figures 1S to 3S in in Additional data file 1). For example, in a comparison of proteomes from budding and fission yeast (Figure 2a), the percentage of similar POFs was typically five-fold less than PDFs, when evaluated with BLAST cutoff thresholds from 10-6 to 10-18. This difference can also be illustrated by comparing the BLAST values that correspond to the point at which 50% of the PDFs or POFs show relatedness with other proteomes (referred to as the 50% similarity point). In the case of budding and fission yeast, these E-values correspond to approximately 1 for POFs and 10-30 for PDFs. Using a 50% similarity point as a standard for comparison, the POFs from fission and budding yeast show 30 orders of magnitude higher divergence than PDFs. This higher divergence of POFs was also corroborated in comparisons of Sc with the nine other proteomes (Sc versus All, Figure 2a; Figures 1S and 2S in Additional data file 1). The same pattern of POFs divergence was seen in a pair-wise comparison of the two insect proteomes, the two plant proteomes, and the two vertebrate proteomes (Figure 2b; Figures 1S and 2S in Additional data file 1).

This relatively high divergence within the group of POFs was also true in a parallel analysis in which only proteins with essential functions were compared. In Figure 2c, the essential POFs and PDFs were compared from proteomes of Ce [9] and Sc [10]. These two organisms have been subjected to systematic deletion or RNAi analyses to define essential genes. In this comparison, essential POFs were still 3 to 9 orders of magnitude more dissimilar as shown by a comparison between the 50% similarity points.

For each proteome we identified the set of unique proteins not found in any of the other nine proteomes analyzed (Tables 3S and 4S in Additional data file 1). We then determined the relative proportion of POFs or PDFs designated as unique (BLAST cutoff of 10-6). As shown in Figure 3a, the relative percentage of POFs designated as unique was always higher than the percentage of PDFs. On average, we found that 60% of POFs (44,236 in total) are species-specific, in contrast to only 7.5% (17,554 in total) of the PDFs.

Figure 3
figure 3

POFs are more likely to be species specific than PDFs. POFs are more likely to be species specific than PDFs among 10 different proteomes (S. cerevisiae, Sc; S. pombe, Sp; A. thaliana, At; O. sativa, Os; D. melanogaster, Dm; A. gambiae, Ag; M. musculus, Mm; R. norvegicus, Rn; H. sapiens, Hs; C. elegans, Ce). (a) Proportion of POFs and PDFs represented in unique protein sets determined among the 10 different proteomes. Specificity of a protein to a particular proteome was determined based on a BLAST e-value cutoff of 10-6. Numbers on top of bars donate the total number of proteins in each group. (b) Relationship trees among the 10 different proteomes shown in (a). Trees were constructed based on PDFs (dashed blue line) or POFs (red line). Proteome analyses, BLAST comparisons and tree construction were performed as described in Materials and methods. Outl., an outlier E.coli genome was used.

Average sequence similarity relationship trees constructed for all 10 proteomes, based on POFs or PDFs (Figure 3b), revealed that the divergence of POFs among the 10 different proteomes was consistently greater than that of PDFs, supporting the contention that POFs account for the majority of phylogenetically specific ORFs (Figure 3a).

Comparative analysis of the human and chimpanzee proteomes

The recent publication of a draft chimpanzee (Pan troglodytes (Pt)) genome [11] provided us with a unique opportunity to compare the similarity of POFs and PDFs encoded in two proteomes that are estimated to have diverged only 5 to 7 million years ago. Compared to the degree of similarity among POFs from human and mouse, the degree of identity among POFs from human and chimpanzee was much higher (Figure 3S in Additional data file 1). To examine what proportion of human-specific proteins are POFs or PDFs we performed a BLAST search (E-value cutoff of 10-6) of all published sequences against the human proteome. Consistent with POFs representing the majority of species-specific proteins (Figure 3a), POFs accounted for all 27 expressed human-specific proteins (Table 5S in Additional data file 1) not observed in the genomes of any other organism.

Relative contribution of POFs to biological functions

To evaluate POFs for their functional relevance, we first compared the percentage of PDFs and POFs that were represented in expressed sequence tag (EST) collections (Figure 4a). The six animal transcriptomes all showed a greater than 95% representation for transcripts encoding POFs. This high representation was only slightly less than that observed for PDFs. Whereas the two plant transcriptomes showed a somewhat larger difference, the POFs still showed a representation of greater than 75%. Thus, POFs are similar to PDFs in their representation as actively expressed mRNAs.

Figure 4
figure 4

Relative contribution of POFs to biological functions. (a) Representation of POFs in EST libraries from different organisms (A. thaliana, At; O. sativa, Os; D. melanogaster, Dm; A. gambiae, Ag; M. musculus, Mm; R. norvegicus, Rn; C. elegans, Ce; H. sapiens, Hs). (b) Representation of POFs in protein-protein interaction networks in Sc and Ce. (c) Percent phenotypic penetrance of PDFs and POFs in Sc. EST libraries, protein-protein interaction data, insertional mutagenesis data, and sequence comparisons were obtained/performed as described in Materials and methods.

To compare the relative contributions of PDFs and POFs to protein-protein interaction networks, we examined the percentage representation of each protein class in global interaction data sets available for the Sc [12] and Ce [13] proteomes (Figure 4b). While the POFs showed a slightly reduced representation compared to PDFs, the relative differences were less than 7%.

To compare the relative phenotypic contribution of both protein groups, we examined the percent representation of corresponding mutant phenotypes from the genome-wide functional analyses conducted for S. cerevisiae (Sc) [10] (Figure 4c). With the exception of a potentially noteworthy two-fold lower representation of POFs in the essential gene categories, PDFs and POFs showed a similar percent contribution to other phenotypic categories. Similar results were found with C. elegans genome-wide functional analysis (Ce) [9] (data not shown).

Together, the findings presented in Figure 4 suggest that POFs, as a group, are not being mis-represented by an unusually high percentage of proteins being incorrectly predicted from inaccurate gene models. Rather, POFs appear comparable to PDFs in their relative contribution to an organism's repertoire of functional proteins.

POFs are typically shorter and contain a higher content of disordered structure

To investigate if there are any structural characteristics other than established motifs and domains that might distinguish POFs from PDFs, the physical properties of the two groups of proteins were examined. Compared to PDFs, POFs as a group are 40% shorter (Figure 5a; ANOVA p < 0.001), have a higher percentage of disordered structure (Figure 5b; ANOVA p < 0.001), a higher content of hydrophilic residues (Figure 5c; ANOVA p < 0.001), a higher content of small amino acids (for example, proline and serine), glutamine, and arginine and a lower content of aliphatic (for example, isoleucine and valine) and aromatic (tyrosine) amino acids, and aspartic acid (Table 6S in Additional data file 1; ANOVA p < 0.01). Therefore, in addition to the absence of established motifs and domains, POFs, on average, have physical characteristics that further distinguish them from PDFs.

Figure 5
figure 5

Distinct differences in several biophysical characteristics between POFs and PDFs. A comparison of PDFs and POFs shows distinct differences in several biophysical characteristics (S. cerevisiae, Sc; S. pombe, Sp; A. thaliana, At; O. sativa, Os; D. melanogaster, Dm; A. gambiae, Ag; M. musculus, Mm; R. norvegicus, Rn; H. sapiens, Hs; C. elegans, Ce). (a) Average length of PDFs and POFs from different species. (b) Structural disorder index (disorder/length) of PDFs and POFs from different species. (c) Hydrophilic index of PDFs and POFs from different species. Proteomes analyses were performed as described in Materials and methods.

Discussion

Assigning a role to proteins with unknown function is a major goal of current and future genomic research. Homology searches have traditionally been used to assign specific domain structures to proteins, thereby classifying them into protein families with putative function [5, 6]. The distinction made in this paper between proteins with defined motifs or domains (PDFs), and those with undefined or obscure features (POFs) (Figure 1), underlined proteins that could not be assigned any known function by homology searches. Interestingly, POFs as a group were found to be shorter, more hydrophilic and more disordered than PDFs (Figure 5).

Our analysis of 10 different proteomes (Figures 1 to 3; Figures 1s to 3s in Additional data file 1) revealed a striking difference in conservation between PDFs and POFs. E-value plots show that this difference is evident whether similarities are measured with stringent or non-stringent criteria. With a minimum E-value similarity threshold of 10-6, a total of 44,236 phylogenetically specific POFs were identified (Figure 3a). In contrast, relatively few (17,544) PDFs were identified as phylogenetically specific. The opposite trend was observed for conserved proteins. Only 3 POF groups were conserved in all 10 proteomes, in contrast to 1,650 PDF groups (Tables 1S and 2S in Additional data file 1). In total, 60% of POFs appear to be phylogenetically restricted, in contrast to only 7.5% of PDFs.

One explanation as to why POFs show a higher degree of species specificity than PDFs is that POFs, in contrast to PDFs, could include a disproportionately higher number of proteins incorrectly predicted from pseudogenes or incorrect gene models. This could result in an artifact in which random sequences or non-functional proteins distort the overall diversity of this group of proteins. However, several lines of evidence presented here (Figures 2, 4 and 5) suggest that this trivial explanation is not the case. For example, even in the relatively well-characterized yeast genomes, the pattern of higher sequence diversification among the POFs is consistent with that seen in the larger more complex genomes. Furthermore, as a group, POFs appear similar to PDFs with respect to mRNA expression, phenotypic penetrance, and involvement in protein-protein interactions.

An additional and noteworthy characteristic supporting the contention that most POFs contribute functional activities is that, on average, POFs show a greater degree of predicted disordered structure (Figure 5). Empirical definitions of disordered structures have been derived from examining regions of proteins that fail to show a consistent or defined structure in a crystallized protein. These regions of disorder show a strong correlation with biochemical studies that suggest their involvement in protein-protein interactions, as well as in providing key regions for regulating a protein's activity via a structural conformation switch [14–18]. Importantly, the disorder prediction software programs used here do not provide high scores to random 'junk' DNA sequences, providing another line of evidence that gene models encoding POFs were not just artifactual predictions derived from 'junk' DNA. Rather, the high levels of disordered structures in POFs support their potential roles in regulatory networks in which protein conformational changes or protein-protein interactions are key.

Together, the above arguments strongly support the view that the average POF is just as likely to have a biological function as a protein with a defined motif or domain. We favor, therefore, a genetic explanation for the unusual diversity within the group of POFs. That explanation is that genes encoding the majority of POFs are arising de novo or diverging at an evolutionary rate much higher than genes encoding PDFs. In support of this possibility, POFs are consistently more divergent among different proteomes (Figure 3b), and are preferentially represented as singletons in the different genomes (Figure 4S in Additional data file 1).

There may be several distinct mechanisms contributing to the genetic diversity of POFs. For example, some POFs may have structures that are highly flexible and can diverge with few structural constraints. By contrast, some POFs may have conserved tertiary structures, but are nevertheless showing rapid divergence in their primary sequence. In either case, a widely conserved motif or distinct domain signature may never be found within the primary sequences of a large subset of the currently defined POFs. Nevertheless, some POFs may ultimately be found to have definable features. One reason that these features currently remain undefined may be related to the sociology of science. In general, scientists have focused their molecular research on relatively few organisms, and devoted most of their resources to in-depth studies of relatively few proteins. Those proteins or pathways are often chosen because of their general relevance to fundamental questions in a broad group of organisms or because these proteins exhibit strong evolutionary conservation and are judged on this basis to have greater intrinsic functional relevance. By contrast, the study of a species-specific protein is often a lonely pursuit. Another source of bias lies in the tendency inherent in classical biochemical methods, which are strongly biased towards the production and characterization of folded, active proteins that have highly ordered structures (for example, PDFs) and for which structural information is more readily obtained [19]. In contrast, disordered proteins (for example, POFs) are less well studied because they lack a readily recognized activity and structural information is more difficult to obtain for these proteins.

Previous work identified a class of proteins termed 'ORFans' that have no significant sequence similarity to any other open reading frame (ORF) and are, therefore, unique to a specific organism [20–22]. In contrast to the definition of POFs that is based on the presence of an observed domain or motif, the definition of an ORFan is based strictly on sequence homology. Thus, ORFans could include POFs as well as PDFs. Indeed, as shown in Figure 3a, POFs and PDFs accounted for 70.4% (42,218) and 29.6% (17,544) of all proteins unique among the 10 analyzed proteomes, respectively. Moreover, the majority of POFs from Mm and Rn were found to be similar (Figure 3a), suggesting that although some overlap exists between ORFans and POFs, homologs of many POFs can be found in similar genomes (Figure 3a). An interesting observation that was recently made for ORFans could also hold true for POFs. It was observed that some ORFans, although demonstrating no sequence homology to any known protein, could fold into a three-dimensional structure that resembled a protein with a known function [22]. In addition to being novel genes unique to an organism or a lineage, some POFs or ORFans could, therefore, be the result of convergent evolution. Thus, they might be distant members of known proteins, with similar functions and three-dimensional structure, but with sequences that have diverged beyond recognition.

Conclusion

The advent of genome sequences has reinvigorated an effort to understand the origins of species specificity. This is a daunting challenge, emphasized by the fact that in the 10 proteomes analyzed here we identified 44,236 phylogenetically specific proteins with undefined or obscure features (POFs). In contrast to PDFs, which have established domains or motifs that can be used to formulate working hypotheses about a protein's function, advancing our understanding of POFs must proceed without such clues. Our analysis here provides an expectation that, on average, 60% of a eukaryote's set of POFs will be highly divergent, and that functional studies will ultimately need to be conducted on a species-specific basis. For example, the human genome encodes 27 proteins that currently cannot be found in genome sequences of any model organism, including the chimpanzee sequence (Table 5S in Additional data file 1). Consistent with expectations from this study, these human-specific proteins are all POFs. Eventually, the function of these unique proteins will need to be studied in humans. Our results support a general expectation that to understand the unique biology of a given organism will ultimately involve understanding the functions of an unexpectedly large number of proteins that have: no defined motifs or domains; are likely to have significant regions of disordered structure; and are restricted to a single species or a closely related phylogenetic branch.

Materials and methods

Protein data and definition of POFs

Gene models and proteomes for Saccharomyces cerevisiae (Sc), Schizosaccharomyces pombe (Sp), Arabidopsis thaliana (At), Oryza sativa (Os), Drosophila melanogaster (Dm), Anopheles gambiae (Ag), Caenorhabditis elegans (Ce), Mus musculus (Mm), Rattus norvegicus (Rn), and Homo sapiens (Hs) were downloaded from the NCBI website [23] and from TAIR [24] on 5 December, 2004. Gene models and proteomes for Pan troglodytes (Pt),Mus musculus (Mm), Rattus norvegicus (Rn), and Homo sapiens (Hs) were also downloaded from Ensembl [25] on 10 September, 2005.

To standardize the classification for which proteins are POFs and which are PDFs, we applied a consistent analysis method to all genes regardless of their current annotation. This analysis method involved an HMMPFAM [26] search against several major signature databases: PFAM [27], TIGRFAM [28], SMART [29], and Superfamily [30]. A protein sequence with a match to one or more of the models in any one of these databases, including domains of unknown function, was flagged as a PDF. Sequences with no matches to any one of the models in any database were flagged as a POF. The definition of POFs used in our work was similar to that used in [9, 31].

BLAST comparisons

BLAST comparisons of PDFs and POFs among different proteomes were performed using TeraBLAST running on an accelerated DeCypher server [32]. The comparisons of PDFs and POFs between each proteome and its respective collection of ESTs or between PDFs and POFs from each proteome and all other genomes translated in all reading frames were accomplished using TBLASTn [33]. ESTs were obtained on 5 December, 2004 from NCBI [23] except for Arabidopsis, which was downloaded at the same time from TAIR [24]. To examine the representation of PDFs and POFs from Sc or Ce in existing phenotypic studies, or existing protein-protein interaction datasets, POFs and PDFs, obtained as described above, were matched to existing datasets [9, 10, 12, 13, 34].

Prediction of protein properties

Prediction of relative disorder for PDFs and POFs was performed with the DisEMBL 1.4 prediction program [35]. Due to the large numbers of proteins that were analyzed, we used DisEMBL locally rather than at the website [36]. To obtain an overall value for the percentage of proteins that were disordered, SAS V9.1 (SAS Institute Inc., Cary, NC, USA) was used to sum the total regions that were predicted to be disordered and to divide it by the length for each protein analyzed. Hydrophilic index and amino acid content were calculated with the ad hoc perl script hydrophil.pl, developed by Garay-Arroyo et al. [37] using the Kyte-Doolittle values for hydrophilicity. SAS was used to perform statistical analysis (descriptive statistics and ANOVA) for the hydrophilic index of POFs and PDFs from hydrophil.pl results (Figure 5c), for variations in amino acid content between POFs and PDFs (Table 6S in Additional data file 1), and for sequence length and relative disorder of POFs and PDFs (Figure 5a, b). Because the average length of the POFs was shorter than that of the PDFs, a length correction was used to eliminate bias in the scoring function of the program.

'All-against-all' comparisons and tree generation

'All-against-all' comparisons used to generate sets of species-specific proteins for both PDFs and POFs from Sc, Sp, At, Os, Dm, Ag, Ce, Mm, Rn and Hs were performed using TeraBLAST running on an accelerated DeCypher server [32], with a cutoff threshold of 10-6. A tree showing the relationships among Sc, Sp, At, Os, Dm, Ag, Ce, Mm, Rn and Hs proteomes was constructed using the reciprocal percentage of the number of genes that the organisms have in common. This tree was constructed using the SAS cluster procedure utilizing the average linkage method, and graphed using the SAS tree procedure [38]. Tree diagrams are discussed in the context of cluster analysis by Hartigan [39], and Everitt [40].

Additional data files

The following additional data are available with the online version of this paper. Additional data file 1 contains supplemental figures and tables. Supplemental Figures 1-1 through 1-10 show the relative similarity among PDFs and POFs in all proteomes studied. Supplemental Figure 2 shows the relative similarity among PDFs and POFs in selected proteomes measured as percentage identity or percentage similarity. Supplemental Figure 3 shows the relative similarity among PDFs and POFs between Hs and Pt, compared to Hs and Mm. Supplemental Figure 4 shows cluster analysis of POFs and PDFs in selected proteomes. Supplemental Table 1 lists common POFs to all proteomes analyzed. Supplemental Table 2 lists common PDFs to all proteomes analyzed. Supplemental Table 3 lists unique PDFs from all proteomes analyzed. Supplemental Table 4 lists unique POFs from all proteomes analyzed. Supplemental Table 5 lists 27 unique Hs proteins with representation in EST databases. Supplemental Table 6 describes the amino acid content of POFs and PDFs from the different proteomes studied.

References

  1. Fay JC, Wu CI: Sequence divergence, functional constraint, and selection in protein evolution. Annu Rev Genomics Hum Genet. 2003, 4: 213-235. 10.1146/annurev.genom.4.020303.162528.

    Article  PubMed  CAS  Google Scholar 

  2. Yang Z: Inference of selection from multiple species alignments. Curr Opin Genet Dev. 2002, 12: 688-694. 10.1016/S0959-437X(02)00348-9.

    Article  PubMed  CAS  Google Scholar 

  3. Robichaux RH, Purugganan MD: Accelerated regulatory gene evolution in an adaptive radiation. Proc Natl Acad Sci USA. 2001, 98: 10208-10213. 10.1073/pnas.181257698.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Liti , Louis EJ: Yeast evolution and comparative genomics. Annu Rev Microbiol. 2005, 59: 135-153. 10.1146/annurev.micro.59.030804.121400.

    Article  PubMed  CAS  Google Scholar 

  5. Orengo CA, Thornton JM: Protein families and their evolution-a structural perspective. Annu Rev Biochem. 2005, 74: 867-900. 10.1146/annurev.biochem.74.082803.133029.

    Article  PubMed  CAS  Google Scholar 

  6. Marsden RL, Lee D, Maibaum M, Yeats C, Orengo CA: Comprehensive genome analysis of 203 genomes provides structural genomics with new insights into protein family space. Nucleic Acids Res. 2006, 34: 1066-1080. 10.1093/nar/gkj494.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  7. Roberts RJ: Identifying protein function - a call for community action. PLoS Biol. 2004, 2: E42-10.1371/journal.pbio.0020042.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Galperin MY, Koonin EV: 'Conserved hypothetical' proteins: prioritization of targets for experimental study. Nucleic Acids Res. 2004, 32: 5452-5463. 10.1093/nar/gkh885.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  9. Kamath RS, Fraser AG, Dong Y, Poulin G, Durbin R, Gotta M, Kanapin A, Le Bot N, Moreno S, Sohrmann M, et al: Systematic functional analysis of the Caenorhabditis elegans genome using RNAi. Nature. 2003, 421: 231-237. 10.1038/nature01278.

    Article  PubMed  CAS  Google Scholar 

  10. Giaever G, Chu AM, Ni L, Connelly C, Riles L, Veronneau S, Dow S, Lucau-Danila A, Anderson K, Andre B, et al: Functional profiling of the Saccharomyces cerevisiae genome. Nature. 2002, 418: 387-391. 10.1038/nature00935.

    Article  PubMed  CAS  Google Scholar 

  11. The Chimpanzee Sequencing and Analysis Consortium: Initial sequence of the chimpanzee genome and comparison with the human genome. Nature. 2005, 437: 69-87. 10.1038/nature04072.

    Article  Google Scholar 

  12. Tong AH, Lesage G, Bader GD, Ding H, Xu H, Xin X, Young J, Berriz GF, Brost RL, Chang M, et al: Global mapping of the yeast genetic interaction network. Science. 2004, 303: 808-813. 10.1126/science.1091317.

    Article  PubMed  CAS  Google Scholar 

  13. Li S, Armstrong CM, Bertin N, Ge H, Milstein S, Boxem M, Vidalain PO, Han JD, Chesneau A, Hao T, et al: A map of theinteractome network of the metazoan C. elegans. Science. 2004, 303: 540-543. 10.1126/science.1091403.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  14. Dunker AK, Lawson JD, Brown CJ, Williams RM, Romero P, Oh JS, Oldfield CJ, Campen AM, Ratliff CM, Hipps KW, et al: Intrinsically disordered protein. J Mol Graph Model. 2001, 19: 26-59. 10.1016/S1093-3263(00)00138-8.

    Article  PubMed  CAS  Google Scholar 

  15. Tompa P, Csermely P: The role of structural disorder in the function of RNA and protein chaperones. FASEB J. 2004, 18: 1169-1175. 10.1096/fj.04-1584rev.

    Article  PubMed  CAS  Google Scholar 

  16. Dunker AK, Brown CJ, Lawson JD, Iakoucheva LM, Obradovic Z: Intrinsic disorder and protein function. Biochemistry. 2002, 41: 6573-6582. 10.1021/bi012159+.

    Article  PubMed  CAS  Google Scholar 

  17. Gunasekaran K, Tsai CJ, Kumar S, Zanuy D, Nussinov R: Extended disordered proteins: targeting function with less scaffold. Trends Biochem Sci. 2003, 28: 81-85. 10.1016/S0968-0004(03)00003-3.

    Article  PubMed  CAS  Google Scholar 

  18. Brown CJ, Takayama S, Campen AM, Vise P, Marshall TW, Oldfield CJ, Williams CJ, Dunker AK: Evolutionary rate heterogeneity in proteins with long disordered regions. J Mol Evol. 2002, 55: 104-110. 10.1007/s00239-001-2309-6.

    Article  PubMed  CAS  Google Scholar 

  19. Dyson HJ, Wright PE: Intrinsically unstructured proteins and their functions. Nat Rev Mol Cell Biol. 2005, 6: 197-208. 10.1038/nrm1589.

    Article  PubMed  CAS  Google Scholar 

  20. Fischer D, Eisenberg D: Finding families for genomic ORFans. Bioinformatics. 1999, 15: 759-762. 10.1093/bioinformatics/15.9.759.

    Article  PubMed  CAS  Google Scholar 

  21. Siew N, Fischer D: Analysis of singleton ORFans in fully sequenced microbial genomes. Proteins. 2003, 53: 241-251. 10.1002/prot.10423.

    Article  PubMed  CAS  Google Scholar 

  22. Siew N, Fischer D: Structural biology sheds light on the puzzle of genomic ORFans. J Mol Biol. 2004, 342: 369-373. 10.1016/j.jmb.2004.06.073.

    Article  PubMed  CAS  Google Scholar 

  23. NCBI Index. [ftp://ftp.ncbi.nlm.nih.gov/]

  24. The Arabidopsis Information Resource. [http://www.arabidopsis.org/]

  25. Ensembl. [http://www.ensembl.org]

  26. HMMER. [http://hmmer.wustl.edu/]

  27. Pfam. [http://www.sanger.ac.uk/Software/Pfam/]

  28. TIGR Protein Families. [http://www.tigr.org/TIGRFAMs/]

  29. SMART. [http://smart.embl-heidelberg.de/]

  30. SUPERFAMILY: Main page. [http://supfam.org/SUPERFAMILY/]

  31. Chothia C, Gough J, Vogel C, Teichmann SA: Evolution of the protein repertoire. Science. 2003, 300: 1701-1703. 10.1126/science.1085371.

    Article  PubMed  CAS  Google Scholar 

  32. TimeLogic. [http://www.timelogic.com]

  33. NCBI BLAST. [http://www.ncbi.nlm.nih.gov/blast/]

  34. Database of Interacting Proteins. [http://dip.doe-mbi.ucla.edu/]

  35. Linding R, Jensen LJ, Diella F, Bork P, Gibson TJ, Russell RB: Protein disorder prediction: implications for structural proteomics. Structure. 2003, 11: 1453-1459. 10.1016/j.str.2003.10.002.

    Article  PubMed  CAS  Google Scholar 

  36. DisEMBL. [http://dis.embl.de]

  37. Garay-Arroyo A, Colmenero-Flores JM, Garciarrubio A, Covarrubias A: Highly hydrophilic proteins in prokaryotes and eukaryotes are common during conditions of water deficit. J Biol Chem. 2000, 275: 5668-5674. 10.1074/jbc.275.8.5668.

    Article  PubMed  CAS  Google Scholar 

  38. SAS 9 Documentation. [http://support.sas.com/documentation/onlinedoc/sas9doc.html]

  39. Hartigan J: Clustering Algorithms. 1975, New York, USA: Wiley

    Google Scholar 

  40. Everitt BS: Cluster Analysis. 1998, London, UK: Edward Arnold

    Google Scholar 

Download references

Acknowledgements

This work was supported in part by grants from the National Science Foundation (2010 Program IBN-0420033 and IBN-0420152) and the Nevada Agricultural Experiment Station, publication no. 03066915. MG acknowledges support from the NIH IDeA Network of Biomedical Research Excellence (INBRE, RR-03-008). The Nevada Bioinformatic Center acknowledges support from the NIH-NCRR Biomedical Research Infrastructure Network (P20 RR016464) and NSF EPSCoR (EPS-0132556) Integrated Approaches to Abiotic Stress Cluster.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Ron Mittler.

Electronic supplementary material

13059_2006_1325_MOESM1_ESM.pdf

Additional data file 1: Supplemental Figures 1-1 through 1-10 show the relative similarity among PDFs and POFs in all proteomes studied. (PDF 903 KB)

Authors’ original submitted files for images

Rights and permissions

Reprints and permissions

About this article

Cite this article

Gollery, M., Harper, J., Cushman, J. et al. What makes species unique? The contribution of proteins with obscure features. Genome Biol 7, R57 (2006). https://doi.org/10.1186/gb-2006-7-7-r57

Download citation

  • Received:

  • Revised:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/gb-2006-7-7-r57

Keywords