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.

A semisupervised model to predict regulatory effects of genetic variants at single nucleotide resolution using massively parallel reporter assays

Prediction of Functional Regulatory SNPs in Monogenic and Complex Disease.

A semi-supervised deep learning approach for predicting the functional effects of genomic non-coding variations

Regsnps: a Strategy for Prioritizing Regulatory Single Nucleotide Substitutions.

Predicting Regulatory Variants with Composite Statistic

Nested epistasis enhancer networks for robust genome regulation

A Probabilistic Model to Predict Clinical Phenotypic Traits from Genome Sequencing

Leveraging functional genomic annotations and genome coverage to improve polygenic prediction of complex traits within and between ancestries

Single-cell multi-ome regression models identify functional and disease-associated enhancers and enable chromatin potential analysis

A Sequence-Based Method to Predict the Impact of Regulatory Variants Using Random Forest.

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

Regulatory Element-Based Prediction Identifies New Susceptibility Regulatory Variants for Osteoporosis.

Powerful mapping of cis-genetic effects on gene expression across diverse populations reveals novel disease-critical genes

Deep learning sequence-based ab initio prediction of variant effects on expression and disease risk

A rapid epistatic mixed-model association analysis by linear retransformations of genomic estimated values.

A comprehensive investigation of statistical and machine learning approaches for predicting complex human diseases on genomic variants

Regbase: Whole Genome Base-Wise Aggregation and Functional Prediction for Human Non-Coding Regulatory Variants

EMO: Predicting Non-coding Mutation-induced Up- and Down-regulation of Risk Gene Expression using Deep Learning

How to improve polygenic prediction from whole-genome sequencing data by leveraging predicted epigenomic features?

In silico generation and augmentation of regulatory variants from massively parallel reporter assay using conditional variational autoencoder

Comparative analysis of models in predicting the effects of SNPs on TF-DNA binding using large-scale in vitro and in vivo data