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

Background

As yet, the genetic abnormalities involved in the exacerbation of Ulcerative colitis (UC) have not been adequately explored based on bioinformatic methods.

Materials and methods

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

Results

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

Conclusion

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

Placental Malaria is Associated with Higher LILRB2 Expression in Monocyte Subsets and Lower Anti-Malarial IgG Antibodies During Infancy

Background

Placental malaria (PM) is associated with a higher susceptibility of infants to Plasmodium falciparum (Pf) malaria. A hypothesis of immune tolerance has been suggested but no clear explanation has been provided so far. Our goal was to investigate the involvement of inhibitory receptors LILRB1 and LILRB2, known to drive immune evasion upon ligation with pathogen and/or host ligands, in PM-induced immune tolerance.

Method

Infants of women with or without PM were enrolled in Allada, southern Benin, and followed-up for 24 months. Antibodies with specificity for five blood stage parasite antigens were quantified by ELISA, and the frequency of immune cell subsets was quantified by flow cytometry. LILRB1 or LILRB2 expression was assessed on cells collected at 18 and 24 months of age.

Findings

Infants born to women with PM had a higher risk of developing symptomatic malaria than those born to women without PM (IRR=1.53, p=0.040), and such infants displayed a lower frequency of non-classical monocytes (OR=0.74, p=0.01) that overexpressed LILRB2 (OR=1.36, p=0.002). Moreover, infants born to women with PM had lower levels of cytophilic IgG and higher levels of IL-10 during active infection.

Interpretation

Modulation of IgG and IL-10 levels could impair monocyte functions (opsonisation/phagocytosis) in infants born to women with PM, possibly contributing to their higher susceptibility to malaria. The long-lasting effect of PM on infants’ monocytes was notable, raising questions about the capacity of ligands such as Rifins or HLA-I molecules to bind to LILRB1 and LILRB2 and to modulate immune responses, and about the reprogramming of neonatal monocytes/macrophages.

Natural LILRB1 D1-D2 variants show frequency differences in populations and bind to HLA class I with various avidities

Leukocyte immunoglobulin-like receptor B1 (LILRB1) is widely expressed on various immune cells and the engagement of LILRB1 to HLA class I and pathogen-derived proteins can modulate the immune response. In the current study, 108 LILRB1 alleles were identified by screening the LILRB1 locus from the 1000 Genomes Phase 3 database. Forty-six alleles that occurred in three or more individuals encode 28 LILRB1 allotypes, and the inferred LILRB1 allotypes were then grouped into 9 LILRB1 D1-D2 variants for further analysis. We found that variants 1, 2, and 3 represent the three most frequent LILRB1 D1-D2 variants and the nine variants show frequency differences in populations. The binding assay demonstrated that variant 1 bound to HLA class I with the highest avidity, and all tested LILRB1 D1-D2 variants bound to HLA-C with lower avidity than to HLA-A and -B. Locus-specific polymorphisms at positions 183, 189, and 268 in HLA class I and dimorphisms in HLA-A (positions 207 and 253) and in HLA-B (position 194) affect their binding to LILRB1. Notably, the electrostatic interaction plays a critical role in the binding of LILRB1 to HLA class I as revealed by electrostatic analysis and by comparison of different binding avidities caused by polymorphisms at positions 72 and 103 of LILRB1. In this paper, we present a comprehensive study of the population genetics and binding abilities of LILRB1. The data will help us better understand the LILRB1-related diversity of the immune system and lay a foundation for functional studies.

B Cells Control Mucosal-Associated Invariant T Cell Responses to Salmonella enterica Serovar Typhi Infection Through the CD85j HLA-G Receptor

Mucosal-associated invariant T (MAIT) cells are an innate-like population of T cells that display a TCR Vα7.2+ CD161+ phenotype and are restricted by the nonclassical MHC-related molecule 1 (MR1). Although B cells control MAIT cell development and function, little is known about the mechanisms underlying their interaction(s). Here, we report, for the first time, that during Salmonella enterica serovar Typhi (S. Typhi) infection, HLA-G expression on B cells downregulates IFN-γ production by MAIT cells. In contrast, blocking HLA-G expression on S. Typhi-infected B cells increases IFN-γ production by MAIT cells. After interacting with MAIT cells, kinetic studies show that B cells upregulate HLA-G expression and downregulate the inhibitory HLA-G receptor CD85j on MAIT cells resulting in their loss. These results provide a new role for HLA-G as a negative feedback loop by which B cells control MAIT cell responses to antigens.

Transcriptomics of type 2 diabetic and healthy human neutrophils

OBJECTIVES

Chronic inflammatory diseases, including diabetes and cardiovascular disease, are heterogeneous and often co-morbid, with increasing global prevalence. Uncontrolled type 2 diabetes (T2D) can result in severe inflammatory complications. As neutrophils are essential to normal and aberrant inflammation, we conducted RNA-seq transcriptomic analyses to investigate the association between neutrophil gene expression and T2D phenotype. As specialized pro-resolving lipid mediators (SPM) act to resolve inflammation, we further surveyed the impact of neutrophil receptor binding SPM resolvin E1 (RvE1) on isolated diabetic and healthy neutrophils.

METHODS

Cell isolation and RNA-seq analysis of neutrophils from N = 11 T2D and N = 7 healthy individuals with available clinical data was conducted. Additionally, cultured neutrophils (N = 3 T2D, N = 3 healthy) were perturbed with increasing RvE1 doses (0 nM, 1 nM, 10 nM, or 100 nM) prior to RNA-seq. Data was evaluated through a bioinformatics pipeline including pathway analysis and post hoc false discovery rate (FDR)-correction.

RESULTS

We observed significant differential expression of 50 genes between T2D and healthy neutrophils (p < 0.05), including decreased T2D gene expression in inflammatory- and lipid-related genes SLC9A4, NECTIN2, and PLPP3 (p < 0.003). RvE1 treatment induced dose-dependent differential gene expression (uncorrected p < 0.05) across groups, including 59 healthy and 216 T2D neutrophil genes. Comparing T2D to healthy neutrophils, 1097 genes were differentially expressed across RvE1 doses, including two significant genes, LILRB5 and AKR1C1, involved in inflammation (p < 0.05).

CONCLUSIONS

The neutrophil transcriptomic database revealed novel chronic inflammatory- and lipid-related genes that were differentially expressed between T2D cells when compared to controls, and cells responded to RvE1 dose-dependently by gene expression changes. Unraveling the mechanisms regulating abnormalities in diabetic neutrophil responses could lead to better diagnostics and therapeutics targeting inflammation and inflammation resolution.

Autophagic Regulation of p62 is Critical for Cancer Therapy

Sequestosome1 (p62/SQSTM 1) is a multidomain protein that interacts with the autophagy machinery as a key adaptor of target cargo. It interacts with phagophores through the LC3-interacting (LIR) domain and with the ubiquitinated protein aggregates through the ubiquitin-associated domain (UBA) domain. It sequesters the target cargo into inclusion bodies by its PB1 domain. This protein is further the central hub that interacts with several key signaling proteins. Emerging evidence implicates p62 in the induction of multiple cellular oncogenic transformations. Indeed, p62 upregulation and/or reduced degradation have been implicated in tumor formation, cancer promotion as well as in resistance to therapy. It has been established that the process of autophagy regulates the levels of p62. Autophagy-dependent apoptotic activity of p62 is recently being reported. It is evident that p62 plays a critical role in both autophagy and apoptosis. Therefore in this review we discuss the role of p62 in autophagy, apoptosis and cancer through its different domains and outline the importance of modulating cellular levels of p62 in cancer therapeutics.