Identification of immune-related key genes in the peripheral blood of ischaemic stroke patients using a weighted gene coexpression network analysis and machine learning

Discovery and validation of molecular patterns and immune characteristics in the peripheral blood of ischemic stroke patients

Background

Stroke is a disease with high morbidity, disability, and mortality. Immune factors play a crucial role in the occurrence of ischemic stroke (IS), but their exact mechanism is not clear. This study aims to identify possible immunological mechanisms by recognizing immune-related biomarkers and evaluating the infiltration pattern of immune cells.

Methods

We downloaded datasets of IS patients from GEO, applied R language to discover differentially expressed genes, and elucidated their biological functions using GO, KEGG analysis, and GSEA analysis. The hub genes were then obtained using two machine learning algorithms (least absolute shrinkage and selection operator (LASSO) and support vector machine-recursive feature elimination (SVM-RFE)) and the immune cell infiltration pattern was revealed by CIBERSORT. Gene-drug target networks and mRNA-miRNA-lncRNA regulatory networks were constructed using Cytoscape. Finally, we used RT-qPCR to validate the hub genes and applied logistic regression methods to build diagnostic models validated with ROC curves.

Results

We screened 188 differentially expressed genes whose functional analysis was enriched to multiple immune-related pathways. Six hub genes (ANTXR2, BAZ2B, C5AR1, PDK4, PPIH, and STK3) were identified using LASSO and SVM-RFE. ANTXR2, BAZ2B, C5AR1, PDK4, and STK3 were positively correlated with neutrophils and gamma delta T cells, and negatively correlated with T follicular helper cells and CD8, while PPIH showed the exact opposite trend. Immune infiltration indicated increased activity of monocytes, macrophages M0, neutrophils, and mast cells, and decreased infiltration of T follicular helper cells and CD8 in the IS group. The ceRNA network consisted of 306 miRNA-mRNA interacting pairs and 285 miRNA-lncRNA interacting pairs. RT-qPCR results indicated that the expression levels of BAZ2B, C5AR1, PDK4, and STK3 were significantly increased in patients with IS. Finally, we developed a diagnostic model based on these four genes. The AUC value of the model was verified to be 0.999 in the training set and 0.940 in the validation set.

Conclusion

Our research explored the immune-related gene expression modules and provided a specific basis for further study of immunomodulatory therapy of IS.

[High expression of UBE2S promotes progression of hepatocellular carcinoma by increasing cancer cell stemness]

OBJECTIVE

To investigate the expression of the ubiquitination enzyme UBE2S in different cell types in hepatocellular carcinoma (HCC) microenvironment and its impact on proliferation and stemness of HCC cells.

METHODS

TCGA and CPTAC database were used to analyze the transcriptional and promoter methylation levels and protein expressions of UBE2S in HCC. Specific expression patterns of UBE2S, intercellular communication and key transcription factors in different cell types were analyzed based on single-cell sequencing data from TISCH website. We further examined UBE2S expressions in clinical samples of HCC tissues, HCC cells and T cells using immunohistochemistry and immunofluorescence staining. We also tested the effects of UBE2S knockdown on stemness of HCC-LM3 and HepG2 cells using clone formation experiments and sphere formation assay.

RESULTS

Analysis based on TCGA database suggested significant overexpression of UBE2S in both paired and non-paired tumor tissues (P < 0.001), and its transcriptional level increased with tumor grades. The methylation level of UBE2S promoter was significantly decreased in HCC (P < 0.001), and its transcription level increased obviously in HCC with TP53 mutation (P < 0.001). Analysis of CPTAC database also demonstrated overexpression of UBE2S protein in HCC tissues (P < 0.001). Three prognostic models suggested that HCC patients with high UBE2S expression had poorer prognosis (P < 0.001). Single-cell sequencing data analysis revealed high expressions of UBE2S in T cells and high intensities of interaction between endothelial cells, epithelial cells and fibroblasts in HCC microenvironment. Immunohistochemistry and immunofluorescence staining demonstrated high UBE2S expressions in clinical samples of HCC tissues, HCC cells and T cells. In HCC-LM3 and HepG2 cells, UBE2S knockdown significantly inhibited cell clone formation and tumor sphere formation (P < 0.05).

CONCLUSION

UBE2S is highly expressed in T cells in HCC microenvironment in close correlation with a poor prognosis. High UBE2S expression promotes the stemness of HCC cells.

Lysosome-related biomarkers in preeclampsia and cancers: Machine learning and bioinformatics analysis

BACKGROUND

Lysosomes serve as regulatory hubs, and play a pivotal role in human diseases. However, the precise functions and mechanisms of action of lysosome-related genes remain unclear in preeclampsia and cancers. This study aimed to identify lysosome-related biomarkers in preeclampsia, and further explore the biomarkers shared between preeclampsia and cancers.

MATERIALS AND METHODS

We obtained GSE60438 and GSE75010 datasets from the Gene Expression Omnibus database, pre-procesed them and merged them into a training cohort. The limma package in R was used to identify the differentially expressed mRNAs between the preeclampsia and normal control groups. Differentially expressed lysosome-related genes were identified by intersecting the differentially expressed mRNAs and lysosome-related genes obtained from Gene Ontology and GSEA databases. Gene Ontology annotations and Kyoto Encyclopedia of Genes and Genomes enrichment analysis were performed using the DAVID database. The CIBERSORT method was used to analyze immune cell infiltration. Weighted gene co-expression analyses and three machine learning algorithm were used to identify lysosome-related diagnostic biomarkers. Lysosome-related diagnostic biomarkers were further validated in the testing cohort GSE25906. Nomogram diagnostic models for preeclampsia were constructed. In addition, pan-cancer analysis of lysosome-related diagnostic biomarkers were identified by was performed using the TIMER, Sangebox and TISIDB databases. Finally, the Drug-Gene Interaction, TheMarker and DSigDB Databases were used for drug-gene interactions analysis.

RESULTS

A total of 11 differentially expressed lysosome-related genes were identified between the preeclampsia and control groups. Three molecular clusters connected to lysosome were identified, and enrichment analysis demonstrated their strong relevance to the development and progression of preeclampsia. Immune infiltration analysis revealed significant immunity heterogeneity among different clusters. GBA, OCRL, TLR7 and HEXB were identified as lysosome-related diagnostic biomarkers with high AUC values, and validated in the testing cohort GSE25906. Nomogram, calibration curve, and decision curve analysis confirmed the accuracy of predicting the occurrence of preeclampsia based on OCRL and HEXB. Pan-cancer analysis showed that GBA, OCRL, TLR7 and HEXB were associated with the prognosis of patients with various tumors and tumor immune cell infiltration. Twelve drugs were identified as potential drugs for the treatment of preeclampsia and cancers.

CONCLUSION

This study identified GBA, OCRL, TLR7 and HEXB as potential lysosome-related diagnostic biomarkers shared between preeclampsia and cancers.

Identification and validation of platelet-related diagnostic markers and potential drug screening in ischemic stroke by integrating comprehensive bioinformatics analysis and machine learning

Background

Ischemic stroke (IS), caused by blood and oxygen deprivation due to cerebral thrombosis, has links to activated and aggregated platelets. Discovering platelet-related biomarkers, developing diagnostic models, and screening antiplatelet drugs are crucial for IS diagnosis and treatment.

Methods and results

Combining and normalizing GSE16561 and GSE22255 datasets identified 1,753 upregulated and 1,187 downregulated genes. Fifty-one genes in the platelet-related module were isolated using weighted gene co-expression network analysis (WGCNA) and other analyses, including 50 upregulated and one downregulated gene. Subsequent enrichment and network analyses resulted in 25 platelet-associated genes and six diagnostic markers for a risk assessment model. This model's area under the ROC curve outperformed single genes, and in the peripheral blood of the high-risk group, immune infiltration indicated a higher proportion of CD4, resting CD4 memory, and activated CD4 memory T cells, along with a lower proportion of CD8 T cells in comparison to the low-risk group. Utilizing the gene expression matrix and the CMap database, we identified two potential drugs for IS. Finally, a rat MACO/R model was used to validate the diagnostic markers' expression and the drugs' predicted anticoagulant effects.

Conclusion

We identified six IS platelet-related biomarkers (APP, THBS1, F13A1, SRC, PPBP, and VCL) for a robust diagnostic model. The drugs alpha-linolenic acid and ciprofibrate have potential antiplatelet effects in IS. This study advances early IS diagnosis and treatment.

Identification of ANXA3 as a biomarker associated with pyroptosis in ischemic stroke

BACKGROUND

Pyroptosis plays an important role in the pathological process of ischemic stroke (IS). However, the exact mechanism of pyroptosis remains unclear. This paper aims to reveal the key molecular markers associated with pyroptosis in IS.

METHODS

We used random forest learning, gene set variation analysis, and Pearson correlation analysis to screen for biomarkers associated with pyroptosis in IS. Middle cerebral artery occlusion/reperfusion (MCAO/R) and oxygen and glucose deprivation/reoxygenation (OGD/R) models were constructed in vitro and in vivo. Cells were transfected with an Annexin A3 silencing (si-ANXA3) plasmid to observe the effects of ANXA3 on OGD/R + lipopolysaccharides (LPS)-induced pyroptosis. qRT‒PCR and western blotting were used to detect the expression of potential biomarkers and pyroptotic pathways.

RESULTS

Samples from a total of 170 IS patients and 109 healthy individuals were obtained from 5 gene expression omnibus databases. Thirty important genes were analyzed by random forest learning from the differentially expressed genes. Then, we investigated the relationship between the above genes and the pyroptosis score, obtaining three potential biomarkers (ANXA3, ANKRD22, ADM). ANXA3 and ADM were upregulated in the MCAO/R model, and the fold difference in ANXA3 expression was greater. Pyroptosis-related factors (NLRP3, NLRC4, AIM2, GSDMD-N, caspase-8, pro-caspase-1, cleaved caspase-1, IL-1β, and IL-18) were upregulated in the MCAO/R model. Silencing ANXA3 alleviated the expression of pyroptosis-related factors (NLRC4, AIM2, GSDMD-N, caspase-8, pro-caspase-1, cleaved caspase-1, and IL-18) induced by OGD/R + LPS or MCAO/R.

CONCLUSION

This study identified ANXA3 as a possible pyroptosis-related gene marker in IS through bioinformatics and experiments. ANXA3 could inhibit pyroptosis through the NLRC4/AIM2 axis.

Investigating the ID3/SLC22A4 as immune-related signatures in ischemic stroke

BACKGROUND

Ischemic stroke (IS) is a fearful disease that can cause a variety of immune events. Nevertheless, precise immune-related mechanisms have yet to be systematically elucidated. This study aimed to identify immune-related signatures using machine learning and to validate them with animal experiments and single cell analysis.

METHODS

In this study, we screened 24 differentially expressed genes (DEGs) while identifying immune-related signatures that may play a key role in IS development through a comprehensive strategy between least absolute shrinkage and selection operation (LASSO) regression, support vector machine (SVM) and immune-related genes. In addition, we explored immune infiltration using the CIBERSORT algorithm. Finally, we performed validation in mouse brain tissue and single cell analysis.

RESULTS

We identified 24 DEGs for follow-up analysis. ID3 and SLC22A4 were finally identified as the better immune-related signatures through a comprehensive strategy among DEGs, LASSO, SVM and immune-related genes. RT-qPCR, western blot, and immunofluorescence revealed a significant decrease in ID3 and a significant increase in SLC22A4 in the middle cerebral artery occlusion group. Single cell analysis revealed that ID3 was mainly concentrated in endothelial_2 cells and SLC22A4 in astrocytes in the MCAO group. A CIBERSORT finds significantly altered levels of immune infiltration in IS patients.

CONCLUSIONS

This study focused on immune-related signatures after stroke and ID3 and SLC22A4 may be new therapeutic targets to promote functional recovery after stroke. Furthermore, the association of ID3 and SLC22A4 with immune cells may be a new direction for post-stroke immunotherapy.

Integrative Analysis of Machine Learning and Molecule Docking Simulations for Ischemic Stroke Diagnosis and Therapy

Due to the narrow therapeutic window and high mortality of ischemic stroke, it is of great significance to investigate its diagnosis and therapy. We employed weighted gene coexpression network analysis (WGCNA) to ascertain gene modules related to stroke and used the maSigPro R package to seek the time-dependent genes in the progression of stroke. Three machine learning algorithms were further employed to identify the feature genes of stroke. A nomogram model was built and applied to evaluate the stroke patients. We analyzed single-cell RNA sequencing (scRNA-seq) data to discern microglia subclusters in ischemic stroke. The RNA velocity, pseudo time, and gene set enrichment analysis (GSEA) were performed to investigate the relationship of microglia subclusters. Connectivity map (CMap) analysis and molecule docking were used to screen a therapeutic agent for stroke. A nomogram model based on the feature genes showed a clinical net benefit and enabled an accurate evaluation of stroke patients. The RNA velocity and pseudo time analysis showed that microglia subcluster 0 would develop toward subcluster 2 within 24 h from stroke onset. The GSEA showed that the function of microglia subcluster 0 was opposite to that of subcluster 2. AZ_628, which screened from CMap analysis, was found to have lower binding energy with Mmp12, Lgals3, Fam20c, Capg, Pkm2, Sdc4, and Itga5 in microglia subcluster 2 and maybe a therapeutic agent for the poor development of microglia subcluster 2 after stroke. Our study presents a nomogram model for stroke diagnosis and provides a potential molecule agent for stroke therapy.

Identification and verification of characteristic differentially expressed ferroptosis-related genes in osteosarcoma using bioinformatics analysis

BACKGROUND

This study identified and verified the characteristic differentially expressed ferroptosis-related genes (CDEFRGs) in osteosarcoma (OS).

METHODS

We extracted ferroptosis-related genes (FRGs), identified differentially expressed FRGs (DEFRGs) in OS, and conducted correlation analysis between DEFRGs. Next, we conducted GO and KEGG analyses to explore the biological functions and pathways of DEFRGs. Furthermore, we used LASSO and SVM-RFE algorithms to screen CDEFRGs, and evaluated its accuracy in diagnosing OS through ROC curves. Then, we demonstrated the molecular function and pathway enrichment of CDEFRGs through GSEA analysis. In addition, we evaluated the differences in immune cell infiltration between OS and NC groups, as well as the correlation between CDEFRGs expressions and immune cell infiltrations. Finally, the expression of CDEFRGs was verified through qRT-PCR, western blotting, and immunohistochemistry experiments.

RESULTS

We identified 51 DEFRGs and the expression relationship between them. GO and KEGG analysis revealed their key functions and important pathways. Based on four CDEFRGs (PEX3, CPEB1, NOX1, and ALOX5), we built the OS diagnostic model, and verified its accuracy. GSEA analysis further revealed the important functions and pathways of CDEFRGs. In addition, there were differences in immune cell infiltration between OS group and NC group, and CDEFRGs showed significant correlation with certain infiltrating immune cells. Finally, we validated the differential expression levels of four CDEFRGs through external experiments.

CONCLUSIONS

This study has shed light on the molecular pathological mechanism of OS and has offered novel perspectives for the early diagnosis and immune-targeted therapy of OS patients.

Single-cell RNA sequencing reveals that VIM and IFITM3 are vital targets of Dengzhan Shengmai capsule to protect against cerebral ischemic injury

ETHNOPHARMACOLOGICAL RELEVANCE

Ischemic stroke is one of the leading causes of mortality, but therapies are limited. Dengzhan Shengmai capsule (DZSM) was included by the Chinese Pharmacopoeia 2020 and has been broadly used for the treatment of ischemic stroke. However, the mechanism of DZSM against ischemic stroke is unclear.

AIM OF THE STUDY

This study used RNA sequencing (RNA-seq) and single-cell RNA sequencing (scRNA-seq) to investigate the mechanism of action of DZSM against ischemic stroke.

MATERIALS AND METHODS

The rats were randomly divided into six groups: the Sham, I/R (water), I/R + DZSM-L (0.1134g/kg), I/R + DZSM-H (0.4536g/kg), I/R + NMDP (20mg/kg), and I/R + Ginaton (20mg/kg). The rats were administrated drugs for 5 days then followed by the ischemic brain injury caused by middle cerebral artery occlusion (MCAO). The neuroprotective effect was assessed by infraction rate, neurological deficit scores, regional cerebral blood flow (rCBF), hematoxylin and eosin (H&E) staining, and Nissl staining. Based on RNA-seq and scRNA-seq, the vital biological processes and core targets of DZSM against cerebral ischemia were revealed. Enzyme-linked immunosorbent assay (ELISA) and immunofluorescence (IF) staining were used to investigate the vital biological processes and core targets of DZSM against ischemic stroke.

RESULTS

Administration of DZSM significantly reduced the infarction rate and Zea Longa score, Garcia JH score, and ameliorated the reduction in rCBF. And alleviated the neuronal damage, such as increased neuronal density level and Nissl bodies density level. RNA-seq analysis revealed that DZSM played important roles in inflammation and apoptosis. ELISA and IF straining validation confirmed that DZSM significantly decreased the expression of IL-6, IL-1β, TNF-α, ICAM-1, IBA-1, MMP9, and Cleaved caspase-3 in MCAO rats. ScRNA-seq analysis identified 8 core targets in neurons including HSPB1, SPP1, MT2A, GFAP, IFITM3, VIM, CRIP1, and GPD1, and VIM and IFITM3 was verified to be decreased by DZSM in neurons.

CONCLUSION

Our study illustrates the neuroprotective effect of DZSM against ischemia stroke, and VIM and IFITM3 were identified as vital targets in neurons of DZSM in protecting against MCAO-induced I/R injury.

m6A regulator-mediated RNA methylation modification patterns are involved in the regulation of the immune microenvironment in ischaemic cardiomyopathy

The role of RNA N6-methyladenosine (m6A) modification in the regulation of the immune microenvironment in ischaemic cardiomyopathy (ICM) remains largely unclear. This study first identified differential m6A regulators between ICM and healthy samples, and then systematically evaluated the effects of m6A modification on the characteristics of the immune microenvironment in ICM, including the infiltration of immune cells, the human leukocyte antigen (HLA) gene, and HALLMARKS pathways. A total of seven key m6A regulators, including WTAP, ZCH3H13, YTHDC1, FMR1, FTO, RBM15 and YTHDF3, were identified using a random forest classifier. A diagnostic nomogram based on these seven key m6A regulators could effectively distinguish patients with ICM from healthy subjects. We further identified two distinct m6A modification patterns (m6A cluster-A and m6A cluster-B) that are mediated by these seven regulators. Meanwhile, we also noted that one m6A regulator, WTAP, was gradually upregulated, while the others were gradually downregulated in the m6A cluster-A vs. m6A cluster-B vs. healthy subjects. In addition, we observed that the degree of infiltration of the activated dendritic cells, macrophages, natural killer (NK) T cells, and type-17 T helper (Th17) cells gradually increased in m6A cluster-A vs. m6A cluster-B vs. healthy subjects. Furthermore, m6A regulators, including FTO, YTHDC1, YTHDF3, FMR1, ZC3H13, and RBM15 were significantly negatively correlated with the above-mentioned immune cells. Additionally, several differential HLA genes and HALLMARKS signalling pathways between the m6A cluster-A and m6A cluster-B groups were also identified. These results suggest that m6A modification plays a key role in the complexity and diversity of the immune microenvironment in ICM, and seven key m6A regulators, including WTAP, ZCH3H13, YTHDC1, FMR1, FTO, RBM15, and YTHDF3, may be novel biomarkers for the accurate diagnosis of ICM. Immunotyping of patients with ICM will help to develop immunotherapy strategies with a higher level of accuracy for patients with a significant immune response.

Identification of predictors for neurological outcome after cardiac arrest in peripheral blood mononuclear cells through integrated bioinformatics analysis and machine learning

Neurological prognostication after cardiac arrest (CA) is important to avoid pursuing futile treatments for poor outcome and inappropriate withdrawal of life-sustaining treatment for good outcome. To predict neurological outcome after CA through biomarkers in peripheral blood mononuclear cells, four datasets were downloaded from the Gene Expression Omnibus database. GSE29546 and GSE74198 were used as training datasets, while GSE92696 and GSE34643 were used as verification datasets. The intersection of differentially expressed genes and hub genes from multiscale embedded gene co-expression network analysis (MEGENA) was utilized in the machine learning screening. Key genes were identified using support vector machine recursive feature elimination (SVM-RFE), least absolute shrinkage and selection operator (LASSO) logistic regression, and random forests (RF). The results were validated using receiver operating characteristic curve analysis. An mRNA-miRNA network was constructed. The distribution of immune cells was evaluated using cell-type identification by estimating relative subsets of RNA transcripts (CIBERSORT). Five biomarkers were identified as predictors for neurological outcome after CA, with an area under the curve (AUC) greater than 0.7: CASP8 and FADD-like apoptosis regulator (CFLAR), human protein kinase X (PRKX), miR-483-5p, let-7a-5p, and let-7c-5p. Interestingly, the combination of CFLAR minus PRKX showed an even higher AUC of 0.814. The mRNA-miRNA network consisted of 30 nodes and 76 edges. Statistical differences were found in immune cell distribution, including neutrophils, NK cells active, NK cells resting, T cells CD4 memory activated, T cells CD4 memory resting, T cells CD8, B cells memory, and mast cells resting between individuals with good and poor neurological outcome after CA. In conclusion, our study identified novel predictors for neurological outcome after CA. Further clinical and laboratory studies are needed to validate our findings.

Machine learning algorithms assisted identification of post-stroke depression associated biological features

Objectives

Post-stroke depression (PSD) is a common and serious psychiatric complication which hinders functional recovery and social participation of stroke patients. Stroke is characterized by dynamic changes in metabolism and hemodynamics, however, there is still a lack of metabolism-associated effective and reliable diagnostic markers and therapeutic targets for PSD. Our study was dedicated to the discovery of metabolism related diagnostic and therapeutic biomarkers for PSD.

Methods

Expression profiles of GSE140275, GSE122709, and GSE180470 were obtained from GEO database. Differentially expressed genes (DEGs) were detected in GSE140275 and GSE122709. Functional enrichment analysis was performed for DEGs in GSE140275. Weighted gene co-expression network analysis (WGCNA) was constructed in GSE122709 to identify key module genes. Moreover, correlation analysis was performed to obtain metabolism related genes. Interaction analysis of key module genes, metabolism related genes, and DEGs in GSE122709 was performed to obtain candidate hub genes. Two machine learning algorithms, least absolute shrinkage and selection operator (LASSO) and random forest, were used to identify signature genes. Expression of signature genes was validated in GSE140275, GSE122709, and GSE180470. Gene set enrichment analysis (GSEA) was applied on signature genes. Based on signature genes, a nomogram model was constructed in our PSD cohort (27 PSD patients vs. 54 controls). ROC curves were performed for the estimation of its diagnostic value. Finally, correlation analysis between expression of signature genes and several clinical traits was performed.

Results

Functional enrichment analysis indicated that DEGs in GSE140275 enriched in metabolism pathway. A total of 8,188 metabolism associated genes were identified by correlation analysis. WGCNA analysis was constructed to obtain 3,471 key module genes. A total of 557 candidate hub genes were identified by interaction analysis. Furthermore, two signature genes (SDHD and FERMT3) were selected using LASSO and random forest analysis. GSEA analysis found that two signature genes had major roles in depression. Subsequently, PSD cohort was collected for constructing a PSD diagnosis. Nomogram model showed good reliability and validity. AUC values of receiver operating characteristic (ROC) curve of SDHD and FERMT3 were 0.896 and 0.964. ROC curves showed that two signature genes played a significant role in diagnosis of PSD. Correlation analysis found that SDHD (r = 0.653, P &lt; 0.001) and FERM3 (r = 0.728, P &lt; 0.001) were positively related to the Hamilton Depression Rating Scale 17-item (HAMD) score.

Conclusion

A total of 557 metabolism associated candidate hub genes were obtained by interaction with DEGs in GSE122709, key modules genes, and metabolism related genes. Based on machine learning algorithms, two signature genes (SDHD and FERMT3) were identified, they were proved to be valuable therapeutic and diagnostic biomarkers for PSD. Early diagnosis and prevention of PSD were made possible by our findings.

Bioinformatic analysis of underlying mechanisms of Kawasaki disease via Weighted Gene Correlation Network Analysis (WGCNA) and the Least Absolute Shrinkage and Selection Operator method (LASSO) regression model

BACKGROUND

Kawasaki disease (KD) is a febrile systemic vasculitis involvingchildren younger than five years old. However, the specific biomarkers and precise mechanisms of this disease are not fully understood, which can delay the best treatment time, hence, this study aimed to detect the potential biomarkers and pathophysiological process of KD through bioinformatic analysis.

METHODS

The Gene Expression Omnibus database (GEO) was the source of the RNA sequencing data from KD patients. Differential expressed genes (DEGs) were screened between KD patients and healthy controls (HCs) with the "limma" R package. Weighted gene correlation network analysis (WGCNA) was performed to discover the most corresponding module and hub genes of KD. The node genes were obtained by the combination of the least absolute shrinkage and selection operator (LASSO) regression model with the top 5 genes from five algorithms in CytoHubba, which were further validated with the receiver operating characteristic curve (ROC curve). CIBERSORTx was employed to discover the constitution of immune cells in KDs and HCs. Functional enrichment analysis was performed to understand the biological implications of the modular genes. Finally, competing endogenous RNAs (ceRNA) networks of node genes were predicted using online databases.

RESULTS

A total of 267 DEGs were analyzed between 153 KD patients and 92 HCs in the training set, spanning two modules according to WGCNA. The turquoise module was identified as the hub module, which was mainly enriched in cell activation involved in immune response, myeloid leukocyte activation, myeloid leukocyte mediated immunity, secretion and leukocyte mediated immunity biological processes; included type II diabetes mellitus, nicotinate and nicotinamide metabolism, O-glycan biosynthesis, glycerolipid and glutathione metabolism pathways. The node genes included ADM, ALPL, HK3, MMP9 and S100A12, and there was good performance in the validation studies. Immune cell infiltration analysis revealed that gamma delta T cells, monocytes, M0 macrophage, activated dendritic cells, activated mast cells and neutrophils were elevated in KD patients. Regarding the ceRNA networks, three intact networks were constructed: NEAT1/NORAD/XIST-hsa-miR-524-5p-ADM, NEAT1/NORAD/XIST-hsa-miR-204-5p-ALPL, NEAT1/NORAD/XIST-hsa-miR-524-5p/hsa-miR-204-5p-MMP9.

CONCLUSION

To conclude, the five-gene signature and three ceRNA networks constructed in our study are of great value in the early diagnosis of KD and might help to elucidate our understanding of KD at the RNA regulatory level.

Identification of Key Biomarkers Associated with Immunogenic Cell Death and Their Regulatory Mechanisms in Severe Acute Pancreatitis Based on WGCNA and Machine Learning

Immunogenic cell death (ICD) is a form of programmed cell death with a strong sense of inflammatory detection, whose powerful situational awareness can cause the reactivation of aberrant immunity. However, the role of ICD in the pathogenesis of severe acute pancreatitis (SAP) has yet to be investigated. This study aims to explore the pivotal genes associated with ICD in SAP and how they relate to immune infiltration and short-chain fatty acids (SCFAs), in order to provide a theoretical foundation for further, in-depth mechanistic studies. We downloaded GSE194331 datasets from the Gene Expression Omnibus (GEO). The use of differentially expressed gene (DEG) analysis; weighted gene co-expression network analysis (WGCNA) and least absolute shrinkage and selection operator (LASSO) regression analysis allowed us to identify a total of three ICD-related hub genes (LY96, BCL2, IFNGR1) in SAP. Furthermore, single sample gene set enrichment analysis (ssGSEA) demonstrated that hub genes are closely associated with the infiltration of specific immune cells, the activation of immune pathways and the metabolism of SCFAs (especially butyrate). These findings were validated through the analysis of gene expression patterns in both clinical patients and rat animal models of SAP. In conclusion, the first concept of ICD in the pathogenesis of SAP was proposed in our study. This has important implications for future investigations into the pro-inflammatory immune mechanisms mediated by damage-associated molecular patterns (DAMPs) in the late stages of SAP.

Identification of anoikis-related genes classification patterns and immune infiltration characterization in ischemic stroke based on machine learning

Introduction

Ischemic stroke (IS) is a type of stroke that leads to high mortality and disability. Anoikis is a form of programmed cell death. When cells detach from the correct extracellular matrix, anoikis disrupts integrin junctions, thus preventing abnormal proliferating cells from growing or attaching to an inappropriate matrix. Although there is growing evidence that anoikis regulates the immune response, which makes a great contribution to the development of IS, the role of anoikis in the pathogenesis of IS is rarely explored.

Methods

First, we downloaded GSE58294 set and GSE16561 set from the NCBI GEO database. And 35 anoikis-related genes (ARGs) were obtained from GSEA website. The CIBERSORT algorithm was used to estimate the relative proportions of 22 infiltrating immune cell types. Next, consensus clustering method was used to classify ischemic stroke samples. In addition, we used least absolute shrinkage and selection operator (LASSO), support vector machine-recursive feature elimination (SVM-RFE) and random forest (RF) algorithms to screen the key ARGs in ischemic stroke. Next, we performed receiver operating characteristics (ROC) analysis to assess the accuracy of each diagnostic gene. At the same time, the nomogram was constructed to diagnose IS by integrating trait genes. Then, we analyzed the correlation between gene expression and immune cell infiltration of the diagnostic genes in the combined database. And gene ontology (GO) and kyoto encyclopedia of genes and genomes (KEGG) analysis were performed on these genes to explore differential signaling pathways and potential functions, as well as the construction and visualization of regulatory networks using NetworkAnalyst and Cytoscape. Finally, we investigated the expression pattern of ARGs in IS patients across age or gender.

Results

Our study comprehensively analyzed the role of ARGs in IS for the first time. We revealed the expression profile of ARGs in IS and the correlation with infiltrating immune cells. And The results of consensus clustering analysis suggested that we can classify IS patients into two clusters. The machine learning analysis screened five signature genes, including AKT1, BRMS1, PTRH2, TFDP1 and TLE1. We also constructed nomogram models based on the five risk genes and evaluated the immune infiltration correlation, gene-miRNA, gene-TF and drug-gene interaction regulatory networks of these signature genes. The expression of ARGs did not differ by sex or age.

Discussion

This study may provide a beneficial reference for further elucidating the pathogenesis of IS, and render new ideas for drug screening, individualized therapy and immunotherapy of IS.

Identification of crucial genes for predicting the risk of atherosclerosis with system lupus erythematosus based on comprehensive bioinformatics analysis and machine learning

BACKGROUND

Systemic lupus erythematosus (SLE) has become a major public health problem over the years, and atherosclerosis (AS) is one of the main complications of SLE associated with serious cardiovascular consequences in this patient population. The present study aimed to identify potential biomarkers for SLE patients with AS.

METHODS

Five microarray datasets (GSE50772, GSE81622, GSE100927, GSE28829, GSE37356) were downloaded from the NCBI Gene Expression Omnibus database. The Limma package was used to identify differentially expressed genes (DEGs) in AS. Weighted gene coexpression network analysis (WGCNA) was used to identify significant module genes associated with SLE. Functional enrichment analysis, protein-protein interaction (PPI) network construction, and machine learning algorithms (least absolute shrinkage and selection operator (Lasso, Support Vector Machine-Recursive Feature Elimination (SVM-RFE), and random forest) were applied to identify hub genes. Subsequently, we generated a nomogram and receiver operating characteristic curve (ROC) for predicting the risk of AS in SLE patients. Finally, immune cell infiltrations were analyzed, and Consensus Cluster Analysis was conducted based on Single Sample Gene Set Enrichment Analysis (ssGSEA) scores.

RESULTS

Five hub genes (SPI1, MMP9, C1QA, CX3CR1, and MNDA) were identified and used to establish a nomogram that yielded a high predictive performance (area under the curve 0.900-0.981). Dysregulated immune cell infiltrations were found in AS, with positive correlations with the five hub genes. Consensus clustering showed that the optimal number of subtypes was 3. Compared to subtypes A and B, subtype C presented higher expression of the five hub genes, immune cell infiltration levels and immune checkpoint expression.

CONCLUSION

Our study systematically identified five candidate hub genes (SPI1, MMP9, C1QA, CX3CR1, MNDA) and established a nomogram that could predict the risk of AS with SLE using various bioinformatic analyses and machine learning algorithms. Our findings provide the foothold for future studies on potential crucial genes for AS in SLE patients. Additionally, the dysregulated immune cell proportions and immune checkpoint expressions in AS with SLE were identified.

Diagnostic model constructed by nine inflammation-related genes for diagnosing ischemic stroke and reflecting the condition of immune-related cells

Background

Ischemic cerebral infarction is the most common type of stroke with high rates of mortality, disability, and recurrence. However, the known diagnostic biomarkers and therapeutic targets for ischemic stroke (IS) are limited. In the current study, we aimed to identify novel inflammation-related biomarkers for IS using machine learning analysis and to explore their relationship with the levels of immune-related cells in whole blood samples.

Methods

Gene expression profiles of healthy controls and patients with IS were download from the Gene Expression Omnibus. Analysis of differentially expressed genes (DEGs) was performed in healthy controls and patients with IS. Single-sample gene set enrichment analysis was performed to calculate inflammation scores, and weighted gene co-expression network analysis was used to analyze genes in significant modules associated with inflammation scores. Key DEGs in significant modules were then analyzed using LASSO regression analysis for constructing a diagnostic model. The effectiveness and specificity of the diagnostic model was verified in healthy controls and patients with IS and with cerebral hemorrhage (CH) using qRT-PCR. The relationship between diagnostic score and the levels of immune-related cells in whole blood were analyzed using Pearson correlations.

Results

A total of 831 DEGs were identified. Both chronic and acute inflammation scores were higher in patients with IS, while 54 DEGs were also clustered in the gene modules associated with chronic and acute inflammation scores. Among them, a total of 9 genes were selected to construct a diagnostic model. Interestingly, RT-qPCR showed that the diagnostic model had better diagnostic value for IS but not for CH. The levels of lymphocytes were lower in blood of patients with IS, while the levels of monocytes and neutrophils were increased. The diagnostic score of the model was negatively associated with the levels of lymphocytes and positively associated with levels of monocytes and neutrophils.

Conclusions

Taken together, the diagnostic model constructed using the inflammation-related genes TNFSF10, ID1, PAQR8, OSR2, PDK4, PEX11B, TNIP1, FFAR2, and JUN exhibited high and specific diagnostic value for IS and reflected the condition of lymphocytes, monocytes, and neutrophils in the blood. The diagnostic model may contribute to the diagnosis of IS.