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.

Hierarchical annotation of eQTLs by H-eQTL enables identification of genes with cell type-divergent regulation

Integrative modeling of eQTLs and cis-regulatory elements suggest mechanisms underlying cell type specificity of eQTLs

Chromatin Potential Identified by Shared Single Cell Profiling of RNA and Chromatin

Structural Basis of Differential Gene Expression at eQTLs Loci from High-Resolution Ensemble Models of 3D Single-Cell Chromatin Conformations

Haplotype-aware modeling of cis-regulatory effects highlights the gaps remaining in eQTL data

A unified framework for cell-type-specific eQTLs prioritization by integrating bulk and scRNA-seq data

A statistical framework for joint eQTL analysis in multiple tissues

Discovering single-cell eQTLs from scRNA-seq data only

Modeling heterogeneity in single-cell perturbation states enhances detection of response eQTLs

Single-cell Eqtlgen Consortium: a Personalized Understanding of Disease

The Single-Cell Eqtlgen Consortium

Cell-type and dynamic state govern genetic regulation of gene expression in heterogeneous differentiating cultures

Epigenome-augmented eQTL-hotspots reveal genome-wide transcriptional programs in 36 human tissues

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

An analytical framework for decoding cell type-specific genetic variation of gene regulation

Research progress on single-cell expression quantitative trait loci

HT-eQTL: Integrative Expression Quantitative Trait Loci Analysis in a Large Number of Human Tissues

Gene Set Enrichment in Eqtl Data Identifies Novel Annotations and Pathway Regulators.

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

Cross-Population Joint Analysis of eQTLs: Fine Mapping and Functional Annotation

HBI: a hierarchical Bayesian interaction model to estimate cell-type-specific methylation quantitative trait loci incorporating priors from cell-sorted bisulfite sequencing data