SNPWay: streamlined SNP-to-function and pathway over-representation analysis

Post-GWAS interpretation frequently requires translating variant lists (e.g., lead SNPs, clumped loci, credible sets, or curated panels) into pathway and functional hypotheses. In practice, obtaining pathway and functional over-representation results from SNP inputs often requires stitching together multiple tools for variant annotation, regulatory annotation, gene identifier handling, and statistical testing. This integration burden can reduce reproducibility and restrict end-to-end analysis to groups with dedicated bioinformatics support.

We present SNPWay, a web server and R package that performs end-to-end SNP-to-function and pathway over-representation analysis in a single standardized workflow. SNPWay accepts rsIDs, VCF files, or hg19/GRCh37 genomic coordinates. It queries Annotation Query (AnnoQ) to obtain SNP-to-gene mappings from ANNOVAR, SnpEff, and VEP under both Ensembl and RefSeq gene models, and incorporates enhancer-gene links via PEREGRINE to augment mappings for noncoding variants. SNPWay aggregates mapped genes into a single, non-redundant, combined gene list and submits it to PANTHER for over-representation testing against the Homo sapiens reference list, returning over-represented pathways and functional categories (e.g., Gene Ontology) with direct links for interactive exploration in PANTHER. SNPWay's modular architecture is designed for extensibility, enabling incorporation of additional analysis methods in future releases. A step-by-step walkthrough is provided in Supplementary Data.