Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2023 Mar 28:13:1144269.
doi: 10.3389/fonc.2023.1144269. eCollection 2023.

Over-expression of RRM2 predicts adverse prognosis correlated with immune infiltrates: A potential biomarker for hepatocellular carcinoma

Affiliations

Over-expression of RRM2 predicts adverse prognosis correlated with immune infiltrates: A potential biomarker for hepatocellular carcinoma

Zhongqiang Qin et al. Front Oncol. .

Abstract

Background: Ribonucleotide reductase regulatory subunit M2 (RRM2) has been reported to be an oncogene in some malignant tumors, such as lung adenocarcinoma, oral squamous cell carcinoma, glioblastoma, and breast cancer. However, the clinical significance of RRM2 in hepatocellular carcinoma has been less studied. The aim of this study was to assess the importance of RRM2 in hepatocellular carcinoma (HCC) based on the Cancer Genome Atlas (TCGA) database.

Methods: The RRM2 expression levels and clinical features were downloaded from the TCGA database. Immunohistochemistry results between tumor tissues and normal tissues were downloaded from the Proteinatlas database. Meanwhile, the expression levels of RRM2 in tumor and paraneoplastic tissues were further verified by qRT-PCR and Western Blotting. Gene Ontology (GO)/Kyoto Encyclopedia of Genes and Genomes (KEGG) and protein-protein-interactions (PPI) network were constructed to analyze RRM2-related downstream molecules. In addition, RRM2 expression-related pathways performed by gene set enrichment analysis (GSEA). Association analysis of RRM2 gene expression and immune infiltration was performed by single-sample GSEA (ssGSEA).

Results: The RRM2 expression level in tumor tissues was higher than normal tissues (P <0.001). The elevated expression of RRM2 in HCC was significantly correlated with T stage (P <0.05), pathologic stage (P <0.05), tumor status (P <0.05), histologic grade (P<0.001), and AFP (P <0.001). HCC with higher RRM2 expression was positively associated with worse OS (overall survival), PFS (progression-free survival), and DSS (disease-specific survival). In the univariate analysis, the expression of RRM2, T stage, M stage, pathologic stage, and tumor status were negatively correlated with OS (P <0.05). Further analysis using multivariate Cox regression showed that tumor status (P<0.01) and RRM2 expression (P<0.05) were independent prognostic factors of OS in HCC. GO/KEGG analysis showed that the critical biological process (chromosome condensation and p53 signaling pathway) might be the possible function mechanism in promoting HCC. Moreover, GSEA showed that several pathways were enriched in RRM2 high-expression samples, including PD-1 signaling, cell cycle, P27 pathway, and T cell receptor signaling pathway. RRM2 was significantly correlated with the infiltration level of CD8 T cells, Cytotoxic cells, DCs, Neutrophils, NK cells, and T helper cells (P <0.05).

Conclusion: Over-expression of RRM2 predict adverse prognosis and is correlated with immune infiltrates in HCC. RRM2 may be a significant molecular biomarker for HCC diagnosis and prognosis.

Keywords: The Cancer Genome Atlas; bioinformatics analysis; biomarker; hepatocellular carcinoma; ribonucleotide reductase regulatory subunit M2.

PubMed Disclaimer

Conflict of interest statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Figures

Figure 1
Figure 1
Relationship between RRM2 expression and hepatocellular carcinoma (HCC). (A) Differential expression of RRM2 between tumor tissues and normal tissues. (B) Differential expression of RRM2 between tumor tissues and matched para-cancerous tissues. (C) The mRNA level of RRM2. (D) The protein level of RRM2 in 5 paired of HCC tissues and adjacent tissues using Western blotting. (E) The results of Immunohistochemistry between the normal tissues and the tumor tissues. (F) The correlation analysis of mRNA expression levels between RRM2 and NCAPG was performed. (G) Diagnostic value of RRM2 expression in HCC. *** represents P <0.001.
Figure 2
Figure 2
Correlation between RRM2 expression and clinical characteristics. Relationship of RRM2 expression with (A) T stage, (B) Pathologic stage, (C) Tumor status, (D) Histologic grade, and (E) AFP level. *represents P <0.05, ** represents P <0.01, *** represents P <0.001.
Figure 3
Figure 3
Increased expression of RRM2 indicated poor prognosis in HCC. Of the 374 cases, patients with high expression of RRM2 had significantly (A) shorter OS (P=0.003), (B) DSS (P=0.003), and (C) PFI (P=0.001); (D) A nomogram that integrate RRM2 and other prognostic factors in HCC from TCGA data; (E) The calibration curve of the nomogram.
Figure 4
Figure 4
Network construction for RRM2 correlated genes in HCC. (A) Top 5 genes of positively or negatively correlated with RRM2 were shown in Heatmap. (B) GO analysis and KEGG pathway reveal the underlying mechanism of RRM2 in the promotion of HCC. (C) The PPI network of RRM2 interaction partners generated by STRING and Cytoscape. The color represents the degree score.
Figure 5
Figure 5
GSEA enrichment analysis results. GSEA results showed that (A) SIGNALING BY RHO FTPASES, (B) GPCR ligand binding, (C) PD1 signaling, (D) Cell cycle, (E) DNA replication, (F) T cell receptor signaling pathway, (G) P27 pathway, (H) RB pathway, (I) SRCRPTP pathway. FDR, false discovery rate; NES, normalized Enrichment Score.
Figure 6
Figure 6
Relationship between RRM2 expression and immune infiltration. (A) Relationship between RRM2 expression and immune cells. (B–M) Further analysis in different immune cell subpopulations in the high- and low-expression groups of RRM2. *represents P <0.05, ** represents P <0.01, *** represents P <0.001. NS, Non Significance.

References

    1. Bresnahan E, Lindblad KE, Ruiz de Galarreta M, Lujambio A. Mouse models of oncoimmunology in hepatocellular carcinoma. Clin Cancer Res (2020) 26(20):5276–86. doi: 10.1158/1078-0432.CCR-19-2923 - DOI - PMC - PubMed
    1. 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. doi: 10.3322/caac.21492 - DOI - PubMed
    1. Huitzil-Melendez FD, Capanu M, O'Reilly EM, Duffy A, Gansukh B, Saltz LL, et al. . Advanced hepatocellular carcinoma: which staging systems best predict prognosis? J Clin Oncol (2010) 28(17):2889–95. doi: 10.1200/JCO.2009.25.9895 - DOI - PMC - PubMed
    1. Amin MB, Greene FL, Edge SB, Compton CC, Gershenwald JE, Brookland RK, et al. . AJCC cancer staging manual. 8th ed. New York: Springer; (2017). - PubMed
    1. Kattan MW, Hess KR, Amin MB, Lu Y, Moons KG, Gershenwald JE, et al. . American Joint committee on cancer acceptance criteria for inclusion of risk models for individualized prognosis in the practice of precision medicine. CA A Cancer J Clin (2016) 66(5):370–4. doi: 10.3322/caac.21339 - DOI - PMC - PubMed