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.

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

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-Asb: A Computational Framework for Identifying Allele-Specific Transcription Factor Binding from ATAC-seq Data

Genome-wide Pqtl Analysis of Protein Expression Regulatory Networks in the Human Liver

Allele-specific transcription factor binding across human brain regions offers mechanistic insight into eQTLs

Systematic Analysis of Binding of Transcription Factors to Noncoding Variants.

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

Global reference mapping of human transcription factor footprints

Liver eQTL meta-analysis illuminates potential molecular mechanisms of cardiometabolic traits

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

Omics-based construction of regulatory variants can be applied to help decipher pig liver-related traits

Contextualising transcription factor binding during embryogenesis using natural sequence variation

Distant regulatory effects of genetic variation in multiple human tissues

A New Liver Expression Quantitative Trait Locus Map from 1,183 Individuals Provides Evidence for Novel Expression Quantitative Trait Loci of Drug Response, Metabolic, and Sex-Biased Phenotypes.

Functional dissection of complex and molecular trait variants at single nucleotide resolution

Abstract 242: Functional Characterization of a Cis-eQTL Locus for Plasma Cholesterol Using CRISPR/Cas Genome Editing in Human Pluripotent Stem Cells

Enhancing Disease Risk Gene Discovery by Integrating Transcription Factor-Linked Trans-located Variants into Transcriptome-Wide Association Analyses

Promoter-interacting expression quantitative trait loci are enriched for functional genetic variants

Fine-mapping causal tissues and genes at disease-associated loci

High-throughput characterization of functional variants highlights heterogeneity and polygenicity underlying lung cancer susceptibility