Abstract:Expression quantitative trait loci (eQTL) analysis is useful for identifying genetic variants correlated with gene expression, however, it cannot distinguish between causal and nearby non-functional variants. Because the majority of disease-associated SNPs are located in regulatory regions, they can impact allele-specific binding (ASB) of transcription factors and result in differential expression of the target gene alleles. In this study, our aim was to identify functional single-nucleotide polymorphisms (SNPs) that alter transcriptional regulation and thus, potentially impact cellular function. Here, we present regSNPs-ASB, a generalized linear model-based approach to identify regulatory SNPs that are located in transcription factor binding sites. The input for this model includes ATAC-seq (assay for transposase-accessible chromatin with high-throughput sequencing) raw read counts from heterozygous loci, where differential transposase-cleavage patterns between two alleles indicate preferential transcription factor binding to one of the alleles. Using regSNPs-ASB, we identified 53 regulatory SNPs in human MCF-7 breast cancer cells and 125 regulatory SNPs in human mesenchymal stem cells (MSC). By integrating the regSNPs-ASB output with RNA-seq experimental data and publicly available chromatin interaction data from MCF-7 cells, we found that these 53 regulatory SNPs were associated with 74 potential target genes and that 32 (43%) of these genes showed significant allele-specific expression. By comparing all of the MCF-7 and MSC regulatory SNPs to the eQTLs in the Genome-Tissue Expression (GTEx) Project database, we found that 30% (16/53) of the regulatory SNPs in MCF-7 and 43% (52/122) of the regulatory SNPs in MSC were also in eQTL regions. The enrichment of regulatory SNPs in eQTLs indicated that many of them are likely responsible for allelic differences in gene expression (chi-square test, p-value < 0.01). In summary, we conclude that regSNPs-ASB is a useful tool for identifying causal variants from ATAC-seq data. This new computational tool will enable efficient prioritization of genetic variants identified as eQTL for further studies to validate their causal regulatory function. Ultimately, identifying causal genetic variants will further our understanding of the underlying molecular mechanisms of disease and the eventual development of potential therapeutic targets.

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

A New Approach for Discovering Functional Links Connecting Non-Coding Regulatory Variants to Gene Targets

Regsnps: a Strategy for Prioritizing Regulatory Single Nucleotide Substitutions.

Systematic Analysis of Binding of Transcription Factors to Noncoding Variants.

Assessing the function of genetic variants in candidate gene association studies

Characterization of non-coding variants associated with transcription factor binding through ATAC-seq-defined footprint QTLs in liver

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

Global reference mapping of human transcription factor footprints

DeFine: Deep Convolutional Neural Networks Accurately Quantify Intensities of Transcription Factor-Dna Binding and Facilitate Evaluation of Functional Non-Coding Variants

Regsnps-Asb: A Computational Framework for Identifying Allele-Specific Transcription Factor Binding from ATAC-seq Data

Understanding the Function of Regulatory DNA Interactions in the Interpretation of Non-Coding GWAS Variants

Functional noncoding SNPs in human endothelial cells fine-map vascular trait associations

Leveraging polygenic functional enrichment to improve GWAS power

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

Incorporating functional annotation information in prioritizing disease associated SNPs from genome wide association studies

A method for scoring the cell type-specific impacts of noncoding variants in personal genomes

Genome-wide enhancer maps link risk variants to disease genes

How chromatin interactions shed light on interpreting non-coding genomic variants: opportunities and future direc-tions

High-throughput characterization of genetic effects on DNA-protein binding and gene transcription

Genomic annotation of disease-associated variants reveals shared functional contexts

Scalable approaches for functional analyses of whole-genome sequencing non-coding variants