Easy-PSAP: an integrated workflow to prioritize pathogenic variants in sequence data from a single individual.

Abstract

Introduction Next-Generation Sequencing data analysis has become an integral part of clinical genetic diagnosis, raising the question of variant prioritization. The Population Sampling Probability (PSAP) method has been developed to tackle the issue of variant prioritization in the exome of a single patient, by leveraging allele frequencies from population databases and a variant pathogenicity score. Methods Here, we present Easy-PSAP, a completely new implementation of the PSAP method comprising two user-friendly and highly adaptable pipelines. Easy-PSAP allows the gene-based recalibration of any in silico pathogenicity prediction score compared to scores of variants seen in the general population, including popular scores like CADD or AlphaMissense. Easy-PSAP can evaluate genetic variants at the scale of a whole exome or a whole genome using information from the latest population and annotation databases. Results Through simulations on synthetic disease exomes, we show that Easy-PSAP is able to rank more than 50% of causal pathogenic variants in the top 10 variants for an autosomal dominant model of transmission and in the top 1 for an autosomal recessive model of transmission. Discussion These findings, along with the accessibility of the pipeline to both researchers and clinicians, make Easy-PSAP a state-of-the-art tool for variant prioritization in Next Generation Sequencing (NGS) data that can continue to evolve as new frameworks and databases become available. Easy-PSAP is implemented in R and bash within an open-source Snakemake framework. It is available on GitHub alongside conda environments containing the required dependencies (https://github.com/msogloblinsky/Easy-PSAP).