Benjamin P. Fairfax

Innate Immune Activity Conditions the Effect of Regulatory Variants upon Monocyte Gene Expression

Introduction Many genetic variants associated with common disease susceptibility occur close to immune-related genes in noncoding DNA, suggestive of a regulatory function. The definition of functional variants and the specific genes that they regulate remains challenging and in many cases is unresolved. We hypothesized that a significant proportion of variants, including those implicated in disease, may show activity in a context-specific manner and therefore only be identifiable upon triggering of immune responses. Context-specific genetic association with differential gene expression in IFN-β signaling. (A) A local association (cis-eQTL) with IFNB1 expression for a single-nucleotide polymorphism (rs2275888) revealed after 2 hours of LPS stimulation of monocytes. (B) This genetic marker shows association with expression of 17 genes on different chromosomes (trans-eQTLs) after 24 hours of LPS stimulation, forming a gene network (C) consistent with the IFN-β signaling cascade. Methods We mapped interindividual variation in gene expression as a quantitative trait, defining expression quantitative trait loci (eQTLs). To investigate the effect of innate immune stimuli on eQTLs, we exposed primary CD14+ human monocytes from 432 European volunteers to the inflammatory proxies interferon-γ (IFN-γ) or differing durations (2 or 24 hours) of lipopolysaccharide (LPS). eQTL mapping was performed on a genome-wide basis with an additive linear model. A subset of 228 individuals with expression data available for all experimental conditions enabled cross-treatment comparisons. Results Stimulation with LPS or IFN-γ resulted in profound effects across monocyte eQTLs, with hundreds of genes and associated pathways demonstrating context-specific eQTLs dependent on the type and duration of stimulus. Context-specific eQTLs frequently intersected established canonical pathways of monocyte signaling and included key nodal genes and effector molecules. These eQTLs are typically more distal to the transcriptional start site and, in some cases, showed reversal of effect between conditions. We also found stimulation reveals novel eQTLs with simultaneous effects involving many genes (trans-eQTLs). Examples included coding polymorphisms in CYP1B1, P2RY11, and IDO2 that modulate activity and develop trans network effects upon stimulation; an LPS-specific IFN-β cytokine network response driven by a cis-eQTL for IFNB1 that was only revealed over time; an interferon regulatory factor 2 (IRF2) transcription factor modulated network up-regulated by IFN-γ involving a cis-eQTL for IRF2; and an IFN-γ–inducible trans gene network involving the transcription factor NFE2L3. We find trans associations to the major histocompatibility complex are dependent on context, paralleling the expression of class II genes. Induced eQTLs were enriched for disease-risk loci with context-specific associations to many putative causal genes, including at ATM, IRF8, and CCR3. Conditional analysis defined additional independent stimulus-specific peaks of association for a given gene. For CARD9 we observed, in addition to a constitutive eQTL informative for a genome-wide association study locus for Crohn’s disease, a stimulus-specific peak eQTL after IFN-γ, defining a further independent signal of disease association. Discussion Interindividual variation in immune responses is accompanied by diverging patterns of gene regulation dependent on underlying genotype. In human monocytes, many regulatory variants display functionality only after pathophysiologically relevant immune stimuli. By considering the cellular and environmental context relevant to disease, it is possible to more extensively resolve functional genetic variants and the specific modulated genes associated with disease. Immune Variation It is difficult to determine the mechanistic consequences of context-dependent genetic variants, some of which may be related to disease (see the Perspective by Gregersen). Two studies now report on the effects of stimulating immunological monocytes and dendritic cells with proteins that can elicit a response to bacterial or viral infection and assess the functional links between genetic variants and profiles of gene expression. M. N. Lee et al. (10.1126/science.1246980) analyzed the expression of more than 400 genes, in dendritic cells from 30 healthy subjects, which revealed how expression quantitative trait loci (eQTLs) affect gene expression within the interferon-β and the Toll-like receptor 3 and 4 pathways. Fairfax et al. (10.1126/science.1246949) performed a genome-wide analysis to show that many eQTLs affected monocyte gene expression in a stimulus- or time-specific manner. Analysis of the transcriptional responses during induced innate immune activity in primary human monocytes is explained. [Also see Perspective by Gregersen] To systematically investigate the impact of immune stimulation upon regulatory variant activity, we exposed primary monocytes from 432 healthy Europeans to interferon-γ (IFN-γ) or differing durations of lipopolysaccharide and mapped expression quantitative trait loci (eQTLs). More than half of cis-eQTLs identified, involving hundreds of genes and associated pathways, are detected specifically in stimulated monocytes. Induced innate immune activity reveals multiple master regulatory trans-eQTLs including the major histocompatibility complex (MHC), coding variants altering enzyme and receptor function, an IFN-β cytokine network showing temporal specificity, and an interferon regulatory factor 2 (IRF2) transcription factor–modulated network. Induced eQTL are significantly enriched for genome-wide association study loci, identifying context-specific associations to putative causal genes including CARD9, ATM, and IRF8. Thus, applying pathophysiologically relevant immune stimuli assists resolution of functional genetic variants.

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

Trans-acting genetic variants have a substantial, albeit poorly characterized, role in the heritable determination of gene expression. Using paired purified primary monocytes and B cells, we identify new predominantly cell type–specific cis and trans expression quantitative trait loci (eQTLs), including multi-locus trans associations to LYZ and KLF4 in monocytes and B cells, respectively. Additionally, we observe a B cell–specific trans association of rs11171739 at 12q13.2, a known autoimmune disease locus, with IP6K2 (P = 5.8 × 10−15), PRIC285 (P = 3.0 × 10−10) and an upstream region of CDKN1A (P = 2 × 10−52), suggesting roles for cell cycle regulation and peroxisome proliferator-activated receptor γ (PPARγ) signaling in autoimmune pathogenesis. We also find that specific human leukocyte antigen (HLA) alleles form trans associations with the expression of AOAH and ARHGAP24 in monocytes but not in B cells. In summary, we show that mapping gene expression in defined primary cell populations identifies new cell type–specific trans-regulated networks and provides insights into the genetic basis of disease susceptibility.

Genetic association analyses implicate aberrant regulation of innate and adaptive immunity genes in the pathogenesis of systemic lupus erythematosus

Systemic lupus erythematosus (SLE) is a genetically complex autoimmune disease characterized by loss of immune tolerance to nuclear and cell surface antigens. Previous genome-wide association studies (GWAS) had modest sample sizes, reducing their scope and reliability. Our study comprised 7,219 cases and 15,991 controls of European ancestry, constituting a new GWAS, a meta-analysis with a published GWAS and a replication study. We have mapped 43 susceptibility loci, including ten new associations. Assisted by dense genome coverage, imputation provided evidence for missense variants underpinning associations in eight genes. Other likely causal genes were established by examining associated alleles for cis-acting eQTL effects in a range of ex vivo immune cells. We found an over-representation (n = 16) of transcription factors among SLE susceptibility genes. This finding supports the view that aberrantly regulated gene expression networks in multiple cell types in both the innate and adaptive immune response contribute to the risk of developing SLE.

Meta-analysis of genome-wide association studies identifies ten loci influencing allergic sensitization

Allergen-specific immunoglobulin E (present in allergic sensitization) has a central role in the pathogenesis of allergic disease. We performed the first large-scale genome-wide association study (GWAS) of allergic sensitization in 5,789 affected individuals and 10,056 controls and followed up the top SNP at each of 26 loci in 6,114 affected individuals and 9,920 controls. We increased the number of susceptibility loci with genome-wide significant association with allergic sensitization from three to ten, including SNPs in or near TLR6, C11orf30, STAT6, SLC25A46, HLA-DQB1, IL1RL1, LPP, MYC, IL2 and HLA-B. All the top SNPs were associated with allergic symptoms in an independent study. Risk-associated variants at these ten loci were estimated to account for at least 25% of allergic sensitization and allergic rhinitis. Understanding the molecular mechanisms underlying these associations may provide new insights into the etiology of allergic disease.

Leprosy and the adaptation of human toll-like receptor 1.

Leprosy is an infectious disease caused by the obligate intracellular pathogen Mycobacterium leprae and remains endemic in many parts of the world. Despite several major studies on susceptibility to leprosy, few genomic loci have been replicated independently. We have conducted an association analysis of more than 1,500 individuals from different case-control and family studies, and observed consistent associations between genetic variants in both TLR1 and the HLA-DRB1/DQA1 regions with susceptibility to leprosy (TLR1 I602S, case-control P = 5.7×10−8, OR = 0.31, 95% CI = 0.20–0.48, and HLA-DQA1 rs1071630, case-control P = 4.9×10−14, OR = 0.43, 95% CI = 0.35–0.54). The effect sizes of these associations suggest that TLR1 and HLA-DRB1/DQA1 are major susceptibility genes in susceptibility to leprosy. Further population differentiation analysis shows that the TLR1 locus is extremely differentiated. The protective dysfunctional 602S allele is rare in Africa but expands to become the dominant allele among individuals of European descent. This supports the hypothesis that this locus may be under selection from mycobacteria or other pathogens that are recognized by TLR1 and its co-receptors. These observations provide insight into the long standing host-pathogen relationship between human and mycobacteria and highlight the key role of the TLR pathway in infectious diseases.

A common haplotype lowers PU.1 expression in myeloid cells and delays onset of Alzheimer's disease

A genome-wide survival analysis of 14,406 Alzheimers disease (AD) cases and 25,849 controls identified eight previously reported AD risk loci and 14 novel loci associated with age at onset. Linkage disequilibrium score regression of 220 cell types implicated the regulation of myeloid gene expression in AD risk. The minor allele of rs1057233 (G), within the previously reported CELF1 AD risk locus, showed association with delayed AD onset and lower expression of SPI1 in monocytes and macrophages. SPI1 encodes PU.1, a transcription factor critical for myeloid cell development and function. AD heritability was enriched within the PU.1 cistrome, implicating a myeloid PU.1 target gene network in AD. Finally, experimentally altered PU.1 levels affected the expression of mouse orthologs of many AD risk genes and the phagocytic activity of mouse microglial cells. Our results suggest that lower SPI1 expression reduces AD risk by regulating myeloid gene expression and cell function.

Genomic modulators of gene expression in human neutrophils

Neutrophils form the most abundant leukocyte subset and are central to many disease processes. Technical challenges in transcriptomic profiling have prohibited genomic approaches to date. Here we map expression quantitative trait loci (eQTL) in peripheral blood CD16+ neutrophils from 101 healthy European adults. We identify cis-eQTL for 3281 neutrophil-expressed genes including many implicated in neutrophil function, with 450 of these not previously observed in myeloid or lymphoid cells. Paired comparison with monocyte eQTL demonstrates nuanced conditioning of genetic regulation of gene expression by cellular context, which relates to cell-type-specific DNA methylation and histone modifications. Neutrophil eQTL are markedly enriched for trait-associated variants particularly autoimmune, allergy and infectious disease. We further demonstrate how eQTL in PADI4 and NOD2 delineate risk variant function in rheumatoid arthritis, leprosy and Crohns disease. Taken together, these data help advance understanding of the genetics of gene expression, neutrophil biology and immune-related diseases.

Distinct Transcriptional and Anti-Mycobacterial Profiles of Peripheral Blood Monocytes Dependent on the Ratio of Monocytes: Lymphocytes.

The ratio of monocytes and lymphocytes (ML ratio) in peripheral blood is associated with tuberculosis and malaria disease risk and cancer and cardiovascular disease outcomes. We studied anti-mycobacterial function and the transcriptome of monocytes in relation to the ML ratio. Mycobacterial growth inhibition assays of whole or sorted blood were performed and mycobacteria were enumerated by liquid culture. Transcriptomes of unstimulated CD14+ monocytes isolated by magnetic bead sorting were characterised by microarray. Transcript expression was tested for association with ML ratio calculated from leucocyte differential counts by linear regression. The ML ratio was associated with mycobacterial growth in vitro (β=2.23, SE 0.91, p=0.02). Using sorted monocytes and lymphocytes, in vivo ML ratio (% variance explained R2=11%, p=0.02) dominated over in vitro ratios (R2=5%, p=0.10) in explaining mycobacterial growth. Expression of 906 genes was associated with the ML ratio and 53 with monocyte count alone. ML-ratio associated genes were enriched for type-I and -II interferon signalling (p=1.2×10-8), and for genes under transcriptional control of IRF1, IRF2, RUNX1, RELA and ESRRB. The ML-ratio-associated gene set was enriched in TB disease (3.11-fold, 95% CI: 2.28-4.19, p=5.7×10-12) and other inflammatory diseases including atopy, HIV, IBD and SLE. The ML ratio is associated with distinct transcriptional and anti-mycobacterial profiles of monocytes that may explain the disease associations of the ML ratio.The ratio of monocytes and lymphocytes (ML ratio) in peripheral blood is associated with tuberculosis and malaria disease risk and cancer and cardiovascular disease outcomes. We studied anti-mycobacterial function and the transcriptome of monocytes in relation to the ML ratio. Mycobacterial growth inhibition assays of whole or sorted blood were performed and mycobacteria were enumerated by liquid culture. Transcriptomes of unstimulated CD14 + monocytes isolated by magnetic bead sorting were characterised by microarray. Transcript expression was tested for association with ML ratio calculated from leucocyte differential counts by linear regression. The ML ratio was associated with mycobacterial growth in vitro (β = 2.23, SE 0.91, p = 0.02). Using sorted monocytes and lymphocytes, in vivo ML ratio (% variance explained R2 = 11%, p = 0.02) dominated over in vitro ratios (R2 = 5%, p = 0.10) in explaining mycobacterial growth. Expression of 906 genes was associated with the ML ratio and 53 with monocyte count alone. ML-ratio associated genes were enriched for type-I and -II interferon signalling (p = 1.2 × 10− 8), and for genes under transcriptional control of IRF1, IRF2, RUNX1, RELA and ESRRB. The ML-ratio-associated gene set was enriched in TB disease (3.11-fold, 95% CI: 2.28–4.19, p = 5.7 × 10− 12) and other inflammatory diseases including atopy, HIV, IBD and SLE. The ML ratio is associated with distinct transcriptional and anti-mycobacterial profiles of monocytes that may explain the disease associations of the ML ratio.

Increased prevalence of sex chromosome aneuploidies in specific language impairment and dyslexia

Sex chromosome aneuploidies increase the risk of spoken or written language disorders but individuals with specific language impairment (SLI) or dyslexia do not routinely undergo cytogenetic analysis. We assess the frequency of sex chromosome aneuploidies in individuals with language impairment or dyslexia.

Pathogenic implications for autoimmune mechanisms derived by comparative eQTL analysis of CD4+ versus CD8+ T cells.

Inappropriate activation or inadequate regulation of CD4+ and CD8+ T cells may contribute to the initiation and progression of multiple autoimmune and inflammatory diseases. Studies on disease-associated genetic polymorphisms have highlighted the importance of biological context for many regulatory variants, which is particularly relevant in understanding the genetic regulation of the immune system and its cellular phenotypes. Here we show cell type-specific regulation of transcript levels of genes associated with several autoimmune diseases in CD4+ and CD8+ T cells including a trans-acting regulatory locus at chr12q13.2 containing the rs1131017 SNP in the RPS26 gene. Most remarkably, we identify a common missense variant in IL27, associated with type 1 diabetes that results in decreased functional activity of the protein and reduced expression levels of downstream IRF1 and STAT1 in CD4+ T cells only. Altogether, our results indicate that eQTL mapping in purified T cells provides novel functional insights into polymorphisms and pathways associated with autoimmune diseases.