An iPSC-derived small intestine-on-chip with self-organizing epithelial, mesenchymal, and neural cells.

Human induced pluripotent stem cell (hiPSC)-derived intestinal organoids are valuable tools for researching developmental biology and personalized therapies, but their closed topology and relative immature state limit applications. Here, we use organ-on-chip technology to develop a hiPSC-derived intestinal barrier with apical and basolateral access in a more physiological in vitro microenvironment. To replicate growth factor gradients along the crypt-villus axis, we locally expose the cells to expansion and differentiation media. In these conditions, intestinal epithelial cells self-organize into villus-like folds with physiological barrier integrity, and myofibroblasts and neurons emerge and form a subepithelial tissue in the bottom channel. The growth factor gradients efficiently balance dividing and mature cell types and induce an intestinal epithelial composition, including absorptive and secretory lineages, resembling the composition of the human small intestine. This well-characterized hiPSC-derived intestine-on-chip system can facilitate personalized studies on physiological processes and therapy development in the human small intestine.

Genetic mapping across autoimmune diseases reveals shared associations and mechanisms.

Autoimmune and inflammatory diseases are polygenic disorders of the immune system. Many genomic loci harbor risk alleles for several diseases, but the limited resolution of genetic mapping prevents determining whether the same allele is responsible, indicating a shared underlying mechanism. Here, using a collection of 129,058 cases and controls across 6 diseases, we show that ~40% of overlapping associations are due to the same allele. We improve fine-mapping resolution for shared alleles twofold by combining cases and controls across diseases, allowing us to identify more expression quantitative trait loci driven by the shared alleles. The patterns indicate widespread sharing of pathogenic mechanisms but not a single global autoimmune mechanism. Our approach can be applied to any set of traits and is particularly valuable as sample collections become depleted.

Phage display sequencing reveals that genetic, environmental, and intrinsic factors influence variation of human antibody epitope repertoire.

Phage-displayed immunoprecipitation sequencing (PhIP-seq) has enabled high-throughput profiling of human antibody repertoires. However, a comprehensive overview of environmental and genetic determinants shaping human adaptive immunity is lacking. In this study, we investigated the effects of genetic, environmental, and intrinsic factors on the variation in human antibody repertoires. We characterized serological antibody repertoires against 344,000 peptides using PhIP-seq libraries from a wide range of microbial and environmental antigens in 1,443 participants from a population cohort. We detected individual-specificity, temporal consistency, and co-housing similarities in antibody repertoires. Genetic analyses showed the involvement of the HLA, IGHV, and FUT2 gene regions in antibody-bound peptide reactivity. Furthermore, we uncovered associations between phenotypic factors (including age, cell counts, sex, smoking behavior, and allergies, among others) and particular antibody-bound peptides. Our results indicate that human antibody epitope repertoires are shaped by both genetics and environmental exposures and highlight specific signatures of distinct phenotypes and genotypes.

iPSC-derived organ-on-a-chip models for personalized human genetics and pharmacogenomics studies.

Genome-wide association studies (GWAS) have now correlated hundreds of genetic variants with complex genetic diseases and drug efficacy. Functional characterization of these factors remains challenging, particularly because of the lack of human model systems. Molecular and nanotechnological advances, in particular the ability to generate patient-specific PSC lines, differentiate them into diverse cell types, and seed and combine them on microfluidic chips, have led to the establishment of organ-on-a-chip (OoC) platforms that recapitulate organ biology. OoC technology thus provides unique personalized platforms for studying the effects of host genetics and environmental factors on organ physiology. In this review we describe the technology and provide examples of how OoCs may be used for disease modeling and pharmacogenetic research.

Gut mucosa dissociation protocols influence cell type proportions and single-cell gene expression levels.

Single-cell RNA sequencing (scRNA-seq) has revolutionized the study of the cellular landscape of organs. Most single-cell protocols require fresh material, which limits sample size per experiment, and consequently, introduces batch effects. This is especially true for samples acquired through complex medical procedures, such as intestinal mucosal biopsies. Moreover, the tissue dissociation procedure required for obtaining single cells is a major source of noise; different dissociation procedures applied to different compartments of the tissue induce artificial gene expression differences between cell subsets. To overcome these challenges, we have developed a one-step dissociation protocol and demonstrated its use on cryopreserved gut mucosal biopsies. Using flow cytometry and scRNA-seq analysis, we compared this one-step dissociation protocol with the current gold standard, two-step collagenase digestion, and an adaptation of a recently published alternative, three-step cold-active Bacillus licheniformus protease digestion. Both cell viability and cell type composition were comparable between the one-step and two-step collagenase dissociation, with the former being more time-efficient. The cold protease digestion resulted in equal cell viability, but better preserves the epithelial cell types. Consequently, to analyze the rarer cell types, such as glial cells, larger total biopsy cell numbers are required as input material. The multi-step protocols affected cell types spanning multiple compartments differently. In summary, we show that cryopreserved gut mucosal biopsies can be used to overcome the logistical challenges and batch effects in large scRNA-seq studies. Furthermore, we demonstrate that using cryopreserved biopsies digested using a one-step collagenase protocol enables large-scale scRNA-seq, FACS, organoid generation and intraepithelial lymphocyte expansion.

Single-Cell RNA Sequencing of Peripheral Blood Mononuclear Cells From Pediatric Coeliac Disease Patients Suggests Potential Pre-Seroconversion Markers.

Celiac Disease (CeD) is a complex immune disorder involving villous atrophy in the small intestine that is triggered by gluten intake. Current CeD diagnosis is based on late-stage pathophysiological parameters such as detection of specific antibodies in blood and histochemical detection of villus atrophy and lymphocyte infiltration in intestinal biopsies. To date, no early onset biomarkers are available that would help prevent widespread villous atrophy and severe symptoms and co-morbidities. To search for novel CeD biomarkers, we used single-cell RNA sequencing (scRNAseq) to investigate PBMC samples from 11 children before and after seroconversion for CeD and 10 control individuals matched for age, sex and HLA-genotype. We generated scRNAseq profiles of 9559 cells and identified the expected major cellular lineages. Cell proportions remained stable across the different timepoints and health conditions, but we observed differences in gene expression profiles in specific cell types when comparing patient samples before and after disease development and comparing patients with controls. Based on the time when transcripts were differentially expressed, we could classify the deregulated genes as biomarkers for active CeD or as potential pre-diagnostic markers. Pathway analysis showed that active CeD biomarkers display a transcriptional profile associated with antigen activation in CD4+ T cells, whereas NK cells express a subset of biomarker genes even before CeD diagnosis. Intersection of biomarker genes with CeD-associated genetic risk loci pinpointed genetic factors that might play a role in CeD onset. Investigation of potential cellular interaction pathways of PBMC cell subpopulations highlighted the importance of TNF pathways in CeD. Altogether, our results pinpoint genes and pathways that are altered prior to and during CeD onset, thereby identifying novel potential biomarkers for CeD diagnosis in blood.

Circulating miRNAs as Potential Biomarkers for Celiac Disease Development.

Background & Aims: Celiac disease (CeD), an immune-mediated disease with enteropathy triggered by gluten, affects ~1% of the general European population. Currently, there are no biomarkers to predict CeD development. MicroRNAs (miRNAs) are short RNAs involved in post-transcriptional gene regulation, and certain disease- and stage-specific miRNA profiles have been found previously. We aimed to investigate whether circulating miRNAs can predict the development of CeD. Methods: Using next-generation miRNA-sequencing, we determined miRNAs in >200 serum samples from 53 participants of the PreventCD study, of whom 33 developed CeD during follow-up. Following study inclusion at 3 months of age, samples were drawn at predefined ages, diagnosis (first anti-transglutaminase antibody (TGA) positivity or diagnostic biopsy) and after the start of a gluten-free diet (GFD). This allowed identification of circulating miRNAs that are deregulated before TGA positivity. For validation of the biomarkers for CeD and GFD response, two additional cohorts were included in subsequent meta-analyses. Additionally, miRNAs were measured in duodenal biopsies in a case-control cohort. Results: 53 circulating miRNAs were increased (27) or decreased (26) in CeD versus controls. We assessed specific trends in these individual miRNAs in the PreventCD cohort by grouping the pre-diagnostic samples of the CeD patients (all had negative TGA) by how close to seroconversion (first sample positive TGA) the samples were taken. 8/53 miRNAs differed significantly between controls and samples taken <1 year before TGA positivity: miR-21-3p, miR-374a-5p, 144-3p, miR-500a-3p, miR-486-3p let-7d-3p, let-7e-5p and miR-3605-3p. 6/26 downregulated miRNAs reconstituted upon GFD, including miR-150-5p/-3p, whereas no upregulated miRNAs were downregulated upon GFD. 15/53 biomarker candidates also differed between CeD biopsies and controls, with a concordant direction, indicating that these circulating miRNAs might originate from the intestine. Conclusions: We identified 53 circulating miRNAs that are potential early biomarkers for CeD, of which several can be detected more than a year before TGA positivity and some start to normalize upon GFD.

A Combined mRNA- and miRNA-Sequencing Approach Reveals miRNAs as Potential Regulators of the Small Intestinal Transcriptome in Celiac Disease.

Celiac disease (CeD) is triggered by gluten and results in inflammation and villous atrophy of the small intestine. We aimed to explore the role of miRNA-mediated deregulation of the transcriptome in CeD. Duodenal biopsies of CeD patients (n = 33) and control subjects (n = 10) were available for miRNA-sequencing, with RNA-sequencing also available for controls (n = 5) and CeD (n = 6). Differential expression analysis was performed to select CeD-associated miRNAs and genes. MiRNA‒target transcript pairs selected from public databases that also displayed a strong negative expression correlation in the current dataset (R < -0.7) were used to construct a CeD miRNA‒target transcript interaction network. The network includes 2030 miRNA‒target transcript interactions, including 423 experimentally validated pairs. Pathway analysis found that interactions are involved in immune-related pathways (e.g., interferon signaling) and metabolic pathways (e.g., lipid metabolism). The network includes 13 genes previously prioritized to be causally deregulated by CeD-associated genomic variants, including STAT1. CeD-associated miRNAs might play a role in promoting inflammation and decreasing lipid metabolism in the small intestine, thereby contributing unbalanced cell turnover in the intestinal crypt. Some CeD-associated miRNAs deregulate genes that are also affected by genomic CeD-risk variants, adding an additional layer of complexity to the deregulated transcriptome in CeD.

Inflammatory Protein Profiles in Plasma of Candidaemia Patients and the Contribution of Host Genetics to Their Variability.

Circulatory inflammatory proteins play a significant role in anti-Candida host immune defence. However, little is known about the genetic variation that contributes to the variability of inflammatory responses in response to C. albicans. To systematically characterize inflammatory responses in Candida infection, we profiled 91 circulatory inflammatory proteins in peripheral blood mononuclear cells (PBMCs) stimulated with C. albicans yeast isolated from 378 individuals of European origin from the 500 Functional Genomics (500FG) cohort of the Human Functional Genomics Project (HFGP) and Lifelines Deep cohort. To identify the genetic factors that determine variation in inflammatory protein responses, we correlated genome-wide single nucleotide polymorphism (SNP) genotypes with protein abundance (protein quantitative trait loci, pQTLs) produced by the Candida-stimulated PBMCs. Furthermore, we investigated whether differences in survival of candidaemia patients can be explained by modulating levels of inflammatory proteins. We identified five genome-wide significant pQTLs that modulate IL-8, MCP-2, MMP-1, and CCL3 in response to C. albicans. In addition, our genetic analysis suggested that GADD45G from rs10114707 locus that reached genome-wide significance could be a potential core gene that regulates a cytokine network upon Candida infection. Last but not least, we observed that a trans-pQTL marked from SNP rs7651677 at chromosome 3 that influences urokinase plasminogen activator (uPA) is strongly associated with patient survival (Psurvival = 3.52 x 10-5, OR 3). Overall, our genetic analysis showed that genetic variation determines the abundance of circulatory proteins in response to Candida infection.

Long Non-Coding RNAs Involved in Progression of Non-Alcoholic Fatty Liver Disease to Steatohepatitis.

Non-alcoholic fatty liver disease (NAFLD) is the most prevalent chronic liver disease and is characterized by different stages varying from benign fat accumulation to non-alcoholic steatohepatitis (NASH) that may progress to cirrhosis and liver cancer. In recent years, a regulatory role of long non-coding RNAs (lncRNAs) in NAFLD has emerged. Therefore, we aimed to characterize the still poorly understood lncRNA contribution to disease progression. Transcriptome analysis in 60 human liver samples with various degrees of NAFLD/NASH was combined with a functional genomics experiment in an in vitro model where we exposed HepG2 cells to free fatty acids (FFA) to induce steatosis, then stimulated them with tumor necrosis factor alpha (TNFα) to mimic inflammation. Bioinformatics analyses provided a functional prediction of novel lncRNAs. We further functionally characterized the involvement of one novel lncRNA in the nuclear-factor-kappa B (NF-κB) signaling pathway by its silencing in Hepatoma G2 (HepG2) cells. We identified 730 protein-coding genes and 18 lncRNAs that responded to FFA/TNFα and associated with human NASH phenotypes with consistent effect direction, with most being linked to inflammation. One novel intergenic lncRNA, designated lncTNF, was 20-fold up-regulated upon TNFα stimulation in HepG2 cells and positively correlated with lobular inflammation in human liver samples. Silencing lncTNF in HepG2 cells reduced NF-κB activity and suppressed expression of the NF-κB target genes A20 and NFKBIA. The lncTNF we identified in the NF-κB signaling pathway may represent a novel target for controlling liver inflammation.

Integration of metabolomics, genomics, and immune phenotypes reveals the causal roles of metabolites in disease.

BACKGROUND: Recent studies highlight the role of metabolites in immune diseases, but it remains unknown how much of this effect is driven by genetic and non-genetic host factors. RESULT: We systematically investigate circulating metabolites in a cohort of 500 healthy subjects (500FG) in whom immune function and activity are deeply measured and whose genetics are profiled. Our data reveal that several major metabolic pathways, including the alanine/glutamate pathway and the arachidonic acid pathway, have a strong impact on cytokine production in response to ex vivo stimulation. We also examine the genetic regulation of metabolites associated with immune phenotypes through genome-wide association analysis and identify 29 significant loci, including eight novel independent loci. Of these, one locus (rs174584-FADS2) associated with arachidonic acid metabolism is causally associated with Crohn's disease, suggesting it is a potential therapeutic target. CONCLUSION: This study provides a comprehensive map of the integration between the blood metabolome and immune phenotypes, reveals novel genetic factors that regulate blood metabolite concentrations, and proposes an integrative approach for identifying new disease treatment targets.

Potential impact of celiac disease genetic risk factors on T cell receptor signaling in gluten-specific CD4+ T cells.

Celiac disease is an auto-immune disease in which an immune response to dietary gluten leads to inflammation and subsequent atrophy of small intestinal villi, causing severe bowel discomfort and malabsorption of nutrients. The major instigating factor for the immune response in celiac disease is the activation of gluten-specific CD4+ T cells expressing T cell receptors that recognize gluten peptides presented in the context of HLA-DQ2 and DQ8. Here we provide an in-depth characterization of 28 gluten-specific T cell clones. We assess their transcriptional and epigenetic response to T cell receptor stimulation and link this to genetic factors associated with celiac disease. Gluten-specific T cells have a distinct transcriptional profile that mostly resembles that of Th1 cells but also express cytokines characteristic of other types of T-helper cells. This transcriptional response appears not to be regulated by changes in chromatin state, but rather by early upregulation of transcription factors and non-coding RNAs that likely orchestrate the subsequent activation of genes that play a role in immune pathways. Finally, integration of chromatin and transcription factor binding profiles suggest that genes activated by T cell receptor stimulation of gluten­specific T cells may be impacted by genetic variation at several genetic loci associated with celiac disease.

Non-classical clinical presentation at diagnosis by male celiac disease patients of older age.

BACKGROUND: . In a biopsy-proven adult celiac disease (CeD) cohort from the Netherlands, male patients were diagnosed with CeD at significantly older ages than female patients. OBJECTIVES: To identify which factors contribute to diagnosis later in life and whether diagnostic delay influences improvement of symptoms after starting a gluten-free diet (GFD). METHODS: . We performed a questionnaire study in 211 CeD patients (67:144, male:female) with median age at diagnosis of 41.8 years (interquartile range: 25-58) and at least Marsh 2 histology. RESULTS: . Classical symptoms (diarrhea, fatigue, abdominal pain and/or weight loss) were more frequent in women than men, but sex was not significantly associated with age at diagnosis. In a multivariate analysis, a non-classical presentation (without any classical symptoms) and a negative family history of CeD were significant predictors of older age at diagnosis (coefficients of 8 and 12 years, respectively). A delay of >3 years between first symptom and diagnosis was associated with slower improvement of symptoms after start of GFD, but not with sex, presentation of classical symptoms or age at diagnosis. CONCLUSION: . Non-classical CeD presentation is more prevalent in men and is associated with a diagnosis of CeD later in life. Recognizing CeD sooner after onset of symptoms is important because a long diagnostic delay is associated with a slower improvement of symptoms after starting a GFD.

Molecular Biomarkers for Celiac Disease: Past, Present and Future.

Celiac disease (CeD) is a complex immune-mediated disorder that is triggered by dietary gluten in genetically predisposed individuals. CeD is characterized by inflammation and villous atrophy of the small intestine, which can lead to gastrointestinal complaints, malnutrition, and malignancies. Currently, diagnosis of CeD relies on serology (antibodies against transglutaminase and endomysium) and small-intestinal biopsies. Since small-intestinal biopsies require invasive upper-endoscopy, and serology cannot predict CeD in an early stage or be used for monitoring disease after initiation of a gluten-free diet, the search for non-invasive biomarkers is ongoing. Here, we summarize current and up-and-coming non-invasive biomarkers that may be able to predict, diagnose, and monitor the progression of CeD. We further discuss how current and emerging techniques, such as (single-cell) transcriptomics and genomics, can be used to uncover the pathophysiology of CeD and identify non-invasive biomarkers.

Deconvolution of bulk blood eQTL effects into immune cell subpopulations.

BACKGROUND: Expression quantitative trait loci (eQTL) studies are used to interpret the function of disease-associated genetic risk factors. To date, most eQTL analyses have been conducted in bulk tissues, such as whole blood and tissue biopsies, which are likely to mask the cell type-context of the eQTL regulatory effects. Although this context can be investigated by generating transcriptional profiles from purified cell subpopulations, current methods to do this are labor-intensive and expensive. We introduce a new method, Decon2, as a framework for estimating cell proportions using expression profiles from bulk blood samples (Decon-cell) followed by deconvolution of cell type eQTLs (Decon-eQTL). RESULTS: The estimated cell proportions from Decon-cell agree with experimental measurements across cohorts (R ≥ 0.77). Using Decon-cell, we could predict the proportions of 34 circulating cell types for 3194 samples from a population-based cohort. Next, we identified 16,362 whole-blood eQTLs and deconvoluted cell type interaction (CTi) eQTLs using the predicted cell proportions from Decon-cell. CTi eQTLs show excellent allelic directional concordance with eQTL (≥ 96-100%) and chromatin mark QTL (≥87-92%) studies that used either purified cell subpopulations or single-cell RNA-seq, outperforming the conventional interaction effect. CONCLUSIONS: Decon2 provides a method to detect cell type interaction effects from bulk blood eQTLs that is useful for pinpointing the most relevant cell type for a given complex disease. Decon2 is available as an R package and Java application (https://github.com/molgenis/systemsgenetics/tree/master/Decon2) and as a web tool (www.molgenis.org/deconvolution).

Tissue alarmins and adaptive cytokine induce dynamic and distinct transcriptional responses in tissue-resident intraepithelial cytotoxic T lymphocytes.

The respective effects of tissue alarmins interleukin (IL)-15 and interferon beta (IFNß), and IL-21 produced by T cells on the reprogramming of cytotoxic T lymphocytes (CTLs) that cause tissue destruction in celiac disease is poorly understood. Transcriptomic and epigenetic profiling of primary intestinal CTLs showed massive and distinct temporal transcriptional changes in response to tissue alarmins, while the impact of IL-21 was limited. Only anti-viral pathways were induced in response to all the three stimuli, albeit with differences in dynamics and strength. Moreover, changes in gene expression were primarily independent of changes in H3K27ac, suggesting that other regulatory mechanisms drive the robust transcriptional response. Finally, we found that IL-15/IFNß/IL-21 transcriptional signatures could be linked to transcriptional alterations in risk loci for complex immune diseases. Together these results provide new insights into molecular mechanisms that fuel the activation of CTLs under conditions that emulate the inflammatory environment in patients with autoimmune diseases.

Systematic Prioritization of Candidate Genes in Disease Loci Identifies TRAFD1 as a Master Regulator of IFNγ Signaling in Celiac Disease.

Celiac disease (CeD) is a complex T cell-mediated enteropathy induced by gluten. Although genome-wide association studies have identified numerous genomic regions associated with CeD, it is difficult to accurately pinpoint which genes in these loci are most likely to cause CeD. We used four different in silico approaches-Mendelian randomization inverse variance weighting, COLOC, LD overlap, and DEPICT-to integrate information gathered from a large transcriptomics dataset. This identified 118 prioritized genes across 50 CeD-associated regions. Co-expression and pathway analysis of these genes indicated an association with adaptive and innate cytokine signaling and T cell activation pathways. Fifty-one of these genes are targets of known drug compounds or likely druggable genes, suggesting that our methods can be used to pinpoint potential therapeutic targets. In addition, we detected 172 gene combinations that were affected by our CeD-prioritized genes in trans. Notably, 41 of these trans-mediated genes appear to be under control of one master regulator, TRAF-type zinc finger domain containing 1 (TRAFD1), and were found to be involved in interferon (IFN)γ signaling and MHC I antigen processing/presentation. Finally, we performed in vitro experiments in a human monocytic cell line that validated the role of TRAFD1 as an immune regulator acting in trans. Our strategy confirmed the role of adaptive immunity in CeD and revealed a genetic link between CeD and IFNγ signaling as well as with MHC I antigen processing, both major players of immune activation and CeD pathogenesis.

Leukocyte-Released Mediators in Response to Both Bacterial and Fungal Infections Trigger IFN Pathways, Independent of IL-1 and TNF-α, in Endothelial Cells.

In sepsis, dysregulated immune responses to infections cause damage to the host. Previous studies have attempted to capture pathogen-induced leukocyte responses. However, the impact of mediators released after pathogen-leukocyte interaction on endothelial cells, and how endothelial cell responses vary depending on the pathogen-type is lacking. Here, we comprehensively characterized the transcriptomic responses of human leukocytes and endothelial cells to Gram negative-bacteria, Gram positive-bacteria, and fungi. We showed that whole pathogen lysates induced strong activation of leukocytes but not endothelial cells. Interestingly, the common response of leukocytes to various pathogens converges on endothelial activation. By exposing endothelial cells to leukocyte-released mediators, we observed a strong activation of endothelial cells at both transcription and protein levels. By adding IL-1RA and TNF-α antibody in leukocyte-released mediators before exposing to endothelial cells, we identified specific roles for IL-1 and TNF-α in driving the most, but not all, endothelial activation. We also showed for the first time, activation of interferon response by endothelial cells in response to leukocyte-released mediators, independently from IL-1 and TNF-α pathways. Our study therefore, not only provides pathogen-dependent transcriptional changes in leukocytes and endothelial cells during infections, but also reveals a role for IFN, together with IL1 and TNFα signaling, in mediating leukocyte-endothelial interaction in infections.

Intestinal Barrier Function in Gluten-Related Disorders.

Gluten-related disorders include distinct disease entities, namely celiac disease, wheat-associated allergy and non-celiac gluten/wheat sensitivity. Despite having in common the contact of the gastrointestinal mucosa with components of wheat and other cereals as a causative factor, these clinical entities have distinct pathophysiological pathways. In celiac disease, a T-cell mediate immune reaction triggered by gluten ingestion is central in the pathogenesis of the enteropathy, while wheat allergy develops as a rapid immunoglobulin E- or non-immunoglobulin E-mediated immune response. In non-celiac wheat sensitivity, classical adaptive immune responses are not involved. Instead, recent research has revealed that an innate immune response to a yet-to-be-defined antigen, as well as the gut microbiota, are pivotal in the development in this disorder. Although impairment of the epithelial barrier has been described in all three clinical conditions, its role as a potential pathogenetic co-factor, specifically in celiac disease and non-celiac wheat sensitivity, is still a matter of investigation. This article gives a short overview of the mucosal barrier of the small intestine, summarizes the aspects of barrier dysfunction observed in all three gluten-related disorders and reviews literature data in favor of a primary involvement of the epithelial barrier in the development of celiac disease and non-celiac wheat sensitivity.

A Genome-Wide Functional Genomics Approach Identifies Susceptibility Pathways to Fungal Bloodstream Infection in Humans.

BACKGROUND: Candidemia, one of the most common causes of fungal bloodstream infection, leads to mortality rates up to 40% in affected patients. Understanding genetic mechanisms for differential susceptibility to candidemia may aid in designing host-directed therapies. METHODS: We performed the first genome-wide association study on candidemia, and we integrated these data with variants that affect cytokines in different cellular systems stimulated with Candida albicans. RESULTS: We observed strong association between candidemia and a variant, rs8028958, that significantly affects the expression levels of PLA2G4B in blood. We found that up to 35% of the susceptibility loci affect in vitro cytokine production in response to Candida. Furthermore, potential causal genes located within these loci are enriched for lipid and arachidonic acid metabolism. Using an independent cohort, we also showed that the numbers of risk alleles at these loci are negatively correlated with reactive oxygen species and interleukin-6 levels in response to Candida. Finally, there was a significant correlation between susceptibility and allelic scores based on 16 independent candidemia-associated single-nucleotide polymorphisms that affect monocyte-derived cytokines, but not with T cell-derived cytokines. CONCLUSIONS: Our results prioritize the disturbed lipid homeostasis and oxidative stress as potential mechanisms that affect monocyte-derived cytokines to influence susceptibility to candidemia.