SNAP predicts effect of mutations on protein function

Many non-synonymous single nucleotide polymorphisms (nsSNPs) in humans are suspected to impact protein function. Here, we present a publicly available server implementation of the method SNAP (screening for non-acceptable polymorphisms) that predicts the functional effects of single amino acid substitutions. SNAP identifies over 80% of the non-neutral mutations at 77% accuracy and over 76% of the neutral mutations at 80% accuracy at its default threshold. Each prediction is associated with a reliability index that correlates with accuracy and thereby enables experimentalists to zoom into the most promising predictions.     

Selective pressures at a codon-level predict deleterious mutations in human disease genes

Deleterious mutations affecting biological function of proteins are constantly being rejected by purifying selection from the gene pool. The non-synonymous/synonymous substitution rate ratio (omega) is a measure of selective pressure on amino acid replacement mutations for protein-coding genes. Different methods have been developed in order to predict non-synonymous changes affecting gene function. However, none has considered the estimation of selective constraints acting on protein residues. Here, we have used codon-based maximum likelihood models in order to estimate the selective pressures on the individual amino acid residues of a well-known model protein: p53. We demonstrate that the number of residues under strong purifying selection in p53 is much higher than those that are strictly conserved during the evolution of the species. In agreement with theoretical expectations, residues that have been noted to be of structural relevance, or in direct association with DNA, were among those showing the highest signals of purifying selection. Conversely, those changing according to a neutral, or nearly neutral mode of evolution, were observed to be irrelevant for protein function. Finally, using more than 40 human disease genes, we demonstrate that residues evolving under strong selective pressures (omega<0.1) are significantly associated (p<0.01) with human disease. We hypothesize that non-synonymous change on amino acids showing omega<0.1 will most likely affect protein function. The application of this evolutionary prediction at a genomic scale will provide an a priori hypothesis of the phenotypic effect of non-synonymous coding single nucleotide polymorphisms (SNPs) in the human genome.     

Proteome-wide analysis of single-nucleotide variations in the N-glycosylation sequon of human genes

N-linked glycosylation is one of the most frequent post-translational modifications of proteins with a profound impact on their biological function. Besides other functions, N-linked glycosylation assists in protein folding, determines protein orientation at the cell surface, or protects proteins from proteases. The N-linked glycans attach to asparagines in the sequence context Asn-X-Ser/Thr, where X is any amino acid except proline. Any variation (e.g. non-synonymous single nucleotide polymorphism or mutation) that abolishes the N-glycosylation sequence motif will lead to the loss of a glycosylation site. On the other hand, variations causing a substitution that creates a new N-glycosylation sequence motif can result in the gain of glycosylation. Although the general importance of glycosylation is well known and acknowledged, the effect of variation on the actual glycoproteome of an organism is still mostly unknown. In this study, we focus on a comprehensive analysis of non-synonymous single nucleotide variations (nsSNV) that lead to either loss or gain of the N-glycosylation motif. We find that 1091 proteins have modified N-glycosylation sequons due to nsSNVs in the genome. Based on analysis of proteins that have a solved 3D structure at the site of variation, we find that 48% of the variations that lead to changes in glycosylation sites occur at the loop and bend regions of the proteins. Pathway and function enrichment analysis show that a significant number of proteins that gained or lost the glycosylation motif are involved in kinase activity, immune response, and blood coagulation. A structure-function analysis of a blood coagulation protein, antithrombin III and a protease, cathepsin D, showcases how a comprehensive study followed by structural analysis can help better understand the functional impact of the nsSNVs.     

Identification and analysis of deleterious human SNPs

We have developed two methods of identifying which non-synonomous single base changes have a deleterious effect on protein function in vivo. One method, described elsewhere, analyzes the effect of the resulting amino acid change on protein stability, utilizing structural information. The other method, introduced here, makes use of the conservation and type of residues observed at a base change position within a protein family. A machine learning technique, the support vector machine, is trained on single amino acid changes that cause monogenic disease, with a control set of amino acid changes fixed between species. Both methods are used to identify deleterious single nucleotide polymorphisms (SNPs) in the human population. After carefully controlling for errors, we find that approximately one quarter of known non-synonymous SNPs are deleterious by these criteria, providing a set of possible contributors to human complex disease traits.     

SNPdbe: constructing an nsSNP functional impacts database

Many existing databases annotate experimentally characterized single nucleotide polymorphisms (SNPs). Each non-synonymous SNP (nsSNP) changes one amino acid in the gene product (single amino acid substitution;SAAS). This change can either affect protein function or be neutral in that respect. Most polymorphisms lack experimental annotation of their functional impact. Here, we introduce SNPdbe-SNP database of effects, with predictions of computationally annotated functional impacts of SNPs. Database entries represent nsSNPs in dbSNP and 1000 Genomes collection, as well as variants from UniProt and PMD. SAASs come from >2600 organisms; 'human' being the most prevalent. The impact of each SAAS on protein function is predicted using the SNAP and SIFT algorithms and augmented with experimentally derived function/structure information and disease associations from PMD, OMIM and UniProt. SNPdbe is consistently updated and easily augmented with new sources of information. The database is available as an MySQL dump and via a web front end that allows searches with any combination of organism names, sequences and mutation IDs. AVAILABILITY: http://www.rostlab.org/services/snpdbe.     

Bison PRNP genotyping and potential association with Brucella spp. seroprevalence

The implication that host cellular prion protein (PrP(C)) may function as a cell surface receptor and/or portal protein for Brucella abortus in mice prompted an evaluation of nucleotide and amino acid variation within exon 3 of the prion protein gene (PRNP) for six US bison populations. A non-synonymous single nucleotide polymorphism (T50C), resulting in the predicted amino acid replacement M17T (Met --> Thr), was identified in each population. To date, no variation (T50; Met) has been detected at the corresponding exon 3 nucleotide and/or amino acid position for domestic cattle. Notably, 80% (20 of 25) of the Yellowstone National Park bison possessing the C/C genotype were Brucella spp. seropositive, representing a significant (P = 0.021) association between seropositivity and the C/C genotypic class. Moreover, significant differences in the distribution of PRNP exon 3 alleles and genotypes were detected between Yellowstone National Park bison and three bison populations that were either founded from seronegative stock or previously subjected to test-and-slaughter management to eradicate brucellosis. Unlike domestic cattle, no indel polymorphisms were detected within the corresponding regions of the putative bison PRNP promoter, intron 1, octapeptide repeat region or 3'-untranslated region for any population examined. This study provides the first evidence of a potential association between nucleotide variation within PRNP exon 3 and the presence of Brucella spp. antibodies in bison, implicating PrP(C) in the natural resistance of bison to brucellosis infection.     

Clinal variation for amino acid polymorphisms at the Pgm locus in Drosophila melanogaster

Clinal variation is common for enzymes in the glycolytic pathway for Drosophila melanogaster and is generally accepted as an adaptive response to different climates. Although the enzyme phosphoglucomutase (PGM) possesses several allozyme polymorphisms, it is unique in that it had been reported to show no clinal variation. Our recent DNA sequence investigation of Pgm found extensive cryptic amino acid polymorphism segregating with the allozyme alleles. In this study, we characterize the geographic variation of Pgm amino acid polymorphisms at the nucleotide level along a latitudinal cline in the eastern United States. A survey of 15 SNPs across the Pgm gene finds significant clinal differentiation for the allozyme polymorphisms as well as for many of the cryptic amino acid polymorphisms. A test of independence shows that pervasive linkage disequilibrium across this gene region can explain many of the amino acid clines. A single Pgm haplotype defined by two amino acid polymorphisms shows the strongest correlation with latitude and the steepest change in allele frequency across the cline. We propose that clinal selection at Pgm may in part explain the extensive amino acid polymorphism at this locus and is consistent with a multilocus response to selection in the glycolytic pathway.     

Unusual signatures of highly adaptable R-loci in closely-related Arabidopsis species

Plant resistance genes (R-genes) evolve rapidly in response to changing environments. What are the most remarkable signatures of fast adaptive genes, besides the commonly revealed rapid divergence and high non-synonymous substitution rate? Here we investigated these changes in five R-loci following recent differentiation between Arabidopsis thaliana and Arabidopsis lyrata. Extreme differences in evolutionary rates were observed: e.g., an overall 5.46-9.83-fold different nucleotide diversity at two R-loci between species, ten-fold higher non-synonymous substitution rates within one species versus the other, significantly different Ka/Ks ratios between species for the same R-gene, and high interspecific divergence at one R-locus. Particularly, we observed an elevated level of trans-specific polymorphism at one R-locus and a differentially maintained presence/absence polymorphism at another. The high frequency of ancestral polymorphisms amongst R-genes suggests that the persistence of some functional variation is an important evolutionary mechanism shaping genetic variation in R-genes, while the variation of presence/absence polymorphisms provides a potential mechanism for malleable activation of adaptive resistance response pathways. The distinct patterns among R-genes suggest that the same R-gene ortholog can be quickly shaped by different evolutionary processes, e.g., purifying selection in one species but positive selection in a closely-related species.     

Origins of functional nucleotide polymorphisms in a major quantitative trait locus, qLTG3-1, controlling low-temperature germinability in rice

qLTG3-1 is a major quantitative trait locus (QTL) controlling tolerance to low-temperature at the seed germination stage (termed low-temperature germinability) in rice using a population derived from the cross between Italica Livorno from Italy and Hayamasari from Japan. Map-based cloning identified that qLTG3-1 encodes a protein of unknown function. The molecular identification of this major QTL could make it possible to identify allelic variation and favorable alleles for rice breeding programs. The present study examined the identification of qLTG3-1 alleles and their distribution among 62 landraces of Asian cultivated rice (Oryza sativa L.) collected from 19 different countries, termed the rice core collection. In the coding region, a single non-synonymous substitution and 3 in-frame insertion/deletion polymorphisms (indels) were detected. The almost completely conserved protein alignment of qLTG3-1 was also identified among 5 Oryza species, suggesting that the function of qLTG3-1 is critical for seed germination or for rice growth by pleiotropic effects of the gene. The functional nucleotide polymorphisms (FNPs), a 71-bp deletion found in Hayamasari and an amino acid substitution found in Nipponbare, was identified in varieties from Japan. These alleles with FNPs might be adapted to rice cultivation in specific local conditions. The present results may contribute to the utilization of favorable alleles of qLTG3-1 for the improvement of low-temperature germinability in rice breeding programs.     

Imputation of missing genotypes: an empirical evaluation of IMPUTE

Background

Polymorphism in genes of regulating enzymes, transporters and receptors of the neurotransmitters of the central nervous system have been associated with altered behaviour, and single nucleotide polymorphisms (SNPs) represent the most frequent type of genetic variation. The serotonin and dopamine signalling systems have a central influence on different behavioural phenotypes, both of invertebrates and vertebrates, and this study was undertaken in order to explore genetic variation that may be associated with variation in behaviour.

Results

Single nucleotide polymorphisms in canine genes related to behaviour were identified by individually sequencing eight dogs (Canis familiaris) of different breeds. Eighteen genes from the dopamine and the serotonin systems were screened, revealing 34 SNPs distributed in 14 of the 18 selected genes. A total of 24,895 bp coding sequence was sequenced yielding an average frequency of one SNP per 732 bp (1/732). A total of 11 non-synonymous SNPs (nsSNPs), which may be involved in alteration of protein function, were detected. Of these 11 nsSNPs, six resulted in a substitution of amino acid residue with concomitant change in structural parameters.

Conclusion

We have identified a number of coding SNPs in behaviour-related genes, several of which change the amino acids of the proteins. Some of the canine SNPs exist in codons that are evolutionary conserved between five compared species, and predictions indicate that they may have a functional effect on the protein. The reported coding SNP frequency of the studied genes falls within the range of SNP frequencies reported earlier in the dog and other mammalian species. Novel SNPs are presented and the results show a significant genetic variation in expressed sequences in this group of genes. The results can contribute to an improved understanding of the genetics of behaviour.     

Novel multi-nucleotide polymorphisms in the human genome characterized by whole genome and exome sequencing

Genomic sequence comparisons between individuals are usually restricted to the analysis of single nucleotide polymorphisms (SNPs). While the interrogation of SNPs is efficient, they are not the only form of divergence between genomes. In this report, we expand the scope of polymorphism detection by investigating the occurrence of double nucleotide polymorphisms (DNPs) and triple nucleotide polymorphisms (TNPs), in which two or three consecutive nucleotides are altered compared to the reference sequence. We have found such DNPs and TNPs throughout two complete genomes and eight exomes. Within exons, these novel polymorphisms are over-represented amongst protein-altering variants; nearly all DNPs and TNPs result in a change in amino acid sequence and, in some cases, two adjacent amino acids are changed. DNPs and TNPs represent a potentially important new source of genetic variation which may underlie human disease and they should be included in future medical genetics studies. As a confirmation of the damaging nature of xNPs, we have identified changes in the exome of a glioblastoma cell line that are important in glioblastoma pathogenesis. We have found a TNP causing a single amino acid change in LAMC2 and a TNP causing a truncation of HUWE1.     

Comparison of the genomic structure and variation in the two human sodium-dependent vitamin C transporters, <Emphasis Type="Italic">SLC23A1</Emphasis> and <Emphasis Type="Italic">SLC23A2</Emphasis>

Vitamin C (L-ascorbic acid) is an essential co-factor for eight mammalian enzymes and quenches reactive oxygen species. Sodium-dependent vitamin C transport is mediated by two transporters, SVCT 1 and SVCT 2, encoded by SLC23A1 and SLC23A2. We characterized the genomic structures of SLC23A1 and SLC23A2, determined the extent of genetic variation and linkage disequilibrium across each gene, analyzed nucleotide diversity to estimate the effect of selective pressure, and compared sequence variation across species. In SLC23A1, the majority of single nucleotide polymorphisms (SNPs) are population-specific in either African Americans or Caucasians, including three of four non-synonymous SNPs. In contrast, most SNPs in SLC23A2 are shared between African Americans and Caucasians, and there are no non-synonymous SNPs in SLC23A2. Our analysis, combined with previous in vitro and in vivo studies, suggests that non-synonymous variation appears to be tolerated in SLC23A1 but not SLC23A2, and that this may be a consequence of different selective pressures following past gene duplication of the sodium-dependent vitamin C transporters. Genetic association studies of these two genes will need to account for the differences in haplotype structure and the population-specific variants. Our data represent a fundamental step toward the application of genetics to refining nutrient recommendations, specifically for vitamin C, and may serve as a paradigm for other vitamins.     

Molecular evolution of bovine Toll-like receptor 2 suggests substitutions of functional relevance

Background  

There is accumulating evidence that polymorphism in Toll-like receptor (TLR) genes might be associated with disease resistance or susceptibility traits in livestock. Polymorphic sites affecting TLR function should exhibit signatures of positive selection, identified as a high ratio of non-synonymous to synonymous nucleotide substitutions (ω). Phylogeny based models of codon substitution based on estimates of ω for each amino acid position can therefore offer a valuable tool to predict sites of functional relevance. We have used this approach to identify such polymorphic sites within the bovine TLR2 genes from ten Bos indicus and Bos taurus cattle breeds. By analysing TLR2 gene phylogeny in a set of mammalian species and a subset of ruminant species we have estimated the selective pressure on individual sites and domains and identified polymorphisms at sites of putative functional importance.     

In vitro characterization of a human calcitonin receptor gene polymorphism

Calcitonin is a 32 amino acid peptide hormone that inhibits bone resorption by stimulating calcitonin receptors (CTR) located on the surfaces of osteoclasts. A polymorphism at nucleotide 1340 of the human calcitonin receptor gene (CALCR) lies within the coding region and has the potential to change the amino acid at codon 447 from leucine to proline. In the present study, we scanned the coding region, portions of the 5'-flanking and 3'-flanking sequences, and the intron-exon boundaries of the human CALCR gene for additional polymorphisms, and determined the frequency of the codon 447 polymorphism in several ethnic groups. Because a leucine to proline change has the potential for significant structural alteration, receptor genes encoding either leucine or proline at residue 447 were transiently expressed in COS-7 cells to determine the binding and functional consequences of this polymorphism. Our complete polymorphism scan of the CALCR gene identified 11 polymorphic sites in the gene and confirmed the presence of the previously identified nucleotide T1340C (codon 447) polymorphism. Ten of the 11 polymorphisms were single nucleotide polymorphisms (SNPs). For the codon 447 polymorphism, the prevalence of the TT genotype (leucine) was 59% in Caucasians, 27% in African-Americans, 0% in Asians, and 20% in Hispanics. The presence of this SNP appears to have no statistically significant difference with the receptor's ability to bind calcitonin or signal when activated with the hormone.     

In silico analysis of single nucleotide polymorphism (SNPs) in human β-globin gene

Single amino acid substitutions in the globin chain are the most common forms of genetic variations that produce hemoglobinopathies--the most widespread inherited disorders worldwide. Several hemoglobinopathies result from homozygosity or compound heterozygosity to beta-globin (HBB) gene mutations, such as that producing sickle cell hemoglobin (HbS), HbC, HbD and HbE. Several of these mutations are deleterious and result in moderate to severe hemolytic anemia, with associated complications, requiring lifelong care and management. Even though many hemoglobinopathies result from single amino acid changes producing similar structural abnormalities, there are functional differences in the generated variants. Using in silico methods, we examined the genetic variations that can alter the expression and function of the HBB gene. Using a sequence homology-based Sorting Intolerant from Tolerant (SIFT) server we have searched for the SNPs, which showed that 200 (80%) non-synonymous polymorphism were found to be deleterious. The structure-based method via PolyPhen server indicated that 135 (40%) non-synonymous polymorphism may modify protein function and structure. The Pupa Suite software showed that the SNPs will have a phenotypic consequence on the structure and function of the altered protein. Structure analysis was performed on the key mutations that occur in the native protein coded by the HBB gene that causes hemoglobinopathies such as: HbC (E→K), HbD (E→Q), HbE (E→K) and HbS (E→V). Atomic Non-Local Environment Assessment (ANOLEA), Yet Another Scientific Artificial Reality Application (YASARA), CHARMM-GUI webserver for macromolecular dynamics and mechanics, and Normal Mode Analysis, Deformation and Refinement (NOMAD-Ref) of Gromacs server were used to perform molecular dynamics simulations and energy minimization calculations on β-Chain residue of the HBB gene before and after mutation. Furthermore, in the native and altered protein models, amino acid residues were determined and secondary structures were observed for solvent accessibility to confirm the protein stability. The functional study in this investigation may be a good model for additional future studies.     

Positive and negative selection on the human genome.

The distinction between deleterious, neutral, and adaptive mutations is a fundamental problem in the study of molecular evolution. Two significant quantities are the fraction of DNA variation in natural populations that is deleterious and destined to be eliminated and the fraction of fixed differences between species driven by positive Darwinian selection. We estimate these quantities using the large number of human genes for which there are polymorphism and divergence data. The fraction of amino acid mutations that is neutral is estimated to be 0.20 from the ratio of common amino acid (A) to synonymous (S) single nucleotide polymorphisms (SNPs) at frequencies of > or =15%. Among the 80% of amino acid mutations that are deleterious at least 20% of them are only slightly deleterious and often attain frequencies of 1-10%. We estimate that these slightly deleterious mutations comprise at least 3% of amino acid SNPs in the average individual or at least 300 per diploid genome. This estimate is not sensitive to human population history. The A/S ratio of fixed differences is greater than that of common SNPs and suggests that a large fraction of protein divergence is adaptive and driven by positive Darwinian selection.     

Molecular anatomy of the cytoplasmic domain of bovine growth hormone receptor, a quantitative trait locus

Quantitative trait loci (QTL) studies have indicated growth hormone receptor (GHR) as a candidate gene affecting cattle milk yield and composition. In order to characterize genetic variation at GHR in cattle, we studied European and East African breeds with different histories of selection, and Bos grunniens, Ovis aries, Sus scrofa, Bison bison and Rangifer tarandus as references. We sequenced most of the cytoplasmic domain (900 bp of exon 10), 89 bp of exon 8, including the putative causative mutation for the QTL effect, and 390 bp of intron 8 for comparison. In the cytoplasmic domain, seven synonymous and seven non-synonymous single nucleotide polymorphisms (SNP) were identified in cattle. Three non-synonymous SNPs were found in sheep and one synonymous SNP in yak, while other studied species were monomorphic. Three major haplotypes were observed, one unique to African breeds, one unique to European breeds and one shared. Bison and yak haplotypes are derivatives of the European haplotype lineage. Most of the exon 10 non-synonymous cattle SNPs appear at phylogenetically highly conserved sites. The polymorphisms in exon 10 cluster around a ruminant-specific tyrosine residue, suggesting that this site may act as an additional signalling domain of GHR in ruminants. Alternative explanations for the persistent polymorphism include balancing selection, hitch-hiking, pleiotropic or sexually antagonistic fitness effects or relaxed functional constraints.     

Population variation in linkage disequilibrium across the COMT gene considering promoter region and coding region variation

Catechol-O-methyl transferase (COMT) catalyzes the first step in one of the major pathways in the degradation of catecholamines. The COMT gene on chromosome 22 has been considered a candidate gene for many neuropsychiatric disorders, in part because an exon 4 single nucleotide polymorphism (SNP) in COMT causes an amino acid substitution associated with significantly altered enzyme activity. This functional variant, detected as an NlaIII restriction site polymorphism (RSP), is polymorphic in populations from around the world. A four-site haplotype spanning 28 kb effectively encompasses COMT. This haplotype is comprised of two novel polymorphisms [a tetranucleotide short tandem repeat polymorphism (STRP) in intron 1 and a HindIII RSP at the 5' end of COMT], the NlaIII site, and another previously published site - a BglI RSP at the 3' end of the gene. Overall linkage disequilibrium (LD) for this haplotype is strong and significant in 32 population samples from around the world. Conditional probabilities indicate that, in spite of moderate to strong disequilibrium in most non-African populations, the NlaIII site, although often used for prediction, would not always be a reliable predictor of allelic variation at the other sites. Because other functional variation might exist, especially regulatory variation, these findings indicate that haplotypes would be more effective indicators of possible involvement of COMT in disease etiology.     

Genome bioinformatic analysis of nonsynonymous SNPs

Background  

Genome-wide association studies of common diseases for common, low penetrance causal variants are underway. A proportion of these will alter protein sequences, the most common of which is the non-synonymous single nucleotide polymorphism (nsSNP). It would be an advantage if the functional effects of an nsSNP on protein structure and function could be predicted, both for the final identification process of a causal variant in a disease-associated chromosome region, and in further functional analyses of the nsSNP and its disease-associated protein.