Dataset Information

Characterization of terminal-ileal and colonic Crohn's disease in treatment-naive paediatric patients based on transcriptomic profile using logistic regression.

ABSTRACT:

Background

Inflammatory bowel disease (IBD) is a chronic and idiopathic inflammatory disorder of the gastrointestinal tract and comprises ulcerative colitis (UC) and Crohn's disease (CD). Crohn's disease can affect any part of the gastrointestinal tract, but mainly the terminal ileum and colon. In the present study, we aimed to characterize terminal-ileal CD (ICD) and colonic CD (CCD) at the molecular level, which might enable a more optimized approach for the clinical care and scientific research of CD.

Methods

We analyzed differentially expressed genes in samples from 23 treatment-naïve paediatric patients with CD and 25 non-IBD controls, and compared the data with previously published RNA-Seq data using multi-statistical tests and confidence intervals. We implemented functional profiling and proposed statistical methods for feature selection using a logistic regression model to identify genes that are highly associated in ICD or CCD. We also validated our final candidate genes in independent paediatric and adult cohorts.

Results

We identified 550 genes specifically expressed in patients with CD compared with those in healthy controls (p < 0.05). Among these DEGs, 240 from patients with CCD were mainly involved in mitochondrial dysfunction, whereas 310 from patients with ICD were enriched in the ileum functions such as digestion, absorption, and metabolism. To choose the most effective gene set, we selected the most powerful genes (p-value ≤ 0.05, accuracy ≥ 0.8, and AUC ≥ 0.8) using logistic regression. Consequently, 33 genes were identified as useful for discriminating CD location; the accuracy and AUC were 0.86 and 0.83, respectively. We then validated the 33 genes with data from another independent paediatric cohort (accuracy = 0.93, AUC = 0.92) and adult cohort (accuracy = 0.88, AUC = 0.72).

Conclusions

In summary, we identified DEGs that are specifically expressed in CCD and ICD compared with those in healthy controls and patients with UC. Based on the feature selection analysis, 33 genes were identified as useful for discriminating CCD and ICD with high accuracy and AUC, for not only paediatric patients but also independent cohorts. We propose that our approach and the final gene set are useful for the molecular classification of patients with CD, and it could be beneficial in treatments based on disease location.

SUBMITTER: Park I

PROVIDER: S-EPMC8185924 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Background and aimsCrohn'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.MethodsWe 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.ResultsIn 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.ConclusionsIleal 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.

Project description:In countries with a high incidence of tuberculosis, the typical clinical features of Crohn's disease (CD) may be covered up after tuberculosis infection, and the identification of atypical Crohn's disease and intestinal tuberculosis (ITB) is still a dilemma for clinicians. Least absolute shrinkage and selection operator (LASSO) regression has been applied to select variables in disease diagnosis. However, its value in discriminating ITB and atypical Crohn's disease remains unknown. A total of 400 patients were enrolled from January 2014 to January 2019 in second Xiangya hospital Central South University.Among them, 57 indicators including clinical manifestations, laboratory results, endoscopic findings, computed tomography enterography features were collected for further analysis. R software version 3.6.1 (glmnet package) was used to perform the LASSO logistic regression analysis. SPSS 20.0 was used to perform Pearson chi-square test and binary logistic regression analysis. In the variable selection step, LASSO regression and Pearson chi-square test were applied to select the most valuable variables as candidates for further logistic regression analysis. Secondly, variables identified from step 1 were applied to construct binary logistic regression analysis. Receiver operating characteristic (ROC) curve analysis was performed on these models to assess the ability and the optimal cutoff value for diagnosis. The area under the ROC curve (AUC), sensitivity, specificity, positive predictive value (PPV), negative predictive value (NPV), accuracy rate, together with their 95% confidence and intervals (CIs) were calculated. MedCalc software (Version 16.8) was applied to analyze the ROC curves of models. 332 patients were eventually enrolled to build a binary logistic regression model to discriminate CD (including comprehensive CD and tuberculosis infected CD) and ITB. However, we did not get a satisfactory diagnostic value via applying the binary logistic regression model of comprehensive CD and ITB to predict tuberculosis infected CD and ITB (accuracy rate:79.2%VS 65.1%). Therefore, we further established a binary logistic regression model to discriminate atypical CD from ITB, based on Pearsonchi-square test (model1) and LASSO regression (model 2). Model 1 showed 89.9% specificity, 65.9% sensitivity, 88.5% PPV, 68.9% NPV, 76.9% diagnostic accuracy, and an AUC value of 0.811, and model 2 showed 80.6% specificity, 84.4% sensitivity, 82.3% PPV, 82.9% NPV, 82.6% diagnostic accuracy, and an AUC value of 0.887. The comparison of AUCs between model1 and model2 was statistically different (P < 0.05). Tuberculosis infection increases the difficulty of discriminating CD from ITB. LASSO regression showed a more efficient ability than Pearson chi-square test based logistic regression on differential diagnosing atypical CD and ITB.

Project description:BackgroundThe development and course of Crohn's disease (CD) is related to both genetic and environmental factors. Smoking has been found to exacerbate the course of CD by increasing the risk of developing fistulas and strictures as well as the need for surgery, possibly because of an interaction between smoking or nicotine on macrophage function and the intestinal microvasculature. Several genes are involved in the pathogenesis of CD, and in this study the gene expression differences of the descending colonic mucosa were investigated in CD (smokers or never smokers) and controls (smokers or never smokers).AimTo identify any difference in gene expression of the descending colonic mucosa between smoking and never-smoking CD patients (and controls) by determining genetic expression profiles from microarray analysis.MethodsFifty-seven specimens were obtained by routine colonoscopy from the included material: CD smokers (n = 28) or never-smokers (n = 14) as compared to fifteen healthy controls (8 smokers and 7 never-smokers). RNA was isolated and gene expression assessed with Affymetrix GeneChip Human Genome U133 Plus 2.0. Data were analyzed by principal component analysis (PCA), Wilcoxon rank sum test and multiple linear regressions. Real-time (RT) PCR was subsequently applied to verify microarray results.ResultsThe PCA analysis showed no intrinsic clustering of smokers versus never-smokers. However, when Wilcoxon rank sum test corrected with Q values were performed, six known genes were significantly expressed differently in the inflamed CD smokers as compared to the inflamed CD never-smokers: ring finger protein 138 (RNF138), metalothionein 2A (MT2A) and six transmembrane epithelial antigen of the prostate 3 (STEAP3), SA hypertension-associated homolog, PGM2L1 and KCNJ2. The subsequent RT-PCR-analyses verified, however, that only RNF138, MT2A and STEAP3 were significantly up-regulated in CD smokers in specimens with inflammatory activity of the descending colon.ConclusionsThe present study demonstrates that the genes, RNF138, MT2A, and STEAP3 are differently expressed in the inflamed descending colon of smoking versus never-smoking CD patients, which might be of relevance for the poorer clinical course among CD smokers. Many gastroenterologists are still not totally aware of the benefits of smoking cessation in relation to CD, and do not put much effort into getting the patients to quit, therefore more information on the negative effects of smoking, seems warranted.