Ziying Zhang1
Peiguo Cao- 1Department of Oncology, Third Xiangya Hospital, Central South University, Changsha, China
- 2Department of Urology, Xiangya Hospital, Central South University, Changsha, China
- 3Department of Thoracic Surgery, Second Xiangya Hospital, Central South University, Changsha, China
Purpose: Our research was aimed to identify the expression, clinical value and biological significance of GINS complex subunit 4 (GINS4) in hepatocellular carcinoma (HCC).
Materials and Methods: GINS4 was initially screened through weighted gene co-expression network analysis (WGCNA). The TCGA, GEO, and TIMER databases were applied for analyzing the GINS4 mRNA expression in HCC. GINS4 protein levels were detected via immunohistochemistry (IHC). Receiver operating characteristic (ROC) curve was applied for estimating the diagnostic significance of GINS4 in HCC. Kaplan-Meier plots, Cox model, and nomogram were used to assess the prognostic performance of GINS4 in HCC. Nomogram validation was conducted through time-dependent ROC and decision curve analysis (DCA). The Wanderer, UALCAN, and DiseaseMeth databases were utilized to identify GINS4 methylation levels in HCC. Genes co-expressed with GINS4 in HCC were estimated through the TCGA, cBioPortal, and GEPIA. GO, KEGG, and GSEA unraveled the possible biological mechanisms of GINS4 in HCC.
Results: WGCNA confirmed that GINS4 was one of hub genes significantly associated with histological grade of HCC. Multiple databases confirmed the significant upregulation of GINS4 in HCC tissues compared with non-tumor controls. IHC analysis of 35 HCC patients demonstrated that overexpressed GINS4 positively correlated with advanced TNM stage and poor pathological differentiation. GINS4 could effectively differentiate HCC cases from healthy individuals, with an AUC of 0.865. Increased GINS4 expression predicted unsatisfactory prognosis in HCC patients, especially in age >60 years, histological grade 1, HBV infection-negative, and occurring relapse subgroup. Nomogram incorporating GINS4 level and TNM stage displayed satisfactory predictive accuracy and clinical utility in predicting HCC prognosis. Upregulated GINS4 exhibited hypomethylated levels in HCC. Functional analysis indicated that GINS4 potentially positively modulated cell cycle and PI3K/AKT/mTOR pathway.
Conclusion: GINS4 is overexpressed in HCC and is correlated with undesirable survival of HCC patients.
Introduction
Liver cancer is the second primary reason for tumor-associated deaths globally, with approximately 841,000 new diagnoses and 782,000 deaths annually (1). Liver cancer kills approximately 383,000 people per year in China, occupying about 51% of liver cancer-related deaths globally (2). HCC, the primary subtype of liver cancer, accounts for 75–85% cases (3, 4). HCC can be triggered by multifarious risk factors, such as chronic hepatitis B virus (HBV) and hepatitis C virus (HCV) infections, aflatoxin exposure, alcoholic abuse, autoimmune hepatitis, and metabolic disorders (3, 5–8). Despite the rapid progression of therapeutic interventions (such as radiofrequency ablation, hepatic resection, hepatic transplantation, transarterial chemoembolization, and stereotactic body radiation), the prognosis of HCC patients is undesirable due to the occurrence of distant metastasis and the increased recurrence (9). Furthermore, the majority of cases are initially diagnosed at advanced HCC owing to the non-specific symptoms in early stage and the deficiency of sensitive diagnostic biomarkers, with a 5-year survival rate of lower than 20% (4, 10). Hence, deep comprehension of the underlying mechanisms concerning HCC progression is required to unravel novel diagnostic and prognostic molecular biomarkers and to develop new effective therapeutic strategies of HCC.
The GINS complex, a heterotetrameric structure composed of four different subunits (Sld5, Psf1, Psf2, and Psf3 from the Japanese go-ichi-ni-san representing 5-1-2-3, also known as GINS4, GINS1, GINS2, and GINS3 in human genome, respectively), can interact with Cdc45 and Mcm2-7 to form the eukaryotic replicative helicase CMG (Cdc45-Mcm helicase-GINS) complex that unties double-stranded DNA prior to moving the replication fork during chromosome duplication (11–13). The GINS complex, without prominent enzymatic activity itself, is pivotal to initiate and elongate chromosome replication through binding to and strengthening the enzymatic function of Mcm helicase (14, 15). GINS4, also known as SLD5, a vital component of GINS complex, exerts a momentous effect on the initiation and prolongation of DNA replication in the G1/S phase cell cycle in eukaryotes (16). GINS4 participates in modulating early embryogenesis in mice and maintaining cell cycle progression and genomic stability in Drosophila (17, 18), indicating its effect on tumorigenesis. Prior studies have demonstrated overexpression of GINS4 in multifarious human cancers tissues and tumor cell lines, including colorectal cancer (CRC) (19, 20), bladder cancer (21), non-small cell lung cancer (NSCLC) (22), gastric cancer (23), and pancreatic cancer (24). Greater expression level of GINS4 in above human tumors is positively correlated with malignant biological properties, such as tumor proliferation, colony forming ability, migration, and invasion as well as epithelial-mesenchymal transition both in vitro and in vivo (19–23). Additionally, survival analysis has revealed that patients with tumor (such as NSCLC, gastric cancer, CRC, and pancreatic cancer) characterized by high GINS4 expression have significantly diminished overall survival (OS) and disease-free survival (DFS) compared with those with lower GINS4 expression (19, 20, 22–24). Thus, GINS4 exerts a vital effect on the malignant progression of tumors and potentially serves as a valuable target for cancer therapy and diagnosis. Nevertheless, no report exists on the role of the GINS4 in HCC so far.
In our report, we investigated the expression, clinical significance, and potential biological functions of GINS4 in HCC based on multiple databases and experiment validation. Initially, the mRNA expression profiles and corresponding clinical information of 371 HCC patients from The Cancer Genome Atlas (TCGA) database and 713 HCC cases from multiple Gene Expression Omnibus (GEO) datasets were analyzed to compare GINS4 mRNA levels between HCC samples and adjacent liver tissues. Meanwhile, GINS4 protein expression was detected through IHC analysis of 35 clinical HCC samples and paired adjacent liver tissues. Secondly, ROC curve evaluated the diagnostic performance of GINS4 and AFP for HCC. The Kaplan-Meier curve, Cox regression models, nomogram, time-dependent ROC curve, and DCA investigated the prognostic performance of GINS4 in HCC. Finally, GINS4 methylation level and its association with clinicopathological factors of HCC were determined via the Wanderer, UALCAN, and human disease methylation (DiseaseMeth) database. Genes co-expressed with GINS4 were identified through TCGA, cBioPortal, and Gene Expression Profiling Interactive Analysis (GEPIA) databases. Multiple bioinformatics analysis methods, including Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Gene Set Enrichment Analysis (GSEA) as well as Pearson correlation analysis, were used to predict the potential mechanism of GINS4 in HCC.
Materials and Methods
Collection of Clinical Specimens
The flowchart of our research was illustrated in Figure 1. From December 2017 and March 2020, a total of 35 paired surgically resected HCC samples and adjacent normal liver specimens were acquired from our hospital, which were used for IHC analysis. All primary HCC individuals had not accepted radiotherapy or chemotherapy before surgery. The present project was approved by the Ethical Committee for Clinical Research of our hospital. All volunteers conferred written informed consent to enroll in this project following receiving a full explanation of the purpose of the project.
Figure 1 Flow chart of the present study. DEGs, differently expressed genes; WGCNA, weighted gene co-expression network analysis; IHC, immunohistochemistry; TIMER, the Tumor Immune Estimation Resource database; AUC, the area under the curve; ROC, the receiver operating characteristic; DCA, Decision curve analysis; GEPIA, Gene Expression Profiling Interactive Analysis.
Public Database-Based Excavation of Data
The mRNA expression profiles extracted from the TCGA database (https://tcga-data.nci.nih.gov/tcga/) were normalized through “DEseq2” package (25). False discovery rate (FDR) <0.05 and |log2 fold change (FC) | > 1.0 were as thresholds for screening differently expressed genes (DEGs) between HCC tissues and adjacent normal liver samples. The heatmap and volcano map of DEGs were constructed by “pheatmap” and “ggplot2” R packages, respectively (5). The corresponding clinicopathologic variables was acquired from the TCGA GDC Data Portal database, including age at diagnosis, gender, ethnicity, body mass index (BMI), AFP level, fibrosis score, Child-Pugh classification, T classification, N classification, M classification, TNM stage, histological grade, tumor status, survival time, and survival status as well as whether had family history, viral hepatitis, hepatic inflammation, residual tumor, vascular invasion, relapse or not. We also used the “limma” package to investigate mRNA expression of DEGs between HCC samples and adjacent noncancer tissues in GEO dataset (http://www.ncbi.nlm.nih.gov/geo/), including GSE14520 (26), GSE25097 (27), GSE54236 (28), and GSE76427 (29).
WGCNA Used for the Screening of GINS4
WGCNA is a novel systematic biology to unravel the association between gene networks modules and clinical phenotype at transcriptome level (30, 31). The dynamic tree cut approach was used to identify module. The modules with high similarity were estimated through cluster analysis. Modules could be merged when their correlation of module eigengene (ME) was higher than 0.95, indicating similar expression profiles among them. Pearson’s correlation analysis assessed the association between MEs and clinicopathological variables, such as gender, AFP level, Child-Pugh classification, TNM stage, histological grade, relapse, and survival status. Generally, the module with the largest absolute of module significance (MS) was selected for subsequent analysis. The gene significance (GS) and module membership (MM) were used to quantify the associations between each gene and external clinical traits in this module (31). Specifically, MM > 0.8 and GS > 0.2 were defined as the thresholds to screen hub genes in above target module. In our report, we utilized the “WGCNA” R package to formulate a co-expression network of 4344 DEGs in 371 HCC patients with corresponding clinical information (32).
Immunohistochemistry
Resected specimens were fixed in 10% formalin, dehydrated, and embedded in paraffin. The paraffin specimens were segmented into 3 μm-thick sections and installed on the glass slide. Initially, the slides were incubated at 60°C for 30 min in a calorstat. Deparaffinization was carried out in xylene and rehydration was then performed in the gradient ethanol. Afterwards, the glass slides were boiled in EDTA solution (pH 8.0) for 5 min to block the endogenous peroxidase activity. Then, the sections were washed with PBS for three times and were then incubated with primary anti-GINS4 antibody (1:100; ab101346, Abcam, UK) at 4°C overnight. A secondary goat anti-rabbit antibody (1:200; ab205718, Abcam, UK) further incubated the slides at 37°C for 30 min. Binding of the primary antibodies was visualized via incubating chromogen diaminobenzidine (DAB, Sigma, UK) for 10 min at 37°C. The sections were counterstained with hematoxylin, dehydrated by a gradient ethanol, followed by xylene, and mounted (20, 23, 33).
The staining of each specimen was evaluated through two independent investigators blinded to the clinicopathological information. The GINS4 expression was considered positive when it was present in the membrane, the cytoplasm, or both. Each specimen was assessed at 200 and 400 magnification. The staining score was assessed according to two parameters: intensity and extension. The percentage of positively stained cells corresponded to five scoring grades: 0, less than 10%; 1, 10 to 25%; 2, 26 to 50%; 3, 51 to 75%; and 4, 76 to 100%. The intensity score was classified as 0, without staining; 1, yellow; 2, yellow-brown; 3, dark brown. The product of intensity and extension was identified as the total staining score which were stratified into three grades: 0 to 3, negative expression; 4 to 6, weakly positive expression; and more than 6, strongly positive expression (20, 23).
Survival Analysis and Establishment of Nomogram
Based on the median GINS4 expression level, a total of 371 HCC cases were stratified into two groups (high versus low expression) in the TCGA database, respectively. The Kaplan-Meier survival curve with Wilcoxon rank sum test was formulated via “survival” R package to evaluate the OS and survival difference between high and low GINS4 expression groups (34). Univariate and multivariate Cox model was formulated for estimating the hazard ratio (HR) with 95% confidence interval (CI). The statistically significant prognostic factors identified by the univariate analysis were further incorporated into the multivariate analysis.
Nomogram is well-acknowledged model to predict long-term prognosis of patients with tumor (35, 36). The “rms” R package was applied for building a prognostic nomogram, thus estimating the probability of the 1-, 3-, and 5-year OS for HCC patients. Discrimination and calibration were used to validate the nomogram. The discrimination of the nomogram was estimated utilizing the concordance index (C-index) via a bootstrap method with 340 resamples. The calibration plot was applied for assessing the consistence between nomogram prediction and practical observation.
Receiver Operating Characteristic Curve and Decision Curve Analysis
The time-dependent ROC curve and its corresponding area under the curve (AUC) value were established through “survivalROC” R package, thus assessing the discriminative accuracy of the predictive nomogram. The AUC value ranges from 0 to 1. The model presents a perfect discrimination when the AUC value is equal to 1. Conversely, the AUC value of 0.5 indicates a random capability to discriminate outcome (37).
DCA is a novel approach to assess the potential clinical net benefit (NB) of prognostic prediction models and to formulate better clinical strategies (38). NB is defined as a pivotal value that sums the benefits (true positives) and subtracts the harms (false positives) (37, 39). It can be plotted for a range of reasonable exchange rates in a decision curve where the potential utility of each decision strategy at each threshold probability is visualized (40). In the present study, we developed DCA by “rmda” R package, thus comparing the clinical utility of nomogram model, TNM stage, and GINS4 expression level (41).
Analysis of GINS4 Methylation in HCC
To explore the mechanism of the dysregulation of GINS4 in HCC, the DiseaseMeth database (http://biobigdata.hrbmu.edu.cn/diseasemeth/) the UALCAN (http://ualcan.path.uab.edu/index.html), and the Wanderer (http://maplab.imppc.org/) database were used to screen for potential methylation sites in the whole sequence of GINS4 DNA and to investigate the correlations between clinicopathologic parameters of HCC patients, GINS4 expression and its methylation values.
Screening of Genes Co-expressed With GINS4
A cluster of lncRNAs, miRNAs, mRNA, and RNA-binding protein (RBP) as well as transcription factors (TF) which interact with GINS4 were identified via the RAID database (RAID v2.0, www.rna-society.org/raid/) (42). The online cBioPortal database (http://www.cbioportal.org/) and GEPIA database (http://gepia.cancer-pku.cn/index.html) were queried to acquire co-expressed genes of GINS4 in HCC. Spearman correlation coefficient and Pearson correlation coefficient (PCC) were used to determine the degree of correlation between GINS4 and its co-expressed genes. R software (version 3.6.3) was used to analyze the HCC transcriptome expression matrix screened from the TCGA database to identify genes co-expressed with GINS4.
Functional Enrichment Analysis
We performed GO and KEGG of genes co-expressed with GINS4 using “clusterProfiler” R package to predict the biological process of GINS4 in HCC (43, 44). HCC with a functional gene set were further determined via GSEA software downloaded from https://www.broadinstitute.org/gsea/, thus acquiring significant biological processes enriched by GINS4 (45). The pathway with a nominal P < 0.05 and FDR < 0.05 was significant (42).
Statistical Analysis
Statistical analysis and graphic production were conducted through R language software (R 3.6.3 version). Chi-square test and Fishers exact test were applied for analyzing the correlation between GINS4 expression and clinicopathological parameters of HCC patients. Wilcoxon rank sum test were applied for comparing the GINS4 expression level in different groups. The ROC curves and corresponding AUC values were applied for determining the diagnostic significance of GINS4 and AFP levels in HCC samples in contrast to the control tissues. It was significant when P value was below 0.05.
Results
GINS4 Is Overexpressed in HCC and Significantly Associated With Clinicopathological Characteristics of HCC Patients
We extracted the mRNA expression profiles from 50 adjacent normal liver samples and 374 HCC tissues in the TCGA database, thus unraveling 4344 DEGs (│log2FC│ > 1, P < 0.05, FDR < 0.05), including 2,019 upregulated DEGs (log2FC > 1, P < 0.05) and 2,325 downregulated DEGs (log2FC < −1, P < 0.05) (Figures 2A, B). Above 4,344 DEGs in the TCGA database were applied for conducting gene co-expression network through WGCNA approach, thus estimating pivotal and candidate mRNAs that modulated histopathological grade in the progression of HCC (46). Based on the standard scale‐free network distribution, the soft threshold power value was set as 5 (Supplementary Figures 1A, B). On the basis of the criterion of dynamic cut tree, the least gene number of every network and the cut-height for the integration of modules was 30 and 0.25, respectively. The correlation of characteristic genes in integrated modules was above 0.95. As revealed in Figure 2C, eight co-expression modules were identified among all genes via the Topological Overlap Matrix (TOM). The gray module indicated a gene set without significant association with any clinical characteristics. The heatmap was applied for estimating a cluster of correlated eigengenes (Figures 2D, E). We further evaluated the associations between MEs and clinical parameters, including gender, AFP level, Child-Pugh classification, T classification, N classification, M classification, TNM stage, histological grade, relapse, and survival status. There was the greatest correlation between the turquoise module and histological grade (r = 0.36, P < 0.0001) (Figure 2F), which was chosen as a module of interest to be further analyzed. GINS4 (GS = 0.3262, MM = 0.8681, P < 0.0001) was one of hub genes significantly associated with histological grade in the turquoise module (Figure 2G), indicating that GINS4 potentially predicts the prognosis of HCC base on histological grade.
Figure 2 Bioinformatic analysis of HCC samples and matched normal adjacent tissues extracted from the TCGA dataset. (A, B) The heatmap and volcano plot of 4,344 DEGs between 374 HCC and 50 normal tissues. (C) Gene dendrogram of DEGs clustering and module screening according to mRNA expression profiles. (D, E) Heatmap of the correlation coefficient in the modules. (F) Correlations of module eigengenes and a series of clinical traits. (G) Analogous scatter plots for the turquoise module. DEGs, differently expressed genes.
Data from the Tumor Immune Estimation Resource (TIMER) database demonstrated that GINS4 mRNA expression was significantly elevated in multiple solid tumors, such as HCC, lung squamous cell carcinoma, gastric cancer, cholangiocarcinoma, and esophageal cancer (Figure 3A). To further identify the GINS4 level in HCC, we extracted the mRNA expression profiles from TCGA and GEO databases. The results revealed that compared with the normal liver tissue samples, GINS4 was significantly overexpressed in the HCC samples from the TCGA database (Figure 3B) and four GEO datasets (including GSE14520, GSE25097, GSE54236, and GSE76427) (Figures 3C–F). Additionally, we undertook IHC of HCC samples and matched liver tissues from 35 cases with primary HCC to further investigate GINS4 protein expression level in HCC. GINS4 protein expression was significantly increased in HCC tissues compared with adjacent normal liver samples (P < 0.01) (Figure 4A). Furthermore, greater GINS4 level was positively related to advanced TNM stage (TNM stage I and II versus TNM stage III and IV, P < 0.05) (Figure 4B) and worse pathological differentiation (Grade 1 and 2 versus Grade 3 and 4, P < 0.01) (Figure 4C). These results highlight that GINS4 is prominently overexpressed at both mRNA and protein levels in HCC and increased GINS4 expression potentially indicates the progression of HCC.
Figure 3 GINS4 mRNA levels are significantly upregulated in HCC. GINS4 mRNA expression in HCC and adjacent normal samples extracted from (A) the TIMER, (B) the TCGA, (C) GSE14520, (D) GSE25097, (E) GSE54236, and (F) GSE76427 database. (*P < 0.05, ***P < 0.001, and ****P < 0.0001).
Figure 4 GINS4 protein is greatly expressed in HCC samples and is positively related to disease stage. IHC analysis revealed GINS4 protein level (A) between HCC tissues and adjacent normal liver tissues, (B) between TNM I/II stage and TNM III/IV stage HCC, and (C) between histological grade 1/2 and histological grade 3/4 HCC. (magnification, 200× and 400×; *P < 0.05 and **P < 0.01). IHC, immunohistochemistry.
Table 1 demonstrated the correlation between GINS4 mRNA levels and clinicopathologic parameters of 371 HCC patients extracted from the TCGA. Statistical analyses suggested that GINS4 expression was significantly correlated with age (P = 0.009), gender (P = 0.001), AFP level (P = 0.049), T classification (P = 0.007), TNM stage (P = 0.011), histologic grade (P < 0.001), the status of residual tumor (P = 0.023), and relapse (P = 0.018) as well as vital status (P = 0.004). Specifically, GINS4 expression was significantly greater in HCC patients that belong to age ≤60 years old (P = 0.011), female (P = 0.010), AFP level >400 (P = 0.007), with residual tumor (P = 0.009), relapse (P = 0.017) (Figures 5A–E). Similarly, GINS4 mRNA level was significantly greater in histologic grade 3 HCC than histologic grade 1 HCC (P < 0.001) (Figure 5F). We also found progressive increase in the GINS4 expression with advanced tumor T classification and TNM stage (P < 0.0001) (Figures 5G, H), highlighting that GINS4 expression was positively related to the progression of HCC. In contrast, there were no correlations between GINS4 mRNA level and other clinicopathological characteristics, including the status of liver fibrosis (P = 0.5), HBV infection (P = 0.73), HCV infection (P = 0.63), vascular invasion (P = 0.26), and tumor status (P = 0.18) (Supplementary Figures 2A–E).
Table 1 Correlation between GINS4 expression and clinicopathological variables in 371 HCC patients extracted from TCGA database.
Figure 5 GINS4 mRNA level is significantly correlated with multiple clinicopathologic parameters in 371 HCC cases from the TCGA database. GINS4 mRNA level in HCC patients stratified by (A) age, (B) gender, (C) AFP level, (D) residual tumor, (E) relapse, (F) histologic grade, (G) T classification, and (H) TNM stage. (*P < 0.05, **P < 0.01, ***P < 0.001, and ****P < 0.0001).
GINS4 Can Effectively Distinguish HCC Patients From Nontumor Individuals
To identify the diagnostic significance of GINS4 in HCC, ROC curve was applied for analyzing the AUC of GINS4 expression stratified by clinical variables of HCC patients in the TCGA database. As revealed in Supplementary Figure 3A, GINS4 could effectively distinguish normal liver samples from HCC samples with an AUC of 0.865 (95% CI = 0.828–0.903). Additionally, the AUC value for the capacity of GINS4 expression level to differentiate HCC samples at TNM I, II, and III stage from adjacent tumor tissues was 0.835 (95% CI = 0.796–0.896), 0.878 (95% CI = 0.822–0.937), and 0.906 (95% CI = 0.840–0.946), respectively (Supplementary Figures 3B–D). Thus, our findings revealed that GINS4 displays a favorable ability to distinguish HCC patients and healthy individuals, even for early-stages HCC.
Furthermore, we formulated ROC curves to differentiate HCC patients from liver cirrhosis cases extracted from the GSE25097 and GSE63898 databases. In the GSE25097 dataset, GINS4 and alpha-fetoprotein (AFP) mRNA expression were both prominently greater in HCC than liver cirrhosis tissues (Figures 6A, B). The AUC of GINS4 (0.832, 95% CI = 0.781–0.882) was higher than that of AFP (0.787, 95% CI = 0.729–0.845) (P = 0.043) (Figure 6C). As for data from the GSE63898 dataset, GINS4 and AFP expression in HCC were remarkably elevated compared with liver cirrhosis (Figures 6D, E). The AUC of 0.708 (95% CI = 0.658–0.758) for GINS4 was significantly higher than that of 0.566 (95% CI = 0.510–0.622) for AFP (P < 0.0001) (Figure 6F). Our findings indicate that GINS4 is endowed with a relatively accurate performance to differentiate HCC from liver cirrhosis.
Figure 6 GINS4 harbors a great performance to differentiate HCC from liver cirrhosis. GINS4 mRNA level between liver cirrhosis and HCC samples from (A) GSE25097 and (D) GSE63898. AFP mRNA level between liver cirrhosis and HCC samples from (B) GSE25097 and (E) GSE63898. ROC curve demonstrates the discriminative efficiency of GINS4 and AFP between liver cirrhosis and HCC samples from (C) GSE25097 and (F) GSE63898. (****P < 0.0001). ROC, the receiver operating characteristic.
We further assessed the discriminative performance of GINS4 between low AFP-expressing HCC individuals and liver cirrhosis patients from the GSE25097 and GSE63898 datasets. GINS4 expression was significantly greater in HCC cases with low AFP expression than liver cirrhosis patients from above two datasets (Figures 7A, D). Conversely, there were no statistically significant difference in AFP level between liver cirrhosis and such HCC in two datasets (Figures 7B, E). ROC curve demonstrated that the AUC value for GINS4 was significantly higher compared with that for AFP in both datasets (0.754, 95% CI = 0.683–0.826 versus 0.575, 95% CI = 0.472–0.678, P = 0.0026 in GSE25097, Figure 7C; 0.654, 95% CI = 0.586–0.723 versus 0.479, 95% CI = 0.413–0.544, P = 0.0002 in GSE63898, Figure 7F). Above results highlight that GINS4 potentially serves as an instrument to screen low AFP-expressing HCC individuals.
Figure 7 GINS4 harbors a satisfactory performance to distinguish low AFP-expressing HCC from liver cirrhosis. GINS4 expression level between liver cirrhosis and low AFP-expressing HCC from (A) GSE25097 and (D) GSE63898. AFP expression level between liver cirrhosis and low AFP-expressing HCC from (B) GSE25097 and (E) GSE63898. ROC curve reveals the capability of GINS4 and AFP in differentiating liver cirrhosis and low AFP-expressing HCC from (C) GSE25097 and (F) GSE63898. (***P < 0.001, ****P < 0.0001). ns, no significance; ROC, the receiver operating characteristic.
Increased GINS4 Expression Predicts Unfavorable Prognosis in HCC Patients
We further estimated the prognostic significance of GINS4 in HCC through the Kaplan‐Meier curve. As revealed in Figure 8A, among all HCC patients, high CINS4- expressing HCC patients exhibited more unfavorable clinical outcome than those with low CINS4 expression (HR = 1.84, 95% CI = 1.21–2.8, P = 0.0038). We further investigated the correlation between GINS4 mRNA level and OS of HCC patients stratified by a variety of clinicopathologic features. Specifically, for HCC patients at histological grade 1, high GINS4-expressing patients were characteristic with worse OS compared with those with low GINS4 level (HR = 2.95, 95% CI = 1.04–8.40, P = 0.033) (Figure 8B). Additionally, high GINS4-expressing HCC patients belonging to age >60 years old (HR = 1.62, 95% CI = 1.03–2.56, P = 0.037), HBV infection-negative (HR = 1.53, 95% CI = 1.03–2.36, P = 0.032), and occurring relapse (HR = 1.63, 95% CI = 1.05–2.51, P = 0.029) subgroups were with significantly diminished OS than low GINS4-expressing HCC (Figures 8C–E).
Figure 8 High GINS4 expression is a significantly adverse prognostic factor in HCC patients. Survival analysis shows that elevated GINS4 expressions predicts poor prognosis in (A) all HCC patients from the TCGA database, (B) HCC patients at histological grade 1, (C) HCC patients with age >60 years old, (D) HBV infection-negative HCC patients, and (E) HCC subgroup with relapse. OS, overall survival.
Furthermore, we formulated univariate and multivariate Cox analyses to estimate the prognostic significance of GINS4 expression in HCC. In the univariate analysis, viral hepatitis infection, vascular invasion, T classification, M classification, TNM stage, tumor status, residual tumor, and GINS4 expression were significantly associated with the prognosis of HCC (P < 0.05) (Table 2). Above parameters were all incorporated into the multivariate Cox analyses. As shown in Table 2, multivariate Cox proportional hazards model revealed that high GINS4 expression (HR = 1.46, 95% CI = 1.01–2.1, P = 0.043) and advanced TNM stage (HR = 1.27, 95% CI = 1.01–1.62 for TNM stage II, P = 0.045; HR = 2.56, 95% CI = 1.66–3.96 for TNM stage III, P < 0.001) were independent unfavorable prognostic factors for the OS of HCC.
Table 2 Univariate and multivariate analysis for OS of 371 HCC patients from the TCGA database.
Development and Validation of Nomogram Model
To estimate the long-term survival of HCC individuals, we included all significant independent prognostic factors identified by the multivariate analyses, thus formulating a user-friendly nomogram with a C-index of 0.724. As revealed in Figure 9A, TNM stage and GINS4 mRNA level made great contributions to clinical outcome of HCC. The calibration plots for the OS probability of 1-year, 3-year, or 5-year in HCC patients showed an optimal consistency between nomogram prediction and practical observation (Figures 9B–D).
Figure 9 Nomogram to predict the long-term OS of 371 HCC patients from the TCGA database. (A) Nomogram model including TNM stage and GINS4 expression level to predict the 1-, 3-, and 5-year OS of HCC patients. The calibration plots of (B) 1-year, (C) 3-year, and (D) 5 year-OS of HCC patients. OS, overall survival.
The time-ROC curves were further employed for nomogram validation. The 1-year OS AUC of the nomogram model, TNM stage, and GINS4 expression level was 0.790 (95% CI = 0.710–0.855), 0.755 (95% CI = 0.673–0.826), and 0.720 (95% CI = 0.657–0.783), respectively (Figure 10A). Similarly, the nomogram model showed the highest 3-year OS AUC of 0.786 (95% CI 0.713–0.856), followed by TNM stage (AUC = 0.772, 95% CI = 0.715–0.842) and GINS4 expression level (AUC = 0.668, 95% CI = 0.601–0.719) (Figure 10B). The AUC at 5 years of the nomogram was 0.774 (95% CI = 0.661–0.877), significantly more discriminative than that of TNM stage being 0.751 (95% CI = 0.661–0.820) and that of GINS4 expression level being 0.676 (95% CI = 0.578–0.728) (Figure 10C). Thus, our nomogram comprising TNM stage and GINS4 expression level displayed a relatively satisfactory predictive accuracy for the long-term prognosis of HCC.
Figure 10 The nomogram model shows a relatively desirable predictive accuracy and clinical utility to predict OS of 371 HCC patients from the TCGA database. The time- dependent ROC demonstrated the AUC values of the nomogram, TNM stage, and GINS4 expression level at (A) 1-year, (B) 3-year, and (C) 5-year OS of HCC patients. (D) DCA for evaluating the clinical utility of nomogram model, TNM stage, and GINS4 expression level in predicting the OS of HCC cases. The x-axis revealed the range of threshold probabilities and the y-axis showed the NB. The gray line farthest left indicated all strategies. The horizontal black line represented none strategy. OS, overall survival; ROC, the receiver operating characteristic; DCA, Decision curve analysis; AUC, the area under the curve; NB, net benefit.
Furthermore, DCA was used to compare the clinical usefulness of nomogram with that of TNM stage and GINS4 expression level based on the threshold probability. Figure 10D revealed that the nomogram model exhibited a greater NB across a wider range of threshold probabilities for predicting long-term OS of HCC patients in the TCGA cohort, followed by TNM stage and GINS4 expression level. Specifically, the patients with OS probability between 0.23 and 0.58 would reap the highest NB if they selected the nomogram model. It also showed that TNM stage indicator would be applicable if OS probability of a patient was within the range of 0.25 to 0.59. Similarly, when OS probability of HCC patients was less than 0.28 or more than 0.45, decisions based on the GINS4 expression level would be meaningless. Therefore, above findings highlight that the nomogram is an excellent predicted evaluation model and it was superior to TNM stage or GINS4 expression level alone.
GINS4 Methylation Level Is Significantly Decreased in HCC Patients
Hypermethylation of CpG sites in promoters frequently results in transcriptional silencing. Conversely, hypomethylation of CpG sites in a gene body generally triggers an enhancive gene expression. A range of tumors are associated with promoter-specific hypomethylation and accompanied gene overexpression (47, 48). We selected the methylation site cg26367730 from the Wanderer database as the most statistically significant candidate site (Supplementary Table 1). As revealed in Figure 11A, GINS4 expression gradually decreased with incremental DNA methylation level in both adjacent normal liver samples and HCC tissues from the Wanderer database (P < 0.05), indicating that there is a potential negative association between the transcript expression of GINS4 and a number of CpG sites. Additionally, data from the UALCAN databases and the DiseaseMeth databases demonstrated that the total methylation value of GINS4 in the HCC samples was significantly decreased than normal liver samples (P < 0.001) (Figures 11B, C). Subsequently, we explored the correlation between GINS4 methylation level and clinicopathologic parameters of HCC patients from the UALCAN database. The lower methylation values of GINS4 in HCC patients were significantly associated with advanced TNM stage (P < 0.01), poorer pathological differentiation (P < 0.01), and lymph node metastasis (P < 0.05) (Figures 11D–F). Thus, above results highlighted that DNA hypomethylation, a primary epigenetic modification, potentially triggers GINS4 overexpression at the transcriptional level, thus exerts crucial effects on carcinogenesis and progression of HCC.
Figure 11 The correlation between GINS4 methylation and clinical characteristics of HCC patients. (A) The association between GINS4 expression and its methylation level of cg26367730 methylation site in HCC samples and matched normal liver tissues from the Wanderer. Total GISN4 methylation in HCC samples was significantly lower compared with the normal liver samples evaluated by (B) UALCAN and (C) DiseaseMeth databases. The association between GINS4 methylation levels and (D) TNM stage, (E) histological grade, and (F) lymph node status of HCC patients downloaded from UALCAN database.
High GINS4 Expression Positively Modulates Cell Cycle and PI3K/AKT/mTOR Signaling Pathway in HCC
As revealed in Figure 12A, GINS4 could interacted with six lncRNAs (TUG1, MALAT1, MANCR, SNHG16, STXBP5-AS1, lincMTX2, and SBF2-AS1) and multiple RBPs (including UPF1, YTHDF2, DGCR8, QK1, FMR1, and SRSF1). We also classified all HCC patients in the TCGA database into high and low GINS4 expression groups in accordance with the median GINS4 level. There was a total of 42 DEGs (│log2FC│>2) between above two groups, including 22 upregulated DEGs (log2FC > 2) and 20 downregulated DEGs (log2FC < −2), which co-expressed with GINS4 in HCC samples (Figure 12B). The top 200 co-expressed genes of GINS4 were obtained from the cBioPortal dataset (Spearman correlation coefficient ≥0.618, P value ≤6.26e-40) (Figure 12C). Meanwhile, the GEPIA database was applied to screen the top 200 genes co-expressed with GINS4 (PCC ≥ 0.62) (Supplementary Table 2). We cross-referenced the co-expressed genes from the above three databases to obtain a total of 41 common GINS4 co-expressed genes (Figure 12D).
Figure 12 High GINS4 expression positively modulates the cell cycle pathway in HCC development. (A) RAID database revealed certain miRNAs, lncRNAs, RBPs, TFs interacting with GINS4. (B) Heatmap of the DEGs co-expressed with GINS4 in HCC samples (│log2FC│> 2). (C) The top 25 genes associated with GINS4 transcript level in HCC based on the cBioPortal database. (D) Venn diagram illustrated the intersection of genes co-expressed with GINS4 in HCC among the TCGA, the cBioPortal, and the GEPIA database. The top 20 signaling pathways in HCC via (E) GO and (F) KEGG analysis. (G) GSEA of the whole DEGs co-expressed with GINS4 in HCC. DEGs, differently expressed genes; TF, transcription factors; NES, normalized enrichment score; RBP, RNA-binding protein.
We further conducted functional analysis of co‐expressed genes to investigate the biological classification of GINS4 in HCC. The top 20 biological process (BP) concerning the significantly enriched GO terms showed that these DEGs were primarily involved in the processes of nuclear division, positive regulation of cell cycle and DNA replication as well as cell cycle G1/S phase transition (Figure 12E), suggesting that GINS4 potentially facilitates HCC growth and proliferation through accelerating G1/S phase transition. Additionally, KEGG pathway analysis demonstrated that enrichment results were significantly correlated with the cell cycle and phosphoinositide 3 kinase (PI3K)- protein kinase B (AKT) signaling pathway (Figure 12F). The GSEA results revealed that the overexpression of GINS4 positively associated with the cell cycle [normalized enrichment score (NES) = 2.112, P < 0.0001], DNA replication (NES = 1.897, P = 0.002), base excision repair (NES = 1.953, P = 0.002), and tight junction (NES = 1.726, P = 0.048) as well as mechanistic target of rapamycin (mTOR) signaling pathway (NES = 1.679, P = 0.009) (Figure 12G).
Pearson correlation analysis further demonstrated that upregulation of GINS4 in HCC was significant positively correlated with the expression of phosphoinositide-3-kinase, catalytic, betapolypeptide (PIK3CB, known as the coding gene of PI3K P110 subunit, R2 = 0.35, P < 0.0001), AKT1 (known as the coding gene of AKT subunit, R2 = 0.23, P < 0.0001) and MTOR (R2 = 0.27, P < 0.0001) as well as CCND1 (known as the coding gene of Cyclin D1 that can promote the G1/S phase transition of mitosis, R2 = 0.16, P < 0.0001) (Supplementary Figures 4A–D). Thus, these studies indicated that GINS4 potentially participates in the regulation of PI3K/AKT/mTOR and cyclin D1 level, thus facilitating the occurrence and progression of HCC.
Discussion
HCC is a highly malignant tumor characterized with unfavorable clinical outcome and extremely high rates of mortality. Therefore, further investigation of HCC oncogenes is conducive to disclose novel and promising prognostic biomarkers and druggable targets, thus improving the clinical outcome of HCC. GINS4, a component of GINS complex, has been demonstrated a series of crucial functions in the biological process, including positive modulating in the initiation and prolongation of DNA replication, accelerating the transition of the cell cycle G1/S phase in eukaryotic cells, conferring protection against DNA damage in in both normal cells and cancer cells (20, 49–52). Significantly increased GINS4 level has been revealed in a series of human cancers, such as CRC (19, 20), NSCLC (22), gastric cancer (23), bladder cancer (21), and pancreatic cancer (24), highlighting the pivotal role of GINS4 in tumorigenesis. Nevertheless, the effect of GINS4 on HCC is relatively indistinct. Herein, our study was designed to identify the expression and the clinical and biology significance of GINS4 in HCC.
In our report, we conducted WGCNA co-expression network and revealed that GINS4 was one of hub DEGs most relevant to histological grade of HCC. GINS4 was overexpressed in HCC samples, and the expression level of GINS4 was significantly positively correlated with TNM stage and histological grade, indicating that GINS4 is an oncogene of HCC. ROC curves also demonstrated that GINS4 expression level could effectively distinguish HCC patients from non-tumor individuals (such as healthy controls and patients with liver cirrhosis). Additionally, the upregulation of GINS4 was associated with poor prognosis of HCC, especially in age >60 years old, histological grade G1, HBV-negative infection, and with recurrence subgroups, suggesting that GINS4 was a potentially independent risk factor affecting OS in HCC patients. The diagnostic and prognostic significance of GINS4 in other human tumors has also been confirmed. For example, the IHC results on tissue microarrays of 106 CRC patients revealed that enhanced GINS4 expression was positively related to advanced T stage, advanced TNM stage, and poor pathological differentiation (20). Additionally, multivariate analysis showed that GINS4 expression level in lung cancer was independent of clinical risk factors, such as gender, smoking, tumor differentiation, and tumor size, whereas it was associated with TNM stage and lymph node metastasis. The Kaplan-Meier curve also presented that high GINS4 expression predicted undesirable prognosis of all lung cancer patients and lung adenocarcinoma cases. Notably, there was no statistically significant correlation between GINS4 level and the survival of patients with lung squamous cell carcinoma (22). Similarly, gastric cancer patients with strongly positive GINS4 staining were characterized with shorter OS and DFS, suggesting that GINS4 may be a promising molecular target in the diagnosis and therapy of gastric cancer (23).
Furthermore, we found that GINS4 potentially positively modulated the cell cycle in HCC through accelerating the transition of mitotic G1/S phase and participated in malignant progression via PI3K/AKT/mTOR pathway based on GO and KEGG analysis. Pearson correlation analysis also demonstrated the significantly positive correlation between GINS4 mRNA and PI3KCB, AKT1, MTOR, and CCND1 transcriptome levels. PI3K/AKT/mTOR pathway is frequently activated in various human cancers, contributing to diversiform oncogenic transformation processes (such as stimulation of proliferation, survival, metabolic reprogramming, metastasis, and inhibition of apoptosis, autophagy, and aging) (53–55). Specifically, GINS4 could directly activate PI3K/AKT and MAPK/ERK pathways, thus accelerating cell proliferation and apoptosis in gastric cancer and CRC (20, 23). AKT is upregulated in 71% of HCC samples, thus accelerating the progressive growth of HCC. The activation of mTOR signaling is also revealed in above 48% of HCC samples and is related to undesirable prognosis in HCC therapy (55). As a pivotal cell cycle regulator, CyclinD1 is essential for accelerating the G1/S phase transition. CCND1, the coding gene of CyclinD1, has also been identified as a candidate proto-oncogene. The amplification and overexpression of CCND1 can alter the progression of the cell cycle and may be involved in the occurrence of tumors (56). Notably, Krüppel-like factor 4 (KLF4) diminishes GINS4 expression through binding to the promoter of GINS4, thus suppressing the development of CRC (20). Lymphoid-specific helicase (LSH) stabilizes and enhances GINS4 expression via binding to 3’UTR region of GINS4, thus facilitating lung cancer development (22). IL-6-induced the upregulation of DNA-methyltransferase (DNMT) inhibits miR-370, leading to high GINS4 expression and tumor growth in bladder cancer (21). Thus, suppression of GINS4 potentially represents a novel strategy to retard tumor development.
In conclusion, GINS4 is upregulated in HCC and high expression of GINS4 is significantly related to shorter survival in HCC patients. GINS4 may positively modulate the cell cycle process of HCC and potentially trigger the tumorigenesis and progression of HCC in a PI3K/AKT/mTOR dependent manner, which needs to be further experimental verification.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Ethics Statement
The studies involving human participants were reviewed and approved by the medical ethics committee of the Third Xiangya Hospital. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
PGC and ZZ designed/planned the study and wrote the paper. ZZ and PC performed experiment operation and computational modeling, and acquired and analyzed data. ZZ, PC, and HX performed imaging analysis. ZZ, PC, HX, and PGC participated in discussion of related data. ZZ drafted the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the National Natural Science Foundation of China (81872473) and the Hunan Province Science and Technology plan (2017SK2052).
Conflict of Interest
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.
Acknowledgments
The authors of this study have no contribution to TCGA and GEO data collection. We would like to thank the TCGA and GEO database for open access.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2021.654185/full#supplementary-material
Supplementary Figure 1 | Analysis of network topology for multiple soft-thresholding powers. (A) The scale-free fit index and (B) the mean connectivity of multiple soft-thresholding powers.
Supplementary Figure 2 | GINS4 mRNA expression in HCC patients stratified by multiple clinicopathologic parameters, including (A) liver fibrosis, (B) HBV infection, (C) HCV infection, (D) vascular invasion, (E) tumor status.
Supplementary Figure 3 | GINS4 exhibits a promising capability to discriminate HCC patients from healthy individuals from the TCGA database. ROC curve for (A) all HCC patients, (B) HCC patients at TNM I stage, (C) HCC patients at TNM II stage, and (D) HCC patients at TNM III stage. ROC, the receiver operating characteristic.
Supplementary Figure 4 | GINS4 transcriptome level is positively correlated with genes associated with PI3K/AKT/mTOR pathway in 371 HCC patients from the TCGA database. Pearson correlation analysis revealed the correlation between GINS4 mRNA level and (A) PIK3CB, (B) AKT1, (C) MTOR, and (D) CCND1.
References
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
2. Dong Z, Chen Y, Yang C, Zhang M, Chen A, Yang J, et al. STAT gene family mRNA expression and prognostic value in hepatocellular carcinoma. Onco Targets Ther (2019) 12:7175–91. doi: 10.2147/ott.s202122
3. Zhang Z, Chen P, Xie H, Cao P. Using circulating tumor DNA as a novel biomarker to screen and diagnose hepatocellular carcinoma: A systematic review and meta-analysis. Cancer Med (2020) 9(4):1349–64. doi: 10.1002/cam4.2799
4. Zhang Z, Tang H, Chen P, Xie H, Tao Y. Demystifying the manipulation of host immunity, metabolism, and extraintestinal tumors by the gut microbiome. Signal Transduct Targ Ther (2019) 4:41. doi: 10.1038/s41392-019-0074-5
5. Liu Y, Yang W, VoPham T, Ma Y, Simon TG, Gao X, et al. Plant-based and animal-based low-carbohydrate diets and risk of hepatocellular carcinoma among US men and women. Hepatology (2020) 73(1):175–85. doi: 10.1002/hep.31251
6. Mendizabal M, Pinero F, Ridruejo E, Wolff FH, Anders M, Reggiardo V, et al. Disease Progression in Patients With Hepatitis C Virus Infection Treated With Direct-Acting Antiviral Agents. Clin Gastroenterol Hepatol (2020) 18(11):2554–63. doi: 10.1016/j.cgh.2020.02.044
7. Juhling F, Hamdane N, Crouchet E, Li S, El Saghire H, Mukherji A, et al. Targeting clinical epigenetic reprogramming for chemoprevention of metabolic and viral hepatocellular carcinoma. Gut (2020) 70(1):157–69. doi: 10.1136/gutjnl-2019-318918
8. Sangiovanni A, Alimenti E, Gattai R, Filomia R, Parente E, Valenti L, et al. Undefined/non-malignant hepatic nodules are associated with early occurrence of HCC in DAA-treated patients with HCV-related cirrhosis. J Hepatol (2020) 73(3):593–602. doi: 10.1016/j.jhep.2020.03.030
9. Couri T, Pillai A. Goals and targets for personalized therapy for HCC. Hepatol Int (2019) 13(2):125–37. doi: 10.1007/s12072-018-9919-1
10. Zhang Z, Xie H, Chen P, Cao P. Development and Identification of a Nomogram Prognostic Model for Patients with Primary Clear Cell Carcinoma of the Liver. Med Sci Monit (2020) 26:e919789. doi: 10.12659/msm.919789
11. Takaya J, Kusunoki S, Ishimi Y. Protein interaction and cellular localization of human CDC45. J Biochem (2013) 153(4):381–8. doi: 10.1093/jb/mvt004
12. Kubota Y, Takase Y, Komori Y, Hashimoto Y, Arata T, Kamimura Y, et al. A novel ring-like complex of Xenopus proteins essential for the initiation of DNA replication. Genes Dev (2003) 17(9):1141–52. doi: 10.1101/gad.1070003
13. Takayama Y, Kamimura Y, Okawa M, Muramatsu S, Sugino A, Araki H. GINS, a novel multiprotein complex required for chromosomal DNA replication in budding yeast. Genes Dev (2003) 17(9):1153–65. doi: 10.1101/gad.1065903
14. Pacek M, Tutter AV, Kubota Y, Takisawa H, Walter JC. Localization of MCM2-7, Cdc45, and GINS to the site of DNA unwinding during eukaryotic DNA replication. Mol Cell (2006) 21(4):581–7. doi: 10.1016/j.molcel.2006.01.030
15. Seo YS, Kang YH. The Human Replicative Helicase, the CMG Complex, as a Target for Anti-cancer Therapy. Front Mol Biosci (2018) 526. doi: 10.3389/fmolb.2018.00026
16. Joshi K, Shah VJ, Maddika S. GINS complex protein Sld5 recruits SIK1 to activate MCM helicase during DNA replication. Cell Signal (2016) 28(12):1852–62. doi: 10.1016/j.cellsig.2016.08.018
17. Labib K, Gambus A. A key role for the GINS complex at DNA replication forks. Trends Cell Biol (2007) 17(6):271–8. doi: 10.1016/j.tcb.2007.04.002
18. Gouge CA, Christensen TW. Drosophila Sld5 is essential for normal cell cycle progression and maintenance of genomic integrity. Biochem Biophys Res Commun (2010) 400(1):145–50. doi: 10.1016/j.bbrc.2010.08.033
19. Wei HB, Wen JZ, Wei B, Han XY, Zhang S. Expression and clinical significance of GINS complex in colorectal cancer. Zhonghua Wei Chang Wai Ke Za Zhi (2011) 14(6):443–7.
20. Rong Z, Luo Z, Zhang J, Li T, Zhu Z, Yu Z, et al. GINS complex subunit 4, a prognostic biomarker and reversely mediated by Kruppel-like factor 4, promotes the growth of colorectal cancer. Cancer Sci (2020) 111(4):1203–17. doi: 10.1111/cas.14341
21. Yamane K, Naito H, Wakabayashi T, Yoshida H, Muramatsu F, Iba T, et al. Regulation of SLD5 gene expression by miR-370 during acute growth of cancer cells. Sci Rep (2016) 6:30941. doi: 10.1038/srep30941
22. Yang R, Liu N, Chen L, Jiang Y, Shi Y, Mao C, et al. LSH interacts with and stabilizes GINS4 transcript that promotes tumourigenesis in non-small cell lung cancer. J Exp Clin Cancer Res (2019) 38(1):280. doi: 10.1186/s13046-019-1276-y
23. Zhu Z, Yu Z, Rong Z, Luo Z, Zhang J, Qiu Z, et al. The novel GINS4 axis promotes gastric cancer growth and progression by activating Rac1 and CDC42. Theranostics (2019) 9(26):8294–311. doi: 10.7150/thno.36256
24. Bu F, Zhu X, Yi X, Luo C, Lin K, Zhu J, et al. Expression Profile of GINS Complex Predicts the Prognosis of Pancreatic Cancer Patients. Onco Targets Ther (2020) 13:11433–44. doi: 10.2147/ott.s275649
25. Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics (2010) 26(1):139–40. doi: 10.1093/bioinformatics/btp616
26. Roessler S, Jia HL, Budhu A, Forgues M, Ye QH, Lee JS, et al. A unique metastasis gene signature enables prediction of tumor relapse in early-stage hepatocellular carcinoma patients. Cancer Res (2010) 70(24):10202–12. doi: 10.1158/0008-5472.can-10-2607
27. Tung EK, Mak CK, Fatima S, Lo RC, Zhao H, Zhang C, et al. Clinicopathological and prognostic significance of serum and tissue Dickkopf-1 levels in human hepatocellular carcinoma. Liver Int (2011) 31(10):1494–504. doi: 10.1111/j.1478-3231.2011.02597.x
28. Villa E, Critelli R, Lei B, Marzocchi G, Camma C, Giannelli G, et al. Neoangiogenesis-related genes are hallmarks of fast-growing hepatocellular carcinomas and worst survival. Results from a prospective study. Gut (2016) 65(5):861–9. doi: 10.1136/gutjnl-2014-308483
29. Grinchuk OV, Yenamandra SP, Iyer R, Singh M, Lee HK, Lim KH, et al. Tumor-adjacent tissue co-expression profile analysis reveals pro-oncogenic ribosomal gene signature for prognosis of resectable hepatocellular carcinoma. Mol Oncol (2018) 12(1):89–113. doi: 10.1002/1878-0261.12153
30. Zhou Q, Su X, Jing G, Ning K. Meta-QC-Chain: comprehensive and fast quality control method for metagenomic data. Genomics Proteomics Bioinf (2014) 12(1):52–6. doi: 10.1016/j.gpb.2014.01.002
31. Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinf (2008) 9:559. doi: 10.1186/1471-2105-9-559
32. Chen L, Yuan L, Qian K, Qian G, Zhu Y, Wu CL, et al. Identification of Biomarkers Associated With Pathological Stage and Prognosis of Clear Cell Renal Cell Carcinoma by Co-expression Network Analysis. Front Physiol (2018) 9:399:399. doi: 10.3389/fphys.2018.00399
33. Fujiya K, Terashima M, Ohshima K, Aizawa D, Sugino T, Serizawa M, et al. MAGEA10 expression is a predictive marker of early hepatic recurrence after curative gastrectomy for gastric and gastroesophageal junction cancer. Gastric Cancer (2020) 9:399. doi: 10.1007/s10120-020-01123-y
34. Nagy A, Lanczky A, Menyhart O, Gyorffy B. Validation of miRNA prognostic power in hepatocellular carcinoma using expression data of independent datasets. Sci Rep (2018) 8(1):9227. doi: 10.1038/s41598-018-27521-y
35. Balachandran VP, Gonen M, Smith JJ, DeMatteo RP. Nomograms in oncology: more than meets the eye. Lancet Oncol (2015) 16(4):e173–80. doi: 10.1016/s1470-2045(14)71116-7
36. Iasonos A, Schrag D, Raj GV, Panageas KS. How to build and interpret a nomogram for cancer prognosis. J Clin Oncol (2008) 26(8):1364–70. doi: 10.1200/jco.2007.12.9791
37. Mo S, Dai W, Xiang W, Li Q, Wang R, Cai G. Predictive factors of synchronous colorectal peritoneal metastases: Development of a nomogram and study of its utilities using decision curve analysis. Int J Surg (2018) 54(Pt A):149–55. doi: 10.1016/j.ijsu.2018.04.051
38. Chen P, Zhang Z, Chen X. Overexpression of PKMYT1 Facilitates Tumor Development and Is Correlated with Poor Prognosis in Clear Cell Renal Cell Carcinoma. Med Sci Monit (2020) 26:e926755. doi: 10.12659/msm.926755
39. Barbieri CE, Cha EK, Chromecki TF, Dunning A, Lotan Y, Svatek RS, et al. Decision curve analysis assessing the clinical benefit of NMP22 in the detection of bladder cancer: secondary analysis of a prospective trial. BJU Int (2012) 109(5):685–90. doi: 10.1111/j.1464-410X.2011.010419.x
40. Fei Z, Qiu X, Li M, Chen C, Li Y, Huang Y. Prognosis viewing for nasopharyngeal carcinoma treated with intensity-modulated radiation therapy: application of nomogram and decision curve analysis. Jpn J Clin Oncol (2020) 50(2):159–68. doi: 10.1093/jjco/hyz165
41. Van Calster B, Wynants L, Verbeek JFM, Verbakel JY, Christodoulou E, Vickers AJ, et al. Reporting and Interpreting Decision Curve Analysis: A Guide for Investigators. Eur Urol (2018) 74(6):796–804. doi: 10.1016/j.eururo.2018.08.038
42. Yi Y, Zhao Y, Li C, Zhang L, Huang H, Li Y, et al. RAID v2.0: an updated resource of RNA-associated interactions across organisms. Nucleic Acids Res (2017) 45(D1):D115–d8. doi: 10.1093/nar/gkw1052
43. Carbon S, Dietze H, Lewis SE, Munoz-Torres MC, Basu S, et al. Expansion of the Gene Ontology knowledgebase and resources. Nucleic Acids Res (2017) 45(D1):D331–d8. doi: 10.1093/nar/gkw1108
44. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res (2000) 28(1):27–30. doi: 10.1093/nar/28.1.27
45. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA (2005) 102(43):15545–50. doi: 10.1073/pnas.0506580102
46. Giulietti M, Occhipinti G, Principato G, Piva F. Identification of candidate miRNA biomarkers for pancreatic ductal adenocarcinoma by weighted gene co-expression network analysis. Cell Oncol (Dordr) (2017) 40(2):181–92. doi: 10.1007/s13402-017-0315-y
47. Morgan AE, Davies TJ, Mc Auley MT. The role of DNA methylation in ageing and cancer. Proc Nutr Soc (2018) 77(4):412–22. doi: 10.1017/s0029665118000150
48. Klutstein M, Nejman D, Greenfield R, Cedar H. DNA Methylation in Cancer and Aging. Cancer Res (2016) 76(12):3446–50. doi: 10.1158/0008-5472.can-15-3278
49. Khemlina G, Ikeda S, Kurzrock R. The biology of Hepatocellular carcinoma: implications for genomic and immune therapies. Mol Cancer (2017) 16(1):149. doi: 10.1186/s12943-017-0712-x
50. Simon AC, Zhou JC, Perera RL, van Deursen F, Evrin C, Ivanova ME, et al. A Ctf4 trimer couples the CMG helicase to DNA polymerase alpha in the eukaryotic replisome. Nature (2014) 510(7504):293–7. doi: 10.1038/nature13234
51. Zhu L, Zhao W, Lu J, Li S, Zhou K, Jiang W, et al. Influenza virus matrix protein M1 interacts with SLD5 to block host cell cycle. Cell Microbiol (2019) 21(8):e13038. doi: 10.1111/cmi.13038
52. Deng L, Wu RA, Sonneville R, Kochenova OV, Labib K, Pellman D, et al. Mitotic CDK Promotes Replisome Disassembly, Fork Breakage, and Complex DNA Rearrangements. Mol Cell (2019) 73(5):915–29.e6. doi: 10.1016/j.molcel.2018.12.021
53. Belinky F, Nativ N, Stelzer G, Zimmerman S, Iny Stein T, Safran M, et al. PathCards: multi-source consolidation of human biological pathways. Database (Oxford) (2015) 2015:1–13. doi: 10.1093/database/bav006
54. Samarin J, Laketa V, Malz M, Roessler S, Stein I, Horwitz E, et al. PI3K/AKT/mTOR-dependent stabilization of oncogenic far-upstream element binding proteins in hepatocellular carcinoma cells. Hepatology (2016) 63(3):813–26. doi: 10.1002/hep.28357
55. Rahmani F, Ziaeemehr A, Shahidsales S, Gharib M, Khazaei M, Ferns GA, et al. Role of regulatory miRNAs of the PI3K/AKT/mTOR signaling in the pathogenesis of hepatocellular carcinoma. J Cell Physiol (2020) 235(5):4146–52. doi: 10.1002/jcp.29333
Keywords: hepatocellular carcinoma, GINS4, expression, prognosis, biological significance
Citation: Zhang Z, Chen P, Xie H and Cao P (2021) Overexpression of GINS4 Is Associated With Tumor Progression and Poor Survival in Hepatocellular Carcinoma. Front. Oncol. 11:654185. doi: 10.3389/fonc.2021.654185
Received: 18 January 2021; Accepted: 08 March 2021;
Published: 25 March 2021.
Edited by:
Shu-Hong Huang, Shandong First Medical University, ChinaReviewed by:
Qingbin Cui, University of Toledo Medical Center, United StatesAnantha Koteswararao Kanugula, University of Massachusetts Medical School, United States
Copyright © 2021 Zhang, Chen, Xie and Cao. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Peiguo Cao, xy3caopg@csu.edu.cn