Abstract:Currently, there are many publicly available Next Generation Sequencing tools developed for variant annotation and classification. However, as modern sequencing technology produces more and more sequencing data, a more efficient analysis program is desired, especially for variant analysis. In this study, we updated SNPAAMapper, a variant annotation pipeline by converting perl codes to python for generating annotation output with an improved computational efficiency and updated information for broader applicability. The new pipeline written in Python can classify variants by region (Coding Sequence, Untranslated Regions, upstream, downstream, intron), predict amino acid change type (missense, nonsense, etc.), and prioritize mutation effects (e.g., synonymous > non-synonymous) while being faster and more efficient. Our new pipeline works in five steps. First, exon annotation files are generated. Next, the exon annotation files are processed, and gene mapping and feature information files are produced. Afterward, the python scrips classify the variants based on genomic regions and predict the amino acid change category. Lastly, another python script prioritizes and ranks the mutation effects of variants to output the result file. The Python version of SNPAAMapper accomplished the overall speed by running most annotation steps in a substantially shorter time. The Python script can classify variants by region in 53 s compared to 166 s for the Perl script in a test sample run on a Latitude 7480 Desktop computer with 8GB RAM and an Intel Core i5-6300 CPU @ 2.4Ghz. Steps of predicting amino acid change type and prioritizing mutation effects of variants were executed within 1 s for both pipelines. SNPAAMapper-Python was developed and tested on the ClinVar database, a NCBI database of information on genomic variation and its relationship to human health. We believe our developed Python version of SNPAAMapper variant annotation pipeline will benefit the community by elucidating the variant consequence and speed up the discovery of causative genetic variants through whole genome/exome sequencing. Source codes, test data files, instructions, and further explanations are available on the web at https://github.com/BaiLab/SNPAAMapper-Python.

NAVIP: Unraveling the Influence of Neighboring Small Sequence Variants on Functional Impact Prediction

Exploring Functional Variant Discovery in Non-Coding Regions with Sinbad

A Comparison on Predicting Functional Impact of Genomic Variants.

DVA: Predicting the Functional Impact of Single Nucleotide Missense Variants

Enhancing missense variant pathogenicity prediction with protein language models using VariPred

Variant effect predictors: a systematic review and practical guide

ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data

VariPred: Enhancing Pathogenicity Prediction of Missense Variants Using Protein Language Models

QAFI: a novel method for quantitative estimation of missense variant impact using protein-specific predictors and ensemble learning

Understanding the characteristic of single nucleotide variants

SNPAAMapper-Python: A highly efficient genome-wide SNP variant analysis pipeline for Next-Generation Sequencing data

SNPeffect 5.0: large-scale structural phenotyping of protein coding variants extracted from next-generation sequencing data using AlphaFold models

Invmap: a Sensitive Mapping Tool for Long Noisy Reads with Inversion Structural Variants.

Prioritization of Nonsynonymous Single Nucleotide Variants for Exome Sequencing Studies Via Integrative Learning on Multiple Genomic Data

Prediction of pathogenic single amino acid substitutions using molecular fragment descriptors

Predicting functional impact of single amino acid polymorphisms by integrating sequence and structural features

TransEFVP: A Two-Stage Approach for the Prediction of Human Pathogenic Variants Based on Protein Sequence Embedding Fusion

Assessment of variant effect predictors unveils variants difficulty as a critical performance indicator

An Integrated Framework for Analysis and Prediction of Impact of Single Nucleotide Polymorphism Associated with Human Diseases

Multimodal learning of noncoding variant effects using genome sequence and chromatin structure

CADD-Splice—improving genome-wide variant effect prediction using deep learning-derived splice scores