The Heterogeneity of Immune Cell Infiltration Landscape and Its Immunotherapeutic Implications in Hepatocellular Carcinoma

Comprehensive analysis of clinical and biological value of ING family genes in liver cancer

BACKGROUND

Liver cancer (LIHC) is a malignant tumor that occurs in the liver and has a high mortality in cancer. The ING family genes were identified as tumor suppressor genes. Dysregulated expression of these genes can lead to cell cycle arrest, senescence and/or apoptosis. ING family genes are promising targets for anticancer therapy. However, their role in LIHC is still not well understood.

AIM

To have a better understanding of the important roles of ING family members in LIHC.

METHODS

A series of bioinformatics approaches (including gene expression analysis, genetic alteration analysis, survival analysis, immune infiltration analysis, prediction of upstream microRNAs (miRNAs) and long noncoding RNAs (lncRNAs) of ING1, and ING1-related gene functional enrichment analysis) was applied to study the expression profile, clinical relationship, prognostic significance and immune infiltration of ING in LIHC. The relationship between ING family genes expression and tumor associated immune checkpoints was investigated in LIHC. The molecular mechanism of ING1 mediated hepatocarcinogenesis was preliminarily discussed.

RESULTS

mRNA/protein expression of different ING family genes in LIHC was analyzed in different databases, showing that ING family genes were highly expressed in LIHC. In 47 samples from 366 LIHC patients, the ING family genes were altered at a rate of 13%. By comprehensively analyzing the expression, clinical pathological parameters and prognostic value of ING family genes, ING1/5 was identified. ING1/5 was related to poor prognosis of LIHC, suggesting that they may play key roles in LIHC tumorigenesis and progression. One of the target miRNAs of ING1 was identified as hsa-miR-214-3p. Two upstream lncRNAs of hsa-miR-214-3p, U91328.1, and HCG17, were identified. At the same time, we found that the expression of ING family genes was correlated with immune cell infiltration and immune checkpoint genes.

CONCLUSION

This study lays a foundation for further research on the potential mechanism and clinical value of ING family genes in the treatment and prognosis of LIHC.

Comprehensive analysis of clinical and biological value of ING family genes in liver cancer

BACKGROUND

Liver cancer (LIHC) is a malignant tumor that occurs in the liver and has a high mortality in cancer. The ING family genes were identified as tumor suppressor genes. Dysregulated expression of these genes can lead to cell cycle arrest, senescence and/or apoptosis. ING family genes are promising targets for anticancer therapy. However, their role in LIHC is still not well understood.

AIM

To have a better understanding of the important roles of ING family members in LIHC.

METHODS

A series of bioinformatics approaches (including gene expression analysis, genetic alteration analysis, survival analysis, immune infiltration analysis, prediction of upstream microRNAs (miRNAs) and long noncoding RNAs (lncRNAs) of ING1, and ING1-related gene functional enrichment analysis) was applied to study the expression profile, clinical relationship, prognostic significance and immune infiltration of ING in LIHC. The relationship between ING family genes expression and tumor associated immune checkpoints was investigated in LIHC. The molecular mechanism of ING1 mediated hepatocarcinogenesis was preliminarily discussed.

RESULTS

mRNA/protein expression of different ING family genes in LIHC was analyzed in different databases, showing that ING family genes were highly expressed in LIHC. In 47 samples from 366 LIHC patients, the ING family genes were altered at a rate of 13%. By comprehensively analyzing the expression, clinical pathological parameters and prognostic value of ING family genes, ING1/5 was identified. ING1/5 was related to poor prognosis of LIHC, suggesting that they may play key roles in LIHC tumorigenesis and progression. One of the target miRNAs of ING1 was identified as hsa-miR-214-3p. Two upstream lncRNAs of hsa-miR-214-3p, U91328.1, and HCG17, were identified. At the same time, we found that the expression of ING family genes was correlated with immune cell infiltration and immune checkpoint genes.

CONCLUSION

This study lays a foundation for further research on the potential mechanism and clinical value of ING family genes in the treatment and prognosis of LIHC.

High LGALS3 expression induced by HCP5/hsa-miR-27b-3p correlates with poor prognosis and tumor immune infiltration in hepatocellular carcinoma

BACKGROUND

Hepatocellular carcinoma (HCC) is widely recognized for its unfavorable prognosis. Increasing evidence has revealed that LGALS3 has an essential function in initiating and developing several malignancies in humans. Nevertheless, thorough analysis of the expression profile, clinical prognosis, pathway prediction, and immune infiltration of LGALS3 has not been fully explored in HCC.

METHODS

In this study, an initial pan-cancer analysis was conducted to investigate the expression and prognosis of LGALS3. Following a comprehensive analysis, which included expression analysis and correlation analysis, noncoding RNAs that contribute to the overexpression of LGALS3 were subsequently identified. This identification was further validated using HCC clinical tissue samples. TIMER2 and GEPIA2 were employed to examine the correlation between LGALS3 and HCP5 with immunological checkpoints, cell chemotaxis, and immune infiltration in HCC. The R program was applied to analyze the expression distribution of immune score in in HCC patients with high and low LGALS3 expression. The expression profiles of immune checkpoints were also analyzed. Use R to perform GSVA analysis in order to explore potential signaling pathways.

RESULTS

First, we conducted pan-cancer analysis for LGALS3 expression level through an in-depth analysis of public databases and found that HCC has a high LGALS3 gene and protein expression level, which were then verified in clinical HCC specimens. Meanwhile, high LGALS3 gene expression is related to malignant progression and poor prognosis of HCC. Univariate and multivariate analyses confirmed that LGALS3 could serve as an independent prognostic marker for HCC. Next, by combining comprehensive analysis and validation on HCC clinical tissue samples, we hypothesize that the HCP5/hsa-miR-27b-3p axis could serve as the most promising LGALS3 regulation mechanism in HCC. KEGG and GO analyses highlighted that the LGALS3-related genes were involved in tumor immunity. Furthermore, we detected a significant positive association between LGALS3 and HCP5 with immunological checkpoints, cell chemotaxis, and immune infiltration. In addition, high LGALS3 expression groups had significantly higher immune cell scores and immune checkpoint expression levels. Finally, GSVA analysis was performed to predict potential signaling pathways linked to LGALS3 and HCP5 in immune evasion and metabolic reprogramming of HCC.

CONCLUSIONS

Our findings indicated that the upregulation of LGALS3 via the HCP5/hsa-miR-27b-3p axis is associated with unfavorable prognosis and increased tumor immune infiltration in HCC.

DCLK1 is Overexpressed and Associated with Immune Cell Infiltration in Hepatocellular Carcinoma

Recent research has shown that Doublecortin-like kinase 1 (DCLK1) is overexpressed in different types of cancer. It has recently been described as a cancer stem cells (CSCs) marker, is associated with carcinogenesis, and positively correlates with infiltration of multiple immune cell types in some cancers. However, studies focused on assessing DCLK1 expression in HCC are limited, and the role of DCLK1 in HCC tumor immunity remains to be determined. In this study, we used a modified model of the resistant hepatocyte (MRHM) to evaluate DCLK1 expression in HCC. Furthermore, DCLK1 expression in HCC was analyzed using TIMER 2.0, UALCAN, GEPIA, GEO, and HPA web-based tools. Correlations between DCLK1 expression and clinicopathological factors in patients were analyzed using the UALCAN web-based tool. Finally, correlations between DCLK1 and immune infiltrates were investigated using the TIMER 2.0 and TISIDB web-based tools. The results showed that DCLK1 is significantly overexpressed during progression of the HCC carcinogenic process in the MRHM. DCLK1 is overexpressed in HCC according to multiple publics web-based tools, and its overexpression is associated with cancer stage. Furthermore, DCLK1 expression was correlated with infiltration levels of multiple immune cells, immunomodulatory factors, immunoinhibitors, MHC molecules, chemokines, receptors, and immune cell-specific markers. These results suggest that DCLK1 is a potential prognostic biomarker that determines cancer progression and correlates with immune cell infiltration in HCC.

SAMD13 serves as a useful prognostic biomarker for hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is the most common form of liver cancer and the 5-year relative overall survival (OS) rate is less than 20%. Since there are no specific symptoms, most patients with HCC are diagnosed in an advanced stage with poor prognosis. Therefore, identifying novel prognostic biomarkers to improve the survival of patients with HCC is urgently needed. In the present study, we attempted to identify SAMD13 (Sterile Alpha Motif Domain-Containing Protein 13) as a novel biomarker associated with the prognosis of HCC using various bioinformatics tools. SAMD13 was found to be highly expressed pan-cancer; however, the SAMD13 expression was significantly correlated with the worst prognosis in HCC. Clinicopathological analysis revealed that SAMD13 upregulation was significantly associated with advanced HCC stage and high-grade tumor type. Simultaneously, high SAMD13 expression resulted in association with various immune markers in the immune cell subsets by TIMER databases and efficacy of immunotherapy. Methylation analysis showed SAMD13 was remarkably associated with prognosis. Furthermore, a six-hub gene signature associated with poor prognosis was correlated with the cell cycle, transcription, and epigenetic regulation and this analysis may support the connection between SAMD13 expression and drug-resistance. Our study illustrated the characteristics of SAMD13 role in patients with HCC using various bioinformatics tools and highlights its potential role as a therapeutic target and promising biomarker for prognosis in HCC.

Pan-cancer analysis identifies RNF43 as a prognostic, therapeutic and immunological biomarker

BACKGROUND

RING finger protein 43 (RNF43), an E3 ubiquitin ligase, is a homologous gene mutated in several cancers. However, the pan-cancer panoramic picture of RNF43 and its predictive value for tumor immune phenotypes and immunotherapeutic efficacy are still largely unclear. Our study aims to clarify the functions of RNF43 in predicting the prognosis, immune signature, and immunotherapeutic efficacy in pan-cancer.

METHODS

By using RNA-seq, mutation, and clinical data from the TCGA database, the expression levels and prognostic significance of RNF43 in pan-cancer were analyzed. The genetic alteration characteristics of RNF43 were displayed by the cBioPortal database. Gene Set Enrichment Analysis (GSEA) was performed to investigate the potential biological functions and signaling pathways modulated by RNF43 in cancers. The relationship of RNF43 expression with immune cell infiltration, and immune modulators expression was interpreted by the ESTIMATE algorithm, CIBERSORT algorithm, and TISIDB database. The correlations between RNF43, microsatellite instability (MSI), and tumor mutation burden (TMB) were also investigated. Furthermore, the predictive value of RNF43 for immunotherapeutic efficacy and drug sensitivity was further illustrated. Besides, immunohistochemistry (IHC) was employed to validate the expression of the RNF43 in different cancer types by our clinical cohorts, including patients with lung cancer, sarcoma, breast cancer, and kidney renal clear cell carcinoma.

RESULTS

The results demonstrated that RNF43 was abnormally expressed in multiple cancers, and RNF43 is a critical prognosis-related factor in several cancers. RNF43 was frequently mutated in several cancers with a high frequency of 4%, and truncating mutation was the most frequent RNF43 mutation type. RNF43 expression was linked to the abundance of several immune cell types, including CD8+ T cells, B cells, and macrophages within the tumor immune microenvironment. Furthermore, RNF43 expression was significantly correlated with the efficacy of anti-PD-1/PD-L1 treatment, and it could predict the sensitivity of various anti-cancer drugs. Finally, IHC explored and validated the different expression levels of RNF43 in different cancers by our clinical samples.

CONCLUSION

Our results first present the expression pattern and the mutation signature of RNF43, highlighting that RNF43 is an important prognostic biomarker in pan-cancer. Furthermore, RNF43 seems to be a critical modulator in the tumor immune microenvironment and can function as a promising biomarker for predicting the immunotherapeutic efficacy of anti-PD-1/PD-L1 treatment, and drug sensitivity in cancer treatment.

Integrated single cell and bulk sequencing analysis identifies tumor reactive CXCR6+ CD8 T cells as a predictor of immune infiltration and immunotherapy outcomes in hepatocellular carcinoma

Background

Various immune cell types in the tumor microenvironment (TME) of hepatocellular carcinoma (HCC) have been identified as important parameters associated with prognosis and responsiveness to immunotherapy. However, how various factors influence immune cell infiltration remains incompletely understood. Hence, we investigated the single cell multi-omics landscape of immune infiltration in HCC, particularly key gene and cell subsets that influence immune infiltration, thus potentially linking the immunotherapy response and immune cell infiltration.

Methods

We grouped patients with HCC according to immune cell infiltration scores calculated by single sample gene set enrichment analysis (ssGSEA). Differential expression analysis, functional enrichment, clinical trait association, gene mutation analysis, tumor immune dysfunction and exclusion (TIDE) and prognostic model construction were used to investigate the immune infiltration landscape through multi-omics. Stepwise regression was further used to identify key genes regulating immune infiltration. Single cell analysis was performed to explore expression patterns of candidate genes and investigate associated cellular populations. Correlation analysis, ROC analysis, Immunotherapy cohorts were used to explore and confirm the role of key gene and cellular population in predicting immune infiltration state and immunotherapy response. Immunohistochemistry and multiplexed fluorescence staining were used to further validated our results.

Results

Patients with HCC were clustered into high and low immune infiltration groups. Mutations of CTNNB1 and TTN were significantly associated with immune infiltration and altered enrichment of cell populations in the TME. TIDE analysis demonstrated that T cell dysfunction and the T cell exclusion score were elevated in the high and low infiltration groups, respectively. Six risk genes and five risk immune cell types were identified and used to construct risk scores and a nomogram model. CXCR6 and LTA, identified by stepwise regression, were highly associated with immune infiltration. Single cell analysis revealed that LTA was expressed primarily in tumor infiltrating T lymphocytes and partial B lymphocytes, whereas CXCR6 was enriched predominantly in T and NK cells. Notably, CXCR6+ CD8 T cells were characterized as tumor enriched cells that may be potential predictors of high immune infiltration and the immune-checkpoint blockade response, and may serve as therapeutic targets.

Conclusion

We constructed a comprehensive single cell and multi-omics landscape of immune infiltration in HCC, and delineated key genes and cellular populations regulating immune infiltration and immunotherapy response, thus providing insights into the mechanisms of immune infiltration and future therapeutic control.

From Reductionistic Approach to Systems Immunology Approach for the Understanding of Tumor Microenvironment

The tumor microenvironment (TME) is a complex and dynamic ecosystem that includes a variety of immune cells mutually interacting with tumor cells, structural/stromal cells, and each other. The immune cells in the TME can have dual functions as pro-tumorigenic and anti-tumorigenic. To understand such paradoxical functions, the reductionistic approach classifies the immune cells into pro- and anti-tumor cells and suggests the therapeutic blockade of the pro-tumor and induction of the anti-tumor immune cells. This strategy has proven to be partially effective in prolonging patients' survival only in a fraction of patients without offering a cancer cure. Recent advances in multi-omics allow taking systems immunology approach. This essay discusses how a systems immunology approach could revolutionize our understanding of the TME by suggesting that internetwork interactions of the immune cell types create distinct collective functions independent of the function of each cellular constituent. Such collective function can be understood by the discovery of the immunological patterns in the TME and may be modulated as a therapeutic means for immunotherapy of cancer.

Chromosome 11q13 amplification correlates with poor response and prognosis to PD-1 blockade in unresectable hepatocellular carcinoma

Background & aims

Little is known about molecular biomarkers that predict the response and prognosis in unresectable hepatocellular carcinoma (HCC) treated with programmed death (PD)-1 inhibitors.

Methods

A total of 62 HCC patients who underwent next-generation sequencing were retrospectively included in our department for this study. Patients with unresectable disease were subjected to systemic therapy. PD-1 inhibitors intervention (PD-1Ab) group and nonPD-1Ab group included 20 and 13 patients, respectively. Primary resistance was defined as initial on-treatment progression or progression with an initial stable disease of less than 6 months.

Results

Chromosome 11q13 amplification (Amp11q13) was the most common copy number variation in our cohort. Fifteen (24.2%) patients harbored Amp11q13 in our dataset. Patients with Amp11q13 showed higher level of Des-γ-carboxy-prothrombin (DCP), tumor number and were more prone to be combined with portal vein tumor thrombosis (PVTT). In the PD-1Ab group, the proportion of progressive disease (PD) in patients with Amp11q13 was significantly higher than that in patients with nonAmp11q13 (100% vs 33.3%, P=0.03). In the nonPD-1Ab group, the proportion of PD in patients with Amp11q13 and nonAmp11q13 had no significant difference (0% vs 11.1%, P>0.99). In the PD-1Ab group, the median progression-free survival (PFS) was 1.5 months in Amp11q13 patients vs 16.2 months in non-Amp11q13 patients (HR, 0.05; 95% CI 0.01-0.45; P = 0.0003). No significant difference was observed in the nonPD-1Ab group. Notably, we found that hyperprogressive disease (HPD) might be associated with Amp11q13. The increased density of Foxp3+ Treg cells in HCC patients with Amp11q13 might be one of potential mechanisms.

Conclusion

HCC patients with Amp11q13 are less likely to benefit from PD-1 blockade therapies. These findings may help guide the use of immunotherapy for HCC in routine clinical practice.

The current understanding of the immune landscape relative to radiotherapy across tumor types

Radiotherapy is part of the standard of care treatment for a great majority of cancer patients. As a result of radiation, both tumor cells and the environment around them are affected directly by radiation, which mainly primes but also might limit the immune response. Multiple immune factors play a role in cancer progression and response to radiotherapy, including the immune tumor microenvironment and systemic immunity referred to as the immune landscape. A heterogeneous tumor microenvironment and the varying patient characteristics complicate the dynamic relationship between radiotherapy and this immune landscape. In this review, we will present the current overview of the immunological landscape in relation to radiotherapy in order to provide insight and encourage research to further improve cancer treatment. An investigation into the impact of radiation therapy on the immune landscape showed in several cancers a common pattern of immunological responses after radiation. Radiation leads to an upsurge in infiltrating T lymphocytes and the expression of programmed death ligand 1 (PD-L1) which can hint at a benefit for the patient when combined with immunotherapy. In spite of this, lymphopenia in the tumor microenvironment of 'cold' tumors or caused by radiation is considered to be an important obstacle to the patient's survival. In several cancers, a rise in the immunosuppressive populations is seen after radiation, mainly pro-tumoral M2 macrophages and myeloid-derived suppressor cells (MDSCs). As a final point, we will highlight how the radiation parameters themselves can influence the immune system and, therefore, be exploited to the advantage of the patient.

Ficolin-2: A potential immune-related therapeutic target with low expression in liver cancer

Objective

This study aimed to investigate the role of ficolin-2 (FCN2) in the development and course of hepatocellular carcinoma (HCC) and to contribute to the evolution of innovative HCC therapeutics.

Methods

Oncomine, GEPIA (Gene Expression Profiling Interactive Analysis), TISIDB (Tumor Immune System Interactions and Drug Bank database), UALCAN (University of Alabama at Birmingham Cancer data analysis portal), UCSC (University of California, Santa Cruz), R package, the Kaplan–Meier technique, Cox regression analysis, LinkedOmics, Pearson’s correlation, and a nomogram were used to investigate the prognostic value of FCN2 in HCC. Co-expressed genes were screened. A protein–protein interaction network was created using the STRING database. Finally, immunohistochemistry was performed to establish the expression of FCN2 in HCC tissues. A pan-cancer study centered on HCC-related molecular analysis was also conducted to look for a link between FCN2 and immune infiltration, immune modulators, and chemokine receptors.

Results

In HCC tissues, the expression of FCN2 was observed to be lower than that in normal tissues. This was connected to the HCC marker alpha-fetoprotein, showing that FCN2 is involved in the development and progression of cancer. FCN2 may act through Staphylococcus aureus infection, lectins, and other pathways. Furthermore, at the immune level, the expression of FCN2 in HCC was associated with some immune cell infiltration, immunomodulators, and chemokine receptors.

Conclusion

FCN2 may be an immune checkpoint inhibitor for HCC, creating a breakthrough in the treatment of HCC.

Development and Verification of a Combined Immune- and Metabolism-Related Prognostic Signature for Hepatocellular Carcinoma

Immune escape and metabolic reprogramming are becoming important characteristics of tumor biology, which play critical roles in tumor initiation and progression. However, the integrative analysis of immune and metabolic characteristics for the tumor microenvironment in hepatocellular carcinoma (HCC) remains unclear. Herein, by univariate and least absolute shrinkage and selection operator (LASSO) Cox regression analyses, a prognostic signature associated with tumor microenvironment was established based on five immune- and metabolism-related genes (IMRGs), which was fully verified and evaluated in both internal and external cohorts. The C-index was superior to previously published HCC signatures, indicating the robustness and reliability of IMRGs prognostic signature. A nomogram was built based on IMRGs prognostic signature and various clinical parameters, such as age and T stage. The AUCs of nomogram at 1-, 3-, and 5-year (AUC = 0.829, 0.749, 0.749) were slightly better than that of IMRGs signature (AUC = 0.809, 0.734, 0.711). The relationship of risk score (RS) with immune checkpoint expressions, immunophenoscore (IPS), as well as microsatellite instability (MSI) together accurately predicted the treatment efficacy. Collectively, the IMRGs signature might have the potential to better predict prognostic risk, evaluate immunotherapy efficacy, and help personalize immunotherapy for HCC patients.

KCNN4 is a Potential Biomarker for Predicting Cancer Prognosis and an Essential Molecule that Remodels Various Components in the Tumor Microenvironment: A Pan-Cancer Study

Objectives: Potassium Calcium-Activated Channel Subfamily N Member 4 (KCNN4) is a member of the KCNN family. Studies have revealed that KCNN4 is implicated in various physiological processes as well as promotes the malignant phenotypes of cancer cells. However, little is known about its associations with survival outcomes across varying cancer types.

Methods: Herein, we systematically explored the prognostic value of KCNN4 in the pan-cancer dataset retrieved from multiple databases. Next, we performed correlation analysis of KCNN4 expression with tumor mutational burden (TMB) and microsatellite instability (MSI), and immune checkpoint genes (ICGs) to assess its potential as a predictor of immunotherapy efficacy. Afterwards, patients were divided into increased-risk group and decreased-risk group based on the contrasting survival outcomes in various cancer types. Furthermore, the underlying mechanisms of the distinctive effects were analyzed using ESTIMATE, CIBERSORT algorithms, and Gene Set Enrichment Analysis (GSEA) analysis.

Results: KCNN4 expression levels were aberrant in transcriptomic and proteomic levels between cancer and normal control tissues in pan-cancer datasets, further survival analysis elucidated that KCNN4 expression was correlated to multiple survival data, and clinical annotations. Besides, KCNN4 expression was correlated to TMB and MSI levels in 14 types and 12 types of pan-cancers, respectively. Meanwhile, different types of cancer have specific tumor-infiltrating immune cell (TICs) profiles.

Conclusions: Our results revealed that KCNN4 could be an essential biomarker for remodeling components in the tumor microenvironment (TME), and a robust indicator for predicting prognosis as well as immunotherapy response in pan-cancer patients.

High Expression of EZH2 Mediated by ncRNAs Correlates with Poor Prognosis and Tumor Immune Infiltration of Hepatocellular Carcinoma

Background: Hepatocellular carcinoma (HCC) is the predominant form of liver cancer and is accompanied by a complex regulatory network. Increasing evidence suggests that an abnormal gene expression of EZH2 is associated with HCC progression. However, the molecular mechanism by which non-coding RNAs (ncRNAs) regulate EZH2 remains elusive. Methods: The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) data were used to perform differential expression analysis and prognostic analysis. We used the Encyclopedia of RNA Interactomes (ENCORI) database to predict candidate miRNAs and lncRNAs that may bind to EZH2. Subsequently, the comprehensive analysis (including expression analysis, correlation analysis, and survival analysis) identified ncRNAs that contribute to EZH2 overexpression. Results: EZH2 was found to be upregulated in the majority of tumor types and associated with a poor prognosis. Hsa-miR-101-3p was identified as a target miRNA of EZH2. Additionally, SNHG6 and MALAT1 were identified as upstream lncRNAs of hsa-miR-101-3p. Meanwhile, correlation analysis revealed that EZH2 expression was significantly associated with the infiltration of several immune cell types in HCC. Conclusion: SNHG6 or MALAT1/hsa-miR-101-3p/EZH2 axis were identified as potential regulatory pathways in the progression of HCC.