• OPEN ACCESS

A Novel mRNA Signature Related to Immunity to Predict Survival and Immunotherapy Response in Hepatocellular Carcinoma

  • Chenhao Zhou1,2,# ,
  • Jialei Weng1,# ,
  • Yuan Gao2,3,#,
  • Chunxiao Liu2,
  • Xiaoqiang Zhu4,
  • Qiang Zhou1,
  • Chia-Wei Li2,
  • Jialei Sun1,
  • Manar Atyah1,
  • Yong Yi1,
  • Qinghai Ye1,
  • Yi Shi5,
  • Qiongzhu Dong2,6,7,
  • Yingbin Liu3,
  • Mien-Chie Hung2,8 and
  • Ning Ren1,6,* 
 Author information
Journal of Clinical and Translational Hepatology 2022;10(5):925-938

DOI: 10.14218/JCTH.2021.00283

Abstract

Background and Aims

Hepatocellular carcinoma (HCC) is the most common primary liver cancer and the incidence and mortality rates are increasing. Given the limited treatments of HCC and promising application of immunotherapy for cancer, we aimed to identify an immune-related prognostic signature that can predict overall survival (OS) rates and immunotherapy response in HCC.

Methods

The initial signature development was conducted using a training dataset from the Cancer Genome Atlas followed by independent internal and external validations from that resource and the Gene Expression Omnibus. A signature based nomogram was generated using multivariate Cox regression analysis. The associations of signature score with tumor immune phenotype and response to immunotherapy were analyzed using single-sample gene set enrichment analysis and tumor immune dysfunction and exclusion algorithm. A cohort from Zhongshan Hospital was employed to verify the predictive robustness of the signature regarding prognostic risk and immunotherapy response.

Results

The prognostic signature, IGSHCC, consisting of 22 immune-related genes, had independent prognostic ability, with training and validation cohorts. Also, IGSHCC stratified HCC patients with different outcomes in subgroups. The prognostic accuracy of IGSHCC was better than three reported prognostic signatures. The IGSHCC-based nomogram had high accuracy and significant clinical benefits in predicting 3- and 5-year OS. IGSHCC reflected distinct immunosuppressive phenotypes in low- and high-score groups. Patients with low IGSHCC scores were more likely than those with high scores to benefit from immunotherapy.

Conclusions

IGSHCC predicted HCC prognosis and response to immunotherapy, and contributed to individualized clinical management.

Keywords

Hepatocellular carcinoma, Gene signature, Immune microenvironment, Prognosis, Immunotherapy

Introduction

Liver cancer is the fourth leading cause of cancer deaths worldwide. Its incidence is increasing most rapidly in the USA, with 42,810 new cases estimated in 2020. Hepatocellular carcinoma (HCC), the most common form of liver cancer, is correlated with chronic infection with hepatitis B or C virus, aflatoxin-contaminated foodstuff consumption, and heavy alcohol intake.1,2 However, due to the universal late diagnosis, most HCC patients are unlikely to be amenable to potentially curative treatments.3 Some multi-kinase inhibitors are currently approved by the U.S. Food and Drug Administration for advanced HCCs. However, use of sorafenib or regorafenib has only improved the median overall survival (OS) duration by fewer than 3 months in advanced HCC patients, whose prognoses remain dismal.4,5 With an increasing number of HCC patients experiencing resistance to targeted therapy, efforts to develop novel therapy for HCC has drawn attention to the HCC microenvironment.6

Recent evidence has highlighted the importance of the tumor immune microenvironment for HCC outcome.7–9 Clinical trials of agents that target immune checkpoints such as programmed death 1 (PD-1, programmed death-ligand 1 (PD-L1), and cytotoxic T-lymphocyte-associate protein 4 (CTLA4) in HCC patients are ongoing, with an overall treatment response rate of less than 20%.10–12 The poor response to immunotherapy demonstrates that current patient-stratification methods are not sufficiently effective. The HCC immune microenvironment remains poorly understood and effective immune biomarkers for the evaluation of HCC prognosis are lacking. Identifying immune-related signatures for HCC patients that can serve to predicting the outcomes of and guide immunotherapy for HCC is needed.

We aimed to add to the knowledge of the association of immune-related genes with HCC prognosis and response to immunotherapy, which will set the groundwork for studies of the identified genes and provide a practical tool to optimize clinical HCC management. We studied 577 HCC cases from one cohort each in two public databases, identifying and validating immune-related gene set signature, IGSHCC as a robust and independent prognostic classifier for HCC with greater predictive accuracy than other prognostic models. We also established an IGSHCC-based nomogram to confirm the clinical applicability of our signature. We discovered a close association of IGSHCC score with tumor immunocyte infiltration status and response to immunotherapy and confirmed the utility of IGSHCC in prognostic risk-stratification and immunotherapy response prediction in another independent HCC cohort receiving anti-PD-1 therapy. Overall, the study finding support IGSHCC as a robust signature with potential clinical significance, given its value in individualized prognostic prediction and immunotherapy guidance for HCC.

Methods

Study populations and gene expression data preprocessing

The gene expression profiles and clinical information for HCC cases were collected from two public microarray datasets, The Cancer Genome Atlas (TCGA)13 and series GSE14520 in the Gene Expression Omnibus (GEO) from the Liver Cancer Institute (LCI) at Zhongshan Hospital, Fudan University.14 Only patients with available follow-up data, gene expression data, and clinical profiles were included. A total of 356 cases from TCGA and 221 cases from the GSE14520 were analyzed in our study.

Level 3 TCGA RNA sequencing data for HCC cases were downloaded and logarithmic-transformed. Gene expression data from GSE14520 were processed on the Affymetrix GeneChip Human Genome U133A 2.0 platform. The difference in the batch effect between the two datasets was adjusted by the ComBat method.15 To ensure consistency in gene names across the cohorts, expression profiles were collapsed from probe level to the corresponding gene symbols based on the annotation platform of each set.

The 356 patients in the TCGA dataset were randomly assigned to training (n=250) and internal validation (n=106) datasets. GSE14520 served as an external validation dataset (LCI validation set 1). Another independent cohort of HCC patient who underwent hepatectomy at Zhongshan Hospital in 2012 (n=56) was retrospectively evaluated as a polymerase chain reaction (PCR) array validation dataset (LCI validation set 2) to further test the reliability of IGSHCC. The clinical characteristics of the patients in the four datasets are summarized in Supplementary Table 1. This study was approved by the ethics committee of Zhongshan Hospital (approval number: B2021-143R), and written informed consent was obtained from all patients. The study was conducted following the ethical standards of the Declaration of Helsinki.

Immune-related gene set definition

A comprehensive list of immune-related genes was downloaded from the ImmPort database.16 The list consists of 1,534 genes that are involved in various immune pathways, including antimicrobial pathways, the T-cell receptor signaling pathway, interleukins, and the tumor necrosis factor family receptor signaling pathways. This list was used as a gene pool from which we screened out key genes to develop our immune-related prognostic signature. Among the genes in the list, 837 whose transcriptional expression levels could be obtained across all platforms, were selected for further analysis.

Development of the immune-related gene signature for HCC

LASSO, a method of regression analysis with high-dimensional factors, has been extended, and is used in Cox proportional hazards regression models for survival analysis.17,18 After application of univariate Cox regression analysis to the TCGA training dataset to identify candidate genes associated with prognosis for HCC (p<0.01) in the immune-related gene list, the LASSO Cox regression model was used with the glmnet package in the R computing language to screen the most powerful prognostic genes from the candidate genes. The best penalty parameter lambda was determined using 10-fold cross-validation. Using the formula below, a weighted mRNA level-based value, IGSHCC, was identified, and normalized to predict prognosis for HCC:

Prognostic score=i=1Nwixi
where N, w, and x represent the number of the selected genes, coefficient value, and gene expression level, respectively. The formula was used to calculate the prognostic scores of HCC patients. The information about and corresponding weight coefficients for 22 identified immune-related genes are listed in Supplementary Table 2. The cutoff value for the low- and high-IGSHCC groups was determined to be the median IGSHCC score in the training set and was used for further validation analysis.

Validation of the immune-related gene signature for HCC

The IGSHCC score was further in the internal and the external validation sets, and patients were assigned to the high- or low-IGSHCC group based on the cutoff value described above. Univariate and multivariate Cox regression analyses were performed to evaluate the independent predictive ability of IGSHCC, age, sex, α-fetoprotein (AFP) level, tumor grade, and TNM stage for OS in the HCC patients. TCGA cases were placed in several subgroups according to TP53 genotype, CTNNB1 genotype, AFP level, age, tumor grade, and TNM stage. The performance of IGSHCC was then evaluated in the subgroups. The OS difference in the high- and low-IGSHCC groups was tested for significance by the Kaplan-Meier method and log-rank test.

Existing HCC prognostic signatures for comparison with IGSHCC

Three other published prognostic signatures were retrospectively collected for comparison with IGSHCC (Supplementary Table 3).19–21 The gene numbers in these signatures ranged from four to 65. The score for each signature was calculated by the corresponding formula, and a receiver operating characteristic (ROC) curve analysis was used to compare the predictive efficiency of different signatures using the pROC package in the R language. The area under the ROC curve (AUC) was calculated for all signatures, and the prognostic signature with the highest AUC value was considered the most efficient one.

Construction and evaluation of the IGSHCC-based nomogram

In cancer care, nomograms are models commonly used to predict prognosis. Nomograms generate an individual numerical probability of a clinical event by integrating diverse variables, fulfilling our desire for biological and clinical integration of models for prognosis prediction and assisting our drive toward personalized cancer medicine.22 Thus, IGSHCC was integrated with age and TNM stage, independent prognostic indicators for HCC identified via multivariate Cox regression analysis of the entire TCGA dataset, to establish a prognostic nomogram for predicting the likelihood of 3- and 5-year OS using the rms package in the R language.

Subsequently, to evaluate the predictive performance and clinical practicability of the nomogram, calibration plots and decision curve analysis (DCA) were performed in the entire TCGA and LCI datasets. Calibration plots were helpful in assessing the consistency of the actual prognoses with the predicted outcomes using the nomogram. In addition, time-dependent ROC analysis was carried out to evaluate the OS-predictive accuracy of the nomogram.

Single-sample gene set enrichment analysis (ssGSEA)

ssGSEA was conducted as described previously.23 Twenty-four immunocyte types that could be discriminated by gene sets specifically overexpressed in each immunocyte using a deconvolution approach were included in this analysis. The immunocytes are involved in both innate immunity (natural killer cells, CD56dim natural killer cells, CD56bright natural killer cells, plasmacytoid dendritic cells [DCs], immature DCs, activated DCs, neutrophils, mast cells, eosinophils, and macrophages) and adaptive immunity (B cells, central memory T cells, effector memory T cells, activated CD8+ T cells, T follicular helper cells, γδ T cells, T helper [Th] cells, type 1 (Th1), Th2, and Th17, and regulatory T cells). The T-cell infiltration score was defined as the average of the standardized ssGSEA values for CD8+ T, central memory T, effector memory T, Th1, Th2, Th17, and regulatory T cells. The cytotoxic cell infiltration score was calculated according to the geometrical mean of perforin (PRF1) and granzyme A (GZMA) expression levels.24 The DC infiltration score and macrophage infiltration score were defined as the standardized ssGSEA values for DC subsets (total, plasmacytoid, immature, and activated) and macrophages, respectively.

Molecular characterization and immunotherapy response prediction

To describe the difference in molecular and signaling pathway characterization between the high- and low-IGSHCC groups, gene set variation analysis25 was applied to our dataset with the C7 (immunological) and H (hallmark) gene sets downloaded from the Molecular Signatures Database.26 For the C7 gene set, the significant pathway sets were identified according to p-values <0.001 and fold-change values of at least 1.85. For the H gene set, the significant pathways were identified according to t-values of at least 2. In addition, the tumor immune dysfunction and exclusion (TIDE) algorithm and subclass mapping were used to predict the response of cancer to immunotherapy, such as anti-PD-1 and anti-CTLA4 therapy.27,28

Real-time PCR array analysis

An RT2 Profiler Custom PCR array (Qiagen, Hilden, Germany) was conducted to simultaneously examine the mRNA expression levels of 22 immune-related genes in our signature. The array was designed and produced by BioTNT and performed following the manufacturer’s instructions. Briefly, total RNA was extracted from HCC samples using TRIzol reagent (Invitrogen, Waltham, MA, USA), and cDNAs were prepared using an iScript cDNA Synthesis Kit (Bio-Rad, Hercules, CA, USA). The PCR primers for the 22 immune-related genes are listed in Supplementary Table 4. Real-time PCR analysis was performed using RT2 SYBR Green qPCR Master Mix (Qiagen) with a ViiA 7 (384-well block) instrument (Thermo Fisher Scientific, Waltham, MA, USA). The relative expression levels of the target genes were calculated using the ΔΔCt method.

Histological and immunofluorescent analyses

HCC tissue sections were prepared and incubated with a rabbit anti-PD-L1 antibody (ab205921; Abcam, Cambridge, UK) and mouse anti-CD8 antibody (ab199016; Abcam). Next, the slides were incubated with goat antirabbit IgG antibody (Alexa Fluor 488; Abcam), goat antimouse IgG antibody (Alexa Fluor 647; Abcam), and 4′,6-diamidino-2-phenylindole (Invitrogen). Images of the sections were acquired using a confocal microscope with 40x objectives and analyzed using image processing software (Olympus Corp., Tokyo, Japan).

Statistical analysis

All statistical analyses were performed in the R language (version 3.6.0) and SPSS software (version 23.0; IBM Corp., Armonk, NY, USA). Continuous variables were reported as means ± standard deviation and were compared using Student t-test. Categorical variables were reported as frequencies and percentages (%) and compared using the Pearson chi-square test. A Cox proportional hazards model was used to conduct standard univariate and multivariate analyses. Prognostic differences between groups or subgroups were analyzed using the Kaplan-Meier method and log-rank tests. For all tests, two-tailed p-values <0.05 were considered significant.

Results

IGSHCC is a powerful prognostic signature with high predictive accuracy

The study design and workflow are shown in Figure 1. Univariate Cox regression analysis identified 78 genes from the immune-related gene list that were significantly associated with the OS of HCC patients (p<0.01) in the TCGA training set as the candidate genes for signature construction (Supplementary Table 5). Subsequently, by applying the LASSO Cox regression model to these genes, we constructed a prognostic signature consisting of 22 immune-related genes (Supplementary Fig. 1A, B). Most of the genes in the signature were cytokines, chemokines, and their receptors. To varying degrees, they played roles in the transforming growth factor-β signaling pathway, antigen processing and presentation, and adaptive immune response (Supplementary Fig. 1C). Detailed information on and corresponding coefficients for the 22 immune-related genes are presented in Supplementary Table 2 and Figure 1D. We calculated the IGSHCC score as the sum of the weighted expression levels for the 22 genes. To guarantee the consistency of our study findings, the cutoff value for high- and low-IGSHCC group was set as the median IGSHCC score in the TCGA training dataset and fixed as 1.291 during further validation processes (Supplementary Fig. 2).

Flow chart of the study design.
Fig. 1  Flow chart of the study design.

The TCGA HCC dataset was randomly divided into training (n=250) and validation (n=106) sets. GSE14520 (n=221) and PCR array (n=56) cohorts served as external validation datasets. Seventy-eight prognosis-related genes were identified in the TCGA training dataset by performing univariate Cox regression analysis, and IGSHCC was constructed using a LASSO Cox regression model. The risk-stratification ability of IGSHCC and its association with tumor immune phenotype and HCC response to immunotherapy were explored. Good predictive performance of the signature was confirmed in several analyses. HCC, hepatocellular carcinoma; LASSO, least absolute shrinkage and selection operator; TCGA, The Cancer Genome Atlas.

Kaplan-Meier survival curve analysis (Fig. 2A) demonstrated an significant prognostic difference in OS between the high- and low-IGSHCC groups in the TCGA training dataset [hazard ratio (HR), 3.87, 95% confidence interval (CI), 2.44–6.14]; p<0.0001). We also confirmed the robust predictive performance of IGSHCC in the internal validation dataset (Fig. 2B), the entire TCGA dataset (Fig. 2C), and the LCI validation dataset (Fig. 2D). Of note, in the LCI validation dataset, the high- and low-IGSHCC groups also had markedly different clinical outcomes (HR: 1.99 95% CI: 1.30–3.05; p=0.0013). Additionally, we compared the expression levels of the 22 immune-related genes in the high- and low-IGSHCC groups in the four datasets (Supplementary Fig. 2). The heatmaps in Supplementary Figure 2 show that the expression levels of genes such as BMP6, GLP1R, and GMFB were higher in the high-IGSHCC group than in the low-IGSHCC group. In contrast, tumors with low IGSHCC scores had much higher expression of genes like CCR8, TNF, and RETN than did those with high scores. Taken together, the results demonstrate that the signature we constructed was a powerful and effective tool for prognostic risk prediction for HCC patients.

Kaplan-Meier analysis of OS in HCC patients classified by IGS<sub>HCC</sub>.
Fig. 2  Kaplan-Meier analysis of OS in HCC patients classified by IGSHCC.

(A) TCGA training dataset, (B) TCGA validation dataset, (C) entire TCGA dataset, and (D) LCI validation dataset. The prognostic difference between the low- and high-IGSHCC groups was determined by log-rank tests. HCC, hepatocellular carcinoma; OS, overall survival; TCGA, The Cancer Genome Atlas.

IGSHCC score is an independent risk factor for HCC prognosis

To confirm the independent prognostic value of our signature, we performed univariate and multivariate Cox regression analyses of the TCGA training dataset (Fig. 3A, B). In addition to IGSHCC, we included several common prognostic factors, like age, sex, AFP level, tumor grade, and TNM stage, in the analysis. The univariate regression model demonstrated that IGSHCC was significantly associated with OS (HR: 4.11, 95% CI: 3.03–5.57; p<0.0001). Furthermore, in the multivariate regression model, IGSHCC remained a strong independent prognostic predictor of OS (HR: 5.51, 95% CI: 3.40–8.96; p<0.0001). The results were consistent with those obtained from the internal validation dataset (Fig. 3C, D), LCI validation set (Fig. 3E, F), and the entire TCGA dataset (Supplementary Fig. 3A, B). In the LCI validation cohort, IGSHCC was still considered a significant prognostic variable (univariate p=0.001 and multivariate p=0.035). The results revealed that the IGSHCC score was a powerful independent risk factor for HCC prognosis.

Univariate and multivariate Cox regression analyses of OS in the HCC patients.
Fig. 3  Univariate and multivariate Cox regression analyses of OS in the HCC patients.

(A and B) TCGA training data set, (C and D) TCGA validation data set, and (E and F) LCI validation dataset based on the IGSHCC score and clinicopathological risk factors, including age, sex, AFP level, tumor grade, TNM stage, and Barcelona Clinic Liver Cancer (BCLC) stage. The forest plots show the HRs for death (solid blue squares) and their 95% CIs (horizontal red lines) for each variable. P-values were calculated by Cox regression hazards analysis. HCC, hepatocellular carcinoma; OS, overall survival; TCGA, The Cancer Genome Atlas.

To determine whether IGSHCC retained its prognostic value for OS in HCC patient subgroups with different molecular or clinical characteristics, we placed the patients in the entire TCGA dataset into different subgroups according to TP53 genotype, CTNN1B genotype, AFP level, age, tumor grade, and TNM stage. IGSHCC succeeded in distinguishing HCC patients with distinct OS durations in all of the subpopulations (Supplementary Fig. 4A–L). Again, these results demonstrated that IGSHCC is a promising predictive signature with high clinical relevance that can be used to further assess the prognostic risk in HCC patients on the basis of existing stratification systems.

IGSHCC outperforms previously published HCC prognostic signatures

The number of prognostic signatures for HCC reported in the literature is increasing. We chose representative signatures to compare their effectiveness with that of IGSHCC. The selection criteria were that the signature was identified in more than 100 HCC patients and validated in an independent cohort. The details of the selected studies are provided in Supplementary Table 3.19–21 Because ROC curves are useful for organizing classifiers and visualizing their performance, we performed ROC analyses to evaluate the predictive accuracy of these selected signatures and IGSHCC regarding OS within the TCGA dataset. As shown in Figure 4A, we drew an ROC curve for each signature according to their sensitivity and specificity in survival prediction. In each curve, we drew a 45-degree line in gray as the negative reference. IGSHCC had a greater AUC (0.67, 95% CI: 0.62–0.72) than those of the other prognostic signatures, indicating that it provided the most accurate prognosis prediction. Although the difference in AUC between the IGSHCC and TNM stage was not significant (0.67 vs. 0.65), the combination of IGSHCC and TNM stage had greater predictive efficacy than did the TNM stage alone (0.72 vs. 0.65; p=0.0046) and age alone (0.72 vs. 0.51; p<0.0001). These results implied that IGSHCC outperforms other published prognostic signatures and, when integrated with the TNM stage, increased the sensitivity and precision of the prediction of OS of HCC patients.

Predictive performance of IGS<sub>HCC</sub> and development of the IGS<sub>HCC</sub>-based nomogram for HCC patients.
Fig. 4  Predictive performance of IGSHCC and development of the IGSHCC-based nomogram for HCC patients.

(A) ROC analyses of the predictive efficacy of TNM stage combined with IGSHCC, IGSHCC alone, TNM stage alone, age alone, and three published prognostic signatures for predicting OS in HCC patients in the entire TCGA dataset. The difference in AUC among the combined models and other factors was tested using the Delong method. (B) Nomogram developed for predicting the probability of 3- and 5-year OS in HCC patients. The point of each variable was determined by drawing a vertical line upward to the point scale. The probability of 3- and 5-year OS was determined by drawing a vertical line from the total point axis downward to the probability axis. (C and D) Calibration curves for the nomogram used to evaluate the consistency of the predicted and actual 5-year OS durations in the TCGA and LCI validation datasets. The 45-degree dashed line represents the ideal consistency. (E and F) ROC curves for the nomogram used to assess its predictive efficacy regarding 3- and 5-year OS. (G and H) DCA of the nomogram (red line) and IGSHCC (blue line) in the entire TCGA and LCI validation datasets. The black line at the bottom and the gray line represent the assumptions that none and all HCC patients survived at 5 years, respectively. The X-axis represents the threshold probability of death, and the Y-axis represents the corresponding net benefit of the model. AUC, area under the curve; HCC, hepatocellular carcinoma; OS, overall survival; ROC, receiver operating characteristic; TCGA, The Cancer Genome Atlas.

An IGSHCC-based nomogram is useful in clinical HCC management

Because IGSHCC proved to be a powerful prognostic signature model for HCC, we further attempted to provide clinicians with a quantitative method of predicting the 3- and 5-year OS rates in HCC patients. Therefore, we generated a nomogram by integrating IGSHCC with age and TNM stage to leverage the complementary value of these independent prognostic factors (Fig. 4B). Calibration curves were consistent with actual observations of OS duration and nomogram-based prediction of 5-year OS likelihood in the TCGA (Fig. 4C) and LCI validation (Fig. 4D) datasets. Similarly, regarding 3-year OS, the bias-corrected lines in calibration curves were still close to the ideal 45-degree line for both datasets (Supplementary Fig. 5A, B), indicating good predictive performance of our nomogram for HCC prognosis. In addition, we plotted ROC curves to evaluate the predictive accuracy of the nomogram. In the entire TCGA dataset, the AUCs for 3- and 5-year OS were both 0.72 (Fig. 4E). We further verified the predictive efficacy of IGSHCC in the LCI validation dataset, with AUCs of 0.71 and 0.69 for 3- and 5-year OS, respectively (Fig. 4F).

DCA can determine whether application of the nomogram to predicting cancer outcomes increases net clinical benefits when compared with the treat-none or treat-all scheme. The results of DCA in the CGA and LCI validation datasets showed that the nomogram provided HCC patients with greater net clinical benefits in predicting 5-year OS likelihood than did IGSHCC alone (Fig. 4G, H). The nomogram also exhibited greater prediction threshold probabilities than IGSHCC. Consistently, the DCA of 3-year OS also demonstrated increased net clinical benefits and threshold probabilities of the nomogram (Supplementary Fig. 5C, D). The results confirmed the clinical practicability of IGSHCC in that, when integrated into a nomogram, it can serve as a quantitative method for predicting OS rates in HCC patients.

IGSHCC can distinguish HCC patients with contrasting tumor immune microenvironments

An increasing number of patients with cancer have benefited from immune checkpoint blockade.7 Considering that we constructed the IGSHCC based on 22 immune-related genes, we wanted to determine whether the signature can provide some clues about tumor immune microenvironments, which would contribute to treatment responder identification and immunotherapy guidance among HCC patients. To that end, we applied ssGSEA to tumor samples from the CGA to deconvolve the relative abundance of 24 immunocyte types based on bulk tumor gene expression profiles. We visualized the abundance of immunocytes in each tumor sample using a heatmap and then categorized tumor immune infiltration into three subgroups (low, moderate, and high) using unsupervised hierarchical clustering (Fig. 5A). We also assessed the associations of immune cell infiltration with other clinical features, such as the TP53 mutation status. HCCs with high IGSHCC scores tended to exhibit relatively low tumor immune cell infiltration statuses and vice versa. This demonstrated that high-IGSHCC patients may have had “cold” tumors and were unlikely to respond to immunotherapy and that low IGSHCC patients may have benefited from immunotherapy. These results demonstrated that IGSHCC reveals, to some extent, the immune phenotypes for HCCs.

IGS<sub>HCC</sub> distinguished HCC patients with contrasting immune microenvironments.
Fig. 5  IGSHCC distinguished HCC patients with contrasting immune microenvironments.

(A) The tumor immune landscape of HCC patients in the TCGA dataset. The relative abundance of 24 immunocytes in each tumor sample was calculated using ssGSEA, and three distinct subgroups (low, moderate, and high infiltration) were defined based on unsupervised hierarchical clustering. IGSHCC, TP53, and CTNNB1 mutation status, AFP level, microvascular invasion (MVI), histological grade, TNM stage, and survival status are annotated in the lower panel. (B-E) Violin plots of the Z-scores for the T-cell infiltration score, cytotoxic cell infiltration score, DC infiltration score, and macrophage infiltration score in the low- and high-IGSHCC groups. The scores were compared using the Student t-test. AFP, α-fetoprotein; HCC, hepatocellular carcinoma; MVI, microvascular invasion; TCGA, The Cancer Genome Atlas.

Furthermore, considering the concomitant infiltration of tumor promoting and suppressing immune cells in the tumor microenvironment, we compared the cytotoxic function of T cells in the low- and high-IGSHCC groups to see if the prognostic signature was correlated with the cytotoxic potential of immunocytes. As expected, the T-cell infiltration score (p<0.0001), cytotoxic cell infiltration score (p=0.0001), and DC infiltration score (p<0.0001) were increased in the low-IGSHCC group (Fig. 5B–D), patients in which tended to have high immunocyte infiltration statuses. In contrast, the macrophage infiltration score (p=0.0015) was higher in HCC patients with high IGSHCC scores than in those with low IGSHCC scores (Fig. 5E). The data demonstrated that HCC patients with low IGSHCC scores may have immune-activated phenotypes, characterized by increased immune cell infiltration and high cytotoxic potential, whereas those with high IGSHCC scores suffer from relatively low levels of immune activation.

Low- and high-IGSHCC groups have different sensitivities to immunotherapy

Because HCC immunotherapy response is closely related to immune activation in the local microenvironment, we then performed gene set variation analysis to explore the association between immune functional condition and IGSHCC score. As depicted in the heatmap in Figure 6A, distinct functional pathways were enriched in the low- and high-IGSHCC groups. Interleukin-6 signaling, interferon-α response, interferon-γ response, and T-cell activation were among the most significantly enriched biological processes in the low-IGSHCC group. Next, we characterized the low- and high-IGSHCC groups regarding the hallmark gene sets. Several tumorigenesis pathways, such as unfolded protein response, G2/M checkpoint, the phosphoinositide 3-kinase/Akt/mammalian target of rapamycin signaling pathway, the transforming growth factor-β signaling pathway, and DNA repair pathways, were enriched in the high-IGSHCC group (Fig. 6B). These results again showed that HCC patients in the low- and high-IGSHCC groups have distinct immune conditions in the tumor microenvironment.

Low- and high-IGS<sub>HCC</sub> groups of HCC patients had different sensitivity to immunotherapy.
Fig. 6  Low- and high-IGSHCC groups of HCC patients had different sensitivity to immunotherapy.

(A) Heatmap of the enriched immune-related and inflammatory pathways in the low- and high-IGSHCC groups in the TCGA dataset. (B) Gene set variation analysis of differentially enriched hallmark gene sets in the low- and high-IGSHCC groups. The vertical dotted lines represent t-value=2. A t-value for the gene set variation analysis score greater than 2 represented significantly enriched pathways in the high-IGSHCC group. (C and D) The association between the IGSHCC score and the possibility of HCC response to immunotherapy calculated using the TIDE algorithm in the entire TCGA and LCI validation datasets. (E) Subclass mapping analysis of putative immunotherapeutic response in the low- and high-IGSHCC groups. HCC, hepatocellular carcinoma; TCGA, The Cancer Genome Atlas; TIDE, tumor immune dysfunction and exclusion.

To determine whether our prognostic signature can predict HCC patient responses to immunotherapy, we used the TIDE algorithm to calculate the probability of tumor response to immunotherapy in the TCGA dataset and LCI validation set. The four-fold contingency tables in Figure 6C, D showed that low-IGSHCC patients were more likely than high-IGSHCC patients to have immunotherapy-sensitive disease (TCGA dataset, p<0.001; LCI validation set, p=0.032). Moreover, we applied subclass mapping to compare the expression profiles for the HCCs in these two groups with a published dataset for 47 patients with melanoma that responded to anti-PD-1 and anti-CTLA4 immunotherapy.29 This demonstrated that compared with high-IGSHCC patients, low-IGSHCC patients were more likely to have responses to anti-PD-1 therapy (Bonferroni corrected p=0.008; Fig. 6E). The findings proved that our prognostic signature can effectively predict sensitivity of HCC to immunotherapy and thus contribute to identification of suitable candidates for anti-PD-1 therapy. IGSHCC score can serve as a useful reference for determining the immunotherapeutic strategy for HCC.

IGSHCC successfully stratifies HCC patients in an independent PCR array cohort

To validate the robust power of IGSHCC in prognostic risk-stratification and immunotherapy response prediction, we collected and analyzed data on an independent cohort of HCC patients in our institution. The detailed clinical characteristics of the HCC patients in this PCR array validation dataset (LCI validation set 2) are shown in Supplementary Table 1. Based on the formula described in Materials and Methods, the IGSHCC score for each individual was calculated, and patients were assigned to the high-IGSHCC (n=33) or low-IGSHCC (n=23) groups according to the cutoff value (Fig. 7A). Kaplan-Meier analysis demonstrated that IGSHCC retained the ability to predict prognosis for HCC patients in this cohort (HR: 2.49, 95% CI, 1.19–5.22; p=0.012; Fig. 7B).

Independent clinical validation of the prognostic risk-stratification and immunotherapy response prediction ability of IGS<sub>HCC</sub>.
Fig. 7  Independent clinical validation of the prognostic risk-stratification and immunotherapy response prediction ability of IGSHCC.

(A) The IGSHCC score and expression of 22 immune-related genes in the HCC patients in the PCR array cohort (n=56). Top panel, the low- and high-IGSHCC groups defined based on the cutoff value for TCGA training dataset. Middle panel, the survival status and duration in the HCC cases in the low- and high-IGSHCC groups. Bottom panel, heatmap of the expression of the 22 immune-related genes in each tumor. (B) Kaplan-Meier analysis of OS in the HCC patients classified according to IGSHCC in the PCR array validation dataset. (C) Representative immunofluorescent images of nuclei (blue), PD-L1 (green), and CD8 (red) in tumor samples obtained from HCC patients who had responses (left panel) or no response (right panel) to anti-PD-1 therapy. (D) The IGSHCC scores in patients with HCC that responded well and did not respond to anti-PD-1 therapy. HCC, hepatocellular carcinoma; OS, overall survival; PCR, polymerase chain reaction; PD-1, programmed death 1; PD-L1, programmed death-ligand 1; PRF1, perforin; TCGA, The Cancer Genome Atlas.

To confirm the potential value of IGSHCC in identifying HCC patients who need immunotherapy, we examined the infiltration of CD8+ T cells and expression of PD-L1 in the tumors in patients who received nivolumab (anti-PD-1 antibody) in our PCR array validation dataset using immunofluorescence. Compared with the nonresponders, nivolumab responders tended to have greater CD8+ T-cell infiltration and lower expression of PD-L1 in tumors (Fig. 7C). Importantly, we found that HCC patients who responded well to treatment with nivolumab had lower IGSHCC scores than did nonresponders (Fig. 7D), which greatly supported our previous conclusion that the signature can help identify suitable HCC candidates for immunotherapy. Collectively, the observations highlighted the clinical significance of IGSHCC in predicting survival and immunotherapy response in HCC patients.

Discussion

The heterogeneity of HCC leads to greatly varied clinical outcomes among patients, and effective prognostic predictors are lacking. Furthermore, although the achievements of immunotherapy have revolutionized cancer treatment and provided significant survival benefits to patients,7,30 accumulating evidence from preclinical and clinical studies demonstrates that traditional immune-related markers such as PD-L1 are neither consistent nor reliable predictors of immunotherapy outcomes.31 In view of that, we designed this study to develop a signature based on immune-related genes that is associated with the tumor immune microenvironment and can provide an individualized prediction of both prognosis and immunotherapy sensitivity for HCC patients. A novel 22-gene signature named IGSHCC was developed, and proved to be a reliable tool to predict HCC survival and immunotherapy response, which was successfully validated in multiple public datasets, as well as an independent clinical cohort of HCC patients in our own institution. In addition, the IGSHCC-based nomogram had satisfying accuracy and sensitivity in predicting the likelihood of 3- and 5-year OS of HCC and provided added clinical benefits to HCC patients. The study results will serve to improve prognostic prediction and personalized therapy in HCC patients by measuring the expression of genes in the IGSHCC panel in tumors.

Based on large public datasets, researchers have developed several prognostic signatures for HCC,19–21 but they have yet to be incorporated into clinical practice because of limitations such as the overfitting risks arising from screening of a large gene pool. In addition, Cox proportional hazards regression analysis is unsuitable for high-dimensional microarray data in which the number of covariates is close to or greater than the number of observations.32 In this study, we used the LASSO Cox regression method, which can select the most valuable variables with robust prognostic value and low correlations to prevent overfitting.17 As expected, IGSHCC had better predictive performance than did other published prognostic models, and there was no reduction in stability and reliability of our signature when applied across independent datasets. In addition, there is no significant gene overlap between our signature and other published models, even if the same construction method was used. The main reason is that the distribution ratios we and other researchers used to randomly divided dataset into training set and internal verification set were different. Therefore, a unified modeling standard and procedure is warranted to promote the consistent use of predictive signatures in clinical practice.

Most of the genes in our prognostic signature were cytokines, chemokines, and their receptors that play key roles in immune and inflammatory response. Several of the genes, including STC2, TNF, CCR8, and BMP6, are associated with human cancers. For example, increased STC2 expression is correlated with aggressive disease and shortened OS in renal cell carcinoma patients.32 Also, BMP6 has been shown to inhibit mast cell recruitment in melanoma, which indicates the role of this gene in modulation of the tumor immune microenvironment.33 Some genes included in the signature may be therapeutic targets in HCC patients, but their biological functions have yet to be fully explored in HCC cases, which is well deserving of further investigation.

Some researchers have proposed that stratifying cancer patients based on the characteristics and quality of tumor immune infiltration is a critical next step in predicting the responsiveness of their disease to immunotherapy.34 In this study, the immune landscape of HCCs depicted by ssGSEA sheds light on the negative correlation between immune infiltration status and IGSHCC score. The complexity and diversity of tumor immune components drove us to explore the association of the IGSHCC score with T-cell activation and cytotoxic function. The decreased T-cell infiltration score, cytotoxic cell infiltration score, and DC infiltration score, and increased macrophage infiltration in patients with high IGSHCC scores demonstrated that their tumors are more likely than those in patients with low scores to be infiltrated-excluded and poorly immunogenic, with a lack of effective antigen-presenting cells and adaptive antitumor cytotoxic T cells. Moreover, several immune-related pathways were enriched in tumors in the low-IGSHCC group, whereas some tumorigenesis pathways were enriched in those in the high-IGSHCC group, further supporting the distinct immune biological processes in low- and high-IGSHCC tumors. In addition, given the close association between IGSHCC score and tumor immune phenotype, we performed an analysis using the TIDE algorithm and subclass mapping that demonstrated that HCC patients with low IGSHCC scores were more likely than those with high scores to benefit from anti-PD-1 immunotherapy.

Most important, we validated the robust prognostic risk-stratification ability and immunotherapy response prediction of IGSHCC in an independent cohort of HCC patients from our institution. External verification of the signature in a real-world cohort strongly demonstrated that prognosis risk and immunotherapy response in HCC patients can be readily evaluated by measuring the expression of 22 immune-related genes in tumor samples using PCR. Thus, our study, for the first time, shows the feasibility of translating IGSHCC into a clinical assay and set the groundwork for future prospective cohort studies to assess the overall translational potential of this signature.

The study had some limitations. First, the retrospective study design and single-center validation require future prospective validation of the results in multicenter cohorts. Second, the patients’ clinical information was incomplete, which impeded comprehensive analyses and decreased statistical reliability and validity. Moreover, the function of some genes in the signature in HCC cases is not clearly understood, so basic experimental studies are warranted to reveal the biological mechanisms underlying our prognostic signature. Finally, although we verified the predictive value of IGSHCC in HCC patients undergoing surgical resection, for major liver cancer patients who cannot receive surgery, whether the expression level of these genes in the tumor can be obtained through minimally invasive methods like needle biopsy and liquid biopsy needs further research to explore whether the signature can be applied to the majority of HCC cases.

In conclusion, IGSHCC is the first immune-related prognostic signature for HCC patients and is a robust and independent biomarker for predicting the risk of OS, even in specific clinical subgroups of patients. Furthermore, this immune-based signature provides a broad view of tumor immune phenotypes, resulting in improved identification of HCC patients who are likely to benefit from immunotherapy. The clinical feasibility of IGSHCC and the IGSHCC-based nomogram should be studied further prospectively in multicenter cohorts.

Supporting information

Supplementary Fig. 1

Development of the immune-related gene signature for HCC.

(A) Partial likelihood deviance of covariates in 78 OS-related candidate genes. (B) LASSO coefficient profiles of the 22 immune-related genes. The coefficient value was chosen using 10-fold cross-validation. (C) The clustered immune function of the 22 immune-related genes included in IGSHCC. (D) The coefficient of each gene in IGSHCC ranked from high to low value. The yellow bars represent the positive coefficients, whereas the blue ones represent the negative coefficients.

(DOCX)

Supplementary Fig. 2

Distribution of IGSHCC scores in HCC patients.

(A) TCGA training data set, (B) TCGA validation set, (C) entire TCGA data set, and (D) LCI validation data set. In each figure, the top panel represents the low- and high-IGSHCC groups defined based on the cutoff value for TCGA training data set, the middle panel represents the survival statuses and durations in the HCC patients in the low- and high-IGSHCC groups, and the bottom panel is a heatmap of the expression of the 22 immune-related genes in IGSHCC in each tumor.

(DOCX)

Supplementary Fig. 3

Cox regression analysis of OS in the HCC patients in the entire TCGA dataset based on the IGSHCC score and clinicopathological risk factors, including age, sex, AFP level, tumor grade, and TNM stage.

(A) Univariate and (B) multivariate. The forest plots show the HRs for death (solid blue squares) and their 95% CIs (horizontal red lines) for each variable.

JCTH-D-21-00283-s003.docx

(DOCX)

Supplementary Fig. 4

Kaplan-Meier analyses of OS in the HCC patients in subgroups with different molecular and clinical characteristics.

Patients in the entire TCGA dataset were stratified into subgroups including (A and B) TP53 genotype, (C and D) CTNNB1 genotype, (E and F) AFP level, (G and H) age, (I and J) tumor grade, and (K and L) TNM stage. The survival difference between the low- and high-IGSHCC patients in these subgroups was analyzed using the log-rank test.

JCTH-D-21-00283-s004.docx

(DOCX)

Supplementary Fig. 5

Evaluation of the predictive accuracy and net benefit of IGSHCC-based nomogram for 3-year OS.

(A and B) Calibration plots for the IGSHCC-based nomogram used to evaluate the consistency of predicted and actual 3-year OS durations in the entire TCGA dataset and LCI validation data set. The 45-degree line represents ideal performance of the nomogram. (C and D) DCA of the nomogram (red line) and IGSHCC (blue line) in the TCGA and LCI validation dataset. The black line at the bottom and gray line represents the assumptions that none and all HCC patients survived at 3 years, respectively. The x-axis represents the threshold probability of death, and the y-axis represents the corresponding net benefit of the model.

JCTH-D-21-00283-s005.docx

(DOCX)

Supplementary Table 1

Clinical characteristics of the HCC patients in the four data sets.

(DOCX)

Supplementary Table 2

Details of the 22 immune-related genes in IGSHCC.

(DOCX)

Supplementary Table 3

Summary of the three reported HCC prognostic signatures.

(DOCX)

Supplementary Table 4

Primer sequences for the 22 immune-related genes in IGSHCC used in the PCR array.

(DOCX)

Supplementary Table 5

Seventy-eight candidate immune-related genes identified by univariate Cox regression analysis.

(DOCX)

Abbreviations

AFP: 

α-fetoprotein

AUC: 

area under the curve

BCLC: 

Barcelona Clinic Liver Cancer

CTLA4: 

cytotoxic T-lymphocyte-associate protein 4

DCA: 

decision curve analysis

DC: 

dendritic cell

GEO: 

Gene Expression Omnibus

GZMA: 

granzyme A

HCC: 

hepatocellular carcinoma

LASSO: 

least absolute shrinkage and selection operator

MVI: 

microvascular invasion

OS: 

overall survival

PCR: 

polymerase chain reaction

PD-1: 

programmed death 1

PD-L1: 

programmed death-ligand 1

PRF1: 

perforin

ROC: 

receiver operating characteristic

ssGSEA: 

single-sample gene set enrichment analysis

TCGA: 

The Cancer Genome Atlas

TIDE: 

tumor immune dysfunction and exclusion

TNM: 

tumor-nodes-metastasis

Declarations

Acknowledgement

We thank the China Scholarship Council for supporting Dr. Chenhao Zhou and Dr. Yuan Gao as a visiting PhD candidate at The University of Texas MD Anderson Cancer Center. We also thank Donald R. Norwood at the Research Medical Library at the University of Texas MD Anderson Cancer Center for editing this manuscript.

Data sharing statement

All data in this study are available from the Cancer Genome Atlas website (http://cancergenome.nih.gov/) and series GSE14520 in the GEO from the Liver Cancer Institute (LCI) at Zhongshan Hospital, Fudan University (https://www.ncbi.nlm.nih.gov/geo/). Other data analyzed in this study are included in this article and its supplementary materials.

Funding

The work was supported by the National Natural Science Foundation of China (82103521, 82073208); the Shanghai Sailing Program (21YF1407500); the China Postdoctoral Science Foundation (2021M690674); the Special Foundation for Science and Technology Basic Research Program (2019FY101103); the Shanghai Shen Kang Hospital Development Center new frontier technology joint project (SHDC12021109).

Conflict of interest

QD has been an editorial board member of Journal of Clinical and Translational Hepatology since 2021. The other authors have no conflicts of interest related to this publication.

Authors’ contributions

Study concept, design, and supervision (NR), analysis and interpretation of data (CZ, JW), drafting of the manuscript (CZ, JW, YG), acquisition of data (CL, XZ, QZ, JS), preparation of figures and tables (CWL, MA), critical revision of the manuscript for important intellectual content (MCH, YS, QD, YL), provision of patient tissue samples (YY, QY).

References

  1. Siegel RL, Miller KD, Jemal A. Cancer statistics, 2020. CA Cancer J Clin 2020;70(1):7-30 View Article PubMed/NCBI
  2. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin 2018;68(6):394-424 View Article PubMed/NCBI
  3. Llovet JM, Zucman-Rossi J, Pikarsky E, Sangro B, Schwartz M, Sherman M, et al. Hepatocellular carcinoma. Nat Rev Dis Primers 2016;2:16018 View Article PubMed/NCBI
  4. Llovet JM, Ricci S, Mazzaferro V, Hilgard P, Gane E, Blanc JF, et al. Sorafenib in advanced hepatocellular carcinoma. N Engl J Med 2008;359(4):378-390 View Article PubMed/NCBI
  5. Bruix J, Qin S, Merle P, Granito A, Huang YH, Bodoky G, et al. Regorafenib for patients with hepatocellular carcinoma who progressed on sorafenib treatment (RESORCE): a randomised, double-blind, placebo-controlled, phase 3 trial. Lancet 2017;389(10064):56-66 View Article PubMed/NCBI
  6. Zhu YJ, Zheng B, Wang HY, Chen L. New knowledge of the mechanisms of sorafenib resistance in liver cancer. Acta Pharmacol Sin 2017;38(5):614-622 View Article PubMed/NCBI
  7. Prieto J, Melero I, Sangro B. Immunological landscape and immunotherapy of hepatocellular carcinoma. Nat Rev Gastroenterol Hepatol 2015;12(12):681-700 View Article PubMed/NCBI
  8. Li H, Li CW, Li X, Ding Q, Guo L, Liu S, et al. MET Inhibitors Promote Liver Tumor Evasion of the Immune Response by Stabilizing PDL1. Gastroenterology 2019;156(6):1849-1861.e1813 View Article PubMed/NCBI
  9. Chan LC, Li CW, Xia W, Hsu JM, Lee HH, Cha JH, et al. IL-6/JAK1 pathway drives PD-L1 Y112 phosphorylation to promote cancer immune evasion. J Clin Invest 2019;129(8):3324-3338 View Article PubMed/NCBI
  10. Zhu AX, Finn RS, Edeline J, Cattan S, Ogasawara S, Palmer D, et al. Pembrolizumab in patients with advanced hepatocellular carcinoma previously treated with sorafenib (KEYNOTE-224): a non-randomised, open-label phase 2 trial. Lancet Oncol 2018;19(7):940-952 View Article PubMed/NCBI
  11. Duffy AG, Ulahannan SV, Makorova-Rusher O, Rahma O, Wedemeyer H, Pratt D, et al. Tremelimumab in combination with ablation in patients with advanced hepatocellular carcinoma. J Hepatol 2017;66(3):545-551 View Article PubMed/NCBI
  12. El-Khoueiry AB, Sangro B, Yau T, Crocenzi TS, Kudo M, Hsu C, et al. Nivolumab in patients with advanced hepatocellular carcinoma (CheckMate 040): an open-label, non-comparative, phase 1/2 dose escalation and expansion trial. Lancet 2017;389(10088):2492-2502 View Article PubMed/NCBI
  13. The Cancer Genome Atlas Research Network. Comprehensive and Integrative Genomic Characterization of Hepatocellular Carcinoma. Cell 2017;169(7):1327-1341.e1323 View Article PubMed/NCBI
  14. Roessler S, Jia HL, Budhu A, Forgues M, Ye QH, Lee JS, et al. A unique metastasis gene signature enables prediction of tumor relapse in early-stage hepatocellular carcinoma patients. Cancer Res 2010;70(24):10202-10212 View Article PubMed/NCBI
  15. Johnson WE, Li C, Rabinovic A. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics 2007;8(1):118-127 View Article PubMed/NCBI
  16. Bhattacharya S, Andorf S, Gomes L, Dunn P, Schaefer H, Pontius J, et al. ImmPort: disseminating data to the public for the future of immunology. Immunol Res 2014;58(2-3):234-239 View Article PubMed/NCBI
  17. Zhu X, Tian X, Sun T, Yu C, Cao Y, Yan T, et al. GeneExpressScore Signature: a robust prognostic and predictive classifier in gastric cancer. Mol Oncol 2018;12(11):1871-1883 View Article PubMed/NCBI
  18. Tibshirani R. The lasso method for variable selection in the Cox model. Stat Med 1997;16(4):385-395 View Article PubMed/NCBI
  19. Kim SM, Leem SH, Chu IS, Park YY, Kim SC, Kim SB, et al. Sixty-five gene-based risk score classifier predicts overall survival in hepatocellular carcinoma. Hepatology 2012;55(5):1443-1452 View Article PubMed/NCBI
  20. Qiao GJ, Chen L, Wu JC, Li ZR. Identification of an eight-gene signature for survival prediction for patients with hepatocellular carcinoma based on integrated bioinformatics analysis. PeerJ 2019;7:e6548 View Article PubMed/NCBI
  21. Zheng Y, Liu Y, Zhao S, Zheng Z, Shen C, An L, et al. Large-scale analysis reveals a novel risk score to predict overall survival in hepatocellular carcinoma. Cancer Manag Res 2018;10:6079-6096 View Article PubMed/NCBI
  22. Balachandran VP, Gonen M, Smith JJ, DeMatteo RP. Nomograms in oncology: more than meets the eye. Lancet Oncol 2015;16(4):e173-180 View Article PubMed/NCBI
  23. Zhang L, Zhao Y, Dai Y, Cheng JN, Gong Z, Feng Y, et al. Immune Landscape of Colorectal Cancer Tumor Microenvironment from Different Primary Tumor Location. Front Immunol 2018;9:1578 View Article PubMed/NCBI
  24. Rooney MS, Shukla SA, Wu CJ, Getz G, Hacohen N. Molecular and genetic properties of tumors associated with local immune cytolytic activity. Cell 2015;160(1-2):48-61 View Article PubMed/NCBI
  25. Hänzelmann S, Castelo R, Guinney J. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinformatics 2013;14:7 View Article PubMed/NCBI
  26. Liberzon A, Birger C, Thorvaldsdóttir H, Ghandi M, Mesirov JP, Tamayo P. The Molecular Signatures Database (MSigDB) hallmark gene set collection. Cell Syst 2015;1(6):417-425 View Article PubMed/NCBI
  27. Jiang P, Gu S, Pan D, Fu J, Sahu A, Hu X, et al. Signatures of T cell dysfunction and exclusion predict cancer immunotherapy response. Nat Med 2018;24(10):1550-1558 View Article PubMed/NCBI
  28. Hoshida Y, Brunet JP, Tamayo P, Golub TR, Mesirov JP. Subclass mapping: identifying common subtypes in independent disease data sets. PLoS One 2007;2(11):e1195 View Article PubMed/NCBI
  29. Roh W, Chen PL, Reuben A, Spencer CN, Prieto PA, Miller JP, et al. Integrated molecular analysis of tumor biopsies on sequential CTLA-4 and PD-1 blockade reveals markers of response and resistance. Sci Transl Med 2017;9(379):eaah3560 View Article PubMed/NCBI
  30. Chen L, Han X. Anti-PD-1/PD-L1 therapy of human cancer: past, present, and future. J Clin Invest 2015;125(9):3384-3391 View Article PubMed/NCBI
  31. Lee HH, Wang YN, Xia W, Chen CH, Rau KM, Ye L, et al. Removal of N-Linked Glycosylation Enhances PD-L1 Detection and Predicts Anti-PD-1/PD-L1 Therapeutic Efficacy. Cancer Cell 2019;36(2):168-178.e164 View Article PubMed/NCBI
  32. Meyer HA, Tölle A, Jung M, Fritzsche FR, Haendler B, Kristiansen I, et al. Identification of stanniocalcin 2 as prognostic marker in renal cell carcinoma. Eur Urol 2009;55(3):669-678 View Article PubMed/NCBI
  33. Stieglitz D, Lamm S, Braig S, Feuerer L, Kuphal S, Dietrich P, et al. BMP6-induced modulation of the tumor micro-milieu. Oncogene 2019;38(5):609-621 View Article PubMed/NCBI
  34. Binnewies M, Roberts EW, Kersten K, Chan V, Fearon DF, Merad M, et al. Understanding the tumor immune microenvironment (TIME) for effective therapy. Nat Med 2018;24(5):541-550 View Article PubMed/NCBI
  • Journal of Clinical and Translational Hepatology
  • pISSN 2225-0719
  • eISSN 2310-8819
  • Copyright © 2022 JCTH. All Rights Reserved.
  • Published by Xia & He Publishing Inc.
  • Address: 14090 Southwest Freeway, Suite 300, Sugar Land, Texas 77478, USA
  • Email: service@xiahepublishing.com