Prostaglandins and calprotectin are genetically and functionally linked to the Inflammatory Bowel Diseases.

BACKGROUND: Genome wide association studies (GWAS) have identified and validated more than 200 genomic loci associated with the inflammatory bowel disease (IBD), although for most the causal gene remains unknown. Given the importance of myeloid cells in IBD pathogenesis, the current study aimed to uncover the role of genes within IBD genetic loci that are endogenously expressed in this cell lineage. METHODS: The open reading frames (ORF) of 42 genes from IBD-associated loci were expressed via lentiviral transfer in the THP-1 model of human monocytes and the impact of each of these on the cell's transcriptome was analyzed using a RNA sequencing-based approach. We used a combination of genetic and pharmacologic approaches to validate our findings in the THP-1 line with further validation in human induced pluripotent stem cell (hiPSC)-derived-monocytes. RESULTS: This functional genomics screen provided evidence that genes in four IBD GWAS loci (PTGIR, ZBTB40, SLC39A11 and NFKB1) are involved in controlling S100A8 and S100A9 gene expression, which encode the two subunits of calprotectin (CP). We demonstrated that increasing PTGIR expression and/or stimulating PTGIR signaling resulted in increased CP expression in THP-1. This was further validated in hiPSC-derived monocytes. Conversely, knocking-down PTGIR endogenous expression and/or inhibiting PTGIR signaling led to decreased CP expression. These analyses were extended to the known IBD gene PTGER4, whereby its specific agonist also led to increased CP expression. Furthermore, we demonstrated that the PTGIR and PTGER4 mediated control of CP expression was dependent on signaling via adenylate cyclase and STAT3. Finally, we demonstrated that LPS-mediated increases in CP expression could be potentiated by agonists of PTGIR and PTGER4, and diminished by their antagonists. CONCLUSION: Our results support a causal role for the PTGIR, PTGER4, ZBTB40, SLC39A11 and NFKB1 genes in IBD, with all five genes regulating the expression of CP in myeloid cells, as well as potential roles for the prostacyclin/prostaglandin biogenesis and signaling pathways in IBD susceptibility and pathogenesis.

Transethnic analysis of the human leukocyte antigen region for ulcerative colitis reveals not only shared but also ethnicity-specific disease associations.

Inflammatory bowel disease (IBD) is a chronic inflammatory disease of the gut. Genetic association studies have identified the highly variable human leukocyte antigen (HLA) region as the strongest susceptibility locus for IBD and specifically DRB1*01:03 as a determining factor for ulcerative colitis (UC). However, for most of the association signal such as delineation could not be made because of tight structures of linkage disequilibrium within the HLA. The aim of this study was therefore to further characterize the HLA signal using a transethnic approach. We performed a comprehensive fine mapping of single HLA alleles in UC in a cohort of 9272 individuals with African American, East Asian, Puerto Rican, Indian and Iranian descent and 40 691 previously analyzed Caucasians, additionally analyzing whole HLA haplotypes. We computationally characterized the binding of associated HLA alleles to human self-peptides and analyzed the physicochemical properties of the HLA proteins and predicted self-peptidomes. Highlighting alleles of the HLA-DRB1*15 group and their correlated HLA-DQ-DR haplotypes, we not only identified consistent associations (regarding effects directions/magnitudes) across different ethnicities but also identified population-specific signals (regarding differences in allele frequencies). We observed that DRB1*01:03 is mostly present in individuals of Western European descent and hardly present in non-Caucasian individuals. We found peptides predicted to bind to risk HLA alleles to be rich in positively charged amino acids. We conclude that the HLA plays an important role for UC susceptibility across different ethnicities. This research further implicates specific features of peptides that are predicted to bind risk and protective HLA proteins.

Fine-mapping inflammatory bowel disease loci to single-variant resolution.

Inflammatory bowel diseases are chronic gastrointestinal inflammatory disorders that affect millions of people worldwide. Genome-wide association studies have identified 200 inflammatory bowel disease-associated loci, but few have been conclusively resolved to specific functional variants. Here we report fine-mapping of 94 inflammatory bowel disease loci using high-density genotyping in 67,852 individuals. We pinpoint 18 associations to a single causal variant with greater than 95% certainty, and an additional 27 associations to a single variant with greater than 50% certainty. These 45 variants are significantly enriched for protein-coding changes (n = 13), direct disruption of transcription-factor binding sites (n = 3), and tissue-specific epigenetic marks (n = 10), with the last category showing enrichment in specific immune cells among associations stronger in Crohn's disease and in gut mucosa among associations stronger in ulcerative colitis. The results of this study suggest that high-resolution fine-mapping in large samples can convert many discoveries from genome-wide association studies into statistically convincing causal variants, providing a powerful substrate for experimental elucidation of disease mechanisms.

Construction and benchmarking of a multi-ethnic reference panel for the imputation of HLA class I and II alleles.

Genotype imputation of the human leukocyte antigen (HLA) region is a cost-effective means to infer classical HLA alleles from inexpensive and dense SNP array data. In the research setting, imputation helps avoid costs for wet lab-based HLA typing and thus renders association analyses of the HLA in large cohorts feasible. Yet, most HLA imputation reference panels target Caucasian ethnicities and multi-ethnic panels are scarce. We compiled a high-quality multi-ethnic reference panel based on genotypes measured with Illumina's Immunochip genotyping array and HLA types established using a high-resolution next generation sequencing approach. Our reference panel includes more than 1,300 samples from Germany, Malta, China, India, Iran, Japan and Korea and samples of African American ancestry for all classical HLA class I and II alleles including HLA-DRB3/4/5. Applying extensive cross-validation, we benchmarked the imputation using the HLA imputation tool HIBAG, our multi-ethnic reference and an independent, previously published data set compiled of subpopulations of the 1000 Genomes project. We achieved average imputation accuracies higher than 0.924 for the commonly studied HLA-A, -B, -C, -DQB1 and -DRB1 genes across all ethnicities. We investigated allele-specific imputation challenges in regard to geographic origin of the samples using sensitivity and specificity measurements as well as allele frequencies and identified HLA alleles that are challenging to impute for each of the populations separately. In conclusion, our new multi-ethnic reference data set allows for high resolution HLA imputation of genotypes at all classical HLA class I and II genes including the HLA-DRB3/4/5 loci based on diverse ancestry populations.

Correction: The Influence of Age and Sex on Genetic Associations with Adult Body Size and Shape: A Large-Scale Genome-Wide Interaction Study.

[This corrects the article DOI: 10.1371/journal.pgen.1005378.].

Comprehensive and Reproducible Untargeted Lipidomic Workflow Using LC-QTOF Validated for Human Plasma Analysis.

The goal of this work was to develop a label-free, comprehensive, and reproducible high-resolution liquid chromatography-mass spectrometry (LC-MS)-based untargeted lipidomic workflow using a single instrument, which could be applied to biomarker discovery in both basic and clinical studies. For this, we have (i) optimized lipid extraction and elution to enhance coverage of polar and nonpolar lipids as well as resolution of their isomers, (ii) ensured MS signal reproducibility and linearity, and (iii) developed a bioinformatic pipeline to correct remaining biases. Workflow validation is reported for 48 replicates of a single human plasma sample: 1124 reproducible LC-MS signals were extracted (median signal intensity RSD = 10%), 50% of which are redundant due to adducts, dimers, in-source fragmentation, contaminations, or positive and negative ion duplicates. From the resulting 578 unique compounds, 428 lipids were identified by MS/MS, including acyl chain composition, of which 394 had RSD < 30% inside their linear intensity range, thereby enabling robust semiquantitation. MS signal intensity spanned 4 orders of magnitude, covering 16 lipid subclasses. Finally, the power of our workflow is illustrated by a proof-of-concept study in which 100 samples from healthy human subjects were analyzed and the data set was investigated using three different statistical testing strategies in order to compare their capacity in identifying the impact of sex and age on circulating lipids.

Genetic Factors Interact With Tobacco Smoke to Modify Risk for Inflammatory Bowel Disease in Humans and Mice.

BACKGROUND & AIMS: The role of tobacco smoke in the etiology of inflammatory bowel disease (IBD) is unclear. We investigated interactions between genes and smoking (gene-smoking interactions) that affect risk for Crohn's disease (CD) and ulcerative colitis (UC) in a case-only study of patients and in mouse models of IBD. METHODS: We used 55 Immunochip-wide datasets that included 19,735 IBD cases (10,856 CD cases and 8879 UC cases) of known smoking status. We performed 3 meta-analyses each for CD, UC, and IBD (CD and UC combined), comparing data for never vs ever smokers, never vs current smokers, and never vs former smokers. We studied the effects of exposure to cigarette smoke in Il10-/- and Nod2-/- mice, as well as in Balb/c mice without disruption of these genes (wild-type mice). Mice were exposed to the smoke of 5 cigarettes per day, 5 days a week, for 8 weeks, in a ventilated smoking chamber, or ambient air (controls). Intestines were collected and analyzed histologically and by reverse transcription polymerase chain reaction. RESULTS: We identified 64 single nucleotide polymorphisms (SNPs) for which the association between the SNP and IBD were modified by smoking behavior (meta-analysis Wald test P < 5.0 × 10-5; heterogeneity Cochrane Q test P > .05). Twenty of these variants were located within the HLA region at 6p21. Analysis of classical HLA alleles (imputed from SNP genotypes) revealed an interaction with smoking. We replicated the interaction of a variant in NOD2 with current smoking in relation to the risk for CD (frameshift variant fs1007insC; rs5743293). We identified 2 variants in the same genomic region (rs2270368 and rs17221417) that interact with smoking in relation to CD risk. Approximately 45% of the SNPs that interact with smoking were in close vicinity (≤1 Mb) to SNPs previously associated with IBD; many were located near or within genes that regulate mucosal barrier function and immune tolerance. Smoking modified the disease risk of some variants in opposite directions for CD vs UC. Exposure of Interleukin 10 (il10)-deficient mice to cigarette smoke accelerated development of colitis and increased expression of interferon gamma in the small intestine compared to wild-type mice exposed to smoke. NOD2-deficient mice exposed to cigarette smoke developed ileitis, characterized by increased expression of interferon gamma, compared to wild-type mice exposed to smoke. CONCLUSIONS: In an analysis of 55 Immunochip-wide datasets, we identified 64 SNPs whose association with risk for IBD is modified by tobacco smoking. Gene-smoking interactions were confirmed in mice with disruption of Il10 and Nod2-variants of these genes have been associated with risk for IBD. Our findings from mice and humans revealed that the effects of smoking on risk for IBD depend on genetic variants.

Host-microbe interactions have shaped the genetic architecture of inflammatory bowel disease.

Crohn's disease and ulcerative colitis, the two common forms of inflammatory bowel disease (IBD), affect over 2.5 million people of European ancestry, with rising prevalence in other populations. Genome-wide association studies and subsequent meta-analyses of these two diseases as separate phenotypes have implicated previously unsuspected mechanisms, such as autophagy, in their pathogenesis and showed that some IBD loci are shared with other inflammatory diseases. Here we expand on the knowledge of relevant pathways by undertaking a meta-analysis of Crohn's disease and ulcerative colitis genome-wide association scans, followed by extensive validation of significant findings, with a combined total of more than 75,000 cases and controls. We identify 71 new associations, for a total of 163 IBD loci, that meet genome-wide significance thresholds. Most loci contribute to both phenotypes, and both directional (consistently favouring one allele over the course of human history) and balancing (favouring the retention of both alleles within populations) selection effects are evident. Many IBD loci are also implicated in other immune-mediated disorders, most notably with ankylosing spondylitis and psoriasis. We also observe considerable overlap between susceptibility loci for IBD and mycobacterial infection. Gene co-expression network analysis emphasizes this relationship, with pathways shared between host responses to mycobacteria and those predisposing to IBD.

The Influence of Age and Sex on Genetic Associations with Adult Body Size and Shape: A Large-Scale Genome-Wide Interaction Study.

Genome-wide association studies (GWAS) have identified more than 100 genetic variants contributing to BMI, a measure of body size, or waist-to-hip ratio (adjusted for BMI, WHRadjBMI), a measure of body shape. Body size and shape change as people grow older and these changes differ substantially between men and women. To systematically screen for age- and/or sex-specific effects of genetic variants on BMI and WHRadjBMI, we performed meta-analyses of 114 studies (up to 320,485 individuals of European descent) with genome-wide chip and/or Metabochip data by the Genetic Investigation of Anthropometric Traits (GIANT) Consortium. Each study tested the association of up to ~2.8M SNPs with BMI and WHRadjBMI in four strata (men ≤50y, men >50y, women ≤50y, women >50y) and summary statistics were combined in stratum-specific meta-analyses. We then screened for variants that showed age-specific effects (G x AGE), sex-specific effects (G x SEX) or age-specific effects that differed between men and women (G x AGE x SEX). For BMI, we identified 15 loci (11 previously established for main effects, four novel) that showed significant (FDR<5%) age-specific effects, of which 11 had larger effects in younger (<50y) than in older adults (≥50y). No sex-dependent effects were identified for BMI. For WHRadjBMI, we identified 44 loci (27 previously established for main effects, 17 novel) with sex-specific effects, of which 28 showed larger effects in women than in men, five showed larger effects in men than in women, and 11 showed opposite effects between sexes. No age-dependent effects were identified for WHRadjBMI. This is the first genome-wide interaction meta-analysis to report convincing evidence of age-dependent genetic effects on BMI. In addition, we confirm the sex-specificity of genetic effects on WHRadjBMI. These results may provide further insights into the biology that underlies weight change with age or the sexually dimorphism of body shape.

Inherited determinants of Crohn's disease and ulcerative colitis phenotypes: a genetic association study.

BACKGROUND: Crohn's disease and ulcerative colitis are the two major forms of inflammatory bowel disease; treatment strategies have historically been determined by this binary categorisation. Genetic studies have identified 163 susceptibility loci for inflammatory bowel disease, mostly shared between Crohn's disease and ulcerative colitis. We undertook the largest genotype association study, to date, in widely used clinical subphenotypes of inflammatory bowel disease with the goal of further understanding the biological relations between diseases. METHODS: This study included patients from 49 centres in 16 countries in Europe, North America, and Australasia. We applied the Montreal classification system of inflammatory bowel disease subphenotypes to 34,819 patients (19,713 with Crohn's disease, 14,683 with ulcerative colitis) genotyped on the Immunochip array. We tested for genotype-phenotype associations across 156,154 genetic variants. We generated genetic risk scores by combining information from all known inflammatory bowel disease associations to summarise the total load of genetic risk for a particular phenotype. We used these risk scores to test the hypothesis that colonic Crohn's disease, ileal Crohn's disease, and ulcerative colitis are all genetically distinct from each other, and to attempt to identify patients with a mismatch between clinical diagnosis and genetic risk profile. FINDINGS: After quality control, the primary analysis included 29,838 patients (16,902 with Crohn's disease, 12,597 with ulcerative colitis). Three loci (NOD2, MHC, and MST1 3p21) were associated with subphenotypes of inflammatory bowel disease, mainly disease location (essentially fixed over time; median follow-up of 10·5 years). Little or no genetic association with disease behaviour (which changed dramatically over time) remained after conditioning on disease location and age at onset. The genetic risk score representing all known risk alleles for inflammatory bowel disease showed strong association with disease subphenotype (p=1·65 × 10(-78)), even after exclusion of NOD2, MHC, and 3p21 (p=9·23 × 10(-18)). Predictive models based on the genetic risk score strongly distinguished colonic from ileal Crohn's disease. Our genetic risk score could also identify a small number of patients with discrepant genetic risk profiles who were significantly more likely to have a revised diagnosis after follow-up (p=6·8 × 10(-4)). INTERPRETATION: Our data support a continuum of disorders within inflammatory bowel disease, much better explained by three groups (ileal Crohn's disease, colonic Crohn's disease, and ulcerative colitis) than by Crohn's disease and ulcerative colitis as currently defined. Disease location is an intrinsic aspect of a patient's disease, in part genetically determined, and the major driver to changes in disease behaviour over time. FUNDING: International Inflammatory Bowel Disease Genetics Consortium members funding sources (see Acknowledgments for full list).

Deep resequencing of GWAS loci identifies rare variants in CARD9, IL23R and RNF186 that are associated with ulcerative colitis.

Genome-wide association studies and follow-up meta-analyses in Crohn's disease (CD) and ulcerative colitis (UC) have recently identified 163 disease-associated loci that meet genome-wide significance for these two inflammatory bowel diseases (IBD). These discoveries have already had a tremendous impact on our understanding of the genetic architecture of these diseases and have directed functional studies that have revealed some of the biological functions that are important to IBD (e.g. autophagy). Nonetheless, these loci can only explain a small proportion of disease variance (~14% in CD and 7.5% in UC), suggesting that not only are additional loci to be found but that the known loci may contain high effect rare risk variants that have gone undetected by GWAS. To test this, we have used a targeted sequencing approach in 200 UC cases and 150 healthy controls (HC), all of French Canadian descent, to study 55 genes in regions associated with UC. We performed follow-up genotyping of 42 rare non-synonymous variants in independent case-control cohorts (totaling 14,435 UC cases and 20,204 HC). Our results confirmed significant association to rare non-synonymous coding variants in both IL23R and CARD9, previously identified from sequencing of CD loci, as well as identified a novel association in RNF186. With the exception of CARD9 (OR = 0.39), the rare non-synonymous variants identified were of moderate effect (OR = 1.49 for RNF186 and OR = 0.79 for IL23R). RNF186 encodes a protein with a RING domain having predicted E3 ubiquitin-protein ligase activity and two transmembrane domains. Importantly, the disease-coding variant is located in the ubiquitin ligase domain. Finally, our results suggest that rare variants in genes identified by genome-wide association in UC are unlikely to contribute significantly to the overall variance for the disease. Rather, these are expected to help focus functional studies of the corresponding disease loci.

A meta-analysis of genome-wide association scans identifies IL18RAP, PTPN2, TAGAP, and PUS10 as shared risk loci for Crohn's disease and celiac disease.

Crohn's disease (CD) and celiac disease (CelD) are chronic intestinal inflammatory diseases, involving genetic and environmental factors in their pathogenesis. The two diseases can co-occur within families, and studies suggest that CelD patients have a higher risk to develop CD than the general population. These observations suggest that CD and CelD may share common genetic risk loci. Two such shared loci, IL18RAP and PTPN2, have already been identified independently in these two diseases. The aim of our study was to explicitly identify shared risk loci for these diseases by combining results from genome-wide association study (GWAS) datasets of CD and CelD. Specifically, GWAS results from CelD (768 cases, 1,422 controls) and CD (3,230 cases, 4,829 controls) were combined in a meta-analysis. Nine independent regions had nominal association p-value <1.0 x 10⁻5 in this meta-analysis and showed evidence of association to the individual diseases in the original scans (p-value < 1 x 10⁻² in CelD and < 1 x 10⁻³ in CD). These include the two previously reported shared loci, IL18RAP and PTPN2, with p-values of 3.37 x 10⁻8 and 6.39 x 10⁻9, respectively, in the meta-analysis. The other seven had not been reported as shared loci and thus were tested in additional CelD (3,149 cases and 4,714 controls) and CD (1,835 cases and 1,669 controls) cohorts. Two of these loci, TAGAP and PUS10, showed significant evidence of replication (Bonferroni corrected p-values <0.0071) in the combined CelD and CD replication cohorts and were firmly established as shared risk loci of genome-wide significance, with overall combined p-values of 1.55 x 10⁻¹° and 1.38 x 10⁻¹¹ respectively. Through a meta-analysis of GWAS data from CD and CelD, we have identified four shared risk loci: PTPN2, IL18RAP, TAGAP, and PUS10. The combined analysis of the two datasets provided the power, lacking in the individual GWAS for single diseases, to detect shared loci with a relatively small effect.

Identification of a sudden cardiac death susceptibility locus at 2q24.2 through genome-wide association in European ancestry individuals.

Sudden cardiac death (SCD) continues to be one of the leading causes of mortality worldwide, with an annual incidence estimated at 250,000-300,000 in the United States and with the vast majority occurring in the setting of coronary disease. We performed a genome-wide association meta-analysis in 1,283 SCD cases and >20,000 control individuals of European ancestry from 5 studies, with follow-up genotyping in up to 3,119 SCD cases and 11,146 controls from 11 European ancestry studies, and identify the BAZ2B locus as associated with SCD (P = 1.8×10(-10)). The risk allele, while ancestral, has a frequency of ~1.4%, suggesting strong negative selection and increases risk for SCD by 1.92-fold per allele (95% CI 1.57-2.34). We also tested the role of 49 SNPs previously implicated in modulating electrocardiographic traits (QRS, QT, and RR intervals). Consistent with epidemiological studies showing increased risk of SCD with prolonged QRS/QT intervals, the interval-prolonging alleles are in aggregate associated with increased risk for SCD (P = 0.006).

Identifying transcript-level differential expression in primary human immune cells.

BACKGROUND: Multipotential hematopoietic stem cells differentiate into a wide variety of immune cells with a diversity of functions, including the ability to respond to a variety of stimuli. Importantly, numerous studies have demonstrated the importance of gene transcription in defining cell identity and functions. While these studies have primarily been performed at the level of the gene, it is known that key immune genes such as CD44 and CD45 generate multiple different transcripts that are differentially expressed across different immune cells, and that encode proteins with different sequences and functions. Prior genomic surveys have shown that the mechanisms for generating diversity in expressed transcripts (alternate splicing, alternate transcription start sites, etc.) are very active in immune cells, but have been lacking in terms of identifying genes with multiple transcripts, that are differentially expressed, and likely to affect cell functions. METHODS: We first identified the set of genes that had at least two transcripts expressed in our RNA sequencing dataset generated from purified populations of neutrophils, monocytes and five lymphocyte populations (B, NK, γδ T, CD4 + T and CD8 + T) from twelve healthy donors. Next, we developed a heuristic approach to identify genes where two or more transcripts have distinct expression patterns across lymphoid and/or myeloid populations. We then focused our annotation and interpretation on differentially expressed transcripts that affect the coding sequence. This process was repeated to identify transcripts that were differentially expressed between monocytes and populations of macrophages and LPS-stimulated macrophages derived from these monocytes in vitro. RESULTS: We found that over 55 % of genes had two or more expressed transcripts, with an average ∼3 transcripts per gene, and that 70 % of these had at least two of the transcripts that encoded proteins with different sequences. As expected, we identified a complex pattern of differential expression for multiple transcripts encoding the CD45 transmembrane protein, but we also found similar evidence for ten other genes (CD300A, FYB1, GPI, LITAF, PSMA1, PTMA, RPL32, SEPTIN9, SH3BP2, SH3KBP1) when comparing the expression patterns of transcripts within myeloid and lymphoid cells. We also identified five genes with differentially expressed transcripts associated with the transition from monocytes to macrophages (FNBP1, KLF6, and SEPTIN9) or between macrophages and LPS-stimulated macrophages (CD44, OAZ2, and SEPTIN9). For the most part, we found that the different transcripts of these genes are expected to impact specific biological functions, for example the different transcripts of SEPTIN9 likely regulate the cytoskeleton in immune cells via their interactions with actins filaments and microtubules. CONCLUSIONS: This analytic approach successfully identified multi-transcript genes that are differentially expressed across immune cells and could be applied to other transcriptomic data. DATA AVAILABILITY STATEMENT: Researchers can request access to the individual-level data from the current study by contacting the Montreal Heart Institute ethics committee at the following institutional email address: cer.icm@icm-mhi.org.

Sex-Dimorphic Analyses Identify Novel and Sex-Specific Genetic Associations in Inflammatory Bowel Disease.

BACKGROUND: Sex is an integral variable often overlooked in complex disease genetics. Differences between sexes have been reported in natural history, disease complications, and age of onset in inflammatory bowel disease (IBD). While association studies have identified >230 IBD loci, there have been a limited number of studies investigating sex differences underlying these genetic associations. METHODS: We report the first investigation of sex-dimorphic associations via meta-analysis of a sex-stratified association study (34 579 IBD cases, 39 125 controls). In addition, we performed chromosome (chr) X-specific analyses, considering models of X inactivation (XCI) and XCI escape. Demographic and clinical characteristics were also compared between sexes. RESULTS: We identified significant differences between sexes for disease location and perianal complication in Crohn's disease and disease extent in ulcerative colitis. We observed genome-wide-significant sex-dimorphic associations (P < 5 × 10-8) at loci not previously reported in large-scale IBD genetic studies, including at chr9q22, CARMIL1, and UBASH3A. We identified variants in known IBD loci, including in chr2p15 and within the major histocompatibility complex on chr6, exhibiting sex-specific patterns of association (P < 5 × 10-7 in one sex only). We identified 3 chrX associations with IBD, including a novel Crohn's disease susceptibility locus at Xp22. CONCLUSIONS: These analyses identified novel IBD loci, in addition to characterizing sex-specific patterns of associations underlying sex-dimorphic associations. By elucidating the role of sex in IBD genetics, our study will help enhance our understanding of the differences between the sexes in IBD biology and underscores a need to move beyond conventional sex-combined analyses to appreciate the genetic architecture of IBD more comprehensively.

Sex-dimorphic meta-analyses of sex-stratified case-control (n = 73 704) regression identified 3 novel inflammatory bowel disease loci reaching genome-wide significance and highlighted chromosome 2 and major histocompatibility complex variants exhibiting sex-specific association. In addition, a novel chromosome X Crohn's disease susceptibility locus was identified.

Serum Lipidomic Screen Identifies Key Metabolites, Pathways, and Disease Classifiers in Crohn's Disease.

BACKGROUND: There is an unmet medical need for biomarkers that capture host and environmental contributions in inflammatory bowel diseases (IBDs). This study aimed at testing the potential of circulating lipids as disease classifiers given their major roles in inflammation. METHODS: We applied a previously validated comprehensive high-resolution liquid chromatography-mass spectrometry-based untargeted lipidomic workflow covering 25 lipid subclasses to serum samples from 100 Crohn's disease (CD) patients and 100 matched control subjects. Findings were replicated and expanded in another 200 CD patients and 200 control subjects. Key metabolites were tested for associations with disease behavior and location, and classification models were built and validated. Their association with disease activity was tested using an independent cohort of 42 CD patients. RESULTS: We identified >70 metabolites with strong association (P < 1 × 10-4, q < 5 × 10-4) to CD. Highly performing classification models (area under the curve > 0.84-0.97) could be built with as few as 5 to 9 different metabolites, representing 6 major correlated lipid clusters. These classifiers included a phosphatidylethanolamine ether (O-16:0/20:4), a sphingomyelin (d18:1/21:0) and a cholesterol ester (14:1), a very long-chain dicarboxylic acid [28:1(OH)] and sitosterol sulfate. These classifiers and correlated lipids indicate a dysregulated metabolism in host cells, notably in peroxisomes, as well as dysbiosis, oxidative stress, compromised inflammation resolution, or intestinal membrane integrity. A subset of these were associated with disease behavior or location. CONCLUSIONS: Untargeted lipidomic analyses uncovered perturbations in the circulating human CD lipidome, likely resulting from multiple pathogenic mechanisms. Models using as few as 5 biomarkers had strong disease classifier characteristics, supporting their potential use in diagnosis or prognosis.

This study reports a comprehensive untargeted lipidomic analysis of 600 serum samples from patients with Crohn's disease and matched control subjects, identified and replicated ~70 metabolites associated with Crohn's disease, and developed highly performing classification models (area under the curve > 0.84-0.97) with as few as 5 metabolites.

IBD-associated G protein-coupled receptor 65 variant compromises signalling and impairs key functions involved in inflammation.

BACKGROUND AND AIMS: Inflammatory bowel diseases (IBD) result in chronic inflammation of the gastrointestinal tract. Genetic studies have shown that the GPR65 gene, as well as its missense coding variant, GPR65*Ile231Leu, is associated with IBD. We aimed to define the signalling and biological pathways downstream of GPR65 activation and evaluate the impact of GPR65*231Leu on these. METHODS: We used HEK 293 cells stably expressing GPR65 and deficient for either Gαs, Gαq/11 or Gα12/13, to define GPR65 signalling pathways, IBD patient biopsies and a panel of human tissues, primary immune cells and cell lines to determine biologic context, and genetic modulation of human THP-1-derived macrophages to examine the impact of GPR65 in bacterial phagocytosis and NLRP3 inflammasome activation. RESULTS: We confirmed that GPR65 signals via the Gαs pathway, leading to cAMP accumulation. GPR65 can also signal via the Gα12/13 pathway leading to formation of stress fibers, actin remodeling and RhoA activation; all impaired by the IBD-associated GPR65*231Leu allele. Gene expression profiling revealed greater expression of GPR65 in biopsies from inflamed compared to non-inflamed tissues from IBD patients or control individuals, potentially explained by infiltration of inflammatory immune cells. Decreased GPR65 expression in THP-1-derived macrophages leads to impaired bacterial phagocytosis, increased NLRP3 inflammasome activation and IL-1ß secretion in response to an inflammatory stimulus. CONCLUSIONS: We demonstrate that GPR65 exerts its effects through Gαs- and Gα12/13-mediated pathways, that the IBD-associated GPR65*231Leu allele has compromised interactions with Gα12/13 and that KD of GPR65 leads to impaired bacterial phagocytosis and increased inflammatory signalling via the NLRP3 inflammasome. This work identifies a target for development of small molecule therapies.

Serum Analyte Profiles Associated With Crohn's Disease and Disease Location.

BACKGROUND: Crohn's disease (CD) can affect any segment of the digestive tract but is most often localized in the ileal, ileocolonic, and colorectal regions of the intestines. It is believed that the chronic inflammation in CD is a result of an imbalance between the epithelial barrier, the immune system, and the intestinal microbiota. The aim of the study was to identify circulating markers associated with CD and/or disease location in CD patients. METHODS: We tested 49 cytokines, chemokines, and growth factors in serum samples from 300 patients with CD and 300 controls. After quality control, analyte levels were tested for association with CD and disease location. RESULTS: We identified 13 analytes that were higher in CD patients relative to healthy controls and that remained significant after conservative Bonferroni correction (P < 0.0015). In particular, CXCL9, CXCL1, and interleukin IL-6 had the greatest effect and were highly significant (P < 5 × 10-7). We also identified 9 analytes that were associated with disease location, with VEGF, IL-12p70, and IL-6 being elevated in patients with colorectal disease (P < 3 × 10-4). CONCLUSIONS: Multiple serum analytes are elevated in CD. These implicate the involvement of multiple cell types from the immune, epithelial, and endothelial systems, suggesting that circulating analytes reflect the inflammatory processes that are ongoing within the gut. Moreover, the identification of distinct profiles according to disease location supports the existence of a biological difference between ileal and colonic CD, consistent with previous genetic and clinical observations.

Functional screen of inflammatory bowel disease genes reveals key epithelial functions.

BACKGROUND: Genetic studies have been tremendously successful in identifying genomic regions associated with a wide variety of phenotypes, although the success of these studies in identifying causal genes, their variants, and their functional impacts has been more limited. METHODS: We identified 145 genes from IBD-associated genomic loci having endogenous expression within the intestinal epithelial cell compartment. We evaluated the impact of lentiviral transfer of the open reading frame (ORF) of these IBD genes into the HT-29 intestinal epithelial cell line via transcriptomic analyses. By comparing the genes in which expression was modulated by each ORF, as well as the functions enriched within these gene lists, we identified ORFs with shared impacts and their putative disease-relevant biological functions. RESULTS: Analysis of the transcriptomic data for cell lines expressing the ORFs for known causal genes such as HNF4a, IFIH1, and SMAD3 identified functions consistent with what is already known for these genes. These analyses also identified two major clusters of genes: Cluster 1 contained the known IBD causal genes IFIH1, SBNO2, NFKB1, and NOD2, as well as genes from other IBD loci (ZFP36L1, IRF1, GIGYF1, OTUD3, AIRE and PITX1), whereas Cluster 2 contained the known causal gene KSR1 and implicated DUSP16 from another IBD locus. Our analyses highlight how multiple IBD gene candidates can impact on epithelial structure and function, including the protection of the mucosa from intestinal microbiota, and demonstrate that DUSP16 acts a regulator of MAPK activity and contributes to mucosal defense, in part via its regulation of the polymeric immunoglobulin receptor, involved in the protection of the intestinal mucosa from enteric microbiota. CONCLUSIONS: This functional screen, based on expressing IBD genes within an appropriate cellular context, in this instance intestinal epithelial cells, resulted in changes to the cell's transcriptome that are relevant to their endogenous biological function(s). This not only helped in identifying likely causal genes within genetic loci but also provided insight into their biological functions. Furthermore, this work has highlighted the central role of intestinal epithelial cells in IBD pathophysiology, providing a scientific rationale for a drug development strategy that targets epithelial functions in addition to the current therapies targeting immune functions.

A transcriptome-based approach to identify functional modules within and across primary human immune cells.

Genome-wide transcriptomic analyses have provided valuable insight into fundamental biology and disease pathophysiology. Many studies have taken advantage of the correlation in the expression patterns of the transcriptome to infer a potential biologic function of uncharacterized genes, and multiple groups have examined the relationship between co-expression, co-regulation, and gene function on a broader scale. Given the unique characteristics of immune cells circulating in the blood, we were interested in determining whether it was possible to identify functional co-expression modules in human immune cells. Specifically, we sequenced the transcriptome of nine immune cell types from peripheral blood cells of healthy donors and, using a combination of global and targeted analyses of genes within co-expression modules, we were able to determine functions for these modules that were cell lineage-specific or shared among multiple cell lineages. In addition, our analyses identified transcription factors likely important for immune cell lineage commitment and/or maintenance.