Development and validation of a novel immune-related prognostic model in hepatocellular carcinoma

Preliminary Analysis of Aging-Related Genes in Intracerebral Hemorrhage by Integration of Bulk and Single-Cell RNA Sequencing Technology

Background

Aging is recognized as the key risk for intracerebral hemorrhage (ICH). The detailed mechanisms of aging in ICH warrant exploration. This study aimed to identify potential aging-related genes associated with ICH.

Methods

ICH-specific aging-related genes were determined by the intersection of differentially expressed genes (DEGs) between perihematomal tissues and corresponding contralateral parts of four patients with ICH (GSE24265) and 349 aging-related genes obtained from the Aging Atlas database. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Set Enrichment Analysis (GSEA) analyses were performed to identify the potential biological functions and pathways in which these ICH-specific aging-related genes may be involved. Then, PPI network was established to identify the hub genes of ICH-specific aging-related genes. Meanwhile, miRNA-mRNA and transcription factor (TF)-mRNA regulatory networks were constructed to further explore the ICH-specific aging-related genes regulation. The relationship between these hub genes and immune infiltration was also further explored. Additional single-cell RNA-seq analysis (scRNA-seq, GSE167593) was used to locate the hub genes in different cell types. Besides, expression levels of the hub genes were validated using clinical samples from our institute and another GEO dataset (GSE206971).

Results

This study identified 24 ICH-specific aging-related genes, including 22 up-regulated and 2 down-regulated genes. The results of GO and KEGG suggested that the ICH-specific aging-related genes mainly enriched in immunity and inflammation-related pathways, suggesting that aging may affect the ich pathogenesis by regulating inflammatory and immune-related pathways.

Conclusion

Our study revealed 24 ICH-specific aging-related genes and their functions highly pertinent to ICH pathogenesis, providing new insights into the impact of aging on ICH.

ANXA5: A Key Regulator of Immune Cell Infiltration in Hepatocellular Carcinoma

BACKGROUND Hepatocellular carcinoma (HCC) poses a significant threat to human life and is the most prevalent form of liver cancer. The intricate interplay between apoptosis, a common form of programmed cell death, and its role in immune regulation stands as a crucial mechanism influencing tumor metastasis. MATERIAL AND METHODS Utilizing HCC samples from the TCGA database and 61 anoikis-related genes (ARGs) sourced from GeneCards, we analyzed the relationship between ARGs and immune cell infiltration in HCC. Subsequently, we identified long non-coding RNAs (lncRNAs) associated with ARGs, using the least absolute shrinkage and selection operator (LASSO) regression analysis to construct a robust prognostic model. The predictive capabilities of the model were then validated through examination in a single-cell dataset. RESULTS Our constructed prognostic model, derived from lncRNAs linked to ARGs, comprised 11 significant lncRNAs: NRAV, MCM3AP-AS1, OTUD6B-AS1, AC026356.1, AC009133.1, DDX11-AS1, AC108463.2, MIR4435-2HG, WARS2-AS1, LINC01094, and HCG18. The risk score assigned to HCC samples demonstrated associations with immune indicators and the infiltration of immune cells. Further, we identified Annexin A5 (ANXA5) as the pivotal gene among ARGs, with it exerting a prominent role in regulating the lncRNA gene signature. Our validation in a single-cell database elucidated the involvement of ANXA5 in immune cell infiltration, specifically in the regulation of mononuclear cells. CONCLUSIONS This study delves into the intricate correlation between ARGs and immune cell infiltration in HCC, culminating in the development of a novel prognostic model reliant on 11 ARGs-associated lncRNAs. Furthermore, our findings highlight ANXA5 as a promising target for immune regulation in HCC, offering new perspectives for immune therapy in the context of HCC.

Development and experimental validation of an M2 macrophage and platelet-associated gene signature to predict prognosis and immunotherapy sensitivity in bladder cancer

Platelets and M2 macrophages both play crucial roles in tumorigenesis, but their relationship and the prognosis value of the relative genes in bladder cancer (BLCA) remain obscure. In the present study, we found that platelets stimulated by BLCA cell lines could promote M2 macrophage polarization, and platelets were significantly associated with the infiltration of M2 macrophages in BLCA samples. Through the bioinformatic analyses, A2M, TGFB3, and MYLK, which were associated with platelets and M2 macrophages, were identified and verified in vitro and then included in the predictive model. A platelet and M2 macrophage-related gene signature was constructed to evaluate the prognosis and immunotherapeutic sensitivity, helping to guide personalized treatment and to disclose the underlying mechanisms.

TuBG1 promotes hepatocellular carcinoma via ATR/P53-apoptosis and cycling pathways

BACKGROUND

As reported, γ-tubulin (TuBG1) is related to the occurrence and development of various types of malignant tumors. However, its role in hepatocellular cancer (HCC) is not clear. The present study was to investigate the relationship between TuBG1 and clinical parameters and survival in HCC patients.

METHODS

The correlation between TuBG1 and clinical parameters and survival in HCC patients was explored by bioinformatics analysis. Immunohistochemistry was used for the verification. The molecular function of TuBG1 was measured using colony formation, scratch assay, trans-well assay and flow cytometry. Gene set enrichment analysis (GSEA) was used to pick up the enriched pathways, followed by investigating the target pathways using Western blotting. The tumor-immune system interactions and drug bank database (TISIDB) was used to evaluate TuBG1 and immunity. Based on the TuBG1-related immune genes, a prognostic model was constructed and was further validated internally and externally.

RESULTS

The bioinformatic analysis found high expressed TuBG1 in HCC tissue, which was confirmed using immunohistochemistry and Western blotting. After silencing the TuBG1 in HCC cell lines, more G1 arrested cells were found, cell proliferation and invasion were inhibited, and apoptosis was promoted. Furthermore, the silence of TuBG1 increased the expressions of Ataxia-Telangiectasia and Rad-3 (ATR), phospho-P38 mitogen-activated protein kinase (P-P38MAPK), phospho-P53 (P-P53), B-cell lymphoma-2 associated X protein (Bax), cleaved caspase 3 and P21; decreased the expressions of B-cell lymphoma-2 (Bcl-2), cyclin D1, cyclin E2, cyclin-dependent kinase 2 (CDK2) and CDK4. The correlation analysis of immunohistochemistry and clinical parameters and survival data revealed that TuBG1 was negatively correlated with the overall survival. The constructed immune prognosis model could effectively evaluate the prognosis.

CONCLUSIONS

The increased expression of TuBG1 in HCC is associated with poor prognosis, which might be involved in the occurrence and development of HCC.

Identification of an immune-related gene signature for predicting prognosis and immunotherapy efficacy in liver cancer via cell-cell communication

BACKGROUND

Liver cancer is one of the deadliest malignant tumors worldwide. Immunotherapy has provided hope to patients with advanced liver cancer, but only a small fraction of patients benefit from this treatment due to individual differences. Identifying immune-related gene signatures in liver cancer patients not only aids physicians in cancer diagnosis but also offers personalized treatment strategies, thereby improving patient survival rates. Although several methods have been developed to predict the prognosis and immunotherapeutic efficacy in patients with liver cancer, the impact of cell-cell interactions in the tumor microenvironment has not been adequately considered.

AIM

To identify immune-related gene signals for predicting liver cancer prognosis and immunotherapy efficacy.

METHODS

Cell grouping and cell-cell communication analysis were performed on single-cell RNA-sequencing data to identify highly active cell groups in immune-related pathways. Highly active immune cells were identified by intersecting the highly active cell groups with B cells and T cells. The significantly differentially expressed genes between highly active immune cells and other cells were subsequently selected as features, and a least absolute shrinkage and selection operator (LASSO) regression model was constructed to screen for diagnostic-related features. Fourteen genes that were selected more than 5 times in 10 LASSO regression experiments were included in a multivariable Cox regression model. Finally, 3 genes (stathmin 1, cofilin 1, and C-C chemokine ligand 5) significantly associated with survival were identified and used to construct an immune-related gene signature.

RESULTS

The immune-related gene signature composed of stathmin 1, cofilin 1, and C-C chemokine ligand 5 was identified through cell-cell communication. The effectiveness of the identified gene signature was validated based on experimental results of predictive immunotherapy response, tumor mutation burden analysis, immune cell infiltration analysis, survival analysis, and expression analysis.

CONCLUSION

The findings suggest that the identified gene signature may contribute to a deeper understanding of the activity patterns of immune cells in the liver tumor microenvironment, providing insights for personalized treatment strategies.

Adverse clinical outcomes and immunosuppressive microenvironment of RHO-GTPase activation pattern in hepatocellular carcinoma

BACKGROUND

Emerging evidence suggests that Rho GTPases play a crucial role in tumorigenesis and metastasis, but their involvement in the tumor microenvironment (TME) and prognosis of hepatocellular carcinoma (HCC) is not well understood.

METHODS

We aim to develop a tumor prognosis prediction system called the Rho GTPases-related gene score (RGPRG score) using Rho GTPase signaling genes and further bioinformatic analyses.

RESULTS

Our work found that HCC patients with a high RGPRG score had significantly worse survival and increased immunosuppressive cell fractions compared to those with a low RGPRG score. Single-cell cohort analysis revealed an immune-active TME in patients with a low RGPRG score, with strengthened communication from T/NK cells to other cells through MIF signaling networks. Targeting these alterations in TME, the patients with high RGPRG score have worse immunotherapeutic outcomes and decreased survival time in the immunotherapy cohort. Moreover, the RGPRG score was found to be correlated with survival in 27 other cancers. In vitro experiments confirmed that knockdown of the key Rho GTPase-signaling biomarker SFN significantly inhibited HCC cell proliferation, invasion, and migration.

CONCLUSIONS

This study provides new insight into the TME features and clinical use of Rho GTPase gene pattern at the bulk-seq and single-cell level, which may contribute to guiding personalized treatment and improving clinical outcome in HCC.

Comprehensive analysis of a tryptophan metabolism-related model in the prognostic prediction and immune status for clear cell renal carcinoma

BACKGROUND

Clear cell renal cell carcinoma (ccRCC) is characterized as one of the most common types of urological cancer with high degrees of malignancy and mortality. Due to the limited effectiveness of existing traditional therapeutic methods and poor prognosis, the treatment and therapy of advanced ccRCC patients remain challenging. Tryptophan metabolism has been widely investigated because it significantly participates in the malignant traits of multiple cancers. The functions and prognostic values of tryptophan metabolism-related genes (TMR) in ccRCC remain virtually obscure.

METHODS

We employed the expression levels of 40 TMR genes to identify the subtypes of ccRCC and explored the clinical characteristics, prognosis, immune features, and immunotherapy response in the subtypes. Then, a model was constructed for the prediction of prognosis based on the differentially expressed genes (DEGs) in the subtypes from the TCGA database and verified using the ICGC database. The prediction performance of this model was confirmed by the receiver operating characteristic (ROC) curves. The relationship of Risk Score with the infiltration of distinct tumor microenvironment cells, the expression profiles of immune checkpoint genes, and the treatment benefits of immunotherapy and chemotherapy drugs were also investigated.

RESULTS

The two subtypes revealed dramatic differences in terms of clinical characteristics, prognosis, immune features, and immunotherapy response. The constructed 6-gene-based model showed that the high Risk Score was significantly connected to poor overall survival (OS) and advanced tumor stages. Furthermore, increased expression of CYP1B1, KMO, and TDO2 was observed in ccRCC tissues at the translation levels, and an unfavorable prognosis for these patients was also found.

CONCLUSION

We identified 2 molecular subtypes of ccRCC based on the expression of TMR genes and constructed a prognosis-related model that may be used as a powerful tool to guide the prediction of ccRCC prognosis and personalized therapy. In addition, CYP1B1, KMO, and TDO2 can be regarded as the risk prognostic genes for ccRCC.

Unleashing the potential of gene signatures as prognostic and predictive tools: A step closer to personalized medicine in hepatocellular carcinoma (HCC)

Hepatocellular carcinoma (HCC) is one of the growing malignancies globally, affecting a myriad of people and causing numerous cancer-related deaths. Despite therapeutic improvements in treatment strategies over the past decades, HCC still remains one of the leading causes of person-years of life lost. Numerous studies have been conducted to assess the characteristics of HCC with the aim of predicting its prognosis and responsiveness to treatment. However, the identified biomarkers have shown limited sensitivity, and the translation of these findings into clinical practice has faced challenges. The development of sequencing techniques has facilitated the exploration of a wide range of genes, leading to the emergence of gene signatures. Although several studies assessed differentially expressed genes in normal and HCC tissues to find the unique gene signature with prognostic value, to date, no study has reviewed the task, and to the best of our knowledge, this review represents the first comprehensive analysis of relevant studies in HCC. Most gene signatures focused on immune-related genes, while others investigated genes related to metabolism, autophagy, and apoptosis. Even though no identical gene signatures were found, NDRG1, SPP1, BIRC5, and NR0B1 were the most extensively studied genes with prognostic value. Finally, despite challenges such as the lack of consistent patterns in gene signatures, we believe that comprehensive analysis of pertinent gene signatures will bring us a step closer to personalized medicine in HCC, where treatment strategies can be tailored to individual patients based on their unique molecular profiles.

Comprehensive analysis of clinical prognosis and biological significance of CNIH4 in cervical cancer

BACKGROUND

Cornichon homolog 4 (CNIH4) belongs to the CNIH family. It functions as an oncogene in many tumors. However, CNIH4's significance in the immune landscape and its predictive potential in cervical cancer (CESC) is unexplored.

METHODS

CNIH4 levels and its effect on the survival of patients with CESC were evaluated using data retrieved from The Cancer Genome Atlas (TCGA). The oncogenic effect of CNIH4 in CESC was determined using small interfering RNA-mediated transfected cell lines and tumorigenesis experiments in animal models.

RESULTS

Higher expression of CNIH4 was found in advanced tumor and pathological stages, as well as lymph node metastasis. CNIH4 expression correlated positively with the infiltration of macrophages M2 and resting dendritic cells into the affected tissue. Additionally, functional enrichment of RNA-sequencing of CNIH4-knocked down CESC cell lines showed the association of CNIH4 to the PI3K-Akt signaling pathway. Single-sample gene set enrichment analysis highlighted several immune pathways that were elevated in the CESC samples with enhanced levels of CNIH4, including Type-I and Type-II IFN-response pathways. The impact of CNIH4 on drug sensitivity was further assessed using the GDSC database. As CNIH4 is linked to the immune landscape in CESC, this study determined a four-gene risk prediction signature utilizing CNIH4-related immunomodulators. The risk score quantified from the prediction signature was an independent predictive indicator in CESC. Receiver operating characteristic curve analysis verified the good predictive ability of the four-gene signature in TCGA-CESC cohort. Thus, the CNIH4-related model showed potential as an auxiliary TNM staging system tool.

CONCLUSION

CNIH4 may be an effective predictive biomarker for patients with cervical cancer, thus providing new ideas and research directions for CESC.

CDK1-SRC Interaction-Dependent Transcriptional Activation of HSP90AB1 Promotes Antitumor Immunity in Hepatocellular Carcinoma

This study aimed to analyze multiomics data and construct a regulatory network involving kinases, transcription factors, and immune genes in hepatocellular carcinoma (HCC) prognosis. The researchers used transcriptomic, proteomic, and clinical data from TCGA and GEO databases to identify immune genes associated with HCC. Statistical analysis, meta-analysis, and protein-protein interaction analyses were performed to identify key immune genes and their relationships. In vitro and in vivo experiments validated the CDK1-SRC-HSP90AB1 network's effects on HCC progression and antitumor immunity. A prognostic risk model was developed using clinicopathological features and immune infiltration. The immune genes LPA, BIRC5, HSP90AB1, ROBO1, and CCL20 were identified as the key prognostic factors. The CDK1-SRC-HSP90AB1 network promoted HCC cell proliferation and migration, with HSP90AB1 being transcriptionally activated by the CDK1-SRC interaction. Manipulating SRC or HSP90AB1 reversed the effects of CDK1 and SRC on HCC. The CDK1-SRC-HSP90AB1 network also influenced HCC tumor formation and antitumor immunity. Overall, this study highlights the importance of the CDK1-SRC-HSP90AB1 network as a crucial immune-regulatory network in the HCC prognosis.

Prediction of liver cancer prognosis based on immune cell marker genes

Introduction

Monitoring the response after treatment of liver cancer and timely adjusting the treatment strategy are crucial to improve the survival rate of liver cancer. At present, the clinical monitoring of liver cancer after treatment is mainly based on serum markers and imaging. Morphological evaluation has limitations, such as the inability to measure small tumors and the poor repeatability of measurement, which is not applicable to cancer evaluation after immunotherapy or targeted treatment. The determination of serum markers is greatly affected by the environment and cannot accurately evaluate the prognosis. With the development of single cell sequencing technology, a large number of immune cell-specific genes have been identified. Immune cells and microenvironment play an important role in the process of prognosis. We speculate that the expression changes of immune cell-specific genes can indicate the process of prognosis.

Method

Therefore, this paper first screened out the immune cell-specific genes related to liver cancer, and then built a deep learning model based on the expression of these genes to predict metastasis and the survival time of liver cancer patients. We verified and compared the model on the data set of 372 patients with liver cancer.

Result

The experiments found that our model is significantly superior to other methods, and can accurately identify whether liver cancer patients have metastasis and predict the survival time of liver cancer patients according to the expression of immune cell-specific genes.

Discussion

We found these immune cell-specific genes participant multiple cancer-related pathways. We fully explored the function of these genes, which would support the development of immunotherapy for liver cancer.

Deciphering the immune heterogeneity dominated by natural killer cells with prognostic and therapeutic implications in hepatocellular carcinoma

Belonging to type 1 innate lymphoid cells (ILC1), natural killer (NK) cells play an important role not only in fighting microbial infections but also in anti-tumor response. Hepatocellular carcinoma (HCC) represents an inflammation-related malignancy and NK cells are enriched in the liver, making them an essential component of the HCC immune microenvironment. In this study, we performed single-cell RNA-sequencing (scRNA-seq) analysis to identify the NK cell marker genes (NKGs) and uncovered 80 prognosis-related ones by the TCGA-LIHC dataset. Based on prognostic NKGs, HCC patients were categorized into two subtypes with distinct clinical outcomes. Subsequently, we conducted LASSO-COX and stepwise regression analysis on prognostic NKGs to establish a five-gene (UBB, CIRBP, GZMH, NUDC, and NCL) prognostic signature-NKscore. Different mutation statuses of the two risk groups stratified by NKscore were comprehensively characterized. Besides, the established NKscore-integrated nomogram presented enhanced predictive performance. Single sample gene set enrichment analysis (ssGSEA) analysis was used to uncover the landscape of the tumor immune microenvironment (TIME) and the high-NKscore risk group was characterized with an immune-exhausted phenotype while the low-NKscore risk group held relatively strong anti-cancer immunity. T cell receptor (TCR) repertoire, tumor inflammation signature (TIS), and Immunophenoscore (IPS) analyses revealed differences in immunotherapy sensitivity between the two NKscore risk groups. Taken together, we developed a novel NK cell-related signature to predict the prognosis and immunotherapy efficacy for HCC patients.

An immune-related signature for optimizing prognosis prediction and treatment decision of hepatocellular carcinoma

BACKGROUND

An immune-related gene signature (IGS) was established for discriminating prognosis, predicting benefit of immunotherapy, and exploring therapeutic options in hepatocellular carcinoma (HCC).

METHODS

Based on Immune-related hub genes and The Cancer Genome Atlas (TCGA) LIHC dataset (n = 363), an immune-related gene signature (IGS) was established by least absolute shrinkage and selection operator (LASSO) analysis. The prognostic significance and clinical implications of IGS were verified in International Cancer Genome Consortium (ICGC) and Chinese HCC (CHCC) cohorts. The molecular and immune characteristics and the benefit of immune checkpoint inhibitor (ICI) therapy in IGS-defined subgroups were analyzed. In addition, by leveraging the Cancer Therapeutics Response Portal (CTRP) and PRISM Repurposing datasets, we determined the potential therapeutic agents for high IGS-risk patients.

RESULTS

The IGS was constructed based on 8 immune-related hub genes with individual coefficients. The IGS risk model could robustly predict the survival of HCC patients in TCGA, ICGC, and CHCC cohorts. Compared with 4 previous established immune genes-based signatures, IGS exhibited superior performance in survival prediction. Additionally, for immunological characteristics and enriched pathways, a low-IGS score was correlated with IL-6/JAK/STAT3 signaling, inflammatory response and interferon α/γ response pathways, low TP53 mutation rate, high infiltration level, and more benefit from ICI therapy. In contrast, high IGS score manifested an immunosuppressive microenvironment and activated aggressive pathways. Finally, by in silico screening potential compounds, vindesine, ispinesib and dasatinib were identified as potential therapeutic agents for high-IGS risk patients.

CONCLUSIONS

This study developed a robust IGS model for survival prediction of HCC patients, providing new insights into integrating tailored risk stratification with precise immunotherapy and screening potentially targeted agents.

Cell Pair Algorithm-Based Immune Infiltrating Cell Signature for Improving Outcomes and Treatment Responses in Patients with Hepatocellular Carcinoma

BACKGROUND

Immune interactions play important roles in the regulation of T cells' cytotoxic function, further impacting the anti-tumor efficacy of immunotherapy. A comprehensive analysis of immune cell types in HCC and immune-cell-related signatures predicting prognosis and monitoring immunotherapy efficacy is still absent.

METHODS

More than 1,300 hepatocellular carcinomas (HCC) patients were collected from public databases and included in the present study. The ssGSEA algorithm was applied to calculate the infiltration level of 28 immunocyte subpopulations. A cell pair algorithm was applied to construct an immune-cell-related prognostic index (ICRPI). Survival analyses were performed to measure the survival difference across ICRPI risk groups. Spearman's correlation analyses were used for the relevance assessment. A Wilcoxon test was used to measure the expression level's differences.

RESULTS

In this study, 28 immune subpopulations were retrieved, and 374 immune cell pairs (ICPs) were established, 38 of which were picked out by the least absolute shrinkage and selection operator (LASSO) algorithm. By using the selected ICPs, the ICRPI was constructed and validated to play crucial roles in survival stratification and dynamic monitoring of immunotherapy effect. We also explored several candidate drugs targeting ICRPI. A composite ICRPI and clinical prognostic index (ICPI) was then constructed, which achieved a more accurate estimation of HCC's survival and is a better choice for prognosis predictions in HCC.

CONCLUSIONS

In conclusion, we constructed and validated ICRPI based on the cell pair algorithm in this study, which might provide some novel insights for increasing the survival estimation and clinical response to immune therapy for individual HCC patients and contribute to the personalized precision immunotherapy strategy of HCC.

A cellular senescence-related classifier based on a tumorigenesis- and immune infiltration-guided strategy can predict prognosis, immunotherapy response, and candidate drugs in hepatocellular carcinoma

BACKGROUND

Cellular senescence plays an irreplaceable role in tumorigenesis, progression, and tumor microenvironment (TME) remodeling. However, to date, there is limited research delineating the landscape of cellular senescence in hepatocellular carcinoma (HCC), and an improved understanding on the interaction of tumor-associated cellular senescence with HCC prognosis, TME, and response to immunotherapy is warrant.

METHODS

Tumorigenic and immune infiltration-associated senescence genes were determined by weighted gene co-expression network analysis (WGCNA) and the Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data (ESTIMATE) algorithm, and subsequently, a prognostic scoring model (named TIS) was constructed using multiple survival analysis algorithms to classify the senescence-related subtypes of HCC. Gene set enrichment analysis (GSEA) and gene set variation analysis (GSVA) were conducted to identify the distinct hallmark pathways between high- and low-risk subtypes. Additionally, we carried out correlation analyses for TIS and clinical traits, senescence-associated secretory phenotype (SASP), immune infiltration and evasion, immune checkpoint factors, drug response, and immunotherapeutic efficacy. External experimental validation was conducted to delineate the association of CPEP3 (a TIS gene) with HCC phenotypes through assays of proliferation, colony formation, and invasion.

RESULTS

A five-gene TIS, composed of NET1, ATP6V0B, MMP1, GTDC1, and CPEB3, was constructed and validated using TCGA and ICGC datasets, respectively, and showed a highly robust and plausible signature for overall survival (OS) prediction of HCC in both training and validation cohorts. Patients in the TIS-high group were accompanied by worse OS, activation of carcinogenetic pathways, infiltration of immunosuppressive cells, exclusion of effector killing cells, overexpression of immunomodulatory genes and SASP, and unsatisfied response to immunotherapy. In response to anticancer drugs, patients in the TIS-high group exhibited enhanced susceptibility to several conventional chemotherapeutic agents (5-fluorouracil, docetaxel, doxorubicin, gemcitabine, and etoposide), as well as several inhibitors of pathways involved in cellular senescence (cell-cycle inhibitors, bromodomain and extraterminal domain family (BET) inhibitors, PI3K-AKT pathway inhibitors, and multikinase inhibitors). Additionally, four putative drugs (palbociclib, JAK3 inhibitor VI, floxuridine, and lestaurtinib) were identified as potential compounds for patients in the TIS-high group. Notably, in vitro functional validation showed that CPEB3 knockdown boosted the phenotypes of proliferation, clonogenicity, and invasion in HCC cells, whereas CPEB3 overexpression attenuated these phenotypes.

CONCLUSIONS

Our study provides comprehensive clues demonstrating the role of novel TIS in predicting HCC prognosis, immunotherapeutic response, and candidate drugs. This work highlights the significance of tumorigenesis- and immune infiltration-related cellular senescence in cancer therapy.

Establishment and analysis of a disease risk prediction model for the systemic lupus erythematosus with random forest

Systemic lupus erythematosus (SLE) is a latent, insidious autoimmune disease, and with the development of gene sequencing in recent years, our study aims to develop a gene-based predictive model to explore the identification of SLE at the genetic level. First, gene expression datasets of SLE whole blood samples were collected from the Gene Expression Omnibus (GEO) database. After the datasets were merged, they were divided into training and validation datasets in the ratio of 7:3, where the SLE samples and healthy samples of the training dataset were 334 and 71, respectively, and the SLE samples and healthy samples of the validation dataset were 143 and 30, respectively. The training dataset was used to build the disease risk prediction model, and the validation dataset was used to verify the model identification ability. We first analyzed differentially expressed genes (DEGs) and then used Lasso and random forest (RF) to screen out six key genes (OAS3, USP18, RTP4, SPATS2L, IFI27 and OAS1), which are essential to distinguish SLE from healthy samples. With six key genes incorporated and five iterations of 10-fold cross-validation performed into the RF model, we finally determined the RF model with optimal mtry. The mean values of area under the curve (AUC) and accuracy of the models were over 0.95. The validation dataset was then used to evaluate the AUC performance and our model had an AUC of 0.948. An external validation dataset (GSE99967) with an AUC of 0.810, an accuracy of 0.836, and a sensitivity of 0.921 was used to assess the model's performance. The external validation dataset (GSE185047) of all SLE patients yielded an SLE sensitivity of up to 0.954. The final high-throughput RF model had a mean value of AUC over 0.9, again showing good results. In conclusion, we identified key genetic biomarkers and successfully developed a novel disease risk prediction model for SLE that can be used as a new SLE disease risk prediction aid and contribute to the identification of SLE.

Hallmark-guided subtypes of hepatocellular carcinoma for the identification of immune-related gene classifiers in the prediction of prognosis, treatment efficacy, and drug candidates

Hepatocellular carcinoma (HCC), accounting for ~90% of all primary liver cancer, is a prevalent malignancy worldwide. The intratumor heterogeneity of its causative etiology, histology, molecular landscape, and immune phenotype makes it difficult to precisely recognize individuals with high mortality risk or tumor-intrinsic treatment resistance, especially immunotherapy. Herein, we comprehensively evaluated the activities of cancer hallmark gene sets and their correlations with the prognosis of HCC patients using gene set variation analysis (GSVA) and identified two HCC subtypes with distinct prognostic outcomes. Based on these subtypes, seven immune-related genes (TMPRSS6, SPP1, S100A9, EPO, BIRC5, PLXNA1, and CDK4) were used to construct a novel prognostic gene signature [hallmark-guided subtypes-based immunologic signature (HGSIS)] via multiple statistical approaches. The HGSIS-integrated nomogram suggested an enhanced predictive performance. Interestingly, oncogenic hallmark pathways were significantly enriched in the high-risk group and positively associated with the risk score. Distinct mutational landscapes and immune profiles were observed between different risk groups. Moreover, immunophenoscore (IPS) and tumor immune dysfunction and exclusion (TIDE) analysis showed different sensitivities of HGSIS risk groups for immune therapy efficacy, and the pRRophetic algorithm indicated distinguishable responses for targeted/chemotherapies in different groups. KIF2C was picked out as the key target concerning HGSIS, and the top 10 small molecules were predicted to bind to the active site of KIF2C via molecular docking, which might be further used for candidate drug discovery of HCC. Taken together, our study offers novel insights for clinically significant subtype recognition, and the proposed signature may be a helpful guide for clinicians to improve the treatment regimens.

Bioinformatic analysis of RNA-seq data from TCGA database reveals prognostic significance of immune-related genes in colon cancer

The tumor immune microenvironment is of crucial importance in cancer progression and anticancer immune responses. Thus, systematic exploration of the expression landscape and prognostic significance of immune-related genes (IRGs) to assist in the prognosis of colon cancer is valuable and significant. The transcriptomic data of 470 colon cancer patients were obtained from The Cancer Genome Atlas database and the differentially expressed genes were analyzed. After an intersection analysis, the hub IRGs were identified and a prognostic index was further developed using multivariable Cox analysis. In addition, the discriminatory ability and prognostic significance of the constructed model were validated and the characteristics of IRGs associated overall survival were analyzed to elucidate the underlying molecular mechanisms. A total of 465 differentially expressed IRGs and 130 survival-associated IRGs were screened. Then, 46 hub IRGs were identified by an intersection analysis. A regulatory network displayed that most of these genes were unfavorable for the prognosis of colon cancer and were regulated by transcription factors. After a least absolute shrinkage and selection operator regression analysis, 14 hub IRGs were ultimately chose to construct a prognostic index. The validation results illustrated that this model could act as an independent indicator to moderately separate colon cancer patients into low- and high-risk groups. This study ascertained the prognostic significance of IRGs in colon cancer and successfully constructed an IRG-based prognostic signature for clinical prediction. Our results provide promising insight for the exploration of diagnostic markers and immunotherapeutic targets in colon cancer.

Construction of Lymph Node Metastasis-Related Prognostic Model and Analysis of Immune Infiltration Mode in Lung Adenocarcinoma

Background

Lung adenocarcinoma (LUAD) is a major cause for global cancer-related deaths. Research reports demonstrate that lymph node metastasis (LNM) is pertinent to the survival rate of LUAD patients, and crux lies in the lack of biomarkers that could distinguish patients with LNM. We aimed to verify the LNM-related prognostic biomarkers in LUAD.

Methods

We firstly accessed the expression data of mRNA from The Cancer Genome Atlas (TCGA) database and then obtained samples with LNM (N+) and without LNM (N-). Differential expression analysis was conducted to acquire differentially expressed genes (DEGs). Univariate-LASSO-multivariate Cox regression analyses were performed on DEGs to build a risk model and obtain optimal genes. Afterwards, effectiveness and independence of risk model were assessed based on TCGA-LUAD and GSE31210 datasets. Moreover, a nomogram was established combining clinical factors and riskscores. Nomogram performance was measured by calibration curves. The infiltration abundance of immune cells was scored with CIBERSORT to explore the differences between high- and low-risk groups. Lastly, gene set enrichment analysis (GSEA) was used to investigate differences in immune features between the two risk groups.

Results

Nine optimal feature genes closely related to LNM in LUAD were identified to construct a risk model. Prognostic ability of the risk model was verified in independent databases. Patients were classified into high- and low-risk groups in accordance with their median riskscores. CIBERSORT score displayed differences in immune cell infiltration like T cells CD4 memory resting between high/low-risk groups. LNM-related genes may also be closely relevant to immune features. Additionally, GSEA indicated that differential genes in the two risk groups were enriched in genes related to immune cells.

Conclusion

This research built a risk model including nine optimal feature genes, which may be potential biomarkers for LUAD.

RNA methylation regulators contribute to poor prognosis of hepatocellular carcinoma associated with the suppression of bile acid metabolism: a multi-omics analysis

RNA methylation has been known to promote the initiation and progression of many types of cancer, including hepatocellular carcinoma (HCC). To fully understand the importance of this post-transcriptional modification in HCC, a thorough investigation that combines different patterns of RNA methylation is urgently needed. In this study, we investigated the regulators of the three most common types of RNA methylation: m6A, N1-methyladenosine (m1A) and 5-methylcytosine (m5C). Based on the genomic and proteomic data, we constructed a classifier consisting of seven RNA methylation regulators. This classifier performed well and robustly predicted the prognosis of HCC patients. By analysis using this classifier, we found that the primary bile acid biosynthesis pathway was mostly downregulated in high-risk HCC patients. Furthermore, we found that the gene expression patterns regulated by several bile acids were similar to those regulated by some well-defined anti-tumor compounds, indicating that bile acid metabolism plays a crucial role in the progression of HCC, and the related metabolites can be used as the potential agents for HCC treatments. Moreover, our study revealed a crosstalk between RNA methylation and bile acid regulators, demonstrating a novel mechanism of the downregulation of bile acid metabolism in HCC and providing new insights into how RNA methylation regulators affect the oncogenesis of HCC.

The DDR-related gene signature with cell cycle checkpoint function predicts prognosis, immune activity, and chemoradiotherapy response in lung adenocarcinoma

BACKGROUND

As a DNA surveillance mechanism, cell cycle checkpoint has recently been discovered to be closely associated with lung adenocarcinoma (LUAD) prognosis. It is also an essential link in the process of DNA damage repair (DDR) that confers resistance to radiotherapy. Whether genes that have both functions play a more crucial role in LUAD prognosis remains unclear.

METHODS

In this study, DDR-related genes with cell cycle checkpoint function (DCGs) were selected to investigate their effects on the prognosis of LUAD. The TCGA-LUAD cohort and two GEO external validation cohorts (GSE31210 and GSE42171) were performed to construct a prognosis model based on the least absolute shrinkage and selection operator (LASSO) regression. Patients were divided into high-risk and low-risk groups based on the model. Subsequently, the multivariate COX regression was used to construct a prognostic nomogram. The ssGSEA, CIBERSORT algorithm, TIMER tool, CMap database, and IC50 of chemotherapeutic agents were used to analyze immune activity and responsiveness to chemoradiotherapy.

RESULTS

4 DCGs were selected as prognostic signatures, and patients in the high-risk group had a lower overall survival (OS). The lower infiltration levels of immune cells and the higher expression levels of immune checkpoints appeared in the high-risk group. The damage repair pathways were upregulated, and chemotherapeutic agent sensitivity was poor in the high-risk group.

CONCLUSIONS

The 4-DCGs signature prognosis model we constructed could predict the survival rate, immune activity, and chemoradiotherapy responsiveness of LUAD patients.

Identification of multi-omics biomarkers and construction of the novel prognostic model for hepatocellular carcinoma

Genome changes play a crucial role in carcinogenesis, and many biomarkers can be used as effective prognostic indicators in various tumors. Although previous studies have constructed many predictive models for hepatocellular carcinoma (HCC) based on molecular signatures, the performance is unsatisfactory. Because multi-omics data can more comprehensively reflect the biological phenomenon of disease, we hope to build a more accurate predictive model by multi-omics analysis. We use the TCGA to identify crucial biomarkers and construct prognostic models through difference analysis, univariate Cox, and LASSO/stepwise Cox analysis. The performances of predictive models were evaluated and validated through survival analysis, Harrell’s concordance index (C-index), receiver operating characteristic (ROC) curve, and decision curve analysis (DCA). Multiple mRNAs, lncRNAs, miRNAs, CNV genes, and SNPs were significantly associated with the prognosis of HCC. We constructed five single-omic models, and the mRNA and lncRNA models showed good performance with c-indexes over 0.70. The multi-omics model presented a robust predictive ability with a c-index over 0.77. This study identified many biomarkers that may help study underlying carcinogenesis mechanisms in HCC. In addition, we constructed multiple single-omic models and an integrated multi-omics model that may provide practical and reliable guides for prognosis assessment.

Prognostic models for outcome prediction in patients with advanced hepatocellular carcinoma treated by systemic therapy: a systematic review and critical appraisal

Objective

To describe and analyze the predictive models of the prognosis of patients with hepatocellular carcinoma (HCC) undergoing systemic treatment.

Design

Systematic review.

Data sources

PubMed and Embase until December 2020 and manually searched references from eligible articles.

Eligibility criteria for study selection

The development, validation, or updating of prognostic models of patients with HCC after systemic treatment.

Results

The systematic search yielded 42 eligible articles: 28 articles described the development of 28 prognostic models of patients with HCC treated with systemic therapy, and 14 articles described the external validation of 32 existing prognostic models of patients with HCC undergoing systemic treatment. Among the 28 prognostic models, six were developed based on genes, of which five were expressed in full equations; the other 22 prognostic models were developed based on common clinical factors. Of the 28 prognostic models, 11 were validated both internally and externally, nine were validated only internally, two were validated only externally, and the remaining six models did not undergo any type of validation. Among the 28 prognostic models, the most common systemic treatment was sorafenib (n = 19); the most prevalent endpoint was overall survival (n = 28); and the most commonly used predictors were alpha-fetoprotein (n = 15), bilirubin (n = 8), albumin (n = 8), Child–Pugh score (n = 8), extrahepatic metastasis (n = 7), and tumor size (n = 7). Further, among 32 externally validated prognostic models, 12 were externally validated > 3 times.

Conclusions

This study describes and analyzes the prognostic models developed and validated for patients with HCC who have undergone systemic treatment. The results show that there are some methodological flaws in the model development process, and that external validation is rarely performed. Future research should focus on validating and updating existing models, and evaluating the effects of these models in clinical practice.

Systematic review registration

PROSPERO CRD42020200187.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-022-09841-5.

A Predictive Model for Prognosis and Therapeutic Response in Hepatocellular Carcinoma Based on a Panel of Three MED8-Related Immunomodulators

The current tumor-node-metastasis (TNM) system is limited in predicting the survival and guiding the treatment of hepatocellular carcinoma (HCC) patients since the TNM system only focuses on the anatomical factors, regardless of the intratumoral molecule heterogeneity. Besides, the landscape of intratumoral immune genes has emerged as a prognostic indicator. The mediator complex subunit 8 (MED8) is a major polymerase regulator and has been described as an oncogene in renal cell carcinoma, but its pathophysiological significance of HCC and its contribution to the prognosis of HCC remain unclear. Here, we aimed to discuss the expression profile and clinical correlation of MED8 in HCC and construct a predictive model based on MED8-related immunomodulators as a supplement to the TNM system. According to our analyses, MED8 was overexpressed in HCC tissues and increased expression of MED8 was an indicator of poor outcome in HCC. The knockdown of MED8 weakened the proliferation, colony forming, and migration of HepG2 and Huh7 cells. Subsequently, a predictive model was identified based on a panel of three MED8-related immunomodulators using The Cancer Genome Atlas (TCGA) database and further validated in International Cancer Genome Consortium (ICGC) database. The combination of the predictive model and the TNM system could improve the performance in predicting the survival of HCC patients. High-risk patients had poor overall survival in TCGA and ICGC databases, as well as in subgroup analysis with early clinicopathology classification. It was also found that high-risk patients had a higher probability of recurrence in TCGA cohort. Furthermore, low-risk score indicated a better response to immunotherapy and drug therapy. This predictive model can be served as a supplement to the TNM system and may have implications in prognosis stratification and therapeutic guidance for HCC.

Identification of a Tumor Immunological Phenotype-Related Gene Signature for Predicting Prognosis, Immunotherapy Efficacy, and Drug Candidates in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is the predominant subtype of primary liver cancer and represents a highly heterogeneous disease, making it hard to predict the prognosis and therapy efficacy. Here, we established a novel tumor immunological phenotype-related gene index (TIPRGPI) consisting of 11 genes by Univariate Cox regression and the least absolute shrinkage and selection operator (LASSO) algorithm to predict HCC prognosis and immunotherapy response. TIPRGPI was validated in multiple datasets and exhibited outstanding performance in predicting the overall survival of HCC. Multivariate analysis verified it as an independent predictor and a TIPRGPI-integrated nomogram was constructed to provide a quantitative tool for clinical practice. Distinct mutation profiles, hallmark pathways, and infiltration of immune cells in tumor microenvironment were shown between the TIPRGPI high and low-risk groups. Notably, significant differences in tumor immunogenicity and tumor immune dysfunction and exclusion (TIDE) were observed between the two risk groups, suggesting a better response to immune checkpoint blockade (ICB) therapy of the low-risk group. Besides, six potential drugs binding to the core target of the TIPRGPI signature were predicted via molecular docking. Taken together, our study shows that the proposed TIPRGPI was a reliable signature to predict the risk classification, immunotherapy response, and drugs candidate with potential application in the clinical decision and treatment of HCC. The novel “TIP genes”-guided strategy for predicting the survival and immunotherapy efficacy, we reported here, might be also applied to more cancers other than HCC.

Systematic Construction and Validation of an Immune-Related Gene-Based Model to Predict Prognosis for Ovarian Cancer

Ovarian cancer (OC) is a malignancy with poor prognosis, stubborn resistance, and frequent recurrence. Recently, it has been widely recognized that immune-related genes (IRGs) have demonstrated their indispensable importance in the occurrence and progression of OC. Given this, this study aimed to identify IRGs with predictive value and build a prognostic model for a more accurate assessment. First, we obtained transcriptome and clinical information of ovarian samples from both TCGA and GTEx databases. After integration, we figured out 10 genes as immune-related prognostic genes (IRPGs) by performing the univariate Cox regression analysis. Subsequently, we established a TF-associated network to investigate its internal mechanism. The prognosis model consisting of 5 IRPGs was constructed later by lasso regression analysis. The comparison of the score with the clinical factors validated its independence and superiority in OC's prognosis. Moreover, the association between the signature and immune cell infiltration demonstrated its ability to image the immune situation of the tumor microenvironment. Finally, the reliability of the risk model was confirmed by the GEO cohort. Together, our study has constructed an independent prognostic model for OC, which may deepen the understanding of the immune microenvironment and help present novel biomarkers or ideas for targeted therapy.

Protein-based prognostic signature for predicting the survival and immunotherapeutic efficiency of endometrial carcinoma

Background

Endometrial cancer (EC) is the most frequent malignancy of the female genital tract worldwide. Our study aimed to construct an effective protein prognostic signature to predict prognosis and immunotherapy responsiveness in patients with endometrial carcinoma.

Methods

Protein expression data, RNA expression profile data and mutation data were obtained from The Cancer Proteome Atlas (TCPA) and The Cancer Genome Atlas (TCGA). Prognosis-related proteins in EC patients were screened by univariate Cox regression analysis. Least absolute shrinkage and selection operator (LASSO) analysis and multivariate Cox regression analysis were performed to establish the protein-based prognostic signature. The CIBERSORT algorithm was used to quantify the proportions of immune cells in a mixed cell population. The Immune Cell Abundance Identifier (ImmuCellAI) and The Cancer Immunome Atlas (TCIA) web tools were used to predict the response to immunochemotherapy. The pRRophetic algorithm was used to estimate the sensitivity of chemotherapeutic and targeted agents.

Results

We constructed a prognostic signature based on 9 prognostic proteins, which could divide patients into high-risk and low-risk groups with distinct prognoses. A novel prognostic nomogram was established based on the prognostic signature and clinicopathological parameters to predict 1, 3 and 5-year overall survival for EC patients. The results obtained with Clinical Proteomic Tumor Analysis Consortium (CPTAC), Human Protein Atlas (HPA) and immunohistochemical (IHC) staining data from EC samples in our hospital supported the predictive ability of these proteins in EC tumors. Next, the CIBERSORT algorithm was used to estimate the proportions of 22 immune cell types. The proportions of CD8 T cells, T follicular helper cells and regulatory T cells were higher in the low-risk group. Moreover, we found that the prognostic signature was positively associated with high tumor mutation burden (TMB) and high microsatellite instability (MSI-H) status in EC patients. Finally, ImmuCellAI and TCIA analyses showed that patients in the low-risk group were more inclined to respond to immunotherapy than patients in the high-risk group. In addition, drug sensitivity analysis indicated that our signature had potential predictive value for chemotherapeutics and targeted therapy.

Conclusion

Our study constructed a novel prognostic protein signature with robust predictive ability for survival and efficiency in predicting the response to immunotherapy, chemotherapy and targeted therapy. This protein signature represents a promising predictor of prognosis and response to cancer treatment in EC patients.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-022-09402-w.

Overweight/obesity-related transcriptomic signature as a correlate of clinical outcome, immune microenvironment, and treatment response in hepatocellular carcinoma

BACKGROUNDS

The pandemic of overweight and obesity (quantified by body mass index (BMI) ≥ 25) has rapidly raised the patient number of non-alcoholic fatty hepatocellular carcinoma (HCC), and several clinical trials have shown that BMI is associated with the prognosis of HCC. However, whether overweight/obesity is an independent prognostic factor is arguable, and the role of overweight/obesity-related metabolisms in the progression of HCC is scarcely known.

MATERIALS AND METHODS

In the present study, clinical information, mRNA expression profile, and genomic data were downloaded from The Cancer Genome Atlas (TCGA) as a training cohort (TCGA-HCC) for the identification of overweight/obesity-related transcriptome. Machine learning and the Cox regression analysis were conducted for the construction of the overweight/obesity-associated gene (OAG) signature. The Kaplan-Meier curve, receiver operating characteristic (ROC) curve, and the Cox regression analysis were performed to assess the prognostic value of the OAG signature, which was further validated in two independent retrospective cohorts from the International Cancer Genome Consortium (ICGC) and Gene Expression Omnibus (GEO). Subsequently, functional enrichment, genomic profiling, and tumor microenvironment (TME) evaluation were utilized to characterize biological activities associated with the OAG signature. GSE109211 and GSE104580 were retrieved to evaluate the underlying response of sorafenib and transcatheter arterial chemoembolization (TACE) treatment, respectively. The Genomics of Drug Sensitivity in Cancer (GDSC) database was employed for the evaluation of chemotherapeutic response.

RESULTS

Overweight/obesity-associated transcriptome was mainly involved in metabolic processes and noticeably and markedly correlated with prognosis and TME of HCC. Afterward, a novel established OAG signature (including 17 genes, namely, GAGE2D, PDE6A, GABRR1, DCAF8L1, DPYSL4, SLC6A3, MMP3, RIBC2, KCNH2, HTRA3, PDX1, ATHL1, PRTG, SHC4, C21orf29, SMIM32, and C1orf133) divided patients into high and low OAG score groups with distinct prognosis (median overall survival (OS): 24.87 vs. 83.51 months, p < 0.0001), and the values of area under ROC curve (AUC) in predicting 1-, 2-, 3-, and 4-year OS were 0.81, 0.80, 0.83, and 0.85, respectively. Moreover, the OAG score was independent of clinical features and also exhibited a good ability for prognosis prediction in the ICGC-LIHC-JP cohort and GSE54236 dataset. Expectedly, the OAG score was also highly correlated with metabolic processes, especially oxidative-related signaling pathways. Furthermore, abundant enrichment of chemokines, receptors, MHC molecules, and other immunomodulators as well as PD-L1/PD-1 expression among patients with high OAG scores indicated that they might have better responses to immunotherapy. However, probably exclusion of T cells from infiltrating tumors resulting in lower infiltration of effective T cells would restrict immunotherapeutic effects. In addition, the OAG score was significantly associated with the response of sorafenib and TACE treatment.

CONCLUSIONS

Overall, this study comprehensively disclosed the relationship between BMI-guided transcriptome and HCC. Moreover, the OAG signature had the potential clinical applications in the future to promote clinical management and precision medicine of HCC.

Cancer Progression Mediated by CAFs Relating to HCC and Identification of Genetic Characteristics Influencing Prognosis

BACKGROUND

Hepatocellular carcinoma (HCC) is one of the most common malignancies, and although there are several treatment options, the overall results are not satisfactory. Cancer-associated fibroblasts (CAFs) can promote cancer progression through various mechanisms.

METHODS

HCC-associated mRNA data were sourced from The Cancer Genome Atlas database (TCGA) and International Cancer Genome Consortium (ICGC) database. First, the differentially expressed CAF-related genes (CAF-DEGs) were acquired by difference analysis and weighted gene coexpression network analysis (WGCNA). Moreover, a CAF-related risk model was built by Cox analysis. Kaplan-Meier (K-M) curves and receiver operating characteristic (ROC) curves were utilized to evaluate the validity of this risk model. Furthermore, enrichment analysis of differentially expressed genes (DEGs) between the high- and low-risk groups was executed to explore the functions relevant to the risk model. Furthermore, this study compared the differences in immune infiltration, immunotherapy, and drug sensitivity between the high- and low-risk groups. Finally, we verified the mRNA expression levels of selected prognostic genes by quantitative real-time polymerase chain reaction (qRT-PCR).

RESULTS

107 CAF-DEGs were identified in the HCC samples, and five prognosis-related genes (ACTA2, IGJ, CTHRC1, CXCL12, and LAMB1) were obtained by Cox analysis and utilized to build a CAF-related risk model. K-M analysis illustrated a low survival in the high-risk group, and ROC curves revealed that the risk model could accurately predict the 1-, 3-, and 5-year overall survival (OS) of HCC patients. In addition, Cox analysis demonstrated that the risk score was an independent prognostic factor. Enrichment analysis illustrated that DEGs between the high- and low-risk groups were related to immune response, amino acid metabolism, and fatty acid metabolism. Furthermore, risk scores were correlated with the tumor microenvironment, CAF scores, and TIDE scores, and CAF-related marker genes were positively correlated with all five model genes. Notably, the risk model was relevant to the sensitivity of chemotherapy drugs. Finally, the results of qRT-PCR demonstrated that the expression levels of 5 model genes were in accordance with the analysis.

CONCLUSION

A CAF-related risk model based on ACTA2, IGJ, CTHRC1, CXCL12, and LAMB1 was built and could be utilized to predict the prognosis and treatment of HCC.

The Use of Machine Learning to Create a Risk Score to Predict Survival in Patients with Hepatocellular Carcinoma: A TCGA Cohort Analysis

INTRODUCTION

Hepatocellular carcinoma (HCC) accounts for approximately 90% of primary liver malignancies and is currently the fourth most common cause of cancer-related death worldwide. Due to varying underlying etiologies, the prognosis of HCC differs greatly among patients. It is important to develop ways to help stratify patients upon initial diagnosis to provide optimal treatment modalities and follow-up plans. The current study uses Artificial Neural Network (ANN) and Classification Tree Analysis (CTA) to create a gene signature score that can help predict survival in patients with HCC.

METHODS

The Cancer Genome Atlas (TCGA-LIHC) was analyzed for differentially expressed genes. Clinicopathological data were obtained from cBioPortal. ANN analysis of the 75 most significant genes predicting disease-free survival (DFS) was performed. Next, CTA results were used for creation of the scoring system. Cox regression was performed to identify the prognostic value of the scoring system.

RESULTS

363 patients diagnosed with HCC were analyzed in this study. ANN provided 15 genes with normalized importance >50%. CTA resulted in a set of three genes (NRM, STAG3, and SNHG20). Patients were then divided in to 4 groups based on the CTA tree cutoff values. The Kaplan-Meier analysis showed significantly reduced DFS in groups 1, 2, and 3 (median DFS: 29.7 months, 16.1 months, and 11.7 months, p < 0.01) compared to group 0 (median not reached). Similar results were observed when overall survival (OS) was analyzed. On multivariate Cox regression, higher scores were associated with significantly shorter DFS (1 point: HR 2.57 (1.38-4.80), 2 points: 3.91 (2.11-7.24), and 3 points: 5.09 (2.70-9.58), p < 0.01).

CONCLUSION

Long-term outcomes of patients with HCC can be predicted using a simplified scoring system based on tumor mRNA gene expression levels. This tool could assist clinicians and researchers in identifying patients at increased risks for recurrence to tailor specific treatment and follow-up strategies for individual patients.

Immunogenomic Landscape and Immune-Related Gene-Based Prognostic Signature in Asian Gastric Cancer

Background

Asians have the highest incidence of gastric cancer (GC), and the prognosis of Asian GC is poor. Furthermore, the therapeutics for Asian GC is limited because of genetic heterogeneity and screening difficulty at the early stage. This study aimed to develop an immune-related gene (IRG)-based prognostic signature and to explore prognosis-related regulatory mechanism and therapeutic target for Asian GC.

Methods

To elucidate the prognostic value of IRGs in Asian GC, a comprehensive analysis of IRG expression profiles and overall survival times in 364 Asian GC patients from the Asian Cancer Research Group (ACRG) and The Cancer Genome Atlas (TCGA) databases was performed, and a novel prognostic index was established. To further explore regulatory prognosis mechanisms and therapeutic targets, a tumor immunogenomic landscape analysis, including stromal and immune subcomponents, cell types, panimmune gene sets, and immunomodulatory genes, was performed.

Result

Our analysis allowed the creation of an optimal risk assessment model, the Asian-specific IRG-based prognostic index (ASIRGPI), which showed a high accuracy in predicting survival in Asian GC. We also developed an ASIRGPI-based nomogram to predict the 3- and 5-year overall survival (OS) of Asian GC patients. The impact of the ASIRGPI on the worse prognosis of Asian GC was possibly related to the stromal component remodeling. Specifically, TGFβ gene sets were significantly associated with the ASIRGPI and worse prognosis. Immunomodulatory gene analysis further revealed that TGFβ1 and EDNRB may be the novel potential therapeutic targets for Asian GC.

Conclusions

As a tumor microenvironment-relevant gene set-based prognostic signature, the ASIRGPI model provides an effective approach for evaluating the prognosis of Asian GC and may even prolong OS by enabling the selection of individualized therapy with the novel targets.

Emerging Roles on Immunological Effect of Indoleamine 2,3-Dioxygenase in Liver Injuries

Indoleamine 2,3-dioxygenase (IDO) is one of the initial rate-limiting enzymes of the kynurenine pathway (KP), which causes immune suppression and induction of T cell anergy. It is associated with the imbalance of immune homeostasis in numerous diseases including cancer, chronic viral infection, allergy, and autoimmune diseases. Recently, IDO has extended its role to liver field. In this review, we summarize the dysregulation and potentials of IDO in the emerging field of liver injuries, as well as current challenges for IDO targets. In particular, we discuss unexpected conclusions against previous work published. IDO is induced by pro-inflammatory cytokines in liver dysfunction and exerts an immunosuppressive effect, whereas the improvement of liver injury may require consideration of multiple factors besides IDO.

Construction of a microenvironment immune gene model for predicting the prognosis of endometrial cancer

Background

Infiltrating immune and stromal cells are important components of the endometrial cancer (EC) microenvironment, which has a significant effect on the biological behavior of EC, suggesting that unique immune-related genes may be associated with the prognosis of EC. However, the association of immune-related genes with the prognosis of EC has not been elucidated. We attempted to identify immune-related genes with potentially prognostic value in EC using The Cancer Genome Atlas database and the relationship between immune microenvironment and EC.

Methods

We analyzed 578 EC samples from TCGA database and used weighted gene co-expression network analysis to screen out immune-related genes. We constructed a protein–protein interaction network and analyzed it using STRING and Cytoscape. Immune-related genes were analyzed through conjoint Cox regression and random forest algorithm analysis were to identify a multi-gene prediction model and stratify low-risk and high-risk groups of EC patients. Based on these data, we constructed a nomogram prediction model to improve prognosis assessment. Evaluation of Immunological, gene mutations and gene enrichment analysis were applied on these groups to quantify additional differences.

Results

Using conjoint Cox regression and random forest algorithm, we found that TRBC2, TRAC, LPXN, and ARHGAP30 were associated with the prognosis of EC and constructed four gene risk models for overall survival and a consistent nomogram. The time-dependent receiver operating characteristic curve analysis revealed that the area under the curve for 1-, 3-, and 5-y overall survival was 0.687, 0.699, and 0.76, respectively. These results were validated using a validation cohort. Immune-related pathways were mostly enriched in the low-risk group, which had higher levels of immune infiltration and immune status.

Conclusion

Our study provides new insights for novel biomarkers and immunotherapy targets in EC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-08935-w.

CDKN2A is a prognostic biomarker and correlated with immune infiltrates in hepatocellular carcinoma

Cyclin dependent kinase inhibitor 2A (CDKN2A) is an essential regulator of immune cell functionality, but the mechanisms whereby it drives immune infiltration in hepatocellular carcinoma (HCC) remain unclear. In the present study, we studied the association with CDKN2A expression and immune invasion with the risk of developing HCC. A totally of 2207 different genes were found between HCC and adjacent liver tissues from TCGA and GEO databases. CDKN2A was highly expressed in HCC and associated with poorer overall survival and disease-free survival. Notably, CDKN2A expression was positively correlated with infiltrating levels into purity, B cell, CD+8 T cell, CD+4 T cell, macrophage, neutrophil, and dendritic cells in HCC. CDKN2A expression showed strong correlations between diverse immune marker sets in HCC. These findings suggest that CDKN2A expression potentially contributes to regulation of tumor-associated macrophages and can be used as a prognostic biomarker for determining prognosis and immune infiltration in HCC.

Genomic Instability of Mutation-Derived Gene Prognostic Signatures for Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is one of the major cancer-related deaths worldwide. Genomic instability is correlated with the prognosis of cancers. A biomarker associated with genomic instability might be effective to predict the prognosis of HCC. In the present study, data of HCC patients from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) databases were used. A total of 370 HCC patients from the TCGA database were randomly classified into a training set and a test set. A prognostic signature of the training set based on nine overall survival (OS)–related genomic instability–derived genes (SLCO2A1, RPS6KA2, EPHB6, SLC2A5, PDZD4, CST2, MARVELD1, MAGEA6, and SEMA6A) was constructed, which was validated in the test and TCGA and ICGC sets. This prognostic signature showed more accurate prediction for prognosis of HCC compared with tumor grade, pathological stage, and four published signatures. Cox multivariate analysis revealed that the risk score could be an independent prognostic factor of HCC. A nomogram that combines pathological stage and risk score performed well compared with an ideal model. Ultimately, paired differential expression profiles of genes in the prognostic signature were validated at mRNA and protein level using HCC and paratumor tissues obtained from our institute. Taken together, we constructed and validated a genomic instability–derived gene prognostic signature, which can help to predict the OS of HCC and help us to explore the potential therapeutic targets of HCC.

Prognostic Signature of Hepatocellular Carcinoma and Analysis of Immune Infiltration Based on m6A-Related lncRNAs

The relationship between m6A-related lncRNAs and prognosis in hepatocellular carcinoma (HCC) is not yet clear. We used Lasso regression to establish a prognostic signature based on m6A-related lncRNAs using a training set from TCGA, and then verified the signature efficacy in a test set. Fluorescence quantitative real-time PCR (qPCR), Survival analysis, clinical risk difference analysis, immune-related analysis, and drug-sensitivity analysis were conducted. The results revealed that 1,651 lncRNAs were differentially expressed in HCC tissues, among which, 163 were m6A-related. Univariate analysis showed that 87 lncRNAs were associated with the overall survival. Six differential m6A-related lncRNAs were validated and selected via Lasso regression to construct a prognostic signature which demonstrated a satisfactory predictive efficacy. In the clinically relevant pathologic stage, histologic grade, and T stage, the risk scores obtained based on this signature showed a statistically significant difference. The high- and low-risk groups exhibited a difference in the tumor immune infiltrating cells, immune checkpoint gene expression, and sensitivity to chemotherapy. In summary, the prognostic signature based on the m6A-related lncRNAs can effectively predict the prognosis of patients and might provide a new vista for the chemotherapy and immunotherapy of HCC.

Construction of a Prognostic Model for Hepatocellular Carcinoma Based on Immunoautophagy-Related Genes and Tumor Microenvironment

Background

The aim of this study was to screen and identify immunoautophagy-related genes (IARGs) in HCC patients and clarify their potential prognostic value in HCC patients.

Methods

Immune-related genes and autophagy-related gene were downloaded from public databases. Cox regression analysis was used to selected several immunoautophagy-related genes to establish a prognostic model, and patients were divided into high- and low-risk groups based on median risk score. We analyzed the overall survival and clinicopathological characteristics between two groups. Meanwhile, internal validation dataset and external ICGC dataset were used to verify robustness of the model. Associations between six immune cells infiltrates and risk score were analyzed.

Results

A prognostic model was established based on CANX and HDAC1. The prognoses of the high-risk group were worse than low-risk group in both TCGA and ICGC datasets. Multivariate Cox regression analysis showed that risk score was an independent prognostic factor for HCC patients. Results showed that the risk score in young group was higher than elderly group. Patients with poorly differentiated tumor may have high risk score and poor survival. The score was positively correlated with immune cells.

Conclusion

Our study shows that immunoautophagy-related genes have potential prognostic value for patients with HCC and may provide new information and direction for targeted therapy.

Xanthine dehydrogenase as a prognostic biomarker related to tumor immunology in hepatocellular carcinoma

Background

Xanthine dehydrogenase (XDH) is a critical enzyme involved in the oxidative metabolism of purines, pterin and aldehydes and a central component of the innate immune system. However, the prognostic value of XDH in predicting tumor-infiltrating lymphocyte abundance, the immune response, and survival in different cancers, including hepatocellular carcinoma (HCC), is still unclear.

Methods

XDH expression was analyzed in multiple databases, including Oncomine, the Tumor Immune Estimation Resource (TIMER), the Kaplan–Meier plotter database, the Gene Expression Profiling Interactive Analysis (GEPIA) database, and The Cancer Genome Atlas (TCGA). XDH-associated transcriptional profiles were detected with an mRNA array, and the levels of infiltrating immune cells were validated by immunohistochemistry (IHC) of HCC tissues. A predictive signature containing multiple XDH-associated immune genes was established using the Cox regression model.

Results

Decreased XDH mRNA expression was detected in human cancers originating from the liver, bladder, breast, colon, bile duct, kidney, and hematolymphoid system. The prognostic potential of XDH mRNA expression was also significant in certain other cancers, including HCC, breast cancer, kidney or bladder carcinoma, gastric cancer, mesothelioma, lung cancer, and ovarian cancer. In HCC, a low XDH mRNA level predicted poorer overall survival, disease-specific survival, disease-free survival, and progression-free survival. The prognostic value of XDH was independent of the clinical features of HCC patients. Indeed, XDH expression in HCC activated several immune-related pathways, including the T cell receptor, PI3K-AKT, and MAPK signaling pathways, which induced a cytotoxic immune response. Importantly, the microenvironment of XDH^high HCC tumors contained abundant infiltrating CD8 + T cells but not exhausted T cells. A risk prediction signature based on multiple XDH-associated immune genes was revealed as an independent predictor in the TCGA liver cancer cohort.

Conclusion

These findings suggest that XDH is a valuable prognostic biomarker in HCC and other cancers and indicate that it may function in tumor immunology. Loss of XDH expression may be an immune evasion mechanism for HCC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-021-02173-7.

MAPK-RAP1A Signaling Enriched in Hepatocellular Carcinoma Is Associated With Favorable Tumor-Infiltrating Immune Cells and Clinical Prognosis

Background

MAPK-RAP1A signaling, which is involved in cancer progression, remains to be defined. Upregulation of MAPK-RAP1A signaling accounts for most cancers that harbor high incident rate, such as non-small cell lung cancer (NSCLC) and pancreatic cancer, especially in hepatocellular carcinoma (HCC). MAPK-RAP1A signaling plays an important function as clinical diagnosis and prognostic value in cancers, and the role of MAPK-RAP1A signaling related with immune infiltration for HCC should be elucidated.

Methods

Microarray data and patient cohort information from The Cancer Genome Atlas (TCGA; n = 425) and International Cancer Genome Consortium (ICGC; n = 405) were selected for validation. The Cox regression and least absolute shrinkage and selection operator (LASSO) were used to construct a clinical prognostic model in this analysis and validation study. We also tested the area under the curve (AUC) of the risk signature that could reflect the status of predictive power by determining model. MAPK-RAP1A signaling is also associated with tumor-infiltrating immune cells (TICs) as well as clinical parameters in HCC. The GSEA and CIBERSORT were used to calculate the proportion of TICs, which should be beneficial for the clinical characteristics (clinical stage, distant metastasis) and positively correlated with the survival of HCC patients.

Results

HCC patients with enrichment of MAPK-RAP1A signaling were associated with clinical characteristics and favorable T cell gamma delta (Vδ T cells), and STMN1, RAP1A, FLT3, HSPA8, ANGPT2, and PGF were used as candidate biomarkers for risk scores of HCC. To determine the molecular mechanism of this signature gene association, Gene Set Enrichment Analysis (GSEA) was proposed. Cytokine-cytokine receptor interaction, TGF-β signaling pathway, and Intestinal immune network for IgA production gene sets were closely related in MAPK-RAP1A gene sets. Thus, we established a novel prognostic prediction of HCC to deepen learning of MAPK-RAP1A signaling pathways.

Conclusion

Our findings demonstrated that HCC patients with enrichment of MAPK-RAP1A signaling were associated with clinical characteristics and favorable T cell gamma delta (Vδ T cells), which may be a novel prognostic prediction of HCC.

Autophagy-Related Gene Pairs Signature for the Prognosis of Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) has been recognized as the third leading cause of cancer-related deaths worldwide. There is increasing evidence that the abnormal expression of autophagy-related genes plays an important role in the occurrence and development of HCC. Therefore, the study of autophagy-related genes can further elucidate the genetic drivers of cancer and provide valuable therapeutic targets for clinical treatment. In this study, we used 232 autophagy-related genes extracted from the Human Autophagy Database (HADb) and Molecular Signatures Database (MSigDB) to construct 1884 autophagy-related gene pairs. On this basis, we developed a prognostic model based on autophagy-related gene pairs using least absolute shrinkage and selection operator (LASSO) Cox regression to evaluate the prognosis of patients after liver cancer resection. We then used 845 liver cancer samples from three different databases to test the reliability of the risk signature through survival analysis, receiver operating characteristic (ROC) curve analysis, univariate and multivariate analysis. To further explore the underlying biological mechanisms, we conducted an enrichment analysis of autophagy-related genes. Finally, we combined the signature with independent prognostic factors to construct a nomogram. Based on the autophagy-related gene pair (ARGP) signature, we can divide patients into high- or low-risk groups. Survival analysis and ROC curve analysis verified the validity of the signature (AUC: 0.786—0.828). Multivariate Cox regression showed that the risk score can be used as an independent predictor of the clinical outcomes of liver cancer patients. Notably, this model has a more accurate predictive effect than most prognostic models for hepatocellular carcinoma. Moreover, our model is a powerful supplement to the HCC staging indicator, and a nomogram comprising both indicators can provide a better prognostic effect. Based on pairs of multiple autophagy-related genes, we proposed a prognostic model for predicting the overall survival rate of HCC patients after surgery, which is a promising prognostic indicator. This study confirms the importance of autophagy in the occurrence and development of HCC, and also provides potential biomarkers for targeted treatments.

Integrative analysis identifies key mRNA biomarkers for diagnosis, prognosis, and therapeutic targets of HCV-associated hepatocellular carcinoma

Hepatitis C virus-associated HCC (HCV-HCC) is a prevalent malignancy worldwide and the molecular mechanisms are still elusive. Here, we screened 240 differentially expressed genes (DEGs) of HCV-HCC from Gene expression omnibus (GEO) and the Cancer Genome Atlas (TCGA), followed by weighted gene coexpression network analysis (WGCNA) to identify the most significant module correlated with the overall survival. 10 hub genes (CCNB1, AURKA, TOP2A, NEK2, CENPF, NUF2, CDKN3, PRC1, ASPM, RACGAP1) were identified by four approaches (Protein-protein interaction networks of the DEGs and of the significant module by WGCNA, and diagnostic and prognostic values), and their abnormal expressions, diagnostic values, and prognostic values were successfully verified. A four hub gene-based prognostic signature was built using the least absolute shrinkage and selection operator (LASSO) algorithm and a multivariate Cox regression model with the ICGC-LIRI-JP cohort (N =112). Kaplan-Meier survival plots (P = 0.0003) and Receiver Operating Characteristic curves (ROC = 0.778) demonstrated the excellent predictive potential for the prognosis of HCV-HCC. Additionally, upstream regulators including transcription factors and miRNAs of hub genes were predicted, and candidate drugs or herbs were identified. These findings provide a firm basis for the exploration of the molecular mechanism and further clinical biomarkers development of HCV-HCC.

Development and validation of epithelial mesenchymal transition-related prognostic model for hepatocellular carcinoma

Epithelial cell transformation (EMT) plays an important role in the pathogenesis and metastasis of hepatocellular carcinoma (HCC). We aimed to establish a genetic risk model to evaluate HCC prognosis based on the expression levels of EMT-related genes. The data of HCC patients were collected from TCGA and ICGC databases. Gene expression differential analysis, univariate analysis, and lasso combined with stepwise Cox regression were used to construct the prognostic model. Kaplan–Meier curve, receiver operating characteristic (ROC) curve, calibration analysis, Harrell’s concordance index (C-index), and decision curve analysis (DCA) were used to evaluate the predictive ability of the risk model or nomogram. GO and KEGG were used to analyze differently expressed EMT genes, or genes that directly or indirectly interact with the risk-associated genes. A 10-gene signature, including TSC2, ACTA2, SLC2A1, PGF, MYCN, PIK3R1, EOMES, BDNF, ZNF746, and TFDP3, was identified. Kaplan–Meier survival analysis showed a significant prognostic difference between high- and low-risk groups of patients. ROC curve analysis showed that the risk score model could effectively predict the 1-, 3-, and 5-year overall survival rates of patients with HCC. The nomogram showed a stronger predictive effect than clinical indicators. C-index, DCA, and calibration analysis demonstrated that the risk score and nomogram had high accuracy. The single sample gene set enrichment analysis results confirmed significant differences in the types of infiltrating immune cells between patients in the high- and low-risk groups. This study established a new prediction model of risk gene signature for predicting prognosis in patients with HCC, and provides a new molecular tool for the clinical evaluation of HCC prognosis.

Development and Validation of a Metabolic-related Prognostic Model for Hepatocellular Carcinoma

BACKGROUND AND AIMS

Growing evidence suggests that metabolic-related genes have a significant impact on the occurrence and development of hepatocellular carcinoma (HCC). However, the prognostic value of metabolic-related genes for HCC has not been fully revealed.

METHODS

mRNA sequencing and clinical data were obtained from The Cancer Genome Atlas and the GTEx Genotype-Tissue Expression comprehensive database. Differentially expressed metabolic-related genes in tumor tissues (n=374) and normal tissues (n=160) were identified by the Wilcoxon test. Time-dependent receiver operating characteristic curve analysis, univariate multivariate Cox regression analysis and Kaplan-Meier survival analysis were used to evaluate the predictive effectiveness and independence of the prognostic model. Two independent cohorts (International Cancer Genome Consortiums and GSE14520) were applied to verify the prognostic model.

RESULTS

Our study included a total of 793 patients with HCC. We constructed a risk score consisting of five metabolic-genes (BDH1, RRM2, CYP2C9, PLA2G7, and TXNRD1). For the overall survival rate, the low-risk group had a considerably higher rate than the high-risk group. Univariate and multivariate Cox regression analyses indicated that the risk score was an independent predictor for the prognosis of HCC.

CONCLUSIONS

We constructed and validated a novel prognostic model, which may provide support for the precise treatment of HCC.

A novel comprehensive immune-related gene signature as a promising survival predictor for the patients with head and neck squamous cell carcinoma

Head and neck squamous cell carcinoma (HNSCC), the most frequent subtype of head and neck cancer, continues to have a poor prognosis with no improvement. The TNM stage is not satisfactory for individualized prognostic assessment and it does not predict response to therapy. In the present study, we downloaded the gene expression profiles from TCGA database to establish a training set and GEO database for a validation set. In the training set, we developed an 10 immune-related genes signature which had superior predictive value compared with TNM stage. A nomogram including clinical characteristics was also constructed for accurate prediction. Furthermore, it was determined that our prognostic signature might act as an independent factor for predicting the survival of HNSCC patients. As for the immune microenvironment, our results showed higher immune checkpoint expression (CLTA-4 and PD-1) in low-risk group which might reflect a positive immunotherapy response. Thus, our signature not only provided a promising biomarker for survival prediction, but might be evaluated as an indicator for personalized immunotherapy in patients with HNSCC.

Development and validation of the prognostic value of the immune-related genes in clear cell renal cell carcinoma

Background

Clear cell renal cell carcinoma (ccRCC) is a highly heterogeneous tumor, resulting a challenge of developing target therapeutics. Not long ago, immune checkpoint blockade regimens combine with tyrosin kinase inhibitors have evolved frontline options in metastatic RCC, which implies arrival of the era of tumor immunotherapy. Studies have demonstrated immune-related genes (IRGs) could characterize tumor milieu and related to patient survival. Nevertheless, the clinical significance of classifier depending on IRGs in ccRCC has not been well established.

Methods

The R package limma, univariate and LASSO cox regression analysis were used to screen the prognostic related IRGs from TCGA database. Multivariate cox regression was utilized to establish a risk prediction model for candidate genes. Quantitative real-time PCR was used to confirm the expression of candidates in clinical samples from our institution. CIBERSORT algorithm and correlation analysis were applied to explore tumor-infiltrating immune cells signature between different risk groups. A clinical nomogram was also developed to predict OS by using the rms R package based on the risk prediction model and other independent risk factors. The ICGC data was used for external validation of either gene risk model or nomogram.

Results

We identified 382 differentially expressed immune related genes. Four unique prognostic IRGs (CRABP2, LTB4R, PTGER1 and TEK) were finally affirmed to associate with tumor survival independently and utilized to establish the risk score model. All candidates’ expression was successfully laboratory confirmed by q-PCR. CIBERSORT analysis implied patients in unfavorable-risk group with high CD8 T cell, regulatory T cell and NK cell infiltration, as well as high expression of PD-1, CTLA4, TNFRSF9, TIGIT and LAG3. A nomogram combined IRGs risk score with age, gender, TNM stage, Fuhrman grade, necrosis was further generated to predict of 3- and 5-year OS, which exhibited superior discriminative power (AUCs were 0.811 and 0.795).

Conclusions

Our study established and validated a survival prognostic model system based on 4 unique immune related genes in ccRCC, which expands knowledge in tumor immune status and provide a potent prediction tool in future.

Identification of Mitochondrial-Related Prognostic Biomarkers Associated With Primary Bile Acid Biosynthesis and Tumor Microenvironment of Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is one of the leading causes of tumor-associated deaths worldwide. Despite great progress in early diagnosis and multidisciplinary tumor management, the long-term prognosis of HCC remains poor. Currently, metabolic reprogramming during tumor development is widely observed to support rapid growth and proliferation of cancer cells, and several metabolic targets that could be used as cancer biomarkers have been identified. The liver and mitochondria are the two centers of human metabolism at the whole organism and cellular levels, respectively. Thus, identification of prognostic biomarkers based on mitochondrial-related genes (Mito-RGs)—the coding-genes of proteins located in the mitochondria—that reflect metabolic changes associated with HCC could lead to better interventions for HCC patients. In the present study, we used HCC data from The Cancer Genome Atlas (TCGA) database to construct a classifier containing 10 Mito-RGs (ACOT7, ADPRHL2, ATAD3A, BSG, FAM72A, PDK3, PDSS1, RAD51C, TOMM34, and TRMU) for predicting the prognosis of HCC by using 10-fold Least Absolute Shrinkage and Selection Operation (LASSO) cross-validation Cox regression. Based on the risk score calculated by the classifier, the samples were divided into high- and low-risk groups. Gene set enrichment analysis (GSEA), gene set variation analysis (GSVA), t-distributed stochastic neighbor embedding (t-SNE), and consensus clusterPlus algorithms were used to identify metabolic pathways that were significantly different between the high- and low-risk groups. We further investigated the relationship between metabolic status and infiltration of immune cells into HCC tumor samples by using the Cell-type Identification By Estimating Relative Subsets Of RNA Transcripts (CIBERSORT) algorithm combined with the Tumor Immune Estimation Resource (TIMER) database. Our results showed that the classifier based on Mito-RGs could act as an independent biomarker for predicting survival of HCC patients. Repression of primary bile acid biosynthesis plays a vital role in the development and poor prognosis of HCC, which provides a potential approach to treatment. Our study revealed cross-talk between bile acid and infiltration of tumors by immune cells, which may provide novel insight into immunotherapy of HCC. Furthermore, our research may provide a novel method for HCC metabolic therapy based on modulation of mitochondrial function.

Identification of a Novel Four-Gene Signature Correlated With the Prognosis of Patients With Hepatocellular Carcinoma: A Comprehensive Analysis

Purpose

Hepatocellular carcinoma (HCC) is a common solid-tumor malignancy with high heterogeneity, and accurate prognostic prediction in HCC remains difficult. This analysis was performed to find a novel prognostic multigene signature.

Methods

The TCGA-LIHC dataset was analyzed for differentially coexpressed genes through weighted gene coexpression network analysis (WGCNA) and differential gene expression analysis. A protein-protein interaction (PPI) network and univariate Cox regression analysis of overall survival (OS) were utilized to identify their prognostic value. Next, we used least absolute shrinkage and selection operator (LASSO) Cox regression to establish a prognostic module. Subsequently, the ICGC-LIRI-JP dataset was applied for further validation. Based on this module, HCC cases were stratified into high-risk and low-risk groups, and differentially expressed genes (DEGs) were identified. Functional enrichment analyses of these DEGs were conducted. Finally, single-sample gene set enrichment analysis (ssGSEA) was performed to explore the correlation between the prognostic signature and immune status.

Results

A total of 393 differentially coexpressed genes were obtained. Forty differentially coexpressed hub genes were identified using the CytoHubba plugin, and 38 of them were closely correlated with OS. Afterward, we established the four-gene prognostic signature with an acceptable accuracy (area under the curve [AUC] of 1-year survival: 0.739). The ICGC-LIRI-JP dataset also supported the acceptable accuracy (AUC of 1-year survival:0.752). Compared with low-risk cohort, HCC cases in the high-risk cohort had shorter OS, higher tumor grades, and higher T stages. The risk scores of this signature still act as independent predictors of OS (P<0.001). Functional enrichment analyses suggest that it was mainly organelle fission and nuclear division that were enriched. Finally, ssGSEA revealed that this signature is strongly associated with the immune status of HCC patients.

Conclusions

The proposed prognostic signature of four differentially coexpressed hub genes has satisfactory prognostic ability, providing important insight into the prediction of HCC prognosis.

Development of a poor-prognostic-mutations derived immune prognostic model for acute myeloid leukemia

Leukemia cell-intrinsic somatic mutations and cytogenetic abnormalities have been used to define risk categories in acute myeloid leukemia (AML). In addition, since the immune microenvironment might influence prognosis and somatic mutations have been demonstrated to modulate the immune microenvironment in AML, there is need for developing and evaluating an immune prognostic model (IPM) derived from mutations associated with poor prognosis. Based on AML cases with intermediate and adverse-cytogenetic risk in the Cancer Genome Atlas (TCGA) database, 64 immune-related differentially expressed genes (DEGs) among patients with RUNX1, TP53, or ASXL1 mutations and patients without these mutations were identified. After Cox proportional hazards analysis, an IPM composed of PYCARD and PEAR1 genes was constructed. IPM defined high-risk (IPM-HR) independently predicted lower 2-year overall survival (OS) rates in both patients with intermediate and adverse-cytogenetic risks and non-M3 patients in the TCGA AML cohort. The poor prognostic impact of IPM-HR on OS was further validated by GSE71014, 37642, and 10358 downloaded from the Gene Expression Omnibus (GEO) database. Furthermore, IPM-HR was remarkably associated with higher proportions of CD8+ T cells and regulatory T cells (Tregs), lower proportions of eosinophils, and higher expression of the checkpoint molecules CTLA-4, PD-1, and LAG3 in the TCGA non-M3 AML cohort. In summary, we developed and validated an IPM derived from mutations related with poor prognosis in AML, which would provide new biomarkers for patient stratification and personalized immunotherapy.

Six Metabolism Related mRNAs Predict the Prognosis of Patients With Hepatocellular Carcinoma

Background: Hepatocellular carcinoma (HCC) is one of the most common aggressive solid malignant tumors and current research regards HCC as a type of metabolic disease. This study aims to establish a metabolism-related mRNA signature model for risk assessment and prognosis prediction in HCC patients. Methods: HCC data were obtained from The Cancer Genome Atlas (TCGA), International Cancer Genome Consortium (ICGC) and Gene Enrichment Analysis (GSEA) website. Least absolute shrinkage and selection operator (LASSO) was used to screen out the candidate mRNAs and calculate the risk coefficient to establish the prognosis model. A high-risk group and low-risk group were separated for further study depending on their median risk score. The reliability of the prediction was evaluated in the validation cohort and the whole cohort. Results: A total of 548 differential mRNAs were identified from HCC samples (n = 374) and normal controls (n = 50), 45 of which were correlated with prognosis. A total of 373 samples met the screening criteria and there were randomly divided into the training cohort (n = 186) and the validation cohort (n = 187). In the training cohort, six metabolism-related mRNAs were used to construct a prognostic model with a LASSO regression model. Based on the risk model, the overall survival rate of the high-risk cohort was significantly lower than that of the low-risk cohort. The results of a time-ROC curve proved that the risk score (AUC = 0.849) had a higher prognostic value than the pathological grade, clinical stage, age or gender. Conclusion: The model constructed by the six metabolism-related mRNAs has a significant value for survival prediction and can be applied to guide the evaluation of HCC and the designation of clinical therapy.