Abstract:A key challenge in the application of whole-genome sequencing (WGS) for clinical diagnostic and research is the high-throughput prioritization of functional variants in the non-coding genome. This challenge is compounded by context-specific genetic modulation of gene expression, and variant-gene mapping depends on the tissues and organ systems affected in a given disease; for instance, a disease affecting the gastrointestinal system would use maps specific to genome regulation in gut-related tissues. While there are large-scale atlases of genome regulation, such as GTEx and NIH Roadmap Epigenomics, the clinical genetics community lacks publicly-available stand-alone software for high-throughput annotation of custom variant data with user-defined tissue-specific epigenetic maps and clinical genetic databases, to prioritize variants for a specific biomedical application. In this work, we provide a simple software pipeline, called SNPnotes, which takes as input variant calls for a patient and prioritizes those using information on clinical relevance from ClinVar, tissue-specific gene regulation from GTEx and disease associations from the NHGRI-EBI GWAS catalogue. This pipeline was developed as part of SVAI Research's "Undiagnosed-1" event for collaborative patient diagnosis. We applied this pipeline to WGS-based variant calls for an individual with a history of gastrointestinal symptoms, using 12 gut-specific eQTL maps and GWAS associations for metabolic diseases, for variant-gene mapping. Out of 6,248,584 SNPs, the pipeline identified 151 high-priority variants, overlapping 129 genes. These top SNPs all have known clinical pathogenicity, modulate gene expression in gut tissues and have genetic associations with metabolic disorders, and serve as starting points for hypotheses about mechanisms driving clinical symptoms. Simple software changes can be made to customize the pipeline for other tissue-specific applications. Future extensions could integrate maps of tissue-specific regulatory elements, higher-order chromatin loops, and mutations affecting splice variants.

VSEAMS: a pipeline for variant set enrichment analysis using summary GWAS data identifies IKZF3, BATF and ESRRA as key transcription factors in type 1 diabetes

From GWAS Variant to Function: A Study of ∼148,000 Variants for Blood Cell Traits

Integrating Regulatory Features Data for Prediction of Functional Disease-Associated SNPs

Functional Variants Identified Efficiently Through An Integrated Transcriptome And Epigenome Analysis

Integrated analysis of variants and pathways in genome-wide association studies using polygenic models of disease

Openness weighted association studies: leveraging personal genome information to prioritize non-coding variants

Variant-to-function dissection of rare non-coding GWAS loci with high impact on blood traits

SNPnotes: high-throughput tissue-specific functional annotation of single nucleotide variants

Integrating Genetic Association, Genetics of Gene Expression, and Single Nucleotide Polymorphism Set Analysis to Identify Susceptibility Loci for Type 2 Diabetes Mellitus.

Identification of novel type 1 and type 2 diabetes genes by co-localization of human islet eQTL and GWAS variants with colocRedRibbon

Genomic annotation of disease-associated variants reveals shared functional contexts

Multivariate Analysis of Genome-Wide Data to Identify Potential Pleiotropic Genes for Type 2 Diabetes, Obesity and Coronary Artery Disease Using MetaCCA

INFIMA leverages multi-omics model organism data to identify effector genes of human GWAS variants

Escan: Scan Regulatory Regions for Aggregate Association Testing Using Whole Genome Sequencing Data

An effector index to predict target genes at GWAS loci

Joint analysis of individual-level and summary-level GWAS data by leveraging pleiotropy

DNA-binding factor footprints and enhancer RNAs identify functional non-coding genetic variants

Addressing the Missing Heritability Problem with the Help of Regulatory Features

GIGSEA: Genotype Imputed Gene Set Enrichment Analysis Using GWAS Summary Level Data

A multilayered post‐GWAS analysis pipeline defines functional variants and target genes for systemic lupus erythematosus (SLE)

Integrating Genome-Wide Association and Eqtls Studies Identifies the Genes and Gene Sets Associated with Diabetes.