The immune system plays a major role in human health and disease, and understanding genetic causes of interindividual variability of immune responses is vital. We isolated monocytes from 134 genotyped individuals, stimulated the cells with three defined microbe-associated molecular patterns (LPS, MDP, and pppdsRNA), and profiled the transcriptome at three time points. After mapping expression quantitative trait loci (eQTL), we identified 417 response eQTLs (reQTLs) with differing effect between the conditions. We characterized the dynamics of genetic regulation on early and late immune response, and observed an enrichment of reQTLs in distal cis-regulatory elements. Response eQTLs are also enriched for recent positive selection with an evolutionary trend towards enhanced immune response. Finally, we uncover novel reQTL effects in multiple GWAS loci, and show a stronger enrichment of response than constant eQTLs in GWAS signals of several autoimmune diseases. This demonstrates the importance of infectious stimuli modifying genetic predisposition to disease.
Main TextAn increasingly popular approach to identify genetic factors affecting interindividual variation in immune response is mapping expression quantitative trait loci (eQTLs) -variants that associate to gene expression -and to identify so-called response eQTLs (reQTLs) where the eQTL effect differs between immune stimuli [1][2][3][4][5][6] .Such genetic variants can impact the transcriptional response to infection, and also represent genetic effects that are modified by the infectious environment via gene-byenvironment interactions. In this study, we create a data set of a large number of . CC-BY-NC-ND 4.0 International license It is made available under a (which was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.The copyright holder for this preprint . http://dx.doi.org/10.1101/116376 doi: bioRxiv preprint first posted online Mar. 13, 2017; -3 -immune stimulus conditions, with monocytes activated with microbial ligands for three different pattern recognition receptor (PRR) families at two different time points, next to the baseline condition.To examine the time course of innate immune responses, we first profiled gene expression in monocytes of five individuals using Human HT-12 v4 Expression BeadChips (Illumina) at six time points after stimulation with three prototypical microbial ligands: Lipopolysaccharide (LPS) was used to activate TLR4, muramyl-dipeptide (MDP) to stimulate NOD2, and 5'-triphosphate RNA (RNA) to activate RIG-I. Hierarchical clustering revealed early differentially expressed (DE) genes at 45 and 90 minutes after stimulation and late DE genes between 3 and 24 hours ( Supplementary Fig. 1). For the full eQTL cohort, we analyzed primary monocytes isolated from 134 healthy male individuals (185 before quality control), which were either untreated (baseline) or stimulated with the same three pathogenderived stimuli, and gene expression was profiled after 90 ...