Immunological Profiling of Paediatric Inflammatory Bowel Disease Using Unsupervised Machine Learning.

OBJECTIVES: The current classification of inflammatory bowel disease (IBD) is based on clinical phenotypes, which is blind to the molecular basis of the disease. The aim of this study was to stratify a treatment-naïve paediatric IBD cohort through specific innate immunity pathway profiling and application of unsupervised machine learning (UML). METHODS: In order to test the molecular integrity of biological pathways implicated in IBD, innate immune responses were assessed at diagnosis in 22 paediatric patients and 10 age-matched controls. Peripheral blood mononuclear cells (PBMCs) were selectively stimulated for assessing the functionality of upstream activation receptors including NOD2, toll-like receptor (TLR) 1-2 and TLR4, and the downstream cytokine responses (IL-10, IL-1ß, IL-6, and TNF-α) using multiplex assays. Cytokine data generated were subjected to hierarchical clustering to assess for patient stratification. RESULTS: Combined immune responses in patients across 12 effector responses were significantly reduced compared with controls (P = 0.003) and driven primarily by "hypofunctional" TLR responses (P values 0.045, 0.010, and 0.018 for TLR4-mediated IL-10, IL-1ß, and TNF-α, respectively; 0.018 and 0.015 for TLR1-2 -mediated IL-10 and IL-1ß). Hierarchical clustering generated 3 distinct clusters of patients and a fourth group of "unclustered" individuals. No relationship was observed between the observed immune clusters and the clinical disease phenotype. CONCLUSIONS: Although a clinically useful outcome was not observed through hierarchical clustering, our study provides a rationale for using an UML approach to stratify patients. The study also highlights the predominance of hypo-inflammatory innate immune responses as a key mechanism in the pathogenesis of IBD.

Deleterious Genetic Variation Across the NOD Signaling Pathway Is Associated With Reduced NFKB Signaling Transcription and Upregulation of Alternative Inflammatory Transcripts in Pediatric Inflammatory Bowel Disease.

BACKGROUND: Inflammatory bowel disease may arise with inadequate immune response to intestinal bacteria. NOD2 is an established gene in Crohn's disease pathogenesis, with deleterious variation associated with reduced NFKB signaling. We hypothesized that deleterious variation across the NOD2 signaling pathway impacts on transcription. METHODS: Treatment-naïve pediatric inflammatory bowel disease patients had ileal biopsies for targeted autoimmune RNA-sequencing and blood for whole exome sequencing collected at diagnostic endoscopy. Utilizing GenePy, a per-individual, per-gene score, genes within the NOD signaling pathway were assigned a quantitative score representing total variant burden. Where multiple genes formed complexes, GenePy scores were summed to create a "complex" score. Normalized transcript expression of 95 genes within this pathway was retrieved. Regression analysis was performed to determine the impact of genomic variation on gene transcription. RESULTS: Thirty-nine patients were included. Limited clustering of patients based on NOD signaling transcripts was related to underlying genomic variation. Patients harboring deleterious variation in NOD2 had reduced NOD2 (ß = -0.702, P = 4.3 × 10-5) and increased NFKBIA (ß = 0.486, P = .001), reflecting reduced NFKB signal activation. Deleterious variation in the NOD2-RIPK2 complex was associated with increased NLRP3 (ß = 0.8, P = 3.1475 × 10-8) and TXN (ß = -0.417, P = 8.4 × 10-5) transcription, components of the NLRP3 inflammasome. Deleterious variation in the TAK1-TAB complex resulted in reduced MAPK14 transcription (ß = -0.677, P = 1.7 × 10-5), a key signal transduction protein in the NOD2 signaling cascade and increased IFNA1 (ß = 0.479, P = .001), indicating reduced transcription of NFKB activators and alternative interferon transcription in these patients. CONCLUSIONS: Data integration identified perturbation of NOD2 signaling transcription correlated with genomic variation. A hypoimmune NFKB signaling transcription response was observed. Alternative inflammatory pathways were activated and may represent therapeutic targets in specific patients.

Ileal Transcriptomic Analysis in Paediatric Crohn's Disease Reveals IL17- and NOD-signalling Expression Signatures in Treatment-naïve Patients and Identifies Epithelial Cells Driving Differentially Expressed Genes.

BACKGROUND AND AIMS: Crohn's disease [CD] arises through host-environment interaction. Abnormal gene expression results from disturbed pathway activation or response to bacteria. We aimed to determine activated pathways and driving cell types in paediatric CD. METHODS: We employed contemporary targeted autoimmune RNA sequencing, in parallel to single-cell sequencing, to ileal tissue derived from paediatric CD and controls. Weighted gene co-expression network analysis [WGCNA] was performed and differentially expressed genes [DEGs] were determined. We integrated clinical data to determine co-expression modules associated with outcomes. RESULTS: In all, 27 treatment-naive CD [TN-CD], 26 established CD patients and 17 controls were included. WGCNA revealed a 31-gene signature characterising TN-CD patients, but not established CD, nor controls. The CSF3R gene is a hub within this module and is key in neutrophil expansion and differentiation. Antimicrobial genes, including S100A12 and the calprotectin subunit S100A9, were significantly upregulated in TN CD compared with controls [p = 2.61 x 10-15 and p = 9.13 x 10-14, respectively] and established CD [both p = 0.0055]. Gene-enrichment analysis confirmed upregulation of the IL17-, NOD- and Oncostatin-M-signalling pathways in TN-CD patients, identified in both WGCNA and DEG analyses. An upregulated gene signature was enriched for transcripts promoting Th17-cell differentiation and correlated with prolonged time to relapse [correlation-coefficient-0.36, p = 0.07]. Single-cell sequencing of TN-CD patients identified specialised epithelial cells driving differential expression of S100A9. Cell groups, determined by single-cell gene expression, demonstrated enrichment of IL17-signalling in monocytes and epithelial cells. CONCLUSIONS: Ileal tissue from treatment-naïve paediatric patients is significantly upregulated for genes driving IL17-, NOD- and Oncostatin-M-signalling. This signal is driven by a distinct subset of epithelial cells expressing antimicrobial gene transcripts.

Expression profile of the matricellular protein periostin in paediatric inflammatory bowel disease.

The precise role of periostin, an extra-cellular matrix protein, in inflammatory bowel disease (IBD) is unclear. Here, we investigated periostin in paediatric IBD including its relationship with disease activity, clinical outcomes, genomic variation and expression in the colonic tissue. Plasma periostin was analysed using ELISA in 144 paediatric patients and 38 controls. Plasma levels were assessed against validated disease activity indices in IBD and clinical outcomes. An immuno-fluorescence for periostin and detailed isoform-expression analysis in the colonic tissue was performed in 23 individuals. We integrated a whole-gene based burden metric 'GenePy' to assess the impact of variation in POSTN and 23 other genes functionally connected to periostin. We found that plasma periostin levels were significantly increased during remission compared to active Crohn's disease. The immuno-fluorescence analysis demonstrated enhanced peri-cryptal ring patterns in patients compared to controls, present throughout inflamed, as well as macroscopically non-inflamed colonic tissue. Interestingly, the pattern of isoforms remained unchanged during bowel inflammation compared to healthy controls. In addition to its role during the inflammatory processes in IBD, periostin may have an additional prominent role in mucosal repair. Additional studies will be necessary to understand its role in the pathogenesis, repair and fibrosis in IBD.

Genetic Sequencing of Pediatric Patients Identifies Mutations in Monogenic Inflammatory Bowel Disease Genes that Translate to Distinct Clinical Phenotypes.

OBJECTIVES: Monogenic inflammatory bowel disease (IBD) comprises rare Mendelian causes of gut inflammation, often presenting in infants with severe and atypical disease. This study aimed to identify clinically relevant variants within 68 monogenic IBD genes in an unselected pediatric IBD cohort. METHODS: Whole exome sequencing was performed on patients with pediatric-onset disease. Variants fulfilling the American College of Medical Genetics criteria as "pathogenic" or "likely pathogenic" were assessed against phenotype at diagnosis and follow-up. Individual patient variants were assessed and processed to generate a per-gene, per-individual, deleteriousness score. RESULTS: Four hundred one patients were included, and the median age of disease-onset was 11.92 years. In total, 11.5% of patients harbored a monogenic variant. TRIM22-related disease was implicated in 5 patients. A pathogenic mutation in the Wiskott-Aldrich syndrome (WAS) gene was confirmed in 2 male children with severe pancolonic inflammation and primary sclerosing cholangitis. In total, 7.3% of patients with Crohn's disease had apparent autosomal recessive, monogenic NOD2-related disease. Compared with non-NOD2 Crohn's disease, these patients had a marked stricturing phenotype (odds ratio 11.52, significant after correction for disease location) and had undergone significantly more intestinal resections (odds ratio 10.75). Variants in ADA, FERMT1, and LRBA did not meet the criteria for monogenic disease in any patients; however, case-control analysis of mutation burden significantly implicated these genes in disease etiology. DISCUSSION: Routine whole exome sequencing in pediatric patients with IBD results in a precise molecular diagnosis for a subset of patients with IBD, providing the opportunity to personalize therapy. NOD2 status informs risk of stricturing disease requiring surgery, allowing clinicians to direct prognosis and intervention.

16S sequencing and functional analysis of the fecal microbiome during treatment of newly diagnosed pediatric inflammatory bowel disease.

The human microbiome is of considerable interest to pediatric inflammatory bowel disease (PIBD) researchers with 1 potential mechanism for disease development being aberrant immune handling of the intestinal bacteria. This study analyses the fecal microbiome through treatment in newly diagnosed PIBD patients and compares to cohabiting siblings where possible. Patients were recruited on clinical suspicion of PIBD before diagnosis. Treatment-naïve fecal samples were collected, with further samples at 2 and 6 weeks into treatment. Samples underwent 16S ribosomal ribonucleic acid (RNA) gene sequencing and short-chain fatty acids (SCFAs) analysis, results were analyzed using quantitative-insights-into-microbial-ecology. Six PIBD patients were included in the cohort: 4 Crohn disease (CD), 1 ulcerative colitis (UC), 1 inflammatory bowel disease (IBD) unclassified, and median age 12.6 (range 10-15.1 years); 3 patients had an unaffected healthy sibling recruited. Microbial diversity (observed species/Chao1/Shannon diversity) was reduced in treatment-naïve patients compared to siblings and patients in remission. Principal coordinate analysis using Bray-Curtis dissimilarity and UniFrac revealed microbial shifts in CD over the treatment course. In treatment-naïve PIBD, there was reduction in functional ability for amino acid metabolism and carbohydrate handling compared to controls (P = .038) and patients in remission (P = .027). Metabolic function returned to normal after remission was achieved. SCFA revealed consistent detection of lactate in treatment-naïve samples. This study adds in-depth 16S rRNA sequencing analysis on a small longitudinal cohort to the literature and includes sibling controls and patients with UC/IBD unclassified. It highlights the initial dysbiosis, reduced diversity, altered functional potential, and subsequent shifts in bacteria from diagnosis over time to remission.

Exome Analysis of Rare and Common Variants within the NOD Signaling Pathway.

Pediatric inflammatory bowel disease (pIBD) is a chronic heterogeneous disorder. This study looks at the burden of common and rare coding mutations within 41 genes comprising the NOD signaling pathway in pIBD patients. 136 pIBD and 106 control samples underwent whole-exome sequencing. We compared the burden of common, rare and private mutation between these two groups using the SKAT-O test. An independent replication cohort of 33 cases and 111 controls was used to validate significant findings. We observed variation in 40 of 41 genes comprising the NOD signaling pathway. Four genes were significantly associated with disease in the discovery cohort (BIRC2 p = 0.004, NFKB1 p = 0.005, NOD2 p = 0.029 and SUGT1 p = 0.047). Statistical significance was replicated for BIRC2 (p = 0.041) and NOD2 (p = 0.045) in an independent validation cohort. A gene based test on the combined discovery and replication cohort confirmed association for BIRC2 (p = 0.030). We successfully applied burden of mutation testing that jointly assesses common and rare variants, identifying two previously implicated genes (NFKB1 and NOD2) and confirmed a possible role in disease risk in a previously unreported gene (BIRC2). The identification of this novel gene provides a wider role for the inhibitor of apoptosis gene family in IBD pathogenesis.

Identification of Variants in Genes Associated with Single-gene Inflammatory Bowel Disease by Whole-exome Sequencing.

BACKGROUND: Most cases of inflammatory bowel disease (IBD) are caused by complex host-environment interaction. There are a number of conditions associated with a single-gene mutation, most cases are very early onset (aged < 6 yr), present with a unique form of disease and often have atypical features. METHODS: Whole-exome data for 147 pediatric patients with IBD were interrogated for a panel of 51 genes associated with monogenic IBD. Observed variation was categorized according to the American College of Medical Genetics (ACMG) guidelines to identify rare, novel, and known variants that might contribute to IBD. RESULTS: Five hundred seventy-four variants were identified across 51 genes. These were categorized in line with ACMG guidance to remove benign variants and to identify "pathogenic" and "likely pathogenic" variants. In 6 patients, we observed 6 pathogenic variants of which CYBA(c.287+2T>C), COL7A1(c.6501+1G>C), LIG4(p.R814X), and XIAP(p.T470S) were known causative mutations, and FERMT1(p.R271Q) and SKIV2L(c.354+5G>A) were novel. In the 3 patients with XIAP, SKIV2L, and FERMT1 variants, individuals' disease features resembled the monogenic phenotype. This was despite apparent heterozygous carriage of pathogenic variation for the latter 2 genes. The XIAP variant was observed in a hemizygous male. CONCLUSIONS: Whole-exome sequencing allows for identification of known and de novo potentially causative mutations in genes associated with monogenic IBD. Although these are rare conditions, it is vital to identify causative mutations early to improve prognosis. We postulate that in a subset of IBD, heterozygous mutations (in genes believed to manifest IBD through autosomal recessive inheritance) may contribute to clinical presentation.

Exome analysis of patients with concurrent pediatric inflammatory bowel disease and autoimmune disease.

BACKGROUND: Pediatric Inflammatory Bowel Disease (PIBD) is a chronic condition seen in genetically predisposed individuals. Genome-wide association studies have implicated >160 genomic loci in IBD with many genes coding for proteins in key immune pathways. This study looks at autoimmune disease burden in patients diagnosed with PIBD and interrogates exome data of a subset of patients. METHODS: Patients were recruited from the Southampton Genetics of PIBD cohort. Clinical diagnosis of autoimmune disease in these individuals was ascertained from medical records. For a subset of patients with PIBD and concurrent asthma, exome data was interrogated to ascertain the burden of pathogenic variants within genes implicated in asthma. Association testing was conducted between cases and population controls using the SKAT-O test. RESULTS: Forty-nine (28.3%) PIBD children (18.49% CD, 8.6% UC, and 21.15% IBDU patients) had a concurrent clinical diagnosis of at least one other autoimmune disorder; asthma was the most prevalent, affecting 16.2% of the PIBD cohort. Rare and common variant association testing revealed 6 significant genes (P < 0.05) before Bonferroni adjustment. Three of these genes were previously implicated in both asthma and IBD (ZPBP2 IL1R1, and IL18R1) and 3 in asthma only (PYHIN1, IL2RB, and GSTP1). CONCLUSIONS: One-third of our cohort had a concurrent autoimmune condition. We observed higher incidence of asthma compared with the overall pediatric prevalence. Despite a small sample size, SKAT-O evaluated a significant burden of rare and common mutations in 6 genes. Variant burden suggests that a systemic immune dysregulation rather than organ-specific could underpin immune dysfunction for a subset of patients.