IBD risk loci are enriched in multigenic regulatory modules encompassing putative causative genes.

GWAS have identified >200 risk loci for Inflammatory Bowel Disease (IBD). The majority of disease associations are known to be driven by regulatory variants. To identify the putative causative genes that are perturbed by these variants, we generate a large transcriptome data set (nine disease-relevant cell types) and identify 23,650 cis-eQTL. We show that these are determined by ∼9720 regulatory modules, of which ∼3000 operate in multiple tissues and ∼970 on multiple genes. We identify regulatory modules that drive the disease association for 63 of the 200 risk loci, and show that these are enriched in multigenic modules. Based on these analyses, we resequence 45 of the corresponding 100 candidate genes in 6600 Crohn disease (CD) cases and 5500 controls, and show with burden tests that they include likely causative genes. Our analyses indicate that ≥10-fold larger sample sizes will be required to demonstrate the causality of individual genes using this approach.

Analysis of Genes Associated With Monogenic Primary Immunodeficiency Identifies Rare Variants in XIAP in Patients With Crohn's Disease.

BACKGROUND & AIMS: A few rare monogenic primary immunodeficiencies (PIDs) are characterized by chronic intestinal inflammation that resembles Crohn's disease (CD). We investigated whether 23 genes associated with 10 of these monogenic disorders contain common, low-frequency, or rare variants that increase risk for CD. METHODS: Common and low frequency variants in 1 Mb loci centered on the candidate genes were analyzed using meta-data corresponding to genotypes of approximately 17,000 patients with CD or without CD (controls) in Europe. The contribution of rare variants was assessed by high-throughput sequencing of 4750 individuals, including 660 early-onset and/or familial cases among the 2390 patients with CD. Variants were expressed from vectors in SW480 or HeLa cells and functions of their products were analyzed in immunofluorescence, luciferase, immunoprecipitation, and immunoblot assays. RESULTS: We reproduced the association of the interleukin 10 locus with CD (P = .007), although none of the significantly associated variants modified the coding sequence of interleukin 10. We found XIAP to be significantly enriched for rare coding mutations in patients with CD vs controls (P = .02). We identified 4 previously unreported missense variants associated with CD. Variants in XIAP cause the PID X-linked lymphoproliferative disease type 2, yet none of the carriers of these variants had all the clinical features of X-linked lymphoproliferative disease type 2. Identified XIAP variants S123N, R233Q, and P257A were associated with an impaired activation of NOD2 signaling after muramyl dipeptide stimulation. CONCLUSIONS: In a systematic analysis of variants in 23 PID-associated genes, we confirmed the association of variants in XIAP with CD. Further screenings for CD-associated variants and analyses of their functions could increase our understanding of the relationship between PID-associated genes and CD pathogenesis.

CCDC88B is required for pathogenesis of inflammatory bowel disease.

Inflammatory bowel disease (IBD) involves interaction between host genetic factors and environmental triggers. CCDC88B maps within one IBD risk locus on human chromosome 11q13. Here we show that CCDC88B protein increases in the colon during intestinal injury, concomitant with an influx of CCDC88B+lymphoid and myeloid cells. Loss of Ccdc88b protects against DSS-induced colitis, with fewer pathological lesions and reduced intestinal inflammation in Ccdc88b-deficient mice. In a T cell transfer model of colitis, Ccdc88b mutant CD4+ T cells do not induce colitis in immunocompromised hosts. Expression of human CCDC88B RNA and protein is higher in IBD patient colons than in control colon tissue. In human CD14+ myeloid cells, CCDC88B is regulated by cis-acting variants. In a cohort of patients with Crohn's disease, CCDC88B expression correlates positively with disease risk. These findings suggest that CCDC88B has a critical function in colon inflammation and the pathogenesis of IBD.Hook-related protein family member CCDC88b is encoded by a locus that has been associated with inflammatory bowel disease. Here the authors show that Ccdc88b inactivation in T cells prevents colitis in a transfer model, and detect high colonic levels of CCDC88b in patients with Crohn disease or ulcerative colitis, identifying that expression correlates with disease risk.

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.

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).

Genome-Wide Copy Number Variation Scan Identifies Complement Component C4 as Novel Susceptibility Gene for Crohn's Disease.

BACKGROUND: The genetic component of Crohn's disease (CD) is well known, with 140 susceptibility loci identified so far. In addition to single nucleotide polymorphisms typically studied in genome-wide scans, copy number variation is responsible for a large proportion of human genetic variation. METHODS: We performed a genome-wide search for copy number variants associated with CD using array comparative genomic hybridization. One of the found regions was validated independently through real-time PCR. Serum levels of the found gene were measured in patients and control subjects. RESULTS: We found copy number differences for the C4S and C4L gene variants of complement component C4 in the central major histocompatibility complex region on chromosome 6p21. Specifically, we saw that CD patients tend to have lower C4L and higher C4S copies than control subjects (P = 5.00 × 10 and P = 9.11 × 10), which was independent of known associated classical HLA I and II alleles (P = 7.68 × 10 and P = 6.29 × 10). Although C4 serum levels were not different between patients and control subjects, the relationship between C4 copy number and serum level was different for patients and control subjects with higher copy numbers leading to higher serum concentrations in control subjects, compared with CD patients (P < 0.001). CONCLUSIONS: C4 is part of the classical activation pathway of the complement system, which is important for (auto)immunity. Low C4L or high C4S copy number, and corresponding effects on C4 serum level, could lead to an exaggerated response against infections, possibly leading to (auto)immune disease.

NF-κB-Independent Role of IKKα/IKKß in Preventing RIPK1 Kinase-Dependent Apoptotic and Necroptotic Cell Death during TNF Signaling.

TNF is a master pro-inflammatory cytokine. Activation of TNFR1 by TNF can result in both RIPK1-independent apoptosis and RIPK1 kinase-dependent apoptosis or necroptosis. These cell death outcomes are regulated by two distinct checkpoints during TNFR1 signaling. TNF-mediated NF-κB-dependent induction of pro-survival or anti-apoptotic molecules is a well-known late checkpoint in the pathway, protecting cells from RIPK1-independent death. On the other hand, the molecular mechanism regulating the contribution of RIPK1 to cell death is far less understood. We demonstrate here that the IKK complex phosphorylates RIPK1 at TNFR1 complex I and protects cells from RIPK1 kinase-dependent death, independent of its function in NF-κB activation. We provide in vitro and in vivo evidence that inhibition of IKKα/IKKß or its upstream activators sensitizes cells to death by inducing RIPK1 kinase-dependent apoptosis or necroptosis. We therefore report on an unexpected, NF-κB-independent role for the IKK complex in protecting cells from RIPK1-dependent death downstream of TNFR1.

DUSP3 Genetic Deletion Confers M2-like Macrophage-Dependent Tolerance to Septic Shock.

DUSP3 is a small dual-specificity protein phosphatase with an unknown physiological function. We report that DUSP3 is strongly expressed in human and mouse monocytes and macrophages, and that its deficiency in mice promotes tolerance to LPS-induced endotoxin shock and to polymicrobial septic shock after cecal ligation and puncture. By using adoptive transfer experiments, we demonstrate that resistance to endotoxin is macrophage dependent and transferable, and that this protection is associated with a striking increase of M2-like macrophages in DUSP3(-/-) mice in both the LPS and cecal ligation and puncture models. We show that the altered response of DUSP3(-/-) mice to sepsis is reflected in decreased TNF production and impaired ERK1/2 activation. Our results demonstrate that DUSP3 plays a key and nonredundant role as a regulator of innate immune responses by mechanisms involving the control of ERK1/2 activation, TNF secretion, and macrophage polarization.

Dual-specificity phosphatase 3 deficiency or inhibition limits platelet activation and arterial thrombosis.

BACKGROUND: A limitation of current antiplatelet therapies is their inability to separate thrombotic events from bleeding occurrences. A better understanding of the molecular mechanisms leading to platelet activation is important for the development of improved therapies. Recently, protein tyrosine phosphatases have emerged as critical regulators of platelet function. METHODS AND RESULTS: This is the first report implicating the dual-specificity phosphatase 3 (DUSP3) in platelet signaling and thrombosis. This phosphatase is highly expressed in human and mouse platelets. Platelets from DUSP3-deficient mice displayed a selective impairment of aggregation and granule secretion mediated by the collagen receptor glycoprotein VI and the C-type lectin-like receptor 2. DUSP3-deficient mice were more resistant to collagen- and epinephrine-induced thromboembolism compared with wild-type mice and showed severely impaired thrombus formation on ferric chloride-induced carotid artery injury. Intriguingly, bleeding times were not altered in DUSP3-deficient mice. At the molecular level, DUSP3 deficiency impaired Syk tyrosine phosphorylation, subsequently reducing phosphorylation of phospholipase Cγ2 and calcium fluxes. To investigate DUSP3 function in human platelets, a novel small-molecule inhibitor of DUSP3 was developed. This compound specifically inhibited collagen- and C-type lectin-like receptor 2-induced human platelet aggregation, thereby phenocopying the effect of DUSP3 deficiency in murine cells. CONCLUSIONS: DUSP3 plays a selective and essential role in collagen- and C-type lectin-like receptor 2-mediated platelet activation and thrombus formation in vivo. Inhibition of DUSP3 may prove therapeutic for arterial thrombosis. This is the first time a protein tyrosine phosphatase, implicated in platelet signaling, has been targeted with a small-molecule drug.

Epilepsy, hippocampal sclerosis and febrile seizures linked by common genetic variation around SCN1A.

Epilepsy comprises several syndromes, amongst the most common being mesial temporal lobe epilepsy with hippocampal sclerosis. Seizures in mesial temporal lobe epilepsy with hippocampal sclerosis are typically drug-resistant, and mesial temporal lobe epilepsy with hippocampal sclerosis is frequently associated with important co-morbidities, mandating the search for better understanding and treatment. The cause of mesial temporal lobe epilepsy with hippocampal sclerosis is unknown, but there is an association with childhood febrile seizures. Several rarer epilepsies featuring febrile seizures are caused by mutations in SCN1A, which encodes a brain-expressed sodium channel subunit targeted by many anti-epileptic drugs. We undertook a genome-wide association study in 1018 people with mesial temporal lobe epilepsy with hippocampal sclerosis and 7552 control subjects, with validation in an independent sample set comprising 959 people with mesial temporal lobe epilepsy with hippocampal sclerosis and 3591 control subjects. To dissect out variants related to a history of febrile seizures, we tested cases with mesial temporal lobe epilepsy with hippocampal sclerosis with (overall n = 757) and without (overall n = 803) a history of febrile seizures. Meta-analysis revealed a genome-wide significant association for mesial temporal lobe epilepsy with hippocampal sclerosis with febrile seizures at the sodium channel gene cluster on chromosome 2q24.3 [rs7587026, within an intron of the SCN1A gene, P = 3.36 × 10(-9), odds ratio (A) = 1.42, 95% confidence interval: 1.26-1.59]. In a cohort of 172 individuals with febrile seizures, who did not develop epilepsy during prospective follow-up to age 13 years, and 6456 controls, no association was found for rs7587026 and febrile seizures. These findings suggest SCN1A involvement in a common epilepsy syndrome, give new direction to biological understanding of mesial temporal lobe epilepsy with hippocampal sclerosis with febrile seizures, and open avenues for investigation of prognostic factors and possible prevention of epilepsy in some children with febrile seizures.

An efficient algorithm to perform multiple testing in epistasis screening.

BACKGROUND: Research in epistasis or gene-gene interaction detection for human complex traits has grown over the last few years. It has been marked by promising methodological developments, improved translation efforts of statistical epistasis to biological epistasis and attempts to integrate different omics information sources into the epistasis screening to enhance power. The quest for gene-gene interactions poses severe multiple-testing problems. In this context, the maxT algorithm is one technique to control the false-positive rate. However, the memory needed by this algorithm rises linearly with the amount of hypothesis tests. Gene-gene interaction studies will require a memory proportional to the squared number of SNPs. A genome-wide epistasis search would therefore require terabytes of memory. Hence, cache problems are likely to occur, increasing the computation time. In this work we present a new version of maxT, requiring an amount of memory independent from the number of genetic effects to be investigated. This algorithm was implemented in C++ in our epistasis screening software MBMDR-3.0.3. We evaluate the new implementation in terms of memory efficiency and speed using simulated data. The software is illustrated on real-life data for Crohn's disease. RESULTS: In the case of a binary (affected/unaffected) trait, the parallel workflow of MBMDR-3.0.3 analyzes all gene-gene interactions with a dataset of 100,000 SNPs typed on 1000 individuals within 4 days and 9 hours, using 999 permutations of the trait to assess statistical significance, on a cluster composed of 10 blades, containing each four Quad-Core AMD Opteron(tm) Processor 2352 2.1 GHz. In the case of a continuous trait, a similar run takes 9 days. Our program found 14 SNP-SNP interactions with a multiple-testing corrected p-value of less than 0.05 on real-life Crohn's disease (CD) data. CONCLUSIONS: Our software is the first implementation of the MB-MDR methodology able to solve large-scale SNP-SNP interactions problems within a few days, without using much memory, while adequately controlling the type I error rates. A new implementation to reach genome-wide epistasis screening is under construction. In the context of Crohn's disease, MBMDR-3.0.3 could identify epistasis involving regions that are well known in the field and could be explained from a biological point of view. This demonstrates the power of our software to find relevant phenotype-genotype higher-order associations.

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.

Overexpression of CD39 in mouse airways promotes bacteria-induced inflammation.

In airways, the ecto-nucleoside triphosphate diphosphohydrolase CD39 plays a central role in the regulation of physiological mucosal nucleotide concentrations and likely contributes to the control of inflammation because accelerated ATP metabolism occurs in chronic inflammatory lung diseases. We sought to determine whether constant elevated CD39 activity in lung epithelia is sufficient to cause inflammation and whether this affects the response to acute LPS or Pseudomonas aeruginosa exposure. We generated transgenic mice overexpressing human CD39 under the control of the airway-specific Clara cell 10-kDa protein gene promoter. Transgenic mice did not develop any spontaneous lung inflammation. However, intratracheal instillation of LPS resulted in accelerated recruitment of neutrophils to the airways of transgenic mice. Macrophage clearance was delayed, and the amounts of CD8(+) T and B cells were augmented. Increased levels of keratinocyte chemoattractant, IL-6, and RANTES were produced in transgenic lungs. Similarly, higher numbers of neutrophils and macrophages were found in the lungs of transgenic mice infected with P. aeruginosa, which correlated with improved bacteria clearance. The transgenic phenotype was partially and differentially restored by coinstillation of P2X(1) or P2X(7) receptor antagonists or of caffeine with LPS. Thus, a chronic increase of epithelial CD39 expression and activity promotes airway inflammation in response to bacterial challenge by enhancing P1 and P2 receptor activation.

A P2X ion channel-triggered NF-kappaB pathway enhances TNF-alpha-induced IL-8 expression in airway epithelial cells.

Extracellular ATP, acting at P2Y and P2X receptors, has recently been shown to contribute to airway inflammation. The aim of our study was to investigate the molecular mechanisms involved in the ATP-dependent regulation of IL-8 production by airway epithelial cells. Treatment of human normal tracheal (NT)-1 cells with ATP or its two analogs, alpha,beta-methylene ATP (alpha,beta-meATP) and 2'- and 3'-O-(4-benzoyl-benzoyl)-ATP (BzATP) activated NF-kappaB through the IkappaB kinase (IKK) complex, a process requiring Ca(2+), calmodulin (CaM), and Ca(2+)/CaM-dependent kinase (CaMK), but independent from phospholipase C. alpha,beta-meATP-induced IKK activation also occurred in the alveolar A549 cell line. Real-time RT-PCR revealed that NT-1 and A549 cells expressed P2X(4), P2X(5),and P2X(6) subtype mRNAs, whereas P2X(7) mRNAs were only detected in NT-1 cells. Polarized human primary nasal epithelial cells expressed all four P2X subtypes. Both alpha,beta-meATP and BzATP caused Ca(2+)-dependent binding of phosphorylated p65 (S536) NF-kappaB subunit to the endogenous IL-8 gene promoter in NT-1 cells. Although these agonists did not induce significant IL-8 gene expression by these cells, they markedly enhanced TNF-alpha-induced NF-kappaB activation, resulting in increased IL-8 expression and release. Application of alpha,beta-meATP or BzATP at the apical side of polarized human primary nasal epithelial cells sufficed to cause CaMK-dependent IL-8 release by these cells. Thus, ATP promotes TNF-alpha-elicited IL-8 expression through P2X ion channel-triggered Ca(2+) entry, leading to CaMK-dependent IKK activation and binding of active p65 to IL-8 gene promoter.

Sensitivity of intestinal fibroblasts to TNF-related apoptosis-inducing ligand-mediated apoptosis in Crohn's disease.

OBJECTIVE: Strictures and fistulas are common complications of Crohn's disease (CD). Collagen deposit and fibroblast proliferation can contribute to their development. Tumour necrosis factor (TNF)-related apoptosis-inducing ligand (TRAIL) binds two pro-apoptotic (TRAIL-R1, TRAIL-R2) and three anti-apoptotic (TRAIL-R3, TRAIL-R4, osteoprotegerin (OPG)) receptors. The aim of this work was to study TRAIL expression and the effects on intestinal fibroblasts (IFs) in CD. MATERIAL AND METHODS: Intestinal samples from 25 CD (with or without fibrostenosing areas) and 38 control patients (with or without inflammation) were used. TRAIL, TRAIL R2 and TRAIL R3 expression in the intestine and in human IFs was studied by real-time reverse transcriptase-polymerase chain reaction (RT-PCR) and immunostaining in IF and intestinal samples. TRAIL-induced IF cell death was studied in the presence or absence of OPG and cytokines. Western blots for poly ADP-ribose polymerase (PARP) and caspase-8 were performed to confirm apoptosis in IFs. RESULTS: Transcripts for TRAIL and its receptors were confirmed in the intestine. Immunostaining showed intestinal expression of TRAIL, TRAIL-R2 and TRAIL-R3 in fibroblasts, immune cells and epithelial cells, mainly in fibrostenosing areas. TRAIL-R3 mRNA expression was lower in IFs from fibrostenosing CD. The sensitivity of IFs to TRAIL-mediated apoptosis was higher in the fibrostenosing areas of CD. The effect of TRAIL was decreased by IL-6 and its soluble receptor and almost completely reversed by OPG in the CD patients involved. CONCLUSIONS: TRAIL is expressed in the intestine and influences fibroblast survival. Variations in TRAIL expression and in TRAIL-mediated apoptosis could be involved in the tissue remodelling associated with CD.

The TNF/ADAM 17 system: implication of an ADAM 17 haplotype in the clinical response to infliximab in Crohn's disease.

Infliximab, a chimeric anti-tumour necrosis factor (TNF)-alpha antibody induces a clinical response in 70% of Crohn's disease patients and the response to infliximab therapy could be partially determined by genetic factors. The implication of both transmembrane and soluble forms of the TNF-alpha in the mechanism of action of infliximab has been demonstrated. The aim of our work was first to perform a complete study of TNF variants role in the response to infliximab in Crohn's disease. Secondly, considering the role of ADAM 17 in TNF-alpha shedding, the ADAM 17 locus was also studied. The response to infliximab was evaluated in 222 Caucasian Crohn's disease patients with a luminal (n=160) or fistulizing (n=62) form of the disease. Clinical and biological response evaluation was based on the Crohn's Disease Activity Index score and C-reactive protein level evolutions, respectively. The entire TNF gene was sequenced on the complete cohort. Twelve single nucleotide polymorphisms spanning the ADAM 17 locus were studied and haplotypes rebuilt. A clinical response was observed in 64% of the patients and biological response in 77.1% of patients. No association was found between the TNF gene and the response to infliximab. One haplotype in the ADAM 17 region was associated with a clinical response to infliximab in CD patients (adjusted P=0.045). In conclusion, our results exclude, with a reasonable power, an implication of the TNF gene in the response to infliximab in Crohn's disease, but reveal a potential role of the ADAM 17 gene in this response.