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.

A generic pipeline for CADD score generation: chickenCADD and turkeyCADD

CADD v1.7: using protein language models, regulatory CNNs and other nucleotide-level scores to improve genome-wide variant predictions

A general framework for estimating the relative pathogenicity of human genetic variants

varCADD: large sets of standing genetic variation enable genome-wide pathogenicity prediction

SNPpath: characterizing cattle SNPs by enriched pathway terms.

DANN: a deep learning approach for annotating the pathogenicity of genetic variants

PhD-SNPg: a webserver and lightweight tool for scoring single nucleotide variants

The GenoPred Pipeline: A Comprehensive and Scalable Pipeline for Polygenic Scoring.

Rapid discrimination between deleterious and benign missense mutations in the CAGI 6 experiment

An integrated pipeline for prediction of Clostridioides difficile infection

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

Cactus: a user-friendly and reproducible ATAC-Seq and mRNA-Seq analysis pipeline for data preprocessing, differential analysis, and enrichment analysis

Computational Assessment of the Regulation-Modulating Potential for Noncoding Variants

Variant analysis pipeline for accurate detection of genomic variants from transcriptome sequencing data

CandiHap: a haplotype analysis toolkit for natural variation study

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

Evaluation of 10 Different Pipelines for Bacterial Single-Nucleotide Variant Detection

CUBAP: an interactive web portal for analyzing codon usage biases across populations

NCBI prokaryotic genome annotation pipeline

Computational Assessment of the Expression-modulating Potential for Non-coding Variants.

Whole Animal Genome Sequencing: user-friendly, rapid, containerized pipelines for processing, variant discovery, and annotation of short-read whole genome sequencing data