Skip Navigation


Nucleic Acids Research Advance Access originally published online on December 17, 2007
Nucleic Acids Research 2008 36(Database issue):D825-D829; doi:10.1093/nar/gkm979
This Article
Right arrow Abstract Freely available
Right arrow Print PDF (1560K) Freely available
Right arrow Screen PDF (267K) Freely available
Right arrow Supplementary Data
Right arrowOA All Versions of this Article:
36/suppl_1/D825    most recent
gkm979v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Commercial Re-use Guidelines
for Open Access NAR Content
Google Scholar
Right arrow Articles by Reumers, J.
Right arrow Articles by Schymkowitz, J.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Reumers, J.
Right arrow Articles by Schymkowitz, J.
Social Bookmarking
Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

Nucleic Acids Research, 2008, Vol. 36, Database issue D825-D829
© 2007 The Author(s)
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/2.0/uk/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article appears in the following Nucleic Acids Research issue: Database issue [View the issue table of contents]

Articles

Joint annotation of coding and non-coding single nucleotide polymorphisms and mutations in the SNPeffect and PupaSuite databases

Joke Reumers1,2, Lucia Conde3, Ignacio Medina3, Sebastian Maurer-Stroh1,2, Joost Van Durme1,2, Joaquin Dopazo3,4, Frederic Rousseau1,2 and Joost Schymkowitz1,2,*

1Switch Laboratory, Department of Applied Biological Sciences, Vrije Universiteit Brussel, 2Switch Laboratory, VIB, Pleinlaan 2, 1050 Brussel, Belgium, 3Bioinformatics Department, Centro de Investigación Príncipe Felipe (CIPF) and 4Functional Genomics Node, INB, CIPF, Valencia 46013, Spain

*To whom correspondence should be addressed. Tel: +32 2 629 14 25; Fax: +32 2 629 19 63; E-mail: joost.schymkowitz{at}vub.ac.be

Correspondence may also be addressed to Frederic Rousseau. Tel: +32 2 629 14 25; Fax: +32 2 629 19 63; E-mail: frederic.rousseau{at}vub.ac.be

Received September 15, 2007. Revised October 17, 2007. Accepted October 18, 2007.


   ABSTRACT
TOP
 ABSTRACT
INTRODUCTION
 ADDITIONS TO SNP ANNOTATION
 DISCUSSION
 SUPPLEMENTARY DATA
 REFERENCES
 
Single nucleotide polymorphisms (SNPs) are, together with copy number variation, the primary source of variation in the human genome. SNPs are associated with altered response to drug treatment, susceptibility to disease and other phenotypic variation. Furthermore, during genetic screens for disease-associated mutations in groups of patients and control individuals, the distinction between disease causing mutation and polymorphism is often unclear. Annotation of the functional and structural implications of single nucleotide changes thus provides valuable information to interpret and guide experiments. The SNPeffect and PupaSuite databases are now synchronized to deliver annotations for both non-coding and coding SNP, as well as annotations for the SwissProt set of human disease mutations. In addition, SNPeffect now contains predictions of Tango2: an improved aggregation detector, and Waltz: a novel predictor of amyloid-forming sequences, as well as improved predictors for regions that are recognized by the Hsp70 family of chaperones. The new PupaSuite version incorporates predictions for SNPs in silencers and miRNAs including their targets, as well as additional methods for predicting SNPs in TFBSs and splice sites. Also predictions for mouse and rat genomes have been added. In addition, a PupaSuite web service has been developed to enable data access, programmatically. The combined database holds annotations for 4 965 073 regulatory as well as 133 505 coding human SNPs and 14 935 disease mutations, and phenotypic descriptions of 43 797 human proteins and is accessible via http://snpeffect.vib.be and http://pupasuite.bioinfo.cipf.es/.


   INTRODUCTION
TOP
 ABSTRACT
 INTRODUCTION
ADDITIONS TO SNP ANNOTATION
 DISCUSSION
 SUPPLEMENTARY DATA
 REFERENCES
 
With the completion of the sequencing of the human genome, much attention has been centered on the study of human genome variability. Single nucleotide polymorphisms (SNPs) are the most common source of human genetic variation, and are undoubtedly a valuable resource for investigating the genetic basis of diseases. SNPs, together with DNA copy number variations (CNVs), have become one of the most actively researched areas of genomics in recent years (1,2).

In this article, we focus on the smaller of these two types of variation, SNPs. Single nucleotide polymorphisms are highly abundant, stable and distributed throughout the genome. This type of variation is associated with diversity in the population, individuality, and although the majority of these variations probably result in neutral phenotypic outcomes, certain polymorphisms can predispose individuals to disease, or influence its severity, progression or individual response to medicine. Viewed at the molecular level, these functional SNPs can affect the human phenotype by interfering on both levels of the protein synthesis machinery: non-coding SNPs may disrupt transcription factor binding sites, splice sites and other functional sites on the transcriptional level, whereas coding SNPs can cause an amino acid change and alter the functional or structural properties of the translated protein. Annotating the way a polymorphism affects an individual's phenotype should therefore focus on both levels: describing the effects on the gene's and the protein's properties. To this end, we merged two databases that each focus on one of these levels: PupaSuite (3) and SNPeffect (4,5).

PupaSuite is a web tool for selecting SNPs with potential phenotypic effect and was originally based on the combined functionality of the PupaSNP (6) and PupasView (7) web tools. Since its release in 2006, PupaSuite was extended with new tools for predicting silencers and enhancers, new conservation measures and predictions on mouse and rat genomes. SNPeffect focuses on the functional annotation of non-synonymous coding SNPs in human proteomes, but now also includes predictions on known disease mutations from the UniProt Knowledge Base.


   ADDITIONS TO SNP ANNOTATION
TOP
 ABSTRACT
 INTRODUCTION
 ADDITIONS TO SNP ANNOTATION
DISCUSSION
 SUPPLEMENTARY DATA
 REFERENCES
 
SNPeffect 3.0 novelties
SNPeffect focuses on the molecular phenotypes of variants, and includes details on structural, functional and cellular effects on the amino acid change of a non-synonymous coding SNP. Destabilizing variations which affect the aggregation behavior of a protein can be potentially disease causing and therefore are of particular interest to researchers. By investigating polymorphisms on a large scale and providing their phenotypic effects together with the ability to filter out specific functional or structural changes, researchers are able to select potentially interesting polymorphisms for further analysis.

Waltz and Tango, two aggregation predictors based on physical properties
The Waltz algorithm combines sequence, physical properties and structural parameters to identify motifs that can nucleate amyloid fiber formation in proteins (Maurer-Stroh et al., submitted for publication). Special emphasis is made to minimize overprediction of amorphous beta aggregation compared to the highly regular cross-beta structure of amyloid fibrils. Creation of new amyloidogenic motifs through nsSNPs is implicated in amyloid deposit diseases. The statistical mechanics algorithm Tango, on the other hand, predicts protein regions that are prone to form amorphous beta aggregates. Included in Tango2 is an improved prediction of mutation effects on aggregation, transmembrane region stability and signal peptide disruption. Both Waltz amylogenic regions and Tango aggregating regions are available as annotations via ProteinDAS servers at the DAS registry at the European Bioinformatics Institute.

Hsp70 chaperone family-binding predictor
A DnaK-binding site predictor was built using a dual method combining sequence and structural information. Experimental DnaK-binding data of 53 non-redundant peptide sequences allowed us to generate a sequence-based position-specific scoring matrix (PSSM) based on logarithm of the odds scores. Following an in silico alanine scan of the substrate peptide in the crystal structure of a DnaK-substrate complex [1dkx, (8)] using the FoldX force field (9), we generated a structure-based PSSM that reflects the individual contribution of certain substrate residue types for DnaK binding. Upon adding the structure-based PSSM with a normalization factor of 0.2 to the sequence-based PSSM, we obtained a DnaK motif predictor that was able to correctly predict 89% of the true positives in the tested peptide set (high sensitivity), with a concurrent amount of only 5.9% false positives for a specific score threshold (high specificity). To assess the robustness of the predictor, we carried out a cross-validation by leaving out each sequence from the learning set together with its close homologs and calculating the rate of repredicting them. This resulted in a prediction accuracy of 72% true positives and 5.9% false positives. The predictor was able to identify an entire known DnaK-binding site in the heat-shock promoter {sigma}32(10).

UniProt disease mutations
The new data source included in SNPeffect, the human dataset available in the UniProt knowledge base (11) version 52.0 (March 2007), enables us to show results on a set of known disease mutations (Table 1). Functional and structural annotations of known disease mutations are of particular interest as they can help direct experimental setup for the elucidation of the molecular mechanism of disease.


View this table:
[in this window]
[in a new window]

  		Table 1. Statistics on putative deleterious SNPs and the affected phenotypic property as annotated by the SNPeffect and PupaSuite databases, and SNPeffect annotation results for known human disease mutations from the Uniprot Knowledge Base release 52.0

 
Additional protein level annotations
Several protein functional annotations were added via the DAS registry at the EBI. A detailed list of annotations provided through a ProteinDAS service is listed in Supplementary Table 2. Additional sequence and structure based tools and databases used to describe variations and proteins in the SNPeffect dataset are listed in Supplementary Table 3.

PupaSuite enhancements
While much attention has been focused on the effects of variation on the amino acid sequence, variations that disrupt gene regulation, expression or splicing can dramatically impact gene function. PupaSuite focuses mainly on the possible effect of these regulatory variations. In this new version of PupaSuite, the database has been updated to analyze the complete set of SNPs cataloged in version 44 of Ensembl (12), which includes dbSNP (13) 126 genotype data and Sanger-caller Celera SNPs. Together with the prediction methods already included, some novel features have been added in this release.

Exonic splicing silencers (ESS)
ESSs are cis-regulatory elements located in coding regions that inhibit the use of adjacent splice sites, often contributing to alternative splicing. Wang et al. (14) described a list of 103 hexamers (the FAS-hex-3 set) identified as ESS candidates by genetic selection; we scanned the exon sequences of all the human genes to identify putative ESSs from Wang's set. SNPs located at these motifs are cataloged as potential SNPs that could disturb the silencer activity. To make prediction more reliable, the tool allows the search to be done in conserved regions.

Transcription factor binding sites
A complementary approach for TFBS identification has been included, which uses the position weight matrices (PWM) deposited in JASPAR (15). JASPAR is an open-access database of annotated, high-quality, matrix-based transcription factor binding site profiles for multicellular eukaryotes. It contains models derived from 111 profiles that were exclusively derived from published collections of experimentally defined TFBSs for multicellular eukaryotes. We use the matrices corresponding to vertebrates to search for TFBSs in the 5 kb upstream region of all the human genes. To this end, we use MatScan (http://genome.imim.es), a program to search binding sites in genomic sequences. Since MatScan does not allow a cutoff to minimize false positives, we also use the Meta program (http://genome.imim.es) to filter the results by searching the coincidences of TFBSs in orthologous genes in mouse.

Prediction of new splice sites
Gene ID (16) is a program to predict genes in genomic sequences, where splice sites, and start and stop codons are predicted and scored along the sequence using PWMs. We use this program to scan the whole genome to find new splice sites and to map SNPs that could have a putative effect in the disruption of these important sites.

MicroRNAs (miRNAs)
miRNAs act as repressors of protein-coding genes by binding to target sites in the 3' UTR of mRNAs. In this release, we scan the genome to find SNPs located at miRNAs. Besides, we use miRanda (14), an algorithm for the detection of potential microRNA target sites in genomic sequences, to localize all the SNPs situated in the region 3' UTR of these targets sites. Both SNPs at miRNAs and SNPs in their target sequences could have an effect in the normal function of these regulatory elements. This effect is measured by the difference of scores among the alleles of the SNPs.

Mouse and rat genomes
Finally, this new release of PupaSuite incorporates the analysis of genetic variations for mouse and rat genomes. Because most of the methods include PWMs for vertebrates, most of the predictions can be extrapolated and the search of regulatory elements can be done in these genomes. Predictions for exonic splicing enhancers and silencers are not extrapolated since the proteins used for building the PWMs correspond exclusively to human proteins. With the inclusion of this information, the tool can aid to better understand the functional diversity in different genomes.

PupaSuite web service
In addition to the updated web page interface, a set of public web services, implemented in Java, have also been developed. These web services constitute a complete and exhaustive API to access all the functional data showed in the web page

An example of how to use the PupaSuite web services is included in the Supplementary Data.

Availability of the databases
The merge of the SNPeffect and PupaSuite databases includes accessibility of the combined data through both portals (at: http://snpeffect.vib.be and http://pupasuite.bioinfo.cipf.es/). PupaSuite data are accessible through the SNPeffect interface and vice versa (Figure 1). Both databases are freely available for academic users. Help pages on both servers provide detailed information on usage of the interfaces. SNPeffect and PupaSuite will be updated with each even version release of Ensembl, which corresponds to a four-monthly update.


Figure 1
View larger version (40K):
[in this window]
[in a new window]
[Download PowerPoint slide]
  		Figure 1. Viewing PupaSuite annotation results in SNPeffect and vice versa.

 

   DISCUSSION
TOP
 ABSTRACT
 INTRODUCTION
 ADDITIONS TO SNP ANNOTATION
 DISCUSSION
SUPPLEMENTARY DATA
 REFERENCES
 
Nowadays, more than 11 million SNPs have been described in databases like dbSNP. Among them, thousands of SNPs can have a direct impact on disease. Recently, different bioinformatics tools have been developed which try to find these putative disruptive polymorphisms. These tools use different information based on sequence, structure, conservation or functional properties to distinguish disease-causing mutations from those that are thought to have a neutral effect. Because of its importance on biomedical research, it would be beneficial to generate bioinformatics tools to extract and merge interesting data coming from these heterogeneous sources, and collect these results in a single database. SNPeffect and PupaSuite are two of the most complete bioinformatics tools for the analysis of SNPs. Both tools integrate different methods for the analysis of SNPs, focusing on different levels of the protein synthesis machinery: (i) SNPeffect focuses on SNPs in gene-coding regions that can lead to changes in the biological properties of the encoded protein, (ii) PupaSuite focuses on SNPs in non-coding gene regulatory regions which may affect gene expression levels and mRNA stability. In this article, we present a joint effort for the integration of both tools, which has lead to the creation of a comprehensive database of putative functional polymorphisms. This way 749 603 regulatory human SNPs and 31 415 non-synonymous coding human SNPs were annotated as putative disruptive polymorphisms, and the same procedure resulted in the suggestion of a molecular mechanism of disease for 3660 known human disease mutations.


   SUPPLEMENTARY DATA
TOP
 ABSTRACT
 INTRODUCTION
 ADDITIONS TO SNP ANNOTATION
 DISCUSSION
 SUPPLEMENTARY DATA
REFERENCES
 
Supplementary Data are available at NAR Online.


   ACKNOWLEDGEMENTS
 
The VIB Switch laboratory was supported by a grant from the Federal Office for Scientific Affairs, Belgium (IUAP P6/43) and the Fund for Scientific Research (FWO Vlaanderen), Flanders. S.M.-S. was supported by a Marie Curie Intra-European fellowship. J.V.D. was supported by a grant from the Fund for Scientific Research, Flanders. L.C. was supported by a fellowship from CeGen (Genoma España), Spain. I.M. was supported by a fellowship from CIBER, Spain. Funding to pay the Open Access publication charges for this article was provided by a grant from the Federal Office of Scientific Affairs, Belgium (IUAP6/43).

Conflict of interest statement. None declared.


   Footnotes
 
The authors wish it to be known that, in their opinion, the first two authors should be regarded as joint First Authors.


   REFERENCES
TOP
 ABSTRACT
 INTRODUCTION
 ADDITIONS TO SNP ANNOTATION
 DISCUSSION
 SUPPLEMENTARY DATA
 REFERENCES
 

  1. Feuk L, Carson AR, Scherer SW. Structural variation in the human genome. Nat. Rev. Genet. (2006) 7:85–97.[Web of Science][Medline]

  2. Stranger BE, Forrest MS, Dunning M, Ingle CE, Beazley C, Thorne N, Redon R, Bird CP, de Grassi A, et al. Relative impact of nucleotide and copy number variation on gene expression phenotypes. Science (2007) 315:848–853.[Abstract/Free Full Text]

  3. Conde L, Vaquerizas JM, Dopazo H, Arbiza L, Reumers J, Rousseau F, Schymkowitz J, Dopazo J. PupaSuite: finding functional single nucleotide polymorphisms for large-scale genotyping purposes. Nucleic Acids Res. (2006) 34:W621–W625.[Abstract/Free Full Text]

  4. Reumers J, Maurer-Stroh S, Schymkowitz J, Rousseau F. SNPeffect v2.0: a new step in investigating the molecular phenotypic effects of human non-synonymous SNPs. Bioinformatics (2006) 22:2183–2185.[Abstract/Free Full Text]

  5. Reumers J, Schymkowitz J, Ferkinghoff-Borg J, Stricher F, Serrano L, Rousseau F. SNPeffect: a database mapping molecular phenotypic effects of human non-synonymous coding SNPs. Nucleic Acids Res. (2005) 33:D527–D532.[Abstract/Free Full Text]

  6. Conde L, Vaquerizas JM, Santoyo J, Al-Shahrour F, Ruiz-Llorente S, Robledo M, Dopazo J. PupaSNP Finder: a web tool for finding SNPs with putative effect at transcriptional level. Nucleic Acids Res. (2004) 32:W242–W248.[Abstract/Free Full Text]

  7. Conde L, Vaquerizas JM, Ferrer-Costa C, de la Cruz X, Orozco M, Dopazo J. PupasView: a visual tool for selecting suitable SNPs, with putative pathological effect in genes, for genotyping purposes. Nucleic Acids Res. (2005) 33:W501–W505.[Abstract/Free Full Text]

  8. Zhu X, Zhao X, Burkholder WF, Gragerov A, Ogata CM, Gottesman ME, Hendrickson WA. Structural analysis of substrate binding by the molecular chaperone DnaK. Science (1996) 272:1606–1614.[Abstract]

  9. Schymkowitz JW, Rousseau F, Martins IC, Ferkinghoff-Borg J, Stricher F, Serrano L. Prediction of water and metal binding sites and their affinities by using the Fold-X force field. Proc. Natl Acad. Sci. USA (2005) 102:10147–10152.[Abstract/Free Full Text]

  10. McCarty JS, Rudiger S, Schonfeld HJ, Schneider-Mergener J, Nakahigashi K, Yura T, Bukau B. Regulatory region C of the E. coli heat shock transcription factor, sigma32, constitutes a DnaK binding site and is conserved among eubacteria. J. Mol. Biol. (1996) 256:829–837.[CrossRef][Web of Science][Medline]

  11. The UniProt Consortium. The Universal Protein Resource (UniProt). Nucleic Acids Res (2007) 35:D193–D197.[Abstract/Free Full Text]

  12. Hubbard TJ, Aken BL, Beal K, Ballester B, Caccamo M, Chen Y, Clarke L, Coates G, Cunningham F, et al. Ensembl 2007. Nucleic Acids Res. (2007) 35:D610–D617.[Abstract/Free Full Text]

  13. Sherry ST, Ward MH, Kholodov M, Baker J, Phan L, Smigielski EM, Sirotkin K. dbSNP: the NCBI database of genetic variation. Nucleic Acids Res. (2001) 29:308–311.[Abstract/Free Full Text]

  14. Wang Z, Rolish ME, Yeo G, Tung V, Mawson M, Burge CB. Systematic identification and analysis of exonic splicing silencers. Cell (2004) 119:831–845.[CrossRef][Web of Science][Medline]

  15. Sandelin A, Alkema W, Engstrom P, Wasserman WW, Lenhard B. JASPAR: an open-access database for eukaryotic transcription factor binding profiles. Nucleic Acids Res. (2004) 32:D91–D94.[Abstract/Free Full Text]

  16. Guigo R. Assembling genes from predicted exons in linear time with dynamic programming. J. Comput. Biol. (1998) 5:681–702.[Web of Science][Medline]


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?


This article has been cited by other articles:


Home page
 Circ Cardiovasc GenetHome page
A. P. Reiner, M. D. Gross, C. S. Carlson, S. J. Bielinski, L. A. Lange, M. Fornage, N. S. Jenny, J. Walston, R. P. Tracy, O. D. Williams, et al.
Common Coding Variants of the HNF1A Gene Are Associated With Multiple Cardiovascular Risk Phenotypes in Community-Based Samples of Younger and Older European-American Adults: The Coronary Artery Risk Development in Young Adults Study and The Cardiovascular Health Study
Circ Cardiovasc Genet, June 1, 2009; 2(3): 244 - 254.
[Abstract] [Full Text] [PDF]

Home page
 Hum Mol GenetHome page
L. Quaye, D. Dafou, S. J. Ramus, H. Song, A. G. Maharaj, M. Notaridou, E. Hogdall, S. K. Kjaer, L. Christensen, C. Hogdall, et al.
Functional complementation studies identify candidate genes and common genetic variants associated with ovarian cancer survival
Hum. Mol. Genet., May 15, 2009; 18(10): 1869 - 1878.
[Abstract] [Full Text] [PDF]

This Article
Right arrow Abstract Freely available
Right arrow Print PDF (1560K) Freely available
Right arrow Screen PDF (267K) Freely available
Right arrow Supplementary Data
Right arrowOA All Versions of this Article:
36/suppl_1/D825    most recent
gkm979v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Commercial Re-use Guidelines
for Open Access NAR Content
Google Scholar
Right arrow Articles by Reumers, J.
Right arrow Articles by Schymkowitz, J.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Reumers, J.
Right arrow Articles by Schymkowitz, J.
Social Bookmarking
Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?