scholarly journals Immune-related miRNA signature identifies prognosis and immune landscape in head and neck squamous cell carcinomas

2020 ◽  
Vol 40 (11) ◽  
Author(s):  
Bo Ma ◽  
Hui Li ◽  
Jia Qiao ◽  
Tao Meng ◽  
Riyue Yu

Abstract Background: Head and neck squamous cell carcinoma (HNSCC) is recognised as an immune active cancer, but little is known about the role of microRNAs (miRNAs) in it. In the present study, we aim to determine a prognostic and immune-related miRNAs signature (IRMS) in HNSCC. Methods: Spearman correlation analysis was used to screen out prognostic immune-related miRNAs based on single-sample gene set enrichment analysis (ssGSEA). Least absolute shrinkage and selection operator (LASSO) Cox regression model was used to establish IRMS in HNSCC. Then, the influence of the IRMS on HNSCC was comprehensively analysed. Results: We obtained 11 prognostic immune-related miRNAs based on ssGSEA. Then an IRMS integrated with six miRNAs was established through LASSO Cox regression analysis. The stratification survival analysis indicated that IRMS was independent from other characteristics and performed favourably in the overall survival (OS) prediction. The function annotation suggested that IRMS was highly associated with the immune-related response biological processes and pathways which are so important for tumorigenesis of HNSCC. Moreover, the nomogram demonstrated that our model was identified as an independent prognostic factor. In addition, we found that IRMS was significantly correlated with the immune infiltration and expression of critical immune checkpoints, indicating that the poor prognosis might be caused partly by immunosuppressive microenvironment. Conclusion: We established a novel IRMS, which exhibited a potent prognostic value and could be representative of immune status in HNSCC.

2022 ◽  
Author(s):  
Feng Liu ◽  
Zewei Tu ◽  
Junzhe Liu ◽  
Xiaoyan Long ◽  
Bing Xiao ◽  
...  

Background: The role of DNAJC10 in cancers have been reported but its function in glioma is not clear. We reveal the prognostic role and underlying functions of DNAJC10 in glioma in this study. Methods: Reverse Transcription and Quantitative Polymerase Chain Reaction (RT-qPCR) was used to quantify the relative DNAJC10 mRNA expression of clinical samples. Protein expressions of clinical samples were tested by Western blot. The overall survival (OS) of glioma patients with different DNAJC10 expression was compared by Kaplan-Meier method (two-sided log-rank test). Single-sample gene set enrichment analysis (ssGSEA) was used to estimate the immune cell infiltrations and immune-related function levels. The independent prognostic role of DNAJC10 was determined by univariate and multivariate Cox regression analysis. The DNAJC10-based nomogram model was established using multivariate Cox regression by R package “rms”. Results: Higher DNAJC10 is observed in gliomas and it’s upregulated in higher grade, IDH-wild, 1p/19q non-codeletion, MGMT unmethylated gliomas. Gliomas with higher DNAJC10 expression present poorer prognosis compared with low-DNAJC10 gliomas. The predictive accuracy of 1/3/5-OS of DNAJC10 is found stable and robust using time-dependent ROC model. Enrichment analysis recognized that T-cell activation and T-cell receptor signaling were enriched in higher DNAJC10 gliomas. Immune/stromal cell infiltrations, tumor mutation burden (TMB), copy Number Alteration (CNA) burden, and immune check-point genes were also positively correlated with DNAJC10 expression in gliomas. DNAJ10-based nomogram model was established and showed strong prognosis-predictive ability. Conclusion: Higher DNAJC10 expression correlates with poor prognosis of glioma and it was a potential prognostic biomarker for glioma.


2021 ◽  
Author(s):  
Pengxiang Li ◽  
Dongchun Qin ◽  
Xuefeng Lv ◽  
Lu Liu ◽  
Mengle Peng

Abstract Background: Cervical cancer (CC) is the most common reproductive neoplasm in women, especially in developing countries. Ferroptosis, a novel type of cell death, and lncRNAs play critical roles in the prognosis of CC patients and antitumor immunity. Methods: A ferroptosis-related lncRNA signature (FRLS) was constructed by LASSO Cox regression analysis. Kaplan-Meier analysis, receiver operating characteristic (ROC) curve analysis, multivariate analysis, and nomogram were used to evaluate and predict the FRLS. Based on the FRLS, immune-related genes, the tumor microenvironment (TME), immune checkpoints, and immunotherapy were investigated. Results: The FRLS was composed of ten lncRNAs and was markedly associated with the overall survival (OS) of CC patients. Gene set enrichment analysis (GSEA) demonstrated that the FRLS was largely associated with immune-related pathways. Weighted gene co-expression network analysis (WGCNA) was performed to analyze immune-related genes and to identify the optimal modules and genes. TLR4 was eventually identified, and its expression was verified in the Gene Expression Omnibus (GEO) database. Then, quantitative real-time PCR (qRT-PCR) was used to validate the results in CC and paracancerous tissues. Besides, our results showed that CD8+ T cells were significantly correlated between the low- and high-risk groups, and it could modulate ferroptosis during tumor immunotherapy. The expression of immune checkpoints was substantially different between the two groups. Additionally, tumor immune dysfunction and exclusion (TIDE) was applied to predict the sensitivity of immune checkpoint inhibitor (ICI) treatment. Conclusion: The FRLS established was significantly associated with prognosis; moreover, the FRLS is a prospective therapeutic target, and combined with immunotherapy, can be used in the treatment of CC patients.


2021 ◽  
Author(s):  
Lili Li ◽  
Rongrong Xie ◽  
Qichun Wei

Abstract Background: Hepatocellular carcinoma (HCC) is one of the leading causes of mortality worldwide. N6-methyladenosine (m6A) methyltransferase, has been proved to act as an oncogene in several human cancers. However, little is known about its relationship with the long non-coding RNAs (lncRNAs) that remains elusive in HCC.Methods: We comprehensively integrated gene expression data acquired from 371 HCC and 50 normal tissues in The Cancer Genome Atlas (TCGA) database. Differentially expressed protein-coding genes (DE-PCGs)/lncRNAs (DE-lncRs) analysis and univariate regression & Kaplan-Meier (K-M) analysis was performed to identify m6A methyltransferase‑related lncRNAs that were related to overall survival (OS). m6A methyltransferase‑related lncRNA signature was constructed using the Least absolute shrinkage and selection operator (LASSO) Cox regression analyses. Furthermore, Cox regression analysis was applied to identify independent prognostic factors in HCC. The signature was validated in an internal validation set. Finally, the correlation analysis between gene signature and immune cells infiltration was also investigated via single-sample Gene Set Enrichment Analysis (ssGSEA) and immunotherapy response was calculated through Tumor Immune Dysfunction and Exclusion (TIDE) algorithm.Results: A total of 21 m6A methyltransferase-related lncRNAs were screened out according to Spearman correlation analysis with the immune score (|R| > 0.3, P < 0.05). We selected 3 prognostic lncRNAs to construct m6A methyltransferase-related lncRNA signature through univariate and LASSO Cox regression analyses. The univariate and multivariate Cox regression analyses demonstrated that the lncRNAs signature was a robust independent prognostic factor in OS prediction with high accuracy. The GSEA also suggested that the m6A methyltransferase-related lncRNAs were involved in the immune-related biological processes and pathways which were very well-known in the context of HCC tumorigenesis. Besides, we found that the lncRNAs signature was strikingly correlated with the tumor microenvironment (TME) immune cells infiltration and expression of critical immune checkpoints. Finally, results from the TIDE analysis revealed that the m6A methyltransferase-related lncRNAs could efficiently predict the clinical response of immunotherapy in HCC.Conclusion: Together, our study screened potential prognostic m6A methyltransferase related lncRNAs and established a novel m6A methyltransferase-based prognostic model of HCC, which not only provides new potential prognostic biomarkers and therapeutic targets but also deepens our understanding of tumor immune microenvironment status and lays a theoretical foundation for immunotherapy.


2021 ◽  
Author(s):  
Hongyang Liu ◽  
Junhu Wan ◽  
Quanling Feng ◽  
Jingyu Li ◽  
Jun Liu ◽  
...  

Abstract Background: Endometrial cancer (EC) is one of the most common types of gynecological cancer. Hypoxia is an important clinical feature and regulates various tumor processes. However, the prognostic value of hypoxia-related lncRNA in EC remains to be further elucidated. Here, we aimed to characterize the molecular features of EC by the development of a classification system based on the expression profile of hypoxia-related lncRNA.Methods: Univariate Cox regression analysis was used to identify hypoxia-related lncRNAs associated with overall survival. The least absolute shrinkage and selection operator (LASSO) Cox regression model was utilized to construct gene signature. Multivariate Cox regression analysis and receiver operating characteristic (ROC) curve analysis were also performed. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEEG) pathway, and Gene Set Enrichment Analysis (GESA) were used to identify hypoxia-related lncRNA pathway. Western blot and real-time PCR were used to detect target gene expression. The cell proliferation was determined by using WST-1 assay.Results: Based on univariate Cox regression analysis, we identified 17 hypoxia-related lncRNAs significantly associated with overall survival. Next, the least absolute shrinkage and selection operator (LASSO) Cox regression model was utilized to construct a multigene signature in the TCGA EC cohort. The risk score was confirmed as an independent predictor for overall survival in multivariate Cox regression analysis and receiver operating characteristic (ROC) curve analysis. Besides, the survival time of EC patients in different risk group was significantly correlated to clinicopathologic factors, such as age, stage and grade. Furthermore, hypoxia-related lncRNA associated with the high-risk group were involved in various aspects of the malignant progression of EC via Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEEG) pathway, and Gene Set Enrichment Analysis (GESA). Besides, using CIBERSORT analysis, we found a different immune cell environment characterization of EC between different cluster and risk group. Moreover, the risk score was closely correlated to immunotherapy response, microsatellite instability and tumor mutation burden (TMB). Finally, we select one hypoxia-related lncRNA SOS1-IT1 to validate its role in hypoxia and EC progression. Interestingly, we found SOS1-IT1 was overexpressed in tumor tissues, and closely correlated with clinicopathological parameters of EC. The expression level of SOS1-IT1 was significantly increased under hypoxia condition. Additionally, the important hypoxia regulatory factor HIF-1α can directly bind SOS1-IT1 promoter region, and affect its expression level. Conclusions: In summary, this study established a new EC classification based on the hypoxia-related lncRNA signature, thereby provide a novel sight to understand the potential mechanism of human EC development.


BMC Cancer ◽  
2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Guichuan Huang ◽  
Jing Zhang ◽  
Ling Gong ◽  
Yi Huang ◽  
Daishun Liu

Abstract Background Lung cancer is one of the most lethal and most prevalent malignant tumors worldwide, and lung squamous cell carcinoma (LUSC) is one of the major histological subtypes. Although numerous biomarkers have been found to be associated with prognosis in LUSC, the prediction effect of a single gene biomarker is insufficient, especially for glycolysis-related genes. Therefore, we aimed to develop a novel glycolysis-related gene signature to predict survival in patients with LUSC. Methods The mRNA expression files and LUSC clinical information were obtained from The Cancer Genome Atlas (TCGA) dataset. Results Based on Gene Set Enrichment Analysis (GSEA), we found 5 glycolysis-related gene sets that were significantly enriched in LUSC tissues. Univariate and multivariate Cox proportional regression models were performed to choose prognostic-related gene signatures. Based on a Cox proportional regression model, a risk score for a three-gene signature (HKDC1, ALDH7A1, and MDH1) was established to divide patients into high-risk and low-risk subgroups. Multivariate Cox regression analysis indicated that the risk score for this three-gene signature can be used as an independent prognostic indicator in LUSC. Additionally, based on the cBioPortal database, the rate of genomic alterations in the HKDC1, ALDH7A1, and MDH1 genes were 1.9, 1.1, and 5% in LUSC patients, respectively. Conclusion A glycolysis-based three-gene signature could serve as a novel biomarker in predicting the prognosis of patients with LUSC and it also provides additional gene targets that can be used to cure LUSC patients.


2020 ◽  
Vol 40 (8) ◽  
Author(s):  
Sihan Chen ◽  
Guodong Cao ◽  
Wei Wu ◽  
Yida Lu ◽  
Xiaobo He ◽  
...  

Abstract Colon adenocarcinoma (COAD) is a malignant gastrointestinal tumor, often occurring in the left colon, which is regulated by glycolysis-related processes. In past studies, multiple genes that influence the prognosis for survival have been discovered through bioinformatics analysis. However, the prediction of disease prognosis using a single gene is not an accurate method. In the present study, a mechanistic model was established to achieve better prediction for the prognosis of COAD. COAD-related data downloaded from The Cancer Genome Atlas (TCGA) were correlated with the glycolysis process using gene set enrichment analysis (GSEA) to determine the glycolysis-related genes that regulate COAD. Using COX regression analysis, glycolysis-related genes associated with the prognosis of COAD were identified, and the genes screened to establish a predictive model. The risk scores of this model were correlated with relevant clinical data to obtain a connection diagram between the model and survival rate, tumor characteristic data, etc. Finally, genes in the model were correlated with cells in the tumor microenvironment, finding that they affected specific immune cells in the model. Seven genes related to glycolysis were identified (PPARGC1A, DLAT, 6PC2, P4HA1, STC2, ANKZF1, and GPC1), which affect the prognosis of patients with COAD and constitute the model for prediction of survival of COAD patients.


PeerJ ◽  
2021 ◽  
Vol 9 ◽  
pp. e11273
Author(s):  
Lei Yang ◽  
Weilong Yin ◽  
Xuechen Liu ◽  
Fangcun Li ◽  
Li Ma ◽  
...  

Background Hepatocellular carcinoma (HCC) is considered to be a malignant tumor with a high incidence and a high mortality. Accurate prognostic models are urgently needed. The present study was aimed at screening the critical genes for prognosis of HCC. Methods The GSE25097, GSE14520, GSE36376 and GSE76427 datasets were obtained from Gene Expression Omnibus (GEO). We used GEO2R to screen differentially expressed genes (DEGs). A protein-protein interaction network of the DEGs was constructed by Cytoscape in order to find hub genes by module analysis. The Metascape was performed to discover biological functions and pathway enrichment of DEGs. MCODE components were calculated to construct a module complex of DEGs. Then, gene set enrichment analysis (GSEA) was used for gene enrichment analysis. ONCOMINE was employed to assess the mRNA expression levels of key genes in HCC, and the survival analysis was conducted using the array from The Cancer Genome Atlas (TCGA) of HCC. Then, the LASSO Cox regression model was performed to establish and identify the prognostic gene signature. We validated the prognostic value of the gene signature in the TCGA cohort. Results We screened out 10 hub genes which were all up-regulated in HCC tissue. They mainly enrich in mitotic cell cycle process. The GSEA results showed that these data sets had good enrichment score and significance in the cell cycle pathway. Each candidate gene may be an indicator of prognostic factors in the development of HCC. However, hub genes expression was weekly associated with overall survival in HCC patients. LASSO Cox regression analysis validated a five-gene signature (including CDC20, CCNB2, NCAPG, ASPM and NUSAP1). These results suggest that five-gene signature model may provide clues for clinical prognostic biomarker of HCC.


2021 ◽  
Vol 12 ◽  
Author(s):  
Zhuomao Mo ◽  
Daiyuan Liu ◽  
Dade Rong ◽  
Shijun Zhang

Background: Generally, hepatocellular carcinoma (HCC) exists in an immunosuppressive microenvironment that promotes tumor evasion. Hypoxia can impact intercellular crosstalk in the tumor microenvironment. This study aimed to explore and elucidate the underlying relationship between hypoxia and immunotherapy in patients with HCC.Methods: HCC genomic and clinicopathological datasets were obtained from The Cancer Genome Atlas (TCGA-LIHC), Gene Expression Omnibus databases (GSE14520) and International Cancer Genome Consortium (ICGC-LIRI). The TCGA-LIHC cases were divided into clusters based on single sample gene set enrichment analysis and hierarchical clustering. After identifying patients with immunosuppressive microenvironment with different hypoxic conditions, correlations between immunological characteristics and hypoxia clusters were investigated. Subsequently, a hypoxia-associated score was established by differential expression, univariable Cox regression, and lasso regression analyses. The score was verified by survival and receiver operating characteristic curve analyses. The GSE14520 cohort was used to validate the findings of immune cell infiltration and immune checkpoints expression, while the ICGC-LIRI cohort was employed to verify the hypoxia-associated score.Results: We identified hypoxic patients with immunosuppressive HCC. This cluster exhibited higher immune cell infiltration and immune checkpoint expression in the TCGA cohort, while similar significant differences were observed in the GEO cohort. The hypoxia-associated score was composed of five genes (ephrin A3, dihydropyrimidinase like 4, solute carrier family 2 member 5, stanniocalcin 2, and lysyl oxidase). In both two cohorts, survival analysis revealed significant differences between the high-risk and low-risk groups. In addition, compared to other clinical parameters, the established score had the highest predictive performance at both 3 and 5 years in two cohorts.Conclusion: This study provides further evidence of the link between hypoxic signals in patients and immunosuppression in HCC. Defining hypoxia-associated HCC subtypes may help reveal potential regulatory mechanisms between hypoxia and the immunosuppressive microenvironment, and our hypoxia-associated score could exhibit potential implications for future predictive models.


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.


Sign in / Sign up

Export Citation Format

Share Document