Abstract

DNA Primase Subunit 1 (PRIM1) is crucial for cancer development and progression. However, there remains a lack of comprehension concerning the clinical implication of PRIM1 in HCC. Here, aberrant expression of PRIM1 was identified in HCC according to available databases. The prognostic value of PRIM1 in patients presenting with HCC was further assessed based on TCGA data. Gene set enrichment analysis (GSEA) was subsequently conducted to investigate the potential function of PRIM1. Additionally, the correlations between tumor-infiltrating immune cells (TIICs) and PRIM1 expression were evaluated. The data from TCGA, GEO, ONCOMINE, and HCCDB databases illustrated that PRIM1 was overexpressed in HCC tissues, compared to normal liver tissues (all ). Kaplan-Meier analysis revealed that high PRIM1 expression in HCC was closely correlated with worse overall survival (). The univariate and multivariate analyses illustrated that PRIM1 expression was an independent novel prognostic indicator in HCC. Additionally, the area under the receiver operating characteristic (AUROC) curve for PRIM1 reached 0.8651, indicating the diagnostic significance of PRIM1 in patients with HCC. GSEA showed that PRIM1 overexpression was significantly enriched in several tumor-related signaling pathways. Besides, TIIC analysis clarified the association between PRIM1 expression and TIICs in HCC. The findings disclose that PRIM1 profoundly implicated in promoting tumorigenesis might work as a desirable biomarker for HCC.

1. Introduction

Hepatocellular carcinoma (HCC) with high mortality rates has become a threat to public health [1]. Although surgery is helpful for prolonging the lifetime of HCC patients, most patients have missed the opportunity for surgical treatment at the time of diagnosis [2]. Patients diagnosed with advanced liver cancer lack curative therapies and have a worse prognosis [3]. Great progress has been made in HCC treatment protocols, including radiofrequency ablation, curative resection, radioembolization, liver transplantation, and systemic targeted therapy [46]. Nonetheless, it remained a huge challenge to discover novel molecular biomarkers for HCC [7].

Primase, a heterodimer of two subunits, is a pivotal enzymatic element in the replication of DNA. DNA primase synthesizes small RNA primers for short DNA fragments generated throughout DNA replication [811]. For this reason, DNA primase, frequently referred to as RNA primase, plays an essential role in RNA polymer synthesis. Primase Subunit 1 (PRIM1), encoding a small, 49 kDa primase subunit, has the catalytic function of an enzyme. Previous studies have implied that PRIM1 closely related to osteosarcoma [12], pancreatic cancer [13], and breast cancer [14] is profoundly implicated in cancer progression [11]. However, the significance of PRIM1 in HCC remains enigmatic.

Recently, a massive number of studies have identified genes associated with the prognosis of HCC [1518]. Until recently, there is still a lack of the accurate and effective markers for the diagnosis and prognostic of HCC [19]. Further studies are still essential to improve the accuracy in diagnosis and prognosis of HCC.

To deepen the comprehension of PRIM1 in HCC, we explored the aberrant expression of PRIM1 based on The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) databases. Additionally, our study also explored the clinical application value of PRIM1 in HCC by Kaplan-Meier survival and Cox regression analyses. As far as we know, we illustrate the role of PRIM1 in HCC for the first time.

2. Material and Methods

2.1. Data Extraction from GEO and TCGA Databases

The RNA-Seq gene expression level 3 Htseq-count data restricted to 371 primary tumor samples and 49 normal samples were selected and downloaded from The Cancer Genome Atlas (TCGA, https://portal.gdc.cancer.gov/) [20]. In this study, we analyzed the pathological data () of HCC collected from TCGA database. Ten patients with histology of fibrolamellar carcinoma () and hepatocholangiocarcinoma (), one patient without survival data, and six patients without gene expression levels were excluded. The study comprised 360 samples of histologically confirmed HCC. Subsequently, the correlations between PRIM1 expression and clinicopathological characteristics of HCC patients were assessed.

GSE25097, GSE6764, GSE14520, GSE45436, GSE55092, and GSE60502 datasets were selected and downloaded from the Gene Expression Omnibus (GEO, https://www.ncbi.nlm.nih.gov/geo/) database [21].

2.2. ONCOMINE Database Analysis

The mRNA levels of PRIM1 in HCC were further elucidated based on the ONCOMINE database (https://www.oncomine.org/), which is a cancer microarray database holding a total of 86,733 samples and 715 gene expression datasets [22]. Genes were screened by fold and , and the top 10% of ranked genes were sorted and considered.

2.3. Integrative Molecular Database of Hepatocellular Carcinoma (HCCDB) Analysis

The Integrative Molecular Database of Hepatocellular Carcinoma (HCCDB, http://lifeome.net/database/hccdb/) database was further applied to annotate the expression patterns of PRIM1 in HCC, which is an integrated database holding a total of 3,917 samples and 15 gene expression datasets [23].

2.4. The Kaplan-Meier Plotter Survival Analysis

Prognostic values of PRIM1 in HCC were further evaluated using the Kaplan-Meier plotter (http://kmplot.com/analysis/) [24]. The log-rank was calculated to display differential survival.

2.5. Gene Set Enrichment Analysis

Gene set enrichment analysis (GSEA, http://software.broadinstitute.org/gsea/) is employed to identify remarkably overrepresented or underrepresented groups of genes [2527]. In the current study, GSEA was performed in terms of the correlation between gene sets and PRIM1 expression, and gene set permutations were carried out 1000 times. Gene sets and the correlative pathways are filtered and ordered by the nominal value and normalized enrichment score (NES).

2.6. Tumor Immune Estimation Resource (TIMER) Database Analysis

Tumor Immune Estimation Resource (TIMER, https://cistrome.shinyapps.io/timer/) database [28] was applied to evaluate the associations between immune cells (B cells, CD4+ T cells, CD8+ T cells, neutrophils, macrophages, and dendritic cells) and PRIM1 expression.

2.7. Statistical Analysis

R software (v.3.6.0; The R Foundation) and GraphPad Prism version 8.0 software (GraphPad Software, Inc.) were applied for statistical analysis and scientific graphing. The correlations of PRIM1 expression and clinicopathological characteristics were evaluated by the Wilcoxon signed-rank test and logistic regression. The Cox regression and Kaplan-Meier analyses were further carried out to elucidate the association between PRIM1 expression and survival data. Furthermore, a receiver operating characteristic curve is plotted to illustrate the underlying diagnostic ability of PRIM1 in HCC. was regarded as statistical significance.

3. Results

3.1. Patient Characteristics

Level 3 mRNA expression and clinical data from 360 primary HCC and 49 normal control tissues were obtained from TCGA database. Data on the clinicopathological features of HCC, such as sex, age, stage, body mass index (BMI), grade, alpha-fetoprotein (AFP) level, platelet, and race, was collected. As listed in Table 1, the patients, median age 60 years (range 16-90), included 243 males and 117 females. Among them, 44.3% of patients were from Asian populations, and other non-Asian populations included White (50.3%), Black or African American (4.9%), and Alaska Native or American Indian (0.6%) populations. There were 6.4% patients with <18.5 kg/m2 body mass index (BMI), 46.5% patients with 18.5-24.99 kg/m2 BMI, 26.6% patients with 25-29.99 kg/m2 BMI, and 20.5% patients with >30 kg/m2 BMI. A total of 167 patients (49.7%) had stage I disease, 81 patients (24.1%) had stage II disease, 84 patients (25.0%) had stage III disease, and only four patients had stage IV disease. Furthermore, there were 14.9% patients with G1 grade, 48.2% of patients with G2 grade, 33.8% of patients with G3 grade, and 3.1% patients with G4 grade. Moreover, this study included 23.3% of patients with AFP and 76.7% of patients with AFP . Additionally, the median platelet was (range ).

3.2. The mRNA Level of PRIM1 Is Upregulated in HCC

As presented in Figure 1(a), PRIM1 overexpression was observed in HCC tissues () based on the gene expression data from TCGA database, compared with that in normal liver tissues (). Moreover, PRIM1 overexpression exhibited in HCC tissues (), compared with that in adjacent liver tissues () (Figure 1(b)). A higher expression of PRIM1 exhibited in HCC tissues than that in the adjacent ones based on six datasets from the GEO database (GSE25097, GSE6764, GSE14520, GSE45436, GSE55092, and GSE60502) (Figures 1(c)1(h)).

>As revealed by an analysis of the Integrative Molecular Database of Hepatocellular Carcinoma (HCCDB) (Figure 2(a)) and ONCOMINE databases (Figure 2(b)), the PRIM1 expression was remarkably enhanced in HCC tissues.

In a word, these data revealed that the mRNA level of PRIM1 was upregulated in HCC, and these results suggested that enhanced expression of PRIM1 might be closely linked with HCC pathogenesis.

3.3. Relationship between the Expression of PRIM1 and Clinicopathological Characteristics in HCC Patients

To further explore the clinical relevance of PRIM1 mRNA expression in HCC, the relationship between PRIM1 expression and the clinicopathological characteristics, such as sex, age, stage, body mass index (BMI), grade, alpha-fetoprotein (AFP) level, platelet, and race, was examined (Table 2). Univariate logistic regression analysis suggested that PRIM1 expression, a categorical dependent variable, was strongly associated with a poor prognosis. A high PRIM1 level in HCC was significantly associated with stage ( [1.07-2.89] for I-II vs. III-IV; ), grade ( [1.47-3.55] for G1-G2 vs. G3-G4; ), and AFP level ( [1.37-4.49] for ≤400 ng/ml vs. >400 ng/ml; ). Additionally, no association between the expression of PRIM1 and other clinicopathological characteristics, such as age (), sex (), BMI (), platelet (), and race (), was found. These results implied that a more advanced stage, high grade, and high AFP level seemed to be attributable to elevated PRIM1 expression.

3.4. The Relationship between PRIM1 Expression and OS in HCC Patients

As displayed in Figure 3(a), Kaplan-Meier analysis was executed based on the data from TCGA database using survival packages in R, which revealed that the overexpression of PRIM1 had an unfavorable OS in patients with HCC (log-rank ). As shown in Figure 3(b), subsequent analysis based on the Kaplan-Meier plotter database was consistent with this result. Moreover, subgroup analysis found that PRIM1 overexpression might be considered a risk factor for the 1-year ( (1.88-6.27), log-rank ), 3-year ( (1.64-3.74), log-rank ), and 5-year ( (1.31-2.72), log-rank ) OS in patients with HCC (Figures 3(c)3(e)).

Additionally, subgroup survival analysis was further executed in various patient populations. PRIM1 overexpression had a remarkable correlation with reduced OS in HCC patients without hepatitis virus infection ( (1.42-3.61), log-rank , Figure 4(a)). Moreover, enhanced PRIM1 expression had a remarkable association with reduced OS in males ( (1.19-2.92), log-rank , Figure 4(b)). Furthermore, PRIM1 overexpression significantly contributed to the poor OS in Asian HCC ( (1.2-4.2), log-rank , Figure 4(c)). Additionally, PRIM1 expression was a risk factor for OS in HCC patients without alcohol consumption ( (1.14-2.91), log-rank , Figure 4(d)).

A high PRIM1 level significantly contributed to worse OS in HCC patients with grade II ( (1-2.8), log-rank ) in Figure 5(a). Additionally, an enhanced expression of PRIM1 had a remarkable association with decreased OS in stage I-II HCC patients ( (1.14-3.03), log-rank , Figure 5(b)) and stage II-III patients ( (1.2-3.12), log-rank , Figure 5(c)).

3.5. PRIM1 Is an Unfavorable Prognostic Factor in HCC

To further explore the prognostic role of PRIM1 in HCC, univariate and multivariate analyses were conducted (Table 3). The univariate analysis showed that elevated PRIM1 expression had a remarkable association with reduced survival (HR: 1.32; 95% CI: 1.06–1.65; ). Moreover, PRIM1 expression was independently associated with overall survival according to multivariate analysis (HR: 1.31; 95% CI: 1.06–1.62; ). Additionally, we found that Asian populations had worse survival compared to that in the non-Asian people (). Overall, these results indicated that elevated PRIM1 expression might have a crucial role in HCC occurrence and development.

To further investigate the potential value of PRIM1 in HCC, the area under the receiver operating characteristic (AUROC) curve for PRIM1 was calculated. As displayed in Figure 6, the value of AUROC reached 0.8651 (), and the results implied that PRIM1 might work as a novel and promising molecular biomarker for HCC.

3.6. GSEA Identifies PRIM1-Related Signaling Pathways

Gene set enrichment analysis (GSEA) was employed to distinguish significant differences between low and high PRIM1 expression datasets. Twelve signaling pathways were significantly enriched according to the normalized enrichment score (NES) (NOM , FDR ). As shown in Figure 7 and Table 4, cell cycle, oocyte meiosis, p53 signaling pathway, Wnt signaling pathway, mTOR signaling pathway, ERBB signaling pathway, phosphatidylinositol signaling system, notch signaling pathway, and RIG-I-like receptor signaling pathway were differentially enriched in PRIM1 high expression phenotype, indicating that PRIM1 expression was closely related to cell growth and death, signal transduction, and immune system.

3.7. Correlation between Tumor-Infiltrating Immune Cells (TIICs) and PRIM1 Expression in HCC

As shown in Figure 8, PRIM1 expression had a positive relationship with the number of B cells (partial correlation, 0.267; ), CD8+ T cells (partial correlation, 0.259; ), neutrophils (partial correlation, 0.181; ), macrophages (partial correlation, 0.199; ), and dendritic cells (partial correlation, 0.305; ), indicating that PRIM1 expression might be closely linked with immunotherapy.

4. Discussion

In the past few years, although therapies of hepatocellular carcinoma (HCC) had achieved significant progress, HCC remained to have high mortality and morbidity rates due to misdiagnosis or delayed diagnosis [29, 30]. It is imperative to identify underlying biomarkers for HCC. Previous results have shown that a particular relationship appears between gene expression and accurate diagnosis. However, it remains a huge challenge to identify a valid biomarker [3133]. As far as we know, the potential prognostic value of PRIM1 in HCC is still confused. Moreover, we describe the features of PRIM1 in HCC based on The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases for the first time, and we centered attention on the underlying impact of PRIM1 on HCC.

High PRIM1 expression was connected with poorly differentiated tumors and poorer survival outcomes in breast cancer [14]. In this study, we initially discovered that PRIM1 expression was significantly increased in HCC using TCGA and GEO databases, which is entirely consistent with the data from the Integrative Molecular Database of Hepatocellular Carcinoma (HCCDB) and ONCOMINE databases. Moreover, we found that HCC patients with high PRIM1 expression showed worse outcomes and lower overall survival. Furthermore, logistic regression showed that PRIM1 expression as a dependent variable was closely related to the poor prognosis. Additionally, multivariate Cox analysis implied that PRIM1 expression was an independent novel prognostic factor for HCC. Collectively, our results indicated that PRIM1 appears to be a desirable prognostic molecular biomarker for HCC. Further studies will be needed to verify the role of PRIM1 in HCC.

Primase Subunit 1 (PRIM1) encodes DNA primase small subunit, which has the catalytic function as enzymes. DNA primase small subunit complexed with polymerase α exerted an import part in DNA replication, indicating that the polα-primase complex might be an underlying target [9, 10]. Previous results have shown that it was necessary to change the conformation of the primer for the initiation and elongation of RNA synthesis [14]. PRIM1 catalyzed RNA synthesis, and it was crucial for the accumulation and stimulation of the DNA damage response. Mutations of PRIM1 induced apoptosis through the ATM-Chk2-p53 pathway [14, 34]. Moreover, PRIM1 might have a synthetically lethal relationship with ATR, indicating that targeting the polα-primase complex might be an effective treatment strategy [13]. However, the underlying molecular mechanism of PRIM1 in HCC remains enigmatic.

A preliminary gene set enrichment analysis (GSEA) was employed to identify PRIM1-related oncogenic pathways, indicating that PRIM1 expression was closely related to cell growth and death, signal transduction, and immune system. Furthermore, the investigation on the correlation between tumor-infiltrating immune cells (TIICs) and PRIM1 expression in HCC showed that PRIM1 expression was positively correlated with the number of B cells, neutrophils, macrophages, and dendritic cells, indicating that PRIM1 expression might be a potential biomarker for immunotherapy in HCC. More studies are needed to elucidate its role in HCC.

5. Conclusions

Our results indicated that PRIM1 might be a desirable molecular biomarker in HCC. Moreover, PRIM1 might be a conceivable therapeutic target for HCC.

Data Availability

Previously reported gene expression data (GSE25097, GSE6764, GSE14520, GSE45436, GSE55092, and GSE60502) were applied to support this study and are available at the GEO database (https://www.ncbi.nlm.nih.gov/geo). TCGA data were obtained from https://portal.gdc.cancer.gov/.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this article.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (No. 81660755) and the Science and Technology Project of Shenzhen of China (No. JCYJ20170307160524377).