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.

Identification of candidate causal cis -regulatory variants underlying electrocardiographic QT interval GWAS loci

Identification of candidate causal -regulatory variants underlying electrocardiographic QT interval GWAS loci

Functional Interrogation Of Enhancer Connectome Prioritizes Candidate Target Genes At Ovarian Cancer Susceptibility Loci

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

Discovering Candidate Genes Regulated by GWAS Signals in Cis and Trans

Genetic Insights into Cardiac Conduction Disorders from Genome-Wide Association Studies

Molecular Mechanisms of Coronary Disease Revealed Using Quantitative Trait Loci for TCF21 Binding, Chromatin Accessibility, and Chromosomal Looping

A Comprehensive cis-eQTL Analysis Revealed Target Genes in Breast Cancer Susceptibility Loci Identified in Genome-wide Association Studies.

Novel Common Variants Associated with Body Mass Index and Coronary Artery Disease Detected Using a Pleiotropic Cfdr Method.

High-resolution genetic mapping of putative causal interactions between regions of open chromatin

High-throughput Screening and CRISPR-Cas9 Modeling of Causal Lipid-associated Expression Quantitative Trait Locus Variants

Distant regulatory effects of genetic variation in multiple human tissues

Dissecting Regulatory Non-Coding Heart Disease GWAS Loci with High-Resolution 3D Chromatin Interactions Reveals Causal Genes with Pathophysiological Relevance to Heart Failure

Functional characterization of QT interval associated enhancer variants identify combined additive effects

Functional Characterization of QT Interval Associatedscn5aenhancer Variants Identify Combined Additive Effects

-eQTL mapping in gene sets identifies network effects of genetic variants

Integrated genomic analysis identifies novel low-frequency <i>cis</i>-regulatory variant rs2279658 associated with VSD risk in Chinese children

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

Impact of disease-associated chromatin accessibility QTLs across immune cell types and contexts

Detailed Regulatory Interaction Map of the Human Heart Facilitates Gene Discovery for Cardiovascular Disease

Functional characterization of QT interval associated SCN5A enhancer variants identify combined additive effects