Ranking insertion, deletion and nonsense mutations based on their effect on genetic information
© Zia and Moses; licensee BioMed Central Ltd. 2011
Received: 24 March 2011
Accepted: 22 July 2011
Published: 22 July 2011
Genetic variations contribute to normal phenotypic differences as well as diseases, and new sequencing technologies are greatly increasing the capacity to identify these variations. Given the large number of variations now being discovered, computational methods to prioritize the functional importance of genetic variations are of growing interest. Thus far, the focus of computational tools has been mainly on the prediction of the effects of amino acid changing single nucleotide polymorphisms (SNPs) and little attention has been paid to indels or nonsense SNPs that result in premature stop codons.
We propose computational methods to rank insertion-deletion mutations in the coding as well as non-coding regions and nonsense mutations. We rank these variations by measuring the extent of their effect on biological function, based on the assumption that evolutionary conservation reflects function. Using sequence data from budding yeast and human, we show that variations which that we predict to have larger effects segregate at significantly lower allele frequencies, and occur less frequently than expected by chance, indicating stronger purifying selection. Furthermore, we find that insertions, deletions and premature stop codons associated with disease in the human have significantly larger predicted effects than those not associated with disease. Interestingly, the large-effect mutations associated with disease show a similar distribution of predicted effects to that expected for completely random mutations.
This demonstrates that the evolutionary conservation context of the sequences that harbour insertions, deletions and nonsense mutations can be used to predict and rank the effects of the mutations.
Genetic variations contribute to normal phenotypic variation . For human, it is estimated that there are more than 10 million SNPs (i.e. 1 in 300 base pairs on average) with an observed minor allele frequency of ≥ 1% in the population . Recent advances in sequencing technologies  have enabled rapid discovery of other types of variations, including mutations expected to have very large effects on protein function such as frame shifting insertions and deletions (indels) and nonsense mutations (mutations that introduce premature stop codons). Amazingly, insertions and deletions are also abundant in the human genome with sizes ranging from single to several million base pairs (bp) [4, 5]. For example, in 179 human genomes there were 1.13 million short indels identified  indicating an estimate of 1 million indels per human genome (1 in 3600 bps on average). Similarly, recent studies have revealed a surprising number (on average, 80-100) of nonsense mutations per genome observed . Consistent with these data from humans, recent studies in model organisms such as yeast have also revealed an abundance of variants that would be expected to have large effects on phenotype [7, 8].
Large-effect mutations (defined here as frame-shift causing indels, nonsense mutations and indels that disrupt highly conserved non-coding DNA) are often assumed to have significant functional impacts, and therefore are likely to cause disease. However, discovery of a large numbers of these variants suggest that many of them might in fact have no (or little) impact on function . A quick survey of dbSNP  reveals that of 8311 nonsense mutations (NM) and 33285 frame-shifting (FS) indels have been identified in humans, of which only 2234 (27%) NMs and 801 (2.4%) FS indels have been clinically studied (i.e. have records in OMIM  and LSDB  and are associated with diseases). Thus, identification of the fraction of these variations that contribute to diseases from those with little effect on function is a great practical challenge [13, 14].
Even in cases where the variations are detected in individuals with diseases, identification of the causative variations is a major challenge. For instance, there are 2729 structural variants (including indels) catalogued in the COSMIC project  that are potentially associated with lung cancer. Similar observations can be made at the level of individual important disease genes: studies of TP53 (mutated in various kinds of cancer) have identified 1256 somatic, 127 cell-line and 36 germline FS indels  as well as 95 somatic, 49 cell-line and 15 germline NMs. Similarly, there are 301 FS indels as well as 160 NMs identified in CFTR (a single mutant gene that is identified to be associated with Cystic Fibrosis ). In all these cases, the identification of causative variations in the list of potential candidates remains a challenge.
As the number of observed mutations increases, it quickly becomes infeasible for researchers to manually assess the impact of each one in laboratory. It is therefore becoming absolutely essential to rank the effect of these variations in terms of their impact to define priority in clinical research as well as to weight their effects in association analyses.
There is a wide range of computational methods that predict the effect of SNP on protein function (for a survey of these methods see  and references therein). Despite existing interest, to the extent of our knowledge, none of these methods are able to deal with indels and NMs. Here we propose to use evolutionary conservation principles to rank the effect of these variations on genetic information.
Evolutionary conservation has been previously used for predicting the effect of SNPs on protein function [19, 7] as well as on non-coding DNA . In SIFT  the conservation of amino-acid residues are measured using protein sequence homology. In this method, a non-synonymous SNP (nsSNP) that substitutes a highly conserved residue (i.e. a residue that is less observed in an alignment of homologue sequences) is predicted to have a more deleterious effect. Similarly, in the case of non-coding DNA, in the so-called LR test , a SNP is predicted to have more deleterious effect if it causes greater change in the rate of evolution of the DNA site that it alters.
We argue that the underlying ideas used in these methods, i.e. using the evolutionary conservation context [20, 21] of the sequence that harbours the mutation for assessing the mutation impact, can be extended to rank the net effect of indels and NMs. In particular, for the protein coding sequences, we propose a scoring scheme that measures the amount of the loss of protein "information content" caused by a NM or FS indel. We expect the variations that interfere with conserved residues of a protein to a larger extent to have more deleterious effects. For the non-coding DNA, we propose to use a likelihood-ratio scoring scheme to measure the conservation of the DNA bases that harbour the indels. We argue that indels that fall in highly conserved DNA sites are expected to have more deleterious effects.
To evaluate our hypothesis, we study the effect of indels and NMs in a population of Saccharomyces cerevisiae (S.cer) yeast . We provide evidence that mutations that disrupt the most highly conserved regions segregate at significantly lower allele frequencies. The paucity of variations in lower allele frequencies suggests that highly deleterious mutations have been removed from the population .
Next, we assess the effects of FS indels as well as NMs on human proteins [10, 6, 11, 12]. We show that variations with no disease association tend to cause less information loss than those associated with disease, suggesting that variations not associated with disease are likely to have less deleterious effects. We further show that NMs that cause higher information loss in human proteins segregate significantly in lower allele frequencies suggesting that not all NMs have the same deleterious effects. We argue that the scoring scheme that quantifies the information loss can be used to rank the effect of mutations in human population.
Large-Effect mutations in protein coding DNA
SIFT  predicts a non-synonymous SNP to be deleterious if it disrupts a highly conserved amino acid residue in a protein sequence, where the conservation of any residue is measured using an alignment of homologous protein sequences.
We propose that the same principle applies to variations that affect more than one amino acid residue. We know that a FS indel changes the translated amino acid residues from its position to the end of the protein, and a NM causes a premature termination of the amino acid sequence. Thus, in order to extend the approach taken in SIFT, we simply add up the effects of each residue affected. In doing so we are assuming that changes in amino acid residues contribute independently to the overall function of the protein, obviously an over-simplification (see Discussion).
Scores used in SIFT are based on a "normalized" transition probability matrix (TPM) that is built using an alignment of protein sequences homologous to the target protein. The TPM is not suitable for adding the effect of mutations on multiple residues because each column is normalized such that for any residue, the most likely transition is normalized to 1. Since the absolute maxima of the observed probabilities are lost in the normalization, the scores cannot be used to compare the residues or, in our case, to add the effect of substitutions for a number of residues together.
Therefore, we modify the scoring method used in SIFT by using the conventional definition of information carried by biological sequences (e.g. ) (for a review see ). According to this definition, protein residues that are important for the species are evolutionarily conserved and therefore have a statistically different distribution compared to the freely evolving residues that are under no selection. Biological information corresponds to the difference between the distribution of conserved and non-conserved (freely evolving) sequence. We quantify this difference as follows.
where, for instance, XiA = 1 if the residue at ith column is amino-acid residue "A". The score S(Xi) shows how likely is that the residue Xi is generated from the distribution fi (compared to the background distribution gi). This score is closely related to the relative information and hence our interpretation of information carried by Xi. The total score of a sequence, S(X), is defined as the sum of information carried by individual residues.
By combining the effects of change in individual residues we assume that residues contribute independently to the information loss (or gain) of the sequence. We normalize the loss of information to the total score, S (of the wild-type (WT) sequence), to obtain a dimensionless quantity that can be used to compare the information loss between protein sequences.
Our predictions are based on the information loss score D which defines a normalized ratio of the information content S of the protein without (i.e. the WT) and with the mutation (mutant). The larger the score D is the greater is the information loss due to the mutation. We use this score to rank the effect of mutations.
examples of mutations in human genes
Mutation at base
Effect of mutation in non-coding DNA
It is widely accepted that the non-coding DNA harbours many functional elements [27, 28] and that variations in these regions can have phenotypic effects and cause disease [29–34]. We use conservation of bases in the non-coding DNA to assess the impact of these variations. While in principle, a similar approach used for computing the information loss in coding regions could also be applied to the non-coding DNA, building sequence alignments of homologous non-coding DNA from distantly related species is infeasible due to the relatively fast rates of evolution.
whereλ syn is the average rate of evolution of synonymous mutations in the coding DNA of the protein and is the maximum likelihood estimate of the rate of evolution of the non-coding DNA site. Here T is the evolutionary tree of the species being used in building the DNA alignments. The LLR measures how fast a DNA site is evolving compared to the synonymous rate of evolution of the protein coding DNA. Slowly evolving sites (i.e. more conserved) give larger LLR values. For these sites, the likelihood (in the LLR numerator) that the site is evolving according to the rate λ* (which is no more than λ syn ) is significantly greater than the likelihood that the site is evolving as fast as the synonymous rate of protein coding region. Alternatively, the DNA sites that evolve at a rate similar to the synonymous sites in the coding region are not conserved and therefore give a lower LLR scores.
The results of our analysis provide evidence that LLR score can also be used to measure the effect of indels that disrupt the conserved DNA sites. A similar approach was previously used in  to assess the effect of SNPs in non-coding as well as the protein-coding regions .
Genome-wide population analysis with mutations in yeast protein coding regions
We sought to test whether our scores for large-effect mutations reflected their functional impact. More deleterious variants are expected to segregate at lower frequencies in the population and occur at lower densities that would be expected of neutral variants . Therefore, in a natural population, we expect mutations with larger predicted effects to segregate at lower frequencies and be found at lower densities than mutations with smaller predicted effects.
In order to test whether our information-based score differentiates FS indels that are likely to be deleterious from those that are likely to be neutral, we computed the fraction of high frequency alleles defined as DAF > 0.2. We compared this fraction between genes with higher score, i.e. D > 0.3, which we considered to have a greater deleterious effects, and genes with score less than 0.3 (which we considered moderate deleterious effects).
We found, for the case of FS indels (see Figure 3c) that the fraction of high frequency alleles for the genes with D > 0.3 was significantly less than the fraction of high frequency alleles for genes with for the genes with D < 0.3 (1/35 versus 7/22 with Fisher's exact test p = 0.004) (due to sparse spectrums, we did not use the KS test for this case). For the NMs, the results were inconclusive because we did not observe mutations with high allele frequency (Figure 3d). The paucity of high frequency alleles is consistent with stronger purifying selection on the mutations we predicted to be more deleterious (using the score D).
Genome-wide population analysis of indels in yeast non-coding DNA
We identified indels that fall within the promoter regions (1000 bp upstream) of genes in a population of 39 strains of S.cer budding yeast . For each identified indel, we computed the conservation of surrounding DNA site (of width 20 bps) using an alignment of non-coding DNA sequences of orthologues genes from 5 closely related yeast species (see Methods and Figure 2).
We sought to test whether indels that fall in highly conserved DNA sites segregate at lower DAFs compared to those that do not disrupt conserved DNA sites (i.e. that fall in relatively non-conserved regions). For this purpose, we define a LLR score threshold THR = 16.6 for dividing indels into two sets: those that fall into highly conserved DNA sites (i.e. sites with LLR > 16.6) and those that fall in less conserved DNA site (i.e. sites with LLR < 16.6). The threshold THR = 16.6 is estimated such that the difference between the selection coefficients of the two classes is maximal (see Methods for detail of using Poisson random fields for computing the optimal threshold).
Figure 3e shows the DAFs for the two sets of indels. The allele frequency spectrum of indels falling into highly conserved DNA sites is shifted towards lower DAFs significantly (KS-test p = 0.0033) suggesting that these mutations are under stronger purifying selection.
To test for the effects of selection on the density of indels, we compared the distribution of DNA sites with respect to their LLR scores with what is expected if the indels were placed randomly. For this purpose, we generated 50 sets of indels distributed randomly in the 1000 bp upstream of all genes in the reference S.cer and computed the LLR scores for the corresponding DNA sites in which they fall. Figure 4e shows that the distribution for DNA sites in our dataset is significantly different than the random dataset (KS test, p < 10-6), such that there is a great enrichment of sites with lower LLR scores. This suggests that indels at highly conserved DNA sites have been removed from the population by purifying selection.
Ranking of mutations in yeast
Our method identifies candidates for new deleterious variations in a pool of genes with mutations. We ranked the mutations in our yeast dataset in terms of their deleterious effects on the corresponding proteins. In the following, we study the top ranked FS indels and the NMs with lowest scores.
highest and lowest D-scores
a) FS indels
Mutation at base
One of the other top ranked indel predictions was in TFB3. A reduction in the function of the gene TFB3 is associated with an increase in resistance to the same drug . Interestingly, we observed the expected effects in the 2 individuals that carry the FS indel (and therefore are predicted to lack the function for this gene). Specifically, the 2 individuals were among the 8 individuals that had REs > -2 (average RE in strains that had genotype data at this locus is -3.06). Because of the small number of individuals that carry the putatively deleterious FS indel alleles (here 2) we were not able to test the significance of these phenotypic observations. However, these examples show the practical uses of the proposed methods.
We further studied the bottom 5 genes with NMs with lowest scores in our dataset. Table 2 shows that while it is possible to observe NMs in essential genes (ranks 74 and 70) , our method predicts that these mutations have no substantial effect on the function. These mutations are located in the C-termini of these genes.
Variations in the human population
To test whether our methods can be applied to the variations in the human population we examined genes with FS indels and NMs reported in dbSNP . We categorized the variations into two classes: variations that have records of diseases association in OMIM  and LSDB  and variations with no such records. We expect that the latter class (i.e., variations with unknown disease association) to contain mutations with less harmful effects.
We then sought to study the allele frequency spectrum of these variations. Heterozygousity information was only available for the NMs with no disease association. We used the heterozygousity to compute the minor-allele frequency (MAF) of the NMs that were used in the following MAF spectrum analysis. We studied the segregation of NMs in the human population by comparing the spectrums of MAF of the NMs that cause greater information loss (i.e. 1101 NMs with D > 0.3) and the NMs with lower effects (i.e. 505 NMs with D < 0.3). As for the non-coding indels above, we computed D = 0.3 as the optimal threshold using the Poisson random fields method (see Methods). We found that the spectrum of MAF of NMs with higher information loss (Figuure 3f) is significantly skewed to left (KS test, p = 0.002) indicating that these variations segregate significantly at lower allele frequencies. Thus, these mutations appear to have more deleterious effects in the population.
To study the effects of selection on mutation' density, we compared the observed distribution of scores with that expected if variations were randomly placed. To do so, we computed the score D for a large number of FS indels as well as NMs placed randomly on the human genes in our dataset. Figures 4c and 4d compare the distribution of these scores for FS indels and NMs, respectively, with variations in our dataset divided into two classes: with and without disease associations. The insignificant different between distributions for variations with disease association compared to randomly generated variations (KS test, p = 0.09 for FS indels and p = 0.09 for NMs) suggests that these variations are distributed randomly throughout human genes. On the other hand, there is a significant difference between the distributions of variations that have no disease association compared to randomly placed variations (KS test, p < 10-6 for FS indels and p < 10-4 for NMs). Similarly, there is also a significant difference between the distributions of variations with and without disease association (KS test p < 10-6 for FS indels and KS test p = 0.0006 for NMs). The significant abundance of mutations with lower deleterious effects in the data with no disease association, or in other words, the paucity of variations with higher information loss scores, indicates that purifying selection had acted on highly deleterious variations.
An example for the application of the method to study human genes
Our proposed methods are useful for practical purposes to sort a huge number of FS indels, NMs, as well as indels in non-coding DNA in terms of their deleterious effect. It is important to note that our methods do not seek to classify variations into deleterious and non-deleterious but rather to rank their effect for further analysis and laboratory experiments.
For the variations in the protein-coding DNA, the proposed score is built upon the principle assumption that the effect of nonsense or FS indel mutations on protein can be computed as sum of effects due to individual residues. This is obviously an over-simplification that is widely accepted in statistically modeling the individual columns of an alignment independently. A more complex method that considered correlations between each residue was also implemented using a profile-HMM based on a generative hidden-Markov model  (data not shown here). The score S were computed as the likelihood of the sequence given the profile-HMM . Similar prediction results were observed, i.e. sequences with mutations in more conserved regions resulted in lower scores.
We observed a strong correlation between the position of the nonsense or FS mutations and the loss of information they cause (Additional file 1, Figure S1). We were not able to demonstrate that the D score outperforms the percentage of the protein that is truncated (the "length lost"). When we compare the distribution of the length lost to the random expectation (Additional file 1, Figure S2) we find that the length-lost score appears to show less deviation from the random expectation than the D score for the human data (Additional file 1, Figure S2c,d). This is consistent with the hypothesis that the D score captures more information than the length lost. While simply considering the number of residues affected provides a reasonable guess at the impact of mutations "on average", there are cases for which the position of the mutation does not reflect its effect on evolutionary conservation (off-diagonal points in Additional file 1, Figure S1). Furthermore, we believe that the D score represents a more principled approach to quantifying the importance of these variants because it directly measures evolutionary information, and because it is consistent with previous approaches to quantify the effects of variants, such as SIFT. However, if multiple sequence alignments are not available, the length lost might also provide a reasonable substitute to quantify the effect of a FS indel or NM.
Identification of causative mutations for diseases remains a challenge even for the case of single genes, let alone in cases where mutations are studied in a network of genes and regulatory elements (e.g. variations affecting genetic pathways). Due to the overwhelming abundance of variations, the information loss score, that captures the evolutionary conservation context of the sequences harbouring mutations, seems to be a good candidate for weighting variations in large-scale association analyses [53, 54].
Information loss in protein coding DNA
To obtain the information loss, D, we compute the scores of the WT and the mutant sequences against a position weight matrix as it is defined above. The mutant proteins with highest scores D are more likely to carry a highly deleterious mutation.
The PWM is built using an alignment of protein sequences homologous to the WT protein. We begin with the WT protein sequence as the query to PSI-BLAST  with adequate number of iterations (2-5) to collect a set of sequences that are locally homologous to the WT protein (similar to SIFT ). The database of protein sequences we used for this study was SWISS-PROT 3.9 . We then remove sequences that are greater than 90% identical to the target sequence from the alignment. This is to maintain a minimum degree of diversity between sequences and to avoid biasing the estimation of the PWM towards closely related species. There are cases where PSI-BLAST did not return sufficient hits to build the alignments. We exclude mutations on genes from our dataset in these cases.
Any column of the PWM, f, is the maximum likelihood estimate of the distribution of amino acid residues observed in the alignment. In ideal cases where there is sufficiently large number of sequences in the alignment, the PWM is simply a matrix with columns equal to the relative frequency of each observed amino acid residue at that column. However, in practice, due to a relatively large number of residues (i.e. 20) compared to the alignment size, there are chances that many residues are not observed and hence their corresponding entry in the PWM is zero. We resolve this issue by considering a minimum number of pseudo-observations for each residue chosen proportionally to the background distribution g. We compute the background distribution, g, by a genome-wide average over all the coding and non-coding regions. The final alignments have always gaps. For each column of the alignment, we uniformly distribute the relative frequency of gaps among all residues.
To validate our prediction method we used complete genome sequences of 39 strains of S.cer as our test data. The data, here referred to as SGRP data , includes a reference lab strain (S 288c) plus 38 other strains from different sources including other labs, pathogenic, baking, wine, food spoilage, natural fermentation, sake, probiotic, and plant isolates and from a wide range of geographic areas including North America, Europe, Malaysia, West Africa and Asia. The data were sequenced using Sanger sequencing on ABI 3730 DNA sequencers as well as Illumina Genome Analyzer. We used genome-wide alignments made by PALAS (see Supplementary data for ). For computation of derived allele frequency of indels and NM, we used a reference Saccharomyces paraduxos (S.par) as the out-group to the S.cer strains.
All DNA base reads as well as gaps in SGRP data have phred quality scores  (according to this standard, each quality character is translated into the probability of error for that read. For instance, a character "A" means that the corresponding read is wrong with the probability of error p = 10-3.2 = 0.00063). We filter out all data reads that have a quality score less than 40 (i.e. those with a probability of error more than 10-4) as well variations that occur in yeast genes annotated as dubious in SGRP. To identify indels with a better quality, we smoothed out the quality scores of the indels with their neighbouring reads (2 each side). We then marked indels with frequency ( < 2) as missing data. Then a simple heuristics algorithm was applied to remove gaps that are due to misalignments of repetitive short elements in genes. It is commonly observed that a FS indel is followed by a stop codon (created due to the frame-shift). We excluded these stop codons from our analysis of NMs and confirm that all our NMs are due to a single nucleotide NMs. In selecting mutant genes, we require that at any event of indel or NM, there is at least 10 high-quality reads available.
In our S.cer dataset, we identified 71 genes with single FS indels and 96 genes with single NM throughout the genome data for which the ancestral state were known. Due to sparse PWMs, we excluded 14 (20%) genes with FS indels and 22 (23%) genes with NMs from the analysis. We then predicted the effect of 57 genes with single frame-shifting indel and 74 genes with a single NM. Analysis of genes with more than one form of variation requires considering mutual effects of variations on function that is the subject of our near future research.
Measuring conservation in yeast non-coding DNA
For each gene, the corresponding alignment was used to estimate the evolutionary tree, T, using Felsenstein peeling algorithm . The estimated tree was then used to compute the synonymous rate of evolution λsyn using CodeML (part of PAML suite of software ) for that gene. In computing the denominator of the LLR, we used evolutionary tree, T, as well as the synonymous rates, λsyn, to compute the likelihood of each DNA site.
We then used these alignments to estimate the rate of evolution of the DNA site, , using BaseML in PAML . In computing the maximum likelihood estimate of the rate, we fixed the structure of the synonymous tree, T, and maximized the likelihood of the alignment by incrementally changing the size of the tree branch lengths while keeping their ratio constant.
The result of the algorithm is summarized as a list of indels with their corresponding DNA sites (width W = 21) ranked by their conservation (LLR score). The sites with highest LLR scores consist of fully conserved DNA bases in all W columns of 5 sequences (Figure 2).
Finding a threshold for deleterious mutations using Poisson random field
We observed that indels that fall in highly conserved DNA sites tend to segregate in lower derived allele frequencies. To test this observation quantitatively, we defined a threshold THR and classified the mutations into two classes by comparing their LLR score with the threshold. Mutations in the class with (LLR > THR) were expected to segregate in comparably lower derived allele frequencies (i.e. the spectrum of their mutations were expected to shift more towards lower allele frequencies [8, 23]).
To measure the difference between the spectrums of allele frequencies for these two classes, we fit two separate Poisson random fields  on each class of mutations and computed corresponding selection coefficients for each class using maximum likelihood estimation .
Each value of the threshold, therefore, results in two classes of mutations with two different selection coefficients. We then repeated the experiment for different values of the threshold to identify an optimal threshold that results in maximal difference between selection coefficients for the two classes.
By dividing the mutations based on their LLR scores and using the optimal threshold we obtained two classes of mutations that have maximal difference in their derived allele frequency spectrum. The maximum and minimum LLR scores in our dataset were 40.6 and 0, respectively. We estimated an optimal threshold of THR = 16.6. We used this threshold to obtain derived allele frequency spectrums as shown in (Figure 3e).
A Similar approach was used to determine the threshold for the information loss score D associated with NMs in the human population. The mutations used in this analysis had no record of disease association and therefore included mutations with lower deleterious effects (there was no genotype information for NMs with records of disease association). We divided these mutations into two classes by comparing their score D with a threshold THR. We then obtained the THR similar to what is explained in the above. We obtained THR = 0.3 for the NMs in the human population (Figure 3f)
Variations in the human population
We used all NMs and FS indels reported in the release 132 of dbSNP  for our human variations dataset. We used NCBI API to fetch the data from dbSNP. There are 33285 (8311) FS indels (NMs) reported in dbSNP from which 801 (2234) FS indels (NMs) have records of diseases association in OMIM  and LSDB . No filtering was done to ensure that the variations were unique, and we noted that ~25% of NMs and FS indels seem to be listed under multiple mutation identifiers.
We were able to compute the information loss score D for 2047 of the 2234 (91%) NMs with disease association records. We randomly chose 4468 (2234 × 2) NMs from the set of 8311 NMs with no records of disease association from which we were able to compute the D score for 3816 (85%) NMs. There was no heterozygocity information available for NMs with disease association. Therefore were able to compute the minor-allele frequency (MAF) of 1606 NMs in the dataset with no record of disease association.
For the case of FS indels, we computed the information loss score D for 696 (86%) of the 801 indels with disease association records. We randomly selected 1602 (801 × 2) indels from the set of 6317 indels with no disease association from which we were able to compute the D score for 1312 (81%) indels. There was no heterozygocity information available for the FS indels.
On average, we excluded 14% of human variants from our analysis due to sparse PWMs.
For computing the minor-allele frequency (MAF) of mutations (q), we used the heterozygocity (h) and solved the equation h = 2q(1-q). Our MAF spectrum analysis is based on the data for the NMs with no disease association record in OMIM and LSDB.
The authors acknowledge useful discussions and comments by Alex Nguyen Ba that enhanced the presented results as well as the manuscript significantly. This research is supported by Canadian Institute for Health Research grant #202372 and an infrastructure grant from the Canadian Foundation for Innovation to AMM.
A simple web application is currently available to compute information loss due to FS indels and NMs on proteins at: (http://www.moseslab.csb.utoronto.ca/indelz.html)
- Feuk L, Carson AR, Scherer SW: Structural variation in the human genome. Nature Reviews Genetics 2006, 7: 85–97.View ArticlePubMedGoogle Scholar
- Gibbs N, et al.: The international HapMap project. Nature 2003, 426: 789–796. 10.1038/nature02168View ArticleGoogle Scholar
- Dalca AV, Brudno M: Genome variation discovery with high-throughput sequencing data. Brief Bioinformmatics 2010, 11(1):3–14. 10.1093/bib/bbp058View ArticleGoogle Scholar
- Barøy T, Misceo D, Frengen E: Structural variation in the human genome contributes to variation of traits. Tidsskr Nor Laegeforen 2008, 128(17):1951–1955.PubMedGoogle Scholar
- Mullaney JM, Mills RE, Pittard WS, Devine SE: Small insertions and deletions in human genomes. Hum Mol Genet 2010. Sep 21 (PMID: 20858594) Sep 21 (PMID: 20858594)Google Scholar
- 1000 Genomes project[http://www.1000genomes.org/]
- Doniger SW, Doniger SW, Kim HS, Swain D, Corcuera D, Williams M, Yang SP, Fay JC: A catalogue of neutral and deleterious polymorphism in yeast. PLOS Genetics 2008., 4(8):Google Scholar
- Liti G, Carter DM, Moses AM, Warringer J, Parts L, James SA, Davey RP, Roberts IN, Burt A, Koufopanou V, Tsai IJ, Bergman CM, Bensasson D, O'Kelly MJ, van Oudenaarden A, Barton DB, Bailes E, Nguyen AN, Jones M, Quail MA, Goodhead I, Sims S, Smith F, Blomberg A, Durbin R, Louis EJ: Population genomics of domestic and wild yeasts. Nature 2009, 458: 337–341. 10.1038/nature07743PubMed CentralView ArticlePubMedGoogle Scholar
- Kimura M: The neutral theory of molecular evolution. Cambridge: Cambridge University Press; 1984.Google Scholar
- dbSNP build 132[http://www.ncbi.nlm.nih.gov/projects/SNP/]
- OMIM (Online Mendelian Inheritance in Man)[http://www.ncbi.nlm.nih.gov/omim]
- LSMD (Locus Specific Mutation Databases)[http://www.hgvs.org/dblist/glsdb.html]
- Pleasance ED, et al.: A comprehensive catalogue of somatic mutations from a human cancer genome. Nature 2010, 463: 191–196. 10.1038/nature08658PubMed CentralView ArticlePubMedGoogle Scholar
- Pleasance ED, et al.: A small-cell lung cancer genome with complex signatures of tobacco exposure. Nature 2010, 463: 184–190. 10.1038/nature08629PubMed CentralView ArticlePubMedGoogle Scholar
- IARC TP53 database[http://www-p53.iarc.fr/]
- The Cystic Fibrosis Genetic Analysis Consortium: Population variation of common Cystic Fibrosis mutations. Human Mutations 1994, 4: 167–177.View ArticleGoogle Scholar
- Ng PC, Henikoff S: Predicting Deleterious Amino Acid Substitutions. Genome Research 2001, 11: 863–874. 10.1101/gr.176601PubMed CentralView ArticlePubMedGoogle Scholar
- Ng PC, Henikoff S: Predicting the effects of amino acid substitutions on protein function. Human Genetics 2006, 7: 61–80.Google Scholar
- DeSalle R, Amato G: The expansion of conservation genetics. Nature Reviews Genetics 2004, 5: 702–712. 10.1038/nrg1425View ArticlePubMedGoogle Scholar
- Allendorf FW, Hohenlohe PA, Luikart G: Genomics and the future of conservation genetics. Nature Reviews Genetics 2010, 11: 697–709.View ArticlePubMedGoogle Scholar
- Schnider TD: Information content of individual genetic sequences. J of Theoretical Biology 1997, 189(4):427–441. 10.1006/jtbi.1997.0540View ArticleGoogle Scholar
- Hartl DL, Clark AG: Principles of Population Genetics. 3rd edition. Sinauer Associates; 1997.Google Scholar
- Moses AM, Sinha S: Regulatory Motif Analysis. In Bioinformatics: Tools and Applications. Edited by: D Edwards et al.. Springer Science+Business Media LLC; 2009:137–163.View ArticleGoogle Scholar
- Rasmussen SA, Friedman JM: NF1 gene and neurofibromatosis 1. American J Epidemiology 2000, 151(1):3340.View ArticleGoogle Scholar
- Konopka B, Paszko Z, Janiec-Jankowska A, Goluda M: Assessment of the quality and frequency of mutations occurrence in PTEN gene in endometrial carcinomas and hyperplasias. Cancer Letters 2002, 178: 43–51. 10.1016/S0304-3835(01)00815-1View ArticlePubMedGoogle Scholar
- Dermitzakis ET, Reymond A, Antonarakis SE: Conserved nongenic sequences- an unexpected feature of mammalian genomes. Nature Reviews Genetics 2005, 6: 151–157.View ArticlePubMedGoogle Scholar
- Drake JA, Bird C, Nemesh J, Thomas DJ, Newton-Cheh C, Reymond A, Excoffier L, Attar H, Antonarakis SE, Dermitzakis ET, Hirschhorn JN: Conserved non-coding sequences are selectively constrained and not mutation cold spots. Nature Genetics 2005, 38: 223–227.View ArticlePubMedGoogle Scholar
- Furney SJ, Higgins DG, Ouzounis CA, López-Bigas N: Structural and functional properties of genes involved in human cancer. BMC Genomics 2006., 7(3):Google Scholar
- Knight JC: Regulatory polymorphisms underlying complex disease traits. J Mol Med 2005, 83: 97–109. 10.1007/s00109-004-0603-7PubMed CentralView ArticlePubMedGoogle Scholar
- Knight JC: Functional implications of genetic variation in non-coding DNA for disease susceptibility and gene regulation. In Clinical. Volume 104. Science (London); 2003:493–501. 10.1042/CS20020304Google Scholar
- Boyadjiev SA, Jabs EW: Online Mendelian Inheritance in Man (OMIM) as a knowledgebase for human developmental disorders. Clinical Genetics 2000, 57: 253–266.View ArticlePubMedGoogle Scholar
- Knight JC: Polymorphisms in tumor necrosis factor and other cytokines as risks for infectious diseases and the septic syndrome. Current Infectious Disease Rep 2001, 3(5):427–439. 10.1007/BF03160478View ArticleGoogle Scholar
- Post SM, Quintás-Cardama A, Pant V, Iwakuma T, Hamir A, Jackson JG, Maccio DR, Bond GL, Johnson DG, Levine AJ, Lozano G: A high-frequency regulatory polymorphism in the p53 pathway accelerates tumor development. Cancer Cell 2010, 18(3):220–230. 10.1016/j.ccr.2010.07.010PubMed CentralView ArticlePubMedGoogle Scholar
- Kellis M, Patterson N, Endrizzi M, Birren B, Lander ES: Sequencing and comparison of yeast species to identify genes and regulatory elements. Nature 2003, 423(6937):241–54. 10.1038/nature01644View ArticlePubMedGoogle Scholar
- Cliften P, Sudarsanam P, Desikan A, Fulton L, Fulton B, Majors J, Waterston R, Cohen BA, Johnston M: Finding functional features in Saccharomyces genomes by phylogenetic footprinting. Science 2003, 301(5629):71–76. 10.1126/science.1084337View ArticlePubMedGoogle Scholar
- Xie X, Lu J, Kulbokas EJ, Golub TR, Mootha V, Lindblad-Toh K, Lander ES, Kellis M: Systematic discovery of regulatory motifs in human promoters and 3' UTRs by comparison of several mammals. Nature 2005, 434(7031):338–345. 10.1038/nature03441PubMed CentralView ArticlePubMedGoogle Scholar
- Ettwiller L, Paten B, Souren M, Loosli F, Wittbrodt J, Birney E: The discovery, positioning and verification of a set of transcription-associated motifs in vertebrates. Genome Biol 2005, 6(12):R104. 10.1186/gb-2005-6-12-r104PubMed CentralView ArticlePubMedGoogle Scholar
- Keightley PD, Kryukov GV, Sunyaev S, Halligan DL, Gaffney DJ: Evolutionary constraints in conserved nongenic sequences of mammals. Genome Res 2005, 15: 1373–1378. 10.1101/gr.3942005PubMed CentralView ArticlePubMedGoogle Scholar
- King DC, Taylor J, Elnitski L, Chiaromonte F, Miller W, Hardison RC: Evaluation of regulatory potential and conservation scores for detecting cis-regulatory modules in aligned mammalian genome sequences. Genome Res 2005, 15: 1051–1060. 10.1101/gr.3642605PubMed CentralView ArticlePubMedGoogle Scholar
- Bustamante CD, Fledel-Alon A, Williamson S, Nielsen R, Hubisz MT, Glanowski S, Tanenbaum DM, White TJ, Sninsky JJ, Hernandez RD, Civello D, Adams MD, Cargill M, Clark AG: Natural selection on protein-coding genes in the human genome. Nature 2005, 437: 1153–1157. 10.1038/nature04240View ArticlePubMedGoogle Scholar
- Fay JC, Wyckoff GJ, Wu CI: Positive and negative selection on the human genome. Genetics 2001, 158: 1227–1234.PubMed CentralPubMedGoogle Scholar
- Chun S, Fay JC: Identification of deleterious mutations within three human genomes. Genome Research 2009, 19: 1553–1561. 10.1101/gr.092619.109PubMed CentralView ArticlePubMedGoogle Scholar
- Yeast essential ORFs[http://www-sequence.stanford.edu/group/yeast_deletion_project/Essential_ORFs.txt]
- Warringer J, Zörgö E, Cubillos FA, Zia A, Gjuvsland A, Simpson JT, Forsmark A, Durbin R, Omholt SW, Louis EJ, Liti G, Moses A, Blomberg A: Trait variation in yeast is defined by population history. , PLoS Genetics 2011, 7(6):e1002111.PubMed CentralView ArticlePubMedGoogle Scholar
- Saccharomyces cerevisiae genome[http://www.yeastgenome.org/]
- Soussi T, Wiman KG: Shaping genetic alterations in human cancer: the p53 mutation paradigm. Cancer Cell 2007, 12(4):303–312. 10.1016/j.ccr.2007.10.001View ArticlePubMedGoogle Scholar
- Zambetti GP: The p53 mutation "gradient effect" and its clinical implications. J Cell Physiol 2007, 213(2):370–373. 10.1002/jcp.21217View ArticlePubMedGoogle Scholar
- Goode EL, et al.: A genome-wide association study identifies susceptibility loci for ovarian cancer at 2q31 and 8q24. Nature Genet 2010, 42: 874–879. 10.1038/ng.668PubMed CentralView ArticlePubMedGoogle Scholar
- Balz V, Scheckenbach K, Götte K, Bockmühl U, Petersen I, Bier H: Is the p53 inactivation frequency in squamous cell carcinomas of the head and neck underestimated? Cancer Research 2003, 63: 1188–1191.PubMedGoogle Scholar
- Sonnhammer EL, Eddy SR, Birney E, Bateman A, Durbin R: Pfam: Multiple sequence alignments and HMM-profiles of protein domains. Nucleic Acids Res 1998, 26(1):320–322. 10.1093/nar/26.1.320PubMed CentralView ArticlePubMedGoogle Scholar
- Durbin R, Eddy SR, Krogh A, Mitchison G: Biological Sequence Analysis: Probabilistic Models of Proteins and Nucleic Acids. Cambridge University Press; 1998.View ArticleGoogle Scholar
- McCarthy MI, Abecasis GR, Cardon LR, Goldstein DB, Little J, Ioannidis JP, Hirschhorn JN: Genome-wide association studies for complex traits: consensus, uncertainty and challenges. Nature Reviews Genetics 2008, 9: 356–369. 10.1038/nrg2344View ArticlePubMedGoogle Scholar
- Hirschhorn JN, Lohmueller K, Byrne E, Hirschhorn K: A comprehensive review of genetic association studies. Genet Med 2002, 4(2):45–61. 10.1097/00125817-200203000-00002View ArticlePubMedGoogle Scholar
- Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ, GappedBLAST and PSI-BLAST: A new generation of protein database search programs. Nucleic Acids Res 1997, 25: 3389–3402. 10.1093/nar/25.17.3389PubMed CentralView ArticlePubMedGoogle Scholar
- Bairoch A, Apweiler R: The SWISS-PROT protein sequence database and its supplement TrEMBL. Nucleic Acid Research 2000, 28: 45–48. 10.1093/nar/28.1.45View ArticleGoogle Scholar
- Phylogenetic Analysis by Maximum Likelihood[http://abacus.gene.ucl.ac.uk/software/paml.html]
- Sawyer SA, Hartl DL: Population genetics of polymorphism and divergence. Genetics 1992, 132: 1161–1176.PubMed CentralPubMedGoogle Scholar
- Bustamante CD, Wakeley J, Sawyer S, Hartl DL: Directional selection and the site-frequency spectrum. Genetics 2001., 159: 1779:1788 1779:1788Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.