Identification of Pathogenesis and Prognoses of Relapse in Myeloma by Bioinformatic Method

Author(s):  
Haoshu Zhong ◽  
Yang Liu ◽  
Jialin Duan ◽  
Xiaomin Chen ◽  
Hao Xiong ◽  
...  

Abstract Background: Multiple myeloma (MM), the second most hematological malignancy, the molecular mechanism and pathogenesis of the relapse of MM is poorly understood. This study aimed to identify novel prognostic model for MM and explore potential mechanism of relapse. Methods: Gene expression data,clinical data(GSE24080) and HTseq-Counts files were downloaded from Gene Expression Omnibus (GEO) and TCGA database. Co-expression modules of genes were built by Weighted Correlation Network Analysis (WGCNA).KEGG and GO enrichment analysis were performed in each module. TATFs (tumor-associated transcription factors) were retrieved from the Cistrome. Twenty-two immune cell compositions was calculated by CIBERSORT algorithm.Univariate and multivariate Cox congression were performed and a predictive model by prognostic genes was constructed,the predictive power of the model was evaluated by Kaplan–Meier curve and time-dependent receiver operating characteristic (ROC) curves. Results: A total of 940 DEGs were identified,and in WGCNA analysis, yellow, brown and sky-blue modules were most associated with clinic traits.The yellow module related with the cell cycle and the brown and sky-blue modules correlated with cytokine and its receptors, where the M2 macrophage fraction is positively correlated with CCL18, CCL2, CCL8, CXCL12 and CCl23 were positively correlated with plasma cells by Cibersort analysis.Prognostic genes were identified and two genes (TPX2,PRAM1) were finally identified to construct a risk model for predicting EFS.

2021 ◽  
Vol 12 ◽  
Author(s):  
Jixin Wang ◽  
Xiangjun Yin ◽  
Yin-Qiang Zhang ◽  
Xuming Ji

Lung adenocarcinoma (LUAD) is a major subtype of lung cancer, the prognosis of patients with which is associated with both lncRNAs and cancer immunity. In this study, we collected gene expression data of 585 LUAD patients from The Cancer Genome Atlas (TCGA) database and 605 subjects from the Gene Expression Omnibus (GEO) database. LUAD patients were divided into high and low immune-cell-infiltrated groups according to the single sample gene set enrichment analysis (ssGSEA) algorithm to identify differentially expressed genes (DEGs). Based on the 49 immune-related DE lncRNAs, a four-lncRNA prognostic signature was constructed by applying least absolute shrinkage and selection operator (LASSO) regression, univariate Cox regression, and stepwise multivariate Cox regression in sequence. Kaplan–Meier curve, ROC analysis, and the testing GEO datasets verified the effectiveness of the signature in predicting overall survival (OS). Univariate Cox regression and multivariate Cox regression suggested that the signature was an independent prognostic factor. The correlation analysis revealed that the infiltration immune cell subtypes were related to these lncRNAs.


2020 ◽  
Author(s):  
Hui Xie ◽  
Xiao-hui Ding ◽  
Ce Yuan ◽  
Jin-jiang Li ◽  
Zhao-yang Li ◽  
...  

Abstract Background: This study aimed to investigate the molecular mechanism of immune cell infiltration during the development of glioblastoma multiforme (GBM), and explore the potential immune cell associated prognostic genes for GBM. Methods: Gene expression data and corresponding clinical data of GBM samples (tumor group) and normal samples (normal group) in TCGA-GBM and GTEx dataset were downloaded. The differentially expression analysis was performed on samples between two groups. Based on tumor immune microenvironment analysis, the immune-related RNAs (lncRNAs and mRNAs) were further explored. Then, functional enrichment analysis, ceRNA network, risk prediction model and prognosis investigation were performed. Finally, the results of survival prognosis of key genes were tested by additional datasets. Results: A total of 4989 up-regulated genes and 2349 down-regulated genes were revealed between tumor group and normal group. M2 macrophages accounted for the largest proportion of tumor infiltrates immune cells in GBM, and was related to the prognosis of GBM patients. Totally 168 mRNAs (KIF18B) and 5 lncRNAs were related to infiltration of M2 Macrophage, of which 25 mRNAs and 5 lncRNAs forms a ceRNA network through 37 miRNAs (eg., miR-6849-3p). These genes were mainly assembled in functions like signal release. A risk model based on 5 mRNAs (such as FOX4 and ELFN2) and lncRNA PR11-161H23.5 was constructed. Verification test showed that all 5 mRNAs were significantly associated with OS prognosis.Conclusions: M2 Macrophage infiltration might participate in tumorigenesis of GBM via RP11-161H23.5-miR6849-3p-KIF18B ceRNA interaction. Furthermore, mRNAs such as FOX4 and ELFN2 might be potential prognostic markers for GBM patients.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Michael Kenn ◽  
Dan Cacsire Castillo-Tong ◽  
Christian F. Singer ◽  
Rudolf Karch ◽  
Michael Cibena ◽  
...  

AbstractCorrectly estimating the hormone receptor status for estrogen (ER) and progesterone (PGR) is crucial for precision therapy of breast cancer. It is known that conventional diagnostics (immunohistochemistry, IHC) yields a significant rate of wrongly diagnosed receptor status. Here we demonstrate how Dempster Shafer decision Theory (DST) enhances diagnostic precision by adding information from gene expression. We downloaded data of 3753 breast cancer patients from Gene Expression Omnibus. Information from IHC and gene expression was fused according to DST, and the clinical criterion for receptor positivity was re-modelled along DST. Receptor status predicted according to DST was compared with conventional assessment via IHC and gene-expression, and deviations were flagged as questionable. The survival of questionable cases turned out significantly worse (Kaplan Meier p < 1%) than for patients with receptor status confirmed by DST, indicating a substantial enhancement of diagnostic precision via DST. This study is not only relevant for precision medicine but also paves the way for introducing decision theory into OMICS data science.


2021 ◽  
Vol 104 (1) ◽  
pp. 003685042199727
Author(s):  
Xinyu Wang ◽  
Jiaojiao Yang ◽  
Xueren Gao

Lung adenocarcinoma (LUAD) is the most common histological type of lung cancer, comprising around 40% of all lung cancer. Until now, the pathogenesis of LUAD has not been fully elucidated. In the current study, we comprehensively analyzed the dysregulated genes in lung adenocarcinoma by mining public datasets. Two sets of gene expression datasets were obtained from the Gene Expression Omnibus (GEO) database. The dysregulated genes were identified by using the GEO2R online tool, and analyzed by R packages, Cytoscape software, STRING, and GPEIA online tools. A total of 275 common dysregulated genes were identified in two independent datasets, including 54 common up-regulated and 221 common down-regulated genes in LUAD. Gene Ontology (GO) enrichment analysis showed that these dysregulated genes were significantly enriched in 258 biological processes (BPs), 27 cellular components (CCs), and 21 molecular functions (MFs). Furthermore, protein-protein interaction (PPI) network analysis showed that PECAM1, ENG, KLF4, CDH5, and VWF were key genes. Survival analysis indicated that the low expression of ENG was associated with poor overall survival (OS) of LUAD patients. The low expression of PECAM1 was associated with poor OS and recurrence-free survival of LUAD patients. The cox regression model developed based on age, tumor stage, ENG, PECAM1 could effectively predict 5-year survival of LUAD patients. This study revealed some key genes, BPs, CCs, and MFs involved in LUAD, which would provide new insights into understanding the pathogenesis of LUAD. In addition, ENG and PECAM1 might serve as promising prognostic markers in LUAD.


2021 ◽  
Vol 9 (Suppl 3) ◽  
pp. A954-A955
Author(s):  
Jacob Kaufman ◽  
Doug Cress ◽  
Theresa Boyle ◽  
David Carbone ◽  
Neal Ready ◽  
...  

BackgroundLKB1 (STK11) is a commonly disrupted tumor suppressor in NSCLC. Its loss promotes an immune exclusion phenotype with evidence of low expression of interferon stimulated genes (ISG) and decreased microenvironment immune infiltration.1 2 Clinically, LKB1 loss induces primary immunotherapy resistance.3 LKB1 is a master regulator of a complex downstream kinase network and has pleiotropic effects on cell biology. Understanding the heterogeneous phenotypes associated with LKB1 loss and their influence on tumor-immune biology will help define and overcome mechanisms of immunotherapy resistance within this subset of lung cancer.MethodsWe applied multi-omic analyses across multiple lung adenocarcinoma datasets2 4–6 (>1000 tumors) to define transcriptional and genetic features enriched in LKB1-deficient lung cancer. Top scoring phenotypes exhibited heterogeneity across LKB1-loss tumors, and were further interrogated to determine association with increased or decreased markers of immune activity. Further, immune cell-types were estimated by Cibersort to identify effects of LKB1 loss on the immune microenvironment. Key conclusions were confirmed by blinded pathology review.ResultsWe show that LKB1 loss significantly affects differentiation patterns, with enrichment of ASCL1-expressing tumors with putative neuroendocrine differentiation. LKB1-deficient neuroendocrine tumors had lower expression of Interferon Stimulated Genes (ISG), MHC1 and MHC2 components, and immune infiltration compared to LKB1-WT and non-neuroendocrine LKB1-deficient tumors (figure 1).The abundances of 22 immune cell types assessed by Cibersort were compared between LKB1-deficient and LKB1-WT tumors. We observe skewing of immune microenvironmental composition by LKB1 loss, with lower abundance of dendritic cells, monocytes, and macrophages, and increased levels of neutrophils and plasma cells (table 1). These trends were most pronounced among tumors with neuroendocrine differentiation, and were concordant across three independent datasets. In a confirmatory subset of 20 tumors, plasma cell abundance was assessed by a blinded pathologist. Pathologist assessment was 100% concordant with Cibersort prediction, and association with LKB1 loss was confirmed (P=0.001).Abstract 909 Figure 1Immune-associated Gene Expression Profiles Affected by Neuroendocrine Differentiation within LKB1-Deficient Lung Adenocarcinomas. Gene expression profiles corresponding to five immune-associated phenotypes are shown with bars indicating average GEP scores for tumors grouped according to LKB1 and neuroendocrine status as indicated. P-values represent results from Student’s T-test between groups as indicated.Abstract 909 Table 1LKB1 Loss Affects Composition of Immune Microenvironment. Values indicate log10 P-values comparing LKB1-loss to LKB1-WT tumors. Positive (red) indicates increased abundance in LKB1 loss. Negative (blue) indicates decreased abundance.ConclusionsWe conclude that tumor differentiation patterns strongly influence the immune microenvironment and immune exclusion characteristics of LKB1-deficient tumors. Neuroendocrine differentiation is associated with the strongest immune exclusion characteristics and should be evaluated clinically for evidence of immunotherapy resistance. A novel observation of increased plasma cell abundance is observed across multiple datasets and confirmed by pathology. Causal mechanisms linking differentiation status to immune activity is not well understood, and the functional role of plasma cells in the immune biology of LKB1-deficient tumors is undefined. These questions warrant further study to inform precision immuno-oncology treatments for these patients.AcknowledgementsThis work was funded by SITC AZ Immunotherapy in Lung Cancer grant (SPS256666) and DOD Lung Cancer Research Program Concept Award (LC180633).ReferencesSkoulidis F, Byers LA, Diao L, et al. Co-occurring genomic alterations define major subsets of KRAS-mutant lung adenocarcinoma with distinct biology, immune profiles, and therapeutic vulnerabilities. Cancer Discov 2015;5:860–77.Schabath MB, Welsh EA, Fulp WJ, et al. Differential association of STK11 and TP53 with KRAS mutation-associated gene expression, proliferation and immune surveillance in lung adenocarcinoma. Oncogene 2016;35:3209–16.Skoulidis F, Goldberg ME, Greenawalt DM, et al. STK11/LKB1 mutations and PD-1 inhibitor resistance in KRAS-mutant lung adenocarcinoma. Cancer Discovery 2018;8:822-835.Cancer Genome Atlas Research Network. Comprehensive molecular profiling of lung adenocarcinoma. Nature 2014;511:543–50.Chitale D, Gong Y, Taylor BS, et al. An integrated genomic analysis of lung cancer reveals loss of DUSP4 in EGFR-mutant tumors. Oncogene 2009;28:2773–83.Shedden K, Taylor JM, Enkemann SA, et al. Gene expression-based survival prediction in lung adenocarcinoma: a multi-site, blinded validation study. Nat Med 2008;14:822–7.


2021 ◽  
Vol 15 ◽  
Author(s):  
Dezhi Shan ◽  
Xing Guo ◽  
Guozheng Yang ◽  
Zheng He ◽  
Rongrong Zhao ◽  
...  

Intracranial aneurysms (IAs) may cause lethal subarachnoid hemorrhage upon rupture, but the molecular mechanisms are poorly understood. The aims of this study were to analyze the transcriptional profiles to explore the functions and regulatory networks of differentially expressed genes (DEGs) in IA rupture by bioinformatics methods and to identify the underlying mechanisms. In this study, 1,471 DEGs were obtained, of which 619 were upregulated and 852 were downregulated. Gene enrichment analysis showed that the DEGs were mainly enriched in the inflammatory response, immune response, neutrophil chemotaxis, and macrophage differentiation. Related pathways include the regulation of actin cytoskeleton, leukocyte transendothelial migration, nuclear factor κB signaling pathway, Toll-like receptor signaling pathway, tumor necrosis factor signaling pathway, and chemokine signaling pathway. The enrichment analysis of 20 hub genes, subnetworks, and significant enrichment modules of weighted gene coexpression network analysis showed that the inflammatory response and immune response had a causal relationship with the rupture of unruptured IAs (UIAs). Next, the CIBERSORT method was used to analyze immune cell infiltration into ruptured IAs (RIAs) and UIAs. Macrophage infiltration into RIAs increased significantly compared with that into UIAs. The result of principal component analysis revealed that there was a difference between RIAs and UIAs in immune cell infiltration. A 4-gene immune-related risk model for IA rupture (IRMIR), containing CXCR4, CXCL3, CX3CL1, and CXCL16, was established using the glmnet package in R software. The receiver operating characteristic value revealed that the model represented an excellent clinical situation for potential application. Enzyme-linked immunosorbent assay was performed and showed that the concentrations of CXCR4 and CXCL3 in serum from RIA patients were significantly higher than those in serum from UIA patients. Finally, a competing endogenous RNA network was constructed to provide a potential explanation for the mechanism of immune cell infiltration into IAs. Our findings highlighted the importance of immune cell infiltration into RIAs, providing a direction for further research.


2021 ◽  
Vol 11 ◽  
Author(s):  
He Ren ◽  
Wanjing Li ◽  
Xin Liu ◽  
Shuliang Li ◽  
Hao Guo ◽  
...  

Hepatocellular carcinoma (HCC) is a common malignant tumor with relatively high malignancy and rapid disease progression. Metabolism-related genes (MRGs) are involved in the pathogenesis of HCC. This study explored potential key MRGs and their effect on T-cell immune function in the tumor immune microenvironment to provide new insight for the treatment of HCC. Of 456 differentially expressed MRGs identified from TCGA database, 21 were screened by MCODE and cytoHubba algorithms. From the key module, GAD1, SPP1, WFS1, GOT2, EHHADH, and APOA1 were selected for validation. The six MRGs were closely correlated with survival outcomes and clinicopathological characteristics in HCC. Receiver operating characteristics analysis and Kaplan-Meier plots showed that these genes had good prognostic value for HCC. Gene set enrichment analysis of the six MRGs indicated that they were associated with HCC development. TIMER and GEPIA databases revealed that WFS1 was significantly positively correlated and EHHADH was negatively correlated with tumor immune cell infiltration and immune checkpoint expression. Finally, quantificational real-time polymerase chain reaction (qRT-PCR) confirmed the expression of WFS1 and EHHADH mRNA in our own patients’ cohort samples and four HCC cell lines. Collectively, the present study identified six potential MRG biomarkers associated with the prognosis and tumor immune infiltration of HCC, thus providing new insight into the pathogenesis and treatment of HCC.


2021 ◽  
Vol 2021 ◽  
pp. 1-11
Author(s):  
Fanyan Meng ◽  
Ningna Du ◽  
Daoming Xu ◽  
Li Kuai ◽  
Lanying Liu ◽  
...  

Ankylosing spondylitis (AS) is an autoimmune disease that mainly affects the spinal joints, sacroiliac joints, and adjacent soft tissues. We conducted bioinformatics analysis to explore the molecular mechanism related to AS pathogenesis and uncover novel potential molecular targets for the treatment of AS. The profiles of GSE25101, containing gene expression data extracted from the blood of 16 AS patients and 16 matched controls, were acquired from the Gene Expression Omnibus (GEO) database. The background correction and standardization were carried out utilizing the transcript per million (TPM) method. After analysis of AS patients and the normal groups, we identified 199 differentially expressed genes (DEGs) with upregulation and 121 DEGs with downregulation by the limma R package. The results of the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway and Gene Ontology (GO) biological process enrichment analysis revealed that the DEGs with upregulation were mainly associated with spliceosome, ribosome, RNA-catabolic process, electron transport chain, etc. And the DEGs with downregulation primarily participated in T cell-associated pathways and processes. After analysis of the protein-protein interaction (PPI) network, our data revealed that the hub genes, comprising MRPL13, MRPL22, LSM3, COX7A2, COX7C, EP300, PTPRC, and CD4, could be the treatment targets in AS. Our data furnish new hints to uncover the features of AS and explore more promising treatment targets towards AS.


2021 ◽  
Author(s):  
Hongpeng Fang ◽  
Zhansen Huang ◽  
Xianzi Zeng ◽  
Jiaming Wan ◽  
Jieying Wu ◽  
...  

Abstract Background As a common malignant cancer of the urinary system, the precise molecular mechanisms of bladder cancer remain to be illuminated. The purpose of this study was to identify core genes with prognostic value as potential oncogenes for the diagnosis, prognosis or novel therapeutic targets of bladder cancer. Methods The gene expression profiles GSE3167 and GSE7476 were available from the Gene Expression Omnibus (GEO) database. Next, PPI network was built to filter the hub gene through the STRING database and Cytoscape software and GEPIA and Kaplan-Meier plotter were implemented. Frequency and type of hub genes and sub groups analysis were performed in cBioportal and ULCAN database. Finally,We used RT-qPCR to confirm our results. Results Totally, 251 DEGs were excavated from two datasets in our study. We only founded high expression of SMC4, TYMS, CCNB1, CKS1B, NUSAP1 and KPNA2 was associated with worse outcomes in bladder cancer patients and no matter from the type of mutation or at the transcriptional level of hub genes, the tumor showed a high form of expression. However, only the expression of SMC4,CCNB1and CKS1B remained changed between the cancer and the normal samples in our results of RT-qPCR. Conclusion In conclusion,These findings indicate that the SMC4,CCNB1 and CKS1B may serve as critical biomarkers in the development and poor prognosis.


Sign in / Sign up

Export Citation Format

Share Document