Characterizing the genetic basis of innate immune response in TLR4-activated human monocytes

A modeling framework for detecting and leveraging node-level information in Bayesian network inference

Bayesian graphical models are powerful tools to infer complex relationships in high dimension, yet are often fraught with computational and statistical challenges. If exploited in a principled way, the increasing information collected alongside the data of primary interest constitutes an opportunity to mitigate these difficulties by guiding the detection of dependence structures. For instance, gene network inference may be informed by the use of publicly available summary statistics on the regulation of genes by genetic variants. Here we present a novel Gaussian graphical modeling framework to identify and leverage information on the centrality of nodes in conditional independence graphs. Specifically, we consider a fully joint hierarchical model to simultaneously infer (i) sparse precision matrices and (ii) the relevance of node-level information for uncovering the sought-after network structure. We encode such information as candidate auxiliary variables using a spike-and-slab submodel on the propensity of nodes to be hubs, which allows hypothesis-free selection and interpretation of a sparse subset of relevant variables. As efficient exploration of large posterior spaces is needed for real-world applications, we develop a variational expectation conditional maximization algorithm that scales inference to hundreds of samples, nodes and auxiliary variables. We illustrate and exploit the advantages of our approach in simulations and in a gene network study which identifies hub genes involved in biological pathways relevant to immune-mediated diseases.

Faecalibacterium prausnitzii promotes intestinal epithelial IL-18 production through activation of the HIF1α pathway

Introduction

Intestinal epithelial cells produce interleukin-18 (IL-18), a key factor in promoting epithelial barrier integrity. Here, we analyzed the potential role of gut bacteria and the hypoxia-inducible factor 1α (HIF1α) pathway in regulating mucosal IL18 expression in inflammatory bowel disease (IBD).

Methods

Mucosal samples from patients with IBD (n = 760) were analyzed for bacterial composition, IL18 levels and HIF1α pathway activation. Wild-type Caco-2 and CRISPR/Cas9-engineered Caco-2-HIF1A-null cells were cocultured with Faecalibacterium prausnitzii in a "Human oxygen-Bacteria anaerobic" in vitro system and analyzed by RNA sequencing.

Results

Mucosal IL18 mRNA levels correlated positively with the abundance of mucosal-associated butyrate-producing bacteria, in particular F. prausnitzii, and with HIF1α pathway activation in patients with IBD. HIF1α-mediated expression of IL18, either by a pharmacological agonist (dimethyloxallyl glycine) or F. prausnitzii, was abrogated in Caco-2-HIF1A-null cells.

Conclusion

Butyrate-producing gut bacteria like F. prausnitzii regulate mucosal IL18 expression in a HIF1α-dependent manner that may aid in mucosal healing in IBD.

Pathogen-specific innate immune response patterns are distinctly affected by genetic diversity

Innate immune responses vary by pathogen and host genetics. We analyze quantitative trait loci (eQTLs) and transcriptomes of monocytes from 215 individuals stimulated by fungal, Gram-negative or Gram-positive bacterial pathogens. We identify conserved monocyte responses to bacterial pathogens and a distinct antifungal response. These include 745 response eQTLs (reQTLs) and corresponding genes with pathogen-specific effects, which we find first in samples of male donors and subsequently confirm for selected reQTLs in females. reQTLs affect predominantly upregulated genes that regulate immune response via e.g., NOD-like, C-type lectin, Toll-like and complement receptor-signaling pathways. Hence, reQTLs provide a functional explanation for individual differences in innate response patterns. Our identified reQTLs are also associated with cancer, autoimmunity, inflammatory and infectious diseases as shown by external genome-wide association studies. Thus, reQTLs help to explain interindividual variation in immune response to infection and provide candidate genes for variants associated with a range of diseases.

Acute systemic myeloid inflammatory and stress response in severe food allergic reactions

INTRODUCTION

Food allergic reactions can be severe and potentially life-threatening and the underlying immunological processes that contribute to the severity of reactions are poorly understood. The aim of this study is to integrate bulk RNA-sequencing of human and mouse peripheral blood mononuclear cells during food allergic reactions and in vivo mouse models of food allergy to identify dysregulated immunological processes associated with severe food allergic reactions.

METHODS

Bulk transcriptomics of whole blood from human and mouse following food allergic reactions combined with integrative differential expressed gene bivariate and module eigengene network analyses to identify the whole blood transcriptome associated with food allergy severity. In vivo validation immune cell and gene expression in mice following IgE-mediated reaction.

RESULTS

Bulk transcriptomics of whole blood from mice with different severity of food allergy identified gene ontology (GO) biological processes associated with innate and inflammatory immune responses, dysregulation of MAPK and NFkB signalling and identified 429 genes that correlated with reaction severity. Utilizing two independent human cohorts, we identified 335 genes that correlated with severity of peanut-induced food allergic reactions. Mapping mouse food allergy severity transcriptome onto the human transcriptome revealed 11 genes significantly dysregulated and correlated with severity. Analyses of whole blood from mice undergoing an IgE-mediated reaction revealed a rapid change in blood leukocytes particularly inflammatory monocytes (Ly6Chi Ly6G- ) and neutrophils that was associated with changes in CLEC4E, CD218A and GPR27 surface expression.

CONCLUSIONS

Collectively, IgE-mediated food allergy severity is associated with a rapid innate inflammatory response associated with acute cellular stress processes and dysregulation of peripheral blood inflammatory myeloid cell frequencies.

Exploration of genotype-by-environment interactions affecting gene expression responses in porcine immune cells

As one of the keys to healthy performance, robustness of farm animals is gaining importance, and with this comes increasing interest in genetic dissection of genotype-by-environment interactions (G×E). Changes in gene expression are among the most sensitive responses conveying adaptation to environmental stimuli. Environmentally responsive regulatory variation thus likely plays a central role in G×E. In the present study, we set out to detect action of environmentally responsive cis-regulatory variation by the analysis of condition-dependent allele specific expression (cd-ASE) in porcine immune cells. For this, we harnessed mRNA-sequencing data of peripheral blood mononuclear cells (PBMCs) stimulated in vitro with lipopolysaccharide, dexamethasone, or their combination. These treatments mimic common challenges such as bacterial infection or stress, and induce vast transcriptome changes. About two thirds of the examined loci showed significant ASE in at least one treatment, and out of those about ten percent exhibited cd-ASE. Most of the ASE variants were not yet reported in the PigGTEx Atlas. Genes showing cd-ASE were enriched in cytokine signaling in immune system and include several key candidates for animal health. In contrast, genes showing no ASE featured cell-cycle related functions. We confirmed LPS-dependent ASE for one of the top candidates, SOD2, which ranks among the major response genes in LPS-stimulated monocytes. The results of the present study demonstrate the potential of in vitro cell models coupled with cd-ASE analysis for the investigation of G×E in farm animals. The identified loci may benefit efforts to unravel the genetic basis of robustness and improvement of health and welfare in pigs.

Using epigenomics to understand cellular responses to environmental influences in diseases

It is a generally accepted model that environmental influences can exert their effects, at least in part, by changing the molecular regulators of transcription that are described as epigenetic. As there is biochemical evidence that some epigenetic regulators of transcription can maintain their states long term and through cell division, an epigenetic model encompasses the idea of maintenance of the effect of an exposure long after it is no longer present. The evidence supporting this model is mostly from the observation of alterations of molecular regulators of transcription following exposures. With the understanding that the interpretation of these associations is more complex than originally recognised, this model may be oversimplistic; therefore, adopting novel perspectives and experimental approaches when examining how environmental exposures are linked to phenotypes may prove worthwhile. In this review, we have chosen to use the example of nonalcoholic fatty liver disease (NAFLD), a common, complex human disease with strong environmental and genetic influences. We describe how epigenomic approaches combined with emerging functional genetic and single-cell genomic techniques are poised to generate new insights into the pathogenesis of environmentally influenced human disease phenotypes exemplified by NAFLD.

IL7 genetic variation and toxicity to immune checkpoint blockade in patients with melanoma

Treatment with immune checkpoint blockade (ICB) frequently triggers immune-related adverse events (irAEs), causing considerable morbidity. In 214 patients receiving ICB for melanoma, we observed increased severe irAE risk in minor allele carriers of rs16906115, intronic to IL7. We found that rs16906115 forms a B cell-specific expression quantitative trait locus (eQTL) to IL7 in patients. Patients carrying the risk allele demonstrate increased pre-treatment B cell IL7 expression, which independently associates with irAE risk, divergent immunoglobulin expression and more B cell receptor mutations. Consistent with the role of IL-7 in T cell development, risk allele carriers have distinct ICB-induced CD8+ T cell subset responses, skewing of T cell clonality and greater proportional repertoire occupancy by large clones. Finally, analysis of TCGA data suggests that risk allele carriers independently have improved melanoma survival. These observations highlight key roles for B cells and IL-7 in both ICB response and toxicity and clinical outcomes in melanoma.

Redefining critical illness

Research and practice in critical care medicine have long been defined by syndromes, which, despite being clinically recognizable entities, are, in fact, loose amalgams of heterogeneous states that may respond differently to therapy. Mounting translational evidence-supported by research on respiratory failure due to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection-suggests that the current syndrome-based framework of critical illness should be reconsidered. Here we discuss recent findings from basic science and clinical research in critical care and explore how these might inform a new conceptual model of critical illness. De-emphasizing syndromes, we focus on the underlying biological changes that underpin critical illness states and that may be amenable to treatment. We hypothesize that such an approach will accelerate critical care research, leading to a richer understanding of the pathobiology of critical illness and of the key determinants of patient outcomes. This, in turn, will support the design of more effective clinical trials and inform a more precise and more effective practice at the bedside.

Functional Characterization of Genetic Variant Effects on Expression

Thousands of common genetic variants in the human population have been associated with disease risk and phenotypic variation by genome-wide association studies (GWAS). However, the majority of GWAS variants fall into noncoding regions of the genome, complicating our understanding of their regulatory functions, and few molecular mechanisms of GWAS variant effects have been clearly elucidated. Here, we set out to review genetic variant effects, focusing on expression quantitative trait loci (eQTLs), including their utility in interpreting GWAS variant mechanisms. We discuss the interrelated challenges and opportunities for eQTL analysis, covering determining causal variants, elucidating molecular mechanisms of action, and understanding context variability. Addressing these questions can enable better functional characterization of disease-associated loci and provide insights into fundamental biological questions of the noncoding genetic regulatory code and its control of gene expression. Expected final online publication date for the Annual Review of Biomedical Data Science, Volume 5 is August 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Mammalian Neuraminidases in Immune-Mediated Diseases: Mucins and Beyond

Mammalian neuraminidases (NEUs), also known as sialidases, are enzymes that cleave off the terminal neuraminic, or sialic, acid resides from the carbohydrate moieties of glycolipids and glycoproteins. A rapidly growing body of literature indicates that in addition to their metabolic functions, NEUs also regulate the activity of their glycoprotein targets. The simple post-translational modification of NEU protein targets-removal of the highly electronegative sialic acid-affects protein folding, alters protein interactions with their ligands, and exposes or covers proteolytic sites. Through such effects, NEUs regulate the downstream processes in which their glycoprotein targets participate. A major target of desialylation by NEUs are mucins (MUCs), and such post-translational modification contributes to regulation of disease processes. In this review, we focus on the regulatory roles of NEU-modified MUCs as coordinators of disease pathogenesis in fibrotic, inflammatory, infectious, and autoimmune diseases. Special attention is placed on the most abundant and best studied NEU1, and its recently discovered important target, mucin-1 (MUC1). The role of the NEU1 - MUC1 axis in disease pathogenesis is discussed, along with regulatory contributions from other MUCs and other pathophysiologically important NEU targets.

17q12-21 risk-variants influence cord blood immune regulation and multitrigger-wheeze

BACKGROUND

Childhood wheeze represents a first symptom of asthma. Early identification of children at risk for wheeze related to 17q12-21 variants and their underlying immunological mechanisms remain unknown. We aimed to assess the influence of 17q12-21 variants and mRNA expression at birth on the development of wheeze.

METHODS

Children were classified as multitrigger/viral/no wheeze until six years of age. The PAULINA/PAULCHEN birth cohorts were genotyped (n = 216; GSA-chip). mRNA expression of 17q21 and innate/adaptive genes was measured (qRT-PCR) in cord blood mononuclear cells. Expression quantitative trait loci (eQTL) and mediation analyses were performed. Genetic variation of 17q12-21 asthma-single nucleotide polymorphisms (SNPs) was summarized as the first principal component (PC1) and used to classify single SNP effects on gene expression as (locus)-dependent/independent eQTL SNPs.

RESULTS

Core region risk variants (IKZF3, ZPBP2, GSDMB, ORMDL3) were associated with multitrigger wheeze (OR: 3.05-5.43) and were locus-dependent eQTL SNPs with higher GSDMA, TLR2, TLR5, and lower TGFB1 expression. Increased risk of multitrigger wheeze with rs9303277 was in part mediated by TLR2 expression. Risk variants distal to the core region were mainly locus-independent eQTL SNPs with decreased CD209, CD86, TRAF6, RORA, and IL-9 expression. Distinct immune signatures in cord blood were associated either with multitrigger wheeze (increased innate genes, e.g., TLR2, IPS1, LY75) or viral wheeze (decreased NF-κB genes, e.g., TNFAIP3 and TNIP2).

CONCLUSION

Locus-dependent eQTL SNPs (core region) associated with increased inflammatory genes (primarily TLR2) at birth and subsequent multitrigger wheeze indicate that early priming and imbalance may be crucial for asthma pathophysiology. Locus-independent eQTL SNPs (mainly distal region, rs1007654) may be involved in the initiation of dendritic cell activation/maturation (TRAF6) and interaction with T cells (CD209, CD86). Identifying potential mechanistic pathways at birth may point to critical key points during early immune development predisposing to asthma.

An autoimmune disease risk variant: A trans master regulatory effect mediated by IRF1 under immune stimulation?

Functional mechanisms remain unknown for most genetic loci associated to complex human traits and diseases. In this study, we first mapped trans-eQTLs in a data set of primary monocytes stimulated with LPS, and discovered that a risk variant for autoimmune disease, rs17622517 in an intron of C5ORF56, affects the expression of the transcription factor IRF1 20 kb away. The cis-regulatory effect specific to IRF1 is active under early immune stimulus, with a large number of trans-eQTL effects across the genome under late LPS response. Using CRISPRi silencing, we showed that perturbation of the SNP locus downregulates IRF1 and causes widespread transcriptional effects. Genome editing by CRISPR had suggestive recapitulation of the LPS-specific trans-eQTL signal and lent support for the rs17622517 site being functional. Our results suggest that this common genetic variant affects inter-individual response to immune stimuli via regulation of IRF1. For this autoimmune GWAS locus, our work provides evidence of the functional variant, demonstrates a condition-specific enhancer effect, identifies IRF1 as the likely causal gene in cis, and indicates that overactivation of the downstream immune-related pathway may be the cellular mechanism increasing disease risk. This work not only provides rare experimental validation of a master-regulatory trans-eQTL, but also demonstrates the power of eQTL mapping to build mechanistic hypotheses amenable for experimental follow-up using the CRISPR toolkit.

Although many genetic loci have been associated to disease, understanding how these variants impact molecular and cellular functions to impact disease risk have been challenging. Here, we first used blood cells from a large number of individuals and stimulated them in the laboratory with a proxy for bacterial infection. We identified that a genetic variant associated to autoimmune diseases also affects the expression of the nearby transcription factor IRF1 gene in early immune response, followed by expression change of other genes in late immune response. We then studied this effect in cell lines, using the CRISPR approach to silence the activity of the genomic element of this variant and cause mutations at that position. We found evidence that this autoimmune disease -associated variant is located in a genomic regulatory element that responds to immune stimulus and affects expression of IRF1 and a complex gene regulatory network. Thus, our characterization of genetic regulatory variation in the human population combined with experimental follow-up suggests a plausible, previously uncharacterized molecular mechanism that may underlie this genetic variant’s effect on immune disease risk.

Dendritic Cells: Neglected Modulators of Peripheral Immune Responses and Neuroinflammation in Mood Disorders?

Affective disorders (AD) including major depressive disorder (MDD) and bipolar disorder (BD) are common mood disorders associated with increased disability and poor health outcomes. Altered immune responses characterized by increased serum levels of pro-inflammatory cytokines and neuroinflammation are common findings in patients with AD and in corresponding animal models. Dendritic cells (DCs) represent a heterogeneous population of myeloid cells that orchestrate innate and adaptive immune responses and self-tolerance. Upon sensing exogenous and endogenous danger signals, mature DCs secrete proinflammatory factors, acquire migratory and antigen presenting capacities and thus contribute to neuroinflammation in trauma, autoimmunity, and neurodegenerative diseases. However, little is known about the involvement of DCs in the pathogenesis of AD. In this review, we summarize the current knowledge on DCs in peripheral immune responses and neuroinflammation in MDD and BD. In addition, we consider the impact of DCs on neuroinflammation and behavior in animal models of AD. Finally, we will discuss therapeutic perspectives targeting DCs and their effector molecules in mood disorders.

[Role of genetic factors in pneumonia and COVID-19]

Die Pneumonie ist die Infektionskrankheit mit der weltweit höchsten Mortalität. Die häufigsten Erreger sind Bakterien, es gibt jedoch auch epidemisch oder pandemisch auftretende virale Lungenentzündungen durch Influenza- oder Coronaviren, wie die aktuelle Pandemie durch das SARS Coronavirus 3766 Fälle (SARS-CoV-2). Wichtige Herausforderungen liegen neben dem Auftreten von Antibiotikaresistenzen und Immunpathologien etwa in der Sepsis in der Betrachtung der Suszeptibilität individueller Patienten: Hier werden vor allen Dingen das Lebensalter, Medikamente und Komorbiditäten betrachtet. Es gibt jedoch auch klare Hinweise für genetische Einflüsse auf das individuelle Risiko, an einer Pneumonie zu erkranken oder einen schweren Verlauf der Erkrankung zu entwickeln. In diesem Beitrag wollen wir die genetischen Einflüsse auf die Pneumonie und ihre klinische Bedeutung darstellen.

Functional and Transcriptional Adaptations of Blood Monocytes Recruited to the Cystic Fibrosis Airway Microenvironment In Vitro

Cystic fibrosis (CF) lung disease is dominated by the recruitment of myeloid cells (neutrophils and monocytes) from the blood which fail to clear the lung of colonizing microbes. In prior in vitro studies, we showed that blood neutrophils migrated through the well-differentiated lung epithelium into the CF airway fluid supernatant (ASN) mimic the dysfunction of CF airway neutrophils in vivo, including decreased bactericidal activity despite an increased metabolism. Here, we hypothesized that, in a similar manner to neutrophils, blood monocytes undergo significant adaptations upon recruitment to CFASN. To test this hypothesis, primary human blood monocytes were transmigrated in our in vitro model into the ASN from healthy control (HC) or CF subjects to mimic in vivo recruitment to normal or CF airways, respectively. Surface phenotype, metabolic and bacterial killing activities, and transcriptomic profile by RNA sequencing were quantified post-transmigration. Unlike neutrophils, monocytes were not metabolically activated, nor did they show broad differences in activation and scavenger receptor expression upon recruitment to the CFASN compared to HCASN. However, monocytes recruited to CFASN showed decreased bactericidal activity. RNASeq analysis showed strong effects of transmigration on monocyte RNA profile, with differences between CFASN and HCASN conditions, notably in immune signaling, including lower expression in the former of the antimicrobial factor ISG15, defensin-like chemokine CXCL11, and nitric oxide-producing enzyme NOS3. While monocytes undergo qualitatively different adaptations from those seen in neutrophils upon recruitment to the CF airway microenvironment, their bactericidal activity is also dysregulated, which could explain why they also fail to protect CF airways from infection.

IRF5 Acts as a Potential Therapeutic Marker in Inflammatory Bowel Diseases

BACKGROUND

Inflammatory bowel diseases (IBDs), including ulcerative colitis (UC) and Crohn's disease (CD), are chronic inflammatory disorders. As is well known, interferon regulatory factor (IRF) 5 is closely associated with the pathogenesis of various inflammatory diseases. But the exact role of IRF5 in IBD remains unclear.

METHODS

In this study, we detected IRF5 expression in peripheral blood mononuclear cells (PBMCs) and inflamed mucosa from IBD patients by immunohistochemistry, western blot, and quantitative real-time polymerase chain reaction. Peripheral blood CD4+ T cells were stimulated with inflammatory cytokines and transfected by lentivirus.

RESULTS

In active IBD patients, the expression of IRF5 in PBMCs and inflamed colonic tissues was obviously increased and significantly associated with disease activity. Ectopic overexpression of IRF5 could promote the differentiation of IBD CD4+ T cells into Th1 and Th17 cells by regulating T-bet and RAR related orphan receptor C, whereas knockdown of IRF5 had the opposite effects. Tumor necrosis factor (TNF)-α upregulated expression of IRF5 in CD4+ T cells, but anti-TNF treatment with infliximab could markedly reduce IRF5 expression in CD4+ T cells and intestinal mucosa of CD patients.

CONCLUSION

Our study reveals a novel mechanism that IRF5 levels are correlated with disease activity in IBD and might function as a possible marker for the management of IBD via regulating Th1 and Th17 immune responses and cytokine production.

The role of toll-like receptor 4 (TLR4) in cancer progression: A possible therapeutic target?

The toll-like receptor (TLR) family consists of vital receptors responsible for pattern recognition in innate immunity, making them the core proteins involved in pathogen detection and eliciting immune responses. The most studied member of this family, TLR4, has been the center of attention regarding its contributory role in many inflammatory diseases including sepsis shock and asthma. Notably, mounting pieces of evidence have proved that this receptor is aberrantly expressed on the tumor cells and the tumor microenvironment in a wide range of cancer types and it is highly associated with the initiation of tumorigenesis as well as tumor progression and drug resistance. Cancer therapy using TLR4 inhibitors has recently drawn scientists' attention, and the promising results of such studies may pave the way for more investigation in the foreseeable future. This review will introduce the key proteins of the TLR4 pathway and how they interact with major growth factors in the tumor microenvironment. Moreover, we will discuss the many aspects of tumor progression affected by the activation of this receptor and provide an overview of the recent therapeutic approaches using various TLR4 antagonists.

An Imperative Need for Further Genetic Studies of Alopecia Areata

Human genetic studies of diseases that are multifactorial and prevalent have generated a wealth of knowledge about the genetic architecture of chronic diseases. Generalizable attributes are shaping the development of models to explain how the human genome influences our health and can be leveraged to improve it. Importantly, both rare and common genetic variants contribute to disease risk and provide complementary information. Although initial genetic studies of alopecia areata have yielded insight with high clinical impact, there remains a number of important unanswered questions pertaining to disease biology and patient care that could be addressed by further genetic investigations.

Multiple sclerosis risk variants regulate gene expression in innate and adaptive immune cells

At least 200 single-nucleotide polymorphisms (SNPs) are associated with multiple sclerosis (MS) risk. A key function that could mediate SNP-encoded MS risk is their regulatory effects on gene expression. We performed microarrays using RNA extracted from purified immune cell types from 73 untreated MS cases and 97 healthy controls and then performed Cis expression quantitative trait loci mapping studies using additive linear models. We describe MS risk expression quantitative trait loci associations for 129 distinct genes. By extending these models to include an interaction term between genotype and phenotype, we identify MS risk SNPs with opposing effects on gene expression in cases compared with controls, namely, rs2256814 MYT1 in CD4 cells (q = 0.05) and rs12087340 RF00136 in monocyte cells (q = 0.04). The rs703842 SNP was also associated with a differential effect size on the expression of the METTL21B gene in CD8 cells of MS cases relative to controls (q = 0.03). Our study provides a detailed map of MS risk loci that function by regulating gene expression in cell types relevant to MS.

A Global-Local Approach for Detecting Hotspots in Multiple-Response Regression

We tackle modelling and inference for variable selection in regression problems with many predictors and many responses. We focus on detecting hotspots, that is, predictors associated with several responses. Such a task is critical in statistical genetics, as hotspot genetic variants shape the architecture of the genome by controlling the expression of many genes and may initiate decisive functional mechanisms underlying disease endpoints. Existing hierarchical regression approaches designed to model hotspots suffer from two limitations: their discrimination of hotspots is sensitive to the choice of top-level scale parameters for the propensity of predictors to be hotspots, and they do not scale to large predictor and response vectors, for example, of dimensions 103-105 in genetic applications. We address these shortcomings by introducing a flexible hierarchical regression framework that is tailored to the detection of hotspots and scalable to the above dimensions. Our proposal implements a fully Bayesian model for hotspots based on the horseshoe shrinkage prior. Its global-local formulation shrinks noise globally and, hence, accommodates the highly sparse nature of genetic analyses while being robust to individual signals, thus leaving the effects of hotspots unshrunk. Inference is carried out using a fast variational algorithm coupled with a novel simulated annealing procedure that allows efficient exploration of multimodal distributions.

CD16+ monocytes give rise to CD103+RALDH2+TCF4+ dendritic cells with unique transcriptional and immunological features

Classical CD16^- vs intermediate/nonclassical CD16⁺ monocytes differ in their homing potential and biological functions, but whether they differentiate into dendritic cells (DCs) with distinct contributions to immunity against bacterial/viral pathogens remains poorly investigated. Here, we employed a systems biology approach to identify clinically relevant differences between CD16⁺ and CD16^- monocyte-derived DCs (MDDCs). Although both CD16⁺ and CD16^- MDDCs acquire classical immature/mature DC markers in vitro, genome-wide transcriptional profiling revealed unique molecular signatures for CD16⁺ MDDCs, including adhesion molecules (ITGAE/CD103), transcription factors (TCF7L2/TCF4), and enzymes (ALDH1A2/RALDH2), whereas CD16^- MDDCs exhibit a CDH1/E-cadherin⁺ phenotype. Of note, lipopolysaccharides (LPS) upregulated distinct transcripts in CD16⁺ (eg, CCL8, SIGLEC1, MIR4439, SCIN, interleukin [IL]-7R, PLTP, tumor necrosis factor [TNF]) and CD16^- MDDCs (eg, MMP10, MMP1, TGM2, IL-1A, TNFRSF11A, lysosomal-associated membrane protein 1, MMP8). Also, unique sets of HIV-modulated genes were identified in the 2 subsets. Further gene set enrichment analysis identified canonical pathways that pointed to "inflammation" as the major feature of CD16⁺ MDDCs at immature stage and on LPS/HIV exposure. Finally, functional validations and meta-analysis comparing the transcriptome of monocyte and MDDC subsets revealed that CD16⁺ vs CD16^- monocytes preserved their superior ability to produce TNF-α and CCL22, as well as other sets of transcripts (eg, TCF4), during differentiation into DC. These results provide evidence that monocyte subsets are transcriptionally imprinted/programmed with specific differentiation fates, with intermediate/nonclassical CD16⁺ monocytes being precursors for pro-inflammatory CD103⁺RALDH2⁺TCF4⁺ DCs that may play key roles in mucosal immunity homeostasis/pathogenesis. Thus, alterations in the CD16⁺ /CD16^- monocyte ratios during pathological conditions may dramatically influence the quality of MDDC-mediated immunity.

Intestinal and Limb Ischemic Preconditioning Provides a Combined Protective Effect in the Late Phase, But not in the Early Phase, Against Intestinal Injury Induced by Intestinal Ischemia-Reperfusion in Rats

Intestinal ischemia/reperfusion (I/R) injury is associated with high morbidity and mortality. This study aimed to compare the protective efficacy of intestinal ischemic preconditioning (IIPC) and limb ischemic preconditioning (LIPC) against intestinal I/R injury and investigate their combined protective effect and the underlying mechanism. Male Sprague-Dawley rats were pretreated with IIPC, LIPC, or IIPC plus LIPC (combined), and intestinal I/R or sham operation was performed. The animals were sacrificed at 2 and 24 h after reperfusion and then blood and tissue samples were harvested for further analyses. In additional groups of animals, a 7-day survival study was conducted. The results showed that ischemic preconditioning (IPC) improved the survival rate and attenuated intestinal edema, injury, and apoptosis. IPC decreased the levels of tumor necrosis factor-α, interleukin -6, malondialdehyde and myeloperoxidase, and increased the activity of superoxide dismutase in serum and intestine after the I/R event. IPC downregulated the expression of Toll-like receptor-4 (TLR4) and nuclear factor-kappa B (NF-κB). The effect of combined pretreatment was better than that of single pretreatment in the late phase (24 h), but not in the early phase (2 h). The study demonstrated that IPC could significantly attenuate intestinal injury induced by intestinal I/R via inhibiting inflammation, oxidative stress, and apoptosis. IIPC and LIPC conferred no synergy in protecting I/R-induced intestinal injury in the early phase, but combined preconditioning had clearly stronger protection in the late phase, which was associated with the inhibition of the activated TLR4/NF-κB signaling pathway. It suggested that LIPC or combined preconditioning could potentially be applied in the clinical settings of surgical patient care.

Multi-Omics Strategies Uncover Host-Pathogen Interactions

With the success of the Human Genome Project, large-scale systemic projects became a reality that enabled rapid development of the systems biology field. Systems biology approaches to host-pathogen interactions have been instrumental in the discovery of some specifics of Gram-negative bacterial recognition, host signal transduction, and immune tolerance. However, further research, particularly using multi-omics approaches, is essential to untangle the genetic, immunologic, (post)transcriptional, (post)translational, and metabolic mechanisms underlying progression from infection to clearance of microbes. The key to understanding host-pathogen interactions lies in acquiring, analyzing, and modeling multimodal data obtained through integrative multi-omics experiments. In this article, we will discuss how multi-omics analyses are adding to our understanding of the molecular basis of host-pathogen interactions and systemic maladaptive immune response of the host to microbes and microbial products.

Identification of trans-eQTLs using mediation analysis with multiple mediators

Background

Mapping expression quantitative trait loci (eQTLs) has provided insight into gene regulation. Compared to cis-eQTLs, the regulatory mechanisms of trans-eQTLs are less known. Previous studies suggest that trans-eQTLs may regulate expression of remote genes by altering the expression of nearby genes. Trans-association has been studied in the mediation analysis with a single mediator. However, prior applications with one mediator are prone to model misspecification due to correlations between genes. Motivated from the observation that trans-eQTLs are more likely to associate with more than one cis-gene than randomly selected SNPs in the GTEx dataset, we developed a computational method to identify trans-eQTLs that are mediated by multiple mediators.

Results

We proposed two hypothesis tests for testing the total mediation effect (TME) and the component-wise mediation effects (CME), respectively. We demonstrated in simulation studies that the type I error rates were controlled in both tests despite model misspecification. The TME test was more powerful than the CME test when the two mediation effects are in the same direction, while the CME test was more powerful than the TME test when the two mediation effects are in opposite direction. Multiple mediator analysis had increased power to detect mediated trans-eQTLs, especially in large samples. In the HapMap3 data, we identified 11 mediated trans-eQTLs that were not detected by the single mediator analysis in the combined samples of African populations. Moreover, the mediated trans-eQTLs in the HapMap3 samples are more likely to be trait-associated SNPs. In terms of computation, although there is no limit in the number of mediators in our model, analysis takes more time when adding additional mediators. In the analysis of the HapMap3 samples, we included at most 5 cis-gene mediators. Majority of the trios we considered have one or two mediators.

Conclusions

Trans-eQTLs are more likely to associate with multiple cis-genes than randomly selected SNPs. Mediation analysis with multiple mediators improves power of identification of mediated trans-eQTLs, especially in large samples.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-2651-6) contains supplementary material, which is available to authorized users.

Transcriptome-wide analysis of filarial extract-primed human monocytes reveal changes in LPS-induced PTX3 expression levels

Filarial nematodes modulate immune responses in their host to enable their survival and mediate protective effects against autoimmunity and allergies. In this study, we examined the immunomodulatory capacity of extracts from the human pathogenic filaria Brugia malayi (BmA) on human monocyte responses in a transcriptome-wide manner to identify associated pathways and diseases. As previous transcriptome studies often observed quiescent responses of innate cells to filariae, the potential of BmA to alter LPS driven responses was investigated by analyzing >47.000 transcripts of monocytes from healthy male volunteers stimulated with BmA, Escherichia coli LPS or a sequential stimulation of both. In comparison to ~2200 differentially expressed genes in LPS-only stimulated monocytes, only a limited number of differentially expressed genes were identified upon BmA priming before LPS re-stimulation with only PTX3↓ reaching statistical significance after correcting for multiple testing. Nominal significant differences were reached for metallothioneins↑, MMP9↑, CXCL5/ENA-78↑, CXCL6/GCP-2↑, TNFRSF21↓, and CCL20/MIP3α↓ and were confirmed by qPCR or ELISA. Flow cytometric analysis of activation markers revealed a reduced LPS-induced expression of HLA-DR and CD86 on BmA-primed monocytes as well as a reduced apoptosis of BmA-stimulated monocytes. While our experimental design does not allow a stringent extrapolation of our results to the development of filarial pathology, several genes that were identified in BmA-primed monocytes had previously been associated with filarial pathology, supporting the need for further research.

Exploring the genetic basis of human population differences in DNA methylation and their causal impact on immune gene regulation

BACKGROUND

DNA methylation is influenced by both environmental and genetic factors and is increasingly thought to affect variation in complex traits and diseases. Yet, the extent of ancestry-related differences in DNA methylation, their genetic determinants, and their respective causal impact on immune gene regulation remain elusive.

RESULTS

We report extensive population differences in DNA methylation between 156 individuals of African and European descent, detected in primary monocytes that are used as a model of a major innate immunity cell type. Most of these differences (~ 70%) are driven by DNA sequence variants nearby CpG sites, which account for ~ 60% of the variance in DNA methylation. We also identify several master regulators of DNA methylation variation in trans, including a regulatory hub nearby the transcription factor-encoding CTCF gene, which contributes markedly to ancestry-related differences in DNA methylation. Furthermore, we establish that variation in DNA methylation is associated with varying gene expression levels following mostly, but not exclusively, a canonical model of negative associations, particularly in enhancer regions. Specifically, we find that DNA methylation highly correlates with transcriptional activity of 811 and 230 genes, at the basal state and upon immune stimulation, respectively. Finally, using a Bayesian approach, we estimate causal mediation effects of DNA methylation on gene expression in ~ 20% of the studied cases, indicating that DNA methylation can play an active role in immune gene regulation.

CONCLUSION

Using a system-level approach, our study reveals substantial ancestry-related differences in DNA methylation and provides evidence for their causal impact on immune gene regulation.

Primary biliary cholangitis: A tale of epigenetically-induced secretory failure?

Primary biliary cholangitis (PBC) is a chronic cholestatic liver disease associated with autoimmune-related destruction of small to medium size intrahepatic bile ducts. The aetiology of PBC is unknown and its pathogenesis remains obscure. Both genetic variants and environmental factors have been linked to increased PBC susceptibility, with other alterations known to cooperate in disease pathobiology. Increasing evidence indicates the presence of epigenetic abnormalities in PBC, particularly alterations of cholangiocellular microRNAs (miRNAs or miRs). This review highlights and discusses the most relevant epigenetic alterations found in patients with PBC, focusing on the role of miR-506 in the promotion of cholestasis and immune activation.

NF-κB/MAPK activation underlies ACVR1-mediated inflammation in human heterotopic ossification

BACKGROUND

Inflammation helps regulate normal growth and tissue repair. Although bone morphogenetic proteins (BMPs) and inflammation are known contributors to abnormal bone formation, how these pathways interact in ossification remains unclear.

METHODS

We examined this potential link in patients with fibrodysplasia ossificans progressiva (FOP), a genetic condition of progressive heterotopic ossification caused by activating mutations in the Activin A type I receptor (ACVR1/ALK2). FOP patients show exquisite sensitivity to trauma, suggesting that BMP pathway activation may alter immune responses. We studied primary blood, monocyte, and macrophage samples from control and FOP subjects using multiplex cytokine, gene expression, and protein analyses; examined CD14+ primary monocyte and macrophage responses to TLR ligands; and assayed BMP, TGF-β activated kinase 1 (TAK1), and NF-κB pathways.

RESULTS

FOP subjects at baseline without clinically evident heterotopic ossification showed increased serum IL-3, IL-7, IL-8, and IL-10. CD14+ primary monocytes treated with the TLR4 activator LPS showed increased CCL5, CCR7, and CXCL10; abnormal cytokine/chemokine secretion; and prolonged activation of the NF-κB pathway. FOP macrophages derived from primary monocytes also showed abnormal cytokine/chemokine secretion, increased TGF-β production, and p38MAPK activation. Surprisingly, SMAD phosphorylation was not significantly changed in the FOP monocytes/macrophages.

CONCLUSIONS

Abnormal ACVR1 activity causes a proinflammatory state via increased NF-κB and p38MAPK activity. Similar changes may contribute to other types of heterotopic ossification, such as in scleroderma and dermatomyositis; after trauma; or with recombinant BMP-induced bone fusion. Our findings suggest that chronic antiinflammatory treatment may be useful for heterotopic ossification.

Are toll-like receptors potential drug targets for atherosclerosis? Evidence from genetic studies to date

Low-density lipoprotein cholesterol lowering, most notably via statin therapy, has successfully reduced the burden of coronary artery disease (CAD) in recent decades. However, the residual risk remaining even after aggressive lipid lowering has renewed interest in alternative targets. Anti-inflammatory drugs are thought to have much potential in this context, but side effects associated with long-term use of conventional anti-inflammatories, such as NSAIDs and glucocorticoids, preclude their use as preventive agents for CAD. Evidence from epidemiological studies and murine models of atherosclerosis suggests that toll-like receptors (TLRs) may have utility as targets for more focused anti-inflammatories, but it remains unclear if this pathway is causally related to CAD in man. Here, we review recent insight into this question gained from genetic studies of cardiovascular risk and innate immune function, focussing on the potential of Mendelian randomisation approaches based on intracellular-signalling pathways to identify and prioritise targets for drug development.

Viewpoint: Toward the Genetic Architecture of Disease Severity in Inflammatory Bowel Diseases

Inflammatory bowel disease (IBD) is characterized by uneven disease courses with various clinical outcomes. A few prognostic markers of disease severity may help stratify patients and identify those who will benefit the most from early aggressive treatment. The concept of disease severity remains too broad and vague, mainly because the definition must embrace several disease mechanisms, mainly inflammation and fibrosis, with various rates of disease progression. The magnitude of inflammation is an obvious key driver of disease severity in IBD that ultimately influence disease behavior. Advances in the genetics underlying disease severity are currently emerging, but attempts to overlap the genetics of disease susceptibility and severity have until now been unsatisfactory, suggesting that the genetic architecture of disease severity may be distinct from the genetics of disease susceptibility. In this review, we report on the current knowledge on disease severity and on the main research venues to decipher the genetic architecture of disease severity.

Baicalin Inhibits Haemophilus Parasuis-Induced High-Mobility Group Box 1 Release during Inflammation

Haemophilus parasuis (H. parasuis) can cause Glässer’s disease in pigs. However, the molecular mechanism of the inflammation response induced by H. parasuis remains unclear. The high-mobility group box 1 (HMGB1) protein is related to the pathogenesis of various infectious pathogens, but little is known about whether H. parasuis can induce the release of HMGB1 in piglet peripheral blood monocytes. Baicalin displays important anti-inflammatory and anti-microbial activities. In the present study, we investigated whether H. parasuis can trigger the secretion of HMGB1 in piglet peripheral blood monocytes and the anti-inflammatory effect of baicalin on the production of HMGB1 in peripheral blood monocytes induced by H. parasuis during the inflammation response. In addition, host cell responses stimulated by H. parasuis were determined with RNA-Seq. The RNA-Seq results showed that H. parasuis infection provokes the expression of cytokines and the activation of numerous pathways. In addition, baicalin significantly reduced the release of HMGB1 in peripheral blood monocytes induced by H. parasuis. Taken together, our study showed that H. parasuis can induce the release of HMGB1 and baicalin can inhibit HMGB1 secretion in an H. parasuis-induced peripheral blood monocytes model, which may provide a new strategy for preventing the inflammatory disorders induced by H. parasuis.

Genetic Variation at Chromosome 2q13 and Its Potential Influence on Endometriosis Susceptibility Through Effects on the IL-1 Family

Endometriosis is characterized by the growth of epithelial and stromal cells outside the uterine cavity. It has a complex etiology and affects ∼10% of reproductive age women. It is accompanied by a chronic inflammatory response with substantial evidence to indicate genetic susceptibility. The causal genes and their pathways leading to endometriosis, however, are still unknown. Recently, genomewide association studies on endometriosis identified 14 genomic risk loci in women of European and Japanese ancestry. It is becoming increasingly clear that these risk regions are intergenic and thus contribute to disease susceptibility through regulatory mechanisms, most likely mediated through regulation of genes within a restricted distance from the risk variants. One endometriosis risk locus has been detected at chromosome 2q13 within an inflammatory-rich region of gene transcripts and thus may play a role in the inflammation component of the disease. We carried out detailed analysis of the genomic region 250 kb on either side of sentinel SNP rs10167914 and identified 21 transcripts which contained 6 interleukin (IL)-1 family genes, 3 previously reported coding genes that have a relationship to inflammation, 4 novel coding, or pseudogenes, and 8 noncoding RNA transcripts. Through an extensive literature search, we examined the roles these genes and their resultant proteins play in endometriosis pathogenesis. The results suggest alteration in the expression the IL-1 family transcripts either alone or as a complex milieu could have a significant influence on endometriosis and should be prioritized for future study on the implications of inflammation on endometriotic lesions.

Diverse mechanisms of IRF5 action in inflammatory responses

Interferon regulatory factor 5 (IRF5) is a key signal-dependent transcription factor in myeloid cells. Its expression is induced by granulocyte-macrophage colony stimulating factor and interferon-gamma. IRF5 protein is further activated in response to stimulation, translocating to the nucleus where it mediates inflammatory responses. IRF5 is capable of both the up-regulation of pro-inflammatory genes and repressing anti-inflammatory mediators, thus polarising macrophages to a pro-inflammatory phenotype. We discuss IRF5 interactions with a wide range of transcriptional regulators that give rise to its diverse effects at the level of chromatin.

Shared genetic effects on chromatin and gene expression indicate a role for enhancer priming in immune response

Regulatory variants are often context specific, modulating gene expression in a subset of possible cellular states. Although these genetic effects can play important roles in disease, the molecular mechanisms underlying context specificity are poorly understood. Here, we identified shared quantitative trait loci (QTLs) for chromatin accessibility and gene expression in human macrophages exposed to IFNγ, Salmonella and IFNγ plus Salmonella. We observed that ~60% of stimulus-specific expression QTLs with a detectable effect on chromatin altered the chromatin accessibility in naive cells, thus suggesting that they perturb enhancer priming. Such variants probably influence binding of cell-type-specific transcription factors, such as PU.1, which can then indirectly alter the binding of stimulus-specific transcription factors, such as NF-κB or STAT2. Thus, although chromatin accessibility assays are powerful for fine-mapping causal regulatory variants, detecting their downstream effects on gene expression will be challenging, requiring profiling of large numbers of stimulated cellular states and time points.

Context-specific effects of genetic variants associated with autoimmune disease

Autoimmune diseases such as rheumatoid arthritis and coeliac disease are typical examples of complex genetic diseases caused by a combination of genetic and non-genetic risk factors. Insight into the genetic risk factors (single nucleotide polymorphisms (SNPs)) has increased since genome-wide association studies (GWAS) became possible in 2007 and, for individual diseases, SNPs can now explain some 15-50% of genetic risk. GWAS have also shown that some 50% of the genetic risk factors for individual autoimmune diseases overlap between different diseases. Thus, shared risk factors may converge to pathways that, when perturbed by genetic variation, predispose to autoimmunity in general. This raises the question of what determines disease specificity, and suggests that identical risk factors may have different effects in various autoimmune diseases. Addressing this question requires translation of genetic risk factors to causal genes and then to molecular and cellular pathways. Since >90% of the genetic risk factors are found in the non-coding part of the genome (i.e. outside the exons of protein-coding genes) and can have an impact on gene regulation, there is an urgent need to better understand the non-coding part of the genome. Here, we will outline the methods being used to unravel the gene regulatory networks perturbed in autoimmune diseases and the importance of doing this in the relevant cell types. We will highlight findings in coeliac disease, which manifests in the small intestine, to demonstrate how cell type and disease context can impact on the consequences of genetic risk factors.

A functional genomics predictive network model identifies regulators of inflammatory bowel disease

A major challenge in inflammatory bowel disease (IBD) is the integration of diverse IBD data sets to construct predictive models of IBD. We present a predictive model of the immune component of IBD that informs causal relationships among loci previously linked to IBD through genome-wide association studies (GWAS) using functional and regulatory annotations that relate to the cells, tissues, and pathophysiology of IBD. Our model consists of individual networks constructed using molecular data generated from intestinal samples isolated from three populations of patients with IBD at different stages of disease. We performed key driver analysis to identify genes predicted to modulate network regulatory states associated with IBD, prioritizing and prospectively validating 12 of the top key drivers experimentally. This validated key driver set not only introduces new regulators of processes central to IBD but also provides the integrated circuits of genetic, molecular, and clinical traits that can be directly queried to interrogate and refine the regulatory framework defining IBD.

Genetic regulatory effects modified by immune activation contribute to autoimmune disease associations

The immune system plays a major role in human health and disease, and understanding genetic causes of interindividual variability of immune responses is vital. Here, we isolate monocytes from 134 genotyped individuals, stimulate these cells with three defined microbe-associated molecular patterns (LPS, MDP, and 5'-ppp-dsRNA), and profile the transcriptomes at three time points. Mapping expression quantitative trait loci (eQTL), we identify 417 response eQTLs (reQTLs) with varying effects between conditions. We characterize the dynamics of genetic regulation on early and late immune response and observe an enrichment of reQTLs in distal cis-regulatory elements. In addition, reQTLs are enriched for recent positive selection with an evolutionary trend towards enhanced immune response. Finally, we uncover reQTL effects in multiple GWAS loci and show a stronger enrichment for response than constant eQTLs in GWAS signals of several autoimmune diseases. This demonstrates the importance of infectious stimuli in modifying genetic predisposition to disease.Insight into the genetic influence on the immune response is important for the understanding of interindividual variability in human pathologies. Here, the authors generate transcriptome data from human blood monocytes stimulated with various immune stimuli and provide a time-resolved response eQTL map.

Genetic Variability as a Regulator of TLR4 and NOD Signaling in Response to Bacterial Driven DNA Damage Response (DDR) and Inflammation: Focus on the Gastrointestinal (GI) Tract

The fundamental role of human Toll-like receptors (TLRs) and NOD-like receptors (NLRs), the two most studied pathogen recognition receptors (PRRs), is the protection against pathogens and excessive tissue injury. Recent evidence supports the association between TLR/NLR gene mutations and susceptibility to inflammatory, autoimmune, and malignant diseases. PRRs also interfere with several cellular processes, such as cell growth, apoptosis, cell proliferation, differentiation, autophagy, angiogenesis, cell motility and migration, and DNA repair mechanisms. We briefly review the impact of TLR4 and NOD1/NOD2 and their genetic variability in the process of inflammation, tumorigenesis and DNA repair, focusing in the gastrointestinal tract. We also review the available data on new therapeutic strategies utilizing TLR/NLR agonists and antagonists for cancer, allergic diseases, viral infections and vaccine development against both infectious diseases and cancer.

Toll-like receptor 4 in glial inflammatory responses to air pollution in vitro and in vivo

Background

Exposure to traffic-related air pollution (TRAP) is associated with accelerated cognitive aging and higher dementia risk in human populations. Rodent brains respond to TRAP with activation of astrocytes and microglia, increased inflammatory cytokines, and neurite atrophy. A role for Toll-like receptor 4 (TLR4) was suggested in mouse TLR4-knockouts, which had attenuated lung macrophage responses to air pollution.

Methods

To further analyze these mechanisms, we examined mixed glial cultures (astrocytes and microglia) for RNA responses to nanoscale particulate matter (nPM; diameter <0.2 μm), a well-characterized nanoscale particulate matter subfraction of TRAP collected from a local freeway (Morgan et al. Environ Health Perspect 2011; 119,1003–1009, 2011). The nPM was compared with responses to the endotoxin lipopolysaccharide (LPS), a classic TLR4 ligand, using Affymetrix whole genome microarray in rats. Expression patterns were analyzed by significance analysis of microarrays (SAM) for fold change and by weighted gene co-expression network analysis (WGCNA) to identify modules of shared responses between nPM and LPS. Finally, we examined TLR4 activation in hippocampal tissue from mice chronically exposed to nPM.

Results

SAM and WGCNA analyses showed strong activation of TLR4 and NF-κB by both nPM and LPS. TLR4 siRNA attenuated TNFα and other inflammatory responses to nPM in vitro, via the MyD88-dependent pathway. In vivo, mice chronically exposed to nPM showed increased TLR4, MyD88, TNFα, and TNFR2 RNA, and decreased NF-κB and TRAF6 RNA TLR4 and NF-κB responses in the hippocampus.

Conclusions

These results show TLR4 activation is integral in brain inflammatory responses to air pollution, and warrant further study of TLR4 in accelerated cognitive aging by air pollution.

Electronic supplementary material

The online version of this article (doi:10.1186/s12974-017-0858-x) contains supplementary material, which is available to authorized users.

Expression QTLs Mapping and Analysis: A Bayesian Perspective

The aim of expression Quantitative Trait Locus (eQTL) mapping is the identification of DNA sequence variants that explain variation in gene expression. Given the recent yield of trait-associated genetic variants identified by large-scale genome-wide association analyses (GWAS), eQTL mapping has become a useful tool to understand the functional context where these variants operate and eventually narrow down functional gene targets for disease. Despite its extensive application to complex (polygenic) traits and disease, the majority of eQTL studies still rely on univariate data modeling strategies, i.e., testing for association of all transcript-marker pairs. However these "one at-a-time" strategies are (1) unable to control the number of false-positives when an intricate Linkage Disequilibrium structure is present and (2) are often underpowered to detect the full spectrum of trans-acting regulatory effects. Here we present our viewpoint on the most recent advances on eQTL mapping approaches, with a focus on Bayesian methodology. We review the advantages of the Bayesian approach over frequentist methods and provide an empirical example of polygenic eQTL mapping to illustrate the different properties of frequentist and Bayesian methods. Finally, we discuss how multivariate eQTL mapping approaches have distinctive features with respect to detection of polygenic effects, accuracy, and interpretability of the results.

Concerted Genetic Function in Blood Traits

The hematopoietic system plays a major role in human health. Two studies by Astle et al. and Chen et al. published in this issue of Cell use genome-wide association and functional genomics approaches to provide deep insights into the role of genetic variants in hematological traits. We discuss these discoveries and future strategies toward completing our understanding of the genetic basis for variation in human traits.

Genome-wide transcriptome induced by nickel in human monocytes

Nickel-containing alloys are frequently used in the biomedical field, although, owing to corrosive processes metal ion leaching is inevitable. Due to nickel ion (Ni(2+)) leaching several adverse effects are described in the literature. However, only a few studies evaluated the genetic profile of Ni(2+) in human cells which is of great importance since nickel-induced effects differ between humans and mice as a result of species-specific receptor variability. Thus, we investigated gene expression induced by Ni(2+)in human monocytes using a transcriptome-wide approach determining new target genes implicated in nickel-induced pathologies. Monocytes were isolated from healthy volunteers of Central European origin using stringent inclusion criteria. Cells were challenged with different Ni(2+) concentrations. Array-based gene expression analysis was performed comprising more than 47,000 transcripts followed by pathway analyses. Transcriptional data were validated by protein and cell surface markers. Ni(2+) significantly influenced the expression of 1385 transcripts in a dose-dependent manner. Apart from known targets (CCL20↑, PTGS2↑, MTs↑, SLCs↑), we identified new candidates implicated in Ni(2+)-elicited processes (various microRNAs↑, INSIG1↑, NAMPT↑, MS4A6A↓, DHRS9↓). Several of these transcripts correspond to immunity, inflammation and were shown to be involved in cellular reactions related to hypersensitivity, cancer, colitis, and encephalitis. Moreover, 459 canonical pathways/signaling, 500 pathologies and 2687 upstream regulators were detected. Protein results validated our findings. To our knowledge, the present systematic transcriptome-wide expression study is the first which explored Ni(2+)-elicited cell responses in human primary monocytes identifying new target genes, pathways and upstream regulators of relevance to diagnostic and therapeutic strategies.

STATEMENT OF SIGNIFICANCE

Nickel is widely applied in the biomedical field, although several adverse effects are documented in the literature due to nickel ion (Ni(2+)) leaching. In humans, allergic reactions like contact dermatitis are the most common adverse effect to Ni(2+), whereas serious concerns relate to possible systemic and carcinogenic activities. Using a systematic genome-wide transcriptional approach in human primary monocytes unveil new target genes, pathways and upstream regulators implicated in nickel-elicited immune response which are of significance to diagnostic and therapeutic strategies. This approach provides new information of how host-derived immune response contributes to the interaction with antigens and supports the interplay between metal ions and systemic diseases.

Gene-based analysis of regulatory variants identifies 4 putative novel asthma risk genes related to nucleotide synthesis and signaling

BACKGROUND

Hundreds of genetic variants are thought to contribute to variation in asthma risk by modulating gene expression. Methods that increase the power of genome-wide association studies (GWASs) to identify risk-associated variants are needed.

OBJECTIVE

We sought to develop a method that aggregates the evidence for association with disease risk across expression quantitative trait loci (eQTLs) of a gene and use this approach to identify asthma risk genes.

METHODS

We developed a gene-based test and software package called EUGENE that (1) is applicable to GWAS summary statistics; (2) considers both cis- and trans-eQTLs; (3) incorporates eQTLs identified in different tissues; and (4) uses simulations to account for multiple testing. We applied this approach to 2 published asthma GWASs (combined n = 46,044) and used mouse studies to provide initial functional insights into 2 genes with novel genetic associations.

RESULTS

We tested the association between asthma and 17,190 genes that were found to have cis- and/or trans-eQTLs across 16 published eQTL studies. At an empirical FDR of 5%, 48 genes were associated with asthma risk. Of these, for 37, the association was driven by eQTLs located in established risk loci for allergic disease, including 6 genes not previously implicated in disease cause (eg, LIMS1, TINF2, and SAFB). The remaining 11 significant genes represent potential novel genetic associations with asthma. The association with 4 of these replicated in an independent GWAS: B4GALT3, USMG5, P2RY13, and P2RY14, which are genes involved in nucleotide synthesis or nucleotide-dependent cell activation. In mouse studies, P2ry13 and P2ry14-purinergic receptors activated by adenosine 5-diphosphate and UDP-sugars, respectively-were upregulated after allergen challenge, notably in airway epithelial cells, eosinophils, and neutrophils. Intranasal exposure with receptor agonists induced the release of IL-33 and subsequent eosinophil infiltration into the lungs.

CONCLUSION

We identified novel associations between asthma and eQTLs for 4 genes related to nucleotide synthesis/signaling and demonstrated the power of gene-based analyses of GWASs.

Targeted genomic analysis reveals widespread autoimmune disease association with regulatory variants in the TNF superfamily cytokine signalling network

BACKGROUND

Tumour necrosis factor (TNF) superfamily cytokines and their receptors regulate diverse immune system functions through a common set of signalling pathways. Genetic variants in and expression of individual TNF superfamily cytokines, receptors and signalling proteins have been associated with autoimmune and inflammatory diseases, but their interconnected biology has been largely unexplored.

METHODS

We took a hypothesis-driven approach using available genome-wide datasets to identify genetic variants regulating gene expression in the TNF superfamily cytokine signalling network and the association of these variants with autoimmune and autoinflammatory disease. Using paired gene expression and genetic data, we identified genetic variants associated with gene expression, expression quantitative trait loci (eQTLs), in four peripheral blood cell subsets. We then examined whether eQTLs were dependent on gene expression level or the presence of active enhancer chromatin marks. Using these eQTLs as genetic markers of the TNF superfamily signalling network, we performed targeted gene set association analysis in eight autoimmune and autoinflammatory disease genome-wide association studies.

RESULTS

Comparison of TNF superfamily network gene expression and regulatory variants across four leucocyte subsets revealed patterns that differed between cell types. eQTLs for genes in this network were not dependent on absolute gene expression levels and were not enriched for chromatin marks of active enhancers. By examining autoimmune disease risk variants among our eQTLs, we found that risk alleles can be associated with either increased or decreased expression of co-stimulatory TNF superfamily cytokines, receptors or downstream signalling molecules. Gene set disease association analysis revealed that eQTLs for genes in the TNF superfamily pathway were associated with six of the eight autoimmune and autoinflammatory diseases examined, demonstrating associations beyond single genome-wide significant hits.

CONCLUSIONS

This systematic analysis of the influence of regulatory genetic variants in the TNF superfamily network reveals widespread and diverse roles for these cytokines in susceptibility to a number of immune-mediated diseases.

Interferon regulatory factor 5 in human autoimmunity and murine models of autoimmune disease

Interferon regulatory factor 5 (IRF5) has been demonstrated as a key transcription factor of the immune system, playing important roles in modulating inflammatory immune responses in numerous cell types including dendritic cells, macrophages, and B cells. As well as driving the expression of type I interferon in antiviral responses, IRF5 is also crucial for driving macrophages toward a proinflammatory phenotype by regulating cytokine and chemokine expression and modulating B-cell maturity and antibody production. This review highlights the functional importance of IRF5 in a disease setting, by discussing polymorphic mutations at the human Irf5 locus that lead to susceptibility to systemic lupus erythematosus, rheumatoid arthritis, and inflammatory bowel disease. In concordance with this, we also discuss lessons in IRF5 functionality learned from murine in vivo models of autoimmune disease and inflammation and hypothesize that modulation of IRF5 activity and expression could provide potential therapeutic benefits in the clinic.

Genome-wide transcriptome induced by Porphyromonas gingivalis LPS supports the notion of host-derived periodontal destruction and its association with systemic diseases

Chronic periodontitis (CP) is a prevalent pathogen-associated inflammatory disorder characterized by the destruction of tooth-supporting tissues, and linked to several systemic diseases. Both the periodontopathogen Porphyromonas gingivalis (Pg), and the genetically determined host immune response, are hypothesized to play a crucial role in this association. To identify new target genes for CP and its associated systemic diseases, we investigated the transcriptome induced by Pg in human monocytes using a genome-wide approach. Monocytes were isolated from healthy male volunteers of European origin and challenged with the Pg virulence factor LPS. Array-based gene expression analysis comprising >47,000 transcripts was performed followed by pathway analyses. Transcriptional data were validated by protein and cell surface markers. LPS Pg challenge led to the significant induction of 902 transcripts. Besides known periodontitis-associated targets, several new candidates were identified (CCL23↑, INDO↑, GBP 1/4↑, CFB↑, ISG20↑, MIR155HG↑, DHRS9↓). Moreover, various transcripts correspond to the host immune response, and have been linked to cancer, atherosclerosis and arthritis, thus highlighting the systemic impact of CP. Protein data of immunological markers validated our results. The present findings expand understanding of Pg elicited immune responses, and indicate new target genes and pathways of relevance to diagnostic and therapeutic strategies.

Transcriptional and Linkage Analyses Identify Loci that Mediate the Differential Macrophage Response to Inflammatory Stimuli and Infection

Macrophages display flexible activation states that range between pro-inflammatory (classical activation) and anti-inflammatory (alternative activation). These macrophage polarization states contribute to a variety of organismal phenotypes such as tissue remodeling and susceptibility to infectious and inflammatory diseases. Several macrophage- or immune-related genes have been shown to modulate infectious and inflammatory disease pathogenesis. However, the potential role that differences in macrophage activation phenotypes play in modulating differences in susceptibility to infectious and inflammatory disease is just emerging. We integrated transcriptional profiling and linkage analyses to determine the genetic basis for the differential murine macrophage response to inflammatory stimuli and to infection with the obligate intracellular parasite Toxoplasma gondii. We show that specific transcriptional programs, defined by distinct genomic loci, modulate macrophage activation phenotypes. In addition, we show that the difference between AJ and C57BL/6J macrophages in controlling Toxoplasma growth after stimulation with interferon gamma and tumor necrosis factor alpha mapped to chromosome 3, proximal to the Guanylate binding protein (Gbp) locus that is known to modulate the murine macrophage response to Toxoplasma. Using an shRNA-knockdown strategy, we show that the transcript levels of an RNA helicase, Ddx1, regulates strain differences in the amount of nitric oxide produced by macrophage after stimulation with interferon gamma and tumor necrosis factor. Our results provide a template for discovering candidate genes that modulate macrophage-mediated complex traits.

Macrophages provide a first line of defense against invading pathogens and play an important role in the initiation and resolution of immune responses. When in contact with pathogens or immune factors, such as cytokines, macrophages assume activation states that range between pro-inflammatory (classical activation) and anti-inflammatory (alternative activation). Even though it is known that macrophages from different individuals are biased towards one of the various activation states, the genetic factors that define individual differences in macrophage activation are not fully understood. Additionally, although macrophages are important in infectious disease pathogenesis, how individual differences in macrophage activation contribute to individual differences in susceptibility to infectious disease is just emerging. We used macrophages from genetically segregating mice to show that discrete transcriptional programs, which are modulated by specific genomic regions, modulate differences in macrophage activation. Murine macrophages differences in controlling Toxoplasma growth mapped to chromosome 3, proximal to the Guanylate binding protein (Gbp) locus that is known to modulate the murine macrophage response to Toxoplasma. Using a shRNA-mediated knockdown approach, we show that the DEAD box polypeptide 1 (Ddx1) modulates nitric oxide production in macrophages stimulated with interferon gamma and tumor necrosis factor. These findings are a step towards the identification of genes that regulate macrophage phenotypes and disease outcome.

Genome-wide analysis of the genetic regulation of gene expression in human neutrophils

Neutrophils are an abundant immune cell type involved in both antimicrobial defence and autoimmunity. The regulation of their gene expression, however, is still largely unknown. Here we report an eQTL study on isolated neutrophils from 114 healthy individuals of Chinese ethnicity, identifying 21,210 eQTLs on 832 unique genes. Unsupervised clustering analysis of these eQTLs confirms their role in inflammatory responses and immunological diseases but also indicates strong involvement in dermatological pathologies. One of the strongest eQTL identified (rs2058660) is also the tagSNP of a linkage block reported to affect leprosy and Crohn's disease in opposite directions. In a functional study, we can link the C allele with low expression of the β-chain of IL18-receptor (IL18RAP). In neutrophils, this results in a reduced responsiveness to IL-18, detected both on the RNA and protein level. Thus, the polymorphic regulation of human neutrophils can impact beneficial as well as pathological inflammatory responses.

Neutrophils are abundant immune cells important for antimicrobial defence and in autoimmunity. Here, by mapping expression quantitative trait loci (eQTL) in neutrophils of Chinese ethnicity from Singapore, Andiappan et al. provide a resource for understanding immune-related trait associated genetic variants.