Journal of Oncology

Journal of Oncology / 2020 / Article

Research Article | Open Access

Volume 2020 |Article ID 8291036 | https://doi.org/10.1155/2020/8291036

Zhuomao Mo, Shuqiao Zhang, Shijun Zhang, "A Novel Signature Based on mTORC1 Pathway in Hepatocellular Carcinoma", Journal of Oncology, vol. 2020, Article ID 8291036, 11 pages, 2020. https://doi.org/10.1155/2020/8291036

A Novel Signature Based on mTORC1 Pathway in Hepatocellular Carcinoma

Academic Editor: San Lin You
Received04 Jun 2020
Revised15 Aug 2020
Accepted01 Sep 2020
Published15 Sep 2020

Abstract

Background. mTORC1 signal pathway plays a role in the initiation and progression of hepatocellular carcinoma (HCC), but no relevant gene signature was developed. This research aimed to explore the potential correlation between the mTORC1 signal pathway and HCC and establish the related gene signature. Methods. HCC cases were retrieved from The Cancer Genome Atlas (TCGA), International Cancer Genome Consortium (ICGC), and Gene Expression Omnibus (GEO) databases. The genes included in mTORC1-associated signature were selected by performing univariate and multivariate Cox regression analyses and lasso regression analysis. The protein expression level of included genes was verified by The Human Protein Altas. Then, the signature was verified by survival analysis and multiple receiver operating characteristic (ROC) curve. Moreover, the correlation between signature and immune cells infiltration was investigated. Furthermore, a nomogram was established and evaluated by C-index and calibration plot. Results. The signature was established with the six genes (ETF1, GSR, SKAP2, HSPD1, CACYBP, and PNP). Three genes (ETF1, GSR, and HSPD1) have verified their protein expression level in HCC. Under the grouping from signature, patients in the high-risk group showed worse survival than those in the low-risk group in both three datasets. The signature was found to be significantly associated with the infiltration of B cells, CD4+ T-cells, CD8+ T-cells, dendritic cells, macrophages, and neutrophils. The univariate and multivariate Cox regression analysis indicated that mTORC1-related signature could be the potential independent prognostic factor in HCC. Finally, the nomogram involving age, gender, stage, and signature has been established and verified. Conclusion. The mTORC1-associated gene signature established and validated in our research could be used as a potential prognostic factor in HCC.

1. Introduction

Hepatocellular carcinoma (HCC) is one of the most prevalent cancers around the world, becoming the second leading causes of tumor-related death [1]. Owing to the high rate of metastasis, HCC patients with advanced stage are usually with a poor prognosis [2]. Although the treatment and biomarkers of HCC have developed, the clinical outcomes of HCC patients are still unsatisfactory [3]. The occurrent and development of HCC involved interactions between genetics, epigenetics, and transcriptomic alterations [4]. Many studies have verified that different biomarkers have certain prognostic value in HCC. For example, Gu et al. found that CCL14 was a potential prognostic biomarker that correlated with tumor immune cell infiltration in HCC [5]. Another study found the strong correlations between PRPF3 expression and prognosis in HCC [6]. However, as a biomarker, a single gene usually has a lower prognostic value than multigene prognostic signature. Therefore, many gene-related signatures have been developed for predicting prognosis of HCC. For instance, Zhang et al. established a gene signature associated with HCC microenvironment and successfully verified it [7]. Other predictive signatures based on immune [8] and glycolysis [9] also play an important role in HCC prognosis.

Generally, the gene researches usually focus on comparing the gene expression between two groups or pay attention to the highly upregulated and downregulated genes. Nevertheless, some genes which showed no significant difference but had important biological function and characteristics were omitted. In view of this, a computational method gene set enrichment analysis (GSEA) determined whether a prior defined set of genes shows statistically significant differences between two biological states [10]. The advantage of GSEA is that it can identify the genes in which expression is based on the trend of overall level. Consequently, in this research, we identified the pathway and gene with GSEA. Then, we constructed the signature based on related genes and verified it, providing the more comprehensive and accurate prognostic model for clinic.

2. Materials and Methods

2.1. Data Collection

The gene expression data with the type of level 3 RNA-seq FPKM dataset and the clinical messages in TCGA website (https://portal.gdc.cancer.gov/) were retrieved. A total of 377 HCC cases have been downloaded and analyzed. As the validated cohorts, GSE76427 datasets were retrieved from Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/), while ICGC-LIRI dataset was retrieved from International Cancer Genome Consortium (ICGC) database (https://icgc.org/).

2.2. Identification of Pathways and Related Genes

After data collection, we extracted the clinical details and generated the expression matrix of all genes in HCC of TCGA dataset. Then, we divided the patients into two groups according to the survival status and employed GSEA to choose the most relevant pathways. The pathways were considered for further analyses if the normalized value was <0.01. After that, we collected the related genes of involved pathway from Molecular Signatures Database (http://software.broadinstitute.org/gsea/msigdb/index.jsp).

2.3. Construction and Verification of Signature

Firstly, we performed the differentially expressed analysis to select the related genes. The “limma” package under R studio software was employed, and a value <0.05 was considered statistically significant. Secondly, we employed the univariate and multivariate Cox regression analysis to choose the prognostic genes among the differentially expressed genes. The genes in this section were eligible for further selection if a value was <0.05. Then, the lasso regression analysis was executed for checking selected genes. In this analysis, a lasso penalty was applied, to simultaneously account for shrinkage and variable selection. The optimal value of the lambda penalty parameter was defined by performing 10 cross-validations. Using the “glmnet” package, the coefficient of each included genes and risk score of each case were calculated. The calculation formula of risk score was as follows:

The cases were divided into two groups (high risk or low risk), according to the risk score median. To explore the time-dependent prognostic value of our gene signature, the survival analysis was performed using the “survival” package in the R studio software. The relationship between signature and other clinical messages was also evaluated and visualized with a heatmap. The protein expression level of included genes of signature was verified by The Human Protein Altas (https://www.proteinatlas.org/). Besides, the multiple receiver operating characteristic (ROC) curve was performed to check the predictive accuracy of risk score. In addition, we investigated the correlation between signature and six different immune cells (B cells, CD4+ T-cells, CD8+ T-cells, dendritic cells, macrophages, and neutrophils). The infiltration data of six immune cells were retrieved from tumor immune estimation resource (https://cistrome.shinyapps.io/timer/). Moreover, the univariate and multivariate Cox regression analyses were performed to verify whether the risk score is an independent prognostic factor.

2.4. Predictive Nomogram Design

A predictive nomogram based on age, gender, stage, and risk score was constructed using the “rms” package and Cox regression model to predict the overall survival (OS) at 1 year, 3 years, and 5 years of HCC patients. Then, we used Harrell’s concordance index (C-index) and calibration plot to evaluate the nomogram.

3. Results

The clinical data details of the patients used in this study are shown in Table 1. Figure 1 shows the screening process and validation of our study. The result of GSEA in Figure 2 showed that a total of 6 pathways were eligible, and the details of pathways are summarized in Table 2. Considering the highest normalized enrichment score of mTORC1 pathways, we chose the 200 relevant genes in this pathway for further analysis. A total of 199 genes have been found in gene expression matrix, and the differentially expressed analysis showed that 160 genes were significantly different in HCC (see Supplementary files 13). After that, the univariate and multivariate Cox regression analyses demonstrated that 15 genes (ETF1, CTSC, GSR, HSPE1, SKAP2, HSPD1, TES, TFRC, ASNS, EPRS, CANX, CACYBP, UNG, TBK1, and PNP) were included with (see Supplementary files 4 and 5). As illustrated in Figures 3(a) and 3(b), the results of lasso regression analysis further confirmed the signature composed of 6 genes (ETF1, GSR, SKAP2, HSPD1, CACYBP, and PNP). And the coefficients of ETF1, GSR, SKAP2, HSPD1, CACYBP, and PNP were 0.03402, 0.00670 0.02556, 0.00181, 0.02034, and 0.00916, respectively.


Clinical characteristicsNumberPercent (%)

TCGA-LIHC (n=377)
Survival statusSurvival24966
Death12834
Age (1 patient missing)≤65 years23562.5
>65 years14137.5
GenderFemale12268
Male25532
Stage (24 patients missing)I17550
II8724.6
III8624.4
IV51
Grade (5 patients missing)G15514
G218048
G312433
G4135
T classification (3 patients missing)T118549
T29526
T38122
T4133

GSE76427 (n=115)
Survival statusSurvival9280
Death2320
Age≤65 years6556.5
>65 years5043.5
GenderFemale2219.1
Male9380.9
StageI5547.8
II3530.4
III2118.3
IV43.5

ICGC-LIRI (n=260)
Survival statusSurvival21482.4
Death4617.6
Age≤65 years9837.7
>65 years16262.3
StageI4015.4
II11745
III8030.8
IV238.8


NameSizeESNESNormalized value

MTORC1 SIGNALING2000.5816941.9822140
UV_RESPONSE_UP1580.5063781.9552580
GLYCOLYSIS2000.5066721.883370
G2M_CHECKPOINT1990.730771.8818150.002088
MYC_TARGETS_V11990.6858111.8604980.004141
E2F_TARGETS2000.7450061.8568940.002079

ES = enrichment score; NES = normalized enrichment score.

Besides, the heatmap of risk score and clinical parameters was shown in Figure 3(c). All the genes included in the signature were highly expressed in the high-risk group and lowly expressed in the low-risk group. The correlation between gene expression of final included genes and clinical parameters is shown in Table 3. Meanwhile, significant difference was found between risk score and stage, grade, T, and M, respectively (Table 3). In terms of protein expression, three genes (ETF1, GSR, and HSPD1) have been found to be highly expressed in HCC tissue according to Figure 4. The survival analysis of TCGA (Figure 5(a)), GSE76427 (Figure 5(b)), and ICGC (Figure 5(c)) both indicated that significant difference was found between two groups (). And the multiple ROC curve plot (Figures 5(d)5(f)) demonstrated that the risk score got the highest predictability among analyzed factors in 1 year (AUC = 0.802), 3 years (AUC = 0.743), and 5 years (AUC = 0.719). In terms of immune cells infiltration, significant difference was found between risk score and B cells, CD4+ T-cell, CD8+ T-cell, dendrites, macrophage, and neutrophil (Figure 6). Furthermore, we performed the univariate and multivariate Cox regression analyses, and the results in Table 4 revealed that the signature can be the independent prognostic factor in HCC ().


IdAgeGenderGradeStageTMN

ETF10.191 (0.849)0.209 (0.834)−1.899 (0.059)−1.797 (0.075)−1.298 (0.197)0.547 (0.636)−1.291 (0.285)
GSR0.782 (0.436)−1.263 (0.208)−2.216 (0.028)−1.886 (0.062)−2.008 (0.047)2.607 (0.093)2.35 (0.049)
SKAP2−0.721 (0.473)1.488 (0.139)1.274 (0.204)−2.457 (0.016)−2.009 (0.047)1.576 (0.241)−0.815 (0.474)
HSPD10.988 (0.325)−0.442 (0.659)−2.609 (0.010)−1.418 (0.158)−1.17 (0.244)3.016 (0.067)−1.426 (0.244)
CACYBP2.79 (0.006)−0.636 (0.526)−2.635 (0.009)−1.547 (0.125)−1.71 (0.090)0.38 (0.738)−0.015 (0.989)
PNP1.614 (0.109)1.766 (0.080)−1.139 (0.256)−2.781 (0.006)−2.819 (0.006)1.759 (0.202)−0.813 (0.475)
Risk score0.94 (0.054)−0.45 (0.654)−2.687 (0.008)−2.603 (0.010)−2.428 (0.017)4.473 (0.009)−1.171 (0.317)

The results in table represent correlation coefficient and value; the bold represents significantly different values.

valueHazard ratio

Univariate Cox regression analysis
Age0.5911.005 (0.987, 1.023)
Gender0.3010.780 (0.487, 1.249)
Grade0.9141.017 (0.746, 1.387)
Stage<0.0011.865 (1.456, 2.388)
T<0.0011.804 (1.434, 2.270)
M0.0233.850 (1.207, 12.281)
N0.3282.022 (0.494, 8.276)
Risk score<0.0011.205 (1.142, 1.272)

Multivariate Cox regression analysis
Age0.0641.019 (0.999, 1.040)
Gender0.8700.957 (0.567, 1.617)
Grade0.7070.935 (0.659, 1.327)
Stage0.8960.935 (0.340, 2.569)
T0.2231.752 (0.711, 4.314)
M0.3171.990 (0.518, 7.654)
N0.5951.686 (0.246, 11.564)
Risk score<0.0011.236 (1.152, 1.328)

Finally, using the TCGA cohort, we built the nomogram based on age, gender, stage, and established signature to predict 1-year, 3-year, and 5-year OS for HCC patients (Figure 7(a)). The results of decision curve analysis (Supplementary file 6) demonstrated that nomogram with age, gender, stage, and established signature showed a higher benefit than other two solo models (stage only or risk score only). The C-index of nomogram was 0.73, and the calibration plot for the probability of survival at 1 year (Figure 7(b)), 3 years (Figure 6(c)), and 5 years (Figure 6(d)) showed good agreement between the prediction by nomogram and real observation. Also, we established the nomograms and calibration plots based on GEO and ICGC cohorts (Supplementary files 7 and 8), and their C-index was 0.70 and 0.76, respectively.

4. Discussion

In the initiation and progression of hepatocellular carcinoma, genetic factors usually play an important role. Meanwhile, mRNA gene signature based on a certain characteristic like glycolysis [11] and immune [12] has been developed for predicting cancer prognosis. In this research, we explored specific function to identify genes by GSEA that could predict the survival of HCC patients. According to our results, six signal pathways were found to be highly related to survival and we established the gene signature with the mTORC1 signal pathway. As we all know, the mTOR pathway is a serine/threonine protein kinase belonging to the PI3K-related kinase family [13], which comprised of two distinct complexes (mTORC1 and mTORC2). With Raptor as its unique and key protein component, mTORC1 plays an important role in cell survival, autophagy, and metabolism [14]. Concerning for the mTOR signal pathway in HCC, it has been found that aberrant mTOR signaling was present in half of the HCC cases [15]. Meanwhile, an intact mTORC1 axis [16] and mTORC2-Akt1 cascade [17] were required for c-Myc-driven hepatocarcinogenesis. Moreover, some research studies [15, 18] provided the theoretical basis of mTOR signaling pathway-oriented targeting treatment for HCC in clinic. Overall, these abovementioned evidences demonstrated that mTOR signal pathway plays an important role in the development of HCC.

In this study, we identified six genes in signature by performing the differentially expressed analysis, univariate Cox regression analysis, and lasso regression analysis. Among our included genes, five genes have been found to be related to HCC from previous studies. Singh et al. found that ETF1, CNOT6, and XRN1 gene in HepG2 cells led to significant alteration in stability of specific mRNAs, and this mechanism may hold novel cancer therapeutic targets [19]. In another research, McLoughlin concluded that GSR, TRXR1, NRF2, and oxidative stress determined hepatocellular carcinoma malignancy [20]. Lee’s et al. study [21] found that HSPD1 was downregulated during early apoptosis of the hepatoma cell mediated by Paeoniae Radix. In terms of CACYBP, it has been verified that CACYBP can promote hepatocellular carcinoma progression in the absence of RNF41-mediated degradation [22]. Moreover, a study [23] found that PNP/fludarabine suicide gene system induced HCC cell apoptosis and inhibited the growth of HCC cells. Although we found no evidence supporting the correlation between SKAP2 and HCC, it has been verified that SKAP2 promotes podosome formation to facilitate tumor-associated macrophage infiltration and metastatic progression [24].

Recently, further investigations have been performed to explore how the mTOR signal transduction mechanisms modulate sensitivity of targeted therapies, angiogenesis, and tumor immunity [25]. The interest in mTOR targeting may improve immune response against cancer and develop new therapeutic strategy. It has been verified that an inflammatory-CCRK circuitry drove mTORC1-dependent metabolic and immunosuppressive reprogramming in obesity-related hepatocellular carcinoma [26]. In another study [27], Tan concluded that Tim-3-mediated PI3K/mTORC1 interference leads to the dysfunction of both tumor-infiltrating conventional natural killer cells and liver-resident natural killer cells. In our research, the results showed that mTORC1 signature significantly associated with B cells, CD4+ T-cell, CD8+ T-cell, dendrites, macrophage, and neutrophil, which indicated that the patients in high-risk group may benefit from immune-targeted therapies and provide a new strategy for immune checkpoint-based targeting.

Being different from the previous prognostic studies in HCC, our predictive model firstly concentrated on the mTORC1 signal pathway. More importantly, the mTORC1 signal pathway was identified by GSEA, which indicated the underlying mechanism between survival of HCC and mTORC1 pathway. Moreover, the validation from three independent datasets and a rigorous screening process enabled the identification of a reliable signature. However, our study has some limitations. First, prognostic signature showed a relatively low diagnostic performance in predicting 5-year OS. It may be attributed to that only 200 associated genes were defined and evaluated for the initiation of the screening process. Second, using a single characteristic (mTORC1 signal pathway) to establish the predictive model is an intrinsic weakness. Indeed, many other mechanisms, such as metabolism [28] and immune [8], have an influence on the development and progression of HCC. Third, the weak correlation between risk score and immune cells infiltration was observed in our research, which may attribute to the calculated method (TIMER) of immune cell infiltration. Furthermore, our signature explored the underlying effect between mTOR signal pathway and HCC, but it is necessary to perform more independent trials and functional experiments to shed light on the mechanism linking them.

5. Conclusion

Our study is the first to identify a novel gene signature related to mTORC1 signal pathway that could be used as a potential prognostic factor in hepatocellular carcinoma.

Abbreviations

HCC:Hepatocellular carcinoma
TCGA:The Cancer Genome Atlas
ICGC:International Cancer Genome Consortium
GEO:Gene Expression Omnibus
ROC:Receiver operating characteristic
GSEA:Gene set enrichment analysis
OS:Overall survival
C-index:Concordance index.

Data Availability

The datasets generated and/or analyzed during the current study are available in the TCGA (https://portal.gdc.cancer.gov/), GEO (https://www.ncbi.nlm.nih.gov/geo/), and ICGC (https://icgc.org/).

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Zhuomao Mo and Shuqiao Zhang contributed equally. Zhuomao Mo and Shijun Zhang designed the manuscript. Zhuomao Mo and Shuqiao Zhang wrote and completed the manuscript. Zhuomao Mo and Shijun Zhang completed the data download and analysis. All the authors approved the final manuscript.

Acknowledgments

This study was funded by the National Natural Science Foundation of China (nos. 81873248 and 81673903).

Supplementary Materials

Supplementary file 1: the gene expression matrix of 199 mTORC1-related genes (TXT 890 kb). Supplementary file 2: the heatmap of 199 mTORC1-related genes (JPG 10489 kb). Supplementary file 3: the volcano plot of 199 mTORC1-related genes (JPG 571 kb). Supplementary file 4: the results of univariate Cox regression analysis for 160 genes (DOC 26.2 kb). Supplementary file 5: the results of multivariate Cox regression analysis for 101 genes (DOC 22.7 kb). Supplementary file 6: the results of decision curve analysis for three different models (JPG 988 kb). Supplementary file 7: the nomogram and calibration plots based on GEO cohort (JPG 608 kb). Supplementary file 8: the nomogram and calibration plots based on ICGC cohort (JPG 653 kb). (Supplementary Materials)

References

  1. J. Ferlay, I. Soerjomataram, R. Dikshit et al., “Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012,” International Journal of Cancer, vol. 136, no. 5, pp. E359–E386, 2015. View at: Publisher Site | Google Scholar
  2. J. Ferlay, H.-R. Shin, F. Bray, D. Forman, C. Mathers, and D. M. Parkin, “Estimates of worldwide burden of cancer in 2008: GLOBOCAN 2008,” International Journal of Cancer, vol. 127, no. 12, pp. 2893–2917, 2010. View at: Publisher Site | Google Scholar
  3. H. B. El-Serag, J. A. Marrero, L. Rudolph, and K. R. Reddy, “Diagnosis and treatment of hepatocellular carcinoma,” Gastroenterology, vol. 134, no. 6, pp. 1752–1763, 2008. View at: Publisher Site | Google Scholar
  4. R. N. Aravalli, E. N. K. Cressman, and C. J. Steer, “Cellular and molecular mechanisms of hepatocellular carcinoma: an update,” Archives of Toxicology, vol. 87, no. 2, pp. 227–247, 2013. View at: Publisher Site | Google Scholar
  5. Y. Gu, X. Li, Y. Bi et al., “CCL14 is a prognostic biomarker and correlates with immune infiltrates in hepatocellular carcinoma,” Aging, vol. 12, no. 1, pp. 784–807, 2020. View at: Publisher Site | Google Scholar
  6. Y. Liu, Y. Yang, Y. Luo et al., “Prognostic potential of PRPF3 in hepatocellular carcinoma,” Aging, vol. 12, no. 1, pp. 912–930, 2020. View at: Publisher Site | Google Scholar
  7. F.-P. Zhang, Y.-P. Huang, W.-X. Luo et al., “Construction of a risk score prognosis model based on hepatocellular carcinoma microenvironment,” World Journal of Gastroenterology, vol. 26, no. 2, pp. 134–153, 2020. View at: Publisher Site | Google Scholar
  8. Y. He, Q. Dang, J. Li et al., “Prediction of hepatocellular carcinoma prognosis based on expression of an immune-related gene set,” Aging, vol. 12, no. 1, pp. 965–977, 2020. View at: Publisher Site | Google Scholar
  9. L. Jiang, L. Zhao, J. Bi et al., “Glycolysis gene expression profilings screen for prognostic risk signature of hepatocellular carcinoma,” Aging, vol. 11, no. 23, pp. 10861–10882, 2019. View at: Publisher Site | Google Scholar
  10. A. Subramanian, P. Tamayo, V. K. Mootha et al., “Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles,” Proceedings of the National Academy of Sciences, vol. 102, no. 43, pp. 15545–15550, 2005. View at: Publisher Site | Google Scholar
  11. L. Mesirov, Z. Zhang, and Z. Yu, “Identification of a novel glycolysis-related gene signature for predicting metastasis and survival in patients with lung adenocarcinoma,” Journal of Translational Medicine, vol. 17, no. 1, p. 423, 2019. View at: Publisher Site | Google Scholar
  12. Y. She, X. Kong, Y. Ge et al., “Immune-related gene signature for predicting the prognosis of head and neck squamous cell carcinoma,” Cancer Cell International, vol. 20, no. 1, p. 22, 2020. View at: Publisher Site | Google Scholar
  13. B. Raught, A.-C. Gingras, and N. Sonenberg, “The target of rapamycin (TOR) proteins,” Proceedings of the National Academy of Sciences, vol. 98, no. 13, pp. 7037–7044, 2001. View at: Publisher Site | Google Scholar
  14. M. Laplante and D. M. Sabatini, “mTOR signaling in growth control and disease,” Cell, vol. 149, no. 2, pp. 274–293, 2012. View at: Publisher Site | Google Scholar
  15. A. Villanueva, D. Y. Chiang, P. Newell et al., “Pivotal role of mTOR signaling in hepatocellular carcinoma,” Gastroenterology, vol. 135, no. 6, pp. 1972–1983, 2008. View at: Publisher Site | Google Scholar
  16. P. Battiston, M. Ge, J. Hu et al., “A functional mammalian target of rapamycin complex 1 signaling is indispensable for c-Myc-driven hepatocarcinogenesis,” Hepatology, vol. 66, no. 1, pp. 167–181, 2017. View at: Publisher Site | Google Scholar
  17. Z. Pes, M. Xu, P. Liu et al., “The mTORC2-Akt1 cascade is crucial for c-Myc to promote hepatocarcinogenesis in mice and humans,” Hepatology, vol. 70, no. 5, pp. 1600–1613, 2019. View at: Publisher Site | Google Scholar
  18. C. Pascale, Z. Ning, A. Wang et al., “USP10 suppresses tumor progression by inhibiting mTOR activation in hepatocellular carcinoma,” Cancer Letters, vol. 436, pp. 139–148, 2018. View at: Publisher Site | Google Scholar
  19. P. Liu, R. S. James, C. J. Mee, and I. Y. Morozov, “mRNA levels are buffered upon knockdown of RNA decay and translation factors via adjustment of transcription rates in human HepG2 cells,” RNA Biology, vol. 16, no. 9, pp. 1147–1155, 2019. View at: Publisher Site | Google Scholar
  20. M. R. McLoughlin, D. J. Orlicky, J. R. Prigge et al., “TrxR1, Gsr, and oxidative stress determine hepatocellular carcinoma malignancy,” Proceedings of the National Academy of Sciences, vol. 116, no. 23, pp. 11408–11417, 2019. View at: Publisher Site | Google Scholar
  21. S. M. Y. Lee, M. L. Y. Li, Y. C. Tse et al., “Paeoniae Radix, a Chinese herbal extract, inhibit hepatoma cells growth by inducing apoptosis in a p53 independent pathway,” Life Sciences, vol. 71, no. 19, pp. 2267–2277, 2002. View at: Publisher Site | Google Scholar
  22. Y.-F. Lian, Y.-L. Huang, Y.-J. Zhang et al., “CACYBP enhances cytoplasmic retention of P27Kip1 to promote hepatocellular carcinoma progression in the absence of RNF41 mediated degradation,” Theranostics, vol. 9, no. 26, pp. 8392–8408, 2019. View at: Publisher Site | Google Scholar
  23. B. Zhang, M. Chen, Y. Zhang, W. Chen, L. Zhang, and L. Chen, “An ultrasonic nanobubble-mediated PNP/fludarabine suicide gene system: a new approach for the treatment of hepatocellular carcinoma,” PLoS One, vol. 13, no. 5, Article ID e0196686, 2018. View at: Publisher Site | Google Scholar
  24. M. Tanaka, S. Shimamura, S. Kuriyama, D. Maeda, A. Goto, and N. Aiba, “SKAP2 promotes podosome formation to facilitate tumor-associated macrophage infiltration and metastatic progression,” Cancer Research, vol. 76, no. 2, pp. 358–369, 2016. View at: Publisher Site | Google Scholar
  25. C. Bazzichetto, F. Conciatori, I. Falcone, and L. Ciuffreda, “Translational landscape of mTOR signaling in integrating cues between cancer and tumor microenvironment,” Advances in Experimental Medicine and Biology, vol. 1223, pp. 69–80, 2020. View at: Publisher Site | Google Scholar
  26. H. Sun, W. Yang, Y. Tian et al., “An inflammatory-CCRK circuitry drives mTORC1-dependent metabolic and immunosuppressive reprogramming in obesity-associated hepatocellular carcinoma,” Nature Communications, vol. 9, no. 1, p. 5214, 2018. View at: Publisher Site | Google Scholar
  27. S. Tan, Y. Xu, Z. Wang et al., “Tim-3 hampers tumor surveillance of liver-resident and conventional NK cells by disrupting PI3K signaling,” Cancer Research, vol. 80, no. 5, pp. 1130–1142, 2020. View at: Publisher Site | Google Scholar
  28. G. M. Liu, W. X. Xie, C. Y. Zhang, and J. W. Xu, “Identification of a four-gene metabolic signature predicting overall survival for hepatocellular carcinoma,” Journal of Cellular Physiology, vol. 235, no. 2, pp. 1624–1636, 2020. View at: Publisher Site | Google Scholar

Copyright © 2020 Zhuomao Mo et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder
Views71
Downloads36
Citations

Related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.