Cell Type Identification By Estimating Relative Subsets Of RNA Transcripts Research Articles

The Prognostic Value and Immunotherapeutic Characteristics of GFPT2 in Pan-cancer.

The purpose of this study is to investigate the underlying relationship of diagnosis and therapy between glutamine-fructose-6-phosphate transaminase 2 (GFPT2) and various cancers. The Cancer Genome Atlas (TCGA) database was utilized to get gene expression RNAseq and clinical data for 33 tumors. The immunotherapeutic cohorts, including GSE35640, GSE78220, GSE67501, GSE181815, and IMvigor210, were derived from the Gene Expression Omnibus database (GEO) and a previously released article. Differential expression analysis of GFPT2 was performed using several clinical factors, and prognostic analysis was performed using Cox proportional hazard regression. In addition, the Cell type Identification By Estimating Relative Subsets Of RNA transcripts (CIBERSORT) and the Estimation of STromal and Immune cells in MAlignant Tumor tissues utilizing Expression data (ESTIMATE) algorithms were used to investigate the connection between GFPT2 and the tumor microenvironment. This approach additionally incorporated dynamic immunological indicators, such as tumor mutational burden (TMB) and microsatellite instability (MSI). In addition, a correlation between GFPT2 expression and the effectiveness of anticancer drugs was plotted for discussion. GFPT2 expression significantly differed in 11 out of 33 cancers. Although the distinct correlation between GFPT2 expression and clinical parameters had no wide distribution in pan-cancer, it demonstrated the potential prognostic validity of gene expression. GFPT2 demonstrated a strong correlation with immune infiltration, immune modulators, and immunerelated biomarkers. Furthermore, a variance analysis demonstrated a significant relationship between GFPT2 and the efficacy of immunotherapy. In addition, GFPT2 was associated with increased sensitivity of drugs such as Olaparib and Lenvatinib and decreased sensitivity of drugs such as Nilotinib. Collectively, GFPT2 is potentially useful as a biomarker for prognostic prediction and immune infiltration in a variety of malignancies ,and could lead to exciting new approaches to personalized oncotherapy.

Landscape of the immune infiltration and identification of molecular diagnostic markers associated with immune cells in patients with kidney transplantation

Rejection seriously affects the success of kidney transplantations. However, the molecular mechanisms underlying this rejection remain unclear. The GSE21374 and GSE36059 datasets were downloaded from the Gene Expression Omnibus (GEO) database. Next, the Cell-type Identification by Estimating Relative Subsets of RNA Transcripts (CIBERSORT) algorithm was used to infer the proportions of 22 immune cells. Moreover, infiltrating immune cell-related genes were identified using weighted gene co-expression network analysis (WGCNA), and enrichment analysis was conducted to observe their biological functions. Extreme Gradient Boosting (XGBoost) and Least Absolute Shrinkage and Selection Operator (LASSO) logistic regression algorithms were used to screen hub genes. Quantitative real-time PCR was conducted to verify the number of immune cells and hub gene expression levels. The rejection and non-rejection groups showed significantly different distributions (P < 0.05) of eight immune cells (B cell memory, Plasma cells, mast cells, follicular helper T cells, T CD8 cells, Macrophages M1, T Cells CD4 memory activated, and gamma delta T cells). Subsequently, CD8A, CRTAM, GBP2, WARS, and VAMP5 were screened as hub genes using the XGBoost and LASSO algorithms and could be used as diagnostic biomarkers. Finally, differential analysis and quantitative real-time PCR suggested that CD8A, CRTAM, GBP2, WARS, and VAMP5 were upregulated in rejection samples compared to non-rejection samples. The present study identified five key infiltrating immune cell-related genes (CD8A, CRTAM, GBP2,WARS, and VAMP5) involved in kidney transplant rejection, which may explain the molecular mechanism of rejection in kidney transplantation development.

Identification and validation of a novel innate lymphoid cell-based signature to predict prognosis and immune response in liver cancer by integrated single-cell RNA analysis and bulk RNA sequencing.

Innate lymphoid cells (ILCs) exert tumor suppressive and tumor promoting effects. However, the prognostic significance of ILC-associated genes remains unclear in hepatocellular carcinoma (HCC). Hence, the aim of this research was to develop an innovative predictive risk classification system using bioinformatics examination. We explored the Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) databases to gather data pertaining to HCC and its clinical details. Significantly different ILC-associated genes were investigated by Seurat analysis. The number of signaling interactions of ILCs with other cells was discovered by CellPhoneDB analysis. ClusterProfiler and Metascape were utilized to perform Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses on ILC genes. In order to identify potential ILC predictors, we utilized univariate Cox regression and least absolute shrinkage and selection operator (LASSO) analyses, subsequently validating these predictors in TCGA and GEO groups. The multi-omics ILC signature model's clinical predictive capabilities, along with drug sensitivity and immune factor relations, were assessed using Cell-type Identification by Estimating Relative Subsets of RNA Transcripts (CIBERSORT) and pRRophetic. We investigated the possible molecular pathways in our predictive ILC signature through the utilization of gene set enrichment analysis (GSEA) and gene set variation analysis (GSVA). Five key genes were screened out by constructing a competing endogenous RNA (ceRNA) network using Cytoscape and their values in clinical indexes were demonstrated. Immunohistochemistry (IHC) in HCC cases confirmed the expression of these genes. ILC cell subsets were identified, and cell-cell communication analysis revealed that the signaling pathways involving ILC cell subsets were mostly abundant in the HCC microenvironment. Subsequently, 270 marker genes of the ILC clusters were subjected to GO and KEGG enrichment analysis. Furthermore, a total of 58 prognostically relevant genes were screened as features for prognostic prediction models. Next, the models were validated and clinically evaluated (P values of Kaplan-Meier survival curves below 0.05). Five key genes (C2, STK4, CALM1, IL7R, and RORA) were further screened by multi omics analysis of immune cell and factor and drug sensitivity and correlation analysis of tumor regulatory genes in liver cancer. Furthermore, the potential clinical value of the 5 key genes was confirmed in HCC patients. Finally, the IHC results confirmed the expression of C2, STK4, CALM1, IL7R, and RORA in HCC. Our experimental results provided preliminary evidence supporting the oncogenic roles of STK4 and CALM1, as well as the tumor-suppressive roles of C2, RORA, and IL7R in HCC. A novel prognostic feature of ILC potentially involved in HCC was discovered. It showed high values in predicting patient overall survival (OS) as well as good differences in immunity and drug sensitivity. Therefore, targeting these ILC signatures may be a potential effective approach in HCC treatment.

Cuproptosis-related signature predicts prognosis and indicates tumor immune infiltration in bladder cancer.

Cuproptosis is a newly identified form of cell death that is dependent on copper (Cu) ions, termed Cu-dependent cytotoxicity. This process is distinct from other forms of cell death such as apoptosis, necrosis, and ferroptosis. The accumulation of copper is known to play a significant role in various biological processes, including angiogenesis (the formation of new blood vessels) and metastasis (the spread of cancer cells to different parts of the body). These processes are crucial for tumor growth and progression, indicating that copper and the cuproptosis-related genes (CPRGs) might be indispensable in the context of cancer development and progression. Given this background, we aimed to explore the relationship between CPRGs and both prognostic predictions and tumor microenvironment (TME) infiltration in bladder cancer (BLCA). For this study, we utilized data from The Cancer Genome Atlas (TCGA) to identify CPRGs and subsequently divided BLCA patients into three distinct molecular clusters based on these genes. To assess the proportions of various immune cell types within the TME, we employed single-sample gene set enrichment analysis (ssGSEA) and the Cell-type Identification by Estimating Relative Subsets of RNA Transcripts (CIBERSORT) method. These computational techniques allowed us to quantify the infiltration of different immune cells, providing insights into the immune landscape of the tumors. Furthermore, we developed a risk score model using CPRGs to predict the survival prospects of BLCA patients. Our analysis identified three molecular clusters of BLCA patients, each exhibiting unique clinical features and patterns of TME infiltration. Among these clusters, cluster 1 was associated with a poor prognosis. Interestingly, this cluster also showed significant infiltration of activated CD4+ (ssGSEA P<0.001) and CD8+ T (ssGSEA P<0.05) cells, which are crucial components of the immune response against tumors. This finding suggests a complex interaction between the immune system and the tumor, where a high presence of T cells does not necessarily correlate with better outcomes. Additionally, our risk score model revealed that the high-risk group, characterized by a specific expression pattern of CPRGs, also had enhanced infiltration of CD4+ and CD8+ T cells. This indicates that the cuproptosis-based risk model has a robust ability to predict patient prognosis and can guide immunotherapy decisions. Our study sheds light on the biological functions of CPRGs within the TME of BLCA and their correlations with clinical parameters and patient prognosis. The identification of distinct molecular clusters with varying prognoses and immune cell infiltrations highlights the heterogeneity of BLCA and underscores the potential of CPRGs as biomarkers for prognosis and therapeutic targets. These findings offer new perspectives for the development of immunotherapeutic strategies in the treatment of BLCA patients, potentially leading to more personalized and effective cancer therapies.

Two-step consensus clustering approach to immune cell infiltration: An integrated exploration and validation of prognostic and immune implications in sarcomas

To conduct a comprehensive investigation of the sarcoma immune cell infiltration (ImmCI) patterns and tumoral microenvironment (TME). We utilized transcriptomic, clinical, and mutation data of sarcoma patients (training cohort) obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) server. Cell-type Identification by Estimating Relative Subsets of RNA Transcripts (CIBERSORT) and Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data (ESTIMATE) algorithms were applied to decipher the immune cell infiltration landscape and TME profiles of sarcomas. An unsupervised clustering method was utilized for classifying ImmCI clusters (initial clustering) and ImmCI-based differentially expressed gene-driven clusters (secondary clustering). Mortality rates and immune checkpoint gene levels was analyzed among the identified clusters. We calculated the ImmCI score through principal component analysis. The tumor immune dysfunction evaluation (TIDE) score was also employed to quantify immunotherapy efficacy between two ImmCI score groups. We further validated the biomarkers for ImmCI and gene-driven clusters via experimental verification and the accuracy of the ImmCI score in predicting survival outcomes and immunotherapy efficacy by external validation cohorts (testing cohort). We demonstrated that ImmCI cluster A and gene-driven cluster A, were beneficial prognostic biomarkers and indicators of immune checkpoint blockade response in sarcomas via in-silico and laboratory experiments. Additionally, the ImmCI score exhibited independent prognostic significance and was predictive of immunotherapy response. Our research underscores the clinical significance of ImmCI scores in identifying sarcoma patients likely to respond to immunotherapy.

Construction of an artificial neural network diagnostic model and investigation of immune cell infiltration characteristics for idiopathic pulmonary fibrosis

BackgroundIdiopathic pulmonary fibrosis (IPF) is a severe lung condition, and finding better ways to diagnose and treat the disease is crucial for improving patient outcomes. Our study sought to develop an artificial neural network (ANN) model for IPF and determine the immune cell types that differed between the IPF and control groups.MethodsFrom the Gene Expression Omnibus (GEO) database, we first obtained IPF microarray datasets. To conduct protein-protein interaction (PPI) networks and enrichment analyses, differentially expressed genes (DEGs) were screened between tissues of patients with IPF and tissues of controls. Afterward, we identified the important feature genes associated with IPF using random forest (RF) analysis, and then constructed and validated a prediction ANN mode. In addition, the proportions of immune cells were quantified using cell-type identification by estimating relative subsets of RNA transcripts (CIBERSORT) analysis, which was performed on microarray datasets based on gene expression profiling.ResultsA total of 11 downregulated and 36 upregulated DEGs were identified. PPI networks and enrichment analyses were carried out; the immune system and extracellular matrix were the subjects of the enrichments. Using RF analysis, the significant feature genes LRRC17, COMP, ASPN, CRTAC1, POSTN, COL3A1, PEBP4, IL13RA2, and CA4 were identified. The nine feature gene scores were integrated into the ANN to develop a diagnostic prediction model. The receiver operating characteristic (ROC) curves demonstrated the strong diagnostic ability of the ANN in predicting IPF in the training and testing sets. An analysis of IPF tissues in comparison to normal tissues revealed a reduction in the infiltration of natural killer cells resting, monocytes, macrophages M0, and neutrophils; conversely, the infiltration of T cells CD4 memory resting, mast cells, and macrophages M0 increased.ConclusionLRRC17, COMP, ASPN, CRTAC1, POSTN, COL3A1, PEBP4, IL13RA2, and CA4 were determined as key feature genes for IPF. The nine feature genes in the ANN model will be extremely important for diagnosing IPF. It may be possible to use differentiated immune cells from IPF samples in comparison to normal samples as targets for immunotherapy in patients with IPF.

Integrated transcriptomic analysis and machine learning for characterizing diagnostic biomarkers and immune cell infiltration in fetal growth restriction.

Fetal growth restriction (FGR) occurs in 10% of pregnancies worldwide. Placenta dysfunction, as one of the most common causes of FGR, is associated with various poor perinatal outcomes. The main objectives of this study were to screen potential diagnostic biomarkers for FGR and to evaluate the function of immune cell infiltration in the process of FGR. Firstly, differential expression genes (DEGs) were identified in two Gene Expression Omnibus (GEO) datasets, and gene set enrichment analysis was performed. Diagnosis-related key genes were identified by using three machine learning algorithms (least absolute shrinkage and selection operator, random forest, and support vector machine model), and the nomogram was then developed. The receiver operating characteristic curve, calibration curve, and decision curve analysis curve were used to verify the validity of the diagnostic model. Using cell-type identification by estimating relative subsets of RNA transcripts (CIBERSORT), the characteristics of immune cell infiltration in placental tissue of FGR were evaluated and the candidate key immune cells of FGR were screened. In addition, this study also validated the diagnostic efficacy of TREM1 in the real world and explored associations between TREM1 and various clinical features. By overlapping the genes selected by three machine learning algorithms, four key genes were identified from 290 DEGs, and the diagnostic model based on the key genes showed good predictive performance (AUC = 0.971). The analysis of immune cell infiltration indicated that a variety of immune cells may be involved in the development of FGR, and nine candidate key immune cells of FGR were screened. Results from real-world data further validated TREM1 as an effective diagnostic biomarker (AUC = 0.894) and TREM1 expression was associated with increased uterine artery PI (UtA-PI) (p-value = 0.029). Four candidate hub genes (SCD, SPINK1, TREM1, and HIST1H2BB) were identified, and the nomogram was constructed for FGR diagnosis. TREM1 was not only associated with a variety of key immune cells but also correlated with increased UtA-PI. The results of this study could provide some new clues for future research on the prediction and treatment of FGR.

Identification of potential diagnostic biomarkers and therapeutic targets in patients with hypoxia pulmonary hypertension

BackgroundPulmonary hypertension is a serious disease. Emerging studies have shown that M2 macrophages play an essential role in pulmonary hypertension; however, their mechanism of action is uncertain. MethodsFour GEO datasets were downloaded. The differentially expressed genes (DEGs) were obtained using the limma package. Simultaneously, the Cell-type Identification By Estimating Relative Subsets Of RNA Transcripts (CIBERSORT) algorithm and weighted gene co-expression network analysis (WGCNA) were used to get the information about M2 macrophage–related modules. Potential key genes were obtained by intersecting DEGs with M2 macrophage–related module genes (M2MRGs), and finally the area under the curve (AUC) was calculated. Rats were exposed to hypoxia condition (10 % O2) for 4 weeks to induce PH. Subsequently, potential key genes with AUC>0.7 were analyzed by quantitative real-time polymerase chain reaction and Western blot using normoxia and hypoxia rat lungs. We knocked down EPHA3 in Raw264.7 cells and detected the protein expression of M2 macrophage markers including arginase 1 (ARG1) and interleukin 10 (IL-10), phospho-protein kinase B (P-Akt), and protein kinase B (Akt) to explore the downstream pathways of EPHA3. ResultsSeven potential hub genes were detected by intersecting M2MRGs and DEGs. Six genes with AUC values above 0.7 were used for further exploration. The expression of EPHA3 mRNA and protein was significantly more upregulated in rats with hypoxia than in rats with normoxia. The expression levels of IL10, ARG1, and P-Akt/Akt decreased after knocking down EPHA3. ConclusionsThis study suggested that the activation of the P-Akt/Akt signaling pathway promoted by EPHA3 played an essential role in the progression of pulmonary hypertension.

Bioinformatics analysis of SH2D4A in glioblastoma multiforme to evaluate immune features and predict prognosis.

Glioblastoma multiforme (GBM) is the most common and aggressive primary brain cancer in adults. This study aimed to obtain data on immune cell infiltration based on public datasets and to examine the prognostic significance of SH2 domain containing 4A (SH2D4A) for GBM. SH2D4A expression in GBM was analyzed using a Tumor Immunity Estimation Resource (TIMER) 2.0 dataset, and a gene expression profile interaction analysis (GEPIA), and the results were validated by quantitative reverse transcription polymerase chain reaction (qRT-PCR). The Chinese Glioma Genome Atlas (CGGA) dataset was used to assess the effect of SH2D4A on GBM patient survival. The SH2D4A co-expression network of the LinkedOmics dataset and GeneMANIA dataset was also investigated. Least absolute shrinkage and selection operator (LASSO) regression models and a nomogram were constructed to assess the prognosis of GBM patients. A Gene Set Enrichment Analysis (GSEA) was performed using The Cancer Genome Atlas (TCGA) dataset to find functional differences. The relationship between SH2D4A expression and tumor-infiltrating immune cells was analyzed using xCELL, the Cell Type Identification by Estimating Relative Subsets of RNA Transcripts (CIBERSORT) algorithm, and the TIMER dataset. We discovered that SH2D4A expression was upregulated in GBM patients, and elevated SH2D4A expression was also substantially correlated with tumor grade. The survival curve analysis and multivariate Cox regression analysis showed that high SH2D4A expression was a significant independent predictor of poor overall survival (OS) in GBM patients. The immunoassay results suggested that altered SH2D4A expression may affect the immune infiltration of GBM tissues and thus the survival outcomes of GBM patients. In addition to being a possible prognostic marker and therapeutic target for GBM, SH2D4A may also accelerate the progression of GBM.

Key genes and immune infiltration patterns and the clinical implications in psoriasis patients.

Psoriasis is an immune-mediated skin disease, closely related to immune regulation. The aim was to understand the pathogenesis of psoriasis further, reveal potential therapeutic targets, and provide new clues for its diagnosis, treatment, and prevention. Expression profiling data were obtained from the Gene Expression Omnibus (GEO) database for skin tissues from healthy population and psoriasis patients. Differentially expressed genes (DEGs) were selected for Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Gene Set Enrichment Analysis (GSEA) analysis separately. Machine learning algorithms were used to obtain characteristic genes closely associated with psoriasis. Receiver operating characteristic (ROC) curve was used to assess the diagnostic value of the characteristic genes for psoriasis. The Cell-type Identification by Estimating Relative Subsets of RNA Transcripts (CIBERSORT) algorithm was used to calculate the proportion of immune cell infiltration. Correlation analysis was used to characterize the connection between gene expression and immune cell, Psoriasis Area and Severity Index (PASI). A total of 254 DEGs were identified in the psoriasis group, including 185 upregulated and 69 downregulated genes. GO was mainly enriched in cytokine-mediated signaling pathway, response to virus, and cytokine activity. KEGG was mainly focused on cytokine-cytokine receptor interaction and IL-17 signaling pathway. GSEA was mainly in chemokine signaling pathway and cytokine-cytokine receptor interaction. The machine learning algorithm screened nine characteristic genes C10orf99, GDA, FCHSD1, C12orf56, S100A7, INA, CHRNA9, IFI44, and CXCL9. In the validation set, the expressions of these nine genes increased in the psoriasis group, and the AUC values were all>0.9, consistent with those of the training set. The immune infiltration results showed increased proportions of macrophages, T cells, and neutrophils in the psoriasis group. The characteristic genes were positively or negatively correlated to varying degrees with T cells and macrophages. Nine characteristic genes were highly expressed in the moderate to severe psoriasis group and positively correlated with PASI scores. High levels of nine characteristic genes C10orf99, GDA, FCHSD1, C12orf56, S100A7, INA, CHRNA9, IFI44, and CXCL9 were risk factors for psoriasis, the differential expression of which was related to the regulation of immune system activity and PASI scores, affecting the proportions of different immune cells and promoting the occurrence and development of psoriasis.

Development of a diagnostic model based on glycolysis-related genes and immune infiltration in intervertebral disc degeneration

BackgroundThe glycolytic pathway and immune response play pivotal roles in the intervertebral disc degeneration (IDD) progression. This study aimed to develop a glycolysis-related diagnostic model and analyze its relationship with the immune response to IDD. MethodsGSE70362, GSE23130, and GSE15227 datasets were collected and merged from the Gene Expression Omnibus, and differential expression analysis was performed. Glycolysis-related differentially expressed genes (GLRDEGs) were identified, and a machine learning-based diagnostic model was constructed and validated, followed by Gene Set Enrichment Analysis (GSEA). Gene Ontology functional enrichment and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were performed, and mRNA-miRNA and mRNA-transcription factor (TF) interaction networks were constructed. Immune infiltration was analyzed using single-sample GSEA (ssGSEA) and cell-type identification by estimating relative subsets of RNA transcripts (CIBERSORT) algorithm between high- and low-risk groups. ResultsIn the combined dataset, samples from 31 patients with IDD and 55 normal controls were analyzed, revealing differential expression of 16 GLRDEGs between the two groups. Using advanced machine learning techniques (LASSO, support vector machine, and random forest algorithms), we identified eight common GLRDEGs (PXK, EIF3D, WSB1, ZNF185, IGFBP3, CKAP4, RPL15, and, SSR1) and developed a diagnostic model, which demonstrated high accuracy in distinguishing IDD from control samples (area under the curve, 0.935). We identified 42 mRNA-miRNA and 33 mRNA-TF interaction pairs. Using the RiskScore from the diagnostic model, the combined dataset was stratified into high- and low-risk groups. SsGSEA revealed significant differences in the infiltration abundances of the four immune cell types between the groups. The CIBERSORT algorithm revealed the strongest correlation between resting natural killer (NK) cells and ZNF185 in the low-risk group and between CD8+ T cells and SSR1 in the high-risk group. ConclusionsOur study reveals a potential interplay between glycolysis-associated genes and immune infiltration in IDD pathogenesis. These findings contribute to our understanding of IDD and may guide development of novel diagnostic markers and therapeutic interventions.

Identification and experimental validation of KMO as a critical immune-associated mitochondrial gene in unstable atherosclerotic plaque

BackgroundThe heightened risk of cardiovascular and cerebrovascular events is associated with the increased instability of atherosclerotic plaques. However, the lack of effective diagnostic biomarkers has impeded the assessment of plaque instability currently. This study was aimed to investigate and identify hub genes associated with unstable plaques through the integration of various bioinformatics tools, providing novel insights into the detection and treatment of this condition.MethodsWeighted Gene Co-expression Network Analysis (WGCNA) combined with two machine learning methods were used to identify hub genes strongly associated with plaque instability. The cell-type identification by estimating relative subsets of RNA transcripts (CIBERSORT) method was utilized to assess immune cell infiltration patterns in atherosclerosis patients. Additionally, Gene Set Variation Analysis (GSVA) was conducted to investigate the potential biological functions, pathways, and mechanisms of hub genes associated with unstable plaques. To further validate the diagnostic efficiency and expression of the hub genes, immunohistochemistry (IHC), quantitative real-time polymerase chain reaction (RT-qPCR), and enzyme-linked immunosorbent assay (ELISA) were performed on collected human carotid plaque and blood samples. Immunofluorescence co-staining was also utilized to confirm the association between hub genes and immune cells, as well as their colocalization with mitochondria.ResultsThe CIBERSORT analysis demonstrated a significant decrease in the infiltration of CD8 T cells and an obvious increase in the infiltration of M0 macrophages in patients with atherosclerosis. Subsequently, two highly relevant modules (blue and green) strongly associated with atherosclerotic plaque instability were identified. Through intersection with mitochondria-related genes, 50 crucial genes were identified. Further analysis employing least absolute shrinkage and selection operator (LASSO) logistic regression and support vector machine recursive feature elimination (SVM-RFE) algorithms revealed six hub genes significantly associated with plaque instability. Among them, NT5DC3, ACADL, SLC25A4, ALDH1B1, and MAOB exhibited positive correlations with CD8 T cells and negative correlations with M0 macrophages, while kynurenine 3-monooxygenas (KMO) demonstrated a positive correlation with M0 macrophages and a negative correlation with CD8 T cells. IHC and RT-qPCR analyses of human carotid plaque samples, as well as ELISA analyses of blood samples, revealed significant upregulation of KMO and MAOB expression, along with decreased ALDH1B1 expression, in both stable and unstable samples compared to the control samples. However, among the three key genes mentioned above, only KMO showed a significant increase in expression in unstable plaque samples compared to stable plaque samples. Furthermore, the expression patterns of KMO in human carotid unstable plaque tissues and cultured mouse macrophage cell lines were assessed using immunofluorescence co-staining techniques. Finally, lentivirus-mediated KMO silencing was successfully transduced into the aortas of high-fat-fed ApoE-/- mice, with results indicating that KMO silencing attenuated plaque formation and promoted plaque stability in ApoE-/- mice.ConclusionsThe results suggest that KMO, a mitochondria-targeted gene associated with macrophage cells, holds promise as a valuable diagnostic biomarker for assessing the instability of atherosclerotic plaques.

Identification of a novel cellular senescence-related lncRNA signature for prognosis and immune response in osteosarcoma.

Cellular senescence, a novel hallmark of cancer, is associated with patient outcomes and tumor immunotherapy. However, at present, there is no systematic study on the use of cellular senescence-related long non-coding RNAs (CSR-lncRNAs) to predict survival in patients with osteosarcoma. In this study, we aimed to identify a CSR-lncRNAs signature and to evaluate its potential use as a survival prognostic marker and predictive tool for immune response of osteosarcoma. We downloaded a cohort of patients with osteosarcoma from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. We performed differential expression and co-expression analyses to identify CSR-lncRNAs. We performed univariate and multivariate Cox regression analyses along with the random forest algorithm to identify lncRNAs significantly correlated with senescence. Subsequently, we assessed the predictive models using survival curves, receiver operating characteristic curves, nomograms, C-index, and decision curve analysis. Based on this model, patients with osteosarcoma were divided into two groups according to their risk scores. Then, using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses, we compared their clinical characteristics to uncover functional differences. We further conducted immune infiltration analyses using estimation of stromal and immune cells in malignant tumor tissues using expression data (ESTIMATE), cell-type identification by estimating relative subsets of rna transcripts (CIBERSORT), and single-sample gene set enrichment analysis for the two groups. We also evaluated the expression of the target genes of immune checkpoint inhibitors (ICIs). We identified six lncRNAs that were significantly correlated with senescence and accordingly established a novel cellular senescence-related lncRNA prognostic signature incorporating these lncRNAs. The nomogram indicated that the risk model was an independent prognostic factor that could predict the survival of patients with osteosarcoma. This model demonstrated high accuracy upon validation. Further analysis revealed that patients with osteosarcoma in the low-risk group exhibited better clinical outcomes and enhanced immune infiltration. The six-CSR-lncRNA prognostic signature effectively predicted survival outcomes and patients in the low-risk group might have improved immune infiltration.

Integrative analysis of bone-formation associated genes and immune cell infiltration in osteoporosis, and the prediction of active ingredients in targeted traditional Chinese medicine

ObjectiveTo explore the differential expression and mechanisms of bone formation-related genes in osteoporosis (OP) leveraging bioinformatics and machine learning methodologies, and to predict the active ingredients of targeted traditional Chinese medicine (TCM) herbs. MethodsThe Gene Expression Omnibus (GEO) and GeneCards databases were employed to conduct a comprehensive screening of genes and disease-associated loci pertinent to the pathogenesis of OP. The R package was utilized as the analytical tool for the identification of differentially expressed genes. Least absolute shrinkage and selection operator (LASSO) logistic regression analysis and support vector machine-recursive feature elimination (SVM-RFE) algorithm were employed in defining the genetic signature specific to OP. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses for the selected pivotal genes were conducted. The cell-type identification by estimating relative subsets of RNA transcripts (CIBERSORT) algorithm was leveraged to examine the infiltration patterns of immune cells, with Spearman’s rank correlation analysis utilized to assess the relationship between the expression levels of the genes and the presence of immune cells. Coremine Medical Database was used to screen out potential TCM herbs for the treatment of OP. Comparative Toxicogenomics Database (CTD) was employed for forecasting the TCM active ingredients targeting the key genes. AutoDock Vina 1.2.2 and GROMACS 2020 softwares were employed to conclude analysis results, facilitating the exploration of binding affinity and conformational dynamics between the TCM active ingredients and their biological targets. ResultsTen genes were identified by intersecting the results from the GEO and GeneCards databases. Through the application of LASSO regression and SVM-RFE algorithm, four pivotal genes were selected: coat protein (CP), kallikrein 3 (KLK3), polymerase γ (POLG), and transient receptor potential vanilloid 4 (TRPV4). GO and KEGG pathway enrichment analyses revealed that these trait genes were predominantly engaged in the regulation of defense response activation, maintenance of cellular metal ion balance, and the production of chemokine ligand 5. These genes were notably associated with signaling pathways such as ferroptosis, porphyrin metabolism, and base excision repair. Immune infiltration analysis showed that key genes were highly correlated with immune cells. Macrophage M0, M1, M2, and resting dendritic cell were significantly different between groups, and there were significant differences between different groups (P < 0.05). The interaction counts of resveratrol, curcumin, and quercetin with KLK3 were 7, 3, and 2, respectively. It shows that the interactions of resveratrol, curcumin, and quercetin with KLK3 were substantial. Molecular docking and molecular dynamics simulations further confirmed the robust binding affinity of these bioactive compounds to the target genes. ConclusionPivotal genes including CP, KLK3, POLG, and TRPV4, exhibited commendable significant prognostic value, and played a crucial role in the diagnostic assessment of OP. Resveratrol, curcumin, and quercetin, natural compounds found in TCM, showed promise in their potential to effectively modulate the bone-forming gene KLK3. This study provides a scientific basis for the interpretation of the pathogenesis of OP and the development of clinical drugs.

Identification of Immune-Related Gene Signature in Schizophrenia.

Schizophrenia (SCZ) is a type of psychiatric disorder characterized by multiple symptoms. Our aim is to decipher the relevant mechanisms of immune-related gene signatures in SCZ. The SCZ dataset and its associated immunoregulatory genes were retrieved using Gene Expression Omnibus (GEO) and single-sample gene set enrichment analysis (ssGSEA). Co-expressed gene modules were determined through weighted gene correlation network analysis (WGCNA). To elucidate the functional characteristics of these clusters, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were used. Additionally, gene set enrichment analysis (GSEA) and Gene Set Variation Analysis (GSVA) were conducted to identify enriched pathways for the immune subgroups. A protein-protein interaction (PPI) network analysis was performed to identify core genes relevant to SCZ. A significantly higher immune score was observed in SCZ compared to control samples. Seven distinct gene modules were identified, with genes highlighted in green selected for further analysis. Using the Cell-type Identification By Estimating Relative Subsets Of RNA Transcripts (CIBERSORT) method, degrees of immune cell adhesion and accumulation related to 22 different immune cell types were calculated. Significantly enriched bioprocesses concerning the immunoregulatory genes with differential expressions included interferon-beta, IgG binding, and response to interferon-gamma, according to GO and KEGG analyses. Eleven hub genes related to immune infiltration emerged as key players among the three top-ranked GO terms. This study underscores the involvement of immunoregulatory reactions in SCZ development. Eleven immune-related genes (IFITM1 (interferon induced transmembrane protein 1), GBP1 (guanylate binding protein 1), BST2 (bone marrow stromal cell antigen 2), IFITM3 (interferon induced transmembrane protein 3), GBP2 (guanylate binding protein 2), CD44 (CD44 molecule), FCER1G (Fc epsilon receptor Ig), HLA-DRA (major histocompatibility complex, class II, DR alpha), FCGR2A (Fc gamma receptor IIa), IFI16 (interferon gamma inducible protein 16), and FCGR3B (Fc gamma receptor IIIb)) were identified as hub genes, representing potential biomarkers and therapeutic targets associated with the immune response in SCZ patients.

Integrated analyses reveal the diagnostic and predictive values of COL5A2 and association with immune environment in Crohn’s disease

The pathogenesis of Crohn’s disease (CD) involves abnormal immune cell infiltration and dysregulated immune response. Therefore, thorough research on immune cell abnormalities in CD is crucial for improved treatment of this disease. Single-cell RNA sequencing (scRNA-seq) and bulk RNA-seq data of CD were obtained from the Gene Expression Omnibus (GEO) database. Cell-type identification by estimating relative subsets of RNA transcripts (CIBERSORT), weighted gene co-expression network analysis (WGCNA), protein–protein interaction (PPI) networks evaluated the proportion of immune infiltrating cells, constructed co-expression network and identified key genes, respectively. Based on the dataset (GSE134809), 15 cell clusters were defined and labeled as different cell types. Among the 11 modules, the yellow module had the closest relationship with plasma cells (cluster 5). Confirmed using RNA sequencing and IHC assay, the expression of COL5A2 in CD samples was higher than that in control samples. Furthermore, the COL5A2 protein expression remarkably decreased in the group of patients who responded to anti–tumor necrosis factor (TNF) treatments, compared to the non-response group. The comprehensive analyses described here provided novel insight into the landscape of CD-associated immune environment. In addition, COL5A2 were identified as potential diagnostic indicators for CD, as well as promising predictive markers for CD patients.

Identification of red blood cell distribution width as a prognostic factor in acute myeloid leukemia

Many prognostic factors have been identified in acute myeloid leukemia (AML). In this study, we investigated novel prognostic biomarkers using machine learning and Cox regression models in a prospective cohort of 591 patients with AML and tried to identify potential therapeutic targets based on transcriptomic data. We found that elevated red blood cell distribution width (RDW) at diagnosis was an adverse prognostic factor for AML, independent of the 2022 European LeukemiaNet (ELN2022) genetic risk. As a continuous variable, higher RDW was associated with shorter overall survival (OS) (hazard ratio [HR] 1.087, 95% confidence interval [CI] 1.036-1.139, p < 0.001) and event-free survival (EFS) (HR 1.078, 95% CI 1.033-1.124, p < 0.001). Elevated RDW returned to normal after consolidation therapy, which indicated that leukemia cells resulted in abnormal RDW. We further investigated the relationship between RDW and transcriptome in another cohort of 191 patients with AML and public datasets using gene set enrichment analysis (GSEA) and cell-type identification by estimating relative subsets of RNA transcripts (CIBERSORT). We found that patients in the high-RDW group were significantly enriched in the positive regulation of erythroid differentiation and inflammation-related pathways. Finally, we identified the inflammation-associated gene IL12RB2 and verified its prognostic relevance with patients with AML in public databases, suggesting it as a potential therapy target.

Unraveling the copper-death connection: Decoding COVID-19‘s immune landscape through advanced bioinformatics and machine learning approaches

ABSTRACT This study aims to analyze Coronavirus Disease 2019 (COVID-19)-associated copper-death genes using the Gene Expression Omnibus (GEO) dataset and machine learning, exploring their immune microenvironment correlation and underlying mechanisms. Utilizing GEO, we analyzed the GSE217948 dataset with control samples. Differential expression analysis identified 16 differentially expressed copper-death genes, and Cell type Identification By Estimating Relative Subsets Of RNA Transcripts (CIBERSORT) quantified immune cell infiltration. Gene classification yielded two copper-death clusters, with Weighted Gene Co-expression Network Analysis (WGCNA) identifying key module genes. Machine learning models (random forest, Support Vector Machine (SVM), Generalized Linear Model (GLM), eXtreme Gradient Boosting (XGBoost)) selected 6 feature genes validated by the GSE213313 dataset. Ferredoxin 1 (FDX1) emerged as the top gene, corroborated by Area Under the Curve (AUC) analysis. Gene Set Enrichment Analysis (GSEA) and Gene Set Variation Analysis (GSVA) revealed enriched pathways in T cell receptor, natural killer cytotoxicity, and Peroxisome Proliferator-Activated Receptor (PPAR). We uncovered differentially expressed copper-death genes and immune infiltration differences, notably CD8 T cells and M0 macrophages. Clustering identified modules with potential implications for COVID-19. Machine learning models effectively predicted COVID-19 risk, with FDX1‘s pivotal role validated. FDX1‘s high expression was associated with immune pathways, suggesting its role in COVID-19 pathogenesis. This comprehensive approach elucidated COVID-19-related copper-death genes, their immune context, and risk prediction potential. FDX1‘s connection to immune pathways offers insights into COVID-19 mechanisms and therapy.

Predicting feature genes correlated with immune infiltration in patients with abdominal aortic aneurysm based on machine learning algorithms

Abdominal aortic aneurysm (AAA) is a condition characterized by a pathological and progressive dilatation of the infrarenal abdominal aorta. The exploration of AAA feature genes is crucial for enhancing the prognosis of AAA patients. Microarray datasets of AAA were downloaded from the Gene Expression Omnibus database. A total of 43 upregulated differentially expressed genes (DEGs) and 32 downregulated DEGs were obtained. Function, pathway, disease, and gene set enrichment analyses were performed, in which enrichments were related to inflammation and immune response. AHR, APLNR, ITGA10 and NR2F6 were defined as feature genes via machine learning algorithms and a validation cohort, which indicated high diagnostic abilities by the receiver operating characteristic curves. The cell-type identification by estimating relative subsets of RNA transcripts (CIBERSORT) method was used to quantify the proportions of immune infiltration in samples of AAA and normal tissues. We have predicted AHR, APLNR, ITGA10 and NR2F6 as feature genes of AAA. CD8 + T cells and M2 macrophages correlated with these genes may be involved in the development of AAA, which have the potential to be developed as risk predictors and immune interventions.

Immunogenic cell death related mRNAs associated signature to predict immunotherapeutic response in osteosarcoma

BackgroundImmunogenic cell death (ICD) is related to cancer prognosis, which has a synergic effect in combination with chemotherapy or immunotherapy. Yet, the relationship between ICD and osteosarcoma remained unclear. Materials and methodsThree osteosarcoma datasets including therapeutically applicable research to generate effective treatments (TARGET), GSE126209 and GSE21257 datasets were included. A protein-protein interaction network was constructed based on ICD-related genes. We performed unsupervised consensus clustering to classify molecular subtypes (clusters). Survival analysis, Estimation of stromal and immune cells in malignant tumour tissues using expression data (ESTIMATE), Cell-type identification by estimating relative subsets of RNA transcripts (CIBERSORT), and differential analysis were employed to characterize the molecular differences between different clusters. Univariate Cox regression analysis was conducted to confirm prognostic genes. Quantitative reverse transcription polymerase chain reaction (qRT-PCR) was used to demonstrate the aberrant expression of ICD-correlated signature genes in osteosarcoma. A series of cellular experiments, including cell counting kit-8 (CCK-8), transwell, and flow cytometry, were used to demonstrate the regulatory role of key genes in the ICD model on the malignant phenotype of osteosarcoma. ResultsThree clusters (cluster1, 2, 3) were constructed and they showed distinct overall survival and immune infiltration. ICD-related genes were highly expressed in cluster1. Moreover, Cluster1 had the best prognosis, high immune score and high expression of human leukocyte antigen (HLA)-related genes. TLR4, LY96, IFNGR1, CD4, and CASP1 were identified as prognostic genes for establishing an ICD-related risk signature. According to the risk signature, two risk groups (high and low risks) showing differential prognosis and response to immunotherapy. The low risks group had a better prognosis but was not sensitive to immunotherapy. Molecular assays verified that prognostic genes were abnormally under-expressed in osteosarcoma. Cellular assays demonstrated that LY96, the most significantly down-regulated gene in osteosarcoma, inhibited the migration, invasion, and proliferation phenotypes of osteosarcoma cells and prolonged the cell cycle. Analysis of oxidative stress related pathway enrichment in tumor microenvironment was conducted by single-sample gene set enrichment analysis (ssGSEA). ConclusionsThis study demonstrated the prognostic significance of ICD-correlated genes in osteosarcoma patients. The five-gene risk signature facilitate prognostic evaluation and prediction of osteosarcoma patients’ response to immunotherapy. The risk signature also offered a possibility for the exploit of novel ICD-related treatment.