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 and 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 (BPs) 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. The present study had deepened the understanding of the TME status of HCC patients and laid a theoretical foundation for the choice of immunotherapy.     

Immune-related long noncoding RNA signature for predicting survival and immune checkpoint blockade in hepatocellular carcinoma

Long noncoding RNAs (lncRNAs) show multiple functions, including immune response. Recently, the immune-related lncRNAs have been reported in some cancers. We first investigated the immune-related lncRNA signature as a potential target in hepatocellular carcinoma (HCC) survival. The training set (n = 368) and the independent external validation cohort (n = 115) were used. Immune genes and lncRNAs coexpression were constructed to identify immune-related lncRNAs. Cox regression analyses were perfumed to establish the immune-related lncRNA signature. Regulatory roles of this signature on cancer pathways and the immunologic features were investigated. The correlation between immune checkpoint inhibitors and this signature was examined. In this study, the immune-related lncRNA signature was identified in HCC, which could stratify patients into high- and low-risk groups. This immune-related lncRNA signature was correlated with disease progression and worse survival and was an independent prognostic biomarker. Our immune-related lncRNA signature was still a powerful tool in predicting survival in each stratum of age, gender, and tumor stage. This signature mediated cell cycle, glycolysis, DNA repair, mammalian target of rapamycin signaling, and immunologic characteristics (i.e., natural killer cells vs. Th1 cells down, etc). This signature was associated with immune cell infiltration (i.e., macrophages M0, Tregs, CD4 memory T cells, and macrophages M1, etc.,) and immune checkpoint blockade (ICB) immunotherapy-related molecules (i.e., PD-L1, PD-L2, and IDO1). Our findings suggested that the immune-related lncRNA signature had an important value for survival prediction and may have the potential to measure the response to ICB immunotherapy. This signature may guide the selection of the immunotherapy for HCC.     

A dual immune signature of CD8+ T cells and MMP9 improves the survival of patients with hepatocellular carcinoma

The 5-year survival of hepatocellular carcinoma (HCC) is difficult due to the high recurrence rate and metastasis. Tumor infiltrating immune cells (TICs) and immune-related genes (IRGs) bring hope to improve survival and treatment of HCC patients. However, there are problems in predicting immune signatures and identifying novel therapeutic targets. In the study, the CIBERSORT algorithm was used to evaluate 22 immune cell infiltration patterns in gene expression omnibus (GEO) and the cancer genome atlas (TCGA) data. Eight immune cells were found to have significant infiltration differences between the tumor and normal groups. The CD8+ T cells immune signature was constructed by least absolute shrinkage and selection operator (LASSO) algorithm. The high infiltration level of CD8+ T cells could significantly improve survival of patients. The weighted gene co-expression network analysis (WGCNA) algorithm identified MMP9 was closely related to the overall survival of HCC patients. K-M survival and tROC analysis confirmed that MMP9 had an excellent prognostic prediction. Cox regression showed that a dual immune signature of CD8+ T cells and MMP9 was independent survival factor in HCC. Therefore, a dual prognostic immune signature could improve the survival of patient and may provide a new strategy for the immunotherapy of HCC.     

Construction and validation of a TP53-associated immune prognostic model for gastric cancer

Increasing evidence indicates that TP53 mutation impacts the patients' prognosis by regulating the gastric cancer (GC) immunophenotype. An immune prognostic signature (IPS) was constructed based on TP53 status. The effects of the IPS on the immune microenvironment of GC were analyzed. We also constructed a nomogram integrating the IPS and other clinical factors. An IPS was constructed in the TCGA cohort and validated in the meta-GEO cohort. TP53 mutation resulted in the downregulation of the immune response in GC. Concretely, high-risk patients were characterized by increased monocyte, macrophage M0 and T cell follicular helper infiltration; increased stromal score, ESTIMATE score and immune score; higher TIM3 and BTLA expression; and decreased dendritic cell and T cell CD4 memory-activated infiltration and tumor purity. The nomogram also showed good predictive performance. These results suggest that the IPS is an effective prognostic indicator for GC patients, which might provide a theoretical foundation for immunotherapy.     

Identification of a prognostic signature of epithelial ovarian cancer based on tumor immune microenvironment exploration

This study aims to develop an immune-related genes (IRGs) prognostic signature to stratify the epithelial ovarian cancer (EOC) patients. We identified 332 up- and 154 down-regulated EOC-specific IRGs. As a result, candidate IRGs were idendified to construct prognostic models respectivy for overall survial and progression-free survival. The risk score was validated as a risk factor for prognosis and was used to built a combined nomogram. According to the IRG-related prognostic model, EOC patients were divided into high- and low- risk group and were further explored their association with tumor immune microenvironment (TME). CIBERSORT algorithm showed higher macrophages M1 cell, T cells follicular helper cell and plasma cells infiltrating levels in the low-risk group. In addition, the low-risk group was found with higher immunophenoscore and distinct mutation signatures compared with the high-risk group. These findings may shed light on the development of novel immune biomarkers and target therapy of EOC.     

Prediction of overall survival and response to immune checkpoint inhibitors: An immune-related signature for gastric cancer

Gastric cancer (GC) is common in East Asia and South and Central America. Most GC patients miss the opportunities for surgery. Despite their therapeutic potential, immune checkpoint inhibitors (ICIs) only work in part of patients with GC. Thus, this study was aimed at constructing a signature for diagnosis, prognosis, and prediction of response to ICIs. A multivariate analysis showed that the 8-immune-related-gene (IRG) signature was an independent prognostic factor of overall survival among GC patients. In the high-risk group of 8IRG signature risk score, the fractions of CD4 T cells, macrophage M2 and monocyte, which is associated with the progression of cancers, were higher. The low-risk group had a higher immunophenoscore, which meant a better response to ICIs.     

Identification of a four-gene metabolic signature predicting overall survival for hepatocellular carcinoma

While hundreds of consistently altered metabolic genes had been identified in hepatocellular carcinoma (HCC), the prognostic role of them remains to be further elucidated. Messenger RNA expression profiles and clinicopathological data were downloaded from The Cancer Genome Atlas—Liver Hepatocellular Carcinoma and GSE14520 data set from the Gene Expression Omnibus database. Univariate Cox regression analysis and lasso Cox regression model established a novel four-gene metabolic signature (including acetyl-CoA acetyltransferase 1, glutamic-oxaloacetic transaminase 2, phosphatidylserine synthase 2, and uridine-cytidine kinase 2) for HCC prognosis prediction. Patients in the high-risk group shown significantly poorer survival than patients in the low-risk group. The signature was significantly correlated with other negative prognostic factors such as higher α-fetoprotein. The signature was found to be an independent prognostic factor for HCC survival. Nomogram including the signature shown some clinical net benefit for overall survival prediction. Furthermore, gene set enrichment analyses revealed several significantly enriched pathways, which might help explain the underlying mechanisms. Our study identified a novel robust four-gene metabolic signature for HCC prognosis prediction. The signature might reflect the dysregulated metabolic microenvironment and provided potential biomarkers for metabolic therapy and treatment response prediction in HCC.     

Profiles of immune infiltration in lung adenocarcinoma and their clinical significant: A gene-expression–based retrospective study

Lung cancer is one of the fatal tumors. The tumor microenvironment plays a key role in regulating tumor progression. To figure out the role of tumor microenvironment in lung adenocarcinoma (LUAD), ESTIMATE algorithm was used to evaluate the immune scores in LUAD. Patients with low immune scores had a worse overall survival (OS) compared with high immune scores. Using RNA-Seq data of 489 patients in The Cancer Genome Atlas (TCGA), differentially expressed genes (DEGs) were identified between high- and low-immune score groups. Based on the DEGs, nine-gene signature was constructed by the least absolute shrinkage and selection operator Cox regression model in TCGA set. The signature demonstrated significant prognostic value in both TCGA and Gene Expression Omnibus database. Multivariate Cox regression analyses indicated that nine-genes signature was an independent prognostic factor. Subgroup analysis also revealed a robust prognostic ability of nine-gene signature. A nomogram with a C-index of 0.722 had a favorable power for predicting 3-, 5-, and 10-year survival for clinical use by integrating nine-gene signature and other clinical features. Co-expression and functional enrichment analysis showed that nine-gene signature was significantly associated with immune response and provided potential profound molecules for revealing the mechanism of tumor initiation and progression. In conclusion, we revealed the significance of immune infiltration and built a novel nine-gene signature as a reliable prognostic factor for patients with LUAD.     

基于肿瘤免疫微环境鉴定IDH突变型弥漫性胶质瘤的预后标志物

异柠檬酸脱氢酶(IDH)突变存在于大多数低级别胶质细胞瘤中,免疫逃逸是肿瘤标志性特征之一,免疫治疗在胶质瘤的治疗中的作用越来越重要。利用生信手段分析TCGA、CGGA、GEO数据集中IDH突变胶质瘤的770个免疫相关基因及其临床相关数据,从而获得每个患者的免疫风险评分(IMRS);结合IMRS和临床信息,筛选出6个差异表达基因(TRAF3、ATG10、BID、TAB1、MAP3K1、RPS6)组成IMRS模型并生成诺莫图对患者预后进行评估,发现低风险组患者的总生存期(OS)较高风险组均明显延长。此外,签名相关免疫细胞浸润分析发现肿瘤相关巨噬细胞浸润评分(TAM)与肿瘤相关T细胞浸润评分(TIS)呈明显的负相关,表明高IMRS富集了促肿瘤免疫浸润,而低IMRS则富集了相对较多的抗肿瘤免疫浸润。     

Typical tumor immune microenvironment status determine prognosis in lung adenocarcinoma

BackgroundImmune cells, vital components of tumor microenvironment, regulate tumor survival and progression. Lung adenocarcinoma (LUAD), the tumor with the highest mortality rate worldwide, reconstitutes tumor immune microenvironment (TIME) to avoid immune destruction. Data have shown that TIME influences LUAD prognosis and predicts immunotherapeutic efficacy. The related information about the role of TIME's characteristics in LUAD is limited.MethodsWe performed unsupervised consensus clustering via machine-learning techniques to identify TIME clusters among 1906 patients and gathered survival data. The characteristics of TIME clusters of LUAD were visualized by multi-omics analysis, pseudo-time dynamic analysis, and enrichment analysis. TIME score model was constructed by principal component analysis. Comprehensive analysis and validation were conducted to test the prognostic efficacy and immunotherapeutic response of TIME score.ResultsTIME clusters (A, B and C) were constructed and exhibited different immune infiltration states. Multi-omics analyses included significant mutated genes (SMG), copy number variation (CNV) and cancer stemness that were significantly different among the three clusters. TIME cluster A had a lower SMG, lower CNV, and lower stemness but a higher immune infiltration level compared to TIME clusters B and C. TIME score showed that patients in low TIME score group had higher overall survival rates, higher immune infiltration level and high expression of immune checkpoints. In validation cohorts, low TIME score subgroup had better drug sensitivity and favorable immunotherapeutic response.ConclusionWe constructed a stable model of LUAD immune microenvironment characteristics that may improve the prognostic accuracy of patients, provide improved explanations of LUAD responses to immunotherapy, and provide new strategies for LUAD treatment.     

Assessment of clinical studies evaluating combinations of immune checkpoint inhibitors with locoregional treatments in solid tumors

In the last decade, immunotherapy with immune checkpoint inhibitors (ICIs) has changed the therapeutic algorithm of cancer patients. ICIs combined with other therapeutic options, such as chemo- and targeted therapies, generate impressive results in cancer patients. Locoregional treatments (LRTs) play an important role in the management of various solid tumors (e.g., hepatocellular carcinoma (HCC), neuroendocrine tumors, etc.), and this therapeutic approach may enhance the activity of the immune response to tumor cells destroying primary tumors and leading to the release of several soluble molecules. This systematic review was performed to identify studies reporting objective response rate (ORR) and survival information in patients with solid tumors treated with ICIs plus LRTs. In the present work, fourteen studies were included, and the majority of them (five studies) enrolled patients with hepatocellular carcinoma (HCC), whereas the others included patients with different diseases. The highest ORRs were seen in HCC (67%, Y-90 RE plus ipilimumab and nivolumab) and melanoma (38%, dendritic cells with mRNA plus ipilimumab) patients. ORRs were not observed in liver metastases from melanoma and colorectal cancer. These data suggest that combination of ICIs and LRTs is feasible and more active in primary tumors (particularly HCC) than metastases with a synergistic effect on antitumor immunity. However, further studies are needed to better select patients, schedules, and setting of treatments.     

Development and validation of a 14-gene signature for prognosis prediction in hepatocellular carcinoma

Worldwide, hepatocellular carcinoma (HCC) remains a crucial medical problem. Precise and concise prognostic models are urgently needed because of the intricate gene variations among liver cancer cells. We conducted this study to identify a prognostic gene signature with biological significance. We applied two algorithms to generate differentially expressed genes (DEGs) between HCC and normal specimens in The Cancer Genome Atlas cohort (training set included) and performed enrichment analyses to expound on their biological significance. A protein-protein interactions network was established based on the STRING online tool. We then used Cytoscape to screen hub genes in crucial modules. A multigene signature was constructed by Cox regression analysis of hub genes to stratify the prognoses of HCC patients in the training set. The prognostic value of the multigene signature was externally validated in two other sets from Gene Expression Omnibus (GSE14520 and GSE76427), and its role in recurrence prediction was also investigated. A total of 2000 DEGs were obtained, including 1542 upregulated genes and 458 downregulated genes. Subsequently, we constructed a 14-gene signature on the basis of 56 hub genes, which was a good predictor of overall survival. The prognostic signature could be replicated in GSE14520 and GSE76427. Moreover, the 14-gene signature could be applied for recurrence prediction in the training set and GSE14520. In summary, the 14-gene signature extracted from hub genes was involved in some of the HCC-related signalling pathways; it not only served as a predictive signature for HCC outcome but could also be used to predict HCC recurrence.     

Profiles of immune-related genes and immune cell infiltration in the tumor microenvironment of diffuse lower-grade gliomas

The tumor microenvironment is highly correlated with tumor occurrence, progress, and prognosis. We aimed to investigate the immune-related gene (IRG) expression and immune infiltration pattern in the tumor microenvironment of lower-grade glioma (LGG). We employed the Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data (ESTIMATE) algorithm to calculate immune and stromal scores and identify prognostic IRG based on The Cancer Genome Atlas data set. The potential molecular functions of these genes were explored with the help of functional enrichment analysis and the protein–protein interaction network. Remarkably, three cohorts that were downloaded from the Chinese Glioma Genome Atlas database were analyzed to further verify the prognostic values of these genes. Moreover, the Tumor IMmune Estimation Resource (TIMER) algorithm was used to estimate the abundance of infiltrating immune cells and explore the immune infiltration pattern in LGG. And unsupervised cluster analysis determined three clusters of the immune infiltration pattern and indicated that CD8⁺ T cells and macrophages were significantly associated with LGG outcomes. Altogether, our study identified a list of prognostic IRGs and provided a perspective to explore the immune infiltration pattern in LGG.     

A novel tailored immune gene pairs signature for overall survival prediction in lower-grade gliomas

Lower-grade gliomas (LGGs) have a good prognosis with a wide range of overall survival (OS) outcomes. An accurate prognostic system can better predict survival time. An RNA-Sequencing (RNA-seq) prognostic signature showed a better predictive power than clinical predictor models. A signature constructed using gene pairs can transcend changes from biological heterogeneity, technical biases, and different measurement platforms. RNA-seq coupled with corresponding clinical information were extracted from The Cancer Genome Atlas (TCGA) and the Chinese Glioma Genome Atlas (CGGA). Immune-related gene pairs (IRGPs) were used to establish a prognostic signature through univariate and multivariate Cox proportional hazards regression. Weighted gene co-expression network analysis (WGCNA) was used to evaluate module eigengenes correlating with immune cell infiltration and to construct gene co-expression networks. Samples in the training and testing cohorts were dichotomized into high- and low-risk groups. Risk score was identified as an independent predictor, and exhibited a closed relationship with prognosis. WGCNA presented a gene set that was positively correlated with age, WHO grade, isocitrate dehydrogenase (IDH) mutation status, 1p/19 codeletion, risk score, and immune cell infiltrations (CD4 T cells, B cells, dendritic cells, and macrophages). A nomogram comprising of age, WHO grade, 1p/19q codeletion, and three gene pairs (BIRC5|SSTR2, BMP2|TNFRSF12A, and NRG3|TGFB2) was established as a tool for predicting OS. The IPGPs signature, which is associated with immune cell infiltration, is a novel tailored tool for individual-level prediction.     

Comprehensive analysis of prognostic immune-related genes in the tumor microenvironment of cutaneous melanoma

Cutaneous malignant melanoma (hereafter called melanoma) is one of the most aggressive cancers with increasing incidence and mortality rates worldwide. In this study, we performed a systematic investigation of the tumor microenvironmental and genetic factors associated with melanoma to identify prognostic biomarkers for melanoma. We calculated the immune and stromal scores of melanoma patients from the Cancer Genome Atlas (TCGA) using the ESTIMATE algorithm and found that they were closely associated with patients’ prognosis. Then the differentially expressed genes were obtained based on the immune and stromal scores, and prognostic immune-related genes further identified. Functional analysis and the protein–protein interaction network further revealed that these genes enriched in many immune-related biological processes. In addition, the abundance of six infiltrating immune cells was analyzed using prognostic immune-related genes by TIMER algorithm. The unsupervised clustering analysis using immune-cell proportions revealed eight clusters with distinct survival patterns, suggesting that dendritic cells were most abundant in the microenvironment and CD8⁺ T cells and neutrophils were significantly related to patients’ prognosis. Finally, we validated these genes in three independent cohorts from the Gene Expression Omnibus database. In conclusion, this study comprehensively analyzed the tumor microenvironment and identified prognostic immune-related biomarkers for melanoma.     

Identification of immune microenvironment subtypes that predicted the prognosis of patients with ovarian cancer

Ovarian cancer (OC) is associated with high mortality rate. However, the correlation between immune microenvironment and prognosis of OC remains unclear. This study aimed to explore prognostic significance of OC tumour microenvironment. The OC data set was selected from the cancer genome atlas (TCGA), and 307 samples were collected. Hierarchical clustering was performed according to the expression of 756 genes. The immune and matrix scores of all immune subtypes were determined, and Kruskal-Wallis test was used to analyse the differences in the immune and matrix scores between OC samples with different immune subtypes. The model for predicting prognosis was constructed based on the expression of immune-related genes. TIDE platform was applied to predict the effect of immunotherapy on patients with OC of different immune subtypes. The 307 OC samples were classified into three immune subtypes A-C. Patients in subtype B had poorer prognosis and lower survival rate. The infiltration of helper T cells and macrophages in microenvironment indicated significant differences between immune subtypes. Enrichment analyses of immune cell molecular pathways showed that JAK–STAT3 pathway changed significantly in subtype B. Furthermore, predictive response to immunotherapy in subtype B was significantly higher than that in subtype A and C. Immune subtyping can be used as an independent predictor of the prognosis of OC patients, which may be related to the infiltration patterns of immune cells in tumour microenvironment. In addition, patients in immune subtype B have superior response to immunotherapy, suggesting that patients in subtype B are suitable for immunotherapy.     

Factors associated with immunotherapy respond and survival in advanced non-small cell lung cancer patients

ObjectivesThis study aimed to explore factors associated with immunotherapy respond and survival in advanced non-small cell lung cancer (aNSCLC) patients treated with immune checkpoint inhibitors (ICIs).MethodsA total of 101 patients with aNSCLC receiving ICIs were included. The association between clinical factors and multiple endpoints including objective response rate (ORR), disease control rate (DCR), overall survival (OS) and progression-free survival (PFS) were investigated by multivariate analyses.ResultsMultivariate logistic analyses revealed that clinical stage, lactate dehydrogena (LDH), and any grade immune-related adverse events (irAEs) were independent predictors of ORR, while LDH and ICIs treatment type were independent predictors of DCR. In Multivariate Cox analysis, Eastern Cooperative Oncology Group performance status (ECOG PS), LDH, albumin (Alb), platelet to lymphocyte ratio (PLR), and any grade irAEs were independent factors for OS. Similarly, clinical stage, LDH, Alb, and any grade irAEs were independent factors for PFS. Pre-treatment prognostic score was established based on clinical stage, ECOG PS, LDH, Alb and PLR to classify patients into three groups: the good group (0–1 score), the intermediate group (2 scores) and the poor group (3-4 scores). The immunotherapy response was significantly different in various prognostic groups. Subset analyses showed pre-treatment prognostic score ≥ 3 tended to have a strong negative impact on survival among patients with programmed cell death-ligand 1 (PD-L1) expression ≥ 50%.ConclusionsPre-treatment prognostic score based on clinical stage, ECOG PS, LDH, Alb and PLR may help to identify aNSCLC patients who may benefit from ICIs.