Genetics of gene expression in primary immune cells identifies cell type-specific master regulators and roles of HLA alleles.

Shared and cell–specific cis and trans associations in B-cells and monocytes(a) Venn diagrams illustrating the number of eSNPs unique to each cell subset and shared between datasets in cis (left panel) and in trans. B-cell data depicted in orange, monocyte in blue. (b) Circos plots for monocyte dataset (left panel) and B-cell dataset. From outside rim inwards: The outermost rim depicts a Manhattan plot for eQTLs from the respective dataset, the second rim depicts relative expression of genes, the third rim depicts an arbitrary selection of genes (constrained due to space) with significant trans-eQTLs (p<1×10−11) and the innermost network depicts ‘spokes’ between nodal eSNPs and their trans-regulated genes. Such nodal eSNPs include rs10784774 which marks a monocyte specific master-regulatory region while rs10816519 identifies a B-cell specific master regulatory region. Red spokes: > 10 eSNPs from this locus map to the trans-eQTLs, blue spokes: >1 and <10, grey spokes: 1 eSNP.

eSNPs shared between cell-types may lead to opposing directional effects on gene expression or associate with expression of different genes in a cell-specific manner(a) Using eSNPs shared between cell types the most significant eSNP per probe is plotted with the fold change in expression this eSNP causes in homozygous form between major and minor alleles (fold-change monocytes y- axis, B-cells x-axis). Whilst the majority of eSNPs shared between datasets cause the same directional change, there are examples of eSNPs that cause opposing directional changes in expression dependent upon cell-type. Only one eSNP is plotted per probe and only eSNPs with examples of >2 individuals homozygous in the minor allele with permuted p<0.001 in both B-cells and monocyte datasets are annotated. (b) rs2223286 is associated with profound directional effects in the expression of SELL dependent upon genotype, with the minor C allele associated with increased expression of SELL in B-cells and reduced expression of SELL in monocytes (pB-cell=4.6×10−11, pmonocyte=1.1×10−22). (c) rs738289 is an example of an eSNP that forms eQTL to differing genes dependent upon cell type. In B-cells this eSNP is strongly associated with the expression of MGAT3 (pB-cell= 9.8×10−26) with no association to SYNGR1 expression; whilst in monocytes this eSNP is significantly associated with the expression of SYNGR1 (pmonocyte = 1.2×10−17) with no association to MGAT3 expression.

rs10784774 marks a monocyte specific cis-eSNP to LYZ and forms a monocyte specific master regulator of multiple genes including CREB1(a & b) rs10784774 is significantly associated with the expression of a probe mapping to the 3′UTR of LYZ, encoding lysozyme, in monocytes only (p 78 B-cell N.S., pmonocyte=8.9×10−78). (c) rs10784774 additionally forms an eSNP to 72 probes mapping across the genome, the most significantly associated being in the genes CREB1(p = 2.0×10−67), DUSP19 (p = 2.2×10−59), and RAB27A (p = 2.0×10−57). (d) Manhattan plot demonstrating chromosomal location of all genes with probes mapping to rs10784774. The blue line represents 1.0×10−11, red line 5.0×10−8.

rs11171739, a SNP at 12q13.2 with strong autoimmune association forms B-cell specific trans-eQTL to 3 separate genes(a) rs11171739 specifically associates with the B-cell trans expression of the genes LAP3P2 (pB-cell=2.0×10−52, pmonocyte=1.7×10−4), IP6K2 (pB-cell =5.8×10−15) and PRIC285 (pB-cell=3.0×10−10). (b) B-cell regional association plots of these genes demonstrating their expression by SNP marker across chromosome 12q13.2, a known autoimmune risk locus. (c) LAP3P2 lies ~1kb 5′ to CDKN1A and is denoted as a pseudogene. We observe a high degree of transcriptional activity in our RNA-seq data from CD19+ primary B-cells, and this is also seen in LCLs (GM12878)24. This region is highly conserved amongst vertebrates.

Imputation of HLA status resolves cell-specific trans-associated gene expression to carriage of specific classical HLA alleles(a) Regional association plots from monocyte and B-cell data demonstrating the monocyte specific trans association of the chromosome 7 gene AOAH, to the class II MHC gene HLA-DRB1. There is no association in B-cell dataset, whilst the peak eSNP from monocyte dataset is rs28366298 (pmonocyte=1.6×10−43). (b) Imputation of class II alleles to 2 digit resolution demonstrates that the C allele of rs28366298 is specific to DRB1 *04, *07 and *09 alleles and these alleles are associated with reduced AOAH expression. For clarity only homozygotes are plotted with two AOAH expression values per individual (corresponding to each DRB allele) are displayed. (c) The number of DRB1 *04/*07/*09 alleles carried by an individual is significantly associated with reduced expression of AOAH (pB-cell <2.2×10−16, one-way ANOVA) and increased expression of ARHGAP24 (pmonocyte=3.0×10−14, one-way ANOVA). (d) Imputation of DQA1 status to 4 digits demonstrates that the G allele of rs35265698 is significantly associated with reduced expression of DEF8 in both B-cells and monocytes (pB-cell=6.2×10−13, pmonocyte=9.0×10−17). This allele is unique to DQA1*0301 as illustrated. For clarity only homozygotes of rs35265698 are plotted with two DEF8 expression values per individual (corresponding to each DQA allele) displayed. (e) Higher numbers of DQA1*0301 alleles possessed by an individual is significantly associated with reduced DEF8 expression (pB-cell=6.2×10−13, pmonocyte <2.2×10−16)

Cell specific cis-eQTL involving SNP markers associated with diseaseGenes showing significant cis-eQTL that involve SNPs reported at genome-wide significance (p < 5×10−8) in the Catalog of Published Genome-Wide Association Studies (www.genome.gov/GWA studies) (accessed 10th September 2011) or proxy SNPs identified for these disease markers from the 1000 Genomes Project (CEU cohort, r2 > 0.8) are shown. (a) List of genes with shared cis-eQTL/GWAS SNPs for UC, SLE and CD grouped by cell type and excluding HLA genes. Specific examples are described further in Supplementary Information. (b) Examples of genes showing eQTL involving SNPs associated with multiple GWAS traits.