Abstract:Background Although high-throughput genotyping arrays have made whole-genome association studies (WGAS) feasible, only a small proportion of SNPs in the human genome are actually surveyed in such studies. In addition, various SNP arrays assay different sets of SNPs, which leads to challenges in comparing results and merging data for meta-analyses. Genome-wide imputation of untyped markers allows us to address these issues in a direct fashion. Methods 384 Caucasian American liver donors were genotyped using Illumina 650Y (Ilmn650Y) arrays, from which we also derived genotypes from the Ilmn317K array. On these data, we compared two imputation methods: MACH and BEAGLE. We imputed 2.5 million HapMap Release22 SNPs, and conducted GWAS on ~40,000 liver mRNA expression traits (eQTL analysis). In addition, 200 Caucasian American and 200 African American subjects were genotyped using the Affymetrix 500 K array plus a custom 164 K fill-in chip. We then imputed the HapMap SNPs and quantified the accuracy by randomly masking observed SNPs. Results MACH and BEAGLE perform similarly with respect to imputation accuracy. The Ilmn650Y results in excellent imputation performance, and it outperforms Affx500K or Ilmn317K sets. For Caucasian Americans, 90% of the HapMap SNPs were imputed at 98% accuracy. As expected, imputation of poorly tagged SNPs (untyped SNPs in weak LD with typed markers) was not as successful. It was more challenging to impute genotypes in the African American population, given (1) shorter LD blocks and (2) admixture with Caucasian populations in this population. To address issue (2), we pooled HapMap CEU and YRI data as an imputation reference set, which greatly improved overall performance. The approximate 40,000 phenotypes scored in these populations provide a path to determine empirically how the power to detect associations is affected by the imputation procedures. That is, at a fixed false discovery rate, the number of cis-eQTL discoveries detected by various methods can be interpreted as their relative statistical power in the GWAS. In this study, we find that imputation offer modest additional power (by 4%) on top of either Ilmn317K or Ilmn650Y, much less than the power gain from Ilmn317K to Ilmn650Y (13%). Conclusion Current algorithms can accurately impute genotypes for untyped markers, which enables researchers to pool data between studies conducted using different SNP sets. While genotyping itself results in a small error rate (e.g. 0.5%), imputing genotypes is surprisingly accurate. We found that dense marker sets (e.g. Ilmn650Y) outperform sparser ones (e.g. Ilmn317K) in terms of imputation yield and accuracy. We also noticed it was harder to impute genotypes for African American samples, partially due to population admixture, although using a pooled reference boosts performance. Interestingly, GWAS carried out using imputed genotypes only slightly increased power on top of assayed SNPs. The reason is likely due to adding more markers via imputation only results in modest gain in genetic coverage, but worsens the multiple testing penalties. Furthermore, cis-eQTL mapping using dense SNP set derived from imputation achieves great resolution, and locate associate peak closer to causal variants than conventional approach.

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

Experimental and Computational Methods for Allelic Imbalance Analysis from Single-Nucleus RNA-seq Data

Probing the limits of cis-acting gene regulation using a model of allelic imbalance quantitative trait loci

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

Defining Loci in Restriction-Based Reduced Representation Genomic Data from Nonmodel Species: Sources of Bias and Diagnostics for Optimal Clustering

Allelic Imbalance (AI) Identifies Novel Tissue-Specific Cis-Regulatory Variation for Human UGT2B15.

Accuracy of genome-wide imputation of untyped markers and impacts on statistical power for association studies

MIXALIME: MIXture models for ALlelic IMbalance Estimation in high-throughput sequencing data

Next-generation DNA sequencing-based assay for measuring allelic expression imbalance (AEI) of candidate neuropsychiatric disorder genes in human brain

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

Imputing Amino Acid Polymorphisms in Human Leukocyte Antigens

Allele-specific binding variants causing ChIP-seq peak height of histone modification are not enriched in expression QTL annotations

Allele-Specific Imbalance Mapping at Human Orthologs of Mouse Susceptibility to Colon Cancer (scc) Loci

A New Approach for Discovering Functional Links Connecting Non-Coding Regulatory Variants to Gene Targets

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

Integrative Analysis of Sequencing and Array Genotype Data for Discovering Disease Associations with Rare Mutations

Abstract B012: Prioritization of ancestry-specific variations in disease data helps to alleviate bias in genomic data

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

A Bayesian Assignment Method for Ambiguous Bisulfite Short Reads

Detecting cell-type-specific allelic expression imbalance by integrative analysis of bulk and single-cell RNA sequencing data

Analyses of Allele-Specific Gene Expression in Highly Divergent Mouse Crosses Identifies Pervasive Allelic Imbalance