scholarly journals Comprehensive Analysis of the Tumor Microenvironment and Ferroptosis-Related Genes Predict Prognosis with Ovarian Cancer

2021 ◽  
Vol 12 ◽  
Author(s):  
Xiao-xue Li ◽  
Li Xiong ◽  
Yu Wen ◽  
Zi-jian Zhang

The early diagnosis of ovarian cancer (OC) is critical to improve the prognosis and prevent recurrence of patients. Nevertheless, there is still a lack of factors which can accurately predict it. In this study, we focused on the interaction of immune infiltration and ferroptosis and selected the ESTIMATE algorithm and 15 ferroptosis-related genes (FRGs) to construct a novel E-FRG scoring model for predicting overall survival of OC patients. The gene expression and corresponding clinical characteristics were obtained from the TCGA dataset (n = 375), GSE18520 (n = 53), and GSE32062 (n = 260). A total of 15 FRGs derived from FerrDb with the immune score and stromal score were identified in the prognostic model by using least absolute shrinkage and selection operator (LASSO)–penalized COX regression analysis. The Kaplan–Meier survival analysis and time-dependent ROC curves performed a powerful prognostic ability of the E-FRG model via multi-validation. Gene Set Enrichment Analysis and Gene Set Variation Analysis elucidate multiple potential pathways between the high and low E-FRG score group. Finally, the proteins of different genes in the model were verified in drug-resistant and non–drug-resistant tumor tissues. The results of this research provide new prospects in the role of immune infiltration and ferroptosis as a helpful tool to predict the outcome of OC patients.

2021 ◽  
Vol 12 ◽  
Author(s):  
Guomin Wu ◽  
Qihao Wang ◽  
Ting Zhu ◽  
Linhai Fu ◽  
Zhupeng Li ◽  
...  

This study aimed to establish a prognostic risk model for lung adenocarcinoma (LUAD). We firstly divided 535 LUAD samples in TCGA-LUAD into high-, medium-, and low-immune infiltration groups by consensus clustering analysis according to immunological competence assessment by single-sample gene set enrichment analysis (ssGSEA). Profile of long non-coding RNAs (lncRNAs) in normal samples and LUAD samples in TCGA was used for a differential expression analysis in the high- and low-immune infiltration groups. A total of 1,570 immune-related differential lncRNAs in LUAD were obtained by intersecting the above results. Afterward, univariate COX regression analysis and multivariate stepwise COX regression analysis were conducted to screen prognosis-related lncRNAs, and an eight-immune-related-lncRNA prognostic signature was finally acquired (AL365181.2, AC012213.4, DRAIC, MRGPRG-AS1, AP002478.1, AC092168.2, FAM30A, and LINC02412). Kaplan–Meier analysis and ROC analysis indicated that the eight-lncRNA-based model was accurate to predict the prognosis of LUAD patients. Simultaneously, univariate COX regression analysis and multivariate COX regression analysis were undertaken on clinical features and risk scores. It was illustrated that the risk score was a prognostic factor independent from clinical features. Moreover, immune data of LUAD in the TIMER database were analyzed. The eight-immune-related-lncRNA prognostic signature was related to the infiltration of B cells, CD4+ T cells, and dendritic cells. GSEA enrichment analysis revealed significant differences in high- and low-risk groups in pathways like pentose phosphate pathway, ubiquitin mediated proteolysis, and P53 signaling pathway. This study helps to treat LUAD patients and explore molecules related to LUAD immune infiltration to deeply understand the specific mechanism.


Author(s):  
Bo Xiao ◽  
Liyan Liu ◽  
Zhuoyuan Chen ◽  
Aoyu Li ◽  
Pingxiao Wang ◽  
...  

Melanoma is the most common cancer of the skin, associated with a worse prognosis and distant metastasis. Epithelial–mesenchymal transition (EMT) is a reversible cellular biological process that plays significant roles in diverse tumor functions, and it is modulated by specific genes and transcription factors. The relevance of EMT-related lncRNAs in melanoma has not been determined. Therefore, RNA expression data and clinical features were collected from the TCGA database (N = 447). Melanoma samples were randomly assigned into the training (315) and testing sets (132). An EMT-related lncRNA signature was constructed via comprehensive analyses of lncRNA expression level and corresponding clinical data. The Kaplan-Meier analysis showed significant differences in overall survival in patients with melanoma in the low and high-risk groups in two sets. Receiver operating characteristic (ROC) curves were used to measure the performance of the model. Cox regression analysis indicated that the risk score was an independent prognostic factor in two sets. Besides, a nomogram was constructed based on the independent variables. Gene Set Enrichment Analysis (GSEA) was applied to evaluate the potential biological functions in the two risk groups. Furthermore, the melanoma microenvironment was evaluated using ESTIMATE and CIBERSORT algorithms in the risk groups. This study indicates that EMT-related lncRNAs can function as potential independent prognostic biomarkers for melanoma survival.


2021 ◽  
Vol 19 (1) ◽  
pp. 169-190
Author(s):  
Peiyuan Li ◽  
◽  
Gangjie Qiao ◽  
Jian Lu ◽  
Wenbin Ji ◽  
...  

<abstract> <p>Plasmacytoma variant translocation 1 (PVT1) is involved in multiple signaling pathways and plays an important regulatory role in a variety of malignant tumors. However, its role in the prognosis and immune invasion of bladder urothelial carcinoma (BLCA) remains unclear. This study investigated the expression of PVT1 in tumor tissue and its relationship with immune invasion, and determined its prognostic role in patients with BLCA. Patients were identified from the cancer genome atlas (TCGA). The enrichment pathway and function of PVT1 were explained by gene ontology (GO) term analysis, gene set enrichment analysis (GSEA) and single-sample gene set enrichment analysis (ssGSEA), and the degree of immune cell infiltration was quantified. Kaplan–Meier analysis and Cox regression were used to analyze the correlation between PVT1 and survival rate. PVT1-high BLCA patients had a lower 10-year disease-specific survival (DSS P &lt; 0.05) and overall survival (OS P &lt; 0.05). Multivariate Cox regression analysis showed that PVT1 (high vs. low) (P = 0.004) was an independent prognostic factor. A nomogram was used to predict the effect of PVT1 on the prognosis. PVT1 plays an important role in the progression and prognosis of BLCA and can be used as a medium biomarker to predict survival after cystectomy.</p> </abstract>


2021 ◽  
Vol 11 ◽  
Author(s):  
Junyu Huo ◽  
Liqun Wu ◽  
Yunjin Zang

BackgroundThe high mutation rate of TP53 in hepatocellular carcinoma (HCC) makes it an attractive potential therapeutic target. However, the mechanism by which TP53 mutation affects the prognosis of HCC is not fully understood.Material and ApproachThis study downloaded a gene expression profile and clinical-related information from The Cancer Genome Atlas (TCGA) database and the international genome consortium (ICGC) database. We used Gene Set Enrichment Analysis (GSEA) to determine the difference in gene expression patterns between HCC samples with wild-type TP53 (n=258) and mutant TP53 (n=116) in the TCGA cohort. We screened prognosis-related genes by univariate Cox regression analysis and Kaplan–Meier (KM) survival analysis. We constructed a six-gene prognostic signature in the TCGA training group (n=184) by Lasso and multivariate Cox regression analysis. To assess the predictive capability and applicability of the signature in HCC, we conducted internal validation, external validation, integrated analysis and subgroup analysis.ResultsA prognostic signature consisting of six genes (EIF2S1, SEC61A1, CDC42EP2, SRM, GRM8, and TBCD) showed good performance in predicting the prognosis of HCC. The area under the curve (AUC) values of the ROC curve of 1-, 2-, and 3-year survival of the model were all greater than 0.7 in each independent cohort (internal testing cohort, n = 181; TCGA cohort, n = 365; ICGC cohort, n = 229; whole cohort, n = 594; subgroup, n = 9). Importantly, by gene set variation analysis (GSVA) and the single sample gene set enrichment analysis (ssGSEA) method, we found three possible causes that may lead to poor prognosis of HCC: high proliferative activity, low metabolic activity and immunosuppression.ConclusionOur study provides a reliable method for the prognostic risk assessment of HCC and has great potential for clinical transformation.


2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Yilin Lin ◽  
Xiaoxian Pan ◽  
Zhihua Chen ◽  
Suyong Lin ◽  
Zhanlong Shen ◽  
...  

Abstract Background Growing evidence has shown that the prognosis for colon cancer depends on changes in microenvironment. The purpose of this study was to elucidate the prognostic value of long noncoding RNAs (lncRNAs) related to immune microenvironment (IM) in colon cancer. Methods Single sample gene set enrichment analysis (ssGSEA) was used to identify the subtypes of colon cancer based on the immune genomes of 29 immune signatures. Cox regression analysis identified a lncRNA signatures associated with immune infiltration. The Tumor Immune Estimation Resource database was used to analyze immune cell content. Results Colon cancer samples were divided into three subtypes by unsupervised cluster analysis. Cox regression analysis identified an immune infiltration-related 5-lncRNA signature. This signature combined with clinical factors can effectively improve the predictive ability for the overall survival (OS) of colon cancer. At the same time, we found that the expression of H19 affects the content of B cells and macrophages in the microenvironment of colon cancer and affects the prognosis of colon cancer. Finally, we constructed the H19 regulatory network and further analyzed the possible mechanisms. We found that knocking down the expression of H19 can significantly inhibit the expression of CCND1 and VEGFA. At the same time, the immunohistochemical assay found that the expression of CCND1 and VEGFA protein was significantly positively correlated with the infiltration of M2 type macrophages. Conclusion The findings may help to formulate clinical strategies and understand the underlying mechanisms of H19 regulation. H19 may be a biomarker for targeted treatment of colon cancer.


Blood ◽  
2020 ◽  
Vol 136 (Supplement 1) ◽  
pp. 50-51
Author(s):  
Guillermo Montalban Bravo ◽  
Faezeh Darbaniyan ◽  
Rashmi Kanagal-Shamanna ◽  
Irene Ganan-Gomez ◽  
Koji Sasaki ◽  
...  

INTRODUCTION: Chronic myelomonocytic leukemia (CMML) is characterized by TET2, SRSF2, ASXL1 and RAS pathway mutations known to induce myelomonocytic bias. We have previously shown that upregulation of KDM6B, a histone demethylase that acts as an innate immune signal activator, leads to myeloid bias and expansion of immature myelomonocytic precursors and LSK cells in mice, resulting in rapid development of a myeloproliferative CMML (MP-CMML) phenotype particularly in cooperation with TET2 deletion. However, the role of genes involved in innate immunity regulation and monocyte differentiation in CMML phenotype and outcomes remains unclear. METHODS: We evaluated bone marrow aspirates from 19 patients with CMML and compared their transcriptomic features to those of healthy individuals obtained from AllCells (Emeryville, CA). CD34+ cells were isolated using the CD34 MicroBead Kit and RNA was isolated using the PicoPure RNA isolation kit. Fastq files were mapped to the human genome (build GRCh38) in TopHat2 using the default options. Differential gene expression analysis was conducted using DESeq2 in R version 3.4.2. Pathway enrichment analysis was performed using gene set enrichment analysis, with the fgsea library in R. Genes were ranked according to their Spearman correlation with the gene of interest, and this ranking was used as the input to fgsea. 10 000 gene permutations were used to calculate statistical significance, and a false discovery corrected p-value of 0.05 was required for statistical significance of a gene set. Cox regression and time ROC curves were used to study the relationship between gene expression and survival. We implemented Kaplan-Meier estimator along with optimum cutoff method to show the survival behavior in high versus low predicted model groups. RESULTS: Compared to healthy controls, a total of 1495 genes had significantly differential expression in CMML (q&lt;0.05, FC&gt;2) including 1271 genes which were significantly upregulated, and 224 which were significantly downregulated in CMML (Figure 1A). Gene set enrichment analysis identified 162 gene sets with differential expression in CMML compared to control (q&lt;0.05). Top upregulated genes were associated with interferon (IFN) alpha and beta signaling, chemokine receptors, IFN-gamma, GPC receptor ligand signaling and genes involved in immunomodulatory interactions between lymphoid and non-lymphoid cells (Figure 1B). Unsupervised clustering of gene expression profiles did not discriminate MP-CMML from myelodysplastic (MD-CMML). However, 20 genes were significantly overexpressed and 16 were significantly downregulated in patients with MP-CMML compared to MD-CMML (q&lt;0.05, FC&gt;2). In addition, 6 gene sets were differentially upregulated and 139 were significantly downregulated in pts with MP-CMML compared to MD-CMML (Figure 1C). To evaluate aberrant monopoiesis in CMML, we compared the expression of genes reported to be involved in regulation of monopoiesis among healthy controls and patients with CMML. A total of 23 genes involved in regulation of monopoiesis were found to be upregulated in CMML (Figure 1D). No significant differences in expression of these genes was found between MP-CMML and MD-CMML. To determine if aberrant expression of genes involved in monopoiesis influenced outcomes of pts with CMML, we developed a prediction model using Cox regression including 18 of these genes. Use of this model with optimum cutoffs allowed segregation of pts into two prognostic subsets with distinct survival outcomes (Figure 1E). Use of ROC curves identified high AUC particularly in pts with prolonged survival (&gt;40 months). CONCLUSIONS: CMML is characterized by upregulation of IFN and chemokine receptor signaling which could represent potential therapeutic targets. Aberrant expression of genes involved in regulation of monopoiesis may influence prognosis in CMML. Figure Disclosures Sasaki: Otsuka: Honoraria; Pfizer Japan: Consultancy; Daiichi Sankyo: Consultancy; Novartis: Consultancy, Research Funding. Kantarjian:Sanofi: Research Funding; Pfizer: Honoraria, Research Funding; Abbvie: Honoraria, Research Funding; Amgen: Honoraria, Research Funding; BMS: Research Funding; Daiichi-Sankyo: Honoraria, Research Funding; Novartis: Honoraria, Research Funding; Immunogen: Research Funding; Actinium: Honoraria, Membership on an entity's Board of Directors or advisory committees; Ascentage: Research Funding; Adaptive biotechnologies: Honoraria; Aptitute Health: Honoraria; BioAscend: Honoraria; Delta Fly: Honoraria; Janssen: Honoraria; Oxford Biomedical: Honoraria; Jazz: Research Funding. Garcia-Manero:AbbVie: Honoraria, Research Funding; Bristol-Myers Squibb: Consultancy, Research Funding; Genentech: Consultancy, Membership on an entity's Board of Directors or advisory committees, Research Funding; Novartis: Research Funding; Onconova: Research Funding; Acceleron Pharmaceuticals: Consultancy, Honoraria; H3 Biomedicine: Research Funding; Jazz Pharmaceuticals: Consultancy; Merck: Research Funding; Amphivena Therapeutics: Research Funding; Helsinn Therapeutics: Consultancy, Honoraria, Research Funding; Astex Pharmaceuticals: Consultancy, Honoraria, Research Funding; Celgene: Consultancy, Honoraria, Research Funding.


2021 ◽  
Author(s):  
Buze Chen ◽  
Xiaoyuan Lu ◽  
Qingmei Zhou ◽  
Qing Chen ◽  
Siyan Zhu ◽  
...  

Abstract Background: The long non-coding RNA (LncRNA) PAXIP1 antisense RNA 1 (PAXIP1-AS1) was found to promote proliferation, migration, EMT, and apoptosis of ovarian cancer (OC) cells in OC cell lines, but the relationship between PAXIP1-AS1 expression and clinical characteristics, prognosis, and immune infiltration of OC patients and its regulatory network are unclear. Methods: QRT-PCR, Kruskal-Wallis test, Wilcoxon sign-rank test, logistic regression, Kaplan-Meier method, Cox regression analysis, Gene set enrichment analysis (GSEA), and immuno-infiltration analysis were used to evaluate the relationship between clinical characteristics and PAXIP1-AS1 expression, prognostic factors, and determine the significant involvement of PAXIP1-AS1 in function. Results: Low PAXIP1-AS1 expression in OC was associated with age (P=0.045), histological grade (P=0.011), and lymphatic invasion (P=0.004). Low PAXIP1-AS1 expression predicted a poorer overall survival (OS) (HR: 0.71; 95% CI: 0.55–0.92; P=0.009), progression free interval (PFS) (HR: 1.776; 95% CI: 1.067–2.955; P=0.001) and disease specific survival (DSS) (HR: 0.67; 95% CI: 0.51–0.89; P=0.006). And PAXIP1-AS1 expression (HR: 0.711; 95% CI: 0.542-0.934; P=0.014) was independently correlated with PFS in OC patients. GSEA demonstrated that neutrophil degranulation, signaling by Interleukins, GPCR-ligand binding, G alpha I signaling events, VEGFAVEGFR-2 signaling pathway, naba secreted factors, Class A 1 Rhodopsin-Like Receptors, PI3K-Akt signaling pathway, and Focal Adhesion-PI3K-Akt-mTOR-signaling pathway were differentially enriched in PAXIP1-AS1 high expression phenotype. PAXIP1-AS1 may inhibit the function of aDC, B cells, CD8 T cells, Cytotoxic cells, DC, iDC, Macrophages, Mast cells, Neutrophils, NK CD56dim cells, T cells, TFH, Tgd, Th1 cells, Th2 cells and Treg. Conclusions: Low expression of PAXIP1-AS1 was significantly associated with poor survival and immune infiltration in OC. PAXIP1-AS1 could be a promising prognosis biomarker for OC.


2021 ◽  
Vol 11 ◽  
Author(s):  
Rui Liu ◽  
Ying Shen ◽  
Jinsong Hu ◽  
Xiaman Wang ◽  
Dong Wu ◽  
...  

BackgroundN6-methyladenosine is the most abundant RNA modification, which plays a prominent role in various biology processes, including tumorigenesis and immune regulation. Multiple myeloma (MM) is the second most frequent hematological malignancy.Materials and MethodsTwenty-two m6A RNA methylation regulators were analyzed between MM patients and normal samples. Kaplan–Meier survival analysis and least absolute shrinkage and selection operator (LASSO) Cox regression analysis were employed to construct the risk signature model. Receiver operation characteristic (ROC) curves were used to verify the prognostic and diagnostic efficiency. Immune infiltration level was evaluated by ESTIMATE algorithm and immune-related single-sample gene set enrichment analysis (ssGSEA).ResultsHigh expression of HNRNPC, HNRNPA2B1, and YTHDF2 and low expression of ZC3H13 were associated with poor survival. Based on these four genes, a prognostic risk signature model was established. Multivariate Cox regression analysis demonstrated that the risk score was an independent prognostic factor of MM. Enrichment analysis showed that cell cycle, immune response, MYC, proteasome, and unfold protein reaction were enriched in high-risk MM patients. Furthermore, patients with higher risk score exhibited lower immune scores and lower immune infiltration level.ConclusionThe m6A-based prognostic risk score accurately and robustly predicts the survival of MM patients and is associated with the immune infiltration level, which complements current prediction models and enhances our cognition of immune infiltration.


2021 ◽  
Vol 11 ◽  
Author(s):  
Sijin Sun ◽  
Kailun Fei ◽  
Guochao Zhang ◽  
Juhong Wang ◽  
Yannan Yang ◽  
...  

For lung adenocarcinoma (LUAD), patients of different stages have strong heterogeneity, and their overall prognosis varies greatly. Thus, exploration of novel biomarkers to better clarify the characteristics of LUAD is urgent. Multi-omics information of LUAD patients were collected form TCGA. Three independent LUAD cohorts were obtained from gene expression omnibus (GEO). A multi-omics correlation analysis of METTL5 was performed in TCGA dataset. To build a METTL5-associated prognostic score (MAPS). Spathial and random forest methods were first applied for feature selection. Then, LASSO was implemented to develop the model in TCGA cohort. The prognostic value of MAPS was validated in three independent GEO datasets. Finally, functional annotation was conducted using gene set enrichment analysis (GSEA) and the abundances of infiltrated immune cells were estimated by ImmuCellAI algorithm. A total of 901 LUAD patients were included. The expression of METTL5 in LUAD was significantly higher than that in normal lung tissue. And high expression of METTL5 indicated poor prognosis in all different stages (P &lt; 0.001, HR = 1.81). Five genes (RAC1, C11of24, METTL5, RCCD1, and SLC7A5) were used to construct MAPS and MAPS was significantly correlated with poor prognosis (P &lt; 0.001, HR = 2.15). Furthermore, multivariate Cox regression analysis suggested MAPS as an independent prognostic factor. Functional enrichment revealed significant association between MAPS and several immune components and pathways. This study provides insights into the potential significance of METTL5 in LUAD and MAPS can serve as a promising biomarker for LUAD.


2021 ◽  
Vol 27 ◽  
Author(s):  
Aoshuang Qi ◽  
Mingyi Ju ◽  
Yinfeng Liu ◽  
Jia Bi ◽  
Qian Wei ◽  
...  

Background: Complex antigen processing and presentation processes are involved in the development and progression of breast cancer (BC). A single biomarker is unlikely to adequately reflect the complex interplay between immune cells and cancer; however, there have been few attempts to find a robust antigen processing and presentation-related signature to predict the survival outcome of BC patients with respect to tumor immunology. Therefore, we aimed to develop an accurate gene signature based on immune-related genes for prognosis prediction of BC.Methods: Information on BC patients was obtained from The Cancer Genome Atlas. Gene set enrichment analysis was used to confirm the gene set related to antigen processing and presentation that contributed to BC. Cox proportional regression, multivariate Cox regression, and stratified analysis were used to identify the prognostic power of the gene signature. Differentially expressed mRNAs between high- and low-risk groups were determined by KEGG analysis.Results: A three-gene signature comprising HSPA5 (heat shock protein family A member 5), PSME2 (proteasome activator subunit 2), and HLA-F (major histocompatibility complex, class I, F) was significantly associated with OS. HSPA5 and PSME2 were protective (hazard ratio (HR) &lt; 1), and HLA-F was risky (HR &gt; 1). Risk score, estrogen receptor (ER), progesterone receptor (PR) and PD-L1 were independent prognostic indicators. KIT and ACACB may have important roles in the mechanism by which the gene signature regulates prognosis of BC.Conclusion: The proposed three-gene signature is a promising biomarker for estimating survival outcomes in BC patients.


Sign in / Sign up

Export Citation Format

Share Document