Dataset Information

Bioinformatics analysis and identification of hub genes and immune-related molecular mechanisms in chronic myeloid leukemia

ABSTRACT:

Background

Chronic myeloid leukemia (CML) is a malignant hyperplastic tumor of the bone marrow originating from pluripotent hematopoietic stem cells. The advent of tyrosine kinase inhibitors (TKIs) has greatly improved the survival rate of patients with CML. However, TKI-resistance leads to the disease recurrence and progression. This study aimed to identify immune-related genes (IRGs) associated with CML progression.

Methods

We extracted the gene’s expression profiles from the Gene Expression Omnibus (GEO). Bioinformatics analysis was used to determine the differentially expressed IRGs of CML and normal peripheral blood mononuclear cells (PBMCs). Functional enrichment and gene set enrichment analysis (GSEA) were used to explore its potential mechanism. Hub genes were identified using Molecular Complex Detection (MCODE) and the CytoHubba plugin. The hub genes’ diagnostic value was evaluated using the receiver operating characteristic (ROC). The relative proportions of infiltrating immune cells in each CML sample were evaluated using CIBERSORT. Quantitative real-time PCR (RT-qPCR) was used to validate the hub gene expression in clinical samples.

Results

A total of 31 differentially expressed IRGs were identified. GO analyses revealed that the modules were typically enriched in the receptor ligand activity, cytokine activity, and endopeptidase activity. KEGG enrichment analysis of IRGs revealed that CML involved Th17 cell differentiation, the NF-kappa B signaling pathway, and cytokine-cytokine receptor interaction. A total of 10 hub genes were selected using the PPI network. GSEA showed that these hub genes were related to the gamma-interferon immune response, inflammatory response, and allograft rejection. ROC curve analysis suggested that six hub genes may be potential biomarkers for CML diagnosis. Further analysis indicated that immune cells were associated with the pathogenesis of CML. The RT-qPCR results showed that proteinase 3 (PRTN3), cathepsin G (CTSG), matrix metalloproteinase 9 (MMP9), resistin (RETN), eosinophil derived neurotoxin (RNase2), eosinophil cationic protein (ECP, RNase3) were significantly elevated in CML patients’ PBMCs compared with healthy controls.

Conclusion

These results improved our understanding of the functional characteristics and immune-related molecular mechanisms involved in CML progression and provided potential diagnostic biomarkers and therapeutic targets.

SUBMITTER: Yao F

PROVIDER: S-EPMC8781323 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:PurposeFerroptosis and immune infiltration are involved in the pathogenesis of chronic obstructive pulmonary disease (COPD). We aim to identify ferroptosis-related hub genes and analyze their association with immune infiltration in COPD through bioinformatics methods.Materials and methodsThe mRNA microarray data of GSE38974 were downloaded from Gene Expression Omnibus to obtain differentially expressed genes (DEGs). The DEGs were intersected with ferroptosis-related genes (FRGs) from FerrDb to obtain differentially expressed FRGs. GO and KEGG enrichment and protein-protein interaction (PPI) analyses of differentially expressed FRGs were conducted in R software and STRING database. The key module and hub genes were screened by Cytoscape software. MiRNAs, transcription factors and signal molecules were predicted in miRNet and NetworkAnalyst. The disease correlation in the Comparative Toxicomics Database (CTD) and the receiver operating characteristic (ROC) curves of hub genes were analyzed. Immune infiltration was evaluated by CIBERSORT algorithm. Spearman correlation analyses were conducted between hub genes and differentially infiltrated immune cells.ResultsFifteen differentially expressed FRGs were identified, which were enriched in some terms involving airway inflammatory responses and structural remodeling. Five hub genes were selected including HIF1A, IL6, PTGS2, CDKN1A and ATM. Inference scores in CTD indicated their association with COPD. Two miRNAs, five transcription factors and one signal molecule were predicted. The combination of hub genes could be a fine diagnostic indicator of COPD (AUC: 0.981, CI: 0.940-1.000). Immune infiltration evaluation showed that monocytes and M0 macrophages were upregulated in COPD lung tissues, while CD8 T cells, activated NK cells, M2 macrophages, resting dendritic cells and resting mast cells were downregulated. The hub genes were significantly associated with differentially infiltrated immune cells.ConclusionWe identified five ferroptosis-related hub genes (HIF1A, IL6, PTGS2, CDKN1A and ATM) in COPD, and found that they may influence the pathogenesis of COPD by regulating ferroptosis and thus affecting infiltrating immune cells.

Project description:Alzheimer's disease (AD) is the leading cause of dementia in aged population. Oxidative stress and neuroinflammation play important roles in the pathogenesis of AD. Investigation of hub genes for the development of potential therapeutic targets and candidate biomarkers is warranted. The differentially expressed genes (DEGs) in AD were screened in GSE48350 dataset. The differentially expressed oxidative stress genes (DEOSGs) were analyzed by intersection of DEGs and oxidative stress-related genes. The immune-related DEOSGs and hub genes were identified by weighted gene co-expression network analysis (WGCNA) and protein-protein interaction (PPI) analysis, respectively. Enrichment analysis was performed by Gene Ontology and Kyoto Encyclopedia of Genes and Genomes. The diagnostic value of hub genes was assessed by receiver operating characteristic analysis and validated in GSE1297. The mRNA expression of diagnostic genes was determined by qRT-PCR analysis. Finally, we constructed the drug, transcription factors (TFs), and microRNA network of the diagnostic genes. A total of 1160 DEGs (259 up-regulated and 901 down-regulated) were screened in GSE48350. Among them 111 DEOSGs were identified in AD. Thereafter, we identified significant difference of infiltrated immune cells (effector memory CD8 T cell, activated B cell, memory B cell, natural killer cell, CD56 bright natural killer cell, natural killer T cell, plasmacytoid dendritic cell, and neutrophil) between AD and control samples. 27 gene modules were obtained through WGCNA and turquoise module was the most relevant module. We obtained 66 immune-related DEOSGs by intersecting turquoise module with the DEOSGs and identified 15 hub genes through PPI analysis. Among them, 9 hub genes (CCK, CNR1, GAD1, GAP43, NEFL, NPY, PENK, SST, and TAC1) were identified with good diagnostic values and verified in GSE1297. qRT-PCR analysis revealed the downregulation of SST, NPY, GAP43, CCK, and PENK and upregulation of NEFL in AD. Finally, we identified 76 therapeutic agents, 152 miRNAs targets, and 91 TFs regulatory networks. Our study identified 9 key genes associated with oxidative stress and immune reaction in AD pathogenesis. The findings may help to provide promising candidate biomarkers and therapeutic targets for AD.

Project description:ObjectiveGastric cancer (GC) is the fourth most common cause of cancer-related deaths in the world. In the current study, we aim to identify the hub genes and uncover the molecular mechanisms of GC.MethodsThe expression profiles of the genes and the miRNAs were extracted from the Gene Expression Omnibus database. The identification of the differentially expressed genes (DEGs), including miRNAs, was performed by the GEO2R. Database for Annotation, Visualization and Integrated Discovery was used to perform GO and KEGG pathway enrichment analysis. The protein-protein interaction (PPI) network and miRNA-gene network were constructed using Cytoscape software. The hub genes were identified by the Molecular Complex Detection (MCODE) plugin, the CytoHubba plugin and miRNA-gene network. Then, the identified genes were verified by Kaplan-Meier plotter database and quantitative real-time PCR (qRT-PCR) in GC tissue samples.ResultsA total of three mRNA expression profiles (GSE13911, GSE79973 and GSE19826) were downloaded from the Gene Expression Omnibus (GEO) database, including 69, 20 and 27cases separately. A total of 120 overlapped upregulated genes and 246 downregulated genes were identified. The majority of the DEGs were enriched in extracellular matrix organization, collagen catabolic process, collagen fibril organization and cell adhesion. In addition, three KEGG pathways were significantly enriched, including ECM-receptor interaction, protein digestion and absorption, and the focal adhesion pathways. In the PPI network, five significant modules were detected, while the genes in the modules were mainly involved in the ECM-receptor interaction and focal adhesion pathways. By combining the results of MCODE, CytoHubba and miRNA-gene network, a total of six hub genes including COL1A2, COL1A1, COL4A1, COL5A2, THBS2 and ITGA5 were chosen. The Kaplan-Meier plotter database confirmed that higher expression levels of these genes were related to lower overall survival, except for COL5A2. Experimental validation showed that the rest of the five genes had the same expression trend as predicted.ConclusionIn conclusion, COL1A2, COL1A1, COL4A1, THBS2 and ITGA5 may be potential biomarkers and therapeutic targets for GC. Moreover, ECM-receptor interaction and focal adhesion pathways play significant roles in the progression of GC.

Project description:Background: Acute Myeloid Leukemia (AML) is a complex and heterogeneous hematologic malignancy. However, the function of prognosis-related signature genes in AML remains unclear. Methods: In the current study, transcriptome sequencing was performed on 15 clinical samples, differentially expressed RNAs were identified using R software. The potential interactions network was constructed by using the common genes between target genes of differentially expressed miRNAs with transcriptome sequencing results. Functional and pathway enrichment analysis was performed to identify candidate gene-mediated aberrant signaling pathways. Hub genes were identified by the cytohubba plugin in Cytoscape software, which then expanded the potential interactions regulatory module for hub genes. TCGA-LAML clinical data were used for the prognostic analysis of the hub genes in the regulatory network, and GVSA analysis was used to identify the immune signature of prognosis-related hub genes. qRT-PCR was used to verify the expression of hub genes in independent clinical samples. Results: We obtained 1,610 differentially expressed lncRNAs, 233 differentially expressed miRNAs, and 2,217 differentially expressed mRNAs from transcriptome sequencing. The potential interactions network is constructed by 12 lncRNAs, 25 miRNAs, and 692 mRNAs. Subsequently, a sub-network including 15 miRNAs as well as 12 lncRNAs was created based on the expanded regulatory modules of 25 key genes. The prognostic analysis results show that CCL5 and lncRNA UCA1 was a significant impact on the prognosis of AML. Besides, we found three potential interactions networks such as lncRNA UCA1/hsa-miR-16-5p/COL4A5, lncRNA UCA1/hsa-miR-16-5p/SPARC, and lncRNA SNORA27/hsa-miR-17-5p/CCL5 may play an important role in AML. Furthermore, the evaluation of the immune infiltration shows that CCL5 is positively correlated with various immune signatures, and lncRNA UCA1 is negatively correlated with the immune signatures. Finally, the result of qRT-PCR showed that CCL5 is down-regulated and lncRNA UCA1 is up-regulated in AML samples separately. Conclusions: In conclusion, we propose that CCL5 and lncRNA UCA1 could be recognized biomarkers for predicting survival prognosis based on constructing competing endogenous RNAs in AML, which will provide us novel insight into developing novel prognostic, diagnostic, and therapeutic for AML.