Unveiling the key genes, environmental toxins, and drug exposures in modulating the severity of ulcerative colitis: a comprehensive analysis

Abstract

BackgroundAs yet, the genetic abnormalities involved in the exacerbation of Ulcerative colitis (UC) have not been adequately explored based on bioinformatic methods.Materials and methodsThe gene microarray data and clinical information were downloaded from Gene Expression Omnibus (GEO) repository. The scale-free gene co-expression networks were constructed by R package “WGCNA”. Gene enrichment analysis was performed via Metascape database. Differential expression analysis was performed using “Limma” R package. The “randomForest” packages in R was used to construct the random forest model. Unsupervised clustering analysis performed by “ConsensusClusterPlus”R package was utilized to identify different subtypes of UC patients. Heat map was established using the R package “pheatmap”. Diagnostic parameter capability was evaluated by ROC curve. The”XSum”packages in R was used to screen out small-molecule drugs for the exacerbation of UC based on cMap database. Molecular docking was performed with Schrodinger molecular docking software.ResultsVia WGCNA, a total 77 high Mayo score-associated genes specific in UC were identified. Subsequently, the 9 gene signatures of the exacerbation of UC was screened out by random forest algorithm and Limma analysis, including BGN,CHST15,CYYR1,GPR137B,GPR4,ITGA5,LILRB1,SLFN11 and ST3GAL2. The ROC curve suggested good predictive performance of the signatures for exacerbation of UC in both the training set and the validation set. We generated a novel genotyping scheme based on the 9 signatures. The percentage of patients achieved remission after 4 weeks intravenous corticosteroids (CS-IV) treatment was higher in cluster C1 than that in cluster C2 (54% vs. 27%, Chi-square test, p=0.02). Energy metabolism-associated signaling pathways were significantly up-regulated in cluster C1, including the oxidative phosphorylation, pentose and glucuronate interconversions and citrate cycle TCA cycle pathways. The cluster C2 had a significant higher level of CD4+ T cells. The”XSum”algorithm revealed that Exisulind has a therapeutic potential for UC. Exisulind showed a good binding affinity for GPR4, ST3GAL2 and LILRB1 protein with the docking glide scores of –7.400 kcal/mol, –7.191 kcal/mol and –6.721 kcal/mol, respectively.We also provided a comprehensive review of the environmental toxins and drug exposures that potentially impact the progression of UC.ConclusionUsing WGCNA and random forest algorithm, we identified 9 gene signatures of the exacerbation of UC. A novel genotyping scheme was constructed to predict the severity of UC and screen UC patients suitable for CS-IV treatment. Subsequently, we identified a small molecule drug (Exisulind) with potential therapeutic effects for UC. Thus, our study provided new ideas and materials for the personalized clinical treatment plans for patients with UC.

Mechanism Exploration on the Immunoregulation of Allogeneic Heart Transplantation Rejection in Rats With Exosome miRNA and Proteins From Overexpressed IDO1 BMSCs

“…The screening criteria were set as|logFC| > 0.5, and P < 0.05. DEPs and DEmiRNAs were shown through volcano map and heatmap, which were plotted by ggplot2 42 and pheatmap 43 packages, respectively.…”

Construction of a Lung Adenocarcinoma Prognostic Model Utilizing Serine and Glycine Metabolism-Related Genes

“…Hence, based on the ssGSEA method, the differences in immune-related cells and immune-related functions between HR and LR groups were explored. Pheatmap R package was used to exhibit the antitumor immune enrichment status of the two groups. Wilcoxon test was carried out to compare the differences in the expression of immune checkpoints between individuals in the two groups.…”

Integrating machine learning algorithms and single-cell analysis to identify gut microbiota-related macrophage biomarkers in atherosclerotic plaques

“…Differential expression analysis was performed on the GSE43292 dataset for atheroma plaque and macroscopically intact tissue samples using the limma package, resulting in 188 upregulated genes in atherosclerosis plaques with criteria of |logFC| >0.8 and adjusted P-value <0.05. WGCNA algorithm was used to calculate gene significance (GS) and module membership (MM), confirming a key atherosclerosis module (blue module with 1115 genes) ( Wang et al., 2023 ). Immune infiltration analysis was conducted using the deconvo_tme function from the IOBR package, and CIBERSORT was applied to WGCNA to calculate immune cell-related modules for M0, M1, and M2 macrophages, which were blue, red, and yellow, respectively ( Wu et al., 2022 ).…”