14-CpG-Based Signature Improves the Prognosis Prediction of Hepatocellular Carcinoma Patients

Jiang, Hong-ye; Ning, Gang; Wang, Yen-sheng; Lv, Wei-biao

doi:https://doi.org/10.1155/2020/9762067

BioMed Research International

On this page

Abstract Introduction Materials and Methods Results Discussion Abbreviations Data Availability Ethical Approval Consent Conflicts of Interest Authors’ Contributions Supplementary Materials References Copyright Related Articles

Research Article | Open Access

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

14-CpG-Based Signature Improves the Prognosis Prediction of Hepatocellular Carcinoma Patients

Hong-ye Jiang,¹Gang Ning,²Yen-sheng Wang,³and Wei-biao Lv¹

Academic Editor: Dimitrios P. Bogdanos

Received26 Aug 2019

Accepted13 Dec 2019

Published08 Jan 2020

Abstract

Background. Epigenetic dysregulation via alteration of DNA methylation often occurs during the development and progression of cancer, including hepatocellular carcinoma (HCC). In the past, many patterns of single-gene DNA methylation have been extensively explored in the context of HCC prognosis prediction. However, the combined model of a mixture of CpGs has rarely been evaluated. In the present study, we aimed to develop and validate a CpG-based signature model for HCC patient prognosis. Methods. Data from methylation profiling of GSE73003, GSE37988, and GSE57958 from the Gene Expression Omnibus (GEO) database and 371 HCC patients from the Cancer Genome Atlas (TCGA) were downloaded. The 371 HCC patients were randomly divided into a development cohort (N = 263) and a validation cohort (N = 108). Two algorithms, least absolute shrinkage and selection operator (LASSO) and robust likelihood-based survival analysis, were used to select the most significant CpGs associated with overall survival (OS) time and were used to develop and validate a methylation-based signature (MSH) for HCC patient prognosis. In addition, the prognostic efficacy of the MSH was compared with that of AJCC TNM classification and other CpG-based MSHs from TCGA. Finally, a nomogram incorporating the MSH and clinicopathologic factors was also developed. Results. Fourteen differential CpGs associated with OS were identified in HCC patients. The MSH, based on these 14 differential CpGs, could effectively divide HCC patients into two distinct subgroups with high risk or low risk of death () in the development cohort (26.35 vs 83.18 months, HR = 3.83, 95% CI: 2.56–5.90, ) and in the validation cohort (40.37 vs 107.03 months, HR = 2.23, 95% CI: 1.22–4.17, ). Univariate analysis showed that the MSH was significantly associated with OS, and the multivariate analysis also showed that the MSH was an independent prognostic factor for the OS of HCC patients in the two cohorts. In addition, stratified survival analysis indicated that the MSH still exhibited good prognostic value in different subgroups classified by AFP, cirrhosis, Child-Pugh A, tumor histologic grade, and AJCC stage. Moreover, time-dependent ROC analysis showed better performance of the MSH in predicting 3-year and 5-year survival of HCC patients than of AJCC stage and other CpG-based signatures from TCGA. The MSH-based nomogram also performed well in predicting 1-year, 3-year, and 5-year OS (C-index: 0.709). Conclusion. The 14-CpG-based signature is significantly associated with OS and may be used as a novel prognostic biomarker for HCC patients.

1. Introduction

Hepatocellular carcinoma (HCC) is predicted to have become the sixth most common cancer and the fourth leading cause of cancer-related death worldwide in 2018. Each year, an estimated 841,000 patients develop HCC, and 782,000 patients die from this disease [1]. Nevertheless, the threat of HCC has not been mitigated, as evidenced by the rapidly increasing incidence of HCC and the high recurrence rate of 50% among early-stage HCC patients after surgery [2, 3]. Late diagnosis and limited treatment options were suggested to account for the high mortality rate in advanced HCC patients [4]. Apart from working to find new treatment methods for this deadly disease, scientists are exploring new models for the early diagnosis and prevention of HCC to improve the prognosis of HCC.

It is well known that cancer genetics, including mutations and single-nucleotide polymorphisms (SNPs), and aberration of epigenetic regulation play important roles in the development and progression of HCC [5–7]. As one of the major epigenetic regulations, DNA methylation is reported to take part in the formation of many malignant tumors, including HCC [8, 9]. Mechanistically, aberrant DNA hypermethylation on the promoter region of CpG islands would result in the silencing of tumor suppressor genes, thus leading to the overexpression of oncogenes [10]. DNA hypermethylation on promoter CpG islands has been observed to be associated with the clinicopathological characteristics and prognosis of HCC patients in previous studies [11–13]. Identifying specific abnormal methylated CpGs may be of promising value for the diagnosis, prognosis, and even treatment of HCC.

The prognostic value of many single-gene DNA methylation patterns for HCC has been extensively explored. However, a combined model that includes assorted CpGs has rarely been evaluated. In the present study, we identified 14 differentially methylated CpGs related to HCC prognosis. We utilized the methylation profiling data of HCC from the Cancer Genome Atlas (TCGA) and developed a methylation-based signature for HCC (MSH) in the development cohort. Next, we validated the model in the validation cohort. Last, we compared the prognostic efficacy of MSH with that of the AJCC classification and other CpG-based MSHs from TCGA [14].

2. Materials and Methods

2.1. Ethics Statement

All data in the study were obtained from online databases, including Gene Expression Omnibus (GEO) and TCGA. Informed consent was obtained from the patients before the study. The study was also approved by the Ethics Committee of the Shunde Hospital of Southern Medical University.

2.2. Methylation Data Collection and Processing from GEO and TCGA

In our present study, DNA methylation profiles between primary HCC tumors and their nontumor counterparts from GSE73003 (including 20 paired tumor and nontumor tissues from Japan), GSE37988 (62 paired tumor and nontumor tissues from Taiwan), and GSE57958 (99 paired tumor and nontumor tissues from Singapore) were first obtained from GEO (https://www.ncbi.nlm.nih.gov/geo/). All three of these datasets were assessed on GPL8490 (Illumina Human Methylation 27 BeadChip). Next, GEO2R, an online software package, was used to identify differential CpGs. The cut-off criterion of differential CpGs was . To find the most significant differentially methylated markers, the top 1,000 CpGs with the lowest values of each dataset were selected (supplementary materials 1, 2, and 3). Finally, an online tool (http://bioinformatics.psb.ugent.be/webtools/Venn/) was used to identify the overlapping CpGs among each of the 1,000 top CpGs of GSE73003, GSE37988, and GSE57958 (Supplementary material 4).

After identifying the most significant differential CpGs between HCC tumors and nontumor tissues, we next verified these CpGs among HCC patients from TCGA (https://cancergenome.nih.gov/). DNA methylation profiling data of 377 HCC patients were downloaded from TCGA. The methylation profiling data were assessed on the GPL13534 platform (Illumina Human Methylation 450 BeadChip), and the methylation level was presented as a β value, which was calculated as the ratio of the intensity of the methylated bead type to the combined locus intensity and ranged from 0 to 1. Subsequently, clinical characteristics including sex, age, BMI, APF, cirrhosis, Child-Pugh stage, adjacent hepatic tissue inflammation, tumor histological grade, surgical margin resection status, AJCC TNM stage, and overall survival (OS) time were also downloaded. Six of 377 HCC patients were excluded because of the absence of OS data. In total, 371 HCC patients with available methylation data and clinical parameters were included in the present study. The clinical parameters of the HCC patients are summarized in Table 1.

2.3. Identification and Selection of HCC Prognosis-Related CpGs

Three hundred seventy-one HCC patients were randomly divided into a development cohort (N = 263) and a training cohort (N = 108) with an allocation of 7 : 3 performed by R software. The development cohort was used to identify key HCC prognosis-related CpGs and develop the MSH. Two different algorithms, least absolute shrinkage and selection operator (LASSO) analysis [15] and robust likelihood-based survival analysis [16, 17], were used to select the most significant methylation markers. Overlapping CpGs between the two selection methods were finally identified as the HCC prognosis-related CpGs.

2.4. Development and Validation of the MSH

After the key prognosis-related CpGs were selected, we next used them to develop the MSH by multivariable Cox regression analysis. With this model, a risk score for each HCC patient was calculated. HCC patients were further classified into a high-risk group and a low-risk group based on the cut-off value of the median risk score. The OS difference between the high-risk patients and low-risk patients was analyzed by Kaplan-Meier analysis. Then, the MSH was validated in the validation cohort. Univariate and multivariate Cox regression analyses were used to further assess the association of MSH with OS in the development and validation cohorts. Furthermore, stratified analysis was also performed to explore the influence of other major clinicopathologic factors (including AFP, cirrhosis, Child-Pugh stage, tumor histologic grade, and AJCC TNM stage) on the prognostic value of MSH in the total cohort by Kaplan–Meier analysis.

2.5. Establishment of a Time-Dependent Receiver Operating Characteristic (ROC) Curve and an MSH-Based Nomogram

To further assess the predictive accuracy and sensitivity of the MSH, time-dependent ROC analysis was performed with HCC patients in the total cohort. The areas under the ROC curve (AUCs) of the MSH for predicting 1-year, 3-year, and 5-year OS were calculated and used for comparisons with other models. Moreover, to make MSH more clinically applicable, an MSH-based nomogram was also developed.

2.6. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) Analysis of the MSH

To explore the biological function and pathways of the MSH, GO and KEGG analyses were conducted. First, the 50 most frequently altered genes related to these 14 genes were downloaded from cBioPortal (http://www.cbioportal.org). The biological function of these 50 genes and the 14 genes were then analyzed by GO and KEGG in the Database for Annotation, Visualization, and Integrated Discovery (DAVID) (https://david.ncifcrf.gov/summary.jsp). The detailed method was described in our previous study [18].

2.7. Statistical Analysis

Statistical analysis was performed with R software (R version 3.5.1) and GraphPad Prism software (version 6). Univariate and multivariate Cox regression analyses were performed with the survival and survminer packages. The robust likelihood-based survival analysis was performed with the survivalROC and rbsurv packages, and the LASSO analysis was conducted with the glmnet and survival packages. Time-dependent ROC analysis was performed with the ROCR and rms packages. The nomogram was constructed with the rms and survival packages and was evaluated by the concordance index and calibration plots. Kaplan–Meier analysis was performed with GraphPad Prism software and was compared with the log-rank test. was considered statistically significant.

3. Results

3.1. Basic Characteristics of the 371 HCC Patients

The flowchart of the present study is shown in Figure 1, and the basic characteristics of the 371 HCC patients are summarized in Table 1. Of the 371 HCC patients, 132 patients (35.6%) died, and 178 patients (25.2%) developed recurrence. The median OS time was 19.78 months (ranging from 0.03 to 120.73 months).

3.2. Selection of Key HCC Prognosis-Related CpGs

Based on the primary filter criteria, 426 differential CpGs between primary HCC tumors and the corresponding nontumor tissue were identified (Figure 2(a)). Then, we validated these selected CpGs in the HCC patients in TCGA and identified 288 CpGs that were detected by two different DNA methylation detection methods (Supplementary material 5). Next, LASSO analysis was used to obtain a set of 21 CpGs (Figures 2(b)–2(d)). A robust likelihood-based survival analysis was also performed and identified a set of 33 CpGs (Table 2). There were 14 overlapping CpGs between the two selection methods (Figure 2(d)), which corresponded to cg00504595 (TNF receptor superfamily member 19, TNFRSF19), cg04711324 (ras-like without CAAX2, RIT2), cg06226384 (calcium voltage-gated channel auxiliary subunit gamma 5, CACNG5), cg07014174 (keratin-associated protein 11-1, KRTAP11-1), cg08668790 (zinc finger protein 154, ZNF154), cg15747595 (TSPY-like 5, TSPYL5), cg16673198 (copine 4, CPNE4), cg18343292 (membrane spanning 4-domains A7, MS4A7), cg18536148 (T-box 4, TBX4), cg21578906 (solute carrier family 5 member 4, SLC5A4), cg23163573 (sulfotransferase family 1C member 2, SULT1C2), cg24432073 (cyclin dependent kinase-like 2, CDKL2), cg24898863 (S100 calcium-binding protein A8, S100A8), and cg26059632 (small proline rich protein 2A, SPRR2A). With the two different algorithms, HCC prognosis-related markers were strictly selected.

(a)

(b)

(c)

(d)

3.3. Construction and Validation of the MSH

To comprehensively explore the association of these 14 selected CpGs with the prognosis of HCC patients, a MSH was built based on the coefficients weighted by multivariable Cox regression analysis in the development cohort (Table 3). The risk score was calculated as follows: risk score = (4.10 ∗ cg00504595) + (3.79 ∗ cg04711324) + (2.83 ∗ cg06226384) + (1.76 ∗ cg07014174) + (0.54 ∗ cg08668790) + (1.08 ∗ cg15747595) + (−2.29 ∗ cg16673198) + (−2.76 ∗ cg18343292) + (0.59 ∗ cg18536148) + (−1.93 ∗ cg21578906) + (−1.78 ∗ cg23163573) + (1.92 ∗ cg24432073) + (−2.28 ∗ cg24898863) + (−0.60 ∗ cg26059632). After the risk score for each patient in the development cohort was calculated, patients with a risk score >1.07 (median score) were assigned to the high-risk group (N = 132), and the other patients were assigned to the low-risk group (N = 131). The methylation levels of cg00504595, cg04711324, cg06226384, cg07014174, cg08668790, cg15747595, cg18536148, and cg24432073 in patients of the high-risk group tended to be higher than those in patients of the low-risk group, while the methylation levels of cg16673198, cg18343292, cg21578906, cg23163573, cg24898863, and cg26059632 tended to be lower in patients of the high-risk group (Figure 3(a)). Moreover, patients in the high-risk group had shorter OS time than those in the low-risk group (median survival time 26.35 vs 83.18 months, HR = 3.83, 95% CI: 2.56–5.90, , Figures 3(b) and 3(c)). To assess the utility and stability of the MSH, verification analysis was performed in the validation cohort (N = 108). Similarly, the MSH also successfully divided 54 patients (50%) into the high-risk group and the other 54 patients into the low-risk group. The OS time of the high-risk group was lower than that of the low-risk group (40.37 vs 107.03 months, HR = 2.23, 95% CI: 1.22–4.17, , Figure 4).

(a)

(b)

(c)

(a)

(b)

(c)

3.4. Prognostic Value of the MSH in HCC Patients

After indicating that the MSH could be used to categorize HCC patients into high-risk (poor OS) and low-risk groups (better OS), we further evaluated the prognostic value of the MSH among HCC patients. Univariate analysis showed that the MSH was significantly associated with OS in the development cohort (HR = 4.3, 95% CI: 2.691–6.871, , Table 4) and the validation cohort (HR = 1.979, 95% CI: 1.019–3.864, , Table 5). Moreover, multivariate analysis also showed that the MSH was an independent prognostic factor for OS in the two cohorts (development cohort: HR = 6.355, 95% CI: 2.524–16, , Table 4; validation cohort, HR = 3.379, 95% CI: 1.054–10.834, , Table 5).

3.5. Stratified Survival Analysis Based on Major Clinicopathological Factors in the Total Cohort

After the MSH was found to be an independent prognostic factor for the OS of HCC patients, we next performed stratified analysis to further explore the prognostic value of the MSH for patients classified by major clinicopathological factors in the total cohort. The number of patients divided into high-risk and low-risk groups and the log-rank tests are shown in Table 6. Our results indicated that the MSH still exhibited good prognostic value in different subgroups classified by AFP, cirrhosis, Child-Pugh A, tumor histologic grade, and AJCC stage (Figure 5), which, to some extent, suggested the greater reliability and general utility of the MSH.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

Figure 5

Stratified survival analysis based on major clinicopathological factors. The 14-CpG-based signature exhibited good prognostic value in different subgroups classified by AFP (a-b), cirrhosis (c-d), Child-Pugh A (e), tumor histologic grade (f-g), and AJCC stage (h-i). (a) Patients with AFP (<25) (N = 166). (b) Patients with AFP (>25) (N = 116). (c) Patients without cirrhosis (N = 136). (d) Patients with cirrhosis (N = 79). (e) Patients with Child-Pugh A (N = 222). (f) Patients with histologic grade G1/G2 (N = 233). (g) Patients with histologic grade G3/G4 (N = 132). (h) Patients with AJCC stage I/II (N = 259). (i) Patients with AJCC stage I/II (N = 88).

3.6. Predictive Value of the MSH for the OS of HCC Patients and Comparison with Other CpG-Based Models Based on TCGA

Time-dependent ROC cure analysis was used to assess the predictive value of the MSH in HCC patients in the total cohort, and this analysis was used to compare the MSH to other CpG-based models based on TCGA. As shown in Figure 6, the AUCs of the MSH for predicting 1-, 3-, and 5-year OS were 0.643, 0.712, and 0.757, respectively, while the AUCs of AJCC stage, which is often used as prognostic model for HCC patients, were 0.657, 0.668, and 0.636, respectively, suggesting that the MSH exhibited a better efficiency in predicting 3- and 5-year OS.

(a)

(b)

(c)

Recently, a five-CpG-based prognostic signature was constructed by Fang et al. on the basis of HCC patients from TCGA [14], and the AUCs of the MSH developed by Fang et al. for predicting 1-, 3-, and 5-year OS were 0.577, 0.587, and 0.603, respectively. Compared to the MSH developed by Fang et al., our 14-CpG-based prognostic signature showed a favorable predictive value in predicting 1-, 3-, and 5-year OS. However, further investigation into external HCC cohorts is needed.

3.7. Development of an MSH-Based Nomogram for OS Prediction in HCC Patients in the Total Cohort

To make MSH more clinically applicable, we developed an MSH-based nomogram to predict the 1-year, 3-year, and 5-year OS of HCC patients in the total cohort (Figure 7(a)). Clinicopathological factors, such as sex, age, AFP, cirrhosis, Child-Pugh stage, tumor histologic grade, AJCC stage, and surgical margin status, were included in the nomogram. The C-index for 1-year, 3-year, and 5-year OS prediction was 0.709, and the calibration plots also exhibited good consistency between the predicted OS and the actual OS (Figures 7(b)–7(d)), suggesting the good predictive value of our MSH-based nomogram.

(a)

(b)

(c)

(d)

3.8. Biological Function and Pathways of the MSH

To explore the biological function and pathways of the MSH, GO and KEGG analyses were performed. Our results showed that the MAPK signaling pathway and neurotrophin signaling pathway were affected by these 14 genes (Figure 8(a)), and as expected, these two common pathways were all reported to play important roles in the development and progression of HCC [19–21], which provided evidence for the rationality and molecular thesis of the MSH. In addition, biological processes such as GO:0006468 (protein phosphorylation), GO:0018105 (peptidyl-serine phosphorylation), and GO:0006351 (transcription, DNA-templated); molecular functions, such as GO:0004674 (protein serine/threonine kinase activity), GO:0004672 (protein kinase activity), and GO:0005524 (ATP binding); and cellular components such as GO:0005634 (nucleus), GO:0005622 (intracellular), and GO:0005856 (cytoskeleton) were also affected by these 14 genes (Figures 8(b)–8(d)).

(a)

(b)

(c)

(d)

4. Discussion

HCC is a highly malignant cancer with poor prognosis. It is still a great challenge to improve the clinical outcome of HCC patients because of the absence of effective prognostic biomarkers or models. In our present study, we aimed to develop and evaluate the prognostic value of methylation-based signatures for HCC patients. Fourteen candidate CpGs related to OS were identified in the development cohort by two distinct algorithms, including LASSO analysis and robust likelihood-based survival analysis. Unlike previous studies that used only one algorithm, we used two algorithms to help minimize the possibility of losing or missing important markers [22]. Subsequently, these 14 CpGs were used to develop an MSH in the development cohort and were validated in an internal validation cohort.

Our results showed that the MSH could effectively divide HCC patients into two distinct subgroups with high risk or low risk of death, suggesting the underlying clinical implications for the management of HCC patients. In addition, MSH was associated with OS and was also an independent prognostic factor for HCC patients. Moreover, stratified analysis also indicated good prognostic value in different subgroups classified by AFP, cirrhosis, Child-Pugh A, tumor histologic grade, and AJCC stage, which, to some extent, suggested the greater reliability and general utility of MSH. With the help of the MSH, high-risk HCC patients can be identified and can receive more intensive surveillance and even active adjuvant treatment to reduce recurrence and improve prognosis. Conversely, HCC patients with low risk may receive less active follow-up and even avoid adverse effects of adjuvant therapies. Therefore, MSH may be a useful method for establishing more individualized follow-up interval schedules and selecting therapeutic strategies for HCC patients after surgery.

The AJCC TNM stage is a well-known useful and common marker for predicting the prognosis of HCC. To further evaluate the predictive value of the MSH, we used time-dependent ROC analysis to compare the prediction efficacy between the MSH and AJCC stage. The prognostic predictive ability of the MSH was stable and good. The AUC for predicting 1-year, 3-year, and 5-year OS increased with increased prediction time (0.643, 0.712, and 0.757, respectively), suggesting the better accuracy of the MSH for long-time survival prediction, which is relatively important for patients at advanced stages. Compared to the AJCC stage (AUC = 0.657), the efficacy of the MSH in predicting 1-year OS was not inferior, but the efficacy of the MSH in predicting 3-year and 5-year OS was superior to that of the AJCC stage (AUC = 0.668, 0.636, respectively), indicating the advantage of MSH in predicting long-time survival. Recently, a five-CpG-based prognostic model was developed by Fang et al. on the basis of TCGA, and the AUCs of the MSH developed by Fang et al. for the prediction of 1-year, 3-year, and 5-year OS were 0.577, 0.587, and 0.603, respectively [14]. Undoubtedly, our 14-CpG-based MSH was better than the five-CpG-based MSH in predicting the OS of HCC patients. Furthermore, we built an MSH-based nomogram to make the MSH more clinically applicable. The C-index and calibration plots exhibited good consistency between the predicted OS and the actual OS, which suggested the accurate prognosis prediction of the MSH-based nomogram.

Fourteen prognosis-related CpGs correspond to TNFRSF19, RIT2, CACNG5, KRTAP11-1, ZNF154, TSPYL5, CPNE4, MS4A7, TBX4, SLC5A4, SULT1C2 CDKL2, S100A8, and SPRR2A. Among all of the CpGs, ZNF154, TSPYL5, CDKL2, and S100A8 have been reported to be associated with HCC. ZNF154, TSPYL5, and CDKL2 were found to be significantly hypermethylated and downregulated in HCC tissues compared to their methylation status in nontumor liver tissues. The methylation of TSPYL5 and CDKL2 could also be used to distinguish HCC tissues from adjacent nontumor tissues [23–27]. Consistent with the above findings, hypermethylation of ZNF154, TSPYL5, and CDKL2 was also found in HCC patients in TCGA. We also found that the methylation levels of these genes in the high-risk patients were higher than those in the low-risk patients, suggesting that ZNF154, TSPYL5, and CDKL2 may play an antitumor role in the development and progression of HCC. In contrast to ZNF154, TSPYL5, and CDKL2, significant hypomethylation of S100A8 was found in HCC tissues when compared to the methylation status of adjacent normal tissues, suggesting shorter OS and progression-free survival (PFS). In addition, overexpression of S100A8 in Huh 7 and MHCC-97H hepatoma cells resulted in increased cell proliferation, migration, and invasion [28]. Furthermore, increased expression of S100A8 and S100A9 promotes the malignant progression of HCC by activating reactive oxygen species (ROS) dependent signaling pathways and inhibiting cell death [29]. In our study, we also found hypomethylated S100A8 in patients with HCC. The methylation level of S100A8 in patients in the high-risk group was lower than that in patients in the low-risk group, indicating the important role of S100A8 in the progression of HCC, which also partly explains why HCC patients with hypomethylation of S100A8 had a shorter OS. Despite the lack of reports about the role of the other 10 genes in HCC, future characterization of these genes may provide new insights into the development and progression of HCC and the discovery of potential novel therapeutic targets for HCC.

Despite our 14-CpG-based signatures showing good performance for the prediction of the prognosis of HCC patients, several limitations of this study should be noted. First, the prognostic value of the MSH was only validated in the internal cohort from TCGA. Other external cohorts with larger sample sizes are still needed to validate our model. Second, although we explored the potential biological functions and pathways of the MSH, more experiments should be conducted to justify the related mechanisms. Finally, as noninvasive “liquid biopsy” has received increasing attention with the potential to revolutionize the diagnosis and treatment of cancer [30], whether these 14 CpGs could be detected in the blood of HCC patients and whether the signature based on these CpGs would still have good prognostic value will need further validation.

In conclusion, we identified 14 differential CpGs that were significantly associated with the OS of HCC patients. The MSH developed by these 14 CpGs showed greater advantage in terms of stability and accuracy in prognosis prediction compared to the predictive ability of AJCC stage and other CpG-based signatures from TCGA. The MSH-based nomogram may help clinicians establish more individualized therapeutic strategies for HCC patients after surgery.

Abbreviations

HCC:	Hepatocellular carcinoma
SNP:	Single-nucleotide polymorphism
TCGA:	The cancer genome atlas
MSH:	Methylation-based signature for HCC
GEO:	Gene Expression Omnibus
OS:	Overall survival
LASSO:	Least absolute shrinkage and selection operator
ROC:	Receiver operating characteristic
AUC:	Area under the ROC curve
GO:	Gene ontology
KEGG:	Kyoto encyclopedia of genes and genomes
DAVID:	Database for annotation, visualization, and integrated discovery
TNFRSF19:	TNF receptor superfamily member 19
RIT2:	Ras-like without CAAX2
CACNG5:	Calcium voltage-gated channel auxiliary subunit gamma 5
KRTAP11-1:	Keratin-associated protein 11–1
ZNF154:	Zinc finger protein 154
TSPYL5:	TSPY-like 5
CPNE4:	Copine 4
MS4A7:	Membrane spanning 4-domains A7
TBX4:	T-box 4
SLC5A4:	Solute carrier family 5 member 4
SULT1C2:	Sulfotransferase family 1C member 2
CDKL2:	Cyclin dependent kinase-like 2
S100A8:	S100 calcium-binding protein A8
SPRR2A:	Small proline rich protein 2A
PFS:	Progression-free survival
ROS:	Reactive oxygen species.

Data Availability

All data in the study were obtained from online databases, including GEO (https://www.ncbi.nlm.nih.gov/geo/) and TCGA (https://cancergenome.nih.gov/).

Ethical Approval

The study was approved by the Ethics Committee of the Shunde Hospital of Southern Medical University.

Informed consent was obtained from the patients before the study.

Conflicts of Interest

All authors declare that there are no conflicts of interest with the contents of this article.

Authors’ Contributions

Hong-ye Jiang and Gang Ning contributed equally to this work.

Supplementary Materials

Supplementary material 1: top 1000 differential CpGs between primary HCC tumors and their corresponding nontumor counterparts of GSE37988. Supplementary material 2: top 1000 differential CpGs between primary HCC tumors and their corresponding nontumor counterparts of GSE57958. Supplementary material 3: top 1000 differential CpGs between primary HCC tumors and their corresponding nontumor counterparts of GSE73003. Supplementary material 4: overlap of each top 1000 differential CpGs from GSE73003, GSE37988, and GSE57958 (426markers). Supplementary material 5: validated differential CpGs of GSE73003, GSE37988, and GSE57958 in HCC patients of TCGA (288 markers). (Supplementary Materials)

References

F. Bray, J. Ferlay, I. Soerjomataram, R. L. Siegel, L. A. Torre, and A. Jemal, “Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries,” CA: A Cancer Journal for Clinicians, vol. 68, no. 6, pp. 394–424, 2018.
View at: Publisher Site | Google Scholar
T. Clark, S. Maximin, J. Meier, S. Pokharel, and P. Bhargava, “Hepatocellular carcinoma: review of epidemiology, screening, imaging diagnosis, response assessment, and treatment,” Current Problems in Diagnostic Radiology, vol. 44, no. 6, pp. 479–486, 2015.
View at: Publisher Site | Google Scholar
P. Tabrizian, G. Jibara, B. Shrager, M. Schwartz, and S. Roayaie, “Recurrence of hepatocellular cancer after resection: patterns, treatments, and prognosis,” Annals of Surgery, vol. 261, no. 5, pp. 947–955, 2015.
View at: Publisher Site | Google Scholar
S. F. Altekruse, J. S. Henley, J. E. Cucinelli, and K. A. McGlynn, “Changing hepatocellular carcinoma incidence and liver cancer mortality rates in the United States,” American Journal of Gastroenterology, vol. 109, no. 4, pp. 542–553, 2014.
View at: Publisher Site | Google Scholar
V. Vedham and M. Verma, “Cancer-associated infectious agents and epigenetic regulation,” Methods in Molecular Biology, vol. 1238, pp. 333–354, 2015.
View at: Publisher Site | Google Scholar
L. Ma, “Epigenetics in hepatocellular carcinoma: an update and future therapy perspectives,” World Journal of Gastroenterology, vol. 20, no. 2, pp. 333–345, 2014.
View at: Publisher Site | Google Scholar
Z. Herceg and A. Paliwal, “Epigenetic mechanisms in hepatocellular carcinoma: how environmental factors influence the epigenome,” Mutation Research/Reviews in Mutation Research, vol. 727, no. 3, pp. 55–61, 2011.
View at: Publisher Site | Google Scholar
B. Wahid, A. Ali, S. Rafique, and M. Idrees, “New insights into the epigenetics of hepatocellular carcinoma,” BioMed Research International, vol. 2017, Article ID 1609575, 16 pages, 2017.
View at: Publisher Site | Google Scholar
M. Nakamura, T. Chiba, K. Kanayama et al., “Epigenetic dysregulation in hepatocellular carcinoma: an up-to-date review,” Hepatology Research, vol. 49, no. 1, pp. 3–13, 2019.
View at: Publisher Site | Google Scholar
N. Yang and P. Sen, “The senescent cell epigenome,” Aging, vol. 10, no. 11, pp. 3590–3609, 2018.
View at: Publisher Site | Google Scholar
X. Wu, X. Yao, Q. Cao et al., “Clinicopathological and prognostic significance of CDH1 hypermethylation in hepatocellular carcinoma: a meta-analysis,” Cancer Management and Research, vol. 11, pp. 857–864, 2019.
View at: Publisher Site | Google Scholar
T. Zhang, G. Guan, T. Chen et al., “Methylation of PCDH19 predicts poor prognosis of hepatocellular carcinoma,” Asia-Pacific Journal of Clinical Oncology, vol. 14, no. 5, pp. e352–e358, 2018.
View at: Publisher Site | Google Scholar
Y. Li, Y. Cai, Y. Li, H. Che, and L. Mao, “Clinical significance and association of GSTP1 hypermethylation with hepatocellular carcinoma: a meta-analysis,” Journal of Cancer Research and Therapeutics, vol. 14, no. 9, pp. S486–S489, 2018.
View at: Publisher Site | Google Scholar
F. Fang, W. Xiaoqing, and S. Tianqiang, “Five-CpG-based prognostic signature for predicting survival in hepatocellular carcinoma patients,” Cancer Biology & Medicine, vol. 15, no. 4, pp. 425–433, 2018.
View at: Publisher Site | Google Scholar
J. Gao, P. W. Kwan, and D. Shi, “Sparse kernel learning with LASSO and Bayesian inference algorithm,” Neural Networks, vol. 23, no. 2, pp. 257–264, 2010.
View at: Publisher Site | Google Scholar
J. Y. Wang and J. J. Tai, “Robust quantitative trait association tests in the parent-offspring triad design: conditional likelihood-based approaches,” Annals of Human Genetics, vol. 73, no. 2, pp. 231–244, 2009.
View at: Publisher Site | Google Scholar
G. Renaud, U. Stenzel, T. Maricic, V. Wiebe, and J. Kelso, “deML: robust demultiplexing of illumina sequences using a likelihood-based approach,” Bioinformatics, vol. 31, no. 5, pp. 770–772, 2015.
View at: Publisher Site | Google Scholar
G. Ning, Y.-L. Huang, L.-M. Zhen et al., “Transcriptional expressions of chromobox 1/2/3/6/8 as independent indicators for survivals in hepatocellular carcinoma patients,” Aging (Albany NY), vol. 10, no. 11, pp. 3450–3473, 2018.
View at: Publisher Site | Google Scholar
X. Zhang, Z.-H. Ye, H.-W. Liang et al., “Down-regulation of miR-146a-5p and its potential targets in hepatocellular carcinoma validated by a TCGA- and GEO-based study,” FEBS Open Bio, vol. 7, no. 4, pp. 504–521, 2017.
View at: Publisher Site | Google Scholar
X. Zhang, W. Tang, G. Chen et al., “An encapsulation of gene signatures for hepatocellular carcinoma, MicroRNA-132 predicted target genes and the corresponding overlaps,” PLoS One, vol. 11, no. 7, Article ID e0159498, 2016.
View at: Publisher Site | Google Scholar
B. Delire and P. Starkel, “The Ras/MAPK pathway and hepatocarcinoma: pathogenesis and therapeutic implications,” European Journal of Clinical Investigation, vol. 45, no. 6, pp. 609–623, 2015.
View at: Publisher Site | Google Scholar
M. Vidyasagar, “Identifying predictive features in drug response using machine learning: opportunities and challenges,” Annual Review of Pharmacology and Toxicology, vol. 55, no. 1, pp. 15–34, 2015.
View at: Publisher Site | Google Scholar
N. Yamada, K. Yasui, O. Dohi et al., “Genome-wide DNA methylation analysis in hepatocellular carcinoma,” Oncology Reports, vol. 35, no. 4, pp. 2228–2236, 2016.
View at: Publisher Site | Google Scholar
X. Qiu, B. Hu, Y. Huang, Y. Deng, X. Wang, and F. Zheng, “Hypermethylation of ACP1, BMP4, and TSPYL5 in hepatocellular carcinoma and their potential clinical significance,” Digestive Diseases and Sciences, vol. 61, no. 1, pp. 149–157, 2016.
View at: Publisher Site | Google Scholar
J. Shen, C. LeFave, I. Sirosh, A. B. Siegel, B. Tycko, and R. M. Santella, “Integrative epigenomic and genomic filtering for methylation markers in hepatocellular carcinomas,” BMC Medical Genomics, vol. 8, no. 1, p. 28, 2015.
View at: Publisher Site | Google Scholar
Y. Zhou, X.-P. Qiu, Z. H. Li et al., “Clinical significance of aberrant cyclin-dependent kinase-like 2 methylation in hepatocellular carcinoma,” Gene, vol. 683, pp. 35–40, 2019.
View at: Publisher Site | Google Scholar
Y. Zheng, Q. Huang, Z. Ding et al., “Genome-wide DNA methylation analysis identifies candidate epigenetic markers and drivers of hepatocellular carcinoma,” Briefings in Bioinformatics, vol. 19, no. 1, pp. 101–108, 2018.
View at: Publisher Site | Google Scholar
K. Liu, Y. Zhang, C. Zhang et al., “Methylation of S100A8 is a promising diagnosis and prognostic marker in hepatocellular carcinoma,” Oncotarget, vol. 7, no. 35, pp. 56798–56810, 2016.
View at: Publisher Site | Google Scholar
J. Nemeth, I. Stein, D. Haag et al., “S100A8 and S100A9 are novel nuclear factor kappa B target genes during malignant progression of murine and human liver carcinogenesis,” Hepatology, vol. 50, no. 4, pp. 1251–1262, 2009.
View at: Publisher Site | Google Scholar
C. Bettegowda, M. Sausen, R. J. Leary et al., “Detection of circulating tumor DNA in early- and late-stage human malignancies,” Science Translational Medicine, vol. 6, no. 224, p. 224ra24, 2014.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2020 Hong-ye Jiang 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.

PDF Download Citation

Download other formats

Order printed copies

Views

573

Downloads

1159

Citations