- Research article
- Open Access
- Open Peer Review
Tissue-infiltrating lymphocytes signature predicts survival in patients with early/intermediate stage hepatocellular carcinoma
BMC Medicinevolume 17, Article number: 106 (2019)
Intratumoral immune infiltrates have manifested a robust prognostic signature in patients with hepatocellular carcinoma (HCC). We hypothesized that a novel tissue-related immune signature (TRIS) could improve the prediction of postoperative survival for patients diagnosed with early/intermediate HCC.
Twenty-eight immune features were immunohistochemically examined on 352 HCC specimens. The LASSO Cox regression model was used to construct a five-feature-based TRIS. The univariate and multivariate Cox analyses were performed. Based on independent predictors, the immune-clinical prognostic index (ICPI) was established. Performance assessment was measured with C-index and compared with seven traditional staging systems. The independent validation cohort (n = 393) was included to validate the model.
By using the LASSO method, the TRIS were constructed on the basis of five immune features, CD3intratumoral (T), CD27T, CD68peritumoral (P), CD103T, and PD1T. Multivariate Cox analysis showed that the TRIS was an independent prognostic predictor. In the training cohort, γ-glutamyl transferase, tumor diameter, tumor differentiation, and TRIS were incorporated into the ICPI. The ICPI presented satisfactory discrimination ability, with C-index values of 0.691 and 0.686 in the training and validation cohorts, respectively. Compared with seven conventional staging systems (C-index, training cohort, 0.548–0.597; validation cohort, 0.519–0.610), the ICPI exhibited better performance for early/intermediate-stage HCCs. Further, the patients were categorized into three subgroups with X-tile software, and the stratified ICPI presented a superior corrected Akaike information criterion and homogeneity in both cohorts.
Our ICPI was a useful and reliable prognostic tool which may offer good individualized prediction capability for HCC patients with early/intermediate stage.
Hepatocellular carcinoma (HCC) is a leading cause of cancer-related death worldwide, accounting for more than 700,000 deaths per year . The prognosis of HCC not only depends on the tumor burden of patients but also on their underlying liver functional reserve  and tumor heterogeneity . Therefore, in order to promote better prognosis and reduce the global burden of this disease, it is crucial to identify new pathological and biological predictors and improve prediction of postsurgical survival for HCC.
While multiple immune components are involved in cancer initiation and progression , the existence of various immune components have been identified in the liver, including Kupffer cells , dendritic cells (DCs) , natural killer (NK) cells, naïve and memory lymphocytes , B cells , T regulatory cells (Tregs) , T follicular helper (Tfh) cells , CD8+ T cells , and CD4+ T cells . Studies have suggested that immune infiltrates are of clinical significance in various types of cancer [13,14,15]. For example, the immune score, which is based on the number of lymphocytic populations in the tumor core and the invasive margin of tumor, was able to indicate outcomes in patients with early-stage colorectal cancer  and has been defined as a new component in the classification of colorectal cancer. In our previous studies, intratumoral neutrophils , margin-infiltrating CD20+ B cells by Shi et al. , and intratumoral balance of regulatory and cytotoxic T cells by Gao et al.  were found to be associated with long-term survival of HCC patients. It was reported that about 25% of HCC samples were found expressing markers of an inflammatory response, characterized with high expression levels of programmed cell death protein 1 (PD1) and CD274, and markers of cytolytic activity . Li et al. have also reported that the expression of PD1 in HCC was able to promote tumor growth independent of adaptive immunity . However, the prognostic potential of infiltrating immune cells in survival prediction for patients with HCC has not yet been assessed comprehensively.
In this study, we used the least absolute shrinkage and selection operator (LASSO) regression model based on 28 immunological features to establish the tissue-related immune signature (TRIS). Then, independent clinicopathological predictors and TRIS were integrated into a novel immune-clinical prognostic index (ICPI). Moreover, we investigated the performance of the ICPI, compared the ICPI with 7 conventional staging systems, and validated the ICPI model in the validation cohort.
Data collection was conducted from all HCC patients who underwent liver resection between April 2005 and September 2008 at the Department of Liver Surgery, Zhongshan Hospital, Shanghai, China. The inclusion criteria were without previous anticancer therapy, absence of any other types of malignancies, complete resection of liver tumors, histopathological confirmation of HCC, and classified as stage 0 or A, or B according to the Barcelona Clinic Liver Cancer (BCLC) staging system. In this study, the patients with BCLC B stage were carefully selected: resectable HCC, adequate liver remnant size after liver resection, no radiological evidence of vascular invasion or extrahepatic metastasis, and liver function status of Child A-B. Patients with hilar or extrahepatic cholangiocarcinoma, tumors of uncertain origin, metastatic liver tumors, combined hepatocellular-cholangiocarcinoma, perioperative mortality, or tumors beyond BCLC stage B were excluded from the study. Recruited patients were divided into 2 cohorts: the training and validation cohorts. The training cohort included patients who received surgery from April 2005 to December 2006. The validation cohort consisted of patients that received surgery from January 2007 to September 2008. The flowchart of patient recruitment and group assignment is presented in Additional file 3: Figure S1.
The study was censored on December 1, 2011. A standardized follow-up protocol was adopted for all patients . The patients attended follow-up visits with computed tomography or abdominal magnetic resonance imaging scans carried out once every 6 months for the first 2 years. Abdominal ultrasound, liver function tests, and serum alpha fetoprotein (AFP) level examinations were performed once every 3 months. The endpoints of the study were overall survival (OS) and recurrence-free survival (RFS). OS was defined as the interval between the date of surgery and the date of patient death or the last follow-up visit. RFS was defined as the time from the date of tumor resection to the date of diagnosis of recurrence/metastasis, or the last follow-up visit.
Tissue microarray (TMA) construction
All HCC tumor specimens were examined independently by two reviewers who were blinded to the clinical characteristics or treatment outcomes of patients. The TMA construction was conducted as previously described . To ensure homogeneity and reproducibility, two representative areas with infiltrating lymphocytes were carefully selected by qualified pathologists in H&E-stained slides of the blocks from the tumor center and peritumor tissues. Then, duplicate cylinders (1 mm in diameter) from two different areas were punched, arrayed, and re-embedded in a recipient block. All final slides were dipped in paraffin for preservation and stored at 4 °C before immunohistochemistry (IHC) assays.
In this study, immune biomarkers were selected for IHC staining based on their close involvement in patient survival and tumor recurrence (CD3, CD8, CD4, CD57, and CD68) [17, 18, 23,24,25,26], local immune response (CD20, CD27, and CD45RO) [8, 27], tumor growth (PD1 and CD66b) [21, 28], and antitumor function (CD14) . CXCR5 was selected because of its participation in the pathogenesis of primary biliary cirrhosis . While the prognostic values of CD45RA  and CD103  in liver cancer still remain unclear, they were also selected due to their presence in tumors. The details regarding IHC, antibodies, and staining conditions are included in the Additional file 1: Supplementary methods.
Evaluation of immunohistochemical staining
To evaluate the tissue-infiltrating immune cells, the tissue sections were screened at low power (× 100), and the 5 most representative and independent fields were selected using the Leica Qwin Plus v3 software to ensure representativeness and homogeneity. The respective areas of the nontumoral and intratumoral regions were then photographed at × 200 magnification. Identical settings were used for each photograph. High-resolution spot images (1360 × 1024) were obtained and stored under TIFF format. All the consecutive images were analyzed using a computer-automated method (Image-pro plus 6.0, Media Cybernetics Inc.) as described elsewhere . The spot and captured spot (× 200) with image software are presented in Additional file 3: Figure S2. The numbers of positively stained cells were recorded, and the mean value was used for statistical analysis. The 5 representative spots showed a good level of homogeneity of stained cell numbers in tumor or non-tumor regions.
Summary statistics were obtained using established methods and presented as percentages or median values. Pearson’s χ2 test or Fisher’s exact test was employed to compare categorical variables. Wilcoxon rank sum test or Student’s t test was used to evaluate continuous variables. The OS was estimated by the Kaplan-Meier method and compared with the log-rank test. Cluster version 3.0 (Michiel de Hoon, Tokyo, Japan) was performed for the hierarchical clustering of multi-immune features . The estimation of the relative fractions of immune cells from tissue expression profiles of HCC was conducted using CIBERSORT . The details regarding CIBERSORT and construction of immune network are included in Additional file 1: Supplementary methods.
LASSO is a broadly used method for regression with high-dimensional predictors . We applied the LASSO Cox analysis to identify significant prognostic immune features and constructed a multi-immune feature (TRIS score) on the basis of OS. The “glment” package was used to do the LASSO Cox analysis. By using the univariate and multivariate Cox proportional hazards regression in the training dataset, we integrated independent prognostic factors into the ICPI model. The nomogram and calibration plots were constructed as previously described . We compared the ICPI model with American Joint Committee on Cancer (AJCC) 7th edition, AJCC 8th edition, Cancer of the Liver Italian Program (CLIP), Barcelona Clínic Liver Cancer (BCLC), Okuda, Japan Integrated Staging (JIS) and Liver Cancer Study Group of Japan (LCSGJ) staging systems based on receiver operating characteristic (ROC) curves. The P value for the c-indices in the 2 models was computed using a bootstrapping method . The rcorrp.cens package in Hmisc was used.
X-tile software was used to generate the optimum cutoff point for continuous variables according to the highest χ2 value defined by the Kaplan-Meier survival analysis and the log-rank test . Based on the ICPI score, HCC patients were categorized into 3 subgroups with X-tile software version 3.6.1 (Yale University School of Medicine, New Haven, CT, USA). To evaluate the performance of stratified ICPI and other staging systems, the corrected Akaike information criterion (AIC) was chosen to present how the staging system correlated with patient survival. And then, Wald’s χ2 test was used to evaluate homogeneity in the same stage within each system, indicating differences in survival among patients . Statistical analyses were performed with R software version 3.1.0 (R Foundation for Statistical Computing, Vienna, Austria). Statistical significance was set at 0.05.
Clinicopathological characteristics of patients
The demographic and clinicopathological data of the 745 HCC patients recruited in this study are shown in Table 1. The patients were divided into the training cohort (n = 352) and the validation cohort (n = 393). The percentages of patients at BCLC stages 0, A, and B in the training cohort were 11.9%, 76.4%, and 11.7%, respectively. Except for albumin and bilirubin, no significant differences were observed between the training and validation cohorts in any other patient, tumor, or operation-related covariates. Despite significant differences existed in the albumin and bilirubin levels between the two cohorts, these two indicators were still within the normal range and would not influence the liver function of HCC patients.
After a median follow-up of 52.2 months (range, 3.0 to 79.3) for the entire study population, 54.8% of patients (408/745) had developed tumor recurrence, and 38.3% (285/745) had died. The 1-, 3-, and 5-year OS rates were 88.9%, 69.7%, and 56.3%, respectively, and the 1-, 3-, and 5-year RFS rates were 73.4%, 54.0%, and 36.5%, respectively.
Immune characteristics of HCC tissues
To investigate the cellular composition of the immune infiltrates in liver cancer, we initially built the CIBERSORT-inferred relative fractions of the different immune cell types with publicly available data (TCGA and 7 GEO datasets) . Among the 8 datasets, the percentage of macrophages was the highest, followed by CD4+ T cells, mast cells, and CD8+ T cells (Fig. 1a). Student’s t test revealed that the percentages of plasma cell, monocyte, CD8+ T cell, and neutrophil contents were decreased in intratumoral tissues, while the percentages of Tfh cells, Tregs, NK cells, and DCs were increased in TCGA and GSE14520 datasets (Fig. 1b). Further, we investigated the coordination of immune cell fractions in TCGA dataset. The correlation analysis was visualized using the unsupervised hierarchical clustering of a correlation matrix of immune cell analysis . Figure 1c shows 2 clusters characterized by immune cells of an exhausted immune response (neutrophilsintratumoral (T), eosinophilsT, and Tregperitumor (P) cells) and an adaptive T cell response (TfhT and TfhP), respectively.
To confirm these results, we evaluated the density of 28 immune features in 2 regions of liver cancer using tissue microarrays: tumor and peritumor. The expression levels of 28 immune features are presented in Additional file 3: Figure S3. Based on the correlation analysis, three major clusters were identified. Two clusters were characterized with functional coordination of T lymphocytes in the intratumoral tissues (CD3T, CD4T, CD27T, and CD45ROT) or the adjacent normal tissues (CD3P, CD45RAP, and CD45ROP). The last cluster was characterized with the functional coordination of peri- and intratumoral CXCR5+ cells (a biomarker of Tfh) , which is consistent with the results of TCGA dataset (TfhT and TfhP) (Fig. 1d). In summary, these data indicated a high degree of functional coordination of specific types of immune cells.
Construction of the TRIS
To select prognostic immune features, we performed the LASSO Cox regression model on the basis of OS . Five immune features were identified in the training cohort: CD3T, CD27T, CD103T, PD1T, and CD68P (Fig. 1e–g). The IHC expression modes of the 5 immune features in HCC tissues are presented in Fig. 1h. And then, we built a personalized TRIS with the coefficients weighted by the LASSO Cox model in the training cohort, where TRIS = (175.6642 × the level of CD68P − 14.7372 × the level of CD3T − 37.5718 × the level of CD27T − 88.1263 × the level of CD103T − 50.5161 × the level of PD1T) × 10−5. The level of immune features was defined as the number of positively stained cells in specific regions of each HCC patient: intra- or peri-tumoral tissues. We then evaluated the correlation of TRIS with clinicopathological variables and observed that the TRIS was statistically associated with sex (P = 0.02), AFP (P = 0.009), tumor diameter (P = 0.02), and tumor numbers (P = 0.046) (Additional file 2: Table S1).
Selection of prognostic predictors
Univariate analyses were performed on age, sex, etiology, TRIS, and clinicopathologic variables to determine their associations with OS in patients with early/intermediate-stage HCC. All variables that were significant were evaluated with multivariate analyses (Table 2). Multivariate analysis revealed that γ-glutamyl transferase (GGT) (hazard ratio (HR), 1.002; 95% CI, 1.000–1.004; P = 0.01), tumor diameter (HR, 1.100; 95% CI, 1.046–1.156; P < 0.001), tumor differentiation (Edmondson-Steiner grade I-II/III-IV) (HR, 1.819; 95% CI, 1.268–2.610; P = 0.001) , and TRIS (HR, 4.606; 95% CI, 1.335–15.890; P = 0.02) were independent prognostic factors of OS in HCC patients.
Establishment of the ICPI
To further improve the accuracy of survival prediction, GGT, TRIS, tumor diameter, and tumor differentiation were integrated. By using the Cox proportional hazards regression model, we then derived an individualized ICPI equation: (0.0889 × GGT + 5 × tumor diameter + 25.4939 × tumor differentiation status + 60.9038 × TRIS − 1.1324). In this formula, Edmondson-Steiner grade I-II was defined as 1, and grade III-IV as 2. And then, we evaluated the predictive value of ICPI with other postoperative variables. The data revealed that the performance of ICPI was superior to that of postoperative variables, including microvascular invasion, lymphoid metastasis, and tumor differentiation (Additional file 2: Table S2).
To compare the predictive power of the ICPI and 7 traditional staging systems, ROC curve analysis was applied. In patients with HBV-associated HCC, our established ICPI achieved significantly improved estimation in survival prediction (C-index, 0.691) when compared with the 7 staging systems in the training cohort (C-index, 0.548–0.597) (Fig. 2a and Additional file 2: Table S3).
Performance of the ICPI in stratifying the risk of patients
To determine the optimal cutoff values of the ICPI, X-tile program  was used in the training cohort (Additional file 3: Figure S4). Then, the patients were stratified into 3 subgroups: score 1, 0 to 58.5; score 2, 58.5 to 86.2; and score 3, ≥ 86.2. In the training cohort, stratification into the 3 subgroups allowed significant distinction for survival outcomes (score 1 vs. 2, P = 0.001; score 2 vs. 3, P < 0.001). Except for Okuda staging system, no significant differences in survival distribution were observed across all stages of BCLC, CLIP, JIS, LCSGJ, and AJCC 7th and 8th edition staging systems (Additional file 3: Figure S5). Further, we assessed the performance of the ICPI model and 7 staging systems with corrected AIC values and homogeneity . Among the eight staging systems, our stratified ICPI revealed the highest homogeneity (43.66) and the lowest AIC value (1594.64) in HBV-related HCC patients in the early/intermediate stage (Table 3), implying that stratified ICPI might be more accurate in predicting postoperative survival.
Validation of constructed ICPI
To substantiate the performance of ICPI, validation analyses were performed in an internal validation cohort (n = 393). In the validation cohort, the C-index value of the ICPI was superior than those of 7 other staging systems [ICPI vs Okuda, 0.686 vs 0.602 (P < 0.001); CLIP, 0.610 (P = 0.01); LCSGJ, 0.553 (P < 0.001); JIS, 0.556 (P < 0.001); AJCC 7th edition, 0.519 (P < 0.001); AJCC 8th edition, 0.573 (P < 0.001); BCLC, 0.544 (P < 0.001), respectively] (Fig. 2 and Additional file 2: Table S3).
Moreover, we performed the Kaplan-Meier survival analysis according to the stratified ICPI and traditional staging systems in the validation cohort. Figure 3 indicates that significant differences in survival distributions were found across all stages of stratified ICPI score 1/2 (P = 0.003), score 2/3 (P < 0.001), Okuda stage I/II (P < 0.001), and AJCC 8th edition stage I/II (P = 0.01) and II/III (P = 0.02). There were no significant survival differences between BCLC stage A/B (P = 0.38), CLIP 0/1 (P = 0.09), 2/3 (P = 0.68), JIS 1/2 (P = 0.14), LCSGJ II/III (P = 0.13), AJCC 7th edition I/II (P = 0.31), and II/III (P = 0.24).
We also assessed corrected AIC values and homogeneity in the validation cohort. The stratified ICPI remains the highest homogeneity (50.94) and lowest AIC (1551.86) (Table 3). In addition, we generated a nomogram to facilitate the clinical application of ICPI (Additional file 3: Figure S6A). The calibration plots achieved acceptable agreement in the validation cohort between the ICPI prediction and the actual observation for 1-, 3-, and 5-year OS (Additional file 3: Figure S6B).
Relationship between the ICPI and local immune status
To investigate the interaction of immune features in the tumor microenvironment, an immune network was constructed on the basis of TMA data. Figure 4a shows a clear separation between the two regions, each exhibiting a distinct and characteristic immune cell pattern, and subnetworks of T cell subpopulations (CD3 cells, CD4 cells, memory T cells, and naïve T cells) in liver tissues. We also observed a correlation between B cells (CD20) and the T cell subset network, implying a functional interaction between B cells and T cells. Interestingly, close correlations between memory T cells (CD45RO) and 5 other immune features, including CD3, CD4, CD27, CD45RA, and CD8, existed in both regions. The local interaction of the immune network suggested the existence of tumor-microenvironment compartments with different compositions, which might influence the activity and mobility of T and B lymphocytes in the tumor progression.
The local immune status plays an essential role in carcinogenesis and in response to cancer therapeutics. Previous reports showed that the intratumoral densities of CD8+, CD57+, and CD45RO+ cells correlated with the local immune status [34, 42, 43]. Then, we investigated whether our constructed ICPI could be a potential indicator of the local immune response. The correlations between the ICPI and intratumoral immune features were assessed. Figure 4 b and c reveal that the levels of CD4+, CD20+, CD45RA+, and CD45RO+cells were inversely associated with individual ICPI. And also, the relationships between the levels of CD4+, CD20+, CD45RA+, and CD45RO+ cells and the TRIS shared similar trends (Additional file 3: Figure S7), implying that the ICPI was a useful prognostic mode and might be correlated with the local immune status.
HCC is a remarkably heterogeneous type of malignancy resulting from the accumulation of epigenetic and genomic alterations and is influenced by tumor-host crosstalk [44, 45]. With inadequate information on the biological behavior and microenvironment characteristics of liver cancer, it may not be accurate enough to simply combine few clinicopathological parameters together for the prediction of survival in HCC patients. In this study, an immune and clinical prognostic index was developed and validated in an independent dataset. The newly developed ICPI was useful in survival prediction for patients with early/intermediate-stage HCC (BCLC stages 0, A, and B).
To identify prognosis-related immune predictors, 28 immune features were evaluated in the training cohort of 352 HCC specimens. In contrast to other studies that used gene expression profiles and flow cytometry, IHC was used for the identification of tumor-infiltrating immune cells, because of its ease of use and clinically applicable. Based on all the examination results, five immune features were selected through LASSO Cox method in our study, including CD3T, CD27T, CD103T, PD1T, and CD68P. This selection is more comprehensive than those previously reported studies in multiple aspects [25, 46]. First, we systematically investigated 28 immunological features based on intra- and peri-tumoral tissues in liver cancer that are not confined to the commonly used CD3, CD4, CD45RO, or CD8. Second, the LASSO Cox regression model was adopted to achieve covariate selection. Also, in contrast to previous studies with relatively smaller sample sizes (from 65 to 348 patients) [25, 26, 46], our study recruited a much larger cohort (n = 745 patients).
The prognosis values of the five immune features selected for TRIS were consistent with previous studies. Our study indicated that a low density of intratumoral CD3+ cells predicted poor prognosis, which is similar to studies on HCC  and gastric cancer . Consistent with the study by Garnelo et al. , high density of CD27-positive cells was associated with superior HCC patients’ survival in our study. Also, we confirmed the interaction between intratumoral B (CD20) cells and T (CD27) cells with immune cell network. In tumor-infiltrating lymphocytes and peripheral blood of HCCs, Zhang et al.  revealed that the prevalence of Tregs expressing high levels of HLA-DR, GITR, and CD103 increased during the progressive stages. In accordance with our previous study , we observed that high infiltration of PD1+ immune cells could predict dismal survival in HCC patients. Furthermore, it has also been reported that Kupffer cells (CD68), the resident macrophages in the liver, play an important role in HCC development and hepatic homeostasis . Similar to the previous study by Zhu et al. , high levels of peritumoral macrophages were correlated with tumor progression and poor prognosis after liver resection, indicating the potential value of peritumoral immune components to serve as a prognosis factor and future target of immunotherapy.
Apart from the selected immune biomarkers of CD3, CD27, CD103, CD68, and PD1, the status and prognostic significance of other immune biomarkers were also evaluated in this study. In HCC, the antitumor effector functions of CD4+ T and CD8+ T cells have been reported , while the role of B lymphocytes in the carcinogenesis of HCC remains controversial . The antitumor function of NK cells has been well-documented, which is gradually diminished with the progression of HCC. In our study, no significant differences in NK cells (CD57) were observed between intra- and peri-tumoral tissues. In our previous study, intratumoral neutrophils, facilitating tumorigenesis and resistance to the tyrosine-kinase inhibitor sorafenib, were reported as a poor indicator for HCC patients . Similar to the studies on colorectal cancer , the intratumoral memory T cells (CD45RO) were also decreased dramatically in our study. In addition, depletion of CD14+ and CD45RA+ cells in the tumor regions was observed. However, the intratumoral CXCR5+ cells did not show any significant changes.
The host immunity against tumors is considered increasingly important in clinical practice and cancer research as the local immune status can influence tumor evolution. Changes in density of intratumoral-infiltrated immune cells during tumor invasion and metastasis may indicate progressive immune escape. In our study, the ICPI was inversely but weakly correlated with CD4+, CD45RA+, CD45RO+, and CD20+cells. Previous studies have reported that the intratumoral immune response was associated with the densities of localized CD3+, CD8+, and CD45RO+ cells . Garnelo et al. reported that tumor-infiltrating CD20+ and CD27+ cells enhanced local immune activation and contributed to a better prognosis for patients with HCC . In addition, CD103+ cells have been proposed as an indicator of the favorable OS in breast and lung cancers [50, 51]. These findings suggest that the combination of selected markers in the ICPI may serve as a potential indicator which represents the localized immune status in liver cancer.
The survival of HCC patients depends on the tumor stage, underlying liver function, and performance status. As a marker of liver injury, GGT was an independent prognostic predictor on multivariate analysis in our study. Recently, other researchers reported that elevated serum GGT was associated with the risk of cancer, especially in liver cancer . Our final ICPI model also integrated tumor diameter and histologic differentiation as reported in previous studies [53, 54]. Based on these clinical predictors, our established ICPI integrated the characteristics of liver function, liver pathology, and host-tumor interaction.
Our study has limitations. First, this study was focused on HBV-associated HCC in the early/intermediate stage. It requires further study whether the ICPI can be applied to patients with advanced-staged disease or HCC with other etiologies requires further studies. Second, the study is retrospective with its inherent defects. Third, the study cohorts were recruited from one single center, and approximately 10% of HCC patients at intermediate stage (BCLC B) underwent surgical resection, which might result in selection bias. In addition, the biological mechanisms remain to be further investigated about how the candidate markers, including CD3, CD27, CD68, CD103, and PD1, are involved in HCC. Finally, the ICPI model was constructed on the basis of pre- and postoperative parameters, which may not be suitable for preoperative decision-making.
In summary, the ICPI improved the accuracy of survival prediction for patients with early/intermediate-stage HCC. Prospective studies are needed to further validate its analytical accuracy in estimating prognosis for individualized management of HCC patients.
Akaike information criterion
American Joint Committee on Cancer
Barcelona Clinic Liver Cancer
Carbohydrate antigen 19-9
Cancer of the Liver Italian Program
Gene Expression Omnibus
Immune-clinical prognostic index
Japan Integrated Staging
Least absolute shrinkage and selection operator
Liver Cancer Study Group of Japan
Receiver operating characteristic curve
The Cancer Genome Atlas
T follicular helper
T regulatory cell
Tissue-related immune signature
Torre LA, Siegel RL, Ward EM, Jemal A. Global cancer incidence and mortality rates and trends--an update. Cancer Epidemiol Biomark Prev. 2016;25(1):16–27.
Johnson PJ, Berhane S, Kagebayashi C, Satomura S, Teng M, Reeves HL, O'Beirne J, Fox R, Skowronska A, Palmer D, et al. Assessment of liver function in patients with hepatocellular carcinoma: a new evidence-based approach-the ALBI grade. J Clin Oncol. 2015;33(6):550–8.
Torrecilla S, Sia D, Harrington AN, Zhang Z, Cabellos L, Cornella H, Moeini A, Camprecios G, Leow WQ, Fiel MI, et al. Trunk mutational events present minimal intra- and inter-tumoral heterogeneity in hepatocellular carcinoma. J Hepatol. 2017;67(6):1222–31.
Ringelhan M, Pfister D, O'Connor T, Pikarsky E, Heikenwalder M. The immunology of hepatocellular carcinoma. Nat Immunol. 2018;19(3):222–32.
Bilzer M, Roggel F, Gerbes AL. Role of Kupffer cells in host defense and liver disease. Liver Int. 2006;26(10):1175–86.
Ouyang FZ, Wu RQ, Wei Y, Liu RX, Yang D, Xiao X, Zheng L, Li B, Lao XM, Kuang DM. Dendritic cell-elicited B-cell activation fosters immune privilege via IL-10 signals in hepatocellular carcinoma. Nat Commun. 2016;7:13453.
Kugler DG, Flomerfelt FA, Costa DL, Laky K, Kamenyeva O, Mittelstadt PR, Gress RE, Rosshart SP, Rehermann B, Ashwell JD, et al. Systemic toxoplasma infection triggers a long-term defect in the generation and function of naive T lymphocytes. J Exp Med. 2016;213(13):3041–56.
Garnelo M, Tan A, Her Z, Yeong J, Lim CJ, Chen J, Lim KH, Weber A, Chow P, Chung A, et al. Interaction between tumour-infiltrating B cells and T cells controls the progression of hepatocellular carcinoma. Gut. 2017;66(2):342–51.
Zheng C, Zheng L, Yoo JK, Guo H, Zhang Y, Guo X, Kang B, Hu R, Huang JY, Zhang Q, et al. Landscape of infiltrating T cells in liver cancer revealed by single-cell sequencing. Cell. 2017;169(7):1342–56 e1316.
Raziorrouh B, Sacher K, Tawar RG, Emmerich F, Neumann-Haefelin C, Baumert TF, Thimme R, Boettler T. Virus-specific CD4+ T cells have functional and phenotypic characteristics of follicular T-helper cells in patients with acute and chronic HCV infections. Gastroenterology. 2016;150(3):696–706.e693.
Huang LR, Wohlleber D, Reisinger F, Jenne CN, Cheng RL, Abdullah Z, Schildberg FA, Odenthal M, Dienes HP, van Rooijen N, et al. Intrahepatic myeloid-cell aggregates enable local proliferation of CD8(+) T cells and successful immunotherapy against chronic viral liver infection. Nat Immunol. 2013;14(6):574–83.
Ma C, Kesarwala AH, Eggert T, Medina-Echeverz J, Kleiner DE, Jin P, Stroncek DF, Terabe M, Kapoor V, ElGindi M, et al. NAFLD causes selective CD4(+) T lymphocyte loss and promotes hepatocarcinogenesis. Nature. 2016;531(7593):253–7.
Fridman WH, Pages F, Sautes-Fridman C, Galon J. The immune contexture in human tumours: impact on clinical outcome. Nat Rev Cancer. 2012;12(4):298–306.
Condamine T, Ramachandran I, Youn JI, Gabrilovich DI. Regulation of tumor metastasis by myeloid-derived suppressor cells. Annu Rev Med. 2015;66:97–110.
Gentles AJ, Newman AM, Liu CL, Bratman SV, Feng W, Kim D, Nair VS, Xu Y, Khuong A, Hoang CD, et al. The prognostic landscape of genes and infiltrating immune cells across human cancers. Nat Med. 2015;21(8):938–45.
Mlecnik B, Tosolini M, Kirilovsky A, Berger A, Bindea G, Meatchi T, Bruneval P, Trajanoski Z, Fridman WH, Pages F, et al. Histopathologic-based prognostic factors of colorectal cancers are associated with the state of the local immune reaction. J Clin Oncol. 2011;29(6):610–8.
Li YW, Qiu SJ, Fan J, Zhou J, Gao Q, Xiao YS, Xu YF. Intratumoral neutrophils: a poor prognostic factor for hepatocellular carcinoma following resection. J Hepatol. 2011;54(3):497–505.
Shi JY, Gao Q, Wang ZC, Zhou J, Wang XY, Min ZH, Shi YH, Shi GM, Ding ZB, Ke AW, et al. Margin-infiltrating CD20(+) B cells display an atypical memory phenotype and correlate with favorable prognosis in hepatocellular carcinoma. Clin Cancer Res. 2013;19(21):5994–6005.
Gao Q, Qiu SJ, Fan J, Zhou J, Wang XY, Xiao YS, Xu Y, Li YW, Tang ZY. Intratumoral balance of regulatory and cytotoxic T cells is associated with prognosis of hepatocellular carcinoma after resection. J Clin Oncol. 2007;25(18):2586–93.
Sia D, Jiao Y, Martinez-Quetglas I, Kuchuk O, Villacorta-Martin C, Castro de Moura M, Putra J, Camprecios G, Bassaganyas L, Akers N, et al. Identification of an immune-specific class of hepatocellular carcinoma, based on molecular features. Gastroenterology. 2017;153(3):812–26.
Li H, Li X, Liu S, Guo L, Zhang B, Zhang J, Ye Q. Programmed cell death-1 (PD-1) checkpoint blockade in combination with a mammalian target of rapamycin inhibitor restrains hepatocellular carcinoma growth induced by hepatoma cell-intrinsic PD-1. Hepatology. 2017;66(6):1920–33.
Zhou J, Sun HC, Wang Z, Cong WM, Wang JH, Zeng MS, Yang JM, Bie P, Liu LX, Wen TF, et al. Guidelines for diagnosis and treatment of primary liver cancer in China (2017 edition). Liver Cancer. 2018;7(3):235–60.
Zhu XD, Zhang JB, Zhuang PY, Zhu HG, Zhang W, Xiong YQ, Wu WZ, Wang L, Tang ZY, Sun HC. High expression of macrophage colony-stimulating factor in peritumoral liver tissue is associated with poor survival after curative resection of hepatocellular carcinoma. J Clin Oncol. 2008;26(16):2707–16.
Fu J, Zhang Z, Zhou L, Qi Z, Xing S, Lv J, Shi J, Fu B, Liu Z, Zhang JY, et al. Impairment of CD4+ cytotoxic T cells predicts poor survival and high recurrence rates in patients with hepatocellular carcinoma. Hepatology. 2013;58(1):139–49.
Gabrielson A, Wu Y, Wang H, Jiang J, Kallakury B, Gatalica Z, Reddy S, Kleiner D, Fishbein T, Johnson L, et al. Intratumoral CD3 and CD8 T-cell densities associated with relapse-free survival in HCC. Cancer Immunol Res. 2016;4(5):419–30.
Gao Q, Zhou J, Wang XY, Qiu SJ, Song K, Huang XW, Sun J, Shi YH, Li BZ, Xiao YS, et al. Infiltrating memory/senescent T cell ratio predicts extrahepatic metastasis of hepatocellular carcinoma. Ann Surg Oncol. 2012;19(2):455–66.
Takata Y, Nakamoto Y, Nakada A, Terashima T, Arihara F, Kitahara M, Kakinoki K, Arai K, Yamashita T, Sakai Y, et al. Frequency of CD45RO+ subset in CD4+CD25(high) regulatory T cells associated with progression of hepatocellular carcinoma. Cancer Lett. 2011;307(2):165–73.
Zhou SL, Zhou ZJ, Hu ZQ, Huang XW, Wang Z, Chen EB, Fan J, Cao Y, Dai Z, Zhou J. Tumor-associated neutrophils recruit macrophages and T-regulatory cells to promote progression of hepatocellular carcinoma and resistance to sorafenib. Gastroenterology. 2016;150(7):1646–1658.e1617.
Asai A, Tsuchimoto Y, Ohama H, Fukunishi S, Tsuda Y, Kobayashi M, Higuchi K, Suzuki F. Host antitumor resistance improved by the macrophage polarization in a chimera model of patients with HCC. Oncoimmunology. 2017;6(4):e1299301.
Wang L, Sun Y, Zhang Z, Jia Y, Zou Z, Ding J, Li Y, Xu X, Jin L, Yang T, et al. CXCR5+ CD4+ T follicular helper cells participate in the pathogenesis of primary biliary cirrhosis. Hepatology. 2015;61(2):627–38.
Yi Y, He HW, Wang JX, Cai XY, Li YW, Zhou J, Cheng YF, Jin JJ, Fan J, Qiu SJ. The functional impairment of HCC-infiltrating gammadelta T cells, partially mediated by regulatory T cells in a TGFbeta- and IL-10-dependent manner. J Hepatol. 2013;58(5):977–83.
Zhang HH, Mei MH, Fei R, Liao WJ, Wang XY, Qin LL, Wang JH, Wei L, Chen HS. Regulatory T cell depletion enhances tumor specific CD8 T-cell responses, elicited by tumor antigen NY-ESO-1b in hepatocellular carcinoma patients, in vitro. Int J Oncol. 2010;36(4):841–8.
Li S, Xu F, Li H, Zhang J, Zhong A, Huang B, Lai M. S100A8(+) stroma cells predict a good prognosis and inhibit aggressiveness in colorectal carcinoma. Oncoimmunology. 2017;6(1):e1260213.
Jiang Y, Zhang Q, Hu Y, Li T, Yu J, Zhao L, Ye G, Deng H, Mou T, Cai S, et al. ImmunoScore signature: a prognostic and predictive tool in gastric cancer. Ann Surg. 2018;267(3):504–13.
Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, Hoang CD, Diehn M, Alizadeh AA. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. 2015;12(5):453–7.
Qiu J, Peng B, Tang Y, Qian Y, Guo P, Li M, Luo J, Chen B, Tang H, Lu C, et al. CpG methylation signature predicts recurrence in early-stage hepatocellular carcinoma: results from a multicenter study. J Clin Oncol. 2017;35(7):734–42.
Li J, Zhou J, Yang PH, Xia Y, Shi YH, Wu D, Lv G, Zheng W, Wang K, Wan XY, et al. Nomograms for survival prediction in patients undergoing liver resection for hepatitis B virus related early stage hepatocellular carcinoma. Eur J Cancer. 2016;62:86–95.
Shim JH, Jun MJ, Han S, Lee YJ, Lee SG, Kim KM, Lim YS, Lee HC. Prognostic nomograms for prediction of recurrence and survival after curative liver resection for hepatocellular carcinoma. Ann Surg. 2015;261(5):939–46.
Camp RL, Dolled-Filhart M, Rimm DL. X-tile: a new bio-informatics tool for biomarker assessment and outcome-based cut-point optimization. Clin Cancer Res. 2004;10(21):7252–9.
Liu PH, Hsu CY, Hsia CY, Lee YH, Su CW, Huang YH, Lee FY, Lin HC, Huo TI. Prognosis of hepatocellular carcinoma: Assessment of eleven staging systems. J Hepatol. 2016;64(3):601–8.
Edmondson HA, Steiner PE. Primary carcinoma of the liver: a study of 100 cases among 48,900 necropsies. Cancer. 1954;7(3):462–503.
Brunner SM, Rubner C, Kesselring R, Martin M, Griesshammer E, Ruemmele P, Stempfl T, Teufel A, Schlitt HJ, Fichtner-Feigl S. Tumor-infiltrating, interleukin-33-producing effector-memory CD8(+) T cells in resected hepatocellular carcinoma prolong patient survival. Hepatology. 2015;61(6):1957–67.
Geissler K, Fornara P, Lautenschlager C, Holzhausen HJ, Seliger B, Riemann D. Immune signature of tumor infiltrating immune cells in renal cancer. Oncoimmunology. 2015;4(1):e985082.
Cancer Genome Atlas Research Network. Electronic address web, Cancer Genome Atlas Research N: Comprehensive and integrative genomic characterization of hepatocellular carcinoma. Cell 2017, 169(7):1327–1341 e1323.
Raggi C, Factor VM, Seo D, Holczbauer A, Gillen MC, Marquardt JU, Andersen JB, Durkin M, Thorgeirsson SS. Epigenetic reprogramming modulates malignant properties of human liver cancer. Hepatology. 2014;59(6):2251–62.
Fu YP, Yi Y, Cai XY, Sun J, Ni XC, He HW, Wang JX, Lu ZF, Huang JL, Cao Y, et al. Overexpression of interleukin-35 associates with hepatocellular carcinoma aggressiveness and recurrence after curative resection. Br J Cancer. 2016;114(7):767–76.
Ma LJ, Feng FL, Dong LQ, Zhang Z, Duan M, Liu LZ, Shi JY, Yang LX, Wang ZC, Zhang S, et al. Clinical significance of PD-1/PD-Ls gene amplification and overexpression in patients with hepatocellular carcinoma. Theranostics. 2018;8(20):5690–702.
Tacke F. Targeting hepatic macrophages to treat liver diseases. J Hepatol. 2017;66(6):1300–12.
Galon J, Costes A, Sanchez-Cabo F, Kirilovsky A, Mlecnik B, Lagorce-Pages C, Tosolini M, Camus M, Berger A, Wind P, et al. Type, density, and location of immune cells within human colorectal tumors predict clinical outcome. Science. 2006;313(5795):1960–4.
Wang ZQ, Milne K, Derocher H, Webb JR, Nelson BH, Watson PH. CD103 and intratumoral immune response in breast cancer. Clin Cancer Res. 2016;22(24):6290–7.
Ganesan AP, Clarke J, Wood O, Garrido-Martin EM, Chee SJ, Mellows T, Samaniego-Castruita D, Singh D, Seumois G, Alzetani A, et al. Tissue-resident memory features are linked to the magnitude of cytotoxic T cell responses in human lung cancer. Nat Immunol. 2017;18(8):940–50.
Mok Y, Son DK, Yun YD, Jee SH, Samet JM. Gamma-glutamyltransferase and cancer risk: the Korean cancer prevention study. Int J Cancer J Int du Cancer. 2016;138(2):311–9.
Minagawa M, Ikai I, Matsuyama Y, Yamaoka Y, Makuuchi M. Staging of hepatocellular carcinoma: assessment of the Japanese TNM and AJCC/UICC TNM systems in a cohort of 13,772 patients in Japan. Ann Surg. 2007;245(6):909–22.
Utsunomiya T, Shimada M, Kudo M, Ichida T, Matsui O, Izumi N, Matsuyama Y, Sakamoto M, Nakashima O, Ku Y, et al. Nationwide study of 4741 patients with non-B non-C hepatocellular carcinoma with special reference to the therapeutic impact. Ann Surg. 2014;259(2):336–45.
We thank Jun Jiang and Wei Le for the bioinformatic analysis, Zun-Song Hu for the statistical analysis, and Hai-Ying Zeng and Dong-Xian Jiang for the pathological analysis.
This work was supported by the grants from the National High Technology Research and Development Program (863 Program) of China (2015AA020401), National Natural Science Foundation of China (No. 81472674, 81530077, 81502486, 81773067 and 81800790), and the Shanghai Committee of Science and Technology, China (No. 16JC1404000). The work was sponsored by the Shanghai Sailing Program (19YF1407800), Shu Guang project of Shanghai Municipal Education Commission and Shanghai Education Development Foundation (13SG04), Shanghai Municipal Science and Technology Major Project (Grant No. 2018SHZDZX05), and the China Postdoctoral Science Foundation (2018 M640343).
Availability of data and materials
The data that support the findings of this study are available from the Department of Liver Surgery in Zhongshan Hospital but restrictions apply to the availability of these data, which were used under license for the current study, and so are not publicly available. Data are however available from the authors upon reasonable request and with permission of the Department of Liver Surgery in Zhongshan Hospital.
Ethics approval and consent to participate
Written informed consent was obtained from each patient. The study was approved by the Institutional Review Board of Zhongshan Hospital and complied with the standards of the Declaration of Helsinki and current ethical guidelines.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary methods. (DOCX 40 kb)
Table S1. The relationship between TRIS and clinicopathological variables. Table S2. Comparison of prognostic performance among postoperative variables and ICPI. Table S3. Comparison of prognostic performance among HCC staging systems and ICPI. Table S4. Components of 7 staging systems for hepatocellular carcinoma. Table S5. Antibody sources and staining conditions. (DOCX 28 kb)
Figure S1. Study flowchart. Figure S2. Digital image analyzed using the image software (Image-pro plus 6.0), with tissue represented in yellow and stained cells represented in red. Figure S3. Immunohistochemistry expression pattern of 14 immune infiltrations in tumor and adjacent liver tissues, including CD3, CD4, CD8, CD14, CD20, CD27, CD45RA, CD45RO, CD57, CD66b, CD68, CD103, CXCR5, and PD1. Figure S4. X-tile plots of ICPI in the training cohort automatically selecting the optimum cut point according to the highest χ2 value defined by the Kaplan-Meier survival analysis and log-rank test. Figure S5. Kaplan-Meier survival curves of the training cohort categorized by different staging systems. Figure S6. (A) Nomogram for predicting the survival probability in HCC patients. (B) Calibration of the predictive models at 1, 3, and 5 years in the derivation and validation cohorts. Figure S7. The correlation between TRIS and the density of intratumoral immune features, including CD4+, CD20+, CD45RO+, and CD45RA+ cells. (DOCX 5565 kb)