High expression of PIMREG predicts poor survival outcomes and is correlated with immune infiltrates in lung adenocarcinoma

Background PIMREG is upregulated in multiple cancer types. However, the potential role of PIMREG in lung adenocarcinoma (LUAD) remains unclear. The present study aimed to explore its clinical significance in LUAD. Methods Using the Cancer Genome Atlas (TCGA) databases, we obtained 513 samples of LUAD and 59 normal samples from the Cancer Genome Atlas (TCGA) databases to analyze the relationship between PIMREG and LUAD. We used t and Chi-square tests to evaluate the level of expression of PIMREG and its clinical implication in LUAD. The prognostic value of PIMREG in LUAD was identified through the Kaplan–Meier method, Cox regression analysis, and nomogram. Gene set enrichment analysis (GSEA) and single-sample gene set enrichment analysis (ssGSEA) were performed to screen biological pathways and analyze the correlation of the immune infiltrating level with the expression of PIMREG in LUAD. Results PIMREG was highly expressed in patients with LUAD. Specifically, the level of PIMREG gradually increased from pathological stage I to IV. Further, we validated the higher expression of PIMREG expressed in LUAD cell lines. Moreover, PIMREG had a high diagnostic value, with an -AUC of 0.955. Kaplan–Meier survival and Cox regression analyses revealed that the high expression of PIMREG was independently associated with poor clinical outcomes. In our prognostic nomogram, the expression of PIMREG implied a significant prognostic value. Gene set enrichment analysis (GSEA) identified that the high expression PIMREG phenotype was involved in the mitotic cell cycle, mRNA splicing, DNA repair, Rho GTPase signaling, TP53 transcriptional regulation, and translation pathways. Next, we also explored the correlation of PIMREG and tumor-immune interactions and found a negative correlation between PIMREG and the immune infiltrating level of T cells, macrophages, B cells, dendritic cells (DCs) , and CD8+ T cells in LUAD. Conclusions High levels of PIMREG correlated with poor prognosis and immune infiltrates in LUAD.

Download Full-text

Identification of Seven Novel Ferroptosis-Related Long Non - Coding RNA Signatures as a Diagnostic Biomarker for Acute Myeloid Leukemia

10.21203/rs.3.rs-502041/v1 ◽

2021 ◽

Author(s):

Zhiyuan Zheng ◽

Wei Wu ◽

Zehang Lin ◽

Shuhan Liu ◽

Qiaoqian Chen ◽

...

Keyword(s):

Acute Myeloid Leukemia ◽

Risk Score ◽

Myeloid Leukemia ◽

Cox Regression ◽

Enrichment Analysis ◽

Gene Set Enrichment Analysis ◽

Gene Set Enrichment ◽

Gene Set ◽

Kaplan Meier ◽

Acute Myeloid

Abstract Background: Ferroptosis is a newly discovered type of programmed cell death that participates in the biological processes of various cancers. However, the mechanism by which ferroptosis modulates acute myeloid leukemia (AML) remains unclear. This study aimed to investigate the role of ferroptosis-related long non-coding RNAs (lncRNAs) in AML and establish a corresponding prognostic model.Methods: RNA-sequencing data and clinicopathological characteristics were obtained from The Cancer Genome Atlas database, and ferroptosis-related genes were obtained from the FerrDb database. The “limma” R package, Cox regression, and the least absolute shrinkage and selection operator were used to determine the ferroptosis-related lncRNA signature with the lowest Akaike information criteria (AIC). The risk score of ferroptosis-related lncRNAs was calculated and patients with AML were divided into high- and low-risk groups based on the median risk score. The Kaplan-Meier curve and Cox regression were used to evaluate the prognostic value of the risk score. Finally, gene set enrichment analysis (GSEA) and single-sample gene set enrichment analysis (ssGSEA) were performed to explore the biological functions of the ferroptosis-related lncRNAs.Results: Seven ferroptosis-related lncRNA signatures were identified in the training group, and Kaplan-Meier and Cox regression analyses confirmed that risk scores were independent prognostic predictors of AML in both the training and validation groups (All P < 0.05). In addition, the area under the curve (AUC) analysis confirmed that the signatures had a good predictive ability for the prognosis of AML. GSEA and ssGSEA showed that the seven ferroptosis-related lncRNAs were related to glutathione metabolism and tumor immunity.Conclusions: In this study, seven novel ferroptosis-related lncRNA signatures (AP001266.2, AC133961.1, AF064858.3, AC007383.2, AC008906.1, AC026771.1, and KIF26B-AS1) were established. These signatures were shown to accurately predict the prognosis of AML, which would provide new insights into strategies for the development of new AML therapies.

Download Full-text

The Prognostic Value of m6A-related LncRNAs in Patients with HNSCC

10.21203/rs.3.rs-778502/v1 ◽

2021 ◽

Author(s):

Liu-qing Zhou ◽

Jie-yu Zhou ◽

Yao Hu

Keyword(s):

Prognostic Value ◽

Risk Model ◽

Cox Regression ◽

Predictive Ability ◽

Enrichment Analysis ◽

The Cancer Genome Atlas ◽

Consensus Clustering ◽

Cox Regression Analysis ◽

Kaplan Meier ◽

Cancer Genome Atlas

Abstract Background: N6-methyladenosine (m6A) modifications play an essential role in tumorigenesis. m6A modifications are known to modulate RNAs, including mRNAs and lncRNAs. However, the prognostic role of m6A-related lncRNAs in head and neck squamous cell carcinoma (HNSCC) is poorly understood.Methods: Based on LASSO Cox regression, enrichment analysis, univariate and multivariate Cox regression analysis, a risk prognostic model, and consensus clustering analysis, we analyzed the 12 m6A-related lncRNAs in HNSCC samples data using the data from The Cancer Genome Atlas (TCGA) database.Results: We found twelve m6A-related lncRNAs in the training cohort and validated in all cohorts by Kaplan-Meier and Cox regression analyses, and revealing their independent prognostic value in HNSCC. Moreover, ROC analysis was conducted, confirming the strong predictive ability of this signature for HNSCC prognosis. GSEA and detailed immune infiltration analyses revealed specific pathways associated with m6A-related lncRNAs.Conclusions: In this study, a novel risk model including twelve genes (SAP30L-AS1, AC022098.1, LINC01475, AC090587.2, AC008115.3, AC015911.3, AL122035.2, AC010226.1, AL513190.1, ZNF32-AS1, AL035587.1 and AL031716.1) was built. It could accurately predict HNSCC prognosis and provide potential prediction outcome and new therapeutic target for HNSCC patients.

Download Full-text

PVT1 is a prognostic marker associated with immune invasion of bladder urothelial carcinoma

Mathematical Biosciences and Engineering ◽

10.3934/mbe.2022009 ◽

2021 ◽

Vol 19 (1) ◽

pp. 169-190

Author(s):

Peiyuan Li ◽

◽

Gangjie Qiao ◽

Jian Lu ◽

Wenbin Ji ◽

...

Keyword(s):

Urothelial Carcinoma ◽

Immune Cell ◽

Cox Regression ◽

Enrichment Analysis ◽

Gene Set Enrichment Analysis ◽

The Cancer Genome Atlas ◽

Gene Set Enrichment ◽

Cox Regression Analysis ◽

Gene Set ◽

Bladder Urothelial Carcinoma

<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 < 0.05) and overall survival (OS P < 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>

Download Full-text

CXCL5 Has Potential to Be a Marker for Hepatocellular Carcinoma Prognosis and Was Correlating With Immune Infiltrates

Frontiers in Oncology ◽

10.3389/fonc.2021.637023 ◽

2021 ◽

Vol 11 ◽

Author(s):

Yuan Nie ◽

Mei-chun Jiang ◽

Cong Liu ◽

Qi Liu ◽

Xuan Zhu

Keyword(s):

Hepatocellular Carcinoma ◽

Cox Regression ◽

Enrichment Analysis ◽

Gene Set Enrichment Analysis ◽

The Cancer Genome Atlas ◽

Hospital Data ◽

Gene Set Enrichment ◽

Local Hospital ◽

High Expression Group ◽

Cancer Genome Atlas

BackgroundsTumor microenvironment (TME) plays a crucial role in the initiation and progression of Hepatocellular Carcinoma (HCC), especially immune infiltrates. However, there is still a challenge in understanding the modulation of the immune and stromal components in TME, especially TME related genes.MethodsThe proportion of tumor-infiltrating immune cells (TICs) and the immune and stromal scores in 374 HCC patients from The Cancer Genome Atlas (TCGA) database were determined using CIBERSORT and ESTIMATE computational methods. The final screened genes were confirmed by the PPI network and univariate Cox regression of the differentially expressed genes based on different immune or stromal scores. The correlation between the expression levels of the final gene interactions and the clinical characteristics was based on TCGA database and local hospital data. Gene set enrichment analysis (GSEA) and the effect of CXCL5 expression on TICs were conducted.ResultsThere were correlations between the expression of CXCL5 and survival of HCC patients and TMN classification both in TCGA database and local hospital data. The immune-related activities were enriched in the high-expression group; however, the metabolic pathways were enriched in the low-expression group. The result of CIBERSORT analyzing had indicated that CXCL5 expression were correlated with the proportion of NK cells activated, macrophages M0, Mast cells resting, Neutrophils.ConclusionsCXCL5 was a potential prognostic marker for HCC and provides clues regarding immune infiltrates, which offers extra insight for therapeutics of HCC, however, more independent cohorts and functional experiments of CXCL5 are warranted.

Download Full-text

The prognostic value of m6A-related LncRNAs in patients with HNSCC: bioinformatics analysis of TCGA database

Scientific Reports ◽

10.1038/s41598-021-04591-z ◽

2022 ◽

Vol 12 (1) ◽

Author(s):

Liu-qing Zhou ◽

Jin-xiong Shen ◽

Jie-yu Zhou ◽

Yao Hu ◽

Hong-jun Xiao

Keyword(s):

Prognostic Value ◽

Risk Model ◽

Cox Regression ◽

Predictive Ability ◽

Enrichment Analysis ◽

The Cancer Genome Atlas ◽

Consensus Clustering ◽

Cox Regression Analysis ◽

Kaplan Meier ◽

Cancer Genome Atlas

AbstractN6-methyladenosine (m6A) modifications play an essential role in tumorigenesis. These modifications modulate RNAs, including mRNAs and lncRNAs. However, the prognostic role of m6A-related lncRNAs in head and neck squamous cell carcinoma (HNSCC) is poorly understood. Based on LASSO Cox regression, enrichment analysis, univariate and multivariate Cox regression analysis, a prognostic risk model, and consensus clustering analysis, we analyzed 12 m6A-related lncRNAs in HNSCC sample data from The Cancer Genome Atlas (TCGA) database. We found 12 m6A-related lncRNAs in the training cohort and validated them in all cohorts by Kaplan–Meier and Cox regression analyses, revealing their independent prognostic value in HNSCC. Moreover, ROC analysis was conducted, confirming the strong predictive ability of this signature for HNSCC survival. GSEA and detailed immune infiltration analyses revealed specific pathways associated with m6A-related lncRNAs. In this study, a novel risk model including twelve genes (SAP30L-AS1, AC022098.1, LINC01475, AC090587.2, AC008115.3, AC015911.3, AL122035.2, AC010226.1, AL513190.1, ZNF32-AS1, AL035587.1 and AL031716.1) was built. It could accurately predict HNSCC outcomes and could provide new therapeutic targets for HNSCC patients.

Download Full-text

Identification of seven novel ferroptosis-related long non-coding RNA signatures as a diagnostic biomarker for acute myeloid leukemia

BMC Medical Genomics ◽

10.1186/s12920-021-01085-9 ◽

2021 ◽

Vol 14 (1) ◽

Author(s):

Zhiyuan Zheng ◽

Wei Wu ◽

Zehang Lin ◽

Shuhan Liu ◽

Qiaoqian Chen ◽

...

Keyword(s):

Acute Myeloid Leukemia ◽

Risk Score ◽

Myeloid Leukemia ◽

Cox Regression ◽

Enrichment Analysis ◽

Gene Set Enrichment Analysis ◽

Gene Set Enrichment ◽

Gene Set ◽

Kaplan Meier ◽

Acute Myeloid

Abstract Background Ferroptosis is a newly discovered type of programmed cell death that participates in the biological processes of various cancers. However, the mechanism by which ferroptosis modulates acute myeloid leukemia (AML) remains unclear. This study aimed to investigate the role of ferroptosis-related long non-coding RNAs (lncRNAs) in AML and establish a corresponding prognostic model. Methods RNA-sequencing data and clinicopathological characteristics were obtained from The Cancer Genome Atlas database, and ferroptosis-related genes were obtained from the FerrDb database. The “limma” R package, Cox regression, and the least absolute shrinkage and selection operator were used to determine the ferroptosis-related lncRNA signature with the lowest Akaike information criteria (AIC). The risk score of ferroptosis-related lncRNAs was calculated and patients with AML were divided into high- and low-risk groups based on the median risk score. The Kaplan–Meier curve and Cox regression were used to evaluate the prognostic value of the risk score. Finally, gene set enrichment analysis (GSEA) and single-sample gene set enrichment analysis (ssGSEA) were performed to explore the biological functions of the ferroptosis-related lncRNAs. Results Seven ferroptosis-related lncRNA signatures were identified in the training group, and Kaplan–Meier and Cox regression analyses confirmed that risk scores were independent prognostic predictors of AML in both the training and validation groups (All P < 0.05). In addition, the area under the curve (AUC) analysis confirmed that the signatures had a good predictive ability for the prognosis of AML. GSEA and ssGSEA showed that the seven ferroptosis-related lncRNAs were related to glutathione metabolism and tumor immunity. Conclusions In this study, seven novel ferroptosis-related lncRNA signatures (AP001266.2, AC133961.1, AF064858.3, AC007383.2, AC008906.1, AC026771.1, and KIF26B-AS1) were established. These signatures were shown to accurately predict the prognosis of AML, which would provide new insights into strategies for the development of new AML therapies.

Download Full-text

Molecular Features of Glioma Determined and Validated Using Combined TCGA and GTEx Data Analyses

Frontiers in Oncology ◽

10.3389/fonc.2021.729137 ◽

2021 ◽

Vol 11 ◽

Author(s):

Zijiang Yang ◽

Weiyi Gong ◽

Ting Zhang ◽

Heng Gao

Keyword(s):

Expression Analysis ◽

Molecular Mechanisms ◽

Enrichment Analysis ◽

Gene Set Enrichment Analysis ◽

The Cancer Genome Atlas ◽

Normal Brain ◽

Gene Set Enrichment ◽

The Core ◽

Kaplan Meier ◽

Cancer Genome Atlas

Gliomas are among the most common intracranial tumors which originated from neuroepithelial cells. Increasing evidence has revealed that long noncoding RNA (lncRNA)-microRNA (miRNA)-mRNA module regulation and tumor-infiltrating immune cells play important regulatory roles in the occurrence and progression of gliomas. However, the precise underlying molecular mechanisms remain largely unknown. Data on gliomas in The Cancer Genome Atlas lack normal control samples; to overcome this limitation, we combined 665 The Cancer Genome Atlas glioma RNA sequence datasets with 188 Genotype-Tissue Expression normal brain RNA sequences to construct an expression matrix profile after normalization. We systematically analyzed the expression of mRNAs, lncRNAs, and miRNAs between gliomas and normal brain tissues. Kaplan–Meier survival analyses were conducted to screen differentially expressed mRNAs, lncRNAs, and miRNAs. A prognostic miRNA-related competitive endogenous RNA network was constructed, and the core subnetworks were filtered using 6 miRNAs, 3 lncRNAs, and 11 mRNAs. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses were performed to investigate the biological functions of significantly dysregulated mRNAs. Co-expression network analysis was performed to analyze and screen the core genes. Furthermore, single-sample Gene Set Enrichment Analysis and immune checkpoint gene expression analysis were performed, as co-expression analysis indicated immune gene dysregulation in glioma. Finally, the expression of representative dysregulated genes was validated in U87 cells at the transcriptional level, establishing a foundation for further research. We identified 7017 mRNAs, 437 lncRNAs, and 9 miRNAs that were differentially expressed in gliomas. Kaplan–Meier survival analysis revealed 5684 mRNAs, 61 lncRNAs, and 7 miRNAs with potential as prognostic signatures in patients with glioma. The hub subnetwork of the competing endogenous RNA network between PART1-hsa-mir-25-SLC12A5/TACC2/BSN/TLN2/ZDHHC8 was screened out. Gene co-expression network, single-sample Gene Set Enrichment Analysis, and immune checkpoint expression analysis demonstrated that tumor-infiltrating immune cells are closely related to gliomas. We identified novel potential biomarkers to predict survival and therapeutic targets for patients with gliomas based on a large-scale sample. Importantly, we filtered pivotal genes that provide valuable information for further exploration of the molecular mechanisms underlying glioma tumorigenesis and progression.

Download Full-text

Development of a Novel Prognostic Signature Based on Antigen Processing and Presentation in Patients with Breast Cancer

Pathology & Oncology Research ◽

10.3389/pore.2021.600727 ◽

2021 ◽

Vol 27 ◽

Author(s):

Aoshuang Qi ◽

Mingyi Ju ◽

Yinfeng Liu ◽

Jia Bi ◽

Qian Wei ◽

...

Keyword(s):

Breast Cancer ◽

Antigen Processing ◽

Cox Regression ◽

Enrichment Analysis ◽

Gene Signature ◽

Gene Set Enrichment Analysis ◽

The Cancer Genome Atlas ◽

Prognostic Indicators ◽

Gene Set ◽

Antigen Processing And Presentation

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) < 1), and HLA-F was risky (HR > 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.

Download Full-text

DPYSL2 as potential diagnostic and prognostic biomarker linked to immune infiltration in lung adenocarcinoma

World Journal of Surgical Oncology ◽

10.1186/s12957-021-02379-z ◽

2021 ◽

Vol 19 (1) ◽

Author(s):

Yang-Jie Wu ◽

Ai-Tao Nai ◽

Gui-Cheng He ◽

Fei Xiao ◽

Zhi-Min Li ◽

...

Keyword(s):

Lung Adenocarcinoma ◽

Roc Curve ◽

Enrichment Analysis ◽

Gene Set Enrichment Analysis ◽

List Type ◽

Gene Set Enrichment ◽

Gene Set ◽

Normal Tissues ◽

Kaplan Meier ◽

Kaplan Meier Plotter

Abstract Background Dihydropyrimidinase like 2 (DPYSL2) has been linked to tumor metastasis. However, the function of DPSY2L in lung adenocarcinoma (LUAD) is yet to be explored. Methods Herein, we assessed DPYSL2 expression in various tumor types via online databases such as Oncomine and Tumor Immune Estimation Resource (TIMER). Further, we verified the low protein and mRNA expressions of DPYSL2 in LUAD via the ULCAN, The TCGA and GEPIA databases. We applied the ROC curve to examine the role of DPYSL2 in diagnosis. The prognostic significance of DPYSL2 was established through the Kaplan–Meier plotter and the Cox analyses (univariate and multivariate). TIMER was used to explore DPYSL2 expression and its connection to immune infiltrated cells. Through Gene Set Enrichment Analysis, the possible mechanism of DPYSL2 in LUAD was investigated. Results In this study, database analysis revealed lower DPYSL2 expression in LUAD than in normal tissues. The ROC curve suggested that expression of DPYSL2 had high diagnostic efficiency in LUAD. The DPYSL2 expression had an association with the survival time of LUAD patients in the Kaplan–Meier plotter and the Cox analyses. The results from TIMER depicted a markedly positive correlation of DPYSL2 expression with immune cells infiltrated in LUAD, such as macrophages, dendritic cells, CD4+ T cells, and neutrophils. Additionally, many gene markers for the immune system had similar positive correlations in the TIMER analysis. In Gene Set Enrichment Analysis, six immune-related signaling pathways were associated with DPYSL2. Conclusions In summary, DPYSL2 is a novel biomarker with diagnostic and prognostic potential for LUAD as well as an immunotherapy target. Highlights Expression of DPYSL2 was considerably lower in LUAD than in normal tissues. Investigation of multiple databases showed a high diagnostic value of DPYSL2 in LUAD. DPYSL2 can independently predict the LUAD outcomes. Immune-related mechanisms may be potential ways for DPYSL2 to play a role in LUAD.

Download Full-text

Identification of Epithelial–Mesenchymal Transition-Related Prognostic lncRNAs Biomarkers Associated With Melanoma Microenvironment

Frontiers in Cell and Developmental Biology ◽

10.3389/fcell.2021.679133 ◽

2021 ◽

Vol 9 ◽

Author(s):

Bo Xiao ◽

Liyan Liu ◽

Zhuoyuan Chen ◽

Aoyu Li ◽

Pingxiao Wang ◽

...

Keyword(s):

Cox Regression ◽

Epithelial Mesenchymal Transition ◽

Enrichment Analysis ◽

Roc Curves ◽

Risk Groups ◽

Gene Set Enrichment Analysis ◽

Gene Set Enrichment ◽

Mesenchymal Transition ◽

Cox Regression Analysis ◽

Kaplan Meier

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.

Download Full-text