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.

Integrative annotation scores of variants for impact on RNA binding protein activities

Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project

A method for scoring the cell type-specific impacts of noncoding variants in personal genomes

Meta-analysis of RNA interaction profiles of RNA-binding protein using the RBPInper tool

Improved indel detection in DNA and RNA via realignment with ABRA2

Leveraging Multiple Genomic Data to Prioritize Disease-Causing Indels from Exome Sequencing Data

RBP-Var: a Database of Functional Variants Involved in Regulation Mediated by RNA-binding Proteins.

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

DeFine: Deep Convolutional Neural Networks Accurately Quantify Intensities of Transcription Factor-Dna Binding and Facilitate Evaluation of Functional Non-Coding Variants

Joint genotypic and phenotypic outcome modeling improves base editing variant effect quantification

A general framework for estimating the relative pathogenicity of human genetic variants

A spectral approach integrating functional genomic annotations for coding and noncoding variants

TIVAN-indel: a computational framework for annotating and predicting non-coding regulatory small insertions and deletions

The Value of Statistical or Bioinformatics Annotation for Rare Variant Association with Quantitative Trait.

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

Assisting the analysis of insertions and deletions using regional allele frequencies

CADD v1.7: using protein language models, regulatory CNNs and other nucleotide-level scores to improve genome-wide variant predictions

ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data

Multicenter integrated analysis of noncoding CRISPRi screens

Rank-in: Enabling Integrative Analysis Across Microarray and RNA-seq for Cancer

NAVIP: Unraveling the Influence of Neighboring Small Sequence Variants on Functional Impact Prediction