SNP annotation

Source: Wikipedia, the free encyclopedia.
SNP annotation
ClassificationBioinformatics
SubclassificationSingle-nucleotide polymorphism
Type of tools usedFunctional annotation tools
Other subjects relatedGenome project, Genomics

Single nucleotide polymorphism annotation (

protein sequences.[1]

Introduction

Directed graph of relationships among SNP prediction webservers and their bioinformatics sources.[2]

Single nucleotide polymorphisms (SNPs) play an important role in

genetic mapping and population genetics.[3] Each SNP represents a nucleotide change between two individuals at a defined location. SNPs are the most common genetic variant found in all individual with one SNP every 100–300 bp in some species.[4] Since there is a massive number of SNPs on the genome
, there is a clear need to prioritize SNPs according to their potential effect in order to expedite genotyping and analysis. [5]

Annotating large numbers of SNPs is a difficult and complex process, which need computational methods to handle such a large dataset. Many tools available have been developed for SNP annotation in different organisms: some of them are optimized for use with

humans), but there are currently few tools available that are species non-specific or support non-model organism data. The majority of SNP annotation tools provide computationally predicted putative deleterious effects of SNPs. These tools examine whether a SNP resides in functional genomic regions such as exons, splice sites, or transcription regulatory sites, and predict the potential corresponding functional effects that the SNP may have using a variety of machine-learning approaches. But the tools and systems that prioritize functionally significant SNPs, suffer from few limitations: First, they examine the putative deleterious effects of SNPs with respect to a single biological function that provide only partial information about the functional significance of SNPs. Second, current systems classify SNPs into deleterious or neutral group.[6]

Many annotation algorithms focus on single nucleotide variants (SNVs), considered more rare than SNPs as defined by their minor allele frequency (MAF).[7][8] As a consequence, training data for the corresponding prediction methods may be different and hence one should be careful to select the appropriate tool for a specific purpose. For the purposes of this article, "SNP" will be used to mean both SNP and SNV, but readers should bear in mind the differences.

SNP annotation

Different type of annotations in genomics

For SNP annotation, many kinds of genetic and genomic information are used. Based on the different features used by each annotation tool, SNP annotation methods may be split roughly into the following categories:

Gene based annotation

Genomic information from surrounding genomic elements is among the most useful information for interpreting the biological function of an observed variant. Information from a known

protein sequence and its function. Gene based annotation is based on the fact that non-synonymous mutations can alter the protein sequence and that splice site mutation may disrupt the transcript splicing pattern.[9]

Knowledge based annotation

Knowledge base annotation is done based on the information of gene attribute, protein function and its

RNA secondary structure, RNA binding protein recognition and miRNA binding activity.[11][12]

Functional annotation

This method mainly identifies variant function based on the information whether the variant loci are in the known functional region that harbor genomic or epigenomic signals. The function of non-coding variants are extensive in terms of the affected genomic region and they involve in almost all processes of gene regulation from transcriptional to post translational level [13]

Transcriptional gene regulation

Transcriptional gene regulation process depends on many spatial and temporal factors in the nucleus such as global or local chromatin states, nucleosome positioning, TF binding, enhancer/promoter activities. Variant that alter the function of any of these biological processes may alter the gene regulation and cause phenotypic abnormality.[14] Genetic variants that located in distal regulatory region can affect the binding motif of TFs, chromatin regulators and other distal transcriptional factors, which disturb the interaction between enhancer/silencer and its target gene.[15]

Alternative splicing

Alternative splicing is one of the most important components that show functional complexity of genome. Modified splicing has significant effect on the phenotype that is relevance to disease or drug metabolism. A change in splicing can be caused by modifying any of the components of the splicing machinery such as splice sites or splice enhancers or silencers.[16] Modification in the alternative splicing site can lead to a different protein form which will show a different function. Humans use an estimated 100,000 different proteins or more, so some genes must be capable of coding for a lot more than just one protein. Alternative splicing occurs more frequently than was previously thought and can be hard to control; genes may produce tens of thousands of different transcripts, necessitating a new gene model for each alternative splice.

RNA processing and post transcriptional regulation

Mutations in the untranslated region (UTR) affect many post-transcriptional regulation. Distinctive structural features are required for many RNA molecules and cis-acting regulatory elements to execute effective functions during gene regulation. SNVs can alter the secondary structure of RNA molecules and then disrupt the proper folding of RNAs, such as tRNA/mRNA/lncRNA folding and miRNA binding recognition regions.[17]

Translation and post translational modifications

Single nucleotide variant can also affect the cis-acting regulatory elements in mRNA’s to inhibit/promote the translation initiation. Change in the synonymous codons region due to mutation may affect the translation efficiency because of codon usage biases. The translation elongation can also be retarded by mutations along the ramp of ribosomal movement. In the post-translational level, genetic variants can contribute to proteostasis and amino acid modifications. However, mechanisms of variant effect in this field are complicated and there are only a few tools available to predict variant’s effect on translation related modifications.[18]

Protein function

Non-synonymous is the variant in exons that change the amino acid sequence encoded by the gene, including single base changes and non frameshift indels. It has been extremely investigated the function of non-synonymous variants on protein and many algorithms have been developed to predict the deleteriousness and pathogenesis of single nucleotide variants (SNVs). Classical bioinformatics tools, such as SIFT, Polyphen and MutationTaster, successfully predict the functional consequence of non-synonymous substitution.[19][20][21][22] PopViz webserver provides a gene-centric approach to visualize the mutation damage prediction scores (CADD, SIFT, PolyPhen-2) or the population genetics (minor allele frequency) versus the amino acid positions of all coding variants of a certain human gene.[23] PopViz is also cross-linked with UniProt database, where the protein domain information can be found, and to then identify the predicted deleterious variants fall into these protein domains on the PopViz plot.[23]

Evolutionary conservation and nature selection

Comparative genomics approaches were used to predict the function-relevant variants under the assumption that the functional genetic locus should be conserved across different species at an extensive phylogenetic distance. On the other hand, some adaptive traits and the population differences are driven by positive selections of advantageous variants, and these genetic mutations are functionally relevant to population specific phenotypes. Functional prediction of variants’ effect in different biological processes is pivotal to pinpoint the molecular mechanism of diseases/traits and direct the experimental validation.[24]

List of available SNP annotation tools

To annotate the vast amounts of available NGS data, currently a large number of SNPs annotation tools are available. Some of them are specific to specific SNPs while others are more general. Some of the available SNPs annotation tools are as follows SNPeff, Ensembl Variant Effect Predictor (VEP), ANNOVAR, FATHMM, PhD-SNP, PolyPhen-2, SuSPect, F-SNP, AnnTools, SeattleSeq, SNPit, SCAN, Snap, SNPs&GO, LS-SNP, Snat, TREAT, TRAMS, Maviant, MutationTaster, SNPdat, Snpranker, NGS – SNP, SVA, VARIANT, SIFT, LIST-S2, PhD-SNP and FAST-SNP. The functions and approaches used in SNPs annotation tools are listed below.

Tools Description External resources use WebsiteURL References
PhyreRisk Maps genetics variants onto experimental and predicted protein structures Variant effect predictor, UniProt, Protein Data Bank, SIFTS, Phyre2 for predicted structures http://phyrerisk.bc.ic.ac.uk/home

[25]

Missense3D Reports structural impact of a missense variant onto PDB and user-supplied protein coordinates. Developed to be applicable to experimental and predicted protein structures Protein Data Bank, Phyre2 for predicted structures http://www.sbg.bio.ic.ac.uk/~missense3d/

[26]

SNPeff SnpEff annotates variants based on their genomic locations and predicts coding effects. Uses an interval forest approach ENSEMBL, UCSC and organism based e.g. FlyBase, WormBase and TAIR http://snpeff.sourceforge.net/SnpEff_manual.html [27]
Ensembl VEP Determines effects of variants (SNPs, insertions, deletions, CNVs or structural variants) on genes, transcripts, proteins and regulatory regions dbSNP, RefSeq, UniProt, COSMIC, PDBe, 1000 Genomes, gnomAD, PubMed https://www.ensembl.org/info/docs/tools/vep/index.html [28]
ANNOVAR This tool is suitable for pinpointing a small subset of functionally important variants. Uses mutation prediction approach for annotation UCSC, RefSeq and Ensembl http://annovar.openbioinformatics.org/ [29]
Jannovar This is a tool and library for genome annotation RefSeq, Ensembl, UCSC, etc. https://github.com/charite/jannovar [30]
PhD-SNP SVM-based method using sequence information retrieved by BLAST algorithm. UniRef90 http://snps.biofold.org/phd-snp/ [31]
PolyPhen-2 Suitable for predicting damaging effects of missense mutations. Uses sequence conservation, structure to model position of amino acid substitution, and SWISS-PROT annotation UniProt http://genetics.bwh.harvard.edu/pph2/ [32]
MutationTaster Suitable for predicting damaging effects of all intragenic mutations (DNA and protein level), including InDels. Ensembl, 1000 Genomes Project, ExAC, UniProt, ClinVar, phyloP, phastCons, nnsplice, polyadq (...) http://www.mutationtaster.org/ [33]
SuSPect An SVM-trained predictor of the damaging effects of missense mutations. Uses sequence conservation, structure and network (interactome) information to model phenotypic effect of amino acid substitution. Accepts VCF file UniProt, PDB, Phyre2 for predicted structures, DOMINE and STRING for interactome http://www.sbg.bio.ic.ac.uk/suspect/index.html [34]
F-SNP Computationally predicts functional SNPs for disease association studies. PolyPhen, SIFT, SNPeffect, SNPs3D, LS-SNP, ESEfinder, RescueESE, ESRSearch, PESX, Ensembl, TFSearch, Consite, GoldenPath, Ensembl, KinasePhos, OGPET, Sulfinator, GoldenPath http://compbio.cs.queensu.ca/F-SNP/ [35]
AnnTools Design to Identify novel and SNP/SNV, INDEL and SV/CNV. AnnTools searches for overlaps with regulatory elements, disease/trait associated loci, known segmental duplications and artifact prone regions dbSNP, UCSC, GATK refGene, GAD, published lists of common structural genomic variation, Database of Genomic Variants, lists of conserved TFBs, miRNA http://anntools.sourceforge.net/ [36]
SNPit Analyses the potential functional significance of SNPs derived from genome wide association studies dbSNP, EntrezGene, UCSC Browser, HGMD, ECR Browser, Haplotter, SIFT -/- [37]
SCAN Uses physical and functional based annotation to categorize according to their position relative to genes and according to linkage disequilibrium (LD) patterns and effects on expression levels -/- http://www.scandb.org/newinterface/about.html [38]
SNAP A neural network-based method for the prediction of the functional effects of non-synonymous SNPs Ensembl, UCSC, Uniprot, UniProt, Pfam, DAS-CBS, MINT, BIND, KEGG, TreeFam http://www.rostlab.org/services/SNAP [39]
SNPs&GO SVM-based method using sequence information, Gene Ontology annotation and when available protein structure. UniRef90, GO, PANTHER, PDB http://snps.biofold.org/snps-and-go/ [40]
LS-SNP Maps nsSNPs onto protein sequences, functional pathways and comparative protein structure models UniProtKB, Genome Browser, dbSNP, PD http://www.salilab.org/LS-SNP [41]
TREAT TREAT is a tool for facile navigation and mining of the variants from both targeted resequencing and whole exome sequencing -/- http://ndc.mayo.edu/mayo/research/biostat/stand-alone-packages.cfm [42]
SNPdat Suitable for species non-specific or support non-model organism data. SNPdat does not require the creation of any local relational databases or pre-processing of any mandatory input files -/- https://code.google.com/p/snpdat/downloads/ [43]
NGS – SNP Annotate SNPs comparing the reference amino acid and the non-reference amino acid to each orthologue Ensembl, NCBI and UniProt http://stothard.afns.ualberta.ca/downloads/NGS-SNP/ [44]
SVA Predicted biological function to variants identified NCBI RefSeq, Ensembl, variation databases, UCSC, HGNC, GO, KEGG, HapMap, 1000 Genomes Project and DG http://www.svaproject.org/ [45]
VARIANT VARIANT increases the information scope outside the coding regions by including all the available information on regulation, DNA structure, conservation, evolutionary pressures, etc. Regulatory variants constitute a recognized, but still unexplored, cause of pathologies dbSNP,1000 genomes, disease-related variants from GWAS, OMIM, COSMIC http://variant.bioinfo.cipf.es/ [46]
SIFT SIFT is a program that predicts whether an amino acid substitution affects protein function. SIFT uses sequence homology to predict whether an amino acid substitution will affect protein function PROT/TrEMBL, or NCBI's http://blocks.fhcrc.org/sift/SIFT.html [47]
LIST-S2 LIST-S2 (Local Identity and Shared Taxa, Species-specific) is based on the assumption that variations observed in closely related species are more significant when assessing conservation compared to those in distantly related species UniProt SwissProt/TrEMBL and NCBI Taxonomy https://gsponerlab.msl.ubc.ca/software/list/ [48][49]
FAST-SNP A web server that allows users to efficiently identify and prioritize high-risk SNPs according to their phenotypic risks and putative functional effects NCBI dbSNP, Ensembl, TFSearch, PolyPhen, ESEfinder, RescueESE, FAS-ESS, SwissProt, UCSC Golden Path, NCBI Blast and HapMap http://fastsnp.ibms.sinica.edu.tw/ [50]
PANTHER PANTHER relate protein sequence evolution to the evolution of specific protein functions and biological roles. The source of protein sequences used to build the protein family trees and used a computer-assisted manual curation step to better define the protein family clusters STKE, KEGG, MetaCyc, FREX and Reactome http://www.pantherdb.org/ [51]
Meta-SNP SVM-based meta predictor including 4 different methods. PhD-SNP, PANTHER, SIFT, SNAP http://snps.biofold.org/meta-snp [52]
PopViz Integrative and interactive gene-centric visualization of population genetics and mutation damage prediction scores of human gene variants gnomAD, Ensembl, UniProt, OMIM, UCSC, CADD, EIGEN, LINSIGHT, SIFT, PolyPhen-2, http://shiva.rockefeller.edu/PopViz/ [23]

Algorithms used in annotation tools

Variant annotation tools use machine learning algorithms to predict variant annotations. Different annotation tools use different algorithms. Common algorithms include:

  • Interval/Random forest-eg.MutPred, SNPeff
  • Neural networks
    -eg.SNAP
  • Support Vector Machines
    -e.g. PhD-SNP, SNPs&GO
  • Bayesian classification-eg.PolyPhen-2

Comparison of variant annotation tools

A large number of variant annotation tools are available for variant annotation. The annotation by different tools does not alway agree amongst each other, as the defined rules for data handling differ between applications. It is frankly impossible to perform a perfect comparison of the available tools. Not all tools have the same input and output nor the same functionality. Below is a table of major annotation tools and their functional area.

Tools Input file Output file SNP INDEL CNV WEB or Program Source
AnnoVar VCF, pileup,

CompleteGenomics, GFF3-SOLiD, SOAPsnp, MAQ, CASAVA

TXT Yes Yes Yes Program [53]
Jannovar VCF VCF Yes Yes Yes Java Program [54]
SNPeff VCF, pileup/TXT VCF, TXT, HTML Yes Yes No Program [27]
Ensembl VEP Ensembl default (coordinates), VCF, variant identifiers, HGVS, SPDI, REST-style regions VCF, VEP, TXT, JSON Yes Yes Yes Web, Perl script, REST API [55]
AnnTools VCF, pileup, TXT VCF Yes Yes No No [56]
SeattleSeq VVCF, MAQ, CASAVA,

GATK BED

VCF, SeattleSeq Yes Yes No Web [57]
VARIANT VCF, GFF2, BED web report, TXT Yes Yes Yes Web [58]

[59]

Application

Different annotations capture diverse aspects of variant function.

rare variants association analysis power of whole exome and whole genome sequencing studies.[61] Some tools have been developed to enable functionally-informed phenotype-genotype association analysis for common and rare variants by incorporating functional annotations in biobank-scale cohorts. [62][63][64][65]

Conclusions

The next generation of SNP annotation webservers can take advantage of the growing amount of data in core bioinformatics resources and use intelligent agents to fetch data from different sources as needed. From a user’s point of view, it is more efficient to submit a set of SNPs and receive results in a single step, which makes meta-servers the most attractive choice. However, if SNP annotation tools deliver heterogeneous data covering sequence, structure, regulation, pathways, etc., they must also provide frameworks for integrating data into a decision algorithms, and quantitative confidence measures so users can assess which data are relevant and which are not.

References

  1. doi:10.1016/S0981-9428(01)01242-6
    .
  2. PMID 19181721
    .
  3. PMID 19327864
    .
  4. ^ N. C. Oraguzie, E.H.A. Rikkerink, S.E. Gardiner, H.N. de Silva (eds.), "Association Mapping in Plants", Springer, 2007
  5. PMID 22247263
    .
  6. Computational Biology
    and Machine Learning Lab, School of Computing, Queen’s University, Kingston, ON, Canada
  7. ^ "Single-nucleotide polymorphism", Wikipedia, 2019-08-12, retrieved 2019-09-03
  8. ^ "Minor allele frequency", Wikipedia, 2019-08-12, retrieved 2019-09-03
  9. ^ M. J. Li, J. Wang, "Current trend of annotating single nucleotide variation in humans – A case study on SNVrap", Elsevier, 2014, pp. 1–9
  10. PMID 19015660
    .
  11. PMID 20808897
    .
  12. PMID 24476892
    .
  13. S2CID 8358824
    .
  14. PMID 24916300
    .
  15. PMID 23540573
    .
  16. PMID 21992029
    .
  17. PMID 23341781
    .
  18. ^ M. J. Li, J. Wang, "Current trend of annotating single nucleotide variation in humans – A case study on SNVrap", Elsevier, 2014, pp. 1–9
  19. ^ J. Wu, R. Jiang, "Prediction of Deleterious Nonsynonymous Single-Nucleotide Polymorphism for Human Diseases", The Scientific World Journal, 2013, 10 pages
  20. PMID 22689647
    .
  21. PMID 20354512
    .
  22. S2CID 26892938
    .
  23. ^
    PMID 30535305
    .
  24. ^ M. J. Li, J. Wang, "Current trend of annotating single nucleotide variation in humans – A case study on SNVrap", Elsevier, 2014, pp. 1–9
  25. PMID 31075275
    .
  26. PMID 30995449
    .
  27. ^
    PMID 22728672
    .
  28. PMID 27268795
    .
  29. PMID 20601685
    .
  30. S2CID 10822001
    .
  31. PMID 16895930
    .
  32. PMID 23315928
    .
  33. S2CID 26892938
    .
  34. PMID 24810707
    .
  35. PMID 17986460
    .
  36. PMID 22257670
    .
  37. PMID 19327864
    .
  38. PMID 19933162
    .
  39. PMID 17526529
    .
  40. S2CID 33900765
    .
  41. PMID 15827081
    .
  42. PMID 22088845
    .
  43. PMID 23390980
    .
  44. PMID 21697123
    .
  45. PMID 21624899
    .
  46. PMID 22693211
    .
  47. PMID 12824425
    .
  48. PMID 30952844
    .
  49. PMID 32352516
    .
  50. PMID 16845089
    .
  51. PMID 17130144
    .
  52. PMID 23819846
    .
  53. PMID 20601685
    .
  54. ^ "charite/jannovar". GitHub. Retrieved 2016-09-25.
  55. PMID 27268795
    .
  56. PMID 22257670
    .
  57. ^ "Input Variation List File for Annotation". SeattleSeq Annotation 151.
  58. PMID 22693211
    .
  59. PMID 23341494
    .
  60. PMID 29288389
    .
  61. PMID 32839606
    .
  62. PMID 29184056
    .
  63. S2CID 243873361
    .
  64. S2CID 253246835
    .
  65. S2CID 255084231
    .
Retrieved from "https://en.wikipedia.org/w/index.php?title=SNP_annotation&oldid=1188719322"