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

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):  
Lili Li ◽  
Rongrong Xie ◽  
Guangrong Lu

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.


2021 ◽  
Vol 2021 ◽  
pp. 1-8
Author(s):  
Honglan Guo ◽  
Qinqiao Fan

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.


2020 ◽  
Author(s):  
Ze-bing Song ◽  
Guo-pei Zhang ◽  
shaoqiang li

Abstract Background: Hepatocellular carcinoma (HCC) is one of the most common malignant tumor in the world which prognosis is poor. Therefore, a precise biomarker is needed to guide treatment and improve prognosis. More and more studies have shown that lncRNAs and immune response are closely related to the prognosis of hepatocellular carcinoma. The aim of this study was to establish a prognostic signature based on immune related lncRNAs for HCC.Methods: Univariate cox regression analysis was performed to identify immune related lncRNAs, which had negative correlation with overall survival (OS) of 370 HCC patients from The Cancer Genome Atlas (TCGA). A prognostic signature based on OS related lncRNAs was identified by using multivariate cox regression analysis. Gene set enrichment analysis (GSEA) and a competing endogenous RNA (ceRNA) network were performed to clarify the potential mechanism of lncRNAs included in prognostic signature. Results: A prognostic signature based on OS related lncRNAs (AC145207.5, AL365203.2, AC009779.2, ZFPM2-AS1, PCAT6, LINC00942) showed moderately in prognosis prediction, and related with pathologic stage (Stage I&II VS Stage III&IV), distant metastasis status (M0 VS M1) and tumor stage (T1-2 VS T3-4). CeRNA network constructed 15 aixs among differentially expressed immune related genes, lncRNAs included in prognostic signature and differentially expressed miRNA. GSEA indicated that these lncRNAs were involved in cancer-related pathways. Conclusion: We constructed a prognostic signature based on immune related lncRNAs which can predict prognosis and guide therapies for HCC.


2020 ◽  
Author(s):  
Andi Ma ◽  
Yukai Sun ◽  
Racheal O. Ogbodu ◽  
Ling Xiao ◽  
Haibing Deng ◽  
...  

Abstract Background: It is well known that long non-coding RNAs (lncRNAs) play a vital role in cancer. We aimed to explore the prognostic value of potential immune-related lncRNAs in hepatocellular carcinoma (HCC). Methods: Validated the established lncRNA signature of 343 patients with HCC from The Cancer Genome Atlas (TCGA) and 81 samples from Gene Expression Omnibus (GEO). Immune-related lncRNAs for HCC prognosis were evaluated using Cox regression and Least Absolute Shrinkage and Selection Operator (LASSO) analyses. LASSO analysis was performed to calculate a risk score formula to explore the difference in overall survival between high- and low-risk groups in TCGA, which was verified using GEO, Gene Ontology (GO), and pathway-enrichment analysis. These analyses were used to identify the function of screened genes and construct a co-expression network of these genes. Results: Using computational difference algorithms and lasso Cox regression analysis, the differentially expressed and survival-related immune-related genes (IRGs) among patients with HCC were established as five novel immune-related lncRNA signatures (AC099850.3, AL031985.3, PRRT3-AS1, AC023157.3, MSC-AS1). Patients in the low‐risk group showed significantly better survival than patients in the high‐risk group ( P = 3.033e−05). The signature identified can be an effective prognostic factor to predict patient survival. The nomogram showed some clinical net benefits predicted by overall survival. In order to explore its underlying mechanism, several methods of enrichment were elucidated using Gene Set Enrichment Analysis. Conclusion: Identifying five immune-related lncRNA signatures has important clinical implications for predicting patient outcome and guiding tailored therapy for patients with HCC with further prospective validation.


2020 ◽  
Vol 10 (1) ◽  
Author(s):  
Qingmiao Shi ◽  
Chen Xue ◽  
Xin Yuan ◽  
Yuting He ◽  
Zujiang Yu

Abstract Hepatocellular carcinoma (HCC) ranks fourth in cancer-related mortality worldwide. N1-methyladenosine (m1A), a methylation modification on RNA, is gaining attention for its role across diverse biological processes. However, m1A-related regulatory genes expression, its relationship with clinical prognosis, and its role in HCC remain unclear. In this study, we utilized The Cancer Genome Atlas-Liver Hepatocellular Carcinoma (TCGA-LIHC) database to investigate alterations within 10 m1A-related regulatory genes and observed a high mutation frequency (23/363). Cox regression analysis and least absolute shrinkage and selection operator were used to explore the association between m1A-related regulatory genes expression and HCC patient survival and identified four regulators that were remarkably associated with HCC patient prognosis. Additionally, an independent cohort from International Cancer Genome Consortium was studied to validate our discoveries and found to be consistent with those in the TCGA dataset. In terms of mechanism, gene set enrichment analysis linked these four genes with various physiological roles in cell division, the MYC pathway, protein metabolism, and mitosis. Kyoto Encyclopedia of Genes and Genomes analysis revealed that PI3K/Akt signaling pathway had potential relevance to m1A-related regulatory genes in HCC. These findings indicate that m1A-related regulatory genes may play crucial roles in regulating HCC progression and be exploited for diagnostic and prognostic purposes.


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):  
Shenglan Cai ◽  
Xingwang Hu ◽  
Ruochan Chen ◽  
Yiya Zhang

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.


2021 ◽  
Vol 7 ◽  
Author(s):  
Xiaoyu Deng ◽  
Qinghua Bi ◽  
Shihan Chen ◽  
Xianhua Chen ◽  
Shuhui Li ◽  
...  

Although great progresses have been made in the diagnosis and treatment of hepatocellular carcinoma (HCC), its prognostic marker remains controversial. In this current study, weighted correlation network analysis and Cox regression analysis showed significant prognostic value of five autophagy-related long non-coding RNAs (AR-lncRNAs) (including TMCC1-AS1, PLBD1-AS1, MKLN1-AS, LINC01063, and CYTOR) for HCC patients from data in The Cancer Genome Atlas. By using them, we constructed a five-AR-lncRNA prognostic signature, which accurately distinguished the high- and low-risk groups of HCC patients. All of the five AR lncRNAs were highly expressed in the high-risk group of HCC patients. This five-AR-lncRNA prognostic signature showed good area under the curve (AUC) value (AUC = 0.751) for the overall survival (OS) prediction in either all HCC patients or HCC patients stratified according to several clinical traits. A prognostic nomogram with this five-AR-lncRNA signature predicted the 3- and 5-year OS outcomes of HCC patients intuitively and accurately (concordance index = 0.745). By parallel comparison, this five-AR-lncRNA signature has better prognosis accuracy than the other three recently published signatures. Furthermore, we discovered the prediction ability of the signature on therapeutic outcomes of HCC patients, including chemotherapy and immunotherapeutic responses. Gene set enrichment analysis and gene mutation analysis revealed that dysregulated cell cycle pathway, purine metabolism, and TP53 mutation may play an important role in determining the OS outcomes of HCC patients in the high-risk group. Collectively, our study suggests a new five-AR-lncRNA prognostic signature for HCC patients.


Sign in / Sign up

Export Citation Format

Share Document