Abstract:Background: Genetic variation can alter transcriptional regulatory activity contributing to variation in complex traits and risk of disease, but identifying individual variants that affect regulatory activity has been challenging. Quantitative sequence-based experiments such as ChIP-seq and DNase-seq can detect sites of allelic imbalance where alleles contribute disproportionately to the overall signal suggesting allelic differences in regulatory activity.Methods: We created an allelic imbalance detection pipeline, AA-ALIGNER, to remove reference mapping biases influencing allelic imbalance detection and evaluate accuracy of allelic imbalance predictions in the absence of complete genotype data. Using the sequence aligner, GSNAP, and varying amounts of genotype information to remove mapping biases we investigated the accuracy of allelic imbalance detection (binomial test) in CREB1 ChIP-seq reads from the GM12878 cell line. Additionally we thoroughly evaluated the influence of experimental and analytical parameters on imbalance detection.Results: Compared to imbalances identified using complete genotypes, using imputed partial sample genotypes, AA-ALIGNER detected >95 % of imbalances with >90 % accuracy. AA-ALIGNER performed nearly as well using common variants when genotypes were unknown. In contrast, predicting additional heterozygous sites and imbalances using the sequence data led to > 50 % false positive rates. We evaluated effects of experimental data characteristics and key analytical parameter settings on imbalance detection. Overall, total base coverage and signal dispersion across the genome most affected our ability to detect imbalances, while parameters such as imbalance significance, imputation quality thresholds, and alignment mismatches had little effect. To assess the biological relevance of imbalance predictions, we used electrophoretic mobility shift assays to functionally test for predicted allelic differences in CREB1 binding in the GM12878 lymphoblast cell line. Six of nine tested variants exhibited allelic differences in binding. Two of these variants, rs2382818 and rs713875, are located within inflammatory bowel disease-associated loci.Conclusions: AA-ALIGNER accurately detects allelic imbalance in quantitative sequence data using partial genotypes or common variants filling a critical methodological gap in these analyses, as full genotypes are rarely available. Importantly, we demonstrate how experimental and analytical features impact imbalance detection providing guidance for similar future studies.

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

Systematic Analysis of Binding of Transcription Factors to Noncoding Variants.

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

Optimized high-throughput screening of non-coding variants identified from genome-wide association studies

Genome-Wide Identification of Allele-Specific Effects on Gene Expression for Single and Multiple Individuals

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

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

Systematic identification of regulatory variants associated with cancer risk

Removing Reference Mapping Biases Using Limited or No Genotype Data Identifies Allelic Differences in Protein Binding at Disease-Associated Loci

Allelic Expression of Deleterious Protein-Coding Variants Across Human Tissues

High-throughput Screening of Prostate Cancer Risk Loci by Single Nucleotide Polymorphisms Sequencing.

Recommendations for Accurate Resolution of Gene and Isoform Allele-Specific Expression in RNA-Seq Data

Uncovering methylation-dependent genetic effects on regulatory element function in diverse genomes

Regsnps: a Strategy for Prioritizing Regulatory Single Nucleotide Substitutions.

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

Discovery of target genes and pathways of blood trait loci using pooled CRISPR screens and single cell RNA sequencing

Polymorphisms Affecting Gene Transcription and Mrna Processing in Pharmacogenetic Candidate Genes: Detection Through Allelic Expression Imbalance in Human Target Tissues

Determining Allele-Specific Protein Expression (ASPE) Using a QconCAT-Based Proteomics Method: a Novel Approach to Identify Cis-acting Genetic Variants

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

Cisase: a Likelihood-Based Method for Detecting Putative Cis-Regulated Allele-Specific Expression in RNA Sequencing Data.

Prioritizing Single-Nucleotide Variations That Potentially Regulate Alternative Splicing.