  Genetic variants in cis-regulatory elements or trans-acting regulators
commonly influence the quantity and spatiotemporal distribution of gene
transcription. Recent interest in expression quantitative trait locus (eQTL)
mapping has paralleled the adoption of genome-wide association studies (GWAS)
for the analysis of complex traits and disease in humans. Under the hypothesis
that many GWAS associations tag non-coding SNPs with small effects, and that
these SNPs exert phenotypic control by modifying gene expression, it has become
common to interpret GWAS associations using eQTL data. To exploit the
mechanistic interpretability of eQTL-GWAS comparisons, an improved
understanding of the genetic architecture and cell type specificity of eQTLs is
required. We address this need by performing an eQTL analysis in four parts:
first we identified eQTLs from eleven studies on seven cell types; next we
quantified cell type specific eQTLs across the studies; then we integrated eQTL
data with cis-regulatory element (CRE) data sets from the ENCODE project;
finally we built a classifier to predict cell type specific eQTLs. Consistent
with prior studies, we demonstrate that allelic heterogeneity is pervasive at
cis-eQTLs and that cis-eQTLs are often cell type specific. Within and between
cell type eQTL replication is associated with eQTL SNP overlap with hundreds of
cell type specific CRE element classes, including enhancer, promoter, and
repressive chromatin marks, regions of open chromatin, and many classes of DNA
binding proteins. Using a random forest classifier including 526 CRE data sets
as features, we successfully predict the cell type specificity of eQTL SNPs in
the absence of gene expression data from the cell type of interest. We
anticipate that such integrative, predictive modeling will improve our ability
to understand the mechanistic basis of human complex phenotypic variation.
