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.

Segpy: a streamlined, user-friendly pipeline for variant segregation analysis

shinyseg: a web application for flexible cosegregation and sensitivity analysis

vSNP: a SNP pipeline for the generation of transparent SNP matrices and phylogenetic trees from whole genome sequencing data sets

Pediatric Cancer Variant Pathogenicity Information Exchange (pecanpie): a Cloud-Based Platform for Curating and Classifying Germline Variants.

Scalable and efficient DNA sequencing analysis on different compute infrastructures aiding variant discovery

The use of segregation analysis in interpretation of sequence variants in SMAD3: A case report

snpQT: flexible, reproducible, and comprehensive quality control and imputation of genomic data

Seqspark: A Complete Analysis Tool for Large-Scale Rare Variant Association Studies Using Whole-Genome and Exome Sequence Data

Rare copy number variant analysis in case-control studies using snp array data: a scalable and automated data analysis pipeline

Onkopipe: A Snakemake Based DNA-Sequencing Pipeline for Clinical Variant Analysis in Precision Medicine

SEQCHIP: a Powerful Method to Integrate Sequence and Genotype Data for the Detection of Rare Variant Associations

Robust and Rapid Algorithms Facilitate Large-Scale Whole Genome Sequencing Downstream Analysis in an Integrative Framework

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

A high-performance computing toolset for relatedness and principal component analysis of SNP data

SPANDx: a genomics pipeline for comparative analysis of large haploid whole genome re-sequencing datasets

PEAS V1.0: a Package for Elementary Analysis of SNP Data

Genopyc: a Python library for investigating the functional effects of genomic variants associated to complex diseases

SEGUL: Ultrafast, memory‐efficient and mobile‐friendly software for manipulating and summarizing phylogenomic datasets

Genomevip: A Cloud Platform For Genomic Variant Discovery And Interpretation

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

PoolParty2: An integrated pipeline for analysing pooled or indexed low-coverage whole-genome sequencing data to discover the genetic basis of diversity