An Immune-Related Signature Predicts Survival in Patients With Lung Adenocarcinoma

Identification of key genes in diabetic nephropathy based on lipid metabolism

Diabetic nephropathy (DN) is a common systemic microvascular complication of diabetes with a high incidence rate. Notably, the disturbance of lipid metabolism is associated with DN progression. The present study aimed to identify lipid metabolism-related hub genes associated with DN for improved diagnosis of DN. The gene expression profile data of DN and healthy samples (GSE142153) were obtained from the Gene Expression Omnibus database, and the lipid metabolism-related genes were obtained from the Molecular Signatures Database. Differentially expressed genes (DEGs) between DN and healthy samples were analyzed. The weighted gene co-expression network analysis (WGCNA) was performed to examine the relationship between genes and clinical traits to identify the key module genes associated with DN. Next, the Venn Diagram R package was used to identify the lipid metabolism-related genes associated with DN and their protein-protein interaction (PPI) network was constructed. Subsequently, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed. The hub genes were identified using machine-learning algorithms. The Gene Set Enrichment Analysis (GSEA) was used to analyze the functions of the hub genes. The present study also investigated the immune infiltration discrepancies between DN and healthy samples, and assessed the correlation between the immune cells and hub genes. Finally, the expression levels of key genes were verified by reverse transcription-quantitative (RT-q)PCR. The present study determined 1,445 DEGs in DN samples. In addition, 694 DN-related genes in MEyellow and MEturquoise modules were identified by WGCNA. Next, the Venn Diagram R package was used to identify 17 lipid metabolism-related genes and to construct a PPI network. GO analysis revealed that these 17 genes were markedly associated with 'phospholipid biosynthetic process' and 'cholesterol biosynthetic process', while the KEGG analysis showed that they were enriched in 'glycerophospholipid metabolism' and 'fatty acid degradation'. In addition, SAMD8 and CYP51A1 were identified through the intersections of two machine-learning algorithms. The results of GSEA revealed that the 'mitochondrial matrix' and 'GTPase activity' were the markedly enriched GO terms in both SAMD8 and CYP51A1. Their KEGG pathways were mainly concentrated in the 'pathways of neurodegeneration-multiple diseases'. Immune infiltration analysis showed that nine types of immune cells had different expression levels in DN (diseased) and healthy samples. Notably, SAMD8 and CYP51A1 were both markedly associated with activated B cells and effector memory CD8 T cells. Finally, RT-qPCR confirmed the high expression of SAMD8 and CYP51A1 in DN. In conclusion, lipid metabolism-related genes SAMD8 and CYP51A1 may play key roles in DN. The present study provides fundamental information on lipid metabolism that may aid the diagnosis and treatment of DN.

Identification of costimulatory molecule signatures for evaluating prognostic risk in non-small cell lung cancer

Background

Non-small cell lung cancer (NSCLC) is a leading cause of cancer-related mortality worldwide. Despite advances in treatment, prognosis remains poor, necessitating the identification of reliable prognostic biomarkers. Costimulatory molecules (CMs) have shown to enhance antitumor immune responses. We aimed to explore their prognostic signals in NSCLC.

Methods

This study is a combination of bioinformatics analysis and laboratory validation. Gene expression profiles from The Cancer Genome Atlas (TCGA), GSE120622, and GSE131907 datasets were collected. NSCLC samples in TCGA were clustered based on CMs using consensus clustering. We used LASSO regression to identify CMs-related signatures and constructed nomogram and risk models. Differences in immune cells and checkpoint expressions between risk models were evaluated. Enrichment analysis was performed for differentially expressed CMs between NSCLC and controls. Key results were validated using qRT-PCR and flow cytometry.

Results

NSCLC samples in TCGA were divided into two clusters based on CMs, with cluster 1 showing poor overall survival. Ten CMs-related signatures were identified using LASSO regression. NSCLC samples in TCGA were stratified into high- and low-risk groups based on the median risk score of these signatures, revealing differences in survival probability, drug sensitivity, immune cell infiltration and checkpoints expression. The area under the ROC curve values (AUC) for EDA, ICOS, PDCD1LG2, and VTCN1 exceeded 0.7 in both datasets and considered as hub genes. Expression of these hub genes was significance in GSE131907 and validated by qRT-PCR. Macrophage M1 and T cell follicular helper showed high correlation with hub genes and were lower in NSCLC than controls detected by flow cytometry.

Conclusion

The identified hub genes can serve as prognostic biomarkers for NSCLC, aiding in treatment decisions and highlighting potential targets for immunotherapy. This study provides new insights into the role of CMs in NSCLC prognosis and suggests future directions for clinical research and therapeutic strategies.

Identification of neutrophil extracellular trap-related biomarkers in non-alcoholic fatty liver disease through machine learning and single-cell analysis

Non-alcoholic Fatty Liver Disease (NAFLD), noted for its widespread prevalence among adults, has become the leading chronic liver condition globally. Simultaneously, the annual disease burden, particularly liver cirrhosis caused by NAFLD, has increased significantly. Neutrophil Extracellular Traps (NETs) play a crucial role in the progression of this disease and are key to the pathogenesis of NAFLD. However, research into the specific roles of NETs-related genes in NAFLD is still a field requiring thorough investigation. Utilizing techniques like AddModuleScore, ssGSEA, and WGCNA, our team conducted gene screening to identify the genes linked to NETs in both single-cell and bulk transcriptomics. Using algorithms including Random Forest, Support Vector Machine, Least Absolute Shrinkage, and Selection Operator, we identified ZFP36L2 and PHLDA1 as key hub genes. The pivotal role of these genes in NAFLD diagnosis was confirmed using the training dataset GSE164760. This study identified 116 genes linked to NETs across single-cell and bulk transcriptomic analyses. These genes demonstrated enrichment in immune and metabolic pathways. Additionally, two NETs-related hub genes, PHLDA1 and ZFP36L2, were selected through machine learning for integration into a prognostic model. These hub genes play roles in inflammatory and metabolic processes. scRNA-seq results showed variations in cellular communication among cells with different expression patterns of these key genes. In conclusion, this study explored the molecular characteristics of NETs-associated genes in NAFLD. It identified two potential biomarkers and analyzed their roles in the hepatic microenvironment. These discoveries could aid in NAFLD diagnosis and management, with the ultimate goal of enhancing patient outcomes.

Identification and Validation of Cytotoxicity-Related Features to Predict Prognostic and Immunotherapy Response in Patients with Clear Cell Renal Cell Carcinoma

Background

Clear cell renal cell carcinoma (ccRCC) is a renal cortical malignancy with a complex pathogenesis. Identifying ideal biomarkers to establish more accurate promising prognostic models is crucial for the survival of kidney cancer patients.

Methods

Seurat R package was used for single-cell RNA-sequencing (scRNA-seq) data filtering, dimensionality reduction, clustering, and differentially expressed genes analysis. Gene coexpression network analysis (WGCNA) was performed to identify the cytotoxicity-related module. The independent cytotoxicity-related risk model was established by the survival R package, and Kaplan-Meier (KM) survival analysis and timeROC with area under the curve (AUC) were employed to confirm the prognosis and effectiveness of the risk model. The risk and prognosis in patients suffering from ccRCC were predicted by establishing a nomogram. A comparison of the level of immune infiltration in different risk groups and subtypes using the CIBERSORT, MCP-counter, and TIMER methods, as well as assessment of drug sensitivity to conventional chemotherapeutic agents in risk groups using the pRRophetic package, was made.

Results

Eleven ccRCC subpopulations were identified by single-cell sequencing data from the GSE224630 dataset. The identified cytotoxicity-related T-cell cluster and module genes defined three cytotoxicity-related molecular subtypes. Six key genes (SOWAHB, SLC16A12, IL20RB, SLC12A8, PLG, and HHLA2) affecting prognosis risk genes were selected for developing a risk model. A nomogram containing the RiskScore and stage revealed that the RiskScore contributed the most and exhibited excellent predicted performance for prognosis in the calibration plots and decision curve analysis (DCA). Notably, high-risk patients with ccRCC demonstrate a poorer prognosis with higher immune infiltration characteristics and TIDE scores, whereas low-risk patients are more likely to benefit from immunotherapy.

Conclusions

A ccRCC survival prognostic model was produced based on the cytotoxicity-related signature, which had important clinical significance and may provide guidance for ccRCC treatment.

Identification of key differentially expressed immune related genes in patients with persistent atrial fibrillation: an integrated bioinformation analysis

OBJECTIVE

We aimed to investigate key differentially expressed immune related genes in persistent atrial fibrillation.

METHODS

Gene expression profiles were downloaded from Gene Expression Omnibus (GEO) using "GEO query" package. "limma" package and "sva" package were used to conduct normalization and eliminate batch effects, respectively. We screened out differentially expressed genes (DEGs) based on "limma" package with the standard of |log fold change (FC)| ≥ 1.5 and false discovery rate (FDR) < 0.05. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of DEGs were performed by "clusterProfler" package. We further applied LASSO to select key DEGs, and intersected key DEGs with immune related genes from ImmPort database. The ROC curve of each DEIRG was constructed to evaluate its diagnostic efficiency for AF.

RESULTS

A total of 103 DEGs we were screened out, of them, 48 genes were down-regulated and 55 genes were up-regulated. Result of functional enrichment analysis show that, most of DEGs were related to immune response, inflammation, and oxidative stress. Ultimately, CYBB, RORB, S100A12, and CHGB were determined as key DEIRGs, each of which displayed a favor efficiency for diagnosing persistent AF.

CONCLUSION

CYBB, RORB, S100A12, and CHGB were identified as key DEIRGs in persistent AF, and future studies are needed to further explore the underlying roles of CYBB, RORB, S100A12, and CHGB in persistent AF.

High-throughput transcriptome profiling indicates ribosomal RNAs to be associated with resistance to immunotherapy in non-small cell lung cancer (NSCLC)

BACKGROUND

Despite the impressive outcomes with immune checkpoint inhibitor (ICI) in non-small cell lung cancer (NSCLC), only a minority of the patients show long-term benefits from ICI. In this study, we used retrospective cohorts of ICI treated patients with NSCLC to discover and validate spatially resolved protein markers associated with resistance to programmed cell death protein-1 (PD-1) axis inhibition.

METHODS

Pretreatment samples from 56 patients with NSCLC treated with ICI were collected and analyzed in a tissue microarray (TMA) format in including four different tumor regions per patient using the GeoMx platform for spatially informed transcriptomics. 34 patients had assessable tissue with tumor compartment in all 4 TMA spots, 22 with leukocyte compartment and 12 with CD68 compartment. The patients' tissue that was not assessable in fourfold redundancy in each compartment was designated as the validation cohort; cytokeratin (CK) (N=22), leukocytes CD45 (N=31), macrophages, CD68 (N=43). The human whole transcriptome, represented by~18,000 individual genes assessed by oligonucleotide-tagged in situ hybridization, was sequenced on the NovaSeq platform to quantify the RNAs present in each region of interest.

RESULTS

54,000 gene variables were generated per case, from them 25,740 were analyzed after removing targets with expression lower than a prespecified frequency. Cox proportional-hazards model analysis was performed for overall and progression-free survival (OS, PFS, respectively). After identifying genes significantly associated with limited survival benefit (HR>1)/progression per spot per patient, we used the intersection of them across the four TMA spots per patient. This resulted in a list of 12 genes in the tumor-cell compartment (RPL13A, GNL3, FAM83A, CYBA, ACSL4, SLC25A6, EPAS1, RPL5, APOL1, HSPD1, RPS4Y1, ADI1). RPL13A, GNL3 in tumor-cell compartment were also significantly associated with OS and PFS, respectively, in the validation cohort (CK: HR, 2.48; p=0.02 and HR, 5.33; p=0.04). In CD45 compartment, secreted frizzled-related protein 2, was associated with OS in the discovery cohort but not in the validation cohort. Similarly, in the CD68 compartment ARHGAP and PNN interacting serine and arginine rich protein were significantly associated with PFS and OS, respectively, in the majority but not all four spots per patient.

CONCLUSION

This work highlights RPL13A and GNL3 as potential indicative biomarkers of resistance to PD-1 axis blockade that might help to improve precision immunotherapy strategies for lung cancer.

T cell-mediated tumor killing based signature to predict the prognosis and immunotherapy for glioblastoma

Despite the significant advancements in cancer treatment brought by immune checkpoint inhibitors (ICIs), their effectiveness in treating glioblastoma (GBM) remains highly dissatisfactory. Immunotherapy relies on the fundamental concept of T cell-mediated tumor killing (TTK). Nevertheless, additional investigation is required to explore its potential in prognostic prediction and regulation of tumor microenvironment (TME) in GBM. TTK sensitivity related genes (referred to as GSTTKs) were obtained from the TISIDB. The training cohort was available from the TCGA-GBM, while the independent validation group was gathered from GEO database. Firstly, we examined differentially expressed GSTTKs (DEGs) with limma package. Afterwards, the prognostic DEGs were identified and the TTK signature was established with univariate and LASSO Cox analyses. Next, we examined the correlation between the TTK signature and outcome of GBM as well as immune phenotypes of TME. Furthermore, the evaluation of TTK signature in predicting the effectiveness of immunotherapy has also been conducted. We successfully developed a TTK signature with an independent predictive value. Patients who had a high score experienced a worse prognosis compared to patients with low scores. The TTK signature showed a strong positive association with the infiltration degree of immunocyte and the presence of various immune checkpoints. Moreover, individuals with a lower score exhibited increased responsiveness to ICIs and experienced improved prognosis. In conclusions, we successfully developed and verified a TTK signature that has the ability to predict the outcome and immune characteristics of GBM. Furthermore, the TTK signature has the potential to direct the personalized immunotherapy for GBM.

Comprehensive analysis reveals that LTBR is a immune-related biomarker for glioma

Glioma is a common malignant brain tumor with great heterogeneity and huge difference in clinical outcomes. Although lymphotoxin (LT) beta receptor (LTBR) has been linked to immune system and response development for decades, the expression and function in glioma have not been investigated. To confirm the expression profile of LTBR, integrated RNA-seq data from glioma and normal brain tissues were analyzed. Functional enrichment analysis, TMEscore analysis, immune infiltration, the correlation of LTBR with immune checkpoints and ferroptosis, and scRNAseq data analysis in gliomas were in turn performed, which pointed out that LTBR was pertinent to immune functions of macrophages in gliomas. In addition, after being trained and validated in the tissue samples of the integrated dataset, an LTBR DNA methylation-based prediction model succeeded to distinguish gliomas from non-gliomas, as well as the grades of glioma. Moreover, by virtue of the candidate LTBR CpG sites, a prognostic risk-score model was finally constructed to guide the chemotherapy, radiotherapy, and immunotherapy for glioma patients. Taken together, LTBR is closely correlated with immune functions in gliomas, and LTBR DNA methylation could serve as a biomarker for diagnosis and prognosis of gliomas.

25-hydroxycholesterol promotes proliferation and metastasis of lung adenocarcinoma cells by regulating ERβ/TNFRSF17 axis

Lung adenocarcinoma is the main type of lung cancer in women. Our previous findings have evidenced that 25-hydroxycholesterol (25-HC) promotes migration and invasion of lung adenocarcinoma cells (LAC), during which LXR as a 25-HC receptor plays an important role. Estrogen receptor beta (ERβ) is a receptor of 27-hydroxycholesterol that is structurally analogous to 25-HC, but its role in the functional actions of 25-HC remained largely unknown. In this study, we demonstrated that 25-HC treatment triggered ERβ expression in LAC. Knockdown of ERβ inhibited 25-HC-mediated proliferation, migration and invasion, and reduced 25-HC-induced LAC metastasis in vivo. Further investigation revealed that ERβ knockdown restrained the expression of TNFRSF17 (BCMA). In vivo experiments also confirmed that ERβ knockdown blocked 25-HC-induced TNFRSF17 expression. TNFRSF17 knockdown also restrained 25-HC-induced proliferation, migration and invasion. Bioinformatic analysis showed that the levels of ERβ and TNFRSF17 were elevated in lung adenocarcinoma, and were closely related to tumor stages and nodal metastasis status. These results suggested that 25-HC promoted the proliferation and metastasis of LAC by regulating ERβ/TNFRSF17 axis.

Development of machine learning-based malignant pericardial effusion-related model in breast cancer: Implications for clinical significance, tumor immune and drug-therapy

Background

Malignant pericardial effusion (MPE) is a common complication of advanced breast cancer (BRCA) and plays an important role in BRCA. This study is aims to construct a prognostic model based on MPE-related genes for predicting the prognosis of breast cancer.

Methods

The BRCA samples are analyzed based on the expression of MPE-related genes by using an unsupervised cluster analysis method. This study processes the data by least absolute shrinkage and selection operator and multivariate Cox analysis, and uses machine learning algorithms to construct BRCA prognostic model and develop web tool.

Results

BRCA patients are classified into three clusters and a BRCA prognostic model is constructed containing 9 MPE-related genes. There are significant differences in signature pathways, immune infiltration, immunotherapy response and drug sensitivity testing between the high and low-risk groups. Of note, a web-based tool (http://wys.helyly.top/cox.html) is developed to predict overall survival as well as drug-therapy response of BRCA patients quickly and conveniently, which can provide a basis for clinicians to formulate individualized treatment plans.

Conclusion

The MPE-related prognostic model developed in this study can be used as an effective tool for predicting the prognosis of BRCA and provides new insights for the diagnosis and treatment of BRCA patients.

Identifying of immune-associated genes for assessing the obesity-associated risk to the offspring in maternal obesity: A bioinformatics and machine learning

BACKGROUND

Perinatal exposure to maternal obesity predisposes offspring to develop obesity later in life. Immune dysregulation in the hypothalamus, the brain center governing energy homeostasis, is pivotal in obesity development. This study aimed to identify key candidate genes associated with the risk of offspring obesity in maternal obesity.

METHODS

We obtained obesity-related datasets from the Gene Expression Omnibus (GEO) database. GSE135830 comprises gene expression data from the hypothalamus of mouse offspring in a maternal obesity model induced by a high-fat diet model (maternal high-fat diet (mHFD) group and maternal chow (mChow) group), while GSE127056 consists of hypothalamus microarray data from young adult mice with obesity (high-fat diet (HFD) and Chow groups). We identified differentially expressed genes (DEGs) and module genes using Limma and weighted gene co-expression network analysis (WGCNA), conducted functional enrichment analysis, and employed a machine learning algorithm (least absolute shrinkage and selection operator (LASSO) regression) to pinpoint candidate hub genes for diagnosing obesity-associated risk in offspring of maternal obesity. We constructed a nomogram receiver operating characteristic (ROC) curve to evaluate the diagnostic value. Additionally, we analyzed immune cell infiltration to investigate immune cell dysregulation in maternal obesity. Furthermore, we verified the expression of the candidate hub genes both in vivo and in vitro.

RESULTS

The GSE135830 dataset revealed 2868 DEGs between the mHFD offspring and the mChow group and 2627 WGCNA module genes related to maternal obesity. The overlap of DEGs and module genes in the offspring with maternal obesity in GSE135830 primarily enriched in neurodevelopment and immune regulation. In the GSE127056 dataset, 133 DEGs were identified in the hypothalamus of HFD-induced adult obese individuals. A total of 13 genes intersected between the GSE127056 adult obesity DEGs and the GSE135830 maternal obesity module genes that were primarily enriched in neurodevelopment and the immune response. Following machine learning, two candidate hub genes were chosen for nomogram construction. Diagnostic value evaluation by ROC analysis determined Sytl4 and Kncn2 as hub genes for maternal obesity in the offspring. A gene regulatory network with transcription factor-miRNA interactions was established. Dysregulated immune cells were observed in the hypothalamus of offspring with maternal obesity. Expression of Sytl4 and Kncn2 was validated in a mouse model of hypothalamic inflammation and a palmitic acid-stimulated microglial inflammation model.

CONCLUSION

Two candidate hub genes (Sytl4 and Kcnc2) were identified and a nomogram was developed to predict obesity risk in offspring with maternal obesity. These findings offer potential diagnostic candidate genes for identifying obesity-associated risks in the offspring of obese mothers.

A diagnostic model based on bioinformatics and machine learning to differentiate bipolar disorder from schizophrenia and major depressive disorder

Bipolar disorder (BD) showed the highest suicide rate of all psychiatric disorders, and its underlying causative genes and effective treatments remain unclear. During diagnosis, BD is often confused with schizophrenia (SC) and major depressive disorder (MDD), due to which patients may receive inadequate or inappropriate treatment, which is detrimental to their prognosis. This study aims to establish a diagnostic model to distinguish BD from SC and MDD in multiple public datasets through bioinformatics and machine learning and to provide new ideas for diagnosing BD in the future. Three brain tissue datasets containing BD, SC, and MDD were chosen from the Gene Expression Omnibus database (GEO), and two peripheral blood datasets were selected for validation. Linear Models for Microarray Data (Limma) analysis was carried out to identify differentially expressed genes (DEGs). Functional enrichment analysis and machine learning were utilized to identify. Least absolute shrinkage and selection operator (LASSO) regression was employed for identifying candidate immune-associated central genes, constructing protein-protein interaction networks (PPI), building artificial neural networks (ANN) for validation, and plotting receiver operating characteristic curve (ROC curve) for differentiating BD from SC and MDD and creating immune cell infiltration to study immune cell dysregulation in the three diseases. RBM10 was obtained as a candidate gene to distinguish BD from SC. Five candidate genes (LYPD1, HMBS, HEBP2, SETD3, and ECM2) were obtained to distinguish BD from MDD. The validation was performed by ANN, and ROC curves were plotted for diagnostic value assessment. The outcomes exhibited the prediction model to have a promising diagnostic value. In the immune infiltration analysis, Naive B, Resting NK, and Activated Mast Cells were found to be substantially different between BD and SC. Naive B and Memory B cells were prominently variant between BD and MDD. In this study, RBM10 was found as a candidate gene to distinguish BD from SC; LYPD1, HMBS, HEBP2, SETD3, and ECM2 serve as five candidate genes to distinguish BD from MDD. The results obtained from the ANN network showed that these candidate genes could perfectly distinguish BD from SC and MDD (76.923% and 81.538%, respectively).

Alveolar Organoids in Lung Disease Modeling

Lung organoids display a tissue-specific functional phenomenon and mimic the features of the original organ. They can reflect the properties of the cells, such as morphology, polarity, proliferation rate, gene expression, and genomic profile. Alveolar type 2 (AT2) cells have a stem cell potential in the adult lung. They produce and secrete pulmonary surfactant and proliferate to restore the epithelium after damage. Therefore, AT2 cells are used to generate alveolar organoids and can recapitulate distal lung structures. Also, AT2 cells in human-induced pluripotent stem cell (iPSC)-derived alveolospheres express surfactant proteins and other factors, indicating their application as suitable models for studying cell-cell interactions. Recently, they have been utilized to define mechanisms of disease development, such as COVID-19, lung cancer, idiopathic pulmonary fibrosis, and chronic obstructive pulmonary disease. In this review, we show lung organoid applications in various pulmonary diseases, drug screening, and personalized medicine. In addition, stem cell-based therapeutics and approaches relevant to lung repair were highlighted. We also described the signaling pathways and epigenetic regulation of lung regeneration. It is critical to identify novel regulators of alveolar organoid generations to promote lung repair in pulmonary diseases.

Novel diagnostic biomarkers of oxidative stress, ferroptosis, immune infiltration characteristics and experimental validation in ischemic stroke

Ischemic stroke (IS) is a prominent type of cerebrovascular disease leading to death and disability in an aging society and is closely related to oxidative stress. Gene expression profiling (GSE222551) was derived from Gene Expression Omnibus (GEO), and 1934 oxidative stress (OS) genes were obtained from the GeneCards database. Subsequently, we identified 149 differentially expressed genes related to OS (DEOSGs). Finally, PTGS2, FOS, and RYR1 were identified as diagnostic markers of IS. Moreover, GSE16561 was used to validate the DEOSGs. Two diagnostic genes (PTGS2 and FOS) were significantly highly expressed, while RYR1 was significantly lowly expressed in the IS group. Remarkably, immune infiltration characteristics of these three genes were analyzed, and we found that PTGS2, FOS, and RYR1 were mainly correlated with Mast cells activated, Neutrophils, and Plasma cells, respectively. Next, we intersected three DEOSGs with the ferroptosis gene set, the findings revealed that only PTGS2 was a differentially expressed gene of ferroptosis. High PTGS2 expression levels in the infarcted cortex of middle cerebral artery occlusion (MCAO) rats were confirmed by immunofluorescence (IF), western blotting (WB), and Immunohistochemistry (IHC). Inhibition of PTGS2 clearly improved the neurological outcome of rats by decreasing infarct volume, neurological problems, and modified neurological severity scores following IS compared with the controls. The protective effect of silencing PTGS2 may be related to anti-oxidative stress and ferroptosis. In conclusion, this work may provide a new perspective for the research of IS, and further research based on PTGS2 may be a breakthrough.

Construction of a predictive model for blood transfusion in patients undergoing total hip arthroplasty and identification of clinical heterogeneity

A precise forecast of the need for blood transfusions (BT) in patients undergoing total hip arthroplasty (THA) is a crucial step toward the implementation of precision medicine. To achieve this goal, we utilized supervised machine learning (SML) techniques to establish a predictive model for BT requirements in THA patients. Additionally, we employed unsupervised machine learning (UML) approaches to identify clinical heterogeneity among these patients. In this study, we recruited 224 patients undergoing THA. To identify factors predictive of BT during the perioperative period of THA, we employed LASSO regression and the random forest (RF) algorithm as part of supervised machine learning (SML). Using logistic regression, we developed a predictive model for BT in THA patients. Furthermore, we utilized unsupervised machine learning (UML) techniques to cluster THA patients who required BT based on similar clinical features. The resulting clusters were subsequently visualized and validated. We constructed a predictive model for THA patients who required BT based on six predictive factors: Age, Body Mass Index (BMI), Hemoglobin (HGB), Platelet (PLT), Bleeding Volume, and Urine Volume. Before surgery, 1 h after surgery, 1 day after surgery, and 1 week after surgery, significant differences were observed in HGB and PLT levels between patients who received BT and those who did not. The predictive model achieved an AUC of 0.899. Employing UML, we identified two distinct clusters with significantly heterogeneous clinical characteristics. Age, BMI, PLT, HGB, bleeding volume, and urine volume were found to be independent predictors of BT requirement in THA patients. The predictive model incorporating these six predictors demonstrated excellent predictive performance. Furthermore, employing UML enabled us to classify a heterogeneous cohort of THA patients who received BT in a meaningful and interpretable manner.

Construction of RNA Methylation Modification-immune-related lncRNA Molecular Subtypes and Prognostic Scoring System in Lung Adenocarcinoma

BACKGROUND

RNA methylation modification is not only intimately interrelated with cancer development and progression but also actively influences immune cell infiltration in the tumor microenvironment (TME). RNA methylation modification genes influence the therapeutic progression of lung adenocarcinoma (LUAD), and mining RNA methylation modification prognosis-related markers in LUAD is crucial for its precise prognosis.

METHODS

RNA-Seq data and Gene sets were collected from online databases or published literature. Genomic variation analysis was conducted by the Maftools package. RNA methylation-immune-related lncRNAs were obtained by Pearson correlation analysis. Then, Consistent clustering analysis was performed to obtain RNA methylation modification- immune molecular subtypes (RMM-I Molecular subtypes) in LUAD based on selected lncRNAs. COX and random survival forest analysis were carried out to construct the RMM-I Score. The receiver operating characteristic (ROC) curve and Kaplan Meier survival analysis were used to assess survival differences. Tumor immune microenvironment was assessed through related gene signatures and CIBERSORT algorithm. In addition, drug sensitivity analysis was executed by the pRRophetic package.

RESULTS

Four RNA methylation modified-immune molecular subtypes (RMM-I1, RMM- I2, RMM-I3, RMM-I4) were presented in LUAD. Patients in RMM-I4 exhibited excellent survival advantages and immune activity. HAVCR2, CD274, and CTLA-4 expression were activated in RMM-I4, which might be heat tumors and a potential beneficial group for immunotherapy. OGFRP1, LINC01116, DLGAP1-AS2, CRNDE, LINC01137, MIR210HG, and CYP1B1-AS1 comprised the RMM-I Score. The RMM-I Score exhibited excellent accuracy in the prognostic assessment of LUAD, as patients with a low RMM- I Score exhibited remarkable survival advantage. Patients with a low RMM-I score might be more sensitive to treatment with Docetaxel, Vinorelbine, Paclitaxel, Cisplatin, and immunotherapy.

CONCLUSION

The RMM-I molecular subtype constituted the novel molecular characteristic subtype of LUAD, which complemented the existing pathological typing. More refined and accurate molecular subtypes provide help to reveal the mechanism of LUAD development. In addition, the RMM-I score offers a reliable tool for accurate prognosis of LUAD.

Identification of neutrophil extracellular traps and crosstalk genes linking inflammatory bowel disease and osteoporosis by integrated bioinformatics analysis and machine learning

Musculoskeletal deficits are among the most common extra-intestinal manifestations and complications of inflammatory bowel disease (IBD). This study aimed to identify crosstalk genes between IBD and osteoporosis (OP) and potential relationships between crosstalk and neutrophil extracellular traps (NETs)-related genes. Three common hub genes from different compared groups are actually the same, namely HDAC6, IL-8, and PPIF. ROC showed that the combined diagnostic value of HDAC6, IL-8, and PPIF was higher than each of the three key hub genes. Immune infiltration results showed that HDAC6 and IL-8 key genes negatively correlated with CD65 bright natural killer cells. USF1 was the common upstream TFs between HDAC6 and PPIF, and MYC was the common upstream TFs between IL-8 and PPIF in RegNetwork. Taken together, this study shows a linked mechanism between IBD and OP via NETs and crosstalk genes. These findings may show light on better diagnosis and treatment of IBD complicated with OP.

Comprehensive analysis reveals key genes and environmental toxin exposures underlying treatment response in ulcerative colitis based on in-silico analysis and Mendelian randomization

BACKGROUND

UC is increasingly prevalent worldwide and represents a significant global disease burden. Although medical therapeutics are employed, they often fall short of being optimal, leaving patients struggling with treatment non-responsiveness and many related complications.

MATERIALS AND METHODS

The study utilized gene microarray data and clinical information from GEO. Gene enrichment and differential expression analyses were conducted using Metascape and Limma, respectively. Lasso Regression Algorithm was constructed using glmnet and heat maps were generated using pheatmap. ROC curves were used to assess diagnostic parameter capability, while XSum was employed to screen for small-molecule drugs exacerbating UC. Molecular docking was carried out using Autodock Vina. The study also performed Mendelian randomization analysis based on TwoSampleMR and used CTD to investigate the relationship between exposure to environmental chemical toxicants and UC therapy responsiveness.

RESULTS

Six genes (ELL2, DAPP1, SAMD9L, CD38, IGSF6, and LYN) were found to be significantly overexpressed in UC patient samples that did not respond to multiple therapies. Lasso analysis identified ELL2 and DAPP1 as key genes influencing UC treatment response. Both genes accurately predicted intestinal inflammation in UC and impacted the immunological infiltration status. Clofibrate showed therapeutic potential for UC by binding to ELL2 and DAPP1 proteins. The study also reviews environmental toxins and drug exposures that could impact UC progression.

CONCLUSIONS

We used microarray technology to identify DAPP1 and ELL2 as key genes that impact UC treatment response and inflammatory progression. Clofibrate was identified as a promising UC treatment. Our review also highlights the impact of environmental toxins on UC treatment response, providing valuable insights for personalized clinical management.

Identification and study of cuproptosis-related genes in prognostic model of multiple myeloma

BACKGROUND

Multiple myeloma (MM) is a highly heterogeneous disease. Cuproptosis is a novel mode of death that is closely associated with several diseases, such as hepatocellular carcinoma. However, its role in MM is unknown.

METHODS

MM transcriptomic and clinical data were obtained from UCSC Xena and gene expression omnibus (GEO) databases. Following MM samples were divided into different subtypes based on the cuproptosis genes, the differentially expressed genes (DEGs) among different subtypes, namely, candidate cuproptosis related genes were analyzed by univariate Cox and least absolute shrinkage and selection operator (LASSO) regression to construct a cuproptosis-related risk model. After the independent prognostic analysis was performed, a nomogram was constructed. Finally, Functional enrichment analysis and immune infiltration analysis were performed in the high- and low-risk groups, potential therapeutic agents were then predicted.

RESULTS

The 784 MM samples in UCSC Xena cohorts were divided into three different subtypes, and 4 out of 346 candidate cuproptosis related genes, namely CDKN2A, BCL3, KCNA3 and TTC14 were used to construct a risk model. Risk score was considered a reliable independent prognostic factor for MM patients. It was investigated that the pathway of cell cycle was significantly enriched in the high-risk group. In addition, immune score, ESTIMATE score and cytolytic activity were significantly different between different risk groups, as well as 13 immune cells such as memory B cells. Nine drugs were predicted in our study.

CONCLUSION

A cuproptosis-related prognostic model was constructed, which may have a potential guiding role in the treatment of MM.

Diagnostic role and immune correlates of programmed cell death-related genes in hepatocellular carcinoma

Programmed cell death (PCD) is thought to have multiple roles in tumors. Here, the roles of PCD-related genes were comprehensively analyzed to evaluate their values in hepatocellular carcinoma (HCC) diagnosis and prognosis. Gene expression and single-cell data of HCC patients, and PCD-related genes were collected from public databases. The diagnostic and prognostic roles of differentially expressed PCD-related genes in HCC were explored by univariate and multivariate Cox regression analyses. Single-cell data were further analyzed for the immune cells and expression of feature genes. Finally, we evaluated the expression of genes by quantitative real-time polymerase chain reaction and Western blot, and the proportion of immune cells was detected by flow cytometry in HCC samples. We obtained 52 differentially expressed PCD-related genes in HCC, based on which the consensus clustering analysis cluster 2 was found to have a worse prognosis than cluster 1. Then 10 feature genes were identified using LASSO analysis, and programmed cell death index (PCDI) was calculated to divided HCC patients into high-PCDI and low-PCDI groups. Worse prognosis was observed in high-PCDI group. Cox regression analysis showed that PCDI is an independent prognostic risk factor for HCC patients. Additionally, SERPINE1 and G6PD of feature genes significantly affect patient survival. Macrophages and Tregs were significantly positively correlated with PCDI. G6PD mainly expressed in macrophages, SERPINE1 mainly expressed in fibroblast. The experimental results confirmed the high expression of SERPINE1 and G6PD in HCC compared with the control, and the infiltration level of macrophages and Treg in HCC was also obviously elevated. PCDI may be a new predictor for the diagnosis of patients with HCC. The association of SERPINE1 and G6PD with the immune environment will provide new clues for HCC therapy.

Identification of MARK2, CCDC71, GATA2, and KLRC3 as candidate diagnostic genes and potential therapeutic targets for repeated implantation failure with antiphospholipid syndrome by integrated bioinformatics analysis and machine learning

Background

Antiphospholipid syndrome (APS) is a group of clinical syndromes of thrombosis or adverse pregnancy outcomes caused by antiphospholipid antibodies, which increase the incidence of in vitro fertilization failure in patients with infertility. However, the common mechanism of repeated implantation failure (RIF) with APS is unclear. This study aimed to search for potential diagnostic genes and potential therapeutic targets for RIF with APS.

Methods

To obtain differentially expressed genes (DEGs), we downloaded the APS and RIF datasets separately from the public Gene Expression Omnibus database and performed differential expression analysis. We then identified the common DEGs of APS and RIF. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were performed, and we then generated protein-protein interaction. Furthermore, immune infiltration was investigated by using the CIBERSORT algorithm on the APS and RIF datasets. LASSO regression analysis was used to screen for candidate diagnostic genes. To evaluate the diagnostic value, we developed a nomogram and validated it with receiver operating characteristic curves, then analyzed these genes in the Comparative Toxicogenomics Database. Finally, the Drug Gene Interaction Database was searched for potential therapeutic drugs, and the interactions between drugs, genes, and immune cells were depicted with a Sankey diagram.

Results

There were 11 common DEGs identified: four downregulated and seven upregulated. The common DEG analysis suggested that an imbalance of immune system-related cells and molecules may be a common feature in the pathophysiology of APS and RIF. Following validation, MARK2, CCDC71, GATA2, and KLRC3 were identified as candidate diagnostic genes. Finally, Acetaminophen and Fasudil were predicted as two candidate drugs.

Conclusion

Four immune-associated candidate diagnostic genes (MARK2, CCDC71, GATA2, and KLRC3) were identified, and a nomogram for RIF with APS diagnosis was developed. Our findings may aid in the investigation of potential biological mechanisms linking APS and RIF, as well as potential targets for diagnosis and treatment.

Identification of the role of immune-related genes in the diagnosis of bipolar disorder with metabolic syndrome through machine learning and comprehensive bioinformatics analysis

Background

Bipolar disorder and metabolic syndrome are both associated with the expression of immune disorders. The current study aims to find the effective diagnostic candidate genes for bipolar affective disorder with metabolic syndrome.

Methods

A validation data set of bipolar disorder and metabolic syndrome was provided by the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were found utilizing the Limma package, followed by weighted gene co-expression network analysis (WGCNA). Further analyses were performed to identify the key immune-related center genes through function enrichment analysis, followed by machine learning-based techniques for the construction of protein-protein interaction (PPI) network and identification of the Least Absolute Shrinkage and Selection Operator (LASSO) and Random Forest (RF). The receiver operating characteristic (ROC) curve was plotted to diagnose bipolar affective disorder with metabolic syndrome. To investigate the immune cell imbalance in bipolar disorder, the infiltration of the immune cells was developed.

Results

There were 2,289 DEGs in bipolar disorder, and 691 module genes in metabolic syndrome were identified. The DEGs of bipolar disorder and metabolic syndrome module genes crossed into 129 genes, so a total of 5 candidate genes were finally selected through machine learning. The ROC curve results-based assessment of the diagnostic value was done. These results suggest that these candidate genes have high diagnostic value.

Conclusion

Potential candidate genes for bipolar disorder with metabolic syndrome were found in 5 candidate genes (AP1G2, C1orf54, DMAC2L, RABEPK and ZFAND5), all of which have diagnostic significance.

Identification of co-diagnostic effect genes for aortic dissection and metabolic syndrome by multiple machine learning algorithms

Aortic dissection (AD) is a life-threatening condition in which the inner layer of the aorta tears. It has been reported that metabolic syndrome (MS) has a close linkage with aortic dissection. However, the inter-relational mechanisms between them were still unclear. This article explored the hub gene signatures and potential molecular mechanisms in AD and MS. We obtained five bulk RNA-seq datasets of AD, one single cell RNA-seq (scRNA-seq) dataset of ascending thoracic aortic aneurysm (ATAA), and one bulk RNA-seq dataset of MS from the gene expression omnibus (GEO) database. Identification of differentially expressed genes (DEGs) and key modules via weighted gene co-expression network analysis (WGCNA), functional enrichment analysis, and machine learning algorithms (Random Forest and LASSO regression) were used to identify hub genes for diagnosing AD with MS. XGBoost further improved the diagnostic performance of the model. The receiver operating characteristic (ROC) and precision-recall (PR) curves were developed to assess the diagnostic value. Then, immune cell infiltration and metabolism-associated pathways analyses were created to investigate immune cell and metabolism-associated pathway dysregulation in AD and MS. Finally, the scRNA-seq dataset was performed to confirm the expression levels of identified hub genes. 406 common DEGs were identified between the merged AD and MS datasets. Functional enrichment analysis revealed these DEGs were enriched for applicable terms of metabolism, cellular processes, organismal systems, and human diseases. Besides, the positively related key modules of AD and MS were mainly enriched in transcription factor binding and inflammatory response. In contrast, the negatively related modules were significantly associated with adaptive immune response and regulation of nuclease activity. Through machine learning, nine genes with common diagnostic effects were found in AD and MS, including MAD2L2, IMP4, PRPF4, CHSY1, SLC20A1, SLC9A1, TIPRL, DPYD, and MAPKAPK2. In the training set, the AUC of the hub gene on RP and RR curves was 1. In the AD verification set, the AUC of the Hub gene on RP and RR curves were 0.946 and 0.955, respectively. In the MS set, the AUC of the Hub gene on RP and RR curves were 0.978 and 0.98, respectively. scRNA-seq analysis revealed that the SLC20A1 was found to be relevant in fatty acid metabolic pathways and expressed in endothelial cells. Our study revealed the common pathogenesis of AD and MS. These common pathways and hub genes might provide new ideas for further mechanism research.

Identification and validation of immune-related hub genes based on machine learning in prostate cancer and AOX1 is an oxidative stress-related biomarker

To investigate potential diagnostic and prognostic biomarkers associated with prostate cancer (PCa), we obtained gene expression data from six datasets in the Gene Expression Omnibus (GEO) database. The datasets included 127 PCa cases and 52 normal controls. We filtered for differentially expressed genes (DEGs) and identified candidate PCa biomarkers using a least absolute shrinkage and selector operation (LASSO) regression model and support vector machine recursive feature elimination (SVM-RFE) analyses. A difference analysis was conducted on these genes in the test group. The discriminating ability of the train group was determined using the area under the receiver operating characteristic curve (AUC) value, with hub genes defined as those having an AUC greater than 85%. The expression levels and diagnostic utility of the biomarkers in PCa were further confirmed in the GSE69223 and GSE71016 datasets. Finally, the invasion of cells per sample was assessed using the CIBERSORT algorithm and the ESTIMATE technique. The possible prostate cancer (PCa) diagnostic biomarkers AOX1, APOC1, ARMCX1, FLRT3, GSTM2, and HPN were identified and validated using the GSE69223 and GSE71016 datasets. Among these biomarkers, AOX1 was found to be associated with oxidative stress and could potentially serve as a prognostic biomarker. Experimental validations showed that AOX1 expression was low in PCa cell lines. Overexpression of AOX1 significantly reduced the proliferation and migration of PCa cells, suggesting that the anti-tumor effect of AOX1 may be attributed to its impact on oxidative stress. Our study employed a comprehensive approach to identify PCa biomarkers and investigate the role of cell infiltration in PCa.

Bioinformatics analysis and prediction of Alzheimer's disease and alcohol dependence based on Ferroptosis-related genes

Background

Alzheimer's disease (AD) is a neurodegenerative disease whose origins have not been universally accepted. Numerous studies have demonstrated the relationship between AD and alcohol dependence; however, few studies have combined the origins of AD, alcohol dependence, and programmed cell death (PCD) to analyze the mechanistic relationship between the development of this pair of diseases. We demonstrated in previous studies the relationship between psychiatric disorders and PCD, and in the same concerning neurodegeneration-related AD, we found an interesting link with the Ferroptosis pathway. In the present study, we explored the bioinformatic interactions between AD, alcohol dependence, and Ferroptosis and tried to elucidate and predict the development of AD from this aspect.

Methods

We selected the Alzheimer's disease dataset GSE118553 and alcohol dependence dataset GSE44456 from the Gene Expression Omnibus (GEO) database. Ferroptosis-related genes were gathered through Gene Set Enrichment Analysis (GSEA), Kyoto Encyclopedia of Genes and Genomes (KEGG), and relevant literature, resulting in a total of 88 related genes. For the AD and alcohol dependence datasets, we conducted Limma analysis to identify differentially expressed genes (DEGs) and performed functional enrichment analysis on the intersection set. Furthermore, we used ferroptosis-related genes and the DEGs to perform machine learning crossover analysis, employing Least Absolute Shrinkage and Selection Operator (LASSO) regression to identify candidate immune-related central genes. This analysis was also used to construct protein-protein interaction networks (PPI) and artificial neural networks (ANN), as well as to plot receiver operating characteristic (ROC) curves for diagnosing AD and alcohol dependence. We analyzed immune cell infiltration to explore the role of immune cell dysregulation in AD. Subsequently, we conducted consensus clustering analysis of AD using three relevant candidate gene models and examined the immune microenvironment and functional pathways between different subgroups. Finally, we generated a network of gene-gene interactions and miRNA-gene interactions using Networkanalyst.

Results

The crossover of AD and alcohol dependence DEG contains 278 genes, and functional enrichment analysis showed that both AD and alcohol dependence were strongly correlated with Ferroptosis, and then crossed them with Ferroptosis-related genes to obtain seven genes. Three candidate genes were finally identified by machine learning to build a diagnostic prediction model. After validation by ANN and PPI analysis, ROC curves were plotted to assess the diagnostic value of AD and alcohol dependence. The results showed a high diagnostic value of the predictive model. In the immune infiltration analysis, functional metabolism and immune microenvironment of AD patients were significantly associated with Ferroptosis. Finally, analysis of target genes and miRNA-gene interaction networks showed that hsa-mir-34a-5p and has-mir-106b-5p could simultaneously regulate the expression of both CYBB and ACSL4.

Conclusion

We obtained a diagnostic prediction model with good effect by comprehensive analysis, and validation of ROC in AD and alcohol dependence data sets showed good diagnostic, predictive value for both AD (AUC 0. 75, CI 0.91-0.60), and alcohol dependence (AUC 0.81, CI 0.95-0.68). In the consensus clustering grouping, we identified variability in the metabolic and immune microenvironment between subgroups as a likely cause of the different prognosis, which was all related to Ferroptosis function. Finally, we discovered that hsa-mir-34a-5p and has-mir-106b-5p could simultaneously regulate the expression of both CYBB and ACSL4.

Identification of common genetic characteristics of rheumatoid arthritis and major depressive disorder by bioinformatics analysis and machine learning

Introduction

Depression is the most common comorbidity of rheumatoid arthritis (RA). In particular, major depressive disorder (MDD) and rheumatoid arthritis share highly overlapping mental and physical manifestations, such as depressed mood, sleep disturbance, fatigue, pain, and worthlessness. This overlap and indistinguishability often lead to the misattribution of physical and mental symptoms of RA patients to depression, and even, the depressive symptoms of MDD patients are ignored when receiving RA treatment. This has serious consequences, since the development of objective diagnostic tools to distinguish psychiatric symptoms from similar symptoms caused by physical diseases is urgent.

Methods

Bioinformatics analysis and machine learning.

Results

The common genetic characteristics of rheumatoid arthritis and major depressive disorder are EAF1, SDCBP and RNF19B.

Discussion

We discovered a connection between RA and MDD through immune infiltration studies: monocyte infiltration. Futhermore, we explored the correlation between the expression of the 3 marker genes and immune cell infiltration using the TIMER 2.0 database. This may help to explain the potential molecular mechanism by which RA and MDD increase the morbidity of each other.

Necroptosis-Related Prognostic Model for Pancreatic Carcinoma Reveals Its Invasion and Metastasis Potential through Hybrid EMT and Immune Escape

Necroptosis, pro-inflammatory programmed necrosis, has been reported to exert momentous roles in pancreatic cancer (PC). Herein, the objective of this study is to construct a necroptosis-related prognostic model for detecting pancreatic cancer. In this study, the intersection between necroptosis-related genes and differentially expressed genes (DEGs) of pancreatic ductal adenocarcinoma (PDAC) was obtained based on GeneCards database, GEO database (GSE28735 and GSE15471), and verified using The Cancer Genome Atlas (TCGA). Next, a prognostic model with Cox and LASSO regression analysis, and divided the patients into high-risk and low-risk groups. Subsequently, the Kaplan-Meier (KM) survival curve and the receiver operating characteristic (ROC) curves were generated to assess the predictive ability of overall survival (OS) of PC patients. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed to predict the potential biofunction and possible mechanical pathways. The EMTome database and an immune analysis were applied to further explore underlying mechanism. Finally, clinical samples of PDAC patients were utilized to verify the expression of model genes via immunohistochemistry (IHC), and the normal human pancreatic ductal cell line, hTERT-HPNE as well as human pancreatic ductal carcinoma cell lines, PANC-1 and PL45, were used to identify the levels of model genes by Western blot (WB) and immunofluorescence (IF) in vitro. The results showed that 13 necroptosis-related DEGs (NRDEGs) were screened based on GEO database, and finally four of five prognostic genes, including KRT7, KRT19, IGF2BP3, CXCL5, were further identified by TCGA to successfully construct a prognostic model. Univariate and multivariate Cox analysis ultimately confirmed that this prognostic model has independent prognostic significance, KM curve suggested that the OS of low-risk group was longer than high-risk group, and the area under receiver (AUC) of ROC for 1, 3, 5 years was 0.733, 0.749 and 0.667, respectively. A GO analysis illustrated that model genes may participate in cell-cell junction, cadherin binding, cell adhesion molecule binding, and neutrophil migration and chemotaxis, while KEGG showed involvement in PI3K-Akt signaling pathway, ECMreceptor interaction, IL-17 signaling pathway, TNF signaling pathway, etc. Moreover, our results showed KRT7 and KRT19 were closely related to EMT markers, and EMTome database manifested that KRT7 and KRT19 are highly expressed in both primary and metastatic pancreatic cancer, declaring that model genes promoted invasion and metastasis potential through EMT. In addition, four model genes were positively correlated with Th2, which has been reported to take part in promoting immune escape, while model genes except CXCL5 were negatively correlated with TFH cells, indicating that model genes may participate in immunity. Additionally, IHC results showed that model genes were higher expressed in PC tissues than that in adjacent tumor tissues, and WB and IF also suggested that model genes were more highly expressed in PANC-1 and PL45 than in hTERT-HPNE. Tracing of a necroptosis-related prognostic model for pancreatic carcinoma reveals its invasion and metastasis potential through EMT and immunity. The construction of this model and the possible mechanism of necroptosis in PDAC was preliminarily explored to provide reliable new biomarkers for the early diagnosis, treatment, and prognosis for pancreatic cancer patients.

Predictive model, miRNA-TF network, related subgroup identification and drug prediction of ischemic stroke complicated with mental disorders based on genes related to gut microbiome

Background

Patients with comorbid schizophrenia, depression, drug use, and multiple psychiatric diagnoses have a greater risk of carotid revascularization following stroke. The gut microbiome (GM) plays a crucial role in the attack of mental illness and IS, which may become an index for the diagnosis of IS. A genomic study of the genetic commonalities between SC and IS, as well as its mediated pathways and immune infiltration, will be conducted to determine how schizophrenia contributes to the high prevalence of IS. According to our study, this could be an indicator of ischemic stroke development.

Methods

We selected two datasets of IS from the Gene Expression Omnibus (GEO), one for training and the other for the verification group. Five genes related to mental disorders and GM were extracted from Gene cards and other databases. Linear models for microarray data (Limma) analysis was utilized to identify differentially expressed genes (DEGs) and perform functional enrichment analysis. It was also used to conduct machine learning exercises such as random forest and regression to identify the best candidate for immune-related central genes. Protein-protein interaction (PPI) network and artificial neural network (ANN) were established for verification. The receiver operating characteristic (ROC) curve was drawn for the diagnosis of IS, and the diagnostic model was verified by qRT-PCR. Further immune cell infiltration analysis was performed to study the IS immune cell imbalance. We also performed consensus clustering (CC) to analyze the expression of candidate models under different subtypes. Finally, miRNA, transcription factors (TFs), and drugs related to candidate genes were collected through the Network analyst online platform.

Results

Through comprehensive analysis, a diagnostic prediction model with good effect was obtained. Both the training group (AUC 0.82, CI 0.93-0.71) and the verification group (AUC 0.81, CI 0.90-0.72) had a good phenotype in the qRT-PCR test. And in verification group 2 we validated between the two groups with and without carotid-related ischemic cerebrovascular events (AUC 0.87, CI 1-0.64). Furthermore, we investigated cytokines in both GSEA and immune infiltration and verified cytokine-related responses by flow cytometry, particularly IL-6, which played an important role in IS occurrence and progression. Therefore, we speculate that mental illness may affect the development of IS in B cells and IL-6 in T cells. MiRNA (hsa-mir-129-2-3p, has-mir-335-5p, and has-mir-16-5p) and TFs (CREB1, FOXL1), which may be related to IS, were obtained.

Conclusion

Through comprehensive analysis, a diagnostic prediction model with good effect was obtained. Both the training group (AUC 0.82, CI 0.93-0.71) and the verification group (AUC 0.81, CI 0.90-0.72) had a good phenotype in the qRT-PCR test. And in verification group 2 we validated between the two groups with and without carotid-related ischemic cerebrovascular events (AUC 0.87, CI 1-0.64). MiRNA (hsa-mir-129-2-3p, has-mir-335-5p, and has-mir-16-5p) and TFs (CREB1, FOXL1), which may be related to IS, were obtained.

Exploring the shared gene signatures of smoking-related osteoporosis and chronic obstructive pulmonary disease using machine learning algorithms

Objectives: Cigarette smoking has been recognized as a predisposing factor for both osteoporosis (OP) and chronic obstructive pulmonary disease (COPD). This study aimed to investigate the shared gene signatures affected by cigarette smoking in OP and COPD through gene expression profiling. Materials and methods: Microarray datasets (GSE11784, GSE13850, GSE10006, and GSE103174) were obtained from Gene Expression Omnibus (GEO) and analyzed for differentially expressed genes (DEGs) and weighted gene co-expression network analysis (WGCNA). Least absolute shrinkage and selection operator (LASSO) regression method and a random forest (RF) machine learning algorithm were used to identify candidate biomarkers. The diagnostic value of the method was assessed using logistic regression and receiver operating characteristic (ROC) curve analysis. Finally, immune cell infiltration was analyzed to identify dysregulated immune cells in cigarette smoking-induced COPD. Results: In the smoking-related OP and COPD datasets, 2858 and 280 DEGs were identified, respectively. WGCNA revealed 982 genes strongly correlated with smoking-related OP, of which 32 overlapped with the hub genes of COPD. Gene Ontology (GO) enrichment analysis showed that the overlapping genes were enriched in the immune system category. Using LASSO regression and RF machine learning, six candidate genes were identified, and a logistic regression model was constructed, which had high diagnostic values for both the training set and external validation datasets. The area under the curves (AUCs) were 0.83 and 0.99, respectively. Immune cell infiltration analysis revealed dysregulation in several immune cells, and six immune-associated genes were identified for smoking-related OP and COPD, namely, mucosa-associated lymphoid tissue lymphoma translocation protein 1 (MALT1), tissue-type plasminogen activator (PLAT), sodium channel 1 subunit alpha (SCNN1A), sine oculis homeobox 3 (SIX3), sperm-associated antigen 9 (SPAG9), and vacuolar protein sorting 35 (VPS35). Conclusion: The findings suggest that immune cell infiltration profiles play a significant role in the shared pathogenesis of smoking-related OP and COPD. The results could provide valuable insights for developing novel therapeutic strategies for managing these disorders, as well as shedding light on their pathogenesis.

Schizophrenia and cell senescence candidate genes screening, machine learning, diagnostic models, and drug prediction

Background

Schizophrenia (SC) is one of the most common psychiatric diseases. Its potential pathogenic genes and effective treatment methods are still unclear. Cell senescence has been confirmed in mental diseases. A link exists between cellular senescence and immunity, and immune-related problems affect suicide rates in individuals suffering from schizophrenia. Therefore, the aims of this study were to identify candidate genes based on cell senescence that can affect the diagnosis and treatment of schizophrenia.

Methods

Two data sets of schizophrenia were provided by the Gene Expression Omnibus (GEO) database, one was taken as training and the other as a validation group. The genes related to cell senescence were obtained from the CellAge database. DEGs were identified using the Limma package and weighted gene co-expression network analysis (WGCNA). The function enrichment analysis was conducted, followed by machine learning-based identification for least absolute shrinking and selection operators (LASSO) regression. Random Forest were used to identify candidate immune-related central genes and establish artificial neural networks for verification of the candidate genes. The receiver operating characteristic curve (ROC curve) was used for the diagnosis of schizophrenia. Immune cell infiltrates were constructed to study immune cell dysregulation in schizophrenia, and relevant drugs with candidate genes were collected from the DrugBank database.

Results

Thirteen co-expression modules were screened for schizophrenia, of which 124 were the most relevant genes.There were 23 intersected genes of schizophrenia (including DEGs and the cellular senescence-related genes), and through machine learning six candidate genes were finally screened out. The diagnostic value was evaluated using the ROC curve data. Based on these results it was confirmed that these candidate genes have high diagnostic value.Two drugs related to candidate genes, Fostamatinib and Ritodine, were collected from the DrugBanks database.

Conclusion

Six potential candidate genes (SFN, KDM5B, MYLK, IRF3, IRF7, and ID1) had been identified, all of which had diagnostic significance. Fostamatinib might be a drug choice for patients with schizophrenia to develop immune thrombocytopenic purpura (ITP) after treatment, providing effective evidence for the pathogenesis and drug treatment of schizophrenia.

Application of machine learning in prediction of bone cement leakage during single-level thoracolumbar percutaneous vertebroplasty

BACKGROUND

In the elderly, osteoporotic vertebral compression fractures (OVCFs) of the thoracolumbar vertebra are common, and percutaneous vertebroplasty (PVP) is a common surgical method after fracture. Machine learning (ML) was used in this study to assist clinicians in preventing bone cement leakage during PVP surgery.

METHODS

The clinical data of 374 patients with thoracolumbar OVCFs who underwent single-level PVP at The First People's Hospital of Chenzhou were chosen. It included 150 patients with bone cement leakage and 224 patients without it. We screened the feature variables using four ML methods and used the intersection to generate the prediction model. In addition, predictive models were used in the validation cohort.

RESULTS

The ML method was used to select five factors to create a Nomogram diagnostic model. The nomogram model's AUC was 0.646667, and its C value was 0.647. The calibration curves revealed a consistent relationship between nomogram predictions and actual probabilities. In 91 randomized samples, the AUC of this nomogram model was 0.7555116.

CONCLUSION

In this study, we invented a prediction model for bone cement leakage in single-segment PVP surgery, which can help doctors in performing better surgery with reduced risk.

WGCNA combined with machine learning algorithms for analyzing key genes and immune cell infiltration in heart failure due to ischemic cardiomyopathy

Background

Ischemic cardiomyopathy (ICM) induced heart failure (HF) is one of the most common causes of death worldwide. This study aimed to find candidate genes for ICM-HF and to identify relevant biomarkers by machine learning (ML).

Methods

The expression data of ICM-HF and normal samples were downloaded from Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) between ICM-HF and normal group were identified. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment and gene ontology (GO) annotation analysis, protein–protein interaction (PPI) network, gene pathway enrichment analysis (GSEA), and single-sample gene set enrichment analysis (ssGSEA) were performed. Weighted gene co-expression network analysis (WGCNA) was applied to screen for disease-associated modules, and relevant genes were derived using four ML algorithms. The diagnostic values of candidate genes were assessed using receiver operating characteristic (ROC) curves. The immune cell infiltration analysis was performed between the ICM-HF and normal group. Validation was performed using another gene set.

Results

A total of 313 DEGs were identified between ICM-HF and normal group of GSE57345, which were mainly enriched in biological processes and pathways related to cell cycle regulation, lipid metabolism pathways, immune response pathways, and intrinsic organelle damage regulation. GSEA results showed positive correlations with pathways such as cholesterol metabolism in the ICM-HF group compared to normal group and lipid metabolism in adipocytes. GSEA results also showed a positive correlation with pathways such as cholesterol metabolism and a negative correlation with pathways such as lipolytic presentation in adipocytes compared to normal group. Combining multiple ML and cytohubba algorithms yielded 11 relevant genes. After validation using the GSE42955 validation sets, the 7 genes obtained by the machine learning algorithm were well verified. The immune cell infiltration analysis showed significant differences in mast cells, plasma cells, naive B cells, and NK cells.

Conclusion

Combined analysis using WGCNA and ML identified coiled-coil-helix-coiled-coil-helix domain containing 4 (CHCHD4), transmembrane protein 53 (TMEM53), acid phosphatase 3 (ACPP), aminoadipate-semialdehyde dehydrogenase (AASDH), purinergic receptor P2Y1 (P2RY1), caspase 3 (CASP3) and aquaporin 7 (AQP7) as potential biomarkers of ICM-HF. ICM-HF may be closely related to pathways such as mitochondrial damage and disorders of lipid metabolism, while the infiltration of multiple immune cells was identified to play a critical role in the progression of the disease.

Machine learning-based predictive models and drug prediction for schizophrenia in multiple programmed cell death patterns

Background

Schizophrenia (SC) is one of the most common mental illnesses. However, the underlying genes that cause it and its effective treatments are unknown. Programmed cell death (PCD) is associated with many immune diseases and plays an important role in schizophrenia, which may be a diagnostic indicator of the disease.

Methods

Two groups as training and validation groups were chosen for schizophrenia datasets from the Gene Expression Omnibus Database (GEO). Furthermore, the PCD-related genes of the 12 patterns were extracted from databases such as KEGG. Limma analysis was performed for differentially expressed genes (DEG) identification and functional enrichment analysis. Machine learning was employed to identify minimum absolute contractions and select operator (LASSO) regression to determine candidate immune-related center genes, construct protein-protein interaction networks (PPI), establish artificial neural networks (ANN), and validate with consensus clustering (CC) analysis, then Receiver operating characteristic curve (ROC curve) was drawn for diagnosis of schizophrenia. Immune cell infiltration was developed to investigate immune cell dysregulation in schizophrenia, and finally, related drugs with candidate genes were collected via the Network analyst online platform.

Results

In schizophrenia, 263 genes were crossed between DEG and PCD-related genes, and machine learning was used to select 42 candidate genes. Ten genes with the most significant differences were selected to establish a diagnostic prediction model by differential expression profiling. It was validated using artificial neural networks (ANN) and consensus clustering (CC), while ROC curves were plotted to assess diagnostic value. According to the findings, the predictive model had a high diagnostic value. Immune infiltration analysis revealed significant differences in Cytotoxic and NK cells in schizophrenia patients. Six candidate gene-related drugs were collected from the Network analyst online platform.

Conclusion

Our study systematically discovered 10 candidate hub genes (DPF2, ATG7, GSK3A, TFDP2, ACVR1, CX3CR1, AP4M1, DEPDC5, NR4A2, and IKBKB). A good diagnostic prediction model was obtained through comprehensive analysis in the training (AUC 0.91, CI 0.95-0.86) and validation group (AUC 0.94, CI 1.00-0.85). Furthermore, drugs that may be useful in the treatment of schizophrenia have been obtained (Valproic Acid, Epigallocatechin gallate).

The role of ferroptosis-related genes in airway epithelial cells of asthmatic patients based on bioinformatics

It has been reported that airway epithelial cells and ferroptosis have certain effect on asthma. However, the action mechanism of ferroptosis-related genes in airway epithelial cells of asthmatic patients is still unclear. Firstly, the study downloaded the GSE43696 training set, GSE63142 validation set and GSE164119 (miRNA) dataset from the gene expression omnibus database. 342 ferroptosis-related genes were downloaded from the ferroptosis database. Moreover, differentially expressed genes (DEGs) between asthma and control samples in the GSE43696 dataset were screened by differential analysis. Consensus clustering analysis was performed on asthma patients to classify clusters, and differential analysis was performed on clusters to obtain inter-cluster DEGs. Asthma-related module was screened by weighted gene co-expression network analysis. Then, DEGs between asthma and control samples, inter-cluster DEGs and asthma-related module were subjected to venn analysis for obtaining candidate genes. The last absolute shrinkage and selection operator and support vector machines were respectively applied to the candidate genes to screen for feature genes, and functional enrichment analysis was performed. Finally, a competition endogenetic RNA network was constructed and drug sensitivity analysis was conducted. There were 438 DEGs (183 up-regulated and 255 down-regulated) between asthma and control samples. 359 inter-cluster DEGs (158 up-regulated and 201 down-regulated) were obtained by screening. Then, the black module was significantly and strongly correlated with asthma. The venn analysis yielded 88 candidate genes. 9 feature genes (NAV3, ITGA10, SYT4, NOX1, SNTG2, RNF182, UPK1B, POSTN, SHISA2) were screened and they were involved in proteasome, dopaminergic synapse etc. Besides, 4 mRNAs, 5 miRNAs, and 2 lncRNAs collectively formed competition endogenetic RNA regulatory network, which included RNF182-hsa-miR-455-3p-LINC00319 and so on. The predicted therapeutic drug network map contained NAV3-bisphenol A and other relationship pairs. The study investigated the potential molecular mechanisms of NAV3, ITGA10, SYT4, NOX1, SNTG2, RNF182, UPK1B, POSTN, SHISA2 in airway epithelial cells of asthmatic patients through bioinformatics analysis, providing a reference for the research of asthma and ferroptosis.

Identification of clinical heterogeneity and construction of a novel subtype predictive model in patients with ankylosing spondylitis: An unsupervised machine learning study

BACKGROUND

Accurate classification of patients with ankylosing spondylitis (AS) is the premise of precision medicine so as to perform different medical interventions for different patient types. AS pathology is closely related to the changes in the immune microenvironment. In this study, we used unsupervised machine learning (UML) to classify patients with AS based on clinical characteristics. We then constructed a novel subtype predictive model for AS based on the clinical classification, after which we investigated the difference in the immune microenvironment to unravel the AS pathogenesis.

METHODS

Overall, 196 patients with AS were enrolled. UML was used to cluster AS patients by similar clinical characteristics. Functional ability, disease status, and grading of radiologic features were assessed to verify the accuracy and heterogeneity of UML clustering. Least Absolute Shrinkage and Selection Operator (LASSO) regression and Random Forest algorithm were used to screen and identify predictive factors for the novel subtype of AS. Logistic regression was also performed to construct a predictive model of this novel subtype. Datasets were downloaded from the Gene Expression Omnibus database to assess immune cell infiltration, and the results were validated using data of routine blood tests from 3671 AS patients and 5720 non-AS patients. The differential expression of Fat Mass and Obesity-Associated Protein (FTO), an m⁶A regulator, between AS patients and healthy control subjects was confirmed using immunohistochemistry.

RESULTS

UML clustering identified two clusters. The clinical characteristics of the two clusters were significantly heterogeneous. For the novel subtype of AS identified in UML clustering, a predictive model was built using three predictive factors, namely, C-reactive protein (CRP), absolute value of neutrophils (NEU), and absolute value of monocytes (MONO). The area under the curve of the predictive model was 0.983. Heterogeneity in the neutrophil and monocyte counts in AS was verified through immune cell infiltration analysis. Data from routine blood tests revealed that NEU and MONO were significantly higher in AS patients than in non-AS patients (p < 0.001). FTO expression was negatively correlated with both NEU and MONO. Immunohistochemistry analysis confirmed the downregulated expression of FTO.

CONCLUSIONS

UML provides an explicable and remarkable classification of a heterogeneous cohort of AS patients. A novel subtype of AS was identified in UML clustering. CRP, NEU, and MONO were the independent predictive factors for the novel subtype of AS. FTO expression was correlated with immune cell infiltration in AS patients.

Clinical diagnosis model of spinal meningiomas based on the surveillance, epidemiology, and end results database

Most spinal meningiomas (SM) are benign lesions of the thoracic spine and are usually treated surgically. This study aimed to explore treatment strategies and construct a nomogram for SM. Data on patients with SM from 2000 to 2019 were extracted from the Surveillance, Epidemiology, and End Results database. First, the distributional properties and characteristics of the patients were descriptively evaluated, and the patients were randomly divided into training and testing groups in a 6:4 ratio. Least absolute shrinkage and selection operator (LASSO) regression was used to screen the survival predictors. Kaplan-Meier curves explained survival probability by different variables. The nomogram was constructed based on the results of LASSO regression. The predictive power of the nomogram was identified using the concordance index, time-receiver operating characteristics, decision curve analysis, and calibration curves. We recruited 1,148 patients with SM. LASSO results for the training group showed that sex (coefficient, 0.004), age (coefficient, 0.034), surgery (coefficient, -0.474), tumor size (coefficient, 0.008), and marital status (coefficient, 0.335) were prognostic factors. The nomogram prognostic model showed good diagnostic ability in both the training and testing groups, with a C-index of 0.726, 95% (0.679, 0.773); 0.827, 95% (0.777, 0.877). The calibration and decision curves suggested that the prognostic model had better diagnostic performance and good clinical benefit. In the training and testing groups, the time-receiver operating characteristic curve showed that SM had moderate diagnostic ability at different times, and the survival rate of the high-risk group was significantly lower than that of the low-risk group (training group: p = 0.0071; testing group: p = 0.00013). Our nomogram prognostic model may have a crucial role in predicting the six-month, one-year, and two-year survival outcomes of patients with SM and may be useful for surgical clinicians to formulate treatment plans.

Extracellular vesicles-encapsulated microRNA in mammalian reproduction: A review

Extracellular vesicles (EVs) are nanoscale cell-derived lipid vesicles that participate in cell-cell communication by delivering cargo, including mRNAs, proteins and non-coding RNAs, to recipient cells. MicroRNA (miRNA), a non-coding RNA typically 22 nucleotides long, is crucial for nearly all developmental and pathophysiological processes in mammals by regulating recipient cells gene expression. Infertility is a worldwide health issue that affects 10-15% of couples during their reproductive years. Although assisted reproductive technology (ART) gives infertility couples hope, the failure of ART is mainly unknown. It is well accepted that EVs-encapsulated miRNAs have a role in different reproductive processes, implying that these EVs-encapsulated miRNAs could optimize ART, improve reproductive rate, and treat infertility. As a result, in this review, we describe the present understanding of EVs-encapsulated miRNAs in reproduction regulation.

Identification of immune-associated genes in diagnosing osteoarthritis with metabolic syndrome by integrated bioinformatics analysis and machine learning

Background

In the pathogenesis of osteoarthritis (OA) and metabolic syndrome (MetS), the immune system plays a particularly important role. The purpose of this study was to find key diagnostic candidate genes in OA patients who also had metabolic syndrome.

Methods

We searched the Gene Expression Omnibus (GEO) database for three OA and one MetS dataset. Limma, weighted gene co-expression network analysis (WGCNA), and machine learning algorithms were used to identify and analyze the immune genes associated with OA and MetS. They were evaluated using nomograms and receiver operating characteristic (ROC) curves, and finally, immune cells dysregulated in OA were investigated using immune infiltration analysis.

Results

After Limma analysis, the integrated OA dataset yielded 2263 DEGs, and the MetS dataset yielded the most relevant module containing 691 genes after WGCNA, with a total of 82 intersections between the two. The immune-related genes were mostly enriched in the enrichment analysis, and the immune infiltration analysis revealed an imbalance in multiple immune cells. Further machine learning screening yielded eight core genes that were evaluated by nomogram and diagnostic value and found to have a high diagnostic value (area under the curve from 0.82 to 0.96).

Conclusion

Eight immune-related core genes were identified (FZD7, IRAK3, KDELR3, PHC2, RHOB, RNF170, SOX13, and ZKSCAN4), and a nomogram for the diagnosis of OA and MetS was established. This research could lead to the identification of potential peripheral blood diagnostic candidate genes for MetS patients who also suffer from OA.

Identification and analysis of hub genes of hypoxia-immunity in type 2 diabetes mellitus

The chronic metabolic disease named type 2 diabetes (T2D) accounts for over 90% of diabetes mellitus. An increasing number of evidences have revealed that hypoxia has a significantly suppressive effect on cell-mediated immunity, as well as the utilization of glucose in diabetics. Therefore, we aimed to screen and identify hypoxia-immune-related hub genes in T2D through bioinformatic analysis. The Gene Expression Omnibus (GEO) database was used to get T2D gene expression profile data in the peripheral blood samples (GSE184050), and hypoxia-related genes were acquired from Molecular Signatures Database (MSigDB). Differentially expressed mRNAs (DEGs) and lncRNAs (DELs) between T2D and normal samples were identified by DeSeq2 package. The clusterProfiler package was used to perform enrichment analyses for the overlapped genes of DEGs and hypoxia-related genes. Further, Hypoxia-related hub genes were discovered using two machine learning algorithms. Next, the compositional patterns of immune and stromal cells in T2D and healthy groups were estimated by using xCell algorithm. Moreover, we used the weighted correlation network analysis (WGCNA) to examine the connection between genes and immune cells to screen immune-related genes. Gene Set Enrichment Analysis (GSEA) was used to investigate the functions of the hypoxia-immune-related hub genes. Finally, two peripheral blood cohorts of T2D (GSE184050 and GSE95849) as well as the quantitative real-time PCR (qRT-PCR) experiments for clicinal peripheral blood samples with T2D were used for verification analyses of hub genes. And meanwhile, a lncRNA-TF-mRNA network was constructed. Following the differentially expressed analysis, 38 out of 3822 DEGs were screened as hypoxia-related DEGs, and 493 DELs were found. These hypoxia-related DEGs were mainly enriched in the GO terms of pyruvate metabolic process, cytoplasmic vesicle lumen and monosaccharide binding, and the KEGG pathways of glycolysis/gluconeogenesis, pentose phosphate pathway and biosynthesis of nucleotide sugars. Moreover, 7 out of hypoxia-related DEGs were identified as hub genes. There were six differentially expressed immune cell types between T2D and healthy samples, which were further used as the clinical traits for WGCNA to identify AMPD3 and IER3 as the hypoxia-immune-related hub genes. The results of the KEGG pathways of genes in high-expression groups of AMPD3 and IER3 were mainly concentrated in glycosaminoglycan degradation and vasopressin-regulated water reabsorption, while the low-expression groups of AMPD3 and IER3 were mainly associated with RNA degradation and nucleotide excision repair. Finally, when compared to normal samples, both the AMPD3 and IER3 were highly expressed in the T2D groups in the GSE184050 and GSE95849 datasets. The result of lncRNA-TF-mRNA regulatory network showed that lncRNAs such as BACH1-IT1 and SNHG15 might induce the expression of the corresponding TFs such as TFAM and THAP12 and upregulate the expression of AMPD3. This study identified AMPD3 and IER3 as hypoxia-immune-related hub genes and potential regulatory mechanism for T2D, which provided a new perspective for elucidating the upstream molecular regulatory mechanism of diabetes mellitus.

Identification of prognostic values of the transcription factor-CpG-gene triplets in lung adenocarcinoma: A narrative review

OBJECTIVE

Abnormal DNA methylation can regulate carcinogenesis in lung adenocarcinoma (LUAD), while transcription factors (TFs) mediate methylation in a site-specific manner to affect downstream transcriptional regulation and tumor progression. Therefore, this study aimed to explore the TF-methylation-gene regulatory relationships that influence LUAD prognosis.

METHODS

Differential analyses of methylation sites and genes were generated by integrating transcriptome and methylome profiles from public databases. Through target gene identification, motif enrichment in the promoter region, and TF prediction, TF-methylation and methylation-gene relation pairs were obtained. Then, the prognostic TF-methylation-gene network was constructed using univariate Cox regression analysis. Prognostic models were constructed based on the key regulatory axes. Finally, Kaplan-Meier curves were created to evaluate the model efficacy and the relationship between candidate genes and prognosis.

RESULTS

A total of 1878 differential expressed genes and 1233 differential methylation sites were screened between LUAD and normal samples. Then 10 TFs were predicted to bind 144 enriched motifs. After integrating TF-methylation and methylation-gene relations, a prognostic TF-methylation-gene network containing 4 TFs, 111 methylation sites, and 177 genes was constructed. In this network, ERG-cg27071152-MTURN and FOXM1-cg19212949-PTPR regulatory axes were selected to construct the prognostic models, which showed robust abilities in predicting 1-, 3-, and 5-year survival probabilities. Finally, ERG and MTURN were downregulated in LUAD samples, whereas FOXM1 and PTPR were upregulated. Their expression levels were related to LUAD prognosis.

CONCLUSION

ERG-cg27071152-MTURN and FOXM1-cg19212949-PTPR regulatory axes were proposed as potential biomarkers for predicting the prognosis of LUAD.

Identification of cancer risk assessment signature in patients with chronic obstructive pulmonary disease and exploration of the potential key genes

It is essential to assess the cancer risk for patients with chronic obstructive pulmonary disease (COPD). Comparing gene expression data from patients with lung cancer (a total of 506 samples) and those with cancer-adjacent normal lung tissues (a total of 370 samples), we generated a qualitative transcriptional signature consisting of 2046 gene pairs. The signature was verified in an evaluation dataset comprising 18 subjects with severe disease and 52 subjects with moderate disease (Wilcoxon rank-sum test; p = 7.33 × 10^-5). Similar results were obtained in other independent datasets. Among the gene pairs in the signature, 326 COPD stage-related gene pairs were identified based on Spearman's rank correlation tests and those gene pairs comprised 368 unique genes. Of these 368 genes, 16 genes were significantly dysregulated in COPD rat model data compared with control data. Some of these genes (Dhx16, Upf2, Notch3, Sec61a1, Dyrk2, and Hmmr) were altered when the COPD rat model was treated with traditional Chinese medicines (TCM), including Bufei Yishen formula, Bufei Jianpi formula, and Yiqi Zishen formula. Overall, the signature could predict the cancer incidence-risk of COPD and the identified key genes might provide guidance regarding both the treatment of COPD using TCM and the prevention of cancer in patients with COPD. KEY MESSAGESA cancer risk assessment signature was identified in patients with COPD.The signature is insensitive to batch effects and is well verified.COPD key genes identified in this study might play a crucial role in TCM treatment and cancer prevention.

Positive regulators of T cell proliferation as biomarkers for predicting prognosis and characterizing the immune landscape in lung adenocarcinoma

Lung adenocarcinoma (LUAD) is the one of the most prevalent and fatal form of malignant tumors worldwide. Recently, immunotherapy is widely used in the treatment of patients with LUAD and has proved to be clinically effective in improve the prognosis of patients. But there still has been a tremendous thrust to further improve the efficacy of immunotherapy in individual patients with LUAD. The suppression of T cells and their effector functions in the tumor microenvironment (TME) of LUAD is one of the primary reasons for the low efficacy of immunotherapy in some patients with LUAD. Therefore, identifying positive regulators of T cell proliferation (TPRs) may offer novel avenues for LUAD immunotherapy. In this study, we comprehensively evaluated the infiltration patterns of TPRs in 1,066 patients with LUAD using unsupervised consensus clustering and identified correlations with genomic and clinicopathological characteristics. Three infiltrating TPR clusters were defined, and a TPR-related risk signature composed of nine TPRs was constructed using least absolute shrinkage and selection operator-Cox regression algorithms to classify the individual TPR infiltration patterns. Cluster 1 exhibited high levels of T cell infiltration and activation of immune-related signaling pathways, whereas cluster 2 was characterized by robust T cell immune infiltration and enrichment of pathways associated with carcinogenic gene sets and tumor immunity. Cluster 3 was characterized as an immune-desert phenotype. Moreover, the TPR signature was confirmed as an independent prognostic biomarker for drug sensitivity in patients with LUAD. In conclusion, the TPR signature may serve as a novel tool for effectively characterizing immune characteristics and evaluating the prognosis of patients with LUAD.

Comprehensive Analysis of GDF10 Methylation Site-Associated Genes as Prognostic Markers for Endometrial Cancer

Growth differentiation factor-10 (GDF10) with its methylation trait has recently been found to play a crucial regulatory and communication role in cancers. This investigation aims to identify GDF10 methylation site-associated genes that are closely associated with endometrial cancer (EC) patients' survival based on normal and UCEC samples from the UCSC Xena database. Our study revealed for the first time that EC exhibited significantly higher levels of GDF10 promoter methylation in comparison with normal tissues. Multiple differentiated methylation sites, which have prognostic value due to their apparent survival differences, were found in the GDF10 promoter region. We performed weighted gene coexpression network analysis (WGCNA) on EC tissues and paraneoplastic tissues while using these differentially methylated sites as phenotypes for selecting the most correlated key modules and their internal genes. To obtain a gene set, the key module genes and differentially expressed genes (DEGs) of EC were intersected. The least absolute shrinkage and selection operator (LASSO) regression along with multivariate Cox regression were performed from the gene set and we screened out the key genes B4GALNT3, DNAJC22, and GREB1. Finally, a prognostic model was validated for effectiveness based on these genes. Additionally, Kaplan-Meier analysis and time-dependent receiver operating characteristics (ROC) were applied to assess and verify the model, and they showed good prognosis prediction. Moreover, the differences in risk scores were statistically significant with age, tumor stage, and grade. They may be related to the immune infiltration of tumors as well. In conclusion, based on the methylation-related genes associated with GDF10, we developed a prognosis model for EC patients. It might provide a fresh view for further research and treatment of EC.

A novel epithelial-mesenchymal transition (EMT)-related gene signature of predictive value for the survival outcomes in lung adenocarcinoma

Lung adenocarcinoma (LUAD) is a remarkably heterogeneous and aggressive disease with dismal prognosis of patients. The identification of promising prognostic biomarkers might enable effective diagnosis and treatment of LUAD. Aberrant activation of epithelial-mesenchymal transition (EMT) is required for LUAD initiation, progression and metastasis. With the purpose of identifying a robust EMT-related gene signature (E-signature) to monitor the survival outcomes of LUAD patients. In The Cancer Genome Atlas (TCGA) database, least absolute shrinkage and selection operator (LASSO) analysis and cox regression analysis were conducted to acquire prognostic and EMT-related genes. A 4 EMT-related and prognostic gene signature, comprising dickkopf-like protein 1 (DKK1), lysyl oxidase-like 2 (LOXL2), matrix Gla protein (MGP) and slit guidance ligand 3 (SLIT3), was identified. By the usage of datum derived from TCGA database and Western blotting analysis, compared with adjacent tissue samples, DKK1 and LOXL2 protein expression in LUAD tissue samples were significantly higher, whereas the trend of MGP and SLIT3 expression were opposite. Concurrent with upregulation of epithelial markers and downregulation of mesenchymal markers, knockdown of DKK1 and LOXL2 impeded the migration and invasion of LUAD cells. Simultaneously, MGP and SLIT3 silencing promoted metastasis and induce EMT of LUAD cells. In the TCGA-LUAD set, receiver operating characteristic (ROC) analysis indicated that our risk model based on the identified E-signature was superior to those reported in literatures. Additionally, the E-signature carried robust prognostic significance. The validity of prediction in the E-signature was validated by the three independent datasets obtained from Gene Expression Omnibus (GEO) database. The probabilistic nomogram including the E-signature, pathological T stage and N stage was constructed and the nomogram demonstrated satisfactory discrimination and calibration. In LUAD patients, the E-signature risk score was associated with T stage, N stage, M stage and TNM stage. GSEA (gene set enrichment analysis) analysis indicated that the E-signature might be linked to the pathways including GLYCOLYSIS, MYC TARGETS, DNA REPAIR and so on. In conclusion, our study explored an innovative EMT based prognostic signature that might serve as a potential target for personalized and precision medicine.

Identification and Experimental Validation of Marker Genes between Diabetes and Alzheimer’s Disease

Currently, Alzheimer's disease (AD) and type 2 diabetes mellitus (T2DM) are widely prevalent in the elderly population, and accumulating evidence implies a strong link between them. For example, patients with T2DM have a higher risk of developing neurocognitive disorders, including AD, but the exact mechanisms are still unclear. This time, by combining bioinformatics analysis and in vivo experimental validation, we attempted to find a common biological link between AD and T2DM. We firstly downloaded the gene expression profiling (AD: GSE122063; T2DM: GSE161355) derived from the temporal cortex. To find the associations, differentially expressed genes (DEGs) of the two datasets were filtered and intersected. Based on them, enrichment analysis was carried out, and the least absolute shrinkage and selection operator (LASSO) logistic regression and support vector machine-recursive feature elimination (SVM-RFE) algorithms were used to identify the specific genes. After verifying in the external dataset and in the samples from the AD and type 2 diabetes animals, the shared targets of the two diseases were finally determined. Based on them, the ceRNA networks were constructed. Besides, the logistic regression and single-sample gene set enrichment analysis (ssGSEA) were performed. As a result, 62 DEGs were totally identified between AD and T2DM, and the enrichment analysis indicated that they were much related to the function of synaptic vesicle and MAPK signaling pathway. Based on the evidence from external dataset and RT-qPCR, CARTPT, EPHA5, and SERPINA3 were identified as the marker genes in both diseases, and their clinical significance and biological functions were further analyzed. In conclusion, discovering and exploring the marker genes that are dysregulated in both 2 diseases could help us better comprehend the intrinsic relationship between T2DM and AD, which may inspire us to develop new strategies for facing the dilemmas of clinical or basic research in cognitive dysfunction.

Development and Validation of a Machine Learning-Based Nomogram for Prediction of Ankylosing Spondylitis

Introduction

Ankylosing spondylitis (AS) is a chronic progressive inflammatory disease of the spine and its affiliated tissues. AS mainly affects the axial bone, sacroiliac joint, hip joint, spinal facet, and adjacent ligaments. We used machine learning (ML) methods to construct diagnostic models based on blood routine examination, liver function test, and kidney function test of patients with AS. This method will help clinicians enhance diagnostic efficiency and allow patients to receive systematic treatment as soon as possible.

Methods

We consecutively screened 348 patients with AS through complete blood routine examination, liver function test, and kidney function test at the First Affiliated Hospital of Guangxi Medical University according to the modified New York criteria (diagnostic criteria for AS). By using random sampling, the patients were randomly divided into training and validation cohorts. The training cohort included 258 patients with AS and 247 patients without AS, and the validation cohort included 90 patients with AS and 113 patients without AS. We used three ML methods (LASSO, random forest, and support vector machine recursive feature elimination) to screen feature variables and then took the intersection to obtain the prediction model. In addition, we used the prediction model on the validation cohort.

Results

Seven factors—erythrocyte sedimentation rate (ESR), red blood cell count (RBC), mean platelet volume (MPV), albumin (ALB), aspartate aminotransferase (AST), and creatinine (Cr)—were selected to construct a nomogram diagnostic model through ML. In the training cohort, the C value and area under the curve (AUC) value of this nomogram was 0.878 and 0.8779462, respectively. The C value and AUC value of the nomogram in the validation cohort was 0.823 and 0.8232055, respectively. Calibration curves in the training and validation cohorts showed satisfactory agreement between nomogram predictions and actual probabilities. The decision curve analysis showed that the nonadherence nomogram was clinically useful when intervention was decided at the nonadherence possibility threshold of 1%.

Conclusion

Our ML model can satisfactorily predict patients with AS. This nomogram can help orthopedic surgeons devise more personalized and rational clinical strategies.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40744-022-00481-6.

AS is a chronic progressive inflammatory disease of the spine and its affiliated tissues. AS starts gradually, and its early symptoms are mild. Some hospitals lack HLA-B27 and related imaging instruments to assist in the diagnosis of AS. There are relatively few studies on liver function and kidney function of patients with AS. We used ML methods to construct diagnostic models. Our model can satisfactorily predict patients with AS. This diagnostic model can help orthopedic surgeons devise more personalized and rational clinical strategies.

Assessing the Prognostic Capability of Immune-Related Gene Scoring Systems in Lung Adenocarcinoma

Background

Lung adenocarcinoma (LUAD) is the commonest of the subtypes of lung cancer histologically. For this study, we intended to analyze the expression profiling of the immune-related genes (IRGs) from an independently available public database and developed a potent signature predictive of patients' prognosis.

Methods

Gene expression profiles and the clinical data of lung adenocarcinoma were gathered from the Gene Expression Omnibus database (GEO) and The Cancer Genome Atlas (TCGA), and the obtained data were split into a training set (n = 226), test set (n = 83), and validation set (n = 400). IRGs were then gathered from the ImmPort database. A prognostic model was constructed by analyzing the training set. Then the GO and KEGG analysis was performed, and a gene correlation prognostic nomogram was constructed. Finally, external validation, such as immune infiltration and immunohistochemistry, was performed.

Results

The 110 genes were significant by univariate Cox regression analysis and randomized survival forest algorithm for the training set and showed a good distinction between the low-risk-score and high-risk-score groups in the training set (P < 0.0001) by screening for four prognosis-related genes (HMOX1, ARRB1, ADM, PDIA3) and validated by the test set GSE30219 (P=0.0025) and TCGA dataset (P=0.00059). Multivariate Cox showed that the four gene signatures were an individual risk factor for LUAD. In addition, the genes in the signatures were externally verified using an online database. In particular, PDIA3 and HMOX1 are essential genes in the prognostic nomogram and play an important role in the model of immune-related genes.

Conclusion

Four immune-related genetic signatures are reliable prognostic indicators for patients with LUAD, providing a relevant theoretical basis and therapeutic rationale for immunotherapy.

Identification of Immune-Associated Genes in Diagnosing Aortic Valve Calcification With Metabolic Syndrome by Integrated Bioinformatics Analysis and Machine Learning

Background

Immune system dysregulation plays a critical role in aortic valve calcification (AVC) and metabolic syndrome (MS) pathogenesis. The study aimed to identify pivotal diagnostic candidate genes for AVC patients with MS.

Methods

We obtained three AVC and one MS dataset from the gene expression omnibus (GEO) database. Identification of differentially expressed genes (DEGs) and module gene via Limma and weighted gene co-expression network analysis (WGCNA), functional enrichment analysis, protein–protein interaction (PPI) network construction, and machine learning algorithms (least absolute shrinkage and selection operator (LASSO) regression and random forest) were used to identify candidate immune-associated hub genes for diagnosing AVC with MS. To assess the diagnostic value, the nomogram and receiver operating characteristic (ROC) curve were developed. Finally, immune cell infiltration was created to investigate immune cell dysregulation in AVC.

Results

The merged AVC dataset included 587 DEGs, and 1,438 module genes were screened out in MS. MS DEGs were primarily enriched in immune regulation. The intersection of DEGs for AVC and module genes for MS was 50, which were mainly enriched in the immune system as well. Following the development of the PPI network, 26 node genes were filtered, and five candidate hub genes were chosen for nomogram building and diagnostic value evaluation after machine learning. The nomogram and all five candidate hub genes had high diagnostic values (area under the curve from 0.732 to 0.982). Various dysregulated immune cells were observed as well.

Conclusion

Five immune-associated candidate hub genes (BEX2, SPRY2, CXCL16, ITGAL, and MORF4L2) were identified, and the nomogram was constructed for AVC with MS diagnosis. Our study could provide potential peripheral blood diagnostic candidate genes for AVC in MS patients.