Email updates

Keep up to date with the latest news and content from BMC Genomics and BioMed Central.

This article is part of the supplement: SNP-SIG 2012: Identification and annotation of SNPs in the context of structure, function, and disease

Open Access Research

A protein domain-centric approach for the comparative analysis of human and yeast phenotypically relevant mutations

Thomas A Peterson1, DoHwan Park2 and Maricel G Kann1*

Author Affiliations

1 Department of Biological Sciences, University of Maryland, Baltimore County, Baltimore, MD, USA

2 Department of Mathematics and Statistics, University of Maryland, Baltimore County, Baltimore, MD, USA

For all author emails, please log on.

BMC Genomics 2013, 14(Suppl 3):S5  doi:10.1186/1471-2164-14-S3-S5


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/bmcgenomics/1471-2164/14/S3/S5


Published:28 May 2013

© 2013 Peterson et al.; licensee 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.

Abstract

Background

The body of disease mutations with known phenotypic relevance continues to increase and is expected to do so even faster with the advent of new experimental techniques such as whole-genome sequencing coupled with disease association studies. However, genomic association studies are limited by the molecular complexity of the phenotype being studied and the population size needed to have adequate statistical power. One way to circumvent this problem, which is critical for the study of rare diseases, is to study the molecular patterns emerging from functional studies of existing disease mutations. Current gene-centric analyses to study mutations in coding regions are limited by their inability to account for the functional modularity of the protein. Previous studies of the functional patterns of known human disease mutations have shown a significant tendency to cluster at protein domain positions, namely position-based domain hotspots of disease mutations. However, the limited number of known disease mutations remains the main factor hindering the advancement of mutation studies at a functional level. In this paper, we address this problem by incorporating mutations known to be disruptive of phenotypes in other species. Focusing on two evolutionarily distant organisms, human and yeast, we describe the first inter-species analysis of mutations of phenotypic relevance at the protein domain level.

Results

The results of this analysis reveal that phenotypic mutations from yeast cluster at specific positions on protein domains, a characteristic previously revealed to be displayed by human disease mutations. We found over one hundred domain hotspots in yeast with approximately 50% in the exact same domain position as known human disease mutations.

Conclusions

We describe an analysis using protein domains as a framework for transferring functional information by studying domain hotspots in human and yeast and relating phenotypic changes in yeast to diseases in human. This first-of-a-kind study of phenotypically relevant yeast mutations in relation to human disease mutations demonstrates the utility of a multi-species analysis for advancing the understanding of the relationship between genetic mutations and phenotypic changes at the organismal level.

Background

The study of human genomic variations, in particular those in protein coding regions, can lead to new hypotheses about the molecular mechanisms of human diseases and might provide critical knowledge about individual response to therapy [1,2]. The advent of large-scale experimental techniques is providing new phenotypic associations for genomic variations [3-5]. However, genomic association studies are limited by the molecular complexity of the phenotype being studied and the cohort size needed to have adequate statistical power. One way to circumvent this problem, which is critical in the study of rare diseases, is to investigate the molecular patterns emerging from functional studies of existing disease mutations. In current large association studies, such as GWAS or upcoming whole-exome and whole-genome sequencing, this is accomplished by aggregating mutations that disrupt the same gene [6,7], pathway [8], or network [9]. In many cases, these molecular variations associated with human diseases have patterns that are similar to those producing a phenotypic change in other species. For example, the comparison between close species has made significant contributions to the biomedical field, such as the use of mice [10] and rats [11] for genetics and drug discovery. In addition, studies across species with longer evolutionary distances to human have many advantages and could bring new perspectives into the study of molecular mechanisms of human phenotypic variations. For instance, the functional analysis of variations in yeast, an organism that can be easily genetically manipulated, has shed light on variations in their human gene orthologs, as shown in McGary et al. [12]. The authors demonstrated the potential of a systematic study of phenotypes produced by variations in human and their orthologs in yeast or other distantly related species, providing novel hypotheses about human diseases, which have already resulted in valuable leads for drug discovery.

The vast majority of studies related to human disease mutations are performed by comparison of whole proteins, which here will be denoted by the genes that encode them. However, these whole-protein approaches are of limited applicability to the study of disease mutations due to the fact that they mostly fail to account for protein modularity. Most proteins contain multiple domains that can be recombined in different arrangements to create proteins with different functions [13-15]. As a consequence, not all protein regions have the same function or produce similar phenotypic changes if disrupted. Thus, the specific location of a particular mutation within the protein could be crucial to understanding the mutation's functional effect. The relevance of studying protein domains in the context of disease was also discussed by Zhong et al. [16] in their study of protein interactions and their relation to diseases. The authors showed that mutations resulting in complete loss of the protein product (removal of a node in the network) could be different from those disrupting only a protein region or domain (edgetic perturbations). Furthermore, Zhong et al. conclude that these edgetic perturbations can cause clinically distinct phenotypes when disrupting different protein domain regions of the same protein. Thus, a domain-centric study of disease mutations has the potential to differentiate among genomic variations by accounting for protein modularity that would have otherwise been grouped together by whole-protein studies.

To capture the disruption of domains by genetic mutations, we have previously created a database to visualize the aggregation patterns of disease mutations at the protein and domain levels for human genomics data (Domain Mapping of Disease Mutations database (DMDM), freely available at http://bioinf.umbc.edu/dmdm/ webcite) [17]. More recently, we have developed a statistical approach, the domain significance score (or DS-Score), for finding significantly mutated positions for individual protein domains [18]. We demonstrated that significant DS-Scores indicate that a mutation at a specific position is highly likely to be a contributor to disease in any protein containing the domain in which the mutations are located. In particular, we have shown that Mendelian disease mutations form clusters at protein domain sites [18]. In addition, results from Yue et al. [19], Nehrt et al. [20], and Peterson et al. [18] have further shown that inherited and somatic cancer mutations cluster at specific sites at the protein domain level. Thus, these studies show how the domain analysis enables the discovery of domain hotspots of mutations with phenotypic relevance by aggregating mutations that share the same domain location but are localized in different genes. However, the discovery of these highly deleterious domain sites by aggregation of mutational data with known phenotypic effect is limited by the availability of such mutational data. As a result, the DS-Score method based on human data has low coverage when analyzing mutations from large-scale sequencing studies. To address this issue, more annotated disease mutation data will need to be incorporated into the analysis, preferably from other species in which the phenotypic effect of putative deleterious mutations could be experimentally tested.

In this paper, we describe the first inter-species analysis of mutations of phenotypic relevance at the protein domain level for human and yeast genomes. We perform the comparison between these species by mapping human and yeast mutations into the corresponding domain sites. Protein domains, such as those defined by CDD [21] and Pfam [22], are protein sequence regions that are highly conserved across distantly related species. For instance, when comparing yeast and human domains, we estimate that 87% of all the protein domains found in yeast are also found in human while, using the Homologene database to compare genes, only 20% of the yeast genes have a human ortholog [23]. Similarly, 58% of the human domains are shared with yeast while only 5% of the human genes have yeast orthologs. Since yeast and human analyses show a significant number of common domains, the protein domain framework facilitates the comparison of a significant number of mutations producing phenotypic changes in both species. Using a domain-centric approach, we show that phenotypically relevant mutations in yeast form hotspots at the protein domain level, and that a significant number of these hotspots map to known human disease mutations. Furthermore, our results show that the feature-based DS-Score, a modification of our statistical method that explicitly incorporates annotation from the protein domain models, was most successful at capturing functional commonalities between human and yeast mutations affecting these organisms' phenotypes.

In summary, the work described in this paper demonstrates that domain-centric, inter-species mutation analyses lead to the identification of new domain sites of relevance to human diseases even when performed among species separated by long evolutionary distances. The patterns of evolutionarily conserved and functional mutations associated with phenotypic changes emerging from this study represent a step towards a new paradigm for the analysis of large-scale genomic studies of human diseases.

Materials and methods

Databases

A human protein database containing 54,372 proteins was created with 33,963 proteins from RefSeq [24] and 20,409 proteins from Swiss-Prot [25] downloaded via NCBI's E-utilities [23]. Since the RefSeq and Swiss-Prot databases contain many redundant protein entries, we selected only one representative protein for each unique Entrez gene ID, either the longest Swiss-Prot protein, or the longest RefSeq protein if no Swiss-Prot protein was listed for the gene ID. A database of 6,717 verified and hypothetical open reading frame yeast reference proteins was downloaded from the Saccharomyces Genome Database (SGD) [26] on September 28th, 2012 (http://downloads.yeastgenome.org/sequence/S288C_reference/orf_protein/ webcite). The Homologene database [23] was downloaded from NCBI's FTP site on September 12th, 2011. A protein domain set was obtained from the Conserved Domain Database (CDD version 2.25) [21], which includes domains from CDD and the SMART [27], COG [28], and Pfam [22] databases, with a total of 23,632 protein domains, 10,925 of which map to at least one human protein, and 7,369 map to at least one yeast protein. Functional feature information was collected for CDD domains from the "cddannot.dat" file located in the CDD FTP directory (ftp://ftp.ncbi.nih.gov/pub/mmdb/cdd webcite), totaling 1,727 unique functional features. The non-overlapping set of human, non-synonymous disease mutations was created from the OMIM [29] and Swiss-Prot variant databases obtained from E-utilities and UniProt's FTP directory (http://www.uniprot.org/docs/humsavar webcite) respectively. The resulting human mutation dataset consists of 32,653 mutations related to human diseases. The set of phenotypic yeast mutations was downloaded from SGD (the phenotype_data.tab database obtained from http://www.yeastgenome.org/download-data/curation/ webcite) and was filtered to exclude records without allelic information and records listing the phenotypic change as "normal," as these records refer to mutations with no phenotypic change. Finally, the yeast mutation database was manually curated to extract single point mutations and to ensure that each mutation record referred to a single point mutation. Mutation records referring to multiple mutations for a single phenotype were separated into multiple records. The final yeast mutational database is comprised of 1,490 unique mutations associated with phenotypic changes and is available upon request.

Mapping mutations to protein domains

Hidden Markov models for protein domains from SMART, COG, CDD, and Pfam were built using multiple sequence alignments from CDD with the hmmerbuild tool (HMMer version 2.3.2) [30]. HMMer's hmmpfam tool was then used with the global option to search for complete domains in human proteins from the RefSeq and Swiss-Prot databases. Protein mutations were distributed to protein domain positions by using HMMer's alignment output for all domains aligning to the protein with an E-Value ≤ 0.001 and by assigning mutations on gap regions of the domain model to the last position before the gap. Each mutation was mapped only to the representative protein for each unique gene in the dataset. The methods for mapping domains to human proteins and disease mutations to their domain positions were previously described for our DMDM tool [17]. After distributing each mutation to all domains that map to the protein position in which the mutation was located, 4,283 human protein domains contained at least one disease mutation and 1,687 yeast protein domains contained at least one phenotypic yeast mutation.

Determining the level of conservation at each domain position

For each column j in a protein domain multiple sequence alignment, we used the AL2CO [31] method to determine the entropies using the following formula

<a onClick="popup('http://www.biomedcentral.com/1471-2164/14/S3/S5/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2164/14/S3/S5/mathml/M1">View MathML</a>

(1)

where p(ai,j) is the frequency of amino acid ai at position j. We then estimated a threshold for identifying highly conserved positions by adding one standard deviation to the average of all AL2CO scores on all domain positions. As a result, a threshold of entropy less than or equal to 0.533 was used to determine conservation. The average entropy of each domain model was determined by estimating the mean of the entropy scores for all positions in the domain model.

Estimating domain significance scores (DS-Scores)

We used a method previously developed by our group [18] to estimate the position-based and feature-based DS-Score for each position in the domain. Let n be the total number of mutations in the domain and let L be the number of possible positions in the domain. The random sample, X1, ... , XL consists of the numbers of mutations aggregated to the domain level from a single organism (Figure 1A) or from any number of organisms (Figure 1B) and let X(1) be the smallest of these Xi , and X(2) be the next order of magnitude, ... , and X(L) is the largest Xi. We used the probability mass function, P(X = x), to test whether the n mutations are randomly distributed into L positions. Finally, we defined our position-based DS-Score by the negative logarithm transformation of the binomial probability of observing a cluster of mutations of a particular size, given the total available positions in a domain and the total number of mutations observed. For consistency, when the number of mutations is equal in m positions, we assigned the same score to all m positions. Where k is the kth order of the number of mutations at domain positions, we define the DS-Score as

<a onClick="popup('http://www.biomedcentral.com/1471-2164/14/S3/S5/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2164/14/S3/S5/mathml/M2">View MathML</a>

thumbnailFigure 1. Aggregating mutations using protein domains to form domain hotspots. (A) Mutations from multiple proteins in a single organism sharing a common protein domain mapped to the domain level. These mutations have different phenotypes but can be shown to form a large cluster at a single position in the domain, indicating similarity. (B) Mutations from multiple organisms mapped to a protein domain. These mutations form a large cluster of mutations at the domain level despite having originated from different species. (C) A domain hotspot is formed using mutations from only organism 1, but mutations aligning to the same domain position from organism 2 do not have enough mutations to form a domain hotspot. However, some association can still be drawn due to the existence of a domain hotspot in organism 1 aligning to a mutation in organism 2.

We defined domain hotspots using three levels of significance from the Fisher's scale of evidence for interpreting p-values as described by Yue et al. [32], namely p-values less than or equal to 0.025, 0.05, and 0.10. Using these three levels of significance, we derived three DS-Score thresholds of greater than or equal to 1.6, 1.3, or 1.0 corresponding to p-values less than or equal to 0.025, 0.05, and 0.10 respectively. In addition to position-based DS-Scores, which are based solely on the mutations found at a specific domain position, feature-based DS-Scores were created by distributing the largest position-based DS-Score for each functional feature to all other positions with the same functional feature annotation (see example in Figure 2). We define these functional feature sites as domain positions that have been manually annotated by CDD as having a known functional role (e.g. a DNA binding site or a flexible hinge region) [33]. Perl and R (http://www.r-project.org webcite) were used to determine and distribute the DS-Scores for each domain position.

thumbnailFigure 2. Visual representation of position-based and feature-based domain hotspots. The structure of the Ras-like protein domain (cd00882) of the human protein Cdc42 (PDB: 1CEE_A) is shown in purple (left). The sequence logo (generated using the WebLogo software [62]) represents a subset of the Ras-like protein domain from positions 1-152 (right). The functional feature residues corresponding to the GTP/Mg++ binding site at domain positions 5-11, 54, 110, 111, 113, and 141-143 are highlighted in orange on the structure (left) and are represented as orange boxes below each domain position in the sequence logo (right). Each of these functional feature protein domain positions will have a different position-based DS-Score (estimated based on number of mutations at each domain position) but equal feature-based DS-Score (estimated based on the maximum number of mutations found in any of the binding sites).

Mutations and protein domain redundancy

Due to large domain superfamily hierarchies and duplicate domains from different sources, a single protein mutation can be mapped to many redundant domains. As a consequence, multiple domain hotspots can be identified on different domains that have originated from the same cluster of mutations. On one hand, to ensure that the DS-Scores and domain hotspots are estimated and identified for all protein domain models, all mutations in our analysis were distributed to all domains that map to the mutated protein position. The results of our analysis for each individual domain model, even for those domains from superfamily hierarchies and disparate sources for which we expect high redundancy, are available upon request. On the other hand, to prevent overestimating the number of hotspots shared between yeast and human, we designed a procedure for the non-redundification of domain hotspots. As a result, the numbers reported for all domain hotspots and multi-species domain hotspots originated from a unique set of mutation clusters. We excluded hotspots that originated from identical sets of mutations using the following method to select a unique representative domain for cluster of mutations. The representative domain is only used for visualization and internal calculations and does not affect the reported results. To select the representative domain for a mutation cluster, all domains were ordered alphanumerically from lowest to highest accession identifier. Preference is given to domains that are listed first within the list and were defined as root in the domain hierarchy. If none of the domains in the list is the root domain of a hierarchy, the representative domain is the first domain in the list that contain known functional annotated sites. In addition, when comparing hotspots using the multi-species DS-Score, the representative domain model is selected only among those that are shared among the species.

Assessing the co-occurrence of human diseases and yeast phenotypes

The significance of overlapping human diseases and yeast phenotypes was calculated using a right-sided Fisher's exact test. The Fisher's test for each possible pair of human diseases and yeast phenotypic changes was estimated using the following values: the number of times the human disease and yeast phenotypic change (H and Y respectively) overlap, the number of times H overlaps with a yeast phenotype that is not Y, the number of times Y overlaps with a human disease that is not H, and the total number of overlaps between yeast and human. Human diseases and yeast phenotypes were considered to overlap if the associated mutations were found to localize at the same position of an identical domain. To avoid overestimation due to domain model redundancy, no protein mutation was counted more than once as overlapping with a single human disease or yeast phenotype.

Results

Distribution of mutations in protein domains

To study the distribution of phenotypically relevant mutations at the protein domain level, all proteins from the considered species, i.e., human and yeast, were aligned to one or multiple domains and their mutations mapped into these domains. First, we mapped all phenotypically relevant mutations from yeast and disease mutations from human to protein domains and analyzed the distribution of the location of these mutations with respect to their functional annotation and conservation over species as measured by the entropy of the domain site. Due to the redundancy and large hierarchies within the protein domain database, we observed that the 1,490 yeast and 32,653 human protein mutations were propagated to 11,016 and 323,840 domain mutations respectively (see Table 1). In Table 1, we also list the total number of domain sites with mutations and their breakdown into functionally annotated and conserved domain sites for yeast and human domains. We found a total of 8,186 domain sites mapped to 11,016 phenotypically relevant mutations in yeast, 3,992 of them with known function and 5,950 were conserved domain sites corresponding to 36% and 45% of the total number of yeast domain mutations. Similarly, we reported 323,840 domain mutations in 130,731 domain sites, with 58,096 (18%) corresponding to mutations in functionally annotated sites and 152,524 (47%) to mutations in conserved domain sites. Using Fisher's exact test, we estimated the enrichment of human and yeast mutations in functionally annotated and conserved sites. The results from our analysis produced p-values close to 0 for functionally annotated conserved sites in both human and yeast, which indicates a tendency for both the phenotypically relevant yeast mutations and the human disease mutations to be located at functional feature and conserved sites.

Table 1. Distribution of yeast and human mutations at functional and conserved sites

Transferring mutational information across species through protein domains

When comparing the similarity between yeast and human genomes at the gene and domain levels, we found that the Homologene database can identify human orthologs for only 20% of the yeast genes, while 87% of protein domains found in yeast are also found in human. Most importantly, we quantified how many of the genetic alterations in the yeast and human mutation databases could possibly be related to the other species using either orthologous genes or common protein domains, as illustrated in Figure 3. We found that 435 (29%) of the yeast mutations and 10,187 (31%) of the human mutations can only be related to the other species using common protein domains and not orthologous genes, i.e., mutations located within a domain that is common between yeast and human but not located on orthologous genes. We also identified 610 (41%) yeast and 2,713 (8%) human mutations that could be related to the other species by means of either the common domain or orthologous gene. A small number of mutations, 68 (5%) from yeast and 310 (1%) from human were related to the other species using orthologous genes and not common protein domains, for example, a mutation outside of a protein domain region on an orthologous gene. The remaining mutations, 377 (25%) from yeast and 19,443 (60%) from human, could not be related between species using either method. This analysis was also performed using the orthologous gene information from OMA [34] and InParanoid [35] showing only minor variations in the results.

thumbnailFigure 3. Comparing orthologous genes and common protein domains between yeast and human. Yeast and human protein mutations that can be related to mutations in other organisms through either orthologous genes from Homologene or through common protein domains. For both yeast (A) and human (B) mutations, each portion of the graph represents the mutations that could not be related to the other organism (orange), the mutations that could be related only through Homologene (yellow), the mutations that could be related using only common protein domains (blue), and the mutations that could be related using either method (purple).

Yeast and human mutations form hotspots at protein domain positions

A statistical approach, the DS-Score, was used to identify clusters of mutations sharing the same domain position. The method can be used to cluster mutations from one (Figure 1A) or multiple species (Figure 1B) and clusters from one species can be associated with mutations from another species (Figure 1C). Three DS-Score thresholds were used in this study to identify significantly relevant clusters of mutations at the domain level, namely 1.6, 1.3, and 1.0, corresponding to p-values less than or equal to 0.025, 0.05 and 0.10 respectively. The number of domain hotspots found at each DS-Score threshold is shown in Figure 4. Most importantly, in this first analysis of phenotypically relevant yeast mutations at the protein domain level, we found that they cluster at specific domain positions forming 101, 114, and 135 position-based domain hotspots at DS-Score thresholds greater than or equal to 1.6, 1.3, and 1.0 respectively. The human mutations also form domain hotspots at these thresholds, resulting in 719, 884, and 1,085 position-based domain hotspots at DS-Score thresholds greater than or equal to 1.6, 1.3, or 1.0 respectively.

thumbnailFigure 4. Yeast and human domain hotspots formed at each DS-Score threshold. Yeast and human position-based and feature-based hotspots formed DS-Score thresholds greater than or equal to 1.6, 1.3, and 1.0.

Using a variation of the DS-Score in which we emphasize the similarity of domain sites with the same functional annotation as depicted in Figure 2, we studied the feature-based domain hotspots in yeast and human. Our results show that yeast mutations on domains formed 791, 869, and 1,022 feature-based domain hotspots when using DS-Score thresholds greater than or equal to 1.6, 1.3, and 1.0 respectively, while human mutations yielded 3,197, 3,446, and 3,968 feature-based domain hotspots for these thresholds.

Distribution of phenotypic and disease mutations at conserved and functional domain sites

The distribution of DS-Score domain hotspots located on functionally annotated and conserved sites at each DS-Score threshold is presented in Table 2. Yeast position-based domain hotspots were located at highly conserved sites 55% of the time for the 1.6 threshold, 56% for the 1.3 threshold, and 50% for the 1.0 threshold. Additionally, the position-based domain hotspots for human were found at conserved sites 41%, 39%, and 40% of the time for thresholds 1.6, 1.3, and 1.0 respectively. Despite the enrichment of domain hotspots at conserved domain sites, the correlation between entropy and DS-Score is extremely low with Pearson correlation coefficients of -0.04 and -0.16 for human and yeast data, respectively (see Figures 5A and 5B). To assess the relationship between the number of domain hotspots per domain with respect to the sequence conservation of the domain model, we plotted the average position domain entropy versus their number of hotspots (see Figures 5C and 5D for yeast and human respectively). We found a correlation coefficient of -0.04 for both yeast and human, indicating that almost no correlation exist between the domain divergency and the number of domain hotspots it contains. We also analyzed how frequently the domain hotspots occurred at functional features sites. Our results show that 31%, 31%, and 33% of the position-based domain hotspots are located within annotated functional feature sites for phenotypically altering yeast mutations and 17%, 16%, and 16% for human disease mutations when using 1.6, 1.3, and 1.0 respectively as thresholds to define the domain position-based domain hotspots. When analyzing feature-based domain hotspots at these positions, our results indicated that 67%, 67%, and 63% of feature-based domain hotspots in yeast were localized at conserved sites using the 1.6, 1.3, and 1.0 thresholds respectively and 58%, 57%, and 57% for human feature-based domain hotspots respectively. Similarly, we found that 91% of the yeast feature-based domain hotspots, estimated using 1.6, 1.3 or 1.0 DS-Score thresholds, are located in functionally annotated sites. Likewise, for the human feature-based domain hotspots, we found that 81%, 78%, and 77% (as defined by the 1.6, 1.3, and 1.0 thresholds respectively) were localized within a functionally annotated domain site.

Table 2. Distribution of yeast and human domain hotspots at functional features and conserved sites

thumbnailFigure 5. Correlation between position-based DS-Score and Conservation (estimated by Entropy) for mutations in yeast and human. The relationship between position-based DS-Score (using 1.3 as a threshold) and the conservation of the mutated domain position, as estimated by Entropy, is depicted for yeast (A) and human (B) domains. Additionally, the average entropy for each domain is compared with the number of domain hotspots (using 1.3 as a threshold) in yeast (C) or human (D). In this plot, the total number of domain hotspots is compared to the average Entropy of all positions in the domain model.

Phenotypically relevant mutations tend to cluster at domain positions in yeast and human

Results from Table 3 on the analysis of domain hotspots for both species show that 103, 118, and 177 feature-based domain hotspots at the 1.6, 1.3, and 1.0 thresholds respectively and one position-based hotspot are common between yeast and human. The remaining yeast and human feature-based and position-based domain hotspots were unique to the organism in which they were found. One of the advantages of using additional species in the domain analysis of mutations is its potential for the identification of new domain sites of phenotypic relevance that become statistically significant when more annotated mutations are introduced. Thus, in addition to considering each organism's mutation datasets individually, our method was used to identify statistically significant clusters of a combined set of yeast and human mutations (as shown in Figure 1B). Table 4 outlines results for position-based and feature-based multi-species domain hotspots. In total, we found 861, 1,078, and 1,331 position-based multi-species domain hotspots for the 1.6, 1.3, and 1.0 threshold respectively. In addition, we found 4,243, 4,809, and 5,439 feature-based multi-species domain hotspots for the respective thresholds. We identified 143, 187, and 238 position-based multi-species domain hotspots, for the 1.6, 1.3, and 1.0 DS-Score thresholds respectively, that had not previously been identified when using human or yeast mutations independently. Similarly, we identified 1,243, 1,533, and 1,672 feature-based multi-species domain hotspots, for each DS-Score threshold level respectively, that had not previously been discovered when considering the human or yeast datasets independently. To illustrate the increase in the number of domain hotspots available for inference, Figures 6A and 6B depict the number of position-based and feature-based domain hotspots respectively that can be obtained using only human data, only yeast data, and by using a combined set of mutations from both organisms.

Table 3. Shared and unique domain hotspots between the yeast and human datasets.

Table 4. Multi-species domain hotspots

thumbnailFigure 6. Increase of domain hotspots by incorporating yeast mutational data. The number of position-based (A) and feature-based (B) domain hotspots at the 1.6, 1.3, and 1.0 threshold. This figure visually represents the number of domain hotspots formed from human mutations (blue), the number of domain hotspots from yeast mutations (orange), and the number of new multi-species domain hotspots that did not reach significance in either the yeast or human analyses individually (purple).

Linking domain hotspots with mutations across organisms

We analyzed the ability of domain hotspots to identify relevant mutations in other organisms. In particular, we focused on testing our ability to use yeast hotspots to identify disease mutations in humans. As depicted in Figure 1C, domain hotspots in one organism can be mapped to mutations in another organism through protein domains. When comparing domain hotspots in yeast to test if they could be mapped to any of the currently known human disease mutations, results in Table 5 show that 54 (53.5%), 56 (49.1%), and 65 (48.8%) of position-based domain hotspots in yeast at the 1.6, 1.3, and 1.0 thresholds can be mapped to at least one human mutation. The Fisher's exact test p-values for these results were 2e-26, 6e-25, and 5e-28 for each respective threshold level. Similarly, for yeast feature-based domain hotspots, we identified 562 (71.1%, p-value ≈ 0), 592 (68.1%, p-value ≈ 0), and 666 (65.2%, p-value ≈ 0) yeast domain hotspots that map to at least one human mutation. In addition, in Additional file 1, we list mappings of human domain hotspots into yeast mutational sites, also found to be statistically significant.

Table 5. Mapping of domain hotspots from yeast to known disease mutations in human.

Additional file 1. This file contains information related to the mapping of hotspots in human to mutations in yeast. This file is in PDF format, and can be viewed using Adobe Reader or similar applications.

Format: PDF Size: 133KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Discussion

Our findings highlight the advantages of using protein domains to transfer information related to genetic mutations across species. We show that protein domain models provide a powerful framework for aggregating known phenotypically relevant mutation data across large evolutionary distances, i.e., from human and yeast. As a model organism, yeast is highly studied, well annotated, and easy to manipulate genetically. Thus, it is advantageous to transfer known information from genetic disruptions in yeast for analyzing human mutations. To infer relationships between mutations in different organisms, most studies use orthologous genes as reference to analyze mutations [36]. However, our analysis shows that yeast and human data share more common protein domains than they do orthologous genes. As a result, we show that mutations in both the yeast and human databases are better mapped across organisms when using shared protein domains than when using orthologous genes. For instance, we found that of the 40% of the human mutations that can be related to yeast, only 9% are through gene orthology while 39% can be related using a protein domain framework with an overlap of 5% of mutations that can be related by either domain or gene comparisons. This suggests that transferring mutational information by common protein domains not only vastly increases the number of mutations that can be transferred but also loses very few mutations that would have otherwise been transferred using only gene orthologs. The latter corresponds to, for instance, human disease mutations in genes for which there is a yeast ortholog but located outside a protein domain (only 1% of the human disease mutations in our analysis). Additionally, the domain approach allows the aggregation of mutations from multiple genes in each organism and the identification of relations between mutations located in non-orthologous genes by their functional annotation, which would normally be missed when analyzing the problem using a gene-centric approach.

Our study of phenotypically relevant mutations using a protein domain framework confirms that both yeast and human mutations show a significant tendency to fall within conserved and annotated functional protein domain sites. This is in agreement with the conclusions by Miller et al.[36]. In their study, the authors analyzed human disease mutations on seven disease-associated genes, cystic fibrosis transmembrane conductance regulator (CFTR), glucose-6-phosphate dehydrogenase (G6PD), neural cell adhesion molecule L1 (L1CAM), phenylalanine hydroxylase (PAH), paired box 6 (PAX6), the X-linked retinoschisis gene, and a gene associated with tuberous sclerosis (TSC2). From the study of mutations in these seven genes and their conservation across 20 organisms, including human, the authors concluded that these mutations are in highly conserved protein positions. Here, we reach similar conclusions, but we estimated conservation based on the protein domain models and not at the gene level. Additionally, our findings at the domain level are consistent with Mooney et al. [27]. The authors conducted a study on a set of 231 human genes with known disease mutations, showing that human disease mutations are statistically more likely to be localized within conserved or functionally relevant positions. To summarize, our domain-centric analysis confirms findings from gene-centric studies about enrichment of human disease mutations with respect to conserved and functionally annotated sites while identifying the same characteristics for phenotypically relevant mutations in yeast.

To analyze and compare yeast and human mutations we used the DS-Score method [18] and identified domain hotspots of human and yeast phenotypically relevant mutations. The DS-Score method was previously developed by our team to study human disease mutations and modified in this work to include mutations from both species resulting in the identification of multi-species domain hotspots. We also adapted the method for a multi-species analysis by removing redundant domain hotspots. As an extreme example of the effect of domain redundancy, a single cluster of mutations in the yeast IRE1 gene was propagated to over a hundred domains within the catalytic protein kinase domain family (cl09925 from the CDD database [21]), resulting in 120 domain hotspots having originated from the same cluster of domains. Similarly, domains from multiple sources (such as an identical domain from CDD and Pfam databases) could yield redundant domain hotspots counts. These redundant domain hotspots are correctly estimated and are of great relevance for the analysis of mutations in the context of individual domains. However, when comparing two species using redundant domain hotspots, if the cluster of mutations in the kinase family happens to be common to both species, we would reach the conclusion that there were 120 additional hotspots in common between yeast and human. To avoid overestimation of clusters of mutations that are aggregated at domain level, we defined domain hotspots as those having originated from a unique cluster of mutations and applied this method to the comparison of position-based and feature-based domain hotspots in both species. Using the catalytic protein kinase domain family as an example, each of the 120 domains in which the hotspot was found will retain this information, but only one representative hotspot, cd00180, which is at the top of the hierarchy in that kinase family from CDD [21], was considered for the final domain hotspot count for each species.

In this first-of-a-kind study of yeast mutations at the domain level, we demonstrate that phenotypically relevant mutations in yeast cluster at the domain level just as human disease mutations do, forming yeast and human domain hotspots that are the focus of this study. The hotspots in yeast present the same patterns as human domain hotspots in terms of enrichment at protein domain sites that are conserved and also in sites with known functional annotation. Neither the yeast nor human DS-Scores were found to correlate with conservation (as measured by entropy of the domain site), making the DS-Score method a complement to other methods for prioritization of mutations with putative phenotypic relevance such as SIFT [37], that use conservation as principal feature for their predictions.

We compared the newly found position-based and feature-based domain hotspots in yeast against those arising from human mutations. Given the limited number of mutations with known phenotypic association for each species, we expect the overlap between hotspots in both species to be low or non-existent. Indeed, the only position-based domain hotspot shared between yeast and human is at position 246 on PKc (cd00180), a kinase domain. Surprisingly, the number of hotspots in yeast and human sharing a common domain and functional annotation is 108, 114 and 177 when estimated using feature-based DS-score greater or equal to 1.6, 1.3, and 1.0 respectively. From this analysis, we concluded that the lower threshold provides a significant increase in the number of common feature-based domain hotspots found at the cost of including hotspots with lower significance. On the other hand, the similarity of mutation clusters between species increases when using feature-based DS-Score, a method designed to exploit the manually curated functional annotation provided by CDD [21]. The feature-based method provides some flexibility when comparing across species since it includes hotspots located in different positions for each species, as long as they are positioned in domain sites from the same domain and with identical functional annotation. For instance, using a feature-based threshold of 1.3, we were able to identify an additional 1,533 feature-based domain hotspots when considering DS-Scores generated from the combination of yeast and human mutations. To further underline the similarities between mutation patterns in both species, we analyzed the functional annotation of the domain hotspots and identified several that were unique to the species in which they were found but were localized on domain sites with identical annotation to sites that formed domain hotspots in the other species. Outlined in Table 6, a threshold greater than or equal to 1.3 was used to highlight examples of functional features that contained domain hotspots in both species. Domain positions annotated with the "Active site" functional feature were found to be identified as position-based domain hotspots in yeast 28 times and in human 32 times. Other functional features contained domain hotspots in both organisms as well. Among others, the "GTP/Mg2+ binding site" functional feature formed three yeast domain hotspots and seven human domain hotspots, and the "Substrate binding site," contained five yeast and 17 human domain hotspots.

Table 6. Selected functional feature sites containing domain hotspots in yeast and human

An example of a domain with several domain hotspots in human and yeast that highlights the advantages of using feature-based domain hotspots is the Ras-like GTPase domain (cd00882 from the CDD [21] database). While no position-based domain hotspots in the Ras-like GTPase domain from yeast and human are located at the same domain position, two hotspots are located at domain sites with the same functional annotation. The GTP/Mg++ binding site (highlighted in orange in the structure of the domain shown in Figure 2) contains position-based domain hotspots at position ten for the yeast mutations and at position five for the human mutations. The yeast hotspot in position ten is an example of mutations related at the domain level, originating from several genes, ARF2, SEC4, GTR1, and NOG1 that may not have been identified without analyzing mutations at the domain level due to the low sequence similarity of NOG1 (i.e. when using BLAST with E-Value < = 10−3). The yeast mutations in position ten of cd00882 in genes ARF2, SEC4, GTR1, and NOG1 were associated with increased sensitivity to cold [38], decreased rates of cytokinesis [39], decreased nutrient uptake [40], and inviability due to malformation of the large ribosomal subunit [41], respectively. In human, this domain position on the ARL6 gene contains a mutation associated with Bardet-Biedl syndrome type 3 [42] but not a domain hotspot. On the other hand, in a different position of the binding site, position five of the GTP/Mg++ binding region of the Ras-like GTPase domain contains a human hotspot that aggregates several positions in human genes that have been heavily studied due to their prominence in human diseases. This domain position corresponds to position 12 of both the human HRAS and KRAS genes, from which many mutations have been implicated in diseases such as Costello syndrome [43-46] and Congenital myopathy [47,48] and have also been found to be mutated frequently in somatic tumor samples from patients with follicular thyroid carcinoma [49], pancreatic carcinoma [50], and Schimmelpenning-Feuerstein-Mims syndrome [51], as well as bladder [52], lung [53], and gastric cancers [54]. While both HRAS and KRAS belong to the same protein family and are thus often implicated in the same studies, domain position five also aligns to position 38 of a gene from a different family, GNAT1, which is not similar in sequence to HRAS (i.e., HRAS-GNAST1 E-value of 0.53 using BLAST [55]) or KRAS (i.e., KRAS-GNAST1 BLAST E-value of 0.42). The GNAT1 mutation has been associated with congenital stationary night blindness [56]. Additionally, other mutations were found in the GTP/Mg++ binding pocket that were not members of position-based domain hotspots in either organism that we were able to identify using our feature-based domain hotspots. These mutations, sharing common functional annotation with position-based domain hotspots in both species, have been associated with autoimmune lymphoproliferative syndrome [57], somatic pilocytic astrocytoma [58], Noonan syndrome [59], and chylomiccron retention disease [60]. Thus, by extrapolating hotspots in human and yeast to common functional feature positions, we were able to identify a common functional disruption of the GTP/Mg++ binding pocket that causes different phenotypes when mutated in different genes sharing the same domain in the same organism as well as across organisms.

Furthermore, the domain-centric approach across species introduced here could also be extended to compare the particular phenotypes across organisms that can be related through mutations clustered with this approach, analogous to the phenotype similarities described by McGary et al. [12] using gene orthology. In Table 7, we show a preliminary analysis of this type of phenotype comparison by highlighting the human disease and yeast phenotypic annotations that most frequently co-occur using a domain-centric analysis. A complete list of all significant co-occurrences of human disease and yeast phenotypic changes can be found in Additional file 2. Interestingly, the most significant disease-phenotype co-occurrence was Wilson's disease, a human genetic disorder in which copper accumulates in tissues, with the yeast phenotype related to "Gain of function; metal resistance: increased." This yeast phenotype was derived from a study that analyzed a mutated copper-transporter gene, which resulted in a gain of function mutation that enabled the control of intracellular levels of cadmium through an enhanced cadmium efflux system [61]. While the relationship between the two phenotypes for the Wilson's disease example seems to be clear, for many of the significant co-occurrence of phenotypes found with our method, the relationship between the phenotypes is not apparent. An in-depth analysis of the molecular mechanisms and existing literature, as well as experimental validation will be needed to test and uncover novel hypothesis about molecular similarities between the two species.

Table 7. Human diseases and yeast phenotypic changes that co-occur at domain sites

Additional file 2. This file contains the human disease and yeast phenotypic changes the co-occur at domain sites. This file is in Microsoft Excel format, and can be viewed using Microsoft Excel or similar applications.

Format: XLSX Size: 16KB Download fileOpen Data

Conclusion and future work

This first-of-a-kind study demonstrates the aggregation of mutations from species spanning large evolutionary distances such as yeast and human. Using the DS-Score method as the framework for the integration of molecular characteristics, such as domain location and functional annotation of phenotypically relevant mutations, we were able to identify common mutation patterns from two distantly related species. The domain-centric approach introduced in this paper provides an ideal framework for the analysis of mutational data across species since the number of mutations that can be related from one species to another is much higher than what could be related through gene orthology. The feature-based method to compare mutations across species introduced here represents a unique way to integrate functional annotation of domains into the statistical analysis of mutations, shown here to be extremely advantageous for capturing similarities between mutations across distantly related species. This analysis also suggests that the approach is useful in relating phenotypes from yeast and human resulting from a particular pattern of mutations, such as being localized at the same domain position or functional site. We plan to perform a detailed analysis of the molecular basis of these related phenotypes. Hypotheses derived from this analysis have great potential for discovering new relationships between pathways and networks in both species. In addition, we plan to extend this study to other species to identify patterns across more closely related species including mouse, and to increase the number of known phenotypically relevant domain hotspots by including all mutational data available for a wide range of organisms.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

MGK designed the study, TAP performed the bioinformatics analysis, MGK and TAP prepared and analysed the data and drafted the manuscript. DP developed the algorithm for computing the mutation frequency significance thresholds. All authors read and approved the final manuscript.

Acknowledgements

Authors would like to acknowledge Amy Voltz for her valuable comments and her help obtaining the yeast data of protein mutations. The National Institutes of Health (NIH) 1K22CA143148 to MGK (PI) and R01LM009722 to MGK (co-investigator) and the American Cancer Society, ACS-IRG grant to MGK (PI) supported this work.

Declarations

The publication costs for this article were funded by the above grant.

This article has been published as part of BMC Genomics Volume 14 Supplement 3, 2013: SNP-SIG 2012: Identification and annotation of SNPs in the context of structure, function, and disease. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcgenomics/supplements/14/S3

References

  1. Kann MG: Advances in translational bioinformatics: computational approaches for the hunting of disease genes.

    Brief Bioinform 2010, 11(1):96-110. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  2. Mooney SD, Krishnan VG, Evani US: Bioinformatic tools for identifying disease gene and SNP candidates.

    Methods Mol Biol 2010, 628:307-19. PubMed Abstract | Publisher Full Text OpenURL

  3. Collins FS, Barker AD: Mapping the cancer genome. Pinpointing the genes involved in cancer will help chart a new course across the complex landscape of human malignancies.

    Sci Am 2007, 296(3):50-7. PubMed Abstract | Publisher Full Text OpenURL

  4. The effect of intensive treatment of diabetes on the development and progression of long-term complications in insulin-dependent diabetes mellitus. The Diabetes Control and Complications Trial Research Group

    N Engl J Med 1993, 329(14):977-86. PubMed Abstract | Publisher Full Text OpenURL

  5. The Age-Related Eye Disease Study: a clinical trial of zinc and antioxidants--Age-Related Eye Disease Study Report No. 2

    J Nutr 2000, 130(5S Suppl):1516S-9S. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  6. Cirulli ET, Goldstein DB: Uncovering the roles of rare variants in common disease through whole-genome sequencing.

    Nature reviews Genetics 2010, 11(6):415-25. PubMed Abstract | Publisher Full Text OpenURL

  7. Li B, Leal SM: Methods for detecting associations with rare variants for common diseases: application to analysis of sequence data.

    American journal of human genetics 2008, 83(3):311-21. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  8. Wang K, Li M, Hakonarson H: Analysing biological pathways in genome-wide association studies.

    Nat Rev Genet 2010, 11(12):843-54. PubMed Abstract | Publisher Full Text OpenURL

  9. Califano A, Butte AJ, Friend S, Ideker T, Schadt E: Leveraging models of cell regulation and GWAS data in integrative network-based association studies.

    Nat Genet 2012, 44(8):841-7. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  10. Waterston RH, Lindblad-Toh K, Birney E, Rogers J, Abril JF, Agarwal P, Agarwala R, Ainscough R, Alexandersson M, An P, et al.: Initial sequencing and comparative analysis of the mouse genome.

    Nature 2002, 420(6915):520-62. PubMed Abstract | Publisher Full Text OpenURL

  11. Gibbs RA, Weinstock GM, Metzker ML, Muzny DM, Sodergren EJ, Scherer S, Scott G, Steffen D, Worley KC, Burch PE, et al.: Genome sequence of the Brown Norway rat yields insights into mammalian evolution.

    Nature 2004, 428(6982):493-521. PubMed Abstract | Publisher Full Text OpenURL

  12. McGary KL, Park TJ, Woods JO, Cha HJ, Wallingford JB, Marcotte EM: Systematic discovery of nonobvious human disease models through orthologous phenotypes.

    Proc Natl Acad Sci USA 2010, 107(14):6544-9. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. Ekman D, Bjorklund AK, Frey-Skott J, Elofsson A: Multi-domain proteins in the three kingdoms of life: orphan domains and other unassigned regions.

    J Mol Biol 2005, 348(1):231-43. PubMed Abstract | Publisher Full Text OpenURL

  14. Doolittle RF: The multiplicity of domains in proteins.

    Annu Rev Biochem 1995, 64:287-314. PubMed Abstract | Publisher Full Text OpenURL

  15. Bornberg-Bauer E, Beaussart F, Kummerfeld SK, Teichmann SA, Weiner J: The evolution of domain arrangements in proteins and interaction networks.

    Cell Mol Life Sci 2005, 62(4):435-45. PubMed Abstract | Publisher Full Text OpenURL

  16. Zhong Q, Simonis N, Li QR, Charloteaux B, Heuze F, Klitgord N, Tam S, Yu H, Venkatesan K, Mou D, et al.: Edgetic perturbation models of human inherited disorders.

    Mol Syst Biol 2009, 5:321. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  17. Peterson TA, Adadey A, Santana-Cruz I, Sun Y, Winder A, Kann MG: DMDM: domain mapping of disease mutations.

    Bioinformatics 2010, 26(19):2458-9. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. Peterson TA, Nehrt NL, Park D, Kann MG: Incorporating molecular and functional context into the analysis and prioritization of human variants associated with cancer.

    J Am Med Inform Assoc 2012, 19(2):275-83. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  19. Yue P, Forrest WF, Kaminker JS, Lohr S, Zhang Z, Cavet G: Inferring the functional effects of mutation through clusters of mutations in homologous proteins.

    Human mutation 2010, 31(3):264-71. PubMed Abstract | Publisher Full Text OpenURL

  20. Nehrt NL, Peterson T, Park D, Kann MG: Domain landscapes of somatic mutations in cancer.

    BMC Genomics 2012, 13(Suppl 4):S9. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  21. Derbyshire MK, Lanczycki CJ, Bryant SH, Marchler-Bauer A: Annotation of functional sites with the Conserved Domain Database.

    Database (Oxford) 2012, 2012:bar058. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Punta M, Coggill PC, Eberhardt RY, Mistry J, Tate J, Boursnell C, Pang N, Forslund K, Ceric G, Clements J, et al.: The Pfam protein families database.

    Nucleic Acids Res 2012, 40(Database):D290-301. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  23. Geer LY, Marchler-Bauer A, Geer RC, Han L, He J, He S, Liu C, Shi W, Bryant SH: The NCBI BioSystems database.

    Nucleic Acids Res 2010, 38(Database):D492-6. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. Pruitt KD, Tatusova T, Maglott DR: NCBI reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins.

    Nucleic Acids Res 2007, 35(Database):D61-5. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. Boeckmann B, Bairoch A, Apweiler R, Blatter MC, Estreicher A, Gasteiger E, Martin MJ, Michoud K, O'Donovan C, Phan I, et al.: The SWISS-PROT protein knowledgebase and its supplement TrEMBL in 2003.

    Nucleic Acids Res 2003, 31(1):365-70. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  26. Cherry JM, Hong EL, Amundsen C, Balakrishnan R, Binkley G, Chan ET, Christie KR, Costanzo MC, Dwight SS, Engel SR, et al.: Saccharomyces Genome Database: the genomics resource of budding yeast.

    Nucleic Acids Res 2012, 40(Database):D700-5. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. Letunic I, Copley RR, Pils B, Pinkert S, Schultz J, Bork P: SMART 5: domains in the context of genomes and networks.

    Nucleic Acids Res 2006, 34(Database):D257-60. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  28. Tatusov RL, Fedorova ND, Jackson JD, Jacobs AR, Kiryutin B, Koonin EV, Krylov DM, Mazumder R, Mekhedov SL, Nikolskaya AN: The COG database: an updated version includes eukaryotes.

    BMC Bioinformatics 2003, 4:41. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  29. McKusick VA: Mendelian Inheritance in Man and its online version, OMIM.

    Am J Hum Genet 2007, 80(4):588-604. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  30. Eddy SR: Hidden Markov models.

    Curr Opin Struct Biol 1996, 6(3):361-365. PubMed Abstract | Publisher Full Text OpenURL

  31. Pei J, Grishin NV: AL2CO: calculation of positional conservation in a protein sequence alignment.

    Bioinformatics 2001, 17(8):700-12. PubMed Abstract | Publisher Full Text OpenURL

  32. Yue P, Forrest WF, Kaminker JS, Lohr S, Zhang Z, Cavet G: Inferring the functional effects of mutation through clusters of mutations in homologous proteins.

    Hum Mutat 2010, 31(3):264-71. PubMed Abstract | Publisher Full Text OpenURL

  33. Marchler-Bauer A, Zheng C, Chitsaz F, Derbyshire MK, Geer LY, Geer RC, Gonzales NR, Gwadz M, Hurwitz DI, Lanczycki CJ, et al.: CDD: conserved domains and protein three-dimensional structure.

    Nucleic Acids Res 2013, 41(D1):D348-52. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  34. Altenhoff AM, Schneider A, Gonnet GH, Dessimoz C: OMA 2011: orthology inference among 1000 complete genomes.

    Nucleic Acids Res 2011, 39(Database):D289-94. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  35. Ostlund G, Schmitt T, Forslund K, Kostler T, Messina DN, Roopra S, Frings O, Sonnhammer EL: InParanoid 7: new algorithms and tools for eukaryotic orthology analysis.

    Nucleic Acids Res 2010, 38(Database):D196-203. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  36. Miller MP, Kumar S: Understanding human disease mutations through the use of interspecific genetic variation.

    Hum Mol Genet 2001, 10(21):2319-28. PubMed Abstract | Publisher Full Text OpenURL

  37. Ng PC, Henikoff S: SIFT: Predicting amino acid changes that affect protein function.

    Nucleic Acids Res 2003, 31(13):3812-4. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  38. Peyroche A, Paris S, Jackson CL: Nucleotide exchange on ARF mediated by yeast Gea1 protein.

    Nature 1996, 384(6608):479-81. PubMed Abstract | Publisher Full Text OpenURL

  39. VerPlank L, Li R: Cell cycle-regulated trafficking of Chs2 controls actomyosin ring stability during cytokinesis.

    Mol Biol Cell 2005, 16(5):2529-43. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  40. Gao M, Kaiser CA: A conserved GTPase-containing complex is required for intracellular sorting of the general amino-acid permease in yeast.

    Nat Cell Biol 2006, 8(7):657-67. PubMed Abstract | Publisher Full Text OpenURL

  41. Fuentes JL, Datta K, Sullivan SM, Walker A, Maddock JR: In vivo functional characterization of the Saccharomyces cerevisiae 60S biogenesis GTPase Nog1.

    Mol Genet Genomics 2007, 278(1):105-23. PubMed Abstract | Publisher Full Text OpenURL

  42. Fan Y, Esmail MA, Ansley SJ, Blacque OE, Boroevich K, Ross AJ, Moore SJ, Badano JL, May-Simera H, Compton DS, et al.: Mutations in a member of the Ras superfamily of small GTP-binding proteins causes Bardet-Biedl syndrome.

    Nat Genet 2004, 36(9):989-93. PubMed Abstract | Publisher Full Text OpenURL

  43. Yoshida R, Fukushima Y, Ohashi H, Asoh M, Fukuyama Y: The Costello syndrome: are nasal papillomata essential?

    Jpn J Hum Genet 1993, 38(4):437-44. PubMed Abstract | Publisher Full Text OpenURL

  44. Smith LP, Podraza J, Proud VK: Polyhydramnios, fetal overgrowth, and macrocephaly: prenatal ultrasound findings of Costello syndrome.

    Am J Med Genet A 2009, 149A(4):779-84. PubMed Abstract | Publisher Full Text OpenURL

  45. Gripp KW, Stabley DL, Nicholson L, Hoffman JD, Sol-Church K: Somatic mosaicism for an HRAS mutation causes Costello syndrome.

    Am J Med Genet A 2006, 140(20):2163-9. PubMed Abstract | Publisher Full Text OpenURL

  46. Sol-Church K, Stabley DL, Demmer LA, Agbulos A, Lin AE, Smoot L, Nicholson L, Gripp KW: Male-to-male transmission of Costello syndrome: G12S HRAS germline mutation inherited from a father with somatic mosaicism.

    Am J Med Genet A 2009, 149A(3):315-21. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  47. van der Burgt I, Kupsky W, Stassou S, Nadroo A, Barroso C, Diem A, Kratz CP, Dvorsky R, Ahmadian MR, Zenker M: Myopathy caused by HRAS germline mutations: implications for disturbed myogenic differentiation in the presence of constitutive HRas activation.

    J Med Genet 2007, 44(7):459-62. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  48. Lo IF, Brewer C, Shannon N, Shorto J, Tang B, Black G, Soo MT, Ng DK, Lam ST, Kerr B: Severe neonatal manifestations of Costello syndrome.

    J Med Genet 2008, 45(3):167-71. PubMed Abstract | Publisher Full Text OpenURL

  49. Dajee M, Lazarov M, Zhang JY, Cai T, Green CL, Russell AJ, Marinkovich MP, Tao S, Lin Q, Kubo Y, et al.: NF-kappaB blockade and oncogenic Ras trigger invasive human epidermal neoplasia.

    Nature 2003, 421(6923):639-43. PubMed Abstract | Publisher Full Text OpenURL

  50. Motojima K, Urano T, Nagata Y, Shiku H, Tsurifune T, Kanematsu T: Detection of point mutations in the Kirsten-ras oncogene provides evidence for the multicentricity of pancreatic carcinoma.

    Ann Surg 1993, 217(2):138-43. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  51. Rijntjes-Jacobs EG, Lopriore E, Steggerda SJ, Kant SG, Walther FJ: Discordance for Schimmelpenning-Feuerstein-Mims syndrome in monochorionic twins supports the concept of a postzygotic mutation.

    Am J Med Genet A 2010, 152A(11):2816-9. PubMed Abstract | Publisher Full Text OpenURL

  52. Di Micco R, Fumagalli M, Cicalese A, Piccinin S, Gasparini P, Luise C, Schurra C, Garre M, Nuciforo PG, Bensimon A, et al.: Oncogene-induced senescence is a DNA damage response triggered by DNA hyper-replication.

    Nature 2006, 444(7119):638-42. PubMed Abstract | Publisher Full Text OpenURL

  53. Santos E, Martin-Zanca D, Reddy EP, Pierotti MA, Della Porta G, Barbacid M: Malignant activation of a K-ras oncogene in lung carcinoma but not in normal tissue of the same patient.

    Science 1984, 223(4637):661-4. PubMed Abstract | Publisher Full Text OpenURL

  54. Haigis KM, Kendall KR, Wang Y, Cheung A, Haigis MC, Glickman JN, Niwa-Kawakita M, Sweet-Cordero A, Sebolt-Leopold J, Shannon KM, et al.: Differential effects of oncogenic K-Ras and N-Ras on proliferation, differentiation and tumor progression in the colon.

    Nat Genet 2008, 40(5):600-8. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  55. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool.

    J Mol Biol 1990, 215(3):403-410. PubMed Abstract OpenURL

  56. Dryja TP, Hahn LB, Reboul T, Arnaud B: Missense mutation in the gene encoding the alpha subunit of rod transducin in the Nougaret form of congenital stationary night blindness.

    Nat Genet 1996, 13(3):358-60. PubMed Abstract | Publisher Full Text OpenURL

  57. Oliveira JB, Bidere N, Niemela JE, Zheng L, Sakai K, Nix CP, Danner RL, Barb J, Munson PJ, Puck JM, et al.: NRAS mutation causes a human autoimmune lymphoproliferative syndrome.

    Proc Natl Acad Sci USA 2007, 104(21):8953-8. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  58. Sharma MK, Zehnbauer BA, Watson MA, Gutmann DH: RAS pathway activation and an oncogenic RAS mutation in sporadic pilocytic astrocytoma.

    Neurology 2005, 65(8):1335-6. PubMed Abstract | Publisher Full Text OpenURL

  59. Schubbert S, Zenker M, Rowe SL, Boll S, Klein C, Bollag G, van der Burgt I, Musante L, Kalscheuer V, Wehner LE, et al.: Germline KRAS mutations cause Noonan syndrome.

    Nat Genet 2006, 38(3):331-6. PubMed Abstract | Publisher Full Text OpenURL

  60. Jones B, Jones EL, Bonney SA, Patel HN, Mensenkamp AR, Eichenbaum-Voline S, Rudling M, Myrdal U, Annesi G, Naik S, et al.: Mutations in a Sar1 GTPase of COPII vesicles are associated with lipid absorption disorders.

    Nat Genet 2003, 34(1):29-31. PubMed Abstract | Publisher Full Text OpenURL

  61. Shiraishi E, Inouhe M, Joho M, Tohoyama H: The cadmium-resistant gene, CAD2, which is a mutated putative copper-transporter gene (PCA1), controls the intracellular cadmium-level in the yeast S. cerevisiae.

    Curr Genet 2000, 37(2):79-86. PubMed Abstract | Publisher Full Text OpenURL

  62. Crooks GE, Hon G, Chandonia JM, Brenner SE: WebLogo: a sequence logo generator.

    Genome Res 2004, 14(6):1188-90. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL