Systematic identification of trans eQTLs as putative drivers of known disease associations.

Trans-eQTL SNPs are enriched for functional elementsWe investigated whether the trans-eQTL SNPs are enriched for certain functional elements. We used the online tools SNPInfo, SNPNexus, and HaploReg that rely upon data from, amongst others, the ENCODE project. (a) We observed that trans-eQTL SNPs are enriched for mapping within miRNA binding sites (b) trans-eQTL SNPs show strong enrichment (as annotated using HaploReg) for enhancer regions that are present in K562 (myeloid) and GM12878 (lymphoid) cell-lines (error bars represent one standard deviation).

Independent trans-eQTL effects emanating from the IKZF1 locusSystemic lupus erythematosus SNP rs4917014 and unlinked mean corpuscular volume SNP rs4917014 both affect expression of IKZF1 in cis. rs12718597 affects 50 trans-genes (mostly involved in hemoglobin metabolism) while rs4917014 affects eight different genes in trans: the rs4917014*T risk allele increases expression of genes involved in type I interferon response. At a somewhat lower significance threshold of FDR 0.28 rs4917014*T decreases complement C1QB expression. Both processes are hallmark features of SLE.

Cholesterol SNP rs174546 affects LDLR in transThe rs174546*T allele is known to be associated with a decrease in serum LDL cholesterol and triglycerides levels. It increases the expression levels of three genes in cis, but also increases gene expression levels of LDLR that encodes the LDL receptor.

For 21 complex traits, pairs of unlinked trait-associated SNPs affect the same downstream genesWe observed that for 21 different traits, there were pairs of unlinked SNPs that have previously been reported to be associated with these traits and which also affect exactly the same downstream genes in trans, whereas this is rarely observed when using an equally sized, but permuted list of trans-eQTLs.

Two unlinked type-1 diabetes risk alleles both increase STAT1 and GBP4 expressionrs3184504*T, a risk allele for type 1 diabetes (chromosome 12), affects the expression of SH2B3 in cis, but also affects the expression levels of fourteen unique genes in trans, including two interferon-γ response genes GBP4 and STAT1. Another unlinked type-1 diabetes risk allele (rs4788084*C on chromosome 16) also increases expression levels of these two interferon-γ response genes, indicating that an elevated interferon-γ response is important in type 1 diabetes.