Identification of m6A methyltransferase-related lncRNA signature for predicting immunotherapy and prognosis in patients with hepatocellular carcinoma

N6-methyladenosine (m6A) methyltransferase has been shown to be an oncogene in a variety of cancers. Nevertheless, the relationship between the long non-coding RNAs (lncRNAs) and hepatocellular carcinoma (HCC) remains elusive. We integrated the gene expression data of 371 HCC and 50 normal tissues from 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 were performed to identify m6A methyltransferase‑related lncRNAs. Three prognostic lncRNAs were selected by univariate and LASSO Cox regression analyses to construct the m6A methyltransferase-related lncRNA signature. Multivariate Cox regression analyses illustrated that this signature was an independent prognostic factor for overall survival (OS) prediction. The Gene Set Enrichment Analysis (GSEA) suggested that the m6A methyltransferase-related lncRNAs were involved in the immune-related biological processes and pathways. Besides, we discovered that the lncRNAs signature was correlated with the tumor microenvironment (TME) and the expression of critical immune checkpoints. Tumor Immune Dysfunction and Exclusion (TIDE) analysis revealed that the lncRNAs could predict the clinical response to immunotherapy. Our study had originated a prognostic signature for HCC based on the potential prognostic m6A methyltransferase-related lncRNAs. This study had deepened the understanding of the TME status of HCC patients and laid a theoretical foundation for the choice of immunotherapy.

Download Full-text

Significance of m 6A Methyltransferase‑Related lncRNAs in the Prognosis and Immunotherapy of Hepatocellular Carcinoma

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

2021 ◽

Author(s):

Lili Li ◽

Rongrong Xie ◽

Qichun Wei

Keyword(s):

Hepatocellular Carcinoma ◽

Correlation Analysis ◽

Immune Cells ◽

Cox Regression ◽

Enrichment Analysis ◽

Gene Set Enrichment Analysis ◽

The Cancer Genome Atlas ◽

Immune Checkpoints ◽

Regression Analyses ◽

Cox Regression Analysis

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.

Download Full-text

Decreased DHRS1 expression is a novel predictor of poor survival in patients with hepatocellular carcinoma

Biomarkers in Medicine ◽

10.2217/bmm-2021-0041 ◽

2021 ◽

Vol 15 (15) ◽

pp. 1319-1331

Author(s):

Li Li ◽

Hui-Jing Situ ◽

Wen-Cheng Ma ◽

Xuan Liu ◽

Lu-Lu Wang

Keyword(s):

Hepatocellular Carcinoma ◽

Cox Regression ◽

Enrichment Analysis ◽

Gene Set Enrichment Analysis ◽

Aberrant Expression ◽

Normal Tissues ◽

Multiple Data ◽

Kaplan Meier ◽

Survival Gene ◽

Mrna And Protein Expression

Aim: To investigate the effect of aberrant expression of DHRS1 on hepatocellular carcinoma (HCC). Materials & methods: Kaplan–Meier and Cox regression analyses were performed to evaluate the correlation between DHRS1 and overall survival. Gene set enrichment analysis was performed to explore the potential function of DHRS1 in HCC. Results: Multiple data analysis revealed that DHRS1 mRNA and protein expression level were remarkably lower in HCC than that in normal tissues. In survival analysis, patients with low DHRS1 expression presented a poorer prognosis, and was an independent risk factor for HCC. Conclusion: Decreased DHRS1 expression may be a potential predictor of poor prognosis in HCC.

Download Full-text

Identification of the HMMR Gene as a Diagnostic and Prognostic Biomarker in Hepatocellular Carcinoma Based on Integrated Bioinformatics Analysis

Evidence-based Complementary and Alternative Medicine ◽

10.1155/2021/5970085 ◽

2021 ◽

Vol 2021 ◽

pp. 1-8

Author(s):

Honglan Guo ◽

Qinqiao Fan

Keyword(s):

Hepatocellular Carcinoma ◽

Roc Curve ◽

Cox Regression ◽

Enrichment Analysis ◽

Cancer Genome ◽

Gene Set Enrichment Analysis ◽

The Cancer Genome Atlas ◽

Roc Curve Analysis ◽

Cox Regression Analysis ◽

Normal Tissues

Background. We aimed to investigate the expression of the hyaluronan-mediated motility receptor (HMMR) gene in hepatocellular carcinoma (HCC) and nonneoplastic tissues and to investigate the diagnostic and prognostic value of HMMR. Method. With the reuse of the publicly available The Cancer Genome Atlas (TCGA) data, 374 HCC patients and 50 nonneoplastic tissues were used to investigate the diagnostic and prognostic values of HMMR genes by receiver operating characteristic (ROC) curve analysis and survival analysis. All patients were divided into low- and high-expression groups based on the median value of HMMR expression level. Univariate and multivariate Cox regression analysis were used to identify prognostic factors. Gene set enrichment analysis (GSEA) was performed to explore the potential mechanism of the HMMR genes involved in HCC. The diagnostic and prognostic values were further validated in an external cohort from the International Cancer Genome Consortium (ICGC). Results. HMMR mRNA expression was significantly elevated in HCC tissues compared with that in normal tissues from both TCGA and the ICGC cohorts (all P values <0.001). Increased HMMR expression was significantly associated with histologic grade, pathological stage, and survival status (all P values <0.05). The area under the ROC curve for HMMR expression in HCC and normal tissues was 0.969 (95% CI: 0.948–0.983) in the TCGA cohort and 0.956 (95% CI: 0.932–0.973) in the ICGC cohort. Patients with high HMMR expression had a poor prognosis than patients with low expression group in both cohorts (all P < 0.001 ). Univariate and multivariate analysis also showed that HMMR is an independent predictor factor associated with overall survival in both cohorts (all P values <0.001). GSEA showed that genes upregulated in the high-HMMR HCC subgroup were mainly significantly enriched in the cell cycle pathway, pathways in cancer, and P53 signaling pathway. Conclusion. HMMR is expressed at high levels in HCC. HMMR overexpression may be an unfavorable prognostic factor for HCC.

Download Full-text

DNM1: A Prognostic Biomarker Associated with Immune Infiltration in Colon Cancer—A Study Based on TCGA Database

BioMed Research International ◽

10.1155/2021/4896106 ◽

2021 ◽

Vol 2021 ◽

pp. 1-9

Author(s):

Mingchao Hu ◽

Jianchun Gu ◽

Wenzhao Su ◽

Zhenjie Zhang ◽

Baosong Zhu ◽

...

Keyword(s):

Colon Cancer ◽

Immune Cell ◽

Cox Regression ◽

Univariate Analysis ◽

Enrichment Analysis ◽

Gene Set Enrichment Analysis ◽

The Cancer Genome Atlas ◽

Receptor Interaction ◽

Regression Analyses ◽

Normal Tissues

Aim. The aim of our work was to determine the utility of DNM1 as a biomarker for the diagnosis and prognosis of colon cancer (CC). Methods. DNM1 expression variations in CC vs. normal tissues were investigated using The Cancer Genome Atlas (TCGA) database. The association of DNM1 expression levels with the clinicopathological variables in CC prognosis was investigated using logistic regression analyses. Independent prognostic factors for CC were evaluated using univariate and multivariate Cox regression analyses. The correlation between DNM1 expression and immune cell infiltration was estimated using single-sample Gene Set Enrichment Analysis (ssGSEA). Results. DNM1 expression in CC tissues was significantly higher than that in normal tissues. High DNM1 expression was significantly correlated with M stage, N stage, perineural invasion and lymphatic invasion and predicted poor prognosis. The univariate analysis highlighted that DNM1 was an independent CC risk factor. Results of ssGSEA showed that DNM1 was linked to several cancer-related pathways, including the neuroactive ligand-receptor interaction, hypertrophic cardiomyopathy, ECM-receptor interaction, dilated cardiomyopathy, and calcium signaling pathway. Moreover, DNM1 expression was positively correlated with the level of infiltration by Neutrophils, Tregs, NK cells, and Macrophages. Conclusion. DNM1 has a significant function and has diagnostic and prognostic potential for CC.

Download Full-text

Hypoxic Characteristic in the Immunosuppressive Microenvironment of Hepatocellular Carcinoma

Frontiers in Immunology ◽

10.3389/fimmu.2021.611058 ◽

2021 ◽

Vol 12 ◽

Author(s):

Zhuomao Mo ◽

Daiyuan Liu ◽

Dade Rong ◽

Shijun Zhang

Keyword(s):

Hepatocellular Carcinoma ◽

Immune Cell ◽

Cox Regression ◽

Cancer Genome ◽

Gene Set Enrichment Analysis ◽

The Cancer Genome Atlas ◽

Cell Infiltration ◽

Immune Checkpoints ◽

Immune Cell Infiltration ◽

Immunosuppressive Microenvironment

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.

Download Full-text

Identification and Validation of an Immune-Related eRNA Prognostic Signature for Hepatocellular Carcinoma

Frontiers in Genetics ◽

10.3389/fgene.2021.657051 ◽

2021 ◽

Vol 12 ◽

Author(s):

Shenglan Cai ◽

Xingwang Hu ◽

Ruochan Chen ◽

Yiya Zhang

Keyword(s):

Hepatocellular Carcinoma ◽

Immune Cells ◽

Cox Regression ◽

Expression Profiles ◽

Gene Expression Profiles ◽

Roc Curves ◽

The Cancer Genome Atlas ◽

Immune Checkpoints ◽

Immune Genes ◽

Normal Tissues

BackgroundEnhancer RNAs (eRNAs) are intergenic long non-coding RNAs (lncRNAs) that participate in the progression of malignancies by targeting tumor-related genes and immune checkpoints. However, the potential role of eRNAs in hepatocellular carcinoma (HCC) is unclear. In this study, we aimed to construct an immune-related eRNA prognostic model that could be used to prospectively assess the prognosis of patients with HCC.MethodsGene expression profiles of patients with HCC were downloaded from The Cancer Genome Atlas (TCGA). The eRNAs co-expressed from immune genes were identified as immune-related eRNAs. Cox regression analyses were applied in a training cohort to construct an immune-related eRNA signature (IReRS), that was subsequently used to analyze a testing cohort and combination of the two cohorts. Kaplan-Meier and receiver operating characteristic (ROC) curves were used to validate the predictive effect in the three cohorts. Gene Set Enrishment Analysis (GSEA) computation was used to identify an IReRS-related signaling pathway. A web-based cell type identification by estimating relative subsets of RNA transcripts (CIBERSORT) computation was used to evaluate the relationship between the IReRS and infiltrating immune cells.ResultsA total of sixty-four immune-related eRNAs (IReRNAs) was identified in HCC, and 14 IReRNAs were associated with overall survival (OS). Five IReRNAs were used for constructing an immune-related eRNA signature (IReRS), which was shown to correlate with poor survival and to be an independent prognostic biomarker for HCC. The GSEA results showed that the IReRS was correlated to cancer-related and immune-related pathways. Moreover, we found that IReRS was correlated to infiltrating immune cells, including CD8+ T cells and M0 macrophages. Finally, differential expressions of the five risk IReRNAs in tumor tissues vs. adjacent normal tissues and their prognostic values were verified, in which the AL445524.1 may function as an oncogene that affects prognosis partly by regulating CD4-CLTA4 related genes.ConclusionOur results suggest that the IReRS could serve as a biomarker for predicting prognosis in patients with HCC. Additionally, it may be correlated to the tumor immune microenvironment and could also be used as a biomarker in immunotherapy for HCC.

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

A Glycolysis-Based Long Non-coding RNA Signature Accurately Predicts Prognosis in Renal Carcinoma Patients

Frontiers in Genetics ◽

10.3389/fgene.2021.638980 ◽

2021 ◽

Vol 12 ◽

Author(s):

Honghao Cao ◽

Hang Tong ◽

Junlong Zhu ◽

Chenchen Xie ◽

Zijia Qin ◽

...

Keyword(s):

Clinical Significance ◽

Risk Score ◽

Cox Regression ◽

Risk Group ◽

Enrichment Analysis ◽

Roc Curves ◽

Gene Set Enrichment Analysis ◽

The Cancer Genome Atlas ◽

Survival Prediction ◽

Kaplan Meier

BackgroundThe prognosis of renal cell carcinoma (RCC) varies greatly among different risk groups, and the traditional indicators have limited effect in the identification of risk grade in patients with RCC. The purpose of our study is to explore a glycolysis-based long non-coding RNAs (lncRNAs) signature and verify its potential clinical significance in prognostic prediction of RCC patients.MethodsIn this study, RNA data and clinical information were downloaded from The Cancer Genome Atlas (TCGA) database. Univariate and multivariate cox regression displayed six significantly related lncRNAs (AC124854.1, AC078778.1, EMX2OS, DLGAP1-AS2, AC084876.1, and AC026401.3) which were utilized in construction of risk score by a formula. The accuracy of risk score was verified by a series of statistical methods such as receiver operating characteristic (ROC) curves, nomogram and Kaplan-Meier curves. Its potential clinical significance was excavated by gene enrichment analysis.ResultsKaplan-Meier curves and ROC curves showed reliability of the risk score to predict the prognosis of RCC patients. Stratification analysis indicated that the risk score was independent predictor compare to other traditional clinical parameters. The clinical nomogram showed highly rigorous with index of 0.73 and precisely predicted 1-, 3-, and 5-year survival time of RCC patients. Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene set enrichment analysis (GSEA) depicted the top ten correlated pathways in both high-risk group and low-risk group. There are 6 lncRNAs and 25 related mRNAs including 36 lncRNA-mRNA links in lncRNA-mRNA co-expression network.ConclusionThis research demonstrated that glycolysis-based lncRNAs possessed an important value in survival prediction of RCC patients, which would be a potential target for future treatment.

Download Full-text

Increased Expression Together with the Gene Regulation Network of NEIL3 Predicts Poor Prognosis and Correlates with Immune Infiltrates in Hepatocellular Carcinoma

10.21203/rs.2.23392/v1 ◽

2020 ◽

Author(s):

Bo Hu ◽

Xiao-Bo Yang ◽

Xinting Sang

Keyword(s):

Hepatocellular Carcinoma ◽

Gene Regulation ◽

Poor Prognosis ◽

Cox Regression ◽

Enrichment Analysis ◽

Gene Set Enrichment Analysis ◽

The Cancer Genome Atlas ◽

Immune Gene ◽

Data Set ◽

Gene Markers

Abstract Background: Abnormal Nei endonuclease VIII-like 3 (NEIL3)expression is associated with carcinogenesis. Methods: We used sequencing data from the Cancer Genome Atlas database, analyzed NEIL3 expression, gene regulation networks and the correlation with immune infiltrates in hepatocellular carcinoma (HCC). Clinicopathologic characteristics associated with overall survival in TCGA patients using Cox regression and the Kaplan-Meier method. Gene Set Enrichment Analysis was performed using TCGA data set. LinkedOmics was used to identify differential gene expression with NEIL3 and to analyze Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathways. Gene enrichment analysis examined target networks of kinases and transcription factors.Correlations between NEIL3 expression and cancer immune infiltrates and immune gene markers were analyzed by TIMER and GEPIA. Results: We found that overexpressed NEIL3 predicted poor prognosis. Functional network analysis suggested that NEIL3 regulates the DNA replication and cell cycle signaling via pathways involving several cancer-related kinases and E2F Transcription Factor 1.NEIL3 was also found to be associated with the infiltration of several immune cells. Conclusions: Our results demonstrate that data mining efficiently reveals information about NEIL3 expression, potential regulatory networks and the relationship with immune infiltration in HCC, laying a foundation for further study of the role of NEIL3 in carcinogenesis.

Download Full-text

A hypoxia-related signature for clinically predicting diagnosis, prognosis and immune microenvironment of hepatocellular carcinoma patients.

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

2020 ◽

Author(s):

Baohui Zhang ◽

Bufu Tang ◽

Jianyao Gao ◽

Jiatong Li ◽

Lingming Kong ◽

...

Keyword(s):

Hepatocellular Carcinoma ◽

Immune Cell ◽

Cox Regression ◽

Enrichment Analysis ◽

High Risk Group ◽

Cancer Genome ◽

Gene Set Enrichment Analysis ◽

The Cancer Genome Atlas ◽

Diagnostic Model ◽

Potential Biomarker

Abstract Background Hypoxia plays an indispensable role in the development of hepatocellular carcinoma (HCC). However, there are few studies on the application of hypoxia molecules in the prognosis predicting of HCC. We aimed to identify the hypoxia-related genes in HCC and construct reliable models for diagnosis, prognosis and recurrence of HCC patients as well as exploring the potential mechanism.Methods Differentially expressed genes (DEGs) analysis was performed using The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database and four clusters were determined by a consistent clustering analysis. Three DEGs closely related to overall survival（OS）were identified using Cox regression and LASSO analysis and the hypoxia-related signature was developed and validated in TCGA and International Cancer Genome Consortium (ICGC) database. Then the Gene Set Enrichment Analysis (GSEA) was performed to explore signaling pathways regulated by the signature and the CIBERSORT was used for estimating the fractions of immune cell types.Results A total of 397 hypoxia-related DEGs were detected and three genes (PDSS1, CDCA8 and SLC7A11) were selected to construct a prognosis, recurrence and diagnosis model. Then patients were divided into high- and low-risk groups. Our hypoxia-related signature was significantly associated with worse prognosis and higher recurrence rate. The diagnostic model also accurately distinguished HCC from normal samples and nodules. Furthermore, the hypoxia-related signature could positively regulate immune response and the high-risk group had higher fractions of macrophages, B memory cells and follicle-helper T cells, and exhibited higher expression of immunocheckpoints such as PD1and PDL1.Conclusions Altogether, our study showed that hypoxia-related signature is a potential biomarker for diagnosis, prognosis and recurrence of HCC, and it provided an immunological perspective for developing personalized therapies.

Download Full-text