SNP@Promoter: a database of human SNPs (Single Nucleotide Polymorphisms) within the putative promoter regions

1471-2105-9-S1-S2 1471-2105 Proceedings SNP@Promoter: a database of human SNPs (Single Nucleotide Polymorphisms) within the putative promoter regions Kim Byoung-Chul bckim@kribb.re.kr Kim Woo-Yeon kimplove@kribb.re.kr Park Daeui daeui@kribb.re.kr Chung Won-Hyong whchung@kribb.re.kr Shin Kwang-sik shinks@kribb.re.kr Bhak Jong jong@kribb.re.kr

Korean BioInformation Center (KOBIC), KRIBB, Daejeon 305-806, Korea

BMC Bioinformatics Asia Pacific Bioinformatics Network (APBioNet) Sixth International Conference on Bioinformatics (InCoB2007) Shoba Ranganathan, Michael Gribskov and Tin Wee Tan Proceedings Sixth International Conference on Bioinformatics (InCoB2007) Hong Kong

27–30 August 2007

http://incob.apbionet.org/ 1471-2105 2008 9 Suppl 1 S2 http://www.biomedcentral.com/1471-2105/9/S1/S2 1831585110.1186/1471-2105-9-S1-S2 13 2 2008 2008 Kim et al; licensee BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Analysis of single nucleotide polymorphism (SNP) is becoming a key research in genomics fields. Many functional analyses of SNPs have been carried out for coding regions and splicing sites that can alter proteins and mRNA splicing. However, SNPs in non-coding regulatory regions can also influence important biological regulation. Presently, there are few databases for SNPs in non-coding regulatory regions.

Description

We identified 488,452 human SNPs in the putative promoter regions that extended from the +5000 bp to -500 bp region of the transcription start sites. Some SNPs occurring in transcription factor (TF) binding sites were also predicted (47,832 SNP; 9.8%). The result is stored in a database: SNP@promoter. Users can search the SNP@Promoter database using three entries: 1) by SNP identifier (rs number from dbSNP), 2) by gene (gene name, gene symbol, refSeq ID), and 3) by disease term. The SNP@Promoter database provides extensive genetic information and graphical views of queried terms.

Conclusion

We present the SNP@Promoter database. It was created in order to predict functional SNPs in putative promoter regions and predicted transcription factor binding sites. SNP@Promoter will help researchers to identify functional SNPs in non-coding regions.

Background

After finishing the Human Genome Project, biologists' interest has shifted to non-repetitive sequence variants in genome, by far the most common of which are single nucleotide polymorphisms (SNPs). For a variation to be considered an SNP, it must occur in at least 1% of the population. SNPs, which make up about 90% of all human genetic variation, occur every 100 to 300 bases along the 3-billion-base human genome 12. It is generally believed that the complete human sequence will reveal at least a million SNPs of coding regions, including introns and promoters. As a general rule, many SNPs have no effect on cell function, but some SNPs are reported to be highly related to diseases or to influence cells' response to a drug. Although more than 99% of human DNA sequences are the same across all populations, some SNPs can have a major impact on how humans respond to diseases; environmental insults such as bacteria, viruses, toxins, and chemicals; and drugs and other therapies. This makes SNPs of great value for biomedical research and for developing pharmaceutical products and for medical diagnostics.

New bioinformatics tools and public SNP resources for SNP studies, specifically for linkage disequilibrium and disease association studies, will form part of the new scientific landscape 3456789. These public SNP resources are possible through the large-scale and high-throughput systems to screen SNPs on many individuals. The challenge is to accomplish this while reducing the cost per genotype and required completion time. The public SNP resources are producing information about SNPs which are related to diseases or that modify biological function. Many functional studies of SNPs were focused on SNPs located in coding regions that can influence phenotype by altering the encoded proteins 910. They can also influence premature termination that can cause nonsense-mediated mRNA decay (NMD) 11. Another function of SNPs is that they affect splice sites which results in alternative splicing 12.

Additionally, there are many SNPs in non-coding regulatory regions. The exact functions of the non-coding regulatory region SNPs are not clear yet. However, some SNPs are predicted to be related to genes by influencing the binding affinity of transcription factors. For example, the G/C polymorphism in the promoter region of the FCGR2B promoter regulates gene expression 13. -783A/G and -1438A/G polymorphisms in the promoter of HTR2A gene regulate gene expression. -783 G allele and -1438 G allele are known to reduce the binding activity of transcription factors 14. However, there are no public resources that provide promoter information of SNPs influencing the non-coding regulatory regions in the human genome. The rSNP_Guide system is the only one that has reported SNPs that are related to potential transcription factor candidates among 41 types of known transcription factor binding sites. 1516. ORegAnno is focused not on SNP information of the regulatory regions in the human genome but on the registration and validation of SNPs from promoters, transcription factor binding sites, and regulatory variation 17.

SNP@Promoter is a large database that contains various types of information on the location and function for putative promoter regions in the human genome for gene regulation study. In particular, SNP@Promoter provides a platform for biologists including disease associated genes, transcription factor binding sites, and a graphic viewer.

Methods and results

We developed an integrated computational system for identifying SNPs in non-coding regulation regions (Fig 1). In this system, we: 1) predicted TF binding sites in putative promoter regions, 2) identified SNPs in the putative promoter regions and selected SNPs within predicted TF binding sites, 3) examined evolutionary conservation of predicted TF binding sites, and 4) integrated a variety of gene annotation information.

Figure 1

Flow chart for identifying SNPs in putative promoter regions

Flow chart for identifying SNPs in putative promoter regions. Cylinders represent databases. Rectangles are computational applications. (a) Putative promoter regions are identified in the human genome sequence. (b) Transcription binding sites are predicted in the putative promoter regions by using TransFac database. (c) SNPs are mapped. (d) Evolution conservation scores are calculated within transcription factor binding sites. (e) The disease association and functional annotation of target genes carried out by using an in-house functional annotation database.

Prediction of TF binding sites in putative promoter region

We identified TF binding sites in the putative promoter regions in the human genome. The promoter region is defined as the sequence of 5 kb upstream to 500 downstream bases of a transcription start site. The annotation information of genes, which is mapped to the genome, was obtained from the NCBI Gene database. To find TF binding sites in the putative promoter regions, we used the MATCH (Matrix Search For Transcription Factor Binding Site) program from the Transfac database (ver. 8.4) 1819. As a result, we predicted 1,497,317 TF binding sites from 28,644 human genes.

Identification of SNPs on predicted TF binding sites

The SNP annotation information was derived from a public SNP database (dbSNP ver. 126). We identified SNPs in putative promoter regions and selected SNPs that are predicted to be within TF binding sites. As a result, we mapped 488,452 SNPs and filtered out 47,832 SNPs within the putative TF binding sites.

Applying a conservation score

Using computational methods for predicting TFBS (TF binding sites) is not optimal due to a high false positive rate. However, recent algorithms have been improved in their reliability in TFBS prediction. Popular algorithms examine well-conserved regulatory sequences by comparing upstream sequences of orthologous genes across species 202122232425262728. Therefore, as an index of reliability for such an approach, we calculated an evolutionary conservation score for all the predicted TF binding sites. Users can see how reliable their predicted TF binding sites are. We used the phastcons16way file derived from UCSC human genome data. This file contains a conservation score from multiple genome alignment data calculated by the phastCons program 29.

Integration with functional annotation

The SNP@Promoter database adopted various gene annotations including pathways (KEGG), gene ontology (GOA), and disease information such as GAD, HGMD, and OMIM. The raw data files were integrated into the SNP@Promoter database based on a gene synonym table from HGNC (HUGO). These annotations provide insight into the effects of SNPs within TF binding sites and help users to characterize target genes regulated by SNPs.

User interface

As shown in Fig. 2(A), a user can search the SNP@Promoter database using three kinds of entries: 1) an SNP identifier (rs number from dbSNP), 2) a gene (Gene name, gene symbol, refSeq ID), or (3) a disease term. When the user submits a gene or a disease term, SNP@Promoter returns a gene list related to queries. In the case of accessing details of the query gene, it shows SNP information, gene information, and transcription factor binding site information of target genes as shown Fig. 2(B). SNP@Promoter provides graphical views of the queried SNPs and genes. Fig. 3 shows a putative promoter region browser.

Figure 2

SNP@Promoter user interface

SNP@Promoter user interface. SNP@Promoter main page. (A) Users can search using three entries: 1) an SNP identifier (rs number from dbSNP), 2) a gene (Gene name, gene symbol, refSeq ID), or 3) a disease term. (B) SNP@Promoter gene retrieval page. The SNP Information table shows identified SNPs within putative promoter region and TF biding sites. The Gene Information table shows various gene annotations including pathways (KEGG), gene ontology (GOA). The Information of Transcription Factor Binding Sites table shows a variety off TF information such as TF start position, upstream position, TF strand, match score, TF binding sequences, conservations score.

Figure 3

A graphic viewer of transcription regulatory region

A graphic viewer of transcription regulatory region. The green bar represents a putative promoter region (5500 bp). The arrows in the green bar show a strand of transcription, orange box is transcription start region, yellow inverted triangles are SNP positions, and purple triangles are predicted transcription binding sites.

Conclusion

SNP@Promoter is a database for functional SNPs within putative promoter regions and predicted TF binding sites. The database provides genetic information and graphical views of queried terms. SNP@Promoter will help researchers to identify functional SNPs in non-coding regions. Users can access the SNP@Promoter at http://variome.net or directly at http://variome.kobic.re.kr/SNPatPromoter.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

BK constructed the database. WYK developed the website and assisted to construction of database. WH and KS helped to develop the website. BK initiated this project and wrote the manuscript. DP assisted the manuscript writing. JB directed the study and helped to draft the manuscript.

Acknowledgements

We thank our colleagues at KOBIC, especially Areum Han. This project was supported by a grant from the KRIBB Research Initiative Program of Korea, by the Korea Science and Engineering Foundation (KOSEF) grant funded by the Korea government (MOST) (No. M10508040002-07N0804-00210), and by the MIC (Ministry of Information and Communication), Korea, under the KADO (Korea Agency Digital Opportunity and Promotion) support program (07-83).

This article has been published as part of BMC Bioinformatics Volume 9 Supplement 1, 2008: Asia Pacific Bioinformatics Network (APBioNet) Sixth International Conference on Bioinformatics (InCoB2007). The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2105/9?issue=S1.

A DNA polymorphism discovery resource for research on human genetic variation Collins FS Brooks LD Chakravarti A Genome Res 1998 8 1229 1231 9872978 The essence of SNPs Brookes AJ Gene 1999 234 177 186 10.1016/S0378-1119(99)00219-X 10395891 FESD: a Functional Element SNPs Database in human Kang HJ Choi KO Kim BD Kim S Kim YJ Nucleic Acids Res 2005 33 D518 522 10.1093/nar/gki082 15608252 540036 PicSNP: a browsable catalog of nonsynonymous single nucleotide polymorphisms in the human genome Chang H Fujita T Biochem Biophys Res Commun 2001 287 288 291 10.1006/bbrc.2001.5576 11549289 SNPper: retrieval and analysis of human SNPs Riva A Kohane IS Bioinformatics 2002 18 1681 1685 10.1093/bioinformatics/18.12.1681 12490454 SNPs3D: candidate gene and SNP selection for association studies Yue P Melamud E Moult J BMC Bioinformatics 2006 7 166 10.1186/1471-2105-7-166 16551372 1435944 SNPeffect: a database mapping molecular phenotypic effects of human non-synonymous coding SNPs Reumers J Schymkowitz J Ferkinghoff-Borg J Stricher F Serrano L Rousseau F Nucleic Acids Res 2005 33 D527 532 10.1093/nar/gki086 15608254 540040 MutDB services: interactive structural analysis of mutation data Dantzer J Moad C Heiland R Mooney S Nucleic Acids Res 2005 33 W311 314 10.1093/nar/gki404 15980479 1160165 SNP@Domain: a web resource of single nucleotide polymorphisms (SNPs) within protein domain structures and sequences Han A Kang HJ Cho Y Lee S Kim YJ Gong S Nucleic Acids Res 2006 1 34 W642 W644 10.1093/nar/gkl323 SIFT: predicting amino acid changes that affect protein function Ng PC Henikoff S Nucleic Acids Res 2003 31 3812 3814 10.1093/nar/gkg509 12824425 168916 SNP2NMD: a database of human single nucleotide polymorphisms causing nonsense-mediated mRNA decay Han A Kim WY Park SM Bioinformatics 2007 23 397 399 10.1093/bioinformatics/btl593 17121775 SNPSplicer: systematic analysis of SNP-dependent splicing in genotyped cDNAs ElSharawy A Manaster C Teuber M Rosenstiel P Kwiatkowski R Huse K Platzer M Becker A Nurnberg P Schreiber S Hampe J Hum Mutat 2006 27 1129 1134 10.1002/humu.20377 16937379 Decreased transcription of the human FCGR2B gene mediated by the -343 G/C promoter polymorphism and association with systemic lupus erythematosus Blank MC Stefanescu RN Masuda E Marti F King PD Redecha PB Wurzburger RJ Peterson MG Tanaka S Pricop L Hum Genet 2005 117 220 227 10.1007/s00439-005-1302-3 15895258 Polymorphisms in the regulatory region of the human serotonin 5-HT2A receptor gene (HTR2A) influence gene expression Myers RL Airey DC Manier DH Shelton RC Sanders-Bush E Biol Psychiatry 2007 61 167 173 10.1016/j.biopsych.2005.12.018 16697352 rSNP_Guide: An integrated database-tools system for studying SNPs and site-directed mutations in transcription factor binding sites Ponomarenko JV Orlova GV Merkulova TI Gorshkova EV Fokin ON Vasiliev GV Frolov AS Ponomarenko MP Hum Mutat 2002 20 239 248 10.1002/humu.10116 12325018 rSNP_Guide, a database system for analysis of transcription factor binding to DNA with variations: application to genome annotation Ponomarenko JV Merkulova TI Orlova GV Fokin ON Gorshkova EV Frolov AS Valuev VP Ponomarenko MP Nucleic Acids Res 2003 31 118 121 10.1093/nar/gkg112 12519962 165559 ORegAnno: an open access database and curation system for literature-derived promoters, transcription factor binding sites and regulatory variation Montgomery SB Griffith OL Sleumer MC Bergman CM Bilenky M Pleasance ED Prychyna Y Zhang X Jones SJ Bioinformatics 2006 22 637 640 10.1093/bioinformatics/btk027 16397004 TRANSFAC: transcriptional regulation, from patterns to profiles Matys V Fricke E Geffers R Goßling E Haubrock M Hehl R Hornischer K Karas D Kel AE Kel-Margoulis OV Kloos DU Land S Lewicki-Potapov B Michael H Munch R Reuter I Rotert S Saxel H Scheer M Thiele S Wingender E Nucleic Acids Res 2003 31 374 378 10.1093/nar/gkg108 12520026 165555 MATCH: A tool for searching transcription factor binding sites in DNA sequences Kel AE Gossling E Reuter I Cheremushkin E Kel-Margoulis OV Wingender E Nucleic Acids Res 2003 31 3576 3579 10.1093/nar/gkg585 12824369 169193 Evolutionary strategies for the elucidation of cis and trans factors that regulate the developmental switching programs of the beta-like globin genes Gumucio DL Shelton DA Zhu W Millinoff D Gray T Bock JH Slightom JL Goodman M Mol Phylogenet Evol 1996 5 18 32 10.1006/mpev.1996.0004 8673285 Long human-mouse sequence alignments reveal novel regulatory elements: A reason to sequence the mouse genome Hardison RC Oeltjen J Miller W Genome Res 1997 7 959 966 9331366 Conserved noncoding sequences are reliable guides to regulatory elements Hardison RC Trends Genet 2000 16 369 372 10.1016/S0168-9525(00)02081-3 10973062 Enrichment of regulatory signals in conserved non-coding genomic sequence Levy S Hannenhalli S Workman C Bioinformatics 2001 17 871 877 10.1093/bioinformatics/17.10.871 11673231 rVista for Comparative Sequence-Based Discovery of Functional Transcription Factor Binding Sites Loots GG Ovcharenko I Pachter L Dubchak I Rubin EM Genome Res 2002 12 832 839 10.1101/gr.225502. Article published online before print in April 2002 11997350 186580 AthaMap web tools for database-assisted identification of combinatorial cis-regulatory elements and the display of highly conserved transcription factor binding sites in Arabidopsis thaliana Steffens NO Galuschka C Schindler M Bulow L Hehl R Nucleic Acids Res 2005 33 W397 W402 10.1093/nar/gki395 15980498 1160156 Cross-species comparison significantly improves genome-wide prediction of cis-regulatory modules in Drosophila Sinha S Schroeder MD Unnerstall U Gaul U Siggia ED BMC Bioinformatics 2004 5 129 10.1186/1471-2105-5-129 15357878 521067 Computational Prediction of cis-Regulatory Modules from Multispecies Alignments Using Galaxy, Table Browser, and GALA Elnitski L King D Hardison RC Methods Mol Biol 2006 338 91 103 16888352 Identifying cis-regulatory modules by combining comparative and compositional analysis of DNA Pierstorff N Bergman CM Wiehe T Bioinformatics 2006 22 2858 64 10.1093/bioinformatics/btl499 17032682 Evolutionarily conserved elements in vertebrate, insect, worm, and yeast genomes Siepel A Bejerano G Pedersen JS Hinrichs A Hou M Rosenbloom K Clawson H Spieth J Hillier LW Richards S Weinstock GM Wilson RK Gibbs RA Kent WJ Miller W Haussler D Genome Res 2005 15 1034 1050 10.1101/gr.3715005 16024819 1182216