1. A framework for pathway knowledge driven prioritization in genome-wide association studies
- Author
-
Shrayashi Biswas, Partha P. Majumder, Samsiddhi Bhattacharjee, and Soumen Pal
- Subjects
Prioritization ,Epidemiology ,Computer science ,Gwas data ,Genome-wide association study ,Computational biology ,Overfitting ,Polymorphism, Single Nucleotide ,Article ,03 medical and health sciences ,Databases, Genetic ,Humans ,Disease process ,Computer Simulation ,Genetics (clinical) ,030304 developmental biology ,Genetic association ,0303 health sciences ,Models, Genetic ,Genome, Human ,030305 genetics & heredity ,Molecular Sequence Annotation ,ComputingMethodologies_PATTERNRECOGNITION ,Diabetes Mellitus, Type 2 ,Scalability ,Type I and type II errors ,Genome-Wide Association Study - Abstract
Many variants with low frequencies or with low to modest effects likely remain unidentified in genome-wide association studies (GWAS) because of stringent genome-wide thresholds for detection. To improve the power of detection, variant prioritization based on their functional annotations and epigenetic landmarks has been used successfully. Here, we propose a novel method of prioritization of a GWAS by exploiting gene-level knowledge (e.g., annotations to pathways and ontologies) and show that it further improves power. Often, disease associated variants are found near genes that are coinvolved in specific biological pathways relevant to disease process. Utilization of this knowledge to conduct a prioritized scan increases the power to detect loci that map to genes clustered in a few specific pathways. We have developed a computationally scalable framework based on penalized logistic regression (termed GKnowMTest-Genomic Knowledge-guided Multiplte Testing) to enable a prioritized pathway-guided GWAS scan with a very large number of gene-level annotations. We demonstrate that the proposed strategy improves overall power and maintains the Type 1 error globally. Our method works on genome-wide summary level data and a user-specified list of pathways (e.g., those extracted from large pathway databases without reference to biology of a specific disease). It automatically reweights the input p values by incorporating the pathway enrichments as "adaptively learned" from the data using a cross-validation technique to avoid overfitting. We used whole-genome simulations and some publicly available GWAS data sets to illustrate the application of our method. The GKnowMTest framework has been implemented as a user-friendly open-source R package.
- Published
- 2020