Skip to main content

Functional analysis of novel SNPs and mutations in human and mouse genomes

Abstract

Background

With the flood of information generated by the new generation of sequencing technologies, more efficient bioinformatics tools are needed for in-depth impact analysis of novel genomic variations. FANS (Functional Analysis of Novel SNPs) was developed to streamline comprehensive but tedious functional analysis steps into a few clicks and to offer a carefully designed presentation of results so researchers can focus more on thinking instead of typing and calculating.

Results

FANS http://fans.ngc.sinica.edu.tw/ harnesses the power of public information databases and powerful tools from six well established websites to enhance the efficiency of analysis of novel variations. FANS can process any point change in any coding region or GT-AG splice site to provide a clear picture of the disease risk of a prioritized variation by classifying splicing and functional alterations into one of nine risk subtypes with five risk levels.

Conclusion

FANS significantly simplifies the analysis operations to a four-step procedure while still covering all major areas of interest to researchers. FANS offers a convenient way to prioritize the variations and select the ones with most functional impact for validation. Additionally, the program offers a distinct improvement in efficiency over manual operations in our benchmark test.

Background

As sequencing technologies continue to advance rapidly and the cost of large-scale genotyping and sequencing falls markedly as a result, there is an increasing need for more efficient bioinformatics tools with which to analyze the flood of information regarding novel variations in human and mouse genomes. For instance, when novel mutations are found in a cancer genome re-sequencing project or when novel SNPs are pinpointed in a genetic study, a bioinformatics tool that can efficiently analyze and prioritize these novel variations for downstream research will surely save a lot of time and labor. Currently, though many well developed web projects like SNP@Promoter [1], F-SNP [2], and Snap [3] provide rich information for functional SNP prediction and integrated SNP annotation, they can only process known SNPs rather than newly found ones. In the analysis of novel variations, the Genomic Mutation Consequence Calculator (GMCC) [4] does offer some help to researchers by allowing human mutation analysis to the synonymous/non-synonymous level. However, to perform more in-depth impact analysis in human and mouse genomes, researchers still need to use a suite of separate online software applications like Rescue-ESE [5] for candidate exonic splicing enhancers, Fas-ESS [6] for candidate exonic splicing silencers, SIFT [7] for effects on protein function, and PolyPhen [8] for structure and function alterations of a human protein. These analyses of functional consequences caused by novel variations require many tedious and error-prone manual steps. To alleviate this problem, FANS (Functional Analysis of Novel SNPs) as an extended application of FASTSNP [9] streamlines information retrieval and functional analysis for multiple novel SNPs as well as mutations, with versatile input options and a graphical presentation of results, and provides a great efficiency improvement over manual operations.

Implementation

Design

According to their location and the kind of genetic alteration they cause, DNA variations are divided into five groups: non-synonymous, synonymous, splice site, intronic, and untranslated. Variations classified into the first three groups are those most likely to have significant functional impact. The first group may substantially affect protein function [10]: a non-synonymous variation leads to an amino acid change or makes a premature stop codon. Variations of the second [11] and the third categories [12] can affect pre-mRNA splicing. For better prioritization, FANS nominates variations as one of five risk types, with nine risk subtypes and five risk levels (Table 1) [13, 14]. Each analyzed variation is assigned one risk type and an accompanying risk level. The impact analysis offered by FANS includes examination of diminished ESE or ESS, altered protein function, protein domain abolition, and GT-AG splice site variation. A very important design feature of FANS is that the analysis is always based on the most up-to-date data retrieved from the source websites NCBI (National Center of Biotechnology Information) [15], Ensembl [16], UCSC BLAT [17], Rescue-ESE [5], Fas-ESS [6], and SIFT [7].

Table 1 Risk types and risk levels used for prioritization of variations

Input

Users of FANS can submit a single variation or a batch query. Currently, FANS offers users the option to search two species, human and mouse. After choosing the species, users need to input a sequence in FASTA format or provide information for each variation including chromosome number, physical position, orientation, and allele change, either through input fields on the web page or by uploading a batch file. For a batch file, each variation occupies a line that must contain the all required information for that variation. An example of a batch file is given online to illustrate the data format accepted by FANS. As all gene information and prediction results are acquired through the Internet in real-time, the maximum allowed number of variations for each batch query is currently 100 to avoid overloading the source websites.

Analysis flow

Figure 1 outlines the analysis flow that is started when a query is submitted. For a sequence query, the FASTA format sequence is sent to UCSC BLAT [17] search. FANS automatically generates the required information about the novel variations from all the returned BLAT results, and users only need to choose a group of variations in one of the BLAT results for the next analysis stage. For each variation, FANS uses the submitted information to search NCBI Map Viewer to retrieve all transcripts of the gene covering the variation. If the variation falls in a non-coding region, it will be checked for its GT-AG splice site risk. Otherwise, the analysis will translate the coding region and then follow either a 'synonymous' or a 'non-synonymous' flow. A synonymous coding variation which affects amino acid sequence is first subjected to ESE and ESS analysis. The exonic DNA fragment of eleven nucleotides centered around the variation is extracted and transmitted to Rescue-ESE [5] and Fas-ESS [6] for ESE and ESS hexamers pattern matching. Any variation found to be covered by an ESE or ESS motif that is likely to diminish exon splicing is subsequently sent to protein domain abolition analysis. At this stage, FANS utilizes protein domain information from NCBI GenPept [18] to check protein domain abolition that may be caused by splicing regulation when the diminished exon splicing changes amino acid sequence in a protein domain.

Figure 1
figure 1

FANS prioritizing and analysis flow path.

A non-synonymous variation that produces a stop codon, resulting in very serious protein structure alteration, is categorized as risk type "Non-sense" with risk level "Very High". Other than the "Non-sense" type, FANS takes advantage of SIFT [7] to see if a particular amino acid replacement may affect protein function. Those variations predicted to affect protein function are categorized as "Mis-sense (non-conservative change)". If no significant functional impact is found despite the substitution of an amino acid, then ESE, ESS, and protein domain abolition analysis will be carried out.

Output

By utilizing retrieved information and analyzed results from six websites, FANS efficiently prioritizes novel variations according to their risk levels in a few seconds with just a few clicks. The integrated results are divided into four parts for easy visualization: Genome View, Gene View, Transcript View, and Variation View (Figure 2).

Figure 2
figure 2

The FANS results page. The output of an analysis of eighteen variations in eight chromosomes.

The first screen, Genome View, shows an overview of the chromosome locations and the risk levels of the queried variations. The higher the risk level, the warmer the color that is used, with red color representing a very high risk level and blue a very low risk. Clicking on a variation label will move from Genome View to Gene View to display gene information of that variation. Gene View also displays the scale, location and all transcripts accommodating the selected variation. A gene selection list is provided for a user to choose the gene he/she is interested in. The transcript accession numbers from NCBI are printed on the left side of the transcript picture and transcript IDs from Ensembl are marked on the other side. Users can link to the source web pages by clicking each ID.

Next, in Transcript View, FANS offers separate transcript tabs for transcript selection (provided that a gene has more than one transcript). Gene structure comprising introns and exons is shown here. Extended exonic regions together with colored vertical lines are drawn for better illustration of variation locations. When the mouse is moved over a line, a pop-up window will show its risk level and risk type. In addition, an upwards-pointing arrow further distinguishes the position of the selected variation.

Finally, Variation View depicts the analysis flow path as well as the collected results of an analyzed variation. Selecting an interesting variation out of the variation list box will bring up its associated results and analysis flow path. FANS colors the flow path that an analysis has gone through according to the final risk level outcome. The highlighted path likewise provides corresponding description pop-up windows when pointed at with the mouse. Moreover, users can download not only all analyzed results in CSV format (by clicking the icon right of Transcript View) but also a 200 bps flanking sequence of any variation listed in Variation View.

Development environment

The following software was used to construct FANS:

  1. i.

    Red Hat Enterprise Linux 4 AS Update 6

  2. ii.

    Eclipse 3.2.2

  3. iii.

    Subversion 1.2.1

  4. iv.

    Tomcat 5.5.20

  5. v.

    RoboSuite 5.5 SP2

  6. vi.

    Java SE Development Kit 1.5.0_11

  7. vii.

    Struts Framework 1.3

  8. viii.

    BioJava 1.5

Java is the core language for integration and data calculation from different websites and for the generation of the final results. BioJava, an open-source project, is used for amino acid translation. RoboSuite handles the submission and extraction of data from websites. All development processes were done on Eclipse and Subversion was adopted as our revision control system.

Results

FANS automates analysis flow and simplifies operations so that researchers can spend more time on thinking instead of mouse clicking and dragging. Without FANS, the manual analysis of just a single variation requires a lot of tedious data extraction and transformation, which is compounded when there are many to analyze. Major processing steps include, for example, the retrieval of intronic and exonic information from every transcript of a gene covering the variation, the calculation of the position for translation if the variation occurs in a coding region, the extraction of an 11-mer sequence for ESE and ESS processing, the capture of protein domain information from NCBI GenPept flat file, the generation of an amino acid change for the SIFT program, and so on. In a manual operation benchmarking test, over an hour (3968 s) and 940 discrete steps were needed for a skilled person to operate all analytic processes over 10 variations in one gene (APOE). Two fifths of the operations were clicks for linking and submission, 170 steps were arithmetic, and others were copy/paste jobs. With FANS, the same analysis took only 118 s and 4 simple steps: selecting a query type, choosing a species, filling in variation information, and clicking the submit button, which dramatically reduced the time required and trimmed the returned results from forty-seven discrete pages to one single well-structured web page (Table 2).

Table 2 A benchmark comparison of analysis of 10 variations by hand and with FANS.

Conclusion

FANS is uniquely designed as a convenient way to perform in-depth integrated functional analysis of novel variations with always up-to-date data. With the application of next-generation sequencing platforms such as Roche GS FLX in genetic association studies or cancer research, many novel variations in patient samples will be unveiled. Under these circumstances, FANS offers a convenient way to prioritize the variations and select the ones with most functional impact for validation. In addition, more of the researchers' time and effort can be dedicated to productive thinking instead of error-prone manual operations. FANS significantly simplifies the analysis operations to a four-step procedure while still covering all major areas that researchers are interested in. Additionally, our benchmark tests showed that the operation time required for a sample analysis fell from more than one hour to less than two minutes. The sequence-based submission allows users to submit simply a human or a mouse genomic sequence containing targeted variations. The top-down design of views and easy data download plus many other intuitive features should make FANS quite user-friendly to researchers.

Our future development of this program will include database localization, an essential step to avoid overloading the websites, and one which should dramatically shorten processing times. The incorporation of more analytic tools such as PolyPhen [8] into our analysis flow is also an important goal in the further development of FANS.

Availability and requirements

Project name: FANS

Project home page: http://fans.ngc.sinica.edu.tw/

Operation system(s): Platform independent

Programming language: Java

Other requirements: Standard web browser with JavaScript enabled

License: The tool is freely available to academic and non-academic users.

Any restrictions to use by non-academics: None

Abbreviations

SNP:

Single Nucleotide Polymorphism

ESE:

Exonic Splicing Enhancer

ESS:

Exonic Splicing Silencer

UCSC:

University of California, Santa Cruz

CSV:

Comma Separated Values.

References

  1. Kim BC, Kim WY, Park D, Chung WH, Shin KS, Bhak J: SNP@Promoter: a database of human SNPs (single nucleotide polymorphisms) within the putative promoter regions. BMC Bioinformatics 2008,9(Suppl 1):S2. 10.1186/1471-2105-9-S1-S2

    Article  PubMed Central  PubMed  Google Scholar 

  2. Lee PH, Shatkay H: F-SNP: computationally predicted functional SNPs for disease association studies. Nucleic Acids Res 2008, 36: D820–824. 10.1093/nar/gkm904

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  3. Li S, Ma L, Li H, Vang S, Hu Y, Bolund L, Wang J: Snap: an integrated SNP annotation platform. Nucleic Acids Res 2007, 35: D707–710. 10.1093/nar/gkl969

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  4. Major JE: Genomic mutation consequence calculator. Bioinformatics 2007, 23: 3091–3092. 10.1093/bioinformatics/btm339

    Article  CAS  PubMed  Google Scholar 

  5. Fairbrother WG, Yeo GW, Yeh R, Goldstein P, Mawson M, Sharp PA, Burge CB: RESCUE-ESE identifies candidate exonic splicing enhancers in vertebrate exons. Nucleic Acids Res 2004, 32: W187–190. 10.1093/nar/gkh393

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  6. 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. 10.1016/j.cell.2004.11.010

    Article  CAS  PubMed  Google Scholar 

  7. Ng PC, Henikoff S: SIFT: predicting amino acid changes that affect protein function. Nucleic Acids Res 2003, 31: 3812–3814. 10.1093/nar/gkg509

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Ramensky V, Bork P, Sunyaev S: Human non-synonymous SNPs: server and survey. Nucleic Acids Res 2002, 30: 3894–3900. 10.1093/nar/gkf493

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Yuan HY, Chiou JJ, Tseng WH, Liu CH, Liu CK, Lin YJ, Wang HH, Yao A, Chen YT, Hsu CN: FASTSNP: an always up-to-date and extendable service for SNP function analysis and prioritization. Nucleic Acids Res 2006, 34: W635–641. 10.1093/nar/gkl236

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  10. DePristo MA, Weinreich DM, Hartl DL: Missense meanderings in sequence space: a biophysical view of protein evolution. Nat Rev Genet 2005, 6: 678–687. 10.1038/nrg1672

    Article  CAS  PubMed  Google Scholar 

  11. Cartegni L, Chew SL, Krainer AR: Listening to silence and understanding nonsense: exonic mutations that affect splicing. Nat Rev Genet 2002, 3: 285–298. 10.1038/nrg775

    Article  CAS  PubMed  Google Scholar 

  12. Shin C, Manley JL: Cell signalling and the control of pre-mRNA splicing. Nat Rev Mol Cell Biol 2004, 5: 727–738. 10.1038/nrm1467

    Article  CAS  PubMed  Google Scholar 

  13. Tabor HK, Risch NJ, Myers RM: Candidate-gene approaches for studying complex genetic traits: practical considerations. Nat Rev Genet 2002, 3: 391–397. 10.1038/nrg796

    Article  CAS  PubMed  Google Scholar 

  14. Yeo G, Hoon S, Venkatesh B, Burge CB: Variation in sequence and organization of splicing regulatory elements in vertebrate genes. Proc Natl Acad Sci USA 2004, 101: 15700–15705. 10.1073/pnas.0404901101

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Wheeler DL, Barrett T, Benson DA, Bryant SH, Canese K, Chetvernin V, Church DM, DiCuccio M, Edgar R, Federhen S, et al.: Database resources of the National Center for Biotechnology Information. Nucleic Acids Res 2007, 35: D5–12. 10.1093/nar/gkl1031

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Curwen V, Eyras E, Andrews TD, Clarke L, Mongin E, Searle SM, Clamp M: The Ensembl automatic gene annotation system. Genome Res 2004, 14: 942–950. 10.1101/gr.1858004

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. Kent WJ: BLAT – the BLAST-like alignment tool. Genome Res 2002, 12: 656–664.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  18. Pruitt KD, Tatusova T, Maglott DR: NCBI Reference Sequence (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res 2005, 33: D501–504. 10.1093/nar/gki025

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The authors thank Dr. Chih-Cheng Chen (IBMS) for his contribution towards the initiation of the project. Special thanks to Dr. Jer-Yuarn Wu (NGC) who provided skillful programmers and to Dr. Harry Wilson for manuscript editing. This work was supported by grants from National Science Council, Taiwan, under NRPGM program and Academia Sinica.

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

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Adam Yao.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

CKL designed the analysis flow and wrote the manuscript. CKL, YHC, CYT, and SCC implemented the program. CKL, YHC, YJL, MFT, and AY had substantial contribution to the user interface design. AY and YTC supervised the FANS project.

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Liu, CK., Chen, YH., Tang, CY. et al. Functional analysis of novel SNPs and mutations in human and mouse genomes. BMC Bioinformatics 9 (Suppl 12), S10 (2008). https://doi.org/10.1186/1471-2105-9-S12-S10

Download citation

  • Published:

  • DOI: https://doi.org/10.1186/1471-2105-9-S12-S10

Keywords