Search

Background: Progress in research on expression profiles in osteoarthritis (OA) has been limited to individual tissues within the joint, such as the synovium, cartilage, or meniscus. This study aimed to comprehensively analyze the common gene expression characteristics of various structures in OA and construct a diagnostic model. Methods: Three datasets were selected: synovium, meniscus, and knee joint cartilage. Modular clustering and differential analysis of genes were used for further functional analyses and the construction of protein networks. Signature genes with the highest diagnostic potential were identified and verified using external gene datasets. The expression of these genes was validated in clinical samples by Real-time (RT)-qPCR and immunohistochemistry (IHC) staining. This study investigated the status of immune cells in OA by examining their infiltration. Results: The merged OA dataset included 438 DEGs clustered into seven modules using WGCNA. The intersection of these DEGs with WGCNA modules identified 190 genes. Using Least Absolute Shrinkage and Selection Operator (LASSO) and Random Forest algorithms, nine signature genes were identified (CDADC1, PPFIBP1, ENO2, NOM1, SLC25A14, METTL2A, LINC01089, L3HYPDH, NPHP3), each demonstrating substantial diagnostic potential (areas under the curve from 0.701 to 0.925). Furthermore, dysregulation of various immune cells has also been observed. Conclusion: CDADC1, PPFIBP1, ENO2, NOM1, SLC25A14, METTL2A, LINC01089, L3HYPDH, NPHP3 demonstrated significant diagnostic efficacy in OA and are involved in immune cell infiltration. [ABSTRACT FROM AUTHOR]

Jianmin Gong,1,2,* Bo Shi,1,3,* Ping Yang,1 Adeel Khan,4 Tao Xiong,2 Zhiyang Li1 1Department of Clinical Laboratory, Nanjing Drum Tower Hospital, Nanjing Drum Tower Hospital Clinical College of Nanjing Medical University, Nanjing, 210000, People’s Republic of China; 2College of Life Science, Yangtze University, Jingzhou, 434025, People’s Republic of China; 3Department of Clinical Laboratory, Nanjing Jiangning Hospital of Chinese Medicine (CM), Nanjing, 211100, People’s Republic of China; 4Department of Biotechnology, University of Science and Technology Bannu, Bannu, 28100, Islamic Republic of Pakistan*These authors contributed equally to this workCorrespondence: Tao Xiong, College of Life Science, Yangtze University, Jingzhou, 434025, People’s Republic of China, Tel +8613872387410, Email xiongtao@yangtzeu.edu.cn Zhiyang Li, Department of Clinical Laboratory, Nanjing Drum Tower Hospital, Nanjing Drum Tower Hospital Clinical College of Nanjing Medical University, Nanjing, 210008, People’s Republic of China, Tel +8618951703765, Email lizhiyang@nju.edu.cnBackground: Hypertrophic cardiomyopathy (HCM) is a dominantly inherited disease associated with sudden immune cell associations that remain unclear. The aim of this study was to comprehensively screen candidate markers associated with HCM and immune cells and explore potential pathogenic pathways.Methods: First, download the GSE32453 dataset to identify differentially expressed genes (DEGs) and perform Gene Ontology and pathway enrichment analysis using DAVID and GSEA. Next, construct protein-protein interaction (PPI) networks using String and Cytoscape to identify hub genes. Afterward, use CIBERSORT to determine the proportion of immune cells attributed to key genes in HCM and conduct ROC analysis based on the external dataset GSE36961 to evaluate their diagnostic value. Finally, validate the expression of key genes in the hypertrophic cardiomyocyte model through qRT-PCR using data from the HPA database.Results: Comprehensive analysis revealed that there were 254 upregulated genes and 181 downregulated genes in HCM. The enrichment study underscored pathways of inflammatory signaling, including MAPK and PI3K-Akt pathways. Pathways abundant in genes associated with HCM encompassed myocardial contraction and NADH dehydrogenase activity. Additionally, the analysis of immune infiltration revealed a notable increase in macrophages, NK cells, and monocytes in the HCM group, showing statistically significant variances in CD4 memory resting T cell infiltration when compared to the healthy control group. Within the validation dataset GSE36961, the Area Under the Curve (AUC) scores for eight crucial genes (FOS, CD86, CD68, BDNF, PIK3R1, PLEK, RAC2, CCL2) each exceeded 0.8. The HPA database revealed the positioning traits and paths of these eight crucial genes in smooth muscle cells, myocardial cells, and fibroblasts. The outcomes of the qRT-PCR were aligned with the sequencing findings.Conclusion: Bioinformatics analysis unveiled pivotal genes, pathways, and immune involvement, illuminating the molecular underpinnings of HCM. These findings suggest promising therapeutic targets for clinical applications.Keywords: hypertrophic cardiomyopathy, bioinformatic, immune cells infiltration, biomarker

Background: Systemic lupus erythematosus (SLE) is a multi-organ chronic autoimmune disease. Inflammatory bowel disease (IBD) is a common chronic inflammatory disease of the gastrointestinal tract. Previous studies have shown that SLE and IBD share common pathogenic pathways and genetic susceptibility, but the specific pathogenic mechanisms remain unclear. Methods: The datasets of SLE and IBD were downloaded from the Gene Expression Omnibus (GEO). Differentially expressed genes (DEGs) were identified using the Limma package. Weighted gene coexpression network analysis (WGCNA) was used to determine co-expression modules related to SLE and IBD. Pathway enrichment was performed using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis for co-driver genes. Using the Least AbsoluteShrinkage and Selection Operator (Lasso) regressionand Support Vector Machine-Recursive Feature Elimination (SVM-RFE), common diagnostic markers for both diseases were further evaluated. Then, we utilizedthe CIBERSORT method to assess the abundance of immune cell infiltration. Finally,we used the single-cell analysis to obtain the location of common diagnostic markers. Results: 71 common driver genes were identified in the SLE and IBD cohorts based on the DEGs and module genes. KEGG and GO enrichment results showed that these genes were closely associated with positive regulation of programmed cell death and inflammatory responses. By using LASSO regression and SVM, five hub genes (KLRF1, GZMK, KLRB1, CD40LG, and IL-7R) were ultimately determined as common diagnostic markers for SLE and IBD. ROC curve analysis also showed good diagnostic performance. The outcomes of immune cell infiltration demonstrated that SLE and IBD shared almost identical immune infiltration patterns. Furthermore, the majority of the hub genes were commonly expressed in NK cells by single-cell analysis. Conclusion: This study demonstrates that SLE and IBD share common diagnostic markers and pathogenic pathways. In addition, SLE and IBD show similar immune cellinfiltration microenvironments which provides newperspectives for future treatment. [ABSTRACT FROM AUTHOR]

Objective: Emerging evidence has shown that gut diseases can regulate the development and function of the immune, metabolic, and nervous systems through dynamic bidirectional communication on the brain-gut axis. However, the specific mechanism of intestinal diseases and vascular dementia (VD) remains unclear. We designed this study especially, to further clarify the connection between VD and inflammatory bowel disease (IBD) from bioinformatics analyses. Methods: We downloaded Gene expression profiles for VD (GSE122063) and IBD (GSE47908, GSE179285) from the Gene Expression Omnibus (GEO) database. Then individual Gene Set Enrichment Analysis (GSEA) was used to confirm the connection between the two diseases respectively. The common differentially expressed genes (coDEGs) were identified, and the STRING database together with Cytoscape software were used to construct proteinprotein interaction (PPI) network and core functional modules. We identified the hub genes by using the Cytohubba plugin. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were applied to identify pathways of coDEGs and hub genes. Subsequently, receiver operating characteristic (ROC) analysis was used to identify the diagnostic ability of these hub genes, and a training dataset was used to verify the expression levels of the hub genes. An alternative single-sample gene set enrichment (ssGSEA) algorithm was used to analyze immune cell infiltration between coDEGs and immune cells. Finally, the correlation between hub genes and immune cells was analyzed. Results: We screened 167 coDEGs. The main articles of coDEGs enrichment analysis focused on immune function. 8 shared hub genes were identified, including PTPRC, ITGB2, CYBB, IL1B, TLR2, CASP1, IL10RA, and BTK. The functional categories of hub genes enrichment analysis were mainly involved in the regulation of immune function and neuroinflammatory response. Compared to the healthy controls, abnormal infiltration of immune cells was found in VD and IBD. We also found the correlation between 8 shared hub genes and immune cells. Conclusions: This study suggests that IBD may be a new risk factor for VD. The 8 hub genes may predict the IBD complicated with VD. Immune-related coDEGS may be related to their association, which requires further research to prove. [ABSTRACT FROM AUTHOR]

Objective: Inflammation plays an important role in the transformation of pulmonary nodules (PNs) from benign to malignant. Prediction of benignancy and malignancy of PNs is still lacking efficacy methods. Although Mayo or Brock model have been widely applied in clinical practices, their application conditions are limited. This study aims to construct a diagnostic model of PNs by machine learning using inflammation-related biological markers (IRBMs). Methods: Inflammatory related genes (IRGs) were first extracted from GSE135304 chip data. Then, differentially expressed genes (DEGs) and infiltrating immune cells were screened between malignant pulmonary nodules (MN) and benign pulmonary nodule (BN). Correlation analysis was performed on DEGs and infiltrating immune cells. Molecular modules of IRGs were identified through Consistency cluster analysis. Subsequently, IRBMs in IRGs modules were filtered through Weighted gene co-expression network analysis (WGCNA). An optimal diagnostic model was established using machine learning methods. Finally, external dataset GSE108375 was used to verify this result. Results: 4 hub IRGs and 3 immune cells showed significantly difference between MN and BN, C1 and C2 module, namely PRTN3, ELANE, NFKB1 and CTLA4, T cells CD4 naïve, NK cells activated and Monocytes. IRBMs were screened from black module and yellowgreen module through WGCNA analysis. The Support vector machines (SVM) was identified as the optimal model with the Area Under Curve (AUC) was 0.753. A nomogram was established based on 5 hub IRBMs, namely HS.137078, KLC3, C13ORF15, STOM and KCTD13. Finally, external dataset GSE108375 verified this result, with the AUC was 0.718. Conclusion: SVM model established by 5 hub IRBMs was able to effectively identify MN or BN. Accumulating inflammation and immune dysfunction were important to the transformation from BN to MN.

Chronic inflammation promotes the development of pancreatic ductal adenocarcinoma (PDAC) and PDAC-related inflammatory tumor microenvironment facilitates tumor growth and metastasis. Thus, we aimed to study the association between inflammatory response and prognosis in patients with PDAC. We conducted the whole transcriptomic sequencing using tissue samples collected from patients diagnosed with PDAC (n = 106) recruited from Shandong Cancer Hospital. We first constructed a prognostic signature using 15 inflammation-related genes in The Cancer Genome Atlas (TCGA) cohort (n = 177) and further validated it in an independent International Cancer Genome Consortium (ICGC) cohort (n = 90) and our in-house cohort. PDAC patients with a higher risk score had poorer overall survival (OS) (P < 0.001; HR, 3.02; 95% CI, 1.94-4.70). The association between the prognostic signature and OS remained significant in the multivariable Cox regression adjusting for age, sex, alcohol exposure, diabetes, and stage (P < 0.001; HR, 2.91; 95% CI, 1.73-4.89). This gene signature also robustly predicted prognosis in the ICGC cohort (P = 0.01; HR, 1.94; 95% CI, 1.14-3.30) and our cohort (P < 0.001; HR, 2.40; 95% CI, 1.45-3.97). Immune subtype C3 (inflammatory) was enriched and CD8+ T cells were higher in patients with a lower risk score (P < 0.05). Furthermore, PDAC patients with higher risk scores were more sensitive to chemotherapy, immunotherapy, and PARP inhibitors (P < 0.05). In sum, we identified a novel gene signature that was associated with inflammatory response for risk stratification, prognosis prediction, and therapy guidance in PDAC patients. Future studies are warranted to validate the clinical utility of the signature. [ABSTRACT FROM AUTHOR]

Inflammatory bowel disease (IBD) is a chronic inflammatory disorder of the gastrointestinal tract which is mediated by the inappropriate immune responses. This study was aimed to identify novel diagnostic biomarkers for diagnosis of IBD and explore the relationship between the diagnostic biomarkers and infiltrated immune cells. GSE38713, GSE53306, and GSE75214 downloaded from the Gene Expression Omnibus (GEO) database were split into training and testing sets. Differentially expressed genes (DEGs) were screened using the "limma" package. Gene Ontology (GO) and KEGG pathway enrichment analysis of DEGs were performed by clusterProfiler package. The LASSO regression and support vector machine recursive feature elimination (SVM-RFE) algorithms were conducted to identify novel diagnostic biomarkers. The receiver operating characteristic (ROC) curve was applied to evaluate the diagnostic value of the candidate biomarkers. The relationship of the candidate biomarkers and infiltrating immune cells in IBD were evaluated by CIBERSOTR. Quantitative Real-Time PCR (qRT-PCR) was applied to measure the expression level of the biomarkers in IBD. A total of 289 dysregulated genes were identified as DEGs in IBD. These DEGs were significantly enriched in chemokine signaling pathway and cytokine–cytokine receptor interaction. RHOU was identified as a critical diagnostic gene in IBD, which was confirmed using ROC curve and qRT-PCR assays. Immune cell infiltration analysis showed that RHOU was correlated with macrophages M2, dendritic cells resting, mast cells resting, T cells CD4 memory resting, macrophages M0, and mast cells activated. Our results imply that RHOU may be a potential diagnostic biomarker for IBD. [ABSTRACT FROM AUTHOR]

Abstract Background Pancreatic cancer (PC) stands out as one of the most formidable malignancies and exhibits an exceptionally unfavorable clinical prognosis due to the absence of well‐defined diagnostic indicators and its tendency to develop resistance to therapeutic interventions. The primary objective of this present study was to identify extracellular matrix (ECM)‐related hub genes (HGs) and their corresponding molecular signatures, with the intent of potentially utilizing them as biomarkers for diagnostic, prognostic, and therapeutic applications. Methods Three microarray datasets were sourced from the NCBI database to acquire upregulated differentially expressed genes (DEGs), while MatrisomeDB was employed for filtering ECM‐related genes. Subsequently, a protein–protein interaction (PPI) network was established using the STRING database. The created network was visually inspected through Cytoscape, and HGs were identified using the CytoHubba plugin tool. Furthermore, enrichment analysis, expression pattern analysis, clinicopathological correlation, survival analysis, immune cell infiltration analysis, and examination of chemical compounds were carried out using Enrichr, GEPIA2, ULCAN, Kaplan Meier plotter, TIMER2.0, and CTD web platforms, respectively. The diagnostic and prognostic significance of HGs was evaluated through the ROC curve analysis. Results Ten genes associated with ECM functions were identified as HGs among 131 DEGs obtained from microarray datasets. Notably, the expression of these HGs exhibited significantly (p 0.75) and prognostic (AUC > 0.64) purposes. Furthermore, we observed a positive correlation between immune cell infiltration and the expression of most HGs. Lastly, our study identified nine compounds with significant interaction profiles that could potentially act as effective chemical agents targeting the identified HGs. Conclusion Taken together, our findings suggest that COL1A1, KRT19, MMP1, COL11A1, SDC1, ITGA2, COL1A2, POSTN, FN1, and COL5A1 hold promise as innovative biomarkers for both the diagnosis and prognosis of PC, and they present as prospective targets for therapeutic interventions aimed at impeding the progression PC.

Abstract Background Current research on autophagy is mainly focused on intervertebral disc tissues and cells, while there is few on human peripheral blood sample. therefore, this study constructed a diagnostic model to identify autophagy-related markers of intervertebral disc degeneration (IVDD). Methods GSE150408 and GSE124272 datasets were acquired from the Gene Expression Omnibus database, and differential expression analysis was performed. The IVDD-autophagy genes were obtained using Weighted Gene Coexpression Network Analysis, and a diagnostic model was constructed and validated, followed by Gene Set Variation Analysis (GSVA) and Gene Set Enrichment Analysis (GSEA). Meanwhile, miRNA–gene and transcription factor–gene interaction networks were constructed. In addition, drug-gene interactions and target genes of methylprednisolone and glucosamine were analyzed. Results A total of 1,776 differentially expressed genes were identified between IVDD and control samples, and the composition of the four immune cell types was significantly different between the IVDD and control samples. The Meturquoise and Mebrown modules were significantly related to immune cells, with significant differences between the control and IVDD samples. A diagnostic model was constructed using five key IVDD-autophagy genes. The area under the curve values of the model in the training and validation datasets were 0.907 and 0.984, respectively. The enrichment scores of the two pathways were significantly different between the IVDD and healthy groups. Eight pathways in the IVDD and healthy groups had significant differences. A total of 16 miRNAs and 3 transcription factors were predicted to be of great value. In total, 84 significantly related drugs were screened for five key IVDD-autophagy genes in the diagnostic model, and three common autophagy-related target genes of methylprednisolone and glucosamine were predicted. Conclusion This study constructs a reliable autophagy-related diagnostic model that is strongly related to the immune microenvironment of IVD. Autophagy-related genes, including PHF23, RAB24, STAT3, TOMM5, and DNAJB9, may participate in IVDD pathogenesis. In addition, methylprednisolone and glucosamine may exert therapeutic effects on IVDD by targeting CTSD, VEGFA, and BAX genes through apoptosis, as well as the sphingolipid and AGE-RAGE signaling pathways in diabetic complications.

Background: Current research on autophagy is mainly focused on intervertebral disc tissues and cells, while there is few on human peripheral blood sample. therefore, this study constructed a diagnostic model to identify autophagy-related markers of intervertebral disc degeneration (IVDD). Methods: GSE150408 and GSE124272 datasets were acquired from the Gene Expression Omnibus database, and differential expression analysis was performed. The IVDD-autophagy genes were obtained using Weighted Gene Coexpression Network Analysis, and a diagnostic model was constructed and validated, followed by Gene Set Variation Analysis (GSVA) and Gene Set Enrichment Analysis (GSEA). Meanwhile, miRNA–gene and transcription factor–gene interaction networks were constructed. In addition, drug-gene interactions and target genes of methylprednisolone and glucosamine were analyzed. Results: A total of 1,776 differentially expressed genes were identified between IVDD and control samples, and the composition of the four immune cell types was significantly different between the IVDD and control samples. The Meturquoise and Mebrown modules were significantly related to immune cells, with significant differences between the control and IVDD samples. A diagnostic model was constructed using five key IVDD-autophagy genes. The area under the curve values of the model in the training and validation datasets were 0.907 and 0.984, respectively. The enrichment scores of the two pathways were significantly different between the IVDD and healthy groups. Eight pathways in the IVDD and healthy groups had significant differences. A total of 16 miRNAs and 3 transcription factors were predicted to be of great value. In total, 84 significantly related drugs were screened for five key IVDD-autophagy genes in the diagnostic model, and three common autophagy-related target genes of methylprednisolone and glucosamine were predicted. Conclusion: This study constructs a reliable autophagy-related diagnostic model that is strongly related to the immune microenvironment of IVD. Autophagy-related genes, including PHF23, RAB24, STAT3, TOMM5, and DNAJB9, may participate in IVDD pathogenesis. In addition, methylprednisolone and glucosamine may exert therapeutic effects on IVDD by targeting CTSD, VEGFA, and BAX genes through apoptosis, as well as the sphingolipid and AGE-RAGE signaling pathways in diabetic complications. Highlights: The autophagy-related diagnostic model closely associated with immunity was constructed based on five key autophagy-related genes in IVDD patients' peripheral blood sample, which could diagnose of IVDD. Methylprednisolone and glucosamine, empirical used clinical therapeutic drugs for IVDD, might exert therapeutic effect by targeting CTSD, VEGFA, and BAX through apoptosis, sphingolipid signaling pathway, and AGE-RAGE signaling pathway in diabetic complications. This study provides new insights for the diagnosis and therapy of IVDD. [ABSTRACT FROM AUTHOR]

Severe Fever with thrombocytopenia syndrome (SFTS) is a highly fatal viral infectious disease that poses a significant threat to public health. Currently, the phase and pathogenesis of SFTS are not well understood, and there are no specific vaccines or effective treatment available. Therefore, it is crucial to identify biomarkers for diagnosing acute SFTS, which has a high mortality rate. In this study, we conducted differentially expressed genes (DEGs) analysis and WGCNA module analysis on the GSE144358 dataset, comparing the acute phase of SFTSV-infected patients with healthy individuals. Through the LASSO–Cox and random forest algorithms, a total of 2128 genes were analyzed, leading to the identification of four genes: ADIPOR1, CENPO, E2F2, and H2AC17. The GSEA analysis of these four genes demonstrated a significant correlation with immune cell function and cell cycle, aligning with the functional enrichment findings of DEGs. Furthermore, we also utilized CIBERSORT to analyze the immune cell infiltration and its correlation with characteristic genes. The results indicate that the combination of ADIPOR1, CENPO, E2F2, and H2AC17 genes has the potential as characteristic genes for diagnosing and studying the acute phase of SFTS virus (SFTSV) infection. [ABSTRACT FROM AUTHOR]

Background: Bladder cancer (BLCA) is a common malignant tumor of urinary system and prognostic biomarkers are needed for better clinical decision-making and patient management. Cancer stem cells (CSCs) are involved in carcinogenesis, development, metastasis and recurrence of BLCA. This study explored the prognostic and predictive value of CSCs-related genes and laid the groundwork for precision treatment development in BLCA. Methods: The mRNA data and corresponding clinical information obtained from TCGA-BLCA cohort was used to discover biomarkers and develop CSCs-related prognostic model, which was further validated in GSE32548 and GSE32894 datasets. In addition, the association between CSCs-related risk score and therapeutic efficacy was analyzed to explore the potential predictive value of the prognostic model. Results: We identified four CSCs-related subtypes and 900 differentially expressed genes (DEGs) among subtypes. Then the CSCs-related prognostic model was built based on 16 CSCs-related DEGs with the most significant prognostic value. Patients in the low-risk group had better overall survival than those in high-risk group (P

Abstract Objective The aim of this study was to construct a model used for the accurate diagnosis of Atopic dermatitis (AD) using pyroptosis related biological markers (PRBMs) through the methods of machine learning. Method The pyroptosis related genes (PRGs) were acquired from molecular signatures database (MSigDB). The chip data of GSE120721, GSE6012, GSE32924, and GSE153007 were downloaded from gene expression omnibus (GEO) database. The data of GSE120721 and GSE6012 were combined as the training group, while the others were served as the testing groups. Subsequently, the expression of PRGs was extracted from the training group and differentially expressed analysis was conducted. CIBERSORT algorithm calculated the immune cells infiltration and differentially expressed analysis was conducted. Consistent cluster analysis divided AD patients into different modules according to the expression levels of PRGs. Then, weighted correlation network analysis (WGCNA) screened the key module. For the key module, we used Random forest (RF), support vector machines (SVM), Extreme Gradient Boosting (XGB), and generalized linear model (GLM) to construct diagnostic models. For the five PRBMs with the highest model importance, we built a nomogram. Finally, the results of the model were validated using GSE32924, and GSE153007 datasets. Results Nine PRGs were significant differences in normal humans and AD patients. Immune cells infiltration showed that the activated CD4+ memory T cells and Dendritic cells (DCs) were significantly higher in AD patients than normal humans, while the activated natural killer (NK) cells and the resting mast cells were significantly lower in AD patients than normal humans. Consistent cluster analysis divided the expressing matrix into 2 modules. Subsequently, WGCNA analysis showed that the turquoise module had a significant difference and high correlation coefficient. Then, the machine model was constructed and the results showed that the XGB model was the optimal model. The nomogram was constructed by using HDAC1, GPALPP1, LGALS3, SLC29A1, and RWDD3 five PRBMs. Finally, the datasets GSE32924 and GSE153007 verified the reliability of this result. Conclusions The XGB model based on five PRBMs can be used for the accurate diagnosis of AD patients.

Abstract Background Dermatomyositis (DM) is an acquired autoimmune disease that can cause damage to various organs, including the heart muscle. However, the mechanisms underlying myocardial injury in DM are not yet fully understood. Methods In this study, we utilized publicly available datasets from the Gene Expression Omnibus (GEO) database to identify hub-genes that are enriched in the immune system process in DM and myocarditis. Weighted gene co-expression network analysis (WGCNA), differentially expressed genes (DEGs) analysis, protein–protein interaction (PPI), and gene ontology (GO) analysis were employed to identify these hub-genes. We then used the CIBERSORT method to analyze immune cell infiltration in skeletal muscle specimens of DM and myocardium specimens of myocarditis respectively. Correlation analysis was performed to investigate the relationship between key genes and infiltrating immune cells. Finally, we predicted regulatory miRNAs of hub-genes through miRNet and validated their expression in online datasets and clinical samples. Results Using integrated bioinformatics analysis, we identified 10 and 5 hub-genes that were enriched in the immune system process in the database of DM and myocarditis respectively. The subsequent intersections between hub-genes were IFIT3, OAS3, ISG15, and RSAD2. We found M2 macrophages increased in DM and myocarditis compared to the healthy control, associating with the expression of IFIT3, OAS3, ISG15, and RSAD2 in DM and myocarditis positively. Gene function enrichment analysis (GSEA) showed that IFIT3, OAS3, ISG15, and RSAD2 were mainly enriched in type I interferon (IFN) signaling pathway, cellular response to type I interferon, and response to type I interferon. Finally, we verified that the expression of miR-146a-5p was significantly higher in the DM with myocardial injury than those without myocardial injury (p = 0.0009). Conclusion Our findings suggest that IFIT3, OAS3, ISG15, and RSAD2 may play crucial roles in the underlying mechanism of myocardial injury in DM. Serum miR-146a-5p could be a potential biomarker for myocardial injury in DM.

Objective: To identify immune-related genes （IRGs） and explore role of immune cells infiltration in systemic sclerosis-related pulmonary arterial hypertension （SSc-PAH） based on bioinformatics analysis. Methods: GSE22356 dataset was downloaded from GEO database and conducted with differentially expressed genes （DEGs） analysis by R software. Human immunity gene dataset was downloaded from ImmPort database, which were intersected with DEGs to obtain differentially expressed immune-related genes （DEIRGs）. GO and KEGG enrichment analysis were performed with DEIRGs. Protein-protein interaction （PPI） network of DEIRGs was built by Cytoscape, and key immune genes were screened according to Degree value. CIBERSORT algorithm was used to evaluate immune cells infiltration between SSc and SSc-PAH patients. Correlation analysis between key DEIRGs and infiltrating immune cells was performed. Results: A total of 56 DEIRGs were detected. Enrichment function demonstrated that DEIRGs were significant in signal transduction, immune response and chemotaxis. Key DEIRGs assessed by PPI network and degree score were FCGR3B, CD28, LYN, LCK. Immune cells infiltration findings indicated that compared with SSc, SSc-PAH patients contained a higher proportion of monocytes and neutrophils, while a lower proportion of naive CD4 T cells, inactive CD4 memory T cells and γδ T cells. Correlation analysis results demonstrated that CD28 and LCK were positively correlated with γδ T cells and negatively correlated with monocytes. FCGR3B was positively correlated with neutrophils and negatively correlated with naive CD4 T cells. LYN was positively correlated with monocytes and negatively correlated with γδ T cells. Conclusion: Key immune-related genes and differences in immune cells infiltration play an essential role in occurrence and progression of SSc-PAH. [ABSTRACT FROM AUTHOR]

Background: Cerebrovascular disease (CVD) related to atherosclerosis and Parkinson's disease (PD) are two prevalent neurological disorders. They share common risk factors and frequently occur together. The aim of this study is to investigate the association between atherosclerosis and PD using genetic databases to gain a comprehensive understanding of underlying biological mechanisms. Methods: The gene expression profiles of atherosclerosis (GSE28829 and GSE100927) and PD (GSE7621 and GSE49036) were downloaded from the Gene Expression Omnibus (GEO) database. After identifying the common differentially expressed genes (DEGs) for these two disorders, we constructed proteinprotein interaction (PPI) networks and functional modules, and further identified hub genes using Least Absolute Shrinkage and Selection Operator (LASSO) regression. The diagnostic effectiveness of these hub genes was evaluated using Receiver Operator Characteristic Curve (ROC) analysis. Furthermore, we used single sample gene set enrichment analysis (ssGSEA) to analyze immune cell infiltration and explored the association of the identified hub genes with infiltrating immune cells through Spearman's rank correlation analysis in R software. Results: A total of 50 shared DEGs, with 36 up-regulated and 14 down-regulated genes, were identified through the intersection of DEGs of atherosclerosis and PD. Using LASSO regression, we identified six hub genes, namely C1QB, CD53, LY96, P2RX7, C3, and TNFSF13B, in the lambda.min model, and CD14, C1QB, CD53, P2RX7, C3, and TNFSF13B in the lambda.1se model. ROC analysis confirmed that both models had good diagnostic efficiency for atherosclerosis datasets GSE28829 (lambda.min AUC = 0.99, lambda.1se AUC = 0.986) and GSE100927 (lambda.min AUC = 0.922, lambda.1se AUC = 0.933), as well as for PD datasets GSE7621 (lambda.min AUC = 0.924, lambda.1se AUC = 0.944) and GSE49036 (lambda.min AUC = 0.894, lambda.1se AUC = 0.881). Furthermore, we found that activated B cells, effector memory CD8 + T cells, and macrophages were the shared correlated types of immune cells in both atherosclerosis and PD. Conclusion: This study provided new sights into shared molecular mechanisms between these two disorders. These common hub genes and infiltrating immune cells offer promising clues for further experimental studies to explore the common pathogenesis of these disorders. [ABSTRACT FROM AUTHOR]

Background: It is well known that patients with systemic lupus erythematosus (SLE) had a high risk of venous thromboembolism (VTE). This study aimed to identify the crosstalk genes between SLE and VTE and explored their clinical value and molecular mechanism initially. Methods: We downloaded microarray datasets of SLE and VTE from the Gene Expression Omnibus (GEO) dataset. Differential expression analysis was applied to identify the crosstalk genes (CGs). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed on the shared genes. The shared diagnostic biomarkers of the two diseases were further screened from CGs using least absolute shrinkage and selection operator (Lasso) regression. Two risk scores for SLE and VTE were constructed separately to predict the likelihood of illness according to the diagnostic biomarkers using a logical regression algorithm. The immune infiltration levels of SEL and VTE were estimated via the CIBERSORT algorithm and the relationship of CGs with immune cell infiltration was investigated. Finally, we explored potential phenotype subgroups in SLE and VTE based on the expression level of CGs through the consensus clustering method and studied immune cell infiltration in different subtypes. Result: A total of 171 CGs were obtained by the intersection of differentially expressed genes (DEGs) between SLE and VTE cohorts. The functional enrichment shown these CGs were mainly related to immune pathways. After screening by lasso regression, we found that three hub CGs (RSAD2, HSP90AB1, and FPR2) were the optimal shared diagnostic biomarkers for SLE and VTE. Based on the expression level of RSAD2 and HSP90AB1, two risk prediction models for SLE and VTE were built by multifactor logistic regression and systemically validated in internal and external validation datasets. The immune infiltration results revealed that CGs were highly correlated with multiple infiltrated immunocytes. Consensus clustering was used to respectively regroup SLE and VTE patients into C1 and C2 clusters based on the CGs expression profile. The levels of immune cell infiltration and immune activation were higher in C1 than in C2 subtypes. Conclusion: In our study, we further screen out diagnostic biomarkers from crosstalk genes SLE and VTE and built two risk scores. Our findings reveal a close relationship between CGs and the immune microenvironment of diseases. This provides clues for further exploring the common mechanism and interaction between the two diseases. [ABSTRACT FROM AUTHOR]

Objective: Evidences show that there may be a link between SLE and COVID-19. The purpose of this study is to screen out the diagnostic biomarkers of systemic lupus erythematosus (SLE) with COVID-19 and explore the possible related mechanisms by the bioinformatics approach. Methods: SLE and COVID-19 datasets were extracted separately from the NCBI Gene Expression Omnibus (GEO) database. The limma package in R was used to obtain the differential genes (DEGs). The protein interaction network information (PPI) and core functional modules were constructed in the STRING database using Cytoscape software. The hub genes were identified by the Cytohubba plugin, and TF-gene together with TF-miRNA regulatory networks were constructed via utilizing the Networkanalyst platform. Subsequently, we generated subject operating characteristic curves (ROC) to verify the diagnostic capabilities of these hub genes to predict the risk of SLE with COVID-19 infection. Finally, a single-sample gene set enrichment (ssGSEA) algorithm was used to analyze immune cell infiltration. Results: A total of 6 common hub genes (CDC6, PLCG1, KIF15, LCK, CDC25C, and RASGRP1) were identified with high diagnostic validity. These gene functional enrichments were mainly involved in cell cycle, and inflammationrelated pathways. Compared to the healthy controls, abnormal infiltration of immune cells was found in SLE and COVID-19, and the proportion of immune cells linked to the 6 hub genes. Conclusion: Our research logically identified 6 candidate hub genes that could predict SLE complicated with COVID-19. This work provides a foothold for further study of potential pathogenesis in SLE and COVID-19. [ABSTRACT FROM AUTHOR]

Objective: Secondary hyperparathyroidism (SHPT) is a frequent complication of chronic kidney disease (CKD) associated with morbidity and mortality. This study aims to identify potential biomarkers that may be used to predict the progression of SHPT and to elucidate the molecular mechanisms of SHPT pathogenesis at the transcriptome level. Methods: We analyzed differentially expressed genes (DEGs) between diffuse and nodular parathyroid hyperplasia of SHPT patients from the GSE75886 dataset, and then verified DEG levels with the GSE83421 data file of primary hyperparathyroidism (PHPT) patients. Candidate gene sets were selected by machine learning screens of differential genes and immune cell infiltration was explored with the CIBERSORT algorithm. RcisTarget was used to predict transcription factors, and Cytoscape was used to construct a lncRNA-miRNA-mRNA network to identify possible molecular mechanisms. Immunohistochemistry (IHC) staining and quantitative real-time polymerase chain reaction (qRT-PCR) were used to verify the expression of screened genes in parathyroid tissues of SHPT patients and animal models. Results: A total of 614 DEGs in GSE75886 were obtained as candidate gene sets for further analysis. Five key genes (USP12, CIDEA, PCOLCE2, CAPZA1, and ACCN2) had significant expression differences between groups and were screened with the best ranking in the machine learning process. These genes were shown to be closely related to immune cell infiltration levels and play important roles in the immune microenvironment. Transcription factor ZBTB6 was identified as the master regulator, alongside multiple other transcription factors. Combined with qPCR and IHC assay of hyperplastic parathyroid tissues from SHPT patients and rats confirm differential expression of USP12, CIDEA, PCOLCE2, CAPZA1, and ACCN2, suggesting that they may play important roles in the proliferation and progression of SHPT. Conclusion: USP12, CIDEA, PCOLCE2, CAPZA1, and ACCN2 have great potential both as biomarkers and as therapeutic targets in the proliferation of SHPT. These findings suggest novel potential targets and future directions for SHPT research. [ABSTRACT FROM AUTHOR]

Abstract Background The incidence and mortality of gastric cancer ranks fifth and fourth worldwide among all malignancies, respectively. Accumulating evidences have revealed the close relationship between mitochondrial dysfunction and the initiation and progression of stomach cancer. However, rare prognostic models for mitochondrial-related gene risk have been built up in stomach cancer. Methods In current study, the expression and prognostic value of mitochondrial-related genes in stomach adenocarcinoma (STAD) patients were systematically analyzed to establish a mitochondrial-related risk model based on available TCGA and GEO databases. The tumor microenvironment (TME), immune cell infiltration, tumor mutation burden, and drug sensitivity of gastric adenocarcinoma patients were also investigated using R language, GraphPad Prism 8 and online databases. Results We established a mitochondrial-related risk prognostic model including NOX4, ALDH3A2, FKBP10 and MAOA and validated its predictive power. This risk model indicated that the immune cell infiltration in high-risk group was significantly different from that in the low-risk group. Besides, the risk score was closely related to TME signature genes and immune checkpoint molecules, suggesting that the immunosuppressive tumor microenvironment might lead to poor prognosis in high-risk groups. Moreover, TIDE analysis demonstrated that combined analysis of risk score and immune score, or stromal score, or microsatellite status could more effectively predict the benefit of immunotherapy in STAD patients with different stratifications. Finally, rapamycin, PD-0325901 and dasatinib were found to be more effective for patients in the high-risk group, whereas AZD7762, CEP-701 and methotrexate were predicted to be more effective for patients in the low-risk group. Conclusions Our results suggest that the mitochondrial-related risk model could be a reliable prognostic biomarker for personalized treatment of STAD patients.

Background: Intervertebral disc degeneration (IDD) is the leading cause of lower back pain, and an overall understanding of the molecular mechanisms related to IDD is still lacking. The purpose of this study was to explore gene signatures and immune cell infiltration related to IDD via bioinformatics analysis.Methods: A total of five expression profiles of mRNA and non-coding RNA were downloaded from the Gene Expression Omnibus (GEO) database. The potentially involved lncRNA/circRNA–miRNA–mRNA networks and protein-protein interaction networks were constructed by miRNet, circBank, STRING, and the Cytoscape database. Gene ontology, Kyoto Encyclopaedia of Genes and Genomes Analysis, Gene Set Enrichment Analysis, Gene Set Variation Analysis, Immune Infiltration Analysis, and Drug-Gene Interaction were used to analyse the top 20 hub genes. RT-qPCR was conducted to confirm the 12 differential expressions of genes both in the nucleus pulposus and annulus fibrosus tissuesResults: There were 346 differentially expressed mRNAs, 12 differentially expressed miRNAs, 883 differentially expressed lncRNAs, and 916 differentially expressed circRNAs in the GEO database. Functional and enrichment analyses revealed hub genes associated with platelet activation, immune responses, focal adhesion, and PI3K-Akt signalling. The apoptotic pathway, the reactive oxygen species pathway, and oxidative phosphorylation play an essential role in IDD. Immune infiltration analysis demonstrated that the Treg cells had significant infiltration, and three levels of immune cells, including dendritic cells, Th2 cells, and tumour-infiltrating lymphocytes, were inhibited in IDD. Drug-gene interaction analysis showed that COL1A1 and COL1A2 were targeted by collagenase clostridium histolyticum, ocriplasmin, and PDGFRA was targeted by 66 drugs or molecular compounds. Finally, 24 cases of IDD tissues and 12 cases of normal disc tissues were collected, and the results of RT-qPCR were consistent with the bioinformatics results.Conclusion: Our data indicated that the 20 hub genes and immune cell infiltration were involved in the pathological process of IDD. In addition, the PDGFRA and two potential drugs were found to be significant in IDD development.

Objective: The aim of this study was to construct a model used for the accurate diagnosis of Atopic dermatitis (AD) using pyroptosis related biological markers (PRBMs) through the methods of machine learning. Method: The pyroptosis related genes (PRGs) were acquired from molecular signatures database (MSigDB). The chip data of GSE120721, GSE6012, GSE32924, and GSE153007 were downloaded from gene expression omnibus (GEO) database. The data of GSE120721 and GSE6012 were combined as the training group, while the others were served as the testing groups. Subsequently, the expression of PRGs was extracted from the training group and differentially expressed analysis was conducted. CIBERSORT algorithm calculated the immune cells infiltration and differentially expressed analysis was conducted. Consistent cluster analysis divided AD patients into different modules according to the expression levels of PRGs. Then, weighted correlation network analysis (WGCNA) screened the key module. For the key module, we used Random forest (RF), support vector machines (SVM), Extreme Gradient Boosting (XGB), and generalized linear model (GLM) to construct diagnostic models. For the five PRBMs with the highest model importance, we built a nomogram. Finally, the results of the model were validated using GSE32924, and GSE153007 datasets. Results: Nine PRGs were significant differences in normal humans and AD patients. Immune cells infiltration showed that the activated CD4+ memory T cells and Dendritic cells (DCs) were significantly higher in AD patients than normal humans, while the activated natural killer (NK) cells and the resting mast cells were significantly lower in AD patients than normal humans. Consistent cluster analysis divided the expressing matrix into 2 modules. Subsequently, WGCNA analysis showed that the turquoise module had a significant difference and high correlation coefficient. Then, the machine model was constructed and the results showed that the XGB model was the optimal model. The nomogram was constructed by using HDAC1, GPALPP1, LGALS3, SLC29A1, and RWDD3 five PRBMs. Finally, the datasets GSE32924 and GSE153007 verified the reliability of this result. Conclusions: The XGB model based on five PRBMs can be used for the accurate diagnosis of AD patients. [ABSTRACT FROM AUTHOR]

Introduction: The similarity between ankylosing spondylitis (AS) and ulcerative colitis (UC) in incidence rate and pathogenesis has been revealed. But the common pathogenesis that explains the relationship between AS and UC is still lacked, and the related genetic research is limited. We purposed to explore shared biomarkers and pathways of AS and UC through integrated bioinformatics. Methods: Gene expression data of AS and UC were obtained in the GEO database. We applied weighted gene co-expression network analysis (WGCNA) to identify AS-related and UC-related co-expression gene modules. Subsequently, machine learning algorithm was used to further screen hub genes. We validated the expression level and diagnostic efficiency of the shared diagnostic gene of AS and UC in external datasets. Gene set enrichment analysis (GSEA) was applied to analyze pathway-level changes between disease group and normal group. Finally, we analyzed the relationship between hub biomarker and immune microenvironment by using the CIBERSORT deconvolution algorithm. Results: 203 genes were obtained by overlapping AS-related gene module and UC-related gene module. Through SVM-RFE algorithm, 19 hub diagnostic genes were selected for AS in GSE25101 and 6 hub diagnostic genes were selected for UC in GSE94648. KCNJ15 was obtained as a common diagnostic gene of AS and UC. The expression of KCNJ15 was validated in independent datasets, and the results showed that KCNJ15 were similarly upregulated in AS samples and UC samples. Besides, ROC analysis also revealed that KCNJ15 had good diagnostic efficacy. The GSEA analysis revealed that oxidative phosphorylation pathway was the shared pathway of AS and UC. In addition, CIBERSORT results revealed the correlation between KCNJ15 gene and immune microenvironment in AS and UC. Conclusion: We have explored a common diagnostic gene KCNJ15 and a shared oxidative phosphorylation pathway of AS and UC through integrated bioinformatics, which may provide a potential diagnostic biomarker and novel insight for studying the mechanism of AS-related UC. [ABSTRACT FROM AUTHOR]

Background: Hepatocellular carcinoma (HCC) is a primary malignant tumor of the liver, and pyroptosis has been identified as a novel cellular program that plays a role in numerous diseases including cancer. However, the functional role of pyroptosis in HCC remains unclear. The purpose of this study is to explore the relationship between the two found hub genes and provide targets for clinical treatment. Methods: The Cancer Genome Atlas (TCGA) database was used to collect the gene data and clinically-related information of patients with HCC. After the differentially expressed genes (DEGs) were identified, they were intersected with the genes related to pyroptosis, and a risk prediction model was established to predict the overall survival (OS). Subsequently, drug sensitivity analysis, Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), Gene Set Enrichment Analysis (GSEA), and Gene Set Variation Analysis (GSVA) was used to analyze the biological characteristics of the DEGs. Different immune cell infiltration and related pathways were analyzed, and hub genes were identified by protein-protein interaction (PPI). Finally, the expression of hub genes was verified by real-time quantitative PCR (qRT- PCR) and immunohistochemistry. Results: We conducted a comprehensive bioinformatics analysis to investigate the molecular mechanisms of pyroptosis in hepatocellular carcinoma (HCC). A total of 8,958 differentially expressed genes were identified, and 37 differentially expressed genes were associated with pyroptosis through intersection. Moreover, we developed an OS model with excellent predictive ability and discovered the differences in biological function, drug sensitivity, and immune microenvironment between high-risk and low-risk groups. Through enrichment analysis, we found that the differentially expressed genes are related to various biological processes. Then, 10 hub genes were identified from protein-protein interaction networks. Finally, midkine (MDK) was screened from the 10 hub genes and further verified by PCR and immunohistochemistry, which revealed its high expression in HCC. Conclusion: We have developed a reliable and consistent predictive model based on the identification of potential hub genes, which can be used to accurately forecast the prognosis of patients, thus providing direction for further clinical research and treatment. [ABSTRACT FROM AUTHOR]

Background: The incidence and mortality of gastric cancer ranks fifth and fourth worldwide among all malignancies, respectively. Accumulating evidences have revealed the close relationship between mitochondrial dysfunction and the initiation and progression of stomach cancer. However, rare prognostic models for mitochondrial-related gene risk have been built up in stomach cancer. Methods: In current study, the expression and prognostic value of mitochondrial-related genes in stomach adenocarcinoma (STAD) patients were systematically analyzed to establish a mitochondrial-related risk model based on available TCGA and GEO databases. The tumor microenvironment (TME), immune cell infiltration, tumor mutation burden, and drug sensitivity of gastric adenocarcinoma patients were also investigated using R language, GraphPad Prism 8 and online databases. Results: We established a mitochondrial-related risk prognostic model including NOX4, ALDH3A2, FKBP10 and MAOA and validated its predictive power. This risk model indicated that the immune cell infiltration in high-risk group was significantly different from that in the low-risk group. Besides, the risk score was closely related to TME signature genes and immune checkpoint molecules, suggesting that the immunosuppressive tumor microenvironment might lead to poor prognosis in high-risk groups. Moreover, TIDE analysis demonstrated that combined analysis of risk score and immune score, or stromal score, or microsatellite status could more effectively predict the benefit of immunotherapy in STAD patients with different stratifications. Finally, rapamycin, PD-0325901 and dasatinib were found to be more effective for patients in the high-risk group, whereas AZD7762, CEP-701 and methotrexate were predicted to be more effective for patients in the low-risk group. Conclusions: Our results suggest that the mitochondrial-related risk model could be a reliable prognostic biomarker for personalized treatment of STAD patients. [ABSTRACT FROM AUTHOR]

Background: The study aimed to identify core biomarkers related to diagnosis and immune microenvironment regulation and explore the immune molecular mechanism of diabetic nephropathy (DN) through bioinformatics analysis. Methods: GSE30529, GSE99325, and GSE104954 were merged with removing batch effects, and different expression genes (DEGs) were screened at a criterion |log2FC| >0.5 and adjusted P <0.05. KEGG, GO, and GSEA analyses were performed. Hub genes were screened by conducting PPI networks and calculating node genes using five algorithms with CytoHubba, followed by LASSO and ROC analysis to accurately identify diagnostic biomarkers. In addition, two different GEO datasets, GSE175759 and GSE47184, and an experiment cohort with 30 controls and 40 DN patients detected by IHC, were used to validate the biomarkers. Moreover, ssGSEA was performed to analyze the immune microenvironment in DN. Wilcoxon test and LASSO regression were used to determine the core immune signatures. The correlation between biomarkers and crucial immune signatures was calculated by Spearman analysis. Finally, cMap was used to explore potential drugs treating renal tubule injury in DN patients. Results: A total of 509 DEGs, including 338 upregulated and 171 downregulated genes, were screened out. "chemokine signaling pathway" and "cell adhesion molecules" were enriched in both GSEA and KEGG analysis. CCR2, CX3CR1, and SELP, especially for the combination model of the three genes, were identified as core biomarkers with high diagnostic capabilities with striking AUC, sensitivity, and specificity in both merged and validated datasets and IHC validation. Immune infiltration analysis showed a notable infiltration advantage for APC co-stimulation, CD8+ T cells, checkpoint, cytolytic activity, macrophages, MHC class I, and parainflammation in the DN group. In addition, the correlation analysis showed that CCR2, CX3CR1, and SELP were strongly and positively correlated with checkpoint, cytolytic activity, macrophages, MHC class I, and parainflammation in the DN group. Finally, dilazep was screened out as an underlying compound for DN analyzed by CMap. Conclusions: CCR2, CX3CR1, and SELP are underlying diagnostic biomarkers for DN, especially in their combination. APC co-stimulation, CD8+ T cells, checkpoint, cytolytic activity, macrophages, MHC class I, and parainflammation may participate in the occurrence and development of DN. At last, dilazep may be a promising drug for treating DN. [ABSTRACT FROM AUTHOR]

Intervertebral disc (IVD) degeneration and its inflammatory microenvironment can result in discogenic pain, which has been shown to stem from the nucleus pulposus (NP). Increasing evidence suggests that mitochondrial related genes are strictly connected to cell functionality and, importantly, it can regulate cell immune activity in response to damaged associated signals. Therefore, identification of mitochondria related genes might offer new diagnostic markers and therapeutic targets for IVD degeneration. In this study, we identified key genes involved in NP tissue immune cell infiltration during IVD degeneration by bioinformatic analysis. The key modules were screened by weighted gene co-expression network analysis (WCGNA). Characteristic genes were identified by random forest analysis. Then gene set enrichment analysis (GSEA) was used to explore the signaling pathways associated with the signature genes. Subsequently, CIBERSORT was used to classify the infiltration of immune cells. Function of the hub gene was confirmed by PCR, Western blotting and ELISA. Finally, we identified MFN2 as a crucial molecule in the process of NP cell pyroptosis and NLRP3 inflammasome activation. We speculate that the increased MFN2 expression in NP tissue along with the infiltration of CD8+ T cells, NK cell and neutrophils play important roles in the pathogenesis of IVD degeneration. [ABSTRACT FROM AUTHOR]

Diabetic peripheral neuropathy (DPN) is one of the most common chronic complications in diabetes. Previous studies have shown that chronic neuroinflammation was associated with DPN. However, further research is needed to investigate the exact immune molecular mechanism underlying the pathogenesis of DPN. Expression profiles were downloaded from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were screened by R software. After functional enrichment analysis of DEGs, a protein–protein interaction (PPI) network analysis was performed. The CIBERSORT algorithm was used to evaluate the infiltration of immune cells in DPN. Next, the least absolute shrinkage and selection operator (LASSO) logistic regression and support vector machine-recursive feature elimination (SVM-RFE) algorithms were applied to identify potential DPN diagnostic markers. Finally, the results were further validated by qRT-PCR. A total of 1308 DEGs were screened in this study. Enrichment analysis identified that DEGs were significantly enriched in immune-related biological functions and pathways. Immune cell infiltration analysis found that M1 and M2 macrophages, monocytes, resting mast cells, resting CD4 memory T cells and follicular helper T cells were involved in the development of DPN. LTBP2 and GPNMB were identified as diagnostic markers of DPN. qRT-PCR results showed that 15 mRNAs, including LTBP2 and GPNMB, were differentially expressed, consistent with the microarray results. In conclusion, LTBP2 and GPNMB can be used as novel candidate molecular diagnostic markers for DPN. Furthermore, the infiltration of immune cells plays an important role in the progression of DPN. [ABSTRACT FROM AUTHOR]

Caveolin-1 (CAV1) is one of the members of the caveolae, and the role of CAV1 in esophageal cancer (ESCA) is not completely clear. In this study, we found that expression of CAV1 was downregulated in ESCA in The Cancer Genome Atlas and the Genotype-Tissue Expression (GTEx) database and we also use immunohistochemistry of tissue microarray for verification. Then, we used bioinformatics methods to investigate the prognostic value of CAV1, influence on immune cell infiltration in tumor microenvironment (TME) and responding to immunotherapy in ESCA. Our result indicated that CAV1 designs an inflamed TME in ESCA based on the evidence that CAV1 positively correlated with immunomodulators, immune score, stomal score, cancer immunity cycles, tumor-infiltrating immune cells, T cell inflamed score, and immune checkpoints. Immunophenoscore, Tumor Immune Dysfunction and Exclusion algorithms, and the mutation analysis show that the downregulated CAV1 expression indicated higher tumor mutation burden and higher rate of response to immune checkpoint inhibitors (ICIs) in the low-expression group. In a word, our study demonstrated the impact of CAV1 to the TME in ESCA and it may be a new target for ESCA immunotherapy. In addition, the expression of CAV1 can predict the clinical response to ICIs, which may provide clinical treatment guidance. [ABSTRACT FROM AUTHOR]

BackgroundMyocardial infarction (MI) ranks among the most prevalent cardiovascular diseases. Insufficient blood flow to the coronary arteries always leads to ischemic necrosis of the cardiac muscle. However, the mechanism of myocardial injury after MI remains unclear. This article aims to explore the potential common genes between mitophagy and MI and to construct a suitable prediction model.MethodsTwo Gene Expression Omnibus (GEO) datasets (GSE62646 and GSE59867) were used to screen the differential expression genes in peripheral blood. SVM, RF, and LASSO algorithm were employed to find MI and mitophagy-related genes. Moreover, DT, KNN, RF, SVM and LR were conducted to build the binary models, and screened the best model to further external validation (GSE61144) and internal validation (10-fold cross validation and Bootstrap), respectively. The performance of various machine learning models was compared. In addition, immune cell infiltration correlation analysis was conducted with MCP-Counter and CIBERSORT.ResultsWe finally identified ATG5, TOMM20, MFN2 transcriptionally differed between MI and stable coronary artery diseases. Both internal and external validation supported that these three genes could accurately predict MI withAUC = 0.914 and 0.930 by logistic regression, respectively. Additionally, functional analysis suggested that monocytes and neutrophils might be involved in mitochondrial autophagy after myocardial infarction.ConclusionThe data showed that the transcritional levels of ATG5, TOMM20 and MFN2 in patients with MI were significantly different from the control group, which might be helpful to further accurately diagnose diseases and have potential application value in clinical practice.

Severe Fever with thrombocytopenia syndrome (SFTS) is a highly fatal viral infectious disease that poses a significant threat to public health. Currently, the phase and pathogenesis of SFTS are not well understood, and there are no specific vaccines or effective treatment available. Therefore, it is crucial to identify biomarkers for diagnosing acute SFTS, which has a high mortality rate. In this study, we conducted differentially expressed genes (DEGs) analysis and WGCNA module analysis on the GSE144358 dataset, comparing the acute phase of SFTSV-infected patients with healthy individuals. Through the LASSO–Cox and random forest algorithms, a total of 2128 genes were analyzed, leading to the identification of four genes: ADIPOR1, CENPO, E2F2, and H2AC17. The GSEA analysis of these four genes demonstrated a significant correlation with immune cell function and cell cycle, aligning with the functional enrichment findings of DEGs. Furthermore, we also utilized CIBERSORT to analyze the immune cell infiltration and its correlation with characteristic genes. The results indicate that the combination of ADIPOR1, CENPO, E2F2, and H2AC17 genes has the potential as characteristic genes for diagnosing and studying the acute phase of SFTS virus (SFTSV) infection.

Abstract Background Epidemiologic studies continue to emphasize that increasing patients with inflammatory bowel disease (IBD) develop to colorectal cancer (CRC). Although the function and mechanisms of IBD-associated genes (IBDGs) in CRC tumorigenesis have been extensively researched, the implications of IBDGs in the prognosis value and tumor immunity of CRC remain unclear. Results In this study, the expression, pathological stages and prognostic value of IBDGs in CRC were systematically analyzed, and 7 prognostic genes including CDH1, CCL11, HLA–DRA, NOS2, NAT2, TIMP1 and TP53 were screened through LASSO–Cox regression analysis. Then, a prognostic signature was established based on the 7 prognostic genes, and the model exhibited a good ability in risk stratification of CRC patients. Subsequent results showed that the genetic alterations of the 7 prognostic genes exhibited more significant and extensive influence on immune cells infiltration in colon adenocarcinoma than that in rectal adenocarcinoma. Meanwhile, immune cells infiltration also showed a significant difference between low-risk group and high-risk group. What’s more, 7 prognostic genes-based risk stratification was associated with microsatellite instability, and its prognostic characteristics were significantly negatively correlated with mismatch repair genes. Conclusions This study provided a promising insight that the 7 IBDGs could be used as valuable biomarkers for prognostic diagnosis and personalized immunotherapy of CRC patients.

Background: Early diagnosis of septic shock in children is critical for prognosis. This study committed to investigate the signature genes and their connection with immune cells in pediatric septic shock. Methods: We screened a dataset of children with septic shock from the GEO database and analyzed differentially expressed genes (DEGs). Functional enrichment analysis was performed for these DEGs. Weighted gene co-expression network analysis (WCGNA) was used to screen the key modules. Least absolute shrinkage and selection operator (LASSO) and random forest analysis were finally applied to identify the signature genes. Then gene set enrichment analysis (GSEA) was exerted to explore the signaling pathways related to the hub genes. And the immune cells infiltration was subsequently classified via using CIBERSORT. Results: A total of 534 DEGs were screened from GSE26440. The data then was clustered into 17 modules via WGCNA, which MEgrey module was significantly related to pediatric septic shock (cor=-0.62, p<0.0001). LASSO and random forest algorithms were applied to select the signature genes, containing UPP1, S100A9, KIF1B, S100A12, SLC26A8. The receiver operating characteristic curve (ROC) of these signature genes was 0.965, 0.977, 0.984, 0.991 and 0.989, respectively, which were verified in the external dataset from GSE13904. GSEA analysis showed these signature genes involve in positively correlated fructose and mannose metabolism and starch and sucrose metabolism signaling pathway. CIBERSORT suggested these signature genes may participate in immune cells infiltration. Conclusion: UPP1, S100A9, KIF1B, S100A12, SLC26A8 emerge remarkable diagnostic performance in pediatric septic shock and involved in immune cells infiltration. [ABSTRACT FROM AUTHOR]

Highlight Rcl1 mRNA expression is down-regulated within HCC tissues and associated with poor prognosis and disease progression. Anti-cancer effects of Rcl1 on HCC were confirmed in vitro. Rcl1 may be a potential tumor suppressor in HCC.

RIO kinase 2 has emerged as a critical kinase for ribosome maturation, and recently it has also been found to play a fundamental role in cancer, being involved in the occurrence and progression of glioblastoma, liver cancer, prostate cancer, non-small cell lung cancer, and acute myeloid leukemia. However, our knowledge in this regard is fragmented and limited and it is difficult to determine the exact role of RIO kinase 2 in tumors. Here, we conducted an integrated pan-cancer analysis comprising 33 cancer-types to determine the function of RIO kinase 2 in malignancies. The results show that RIO kinase 2 is highly expressed in all types of cancer and is significantly associated with tumor survival, metastasis, and immune cell infiltration. Moreover, RIO kinase 2 alteration via DNA methylation, and protein phosphorylation are involved in tumorigenesis. In summary, RIO kinase two serves as a promising target for the identification of cancer and increases our understanding of tumorigenesis and cancer progression and enhancing the ultimate goal of improved treatment for these diseases.

Neuroinflammation is an important mediator of secondary injury pathogenesis that exerts dual beneficial and detrimental effects on pathophysiology of the central nervous system (CNS) after traumatic brain injury (TBI). Fluvoxamine is a serotonin selective reuptake inhibitor (SSRI) and has been reported to have the anti-inflammatory properties. However, the mechanisms and therapeutic effects of fluvoxamine in neuroinflammation after TBI have not be defined. In this study, we showed that fluvoxamine inhibited peripheral immune cell infiltration and glia activation at 3 days in mice subjected to TBI. Fluvoxamine treatment promoted microglial/macrophage phenotypic transformation from pro-inflammatory M1-phenotype to anti-inflammatory M2-phenotype in in vivo and in vitro experiments. In addition, fluvoxamine treatment attenuated neuronal apoptosis, blood–brain barrier (BBB) disruption, cerebrovascular damage, and post-traumatic edema formation, thereby improving neurological function of mice subjected to TBI. These findings support the clinical evaluation of fluvoxamine as a neuroprotective therapy for TBI. [ABSTRACT FROM AUTHOR]

Background: Although B‐acute lymphoblastic leukemia (B‐ALL) patients' survival has been improved dramatically, some cases still relapse. This study aimed to explore the prognosis‐related novel differentially expressed genes (DEGs) for predicting the overall survival (OS) of children and young adults (CAYAs) with B‐ALL and analyze the immune‐related factors contributing to poor prognosis. Methods: GSE48558 and GSE79533 from Gene Expression Omnibus (GEO) and clinical sample information and mRNA‐seq from Therapeutically Applicable Research to Generate Effective Treatments (TARGET) database were retrieved. Prognosis‐related key genes were enrolled to build a Cox proportional model using multivariate Cox regression. Five‐year OS of patients, clinical characteristic relevance and clinical independence were assessed based on the model. The mRNA levels of prognosis‐related genes were validated in our samples and the difference of immune cells composition between high‐risk and low‐risk patients were compared. Results: One hundred and twelve DEGs between normal B cells and B‐ALL cells were identified based on GSE datasets. They were mainly participated in protein binding and HIF‐1 signaling pathway. One hundred and eighty‐nine clinical samples were enrolled in the study, both Kaplan–Meier (KM) analysis and univariate Cox regression analysis showed that CYBB, BCL2A1, IFI30, and EFNB1 were associated with prognosis, CYBB, BCL2A1, and EFNB1 were used to construct prognostic risk model. Moreover, compared to clinical indicators, the three‐gene signature was an independent prognostic factor for CAYAs with B‐ALL. Finally, the mRNA levels of CYBB, BCL2A1, and EFNB1 were significantly lower in B‐ALL group as compared to controls. The high‐risk group had a significantly higher percentage of infiltrated immune cells. Conclusion: We constructed a novel three‐gene signature with independent prognostic factor for predicting 5‐year OS of CAYAs with B‐ALL. Additionally, we discovered the difference of immune cells composition between high‐risk and low‐risk groups. This study may help to customize individual treatment and improve prognosis of CAYAs with B‐ALL. [ABSTRACT FROM AUTHOR]

Objective: Inflammation plays an important role in the transformation of pulmonary nodules (PNs) from benign to malignant. Prediction of benignancy and malignancy of PNs is still lacking efficacy methods. Although Mayo or Brock model have been widely applied in clinical practices, their application conditions are limited. This study aims to construct a diagnostic model of PNs by machine learning using inflammation-related biological markers (IRBMs)., Methods: Inflammatory related genes (IRGs) were first extracted from GSE135304 chip data. Then, differentially expressed genes (DEGs) and infiltrating immune cells were screened between malignant pulmonary nodules (MN) and benign pulmonary nodule (BN). Correlation analysis was performed on DEGs and infiltrating immune cells. Molecular modules of IRGs were identified through Consistency cluster analysis. Subsequently, IRBMs in IRGs modules were filtered through Weighted gene co-expression network analysis (WGCNA). An optimal diagnostic model was established using machine learning methods. Finally, external dataset GSE108375 was used to verify this result., Results: 4 hub IRGs and 3 immune cells showed significantly difference between MN and BN, C1 and C2 module, namely PRTN3, ELANE, NFKB1 and CTLA4, T cells CD4 naïve, NK cells activated and Monocytes. IRBMs were screened from black module and yellowgreen module through WGCNA analysis. The Support vector machines (SVM) was identified as the optimal model with the Area Under Curve (AUC) was 0.753. A nomogram was established based on 5 hub IRBMs, namely HS.137078, KLC3, C13ORF15, STOM and KCTD13. Finally, external dataset GSE108375 verified this result, with the AUC was 0.718., Conclusion: SVM model established by 5 hub IRBMs was able to effectively identify MN or BN. Accumulating inflammation and immune dysfunction were important to the transformation from BN to MN., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2024 The Authors. Published by Elsevier Ltd.)

Shengyi Long,1,&ast; Fang Peng,1,&ast; Baohui Song,1,&ast; Liang Wang,2 Jun Chen,2 Bingbing Shang1 1The Second Hospital of Dalian Medical University, Dalian, Liaoning Province, People’s Republic of China; 2Laboratory Animal Center, Dalian Medical University, Dalian, Liaoning Province, People’s Republic of China&ast;These authors contributed equally to this workCorrespondence: Jun Chen Email chenjun@dmu.edu.cnBingbing Shang Email shangbingbing@dmu.edu.cnBackground: Hepatocellular carcinoma (HCC) is one of the most serious malignancies. The main features of HCC are vascular invasion and drug resistance. Ferroptosis is a novel cell program that is involved in several diseases, such as cancer. Heat shock protein beta 1 (HSPB1) is a major component of heat shock proteins. A recent study showed that HSPB1 could be a new therapeutic target for colorectal cancer with 5-fluorouracil-acquired resistance. However, the functional role of HSPB1 in HCC remains unclear.Aim: The aim of this study is to clarify HSPB1 expression in HCC and its potential therapeutic and prognostic value.Methods: We collected data on HSPB1 expression levels in HCC and normal liver tissues from The Cancer Genome Atlas and Gene Expression Omnibus databases. We then validated it using immunohistochemistry (IHC). Receiver operating characteristic and Kaplan–Meier survival curves were used to investigate the role of HSPB1 in the prognosis analysis of HCC. Further, we used the online Search Tool for the Retrieval of Interacting Genes/Proteins website, Gene Ontology, and Kyoto Encyclopedia of Genes and Genomes to conduct enrichment analysis and identify the predictive signaling pathways. Meanwhile, we used the TIMER and GSVA package of R (v3.6.3) to analyze the association between HSPB1 and immunocyte infiltration.Results: Compared to normal tissues, there was differential expression of HSPB1 in pan-cancers. HSPB1 expression was higher in HCC tissues than in normal tissues (p< 0.05). There was an evident significant difference between HSPB1 mRNA levels and histologic grade, vascular invasion, and alpha-fetoprotein level (all p values< 0.05). Univariate analysis indicated that HCC patients with high HSPB1 levels had shorter overall survival rates than those with low HSPB1 levels (p< 0.05). MAPK14, HSPA8, MAPKAPK3, MAPKAPK5, and MAPKAPK2 are essential proteins that interact with HSPB1. There was a significant correlation between HSPB1 expression levels and immune cell infiltration, including CD4+ T cells (r=0.203, p< 0.05).Conclusion: High HSPB1 expression is closely associated with a worse prognosis in HCC patients, and HSPB1 may be a target of immunotherapy in HCC.Keywords: hepatocellular carcinoma, ferroptosis, HSPB1, immune cells infiltration, bioinformatics analysis

Pancreatic ductal adenocarcinoma (PDAC) is the most aggressive digestive malignancy due to frequent late-stage diagnosis, rapid progression and resistance to therapy. With increasing PDAC incidence worldwide, there is an urgent need for new prognostic biomarkers and therapy targets. Recently, RNA methylation has emerged as a new tumorigenic mechanism in different cancers. 5-methylcytosine (m5C) is one of the most frequent RNA modifications and occurs on a variety of RNA species including mRNA, thereby regulating gene expression. Here we investigated the prognostic role of m5C-regulator-associated transcriptional signature in PDAC. We evaluated m5C-regulator status and expression in 239 PDAC samples from publicly available datasets. We used unsupervised consensus clustering analyses to classify PDACs based on m5C-regulator expression. From the resulting signature of differentially expressed genes (DEGs), we selected prognosis-relevant DEGs to stratify patients and build a scoring signature (m5C-score) through LASSO and multivariate Cox regression analyses. The m5C-score represented a highly significant independent prognostic marker. A high m5C-score correlated with poor prognosis in different PDAC cohorts, and was associated with the squamous/basal subtype as well as activated cancer-related pathways including Ras, MAPK and PI3K pathways. Furthermore, the m5C-score correlated with sensitivity to pathway-specific inhibitors of PARP, EGFR, AKT, HER2 and mTOR. Tumors with high m5C-score were characterized by overall immune exclusion, low CD8+ T-cell infiltration, and higher PD-L1 expression. Overall, the m5C-score represented a robust predictor of prognosis and therapy response in PDAC, which was associated with unfavorable molecular subtypes and immune microenvironment.

Abstract Background Although B‐acute lymphoblastic leukemia (B‐ALL) patients' survival has been improved dramatically, some cases still relapse. This study aimed to explore the prognosis‐related novel differentially expressed genes (DEGs) for predicting the overall survival (OS) of children and young adults (CAYAs) with B‐ALL and analyze the immune‐related factors contributing to poor prognosis. Methods GSE48558 and GSE79533 from Gene Expression Omnibus (GEO) and clinical sample information and mRNA‐seq from Therapeutically Applicable Research to Generate Effective Treatments (TARGET) database were retrieved. Prognosis‐related key genes were enrolled to build a Cox proportional model using multivariate Cox regression. Five‐year OS of patients, clinical characteristic relevance and clinical independence were assessed based on the model. The mRNA levels of prognosis‐related genes were validated in our samples and the difference of immune cells composition between high‐risk and low‐risk patients were compared. Results One hundred and twelve DEGs between normal B cells and B‐ALL cells were identified based on GSE datasets. They were mainly participated in protein binding and HIF‐1 signaling pathway. One hundred and eighty‐nine clinical samples were enrolled in the study, both Kaplan–Meier (KM) analysis and univariate Cox regression analysis showed that CYBB, BCL2A1, IFI30, and EFNB1 were associated with prognosis, CYBB, BCL2A1, and EFNB1 were used to construct prognostic risk model. Moreover, compared to clinical indicators, the three‐gene signature was an independent prognostic factor for CAYAs with B‐ALL. Finally, the mRNA levels of CYBB, BCL2A1, and EFNB1 were significantly lower in B‐ALL group as compared to controls. The high‐risk group had a significantly higher percentage of infiltrated immune cells. Conclusion We constructed a novel three‐gene signature with independent prognostic factor for predicting 5‐year OS of CAYAs with B‐ALL. Additionally, we discovered the difference of immune cells composition between high‐risk and low‐risk groups. This study may help to customize individual treatment and improve prognosis of CAYAs with B‐ALL.

Background: Atopic dermatitis (AD) is a dermatological disorder characterized by symptoms such as chronically inflamed skin and frequently intolerable itching. The mechanism underlying AD development is still unclear. Our study aims to identify the diagnostic and therapeutic biomarkers for AD and provide insight into immune mechanisms at the molecular level through bioinformatics analysis.Methods: The GSE6012, GSE32924, and GSE36842 gene expression profiles were obtained for analysis from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) were segregated using the “Batch correction” and “RobustRankAggreg” methods. Weighted gene co-expression network analysis (WGCNA) was performed to screen for module genes with AD traits. Then, common DEGs (co-DEGs) were screened out via combined differential expression analysis and WGCNA. Functional enrichment analysis was performed for these co-DEGs using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG), followed by protein-protein interaction network analysis. Candidate hub genes were identified using the “cytoHubba” plugin in Cytoscape, and their value for AD diagnosis was validated using receiver operating characteristic curve analysis in the external database GSE120721. Immunohistochemical staining was performed for further validation. The CIBERSORT algorithm was used to evaluate skin samples obtained from healthy controls (HCs) and lesions of AD patients, to determine the extent of immune cell infiltration. The association between the identified hub genes and significant differential immune cells was analyzed using Pearson correlation analysis.Results: A total of 259 DEGs were acquired from the intersection of DEGs obtained by the two independent procedures, and 331 AD-trait module genes were separated out from the blue module via WGCNA analysis. Then, 169 co-DEGs arising from the intersection of the 259 DEGs and the 331 AD-trait module genes were obtained. We found that co-DEGs were significantly enhanced in the type I interferon and IL-17 signal transduction pathways. Thirteen potential hub genes were identified using Cytoscape. Five hub genes (CCR7, CXCL10, IRF7, MMP1, and RRM2) were identified after screening via external dataset validation and immunohistochemical analysis. We also identified four significant differential immune cells, i.e., activated dendritic cells, plasma cells, resting mast cells, and CD4+ naïve T cells, between AD patients and HCs. Moreover, the relationship between the identified hub genes and significant differential immune cells was analyzed. The results showed that the CCR7 expression level was positively correlated with the number of CD4+ naïve T cells (R = 0.42, p = 0.011).Conclusion: CCR7, CXCL10, IRF7, MMP1, and RRM2 could be potential diagnostic and therapeutic biomarkers for AD. CCR7 expression level was positively correlated with the number of CD4+ naïve T cells in AD. These findings need to be corroborated in future studies.

Rheumatoid arthritis(RA) is the most common inflammatory arthritis, and a significant cause of morbidity and mortality. RA patients' synovial inflammation contains a variety of genes and signalling pathways that are poorly understood. It was the goal of this research to discover the major biomarkers related to the course of RA and how they connect to immune cell infiltration. The Gene Expression Omnibus was used to download gene microarray data. Differential expression analysis, weighted gene co-expression network analysis (WGCNA), and least absolute shrinkage and selection operator (LASSO) regression were used to identify hub markers for RA. Single-sample GSEA was used to examine the infiltration levels of 28 immune cells and their connection to hub gene markers. The hub genes' expression in RA-HFLS and HFLS cells was verified by RT-PCR. The CCK-8 assay was applied to determine the roles of hub genes in RA. In this study, we identified 21 differentially expressed genes (DEGs) in RA. WGCNA yielded two co-expression modules, one of which exhibited the strongest connection with RA. Using a combination of differential genes, a total of 6 intersecting genes was discovered. Six hub genes were identified as possible biomarkers for RA after a lasso analysis was performed on the data. Three hub genes, CKS2, CSTA, and LY96, were found to have high diagnostic value using ROC curve analysis. They were shown to be closely related to the concentrations of several immune cells. RT-PCR confirmed that the expressions of CKS2, CSTA and LY96 were distinctly upregulated in RA‐HFLS cells compared with HFLS cells. More importantly, knockdown of CKS2 suppressed the proliferation of RA‐HFLS cells. Overall, to help diagnose and treat RA, it's expected that CKS2, CSTA, and LY96 will be available, and the aforementioned infiltration of immune cells may have a significant impact on the onset and progression of the disease. [ABSTRACT FROM AUTHOR]

Background: Epidemiologic studies continue to emphasize that increasing patients with inflammatory bowel disease (IBD) develop to colorectal cancer (CRC). Although the function and mechanisms of IBD-associated genes (IBDGs) in CRC tumorigenesis have been extensively researched, the implications of IBDGs in the prognosis value and tumor immunity of CRC remain unclear. Results: In this study, the expression, pathological stages and prognostic value of IBDGs in CRC were systematically analyzed, and 7 prognostic genes including CDH1, CCL11, HLA–DRA, NOS2, NAT2, TIMP1 and TP53 were screened through LASSO–Cox regression analysis. Then, a prognostic signature was established based on the 7 prognostic genes, and the model exhibited a good ability in risk stratification of CRC patients. Subsequent results showed that the genetic alterations of the 7 prognostic genes exhibited more significant and extensive influence on immune cells infiltration in colon adenocarcinoma than that in rectal adenocarcinoma. Meanwhile, immune cells infiltration also showed a significant difference between low-risk group and high-risk group. What's more, 7 prognostic genes-based risk stratification was associated with microsatellite instability, and its prognostic characteristics were significantly negatively correlated with mismatch repair genes. Conclusions: This study provided a promising insight that the 7 IBDGs could be used as valuable biomarkers for prognostic diagnosis and personalized immunotherapy of CRC patients. [ABSTRACT FROM AUTHOR]

Rheumatoid arthritis(RA) is the most common inflammatory arthritis, and a significant cause of morbidity and mortality. RA patients’ synovial inflammation contains a variety of genes and signalling pathways that are poorly understood. It was the goal of this research to discover the major biomarkers related to the course of RA and how they connect to immune cell infiltration. The Gene Expression Omnibus was used to download gene microarray data. Differential expression analysis, weighted gene co-expression network analysis (WGCNA), and least absolute shrinkage and selection operator (LASSO) regression were used to identify hub markers for RA. Single-sample GSEA was used to examine the infiltration levels of 28 immune cells and their connection to hub gene markers. The hub genes’ expression in RA-HFLS and HFLS cells was verified by RT-PCR. The CCK-8 assay was applied to determine the roles of hub genes in RA. In this study, we identified 21 differentially expressed genes (DEGs) in RA. WGCNA yielded two co-expression modules, one of which exhibited the strongest connection with RA. Using a combination of differential genes, a total of 6 intersecting genes was discovered. Six hub genes were identified as possible biomarkers for RA after a lasso analysis was performed on the data. Three hub genes, CKS2, CSTA, and LY96, were found to have high diagnostic value using ROC curve analysis. They were shown to be closely related to the concentrations of several immune cells. RT-PCR confirmed that the expressions of CKS2, CSTA and LY96 were distinctly upregulated in RA‐HFLS cells compared with HFLS cells. More importantly, knockdown of CKS2 suppressed the proliferation of RA‐HFLS cells. Overall, to help diagnose and treat RA, it’s expected that CKS2, CSTA, and LY96 will be available, and the aforementioned infiltration of immune cells may have a significant impact on the onset and progression of the disease.

Abstract Background Hepatocellular carcinoma (HCC) is one of the most prevalent and inflammation-associated cancers. The tumor microenvironment (TME) plays an essential role in HCC development and metastasis, leading to poor prognosis. The overall TME immune cells infiltration characterizations mediated by immune-related genes (IRGs) remain unclear. In this study, we aimed to investigate whether immune-related genes could be indicators for the prognosis of HCC patients and TME cell infiltration characterization as well as responses to immunotherapy. Methods We obtained differentially expressed immune-related genes (DE IRGs) between normal liver tissues and liver cancer tissues from The Cancer Genome Atlas (TCGA) database. To identify the prognostic genes and establish an immune risk signature, we performed univariable Cox regression survival analysis and the Least Absolute Shrinkage and Selector Operation (LASSO) regression based on the DE IRGs by robust rank aggregation method. Cox regression analysis was used to identify independent prognostic factors in HCC. We estimated the immune cell infiltration in TME via CIBERSORT and immunotherapy response through TIDE algorithm. Results We constructed an immune signature and validated its predictive capability. The immune signature included 7 differentially expressed IRGs: BIRC5, CACYBP, NR0B1, RAET1E, S100A8, SPINK5, and SPP1. The univariate and multivariate cox analysis showed that the 7-IRGs signature was a robust independent prognostic factor in the overall survival of HCC patients. The 7-IRG signature was associated with some clinical features, including gender, vascular invasion, histological grade, clinical stage, T stage. We also found that the 7-IRG signature could reflect the infiltration characterization of different immunocytes in the tumor microenvironment (TME) and had a good correlation with immune checkpoint molecules, revealing that the poor prognosis might be partly due to immunosuppressive TME. The Tumour Immune Dysfunction and Exclusion (TIDE) analysis data showed that the 7-IRG signature had great potential for indicating the immunotherapy response in HCC patients. The mutation analysis demonstrated a significant difference in the tumor mutation burden (TMB) between the high- and low-risk groups, partially explaining this signature's predictive value. Conclusion In a word, we constructed and validated a novel, immune-related prognostic signature for HCC patients. This signature could effectively indicate HCC patients' survival and immunotherapy response. And it might act as potential immunotherapeutic targets for HCC patients.

Background: RNA 3'-terminal phosphate cyclase-like protein (Rcl1) is involved in pre-rRNA processing, but its implication in cancers remains unclear. Methods: RCL1 expressions in 21 malignancies was examinated through GEPIA website portal. Clinical implication data related to RCL1 level in Hepatocellular Carcinoma (HCC) samples were downloaded through TCGA, ICGC, GEO databases. Survival analysis and gene function enrichment analyses were performed through R software. The correlation between RCL1 expression and tumor immune infiltration was assessed via the TIMER2.0 database. The effects of Rcl1 overexpression or knockdown on cell growth and metastasis was evaluated by CCK8, transwell, and cell cycle assays. Results: RCL1 expression is commonly down-regulated in HCC. The lower expression of RCL1 is associated with higher tumor stage, higher AFP level, vascular invasion, and poor prognosis. RCL1 expression has a significant correlation with immune cells infiltration in HCC, especially myeloid-derived suppressor cell (MDSC). Moreover, it was further identified that Rcl1 expression was reduced in HCC cell lines and negatively correlated with invasion of HCC cell lines. Immunofluorescence (IF) analysis revealed that the level of Rcl1 expression in the cytoplasm of HCC cells is significantly lower than that in the cytoplasm of L-02 cell. Moreover, both gain- and loss-of-function studies demonstrated that Rcl1 inhibited the growth and metastasis of HCC cells and regulated cell cycle progression in vitro. Conclusions: Rcl1 may serve as a novel tumor suppressor in HCC, and its biological effect needs further study. Highlight: Rcl1 mRNA expression is down-regulated within HCC tissues and associated with poor prognosis and disease progression. Anti-cancer effects of Rcl1 on HCC were confirmed in vitro. Rcl1 may be a potential tumor suppressor in HCC. [ABSTRACT FROM AUTHOR]

Leucine-rich repeat kinase2 (LRRK2) influences the host immune responses and correlates with the pathogenesis of inflammation, cancer as well as Parkinson' Disease. Herein, we explored the oncogenic role of LRRK2 at pan-cancer level and validated the analysis by single cell RNA-sequencing and in-vitro experiments. As a result, LRRK2 significantly correlated with the survival events. Specifically, LRRK2 increased the risk of Low-Grade Glioma whereas improved the survival probability of patients with Skin Cutaneous Melanoma. Gene set enrichment analysis demonstrated the involvement of LRRK2 in the host immune responses. Within the tumor microenvironment, LRRK2 was positively associated with the recruitment of macrophages. Furthermore, scRNA-seq and co-culture experiments demonstrated that LRRK2 deficiency impaired macrophage functions, and influenced the neoplastic progression in a cancer type-specific manner. Therefore, the present study provided a therapeutic strategy for LGG based on the interference with LRRK2 expression and activity to prevent macrophage recruitment and promote tumor eradication. • LRRK2 correlated with the recruitment of tumor-associated macrophages. • LRRK2 influenced the cross-talking within the tumor microenvironment in a cancer type-dependent manner. • LRRK2 deficiency facilitated M2 transition of macrophages and promoted the neoplastic virtues of cancer cells at the early stage of tumor progression. [ABSTRACT FROM AUTHOR]