- Letter to the Editor
- Open Access
Integrated proteogenomic analysis revealed the metabolic heterogeneity in noncancerous liver tissues of patients with hepatocellular carcinoma
Journal of Hematology & Oncology volume 14, Article number: 205 (2021)
Understanding the adjacent liver microenvironment of hepatocellular carcinoma (HCC) with possible metastasis tendency might provide a strategy for risk classification of patients and potential therapies by converting the unique metastasis-inclined microenvironment to a metastasis-averse one. In this study, we performed an integrated proteogenomic analysis to have a comprehensive view on the heterogeneity of hepatic microenvironment contributing to HCC metastasis. Pairing mRNA-protein analysis revealed consistent and discordant mRNA-protein expressions in metabolism regulations and cancer-related pathways, respectively. Proteomic profiling identified three subgroups associated with the recurrence-free survival of patients. These proteomic subgroups demonstrated distinct features in metabolic reprogramming, which was potentially modified by epigenetic alterations. This study raises the point of metabolic heterogeneity in HCC noncancerous tissues and may offer a new perspective on HCC treatment.
To the Editor,
Hepatocellular carcinoma (HCC) is an aggressive malignancy characterized by high rate of relapse after surgery due to intrahepatic metastasis. Nevertheless, there are still patients who undergo surgical resection having favorable outcomes with no detectable recurrence during follow-up, suggesting the heterogeneity in local microenvironment of noncancerous sites dictating the ability of a tumor to metastasize [1, 2]. On the other hand, HCC is usually present in the background of inflamed fibrotic and/or cirrhotic liver with extensive lymphocyte infiltration caused by chronic hepatitis, indicating that the metastatic propensity of HCC might be influenced by local tissue microenvironment of the host . Indeed, several cellular and acellular components in liver microenvironment have been highlighted in the development or potential regression of HCC hepatic metastasis . Based on the above reasons, understanding the liver microenvironment with possible metastasis tendency might provide a strategy for risk classification of patients and potential therapies of metastatic HCC by converting the unique metastasis-inclined phenotype (MI) to a metastasis-averse (MA) one.
To comprehensively understand the heterogeneity of hepatic microenvironment contributing to HCC metastasis, we performed a proteogenomic analyses integrating epigenomic, transcriptomic and proteomic data. A total of 77 noncancerous liver tissues from 77 HCC patients were selected in current study (Fig. 1a). After pathological inspection and quality controls during library preparation, the availability of methylation, RNA-sequencing (RNA-seq) and proteome data for each HCC candidate is listed in Fig. 1a (lower panel), in which > 96% patients have all of the three omics data.
Initially, 6,350 proteins were identified by multiplexed tandem mass tag (TMT) mass spectrometry (MS) with an average of 5,331 in each sample. After removing proteins with missing value in > 50% of the samples, 5,445 proteins with imputed value were selected for the following analyses. Pairing transcriptomic and proteomic data from 74 patients resulted in 3,834 mRNA-protein pairs, which showed an overall correlation close to 0 with only 16.1% (617/3,834) significant positive correlations (multiple-test adjusted P < 0.01, Fig. 1b). This fact highly suggested a large amount of post-transcriptional modifications in noncancerous liver tissues. Table S1 (Additional file 3) listed the top 20 genes with the highest mRNA-protein correlations, among which a larger number of metabolism-related genes were identified, such as FN3KRP, GSTT1, AKR1C2, DPYS, SULT2A1 and GSTZ1. Strikingly, genes involved in various metabolic processes had the strongest positive mRNA-protein correlations (Fig. 1c; Additional file 1: Table S2), which was consistent with the findings in HCC tumor tissues . Notably, mRNA-protein pairs with significant discordant expression were involved in various cancer-related pathways (Additional file 1: Fig. S1 and Table S2), which might be owing to the high level of somatic mutation burden in cirrhotic liver , potentially leading to the aberrant post-transcriptional regulation of these genes . Despite that co-expression network analysis identified a major functional module consisting of multiple metabolic genes in RNA-seq network, no enrichment module was detected in proteomic network (Fig. 1d). These facts highlighted the necessity of using proteome data to reflect the authentic expression profile of liver microenvironment.
In order to evaluate the heterogeneity and homogeneity among noncancerous samples, we then performed unsupervised clustering based on top 25% most varied proteins among the included samples (Additional file 1: supplementary methods and Fig. S2) and identified three subgroups among the 77 samples (Fig. 2a). While S1 and S2 patients have similar expression patterns, there is a perceptible difference between S3 and each of them (Kruskal–Wallis test). Clinicopathologic parameters such as higher AFP level, tumor thrombus and advanced TNM stages were more prominent in S3 than the other two subgroups. Through manual curation (full list in Table S3, Additional file 3), we found that each subgroup was characterized by various level of metabolism-related proteins. For example, fructose-bisphosphate aldolase C, which is involved in numerous metabolism pathways , was ranked in the top 20 proteins in Fig. 2a (Additional file 3: Table S4). To explore this inter-sample heterogeneity, we quantified the activity of a metabolic pathways in each subgroup using a pathway activity score algorithm as previously described . Indeed, significant metabolic heterogeneity was found among these three subgroups, each of which demonstrated a unique profile of metabolic activity (Fig. 2b). By dividing patients into two groups based on the occurrence of early relapse (within 2 years after surgery), we found that the proportion of early relapse samples in S2 subgroups was significantly different from both S1 and S3 subgroups (Fig. 2c, left panel). Moreover, the proteomic subgroups significantly had different recurrence-free survival within 2 years after surgery (Fig. 2c, right panel). Multivariable analysis authenticated the role of proteomic subgroups as independent prognosticators after adjusted for clinicopathologic factors. After stratifying patients according to TNM stage, however, proteomic subgroups showed no correlation with patient prognosis (Fig. 2d, right panel), supporting the effect from tumor tissue on the molecular features of corresponding noncancerous liver tissue. Nevertheless, the proportion of early relapse samples in S2 subgroups was still significantly different from both S1 and S3 subgroups regardless of tumor stages (Fig. 2d, left panel). Meanwhile, the abundance of p62 protein, high expression of which in non-tumor human liver has been proved to be able to predict rapid HCC recurrence after curative treatment , was also compared among these three subgroups. The result showed distributed difference of p62 abundance among proteomic subgroups (Fig. 2e), which further validated the findings in our survival analyses. Given that tumor thrombus is one of the most important recurrence-related factors, we compared proteomic profiles between samples with or without tumor thrombus, which revealed 447 differentially expressed proteins (Fig. 2f), among which 31 metabolism-related proteins and p62 protein were identified. Moreover, Gene Set Enrichment Analysis (GSEA) also showed that the differentially expressed proteins were enriched in metabolism-related pathways (Fig. 2g). Based on the above, a crucial dysregulation of metabolism was also observed in noncancerous tissues, which was similar with the situation in tumor tissues .
Although it was not the main purpose of this study, we also conducted a supervised analysis to identify representative prognostic protein markers (Additional file 2: Fig. S3a), which resulted in two proteins, GNG7 (G Protein Subunit Gamma 7) and MCIDAS (Multiciliate Differentiation And DNA Synthesis Associated Cell Cycle Protein), showing differential expression between MI and MA samples (defined by Budhu et al. , see Supplementary methods in Additional file 1). Using median as the cutoff (Additional file 1: Supplementary methods), both two proteins demonstrated significant predicting values on the recurrence-free survival, which was further confirmed by multivariable analysis (Additional file 2: Fig. S3b). Regarding the potential function of these two markers, tumors with high GNG7 showed significant upregulation of various pathways relevant to cell signaling and metabolism (Additional file 2: Fig. S3c), while those with high MCIDAS were featured with upregulation of pathways related to tumor progressions, cell cycle, immune response and metabolism (Additional file 2: Fig. S3d).
Intrahepatic metastasis of HCC is accompanied by alterations in the immune status of the tumor-surrounding issue , suggesting the potential relationship between HCC recurrence and liver immune microenvironment. To further depict the microenvironment of the noncancerous livers, we enumerated cell subsets in each proteomic subgroup from transcriptomes. While S1 and S2 samples demonstrated similar enrichment for each cell type, perceptible difference was also found in the proportion of cell component between S3 and S1/S2 samples (Additional file 2: Fig. S4a). Since no difference in stroma scores was observed among the three subgroups (Additional file 2: Fig. S4c), the difference in microenvironment score (Additional file 2: Fig. S4d) was mainly contributed by the heterogeneity in immune cell infiltration (Additional file 2: Fig. S4b).
At the same time, we also evaluated the epigenetic regulation of gene expression by using methylation data. A total of 10,9624 differentially methylated regions (DMRs) across all autosomes were detected among the three subgroups (Additional file 2: Fig. S5a), the majority of which are located within intergenic regions (Additional file 2: Fig. S5b). Since hypermethylated DMRs are expected to be associated with decreased gene expression, we analyzed the potential effect of DMRs on mRNA and protein expressions using a Gaussian mixture model, which identified 1,664 proteins significantly attenuated by DMRs regions (Additional file 2: Fig. S5c), corresponding to 43.4% of all the 3,834 genes analyzed. Intriguingly, these attenuated proteins were mainly enriched in various metabolism pathways (Additional file 2: Fig. S5d), which highly suggested that expressions of these metabolic proteins can be regulated by the methylation status of the corresponding genes, which led to the metabolic heterogeneity among the proteomic subgroups. Fig. S6 (Additional file 2) demonstrated the overlaps of results from RNA-seq, methylation profiling and proteome analysis.
In summary, we proved the heterogeneity of HCC noncancerous tissues was caused by dysregulation of metabolism, which was potentially modified by epigenetic alterations and significantly correlated different survival outcomes. These findings revealed the potential of using metabolic inhibitors to prevent early recurrence of HCC after surgery.
Availability of data and materials
All the datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Tandem mass tag
Gene Set Enrichment Analysis
Differentially methylated region
Bernards R, Weinberg RA. Metastasis genes: a progression puzzle. Nature. 2002;418(6900):823.
Fidler IJ. The organ microenvironment and cancer metastasis. Differentiation. 2002;70(9):498–505.
Budhu A, Forgues M, Ye QH, Jia HL, He P, Zanetti KA, et al. Prediction of venous metastases, recurrence, and prognosis in hepatocellular carcinoma based on a unique immune response signature of the liver microenvironment. Cancer Cell. 2006;10(2):99–111.
Williamson T, Sultanpuram N, Sendi H. The role of liver microenvironment in hepatic metastasis. Clin Transl Med. 2019;8(1):21.
Gao Q, Zhu H, Dong L, Shi W, Chen R, Song Z, et al. Integrated proteogenomic characterization of hbv-related hepatocellular carcinoma. Cell. 2019;179(2):561-77.e22.
Brunner SF, Roberts ND, Wylie LA, Moore L, Aitken SJ, Davies SE, et al. Somatic mutations and clonal dynamics in healthy and cirrhotic human liver. Nature. 2019;574(7779):538–42.
Jewer M, Findlay SD, Postovit L-M. Post-transcriptional regulation in cancer progression : microenvironmental control of alternative splicing and translation. J Cell Commun Signal. 2012;6(4):233–48.
Arakaki TL, Pezza JA, Cronin MA, Hopkins CE, Zimmer DB, Tolan DR, et al. Structure of human brain fructose 1,6-(bis)phosphate aldolase: linking isozyme structure with function. Protein Sci. 2004;13(12):3077–84.
Xiao Z, Dai Z, Locasale JW. Metabolic landscape of the tumor microenvironment at single cell resolution. Nat Commun. 2019;10(1):3763.
Umemura A, He F, Taniguchi K, Nakagawa H, Yamachika S, Font-Burgada J, et al. p62, upregulated during preneoplasia, induces hepatocellular carcinogenesis by maintaining survival of stressed HCC-initiating cells. Cancer Cell. 2016;29(6):935–48.
We are most grateful for West China Biobanks, Department of Clinical Research Management, West China Hospital, Sichuan University for their support of human tissue samples. We thank Professor Jiajie Peng and Core Facility of West China Hospital for their technique support.
This work was supported by grants from the National Key Technologies R&D Program (2018YFC1106800), the Natural Science Foundation of China (82002572, 82002967, 81972747, and 81872004), the fellowship of China National Postdoctoral Program for Innovation Talents (BX20200225, BX20200227), the 1.3.5 project for disciplines of excellence, West China Hospital, Sichuan University (ZYJC18008).
Ethics approval and consent to participate
Informed consensus was obtained from all patients. This study was authorized by The Ethics Committee of West China Hospital.
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.
About this article
Cite this article
Liao, H., Du, J., Wang, H. et al. Integrated proteogenomic analysis revealed the metabolic heterogeneity in noncancerous liver tissues of patients with hepatocellular carcinoma. J Hematol Oncol 14, 205 (2021). https://doi.org/10.1186/s13045-021-01195-y
- Hepatocellular carcinoma
- Liver microenvironment