- Research
- Open access
- Published:
Screening for aberrantly methylated and differentially expressed genes in nonalcoholic fatty liver disease of hepatocellular carcinoma patients with cirrhosis
World Journal of Surgical Oncology volume 20, Article number: 364 (2022)
Abstract
Background
Nonalcoholic fatty liver disease (NAFLD) as the leading chronic liver disease worldwide causes hepatic fibrosis, cirrhosis and hepatocellular carcinoma (HCC). The aim of this study was to find potential aberrantly methylated and differentially expressed genes in NAFLD of HCC patients with cirrhosis.
Methods
DNA methylation data, mRNA expression data, and the corresponding clinical information of HCC were downloaded from the Cancer Genome Atlas (TCGA, tissue sample) database. HCC patients with cirrhosis were divided into two groups according to the presence of NAFLD. The differentially expressed genes (DEGs) and differentially methylated genes (DMGs) were obtained.
Results
By overlapping 79 up-regulated genes and 1020 hypomethylated genes, we obtained 5 hypomethylated-highly expressed genes (HypoHGs). By overlapping 365 down-regulated genes and 481 hypermethylated genes, we identified 13 hypermethylated-lowly expressed genes (Hyper-LGs). Survival analysis of these 18 MDEGs indicated that the expression of DGKK and HOXD9 was significantly correlated with the overall survival time of NAFLD patients.
Conclusions
We identified several candidate genes whose expressions were regulated by DNA methylation of NAFLD of HCC with cirrhosis, which may provide a new field in understanding the clinical pathological mechanism of NAFLD of HCC with cirrhosis.
Background
Nonalcoholic fatty liver disease (NAFLD) is emerging as the leading chronic liver disease worldwide [18]. NAFLD, including nonalcoholic steatohepatitis (NASH), causes hepatic fibrosis, cirrhosis, and hepatocellular carcinoma (HCC) [9]. A population-based cohort (n = 8,892) study indicated that NAFLD patients had significantly increased overall cancer incidence, especially HCC, compared to the controls, adding evidence of cancer risk in patients with NAFLD [28]. The characteristics of patients with HCC, secondary to NAFLD, are older age, large tumors due to late diagnosis, often without cirrhosis, and high prevalence of the metabolic syndrome components, leading to an increased mortality rate [2].
HCC, the most prevalent form of liver cancer, is one of the most frequently occurring malignancy around the world [30]. HCC is an insidious tumor that is often diagnosed in the later stage of life [25]. Moreover, the HCC recurrence occurred in patients with hepatitis D virus when hepatitis B virus recurrence developed [3]. It is reported that Hepatitis B and C virus-related chronic liver disease have exceptionally high risk of HCC development [15, 29]. The epidemiologic data also support that HCC incidence is highest in regions with the higher incidence of HBV. In addition, cirrhosis and smoking are important risk factors for HCC. NAFLD has become one of the leading etiologies for HCC [24]. It is reported that NAFLD-related HCC tends to occur in older individuals and tends to be diagnosed at a later stage [38]. In addition, lifestyle, social conditions, and ethnicity may contribute to the incidence of NAFLD-related HCC [24]. Given the high incidence of NAFLD and its close correlation with the occurrence of HCC, it is of great significance to reveal the pathogenesis of NAFLD and the possible mechanism of its transformation to HCC.
DNA methylation has emerged as an important epigenetic modification and plays a key role in the regulation of gene expression and genome stability [5]. Alterations in DNA methylation in promoter regions has been reported to contribute to the occurrence, development, and prognosis of multiple cancers [7]. Accumulating evidence indicates that the dynamic patterns of DNA methylation are closely associated with the development, diagnosis, and prognosis of liver cancer [11, 22]. It has been reported that DNA methylation silenced BCLB gene expression participates in the progression of HCC, indicating its therapeutic implications for HCC [16]. In alcohol-related HCC, retinol metabolism genes and serine hydroxymethyltransferase 1 are epigenetically regulated via promoter DNA methylation [33]. Besides, Kuo et al. found that IRAK3 methylation was associated with the tumor stage and poor prognosis of HCC patients [10]. Therefore, DNA methylation alteration may play a coordinating role in promoting the carcinogenesis and progression of liver cancer.
In the present study, we aimed to evaluate methylation changes specific to NAFLD of HCC with cirrhosis that could be used as tools in the clinical setting for prognostic assessment of patients. To achieve this goal, we used publicly available microarray data to determine the differentially expressed genes (DEGs) and differentially methylated genes (DMGs) between non-NAFLD and NAFLD HCC patients with cirrhosis. Methylated differentially expressed genes (MDEGs) were obtained by overlapping DEGs and DMGs. Comprehensive and advanced bioinformatics analysis of existing microarray data can reveal more reliable and precise disease-related results.
Materials and methods
Data collection
The Cancer Genome Atlas (TCGA) database contains clinical information from tumor tissues, normal tissue, or blood samples from hundreds of patients with specific cancers on a large scale, which is used for comprehensive genomic data analysis and integration analysis. Currently, the TCGA database covers global information with more than 30 kinds of cancers, including genome variation, gene expression, copy number, genotypes, DNA methylation, and exon sorting. For HCC, a total of 377 patients with HCC were included in the TCGA database, including clinical data of 377 patients, RNA sequencing data of 371 patients, and methylation array data of 377 patients. The DNA methylation data, mRNA expression data, and the corresponding clinical information of HCC were downloaded from the TCGA database (involving tissue sample). According to the fibrosis score, 142 HCC patients with cirrhosis were selected. Then, 142 HCC patients with cirrhosis were divided into two groups according to the presence of NAFLD samples with both mRNA expression data and DNA methylation data included. Consequently, 126 non-NAFLD and 12 NAFLD HCC patients with cirrhosis were included in this study. Clinical information of patients included in this study is indicated in Table 1. Student’s t-test was performed for continuous variable, and chi-square test was performed for categorical variable.
Differential analysis of genes
Firstly, the difficultly detected mRNAs with read count value = 0 in more than 20% samples were filtered and deleted. Based on the read count of each sample, the DEGs in NAFLD compared to non-NAFLD were determined by DESeq2 with p-value <0.05 and |log2 fold change (FC)| > 1. With R package “pheatmap”, hierarchical clustering analysis of top 100 DEGs was conducted. ClusterProfiler (version 3.10.1) was applied to perform GO and KEGG pathway enrichment analysis of DEGs. The threshold was set at p-value <0.05.
Differential analysis of DNA methylation
COHCAP package in R was used to analyze differentially methylated sites between NAFLD and non-NAFLD. Methylated sites with a β value = N/A in more than 20% samples were filtered and deleted. The threshold of differentially methylated sites was set as |Δβ| > 0.1 and p-value <0.05. Then, DMGs and MDEGs were obtained as well. Overlapping down-regulated and hypermethylation genes were identified as hypermethylated-lowly expressed genes (Hyper-LGs). Similarly, overlapping up-regulated and hypomethylation genes were considered hypomethylated-highly expressed genes (HypoHGs).
Proteinprotein interaction (PPI) network and survival analysis
A PPI network of Hyper-LGs and Hyper-LGs was built using the STRING database. Consequently, the PPI network was visualized by Cytoscape. In addition, to further investigate the prognostic value of MDEGs, the clinical data of samples were downloaded from TCGA to obtain survival information. Survival and SurvMiner packages were used for survival analysis after the integration of survival information and expression matrix.
Expression validation of Hyper-LGs and HypoHGs by RT-PCR
To validate the expression of Hyper-LGs and HypoHGs, in vitro RT-PCR was performed. A total of 7 NAFLD of HCC patients with cirrhosis and 6 non-NAFLD of HCC patients with cirrhosis were enrolled in this study. The blood samples from the above individuals were collected for RT-PCR. GAPDH and ACTB were used as internal reference. The study was approved by the ethics committee of the Second People’s Hospital of Liaocheng (2022-44). In addition, all individuals provided the informed consent of the patients and their families.
Results
Identification of DEGs
With p-value <0.05 and |log2FC| > 1, a total of 444 DEGs were identified in NAFLD, of which 79 genes (17.8%) were up-regulated while others were down-regulated. The heatmap of top 100 DEGs is shown in Fig. 1. Among them, PTCRA and MUC5B were the most up-regulated and down-regulated DEGs (Table 2). GO analysis indicated that DEGs were significantly enriched in detoxification of copper ion (p-value = 8.09E−09), glutamatergic synapse (p-value = 4.39E−04), acetylcholine receptor regulator activity (p-value = 1.25E−04), and neurotransmitter receptor regulator activity (p-value = 1.25E−04) (Fig. 2A-C). According to the KEGG pathway enrichment analysis, several pathways, including mineral absorption (p-value = 3.28E−08), protein digestion and absorption (p-value = 2.80E−04), and neuroactive ligand-receptor interaction (p-value = 3.82E−04) were significantly enriched (Fig. 2D).
Identification of DMGs
Following the pre-treatment of methylation data without a β value, 395,042 sites were obtained. With the screening criteria of |Δβ| >0.1 and p-value <0.05, a total of 3082 differentially methylated sites and 1501 DMGs (481 hypermethylated genes and 1020 hypomethylated genes) were obtained. The Manhattan plot of these differentially methylated sites is presented in Fig. 3. By overlapping 79 up-regulated genes and 1020 hypomethylated genes, we obtained 5 HypoHGs (HOXD9, RAI2, ADPRHL1, C12orf42 and PCDHB16). By overlapping 365 down-regulated genes and 481 hypermethylated genes, we identified 13 Hyper-LGs (EPCAM, GNG4, SLFN13, USH1C, SPINT1, SLC39A4, LYZ, SPARCL1, DGKK, WNK2, DNAH9, STRA8 and ST8SIA3).
PPI network and survival analysis
With STRING online database and Cytoscape software, the PPI networks were constructed. PPI networks included 18 proteins, which consisted of two subgroups: 6 proteins with strong connections with others and 12 separated proteins (Fig. 4). EPCAM and SPARCL1 were two hub genes in PPI network. Survival analysis of these 18 MDEGs was performed, and the results indicated that the expression of DGKK (p-value = 0.038) and HOXD9 (p-value <0.001) was significantly correlated with the overall survival time of NAFLD patients (Fig. 5).
Expression validation of Hyper-LGs and HypoHGs by RT-PCR
To validate the expression of 5 Hyper-LGs (GNG4, EPCAM, SPARCL1, DGKK, and SLC39A4) and 1 HypoHGs (HOXD9), in vitro RT-PCR was performed in blood samples from 7 NAFLD of HCC patients with cirrhosis and 6 non-NAFLD of HCC patients with cirrhosis (Fig. 6). The clinical information of these individuals is listed in Table 3. Compared to non-NAFLD of HCC patients with cirrhosis, GNG4, EPCAM, SPARCL1, DGKK, and SLC39A4 were down-regulated and HOXD9 was up-regulated in NAFLD of HCC patients with cirrhosis. The expression trend of these genes was consistent with the bioinformatics analysis.
Discussion
It has been widely recognized that aberrant DNA methylation is significantly associated with HCC. In the present study, using the TCGA database, we searched for the candidates of cancer-related genes whose expressions are regulated by DNA methylation of NAFLD of HCC with cirrhosis.
High levels of GNG4 were reported in primary gastric cancer tissues as well as liver metastatic lesions, which were associated with short overall survival (OS) and the likelihood of liver recurrence [32]. It has been reported that GNG4 is up-regulated in colorectal carcinoma and liver metastases from colorectal carcinoma tissues, which is related to the OS and tumor-free survival of colorectal carcinoma patients [41]. Increased GNG4 expression is related to the poor prognosis and hypoxic microenvironment in lung adenocarcinoma [42]. Pal et al. reported that the promoter region of GNG4 was significantly hypermethylated and that its transcript level was significantly down-regulated in glioblastoma and renal cell carcinoma [23]. Mao et al. indicated that GNG4 was hypermethylated and its mRNA expression was significantly decreased in breast cancer [20]. In this study, GNG4 was significantly hypermethylated and one of the top 10 significantly down-regulated genes in NAFLD HCC patients with cirrhosis, suggesting that GNG4 may be a potential clinical therapeutic target for NAFLD HCC patients with cirrhosis by inhibiting tumor metastasis.
Epithelial cell adhesion molecule (EpCAM) is a type I transmembrane glycoprotein, acting as a Ca2+-independent homophilic cell adhesion molecule [13]. EpCAM has been reported to be involved in malignant proliferation, invasion, metastasis, and tumor recurrence [1]. Elevated EpCAM has been detected in various human tumors, including HCC [37]. The expression of EpCAM was significantly associated with inflammation in HBV infection, which serves as an early biomarker for HCC [4]. EpCAM has been demonstrated to be one of the targets of chemoresistance in human hepatocellular carcinoma cell lines [12]. At present, EPCAM was a hub gene in PPI network and one of the top 10 DEGs, as well as Hyper-LG, in NAFLD HCC patients with cirrhosis, which indicated that EPCAM may exert a momentous role in inflammatory response and tumorigenesis in NAFLD of HCC with cirrhosis. Maybe, EPCAM can be considered as a target of drug in the treatment of NAFLD of HCC with cirrhosis.
SPARCL1 has been reported to be expressed in confluent endothelial cells and is one of the signature genes for tumor angiogenesis [34]. Zhang et al. indicated that SPARCL1 was a prognostic biomarker in colorectal cancer and likely played a more significant role in the metastasis of primary colorectal cancer cells to normal liver tissues [39]. Liu et al. demonstrated that SPARCL1 was highly up-regulated in adipose tissue and played a role in exacerbating NASH progression in a mouse model of NASH [14]. Gao et al. suggested that SPARCL1 with AUC greater than 90% could be used as a diagnostic biomarker for liver cancer [6]. In this analysis, SPARCL1 was a hub gene in PPI network and Hyper-LG in NAFLD HCC patients with cirrhosis, which reminds us to focus on the role of tumor angiogenesis and metastasis of SPARCL1 in NAFLD of HCC with cirrhosis.
The DGKK gene (OMIM *300837), located on chromosome Xp11.22, encodes the diacylglycerol kinase kappa [8]. This enzyme is involved in the down-regulation of diacylglycerol signaling since it phosphorylates diacylglycerol, converting it to phosphatidic acid [27]. Genetic variants in DGKK have been strongly associated with risk for hypospadias [26]. In addition, up-regulated DGKK proteins are detected in HCC tumor tissue samples from mice treated with high-dose ascorbate [40]. Apart from this study, DGKK has rarely been reported in HCC. It is noted that DGKK was significantly correlated with overall survival time of NAFLD patients, which indicated that DGKK could be regarded as a potential prognostic marker molecule for NAFLD of HCC with cirrhosis patients.
Numerous studies have reported the highly expressed HOXD9 in HCC. It has been demonstrated that HOXD9 was strongly expressed and functioned as an oncogene to promote epithelial-mesenchymal transition and cancer metastasis in HCC [19]. Long et al. indicated that high expression levels of HOXD9 were relevant to a poor prognosis in HCC patients [17]. Over expressed HOXD9 was detected in HCC patients with microvascular invasion compared to patients without microvascular invasion and associated with poorer prognosis [35]. In addition, high level of HOXD9 has been closely linked to metastasis rate and poor prognosis in cervical cancer patients [36]. Previous studies have revealed that HOXD9 promoter methylation is higher in tumors than in healthy tissue and that DNA methylation levels correlate with the expression of HOXD9 mRNA and protein in malignant melanoma and glioma [21, 31]. Consistent with previous studies, HOXD9 was significantly associated with prognosis of NAFLD patients in this analysis, indicating that HOXD9 may be an effective clinical therapeutic target for NAFLD of HCC with cirrhosis.
Conclusions
In conclusion, 5 HypoHGs (HOXD9, RAI2, ADPRHL1, C12orf42, and PCDHB16) and 13 Hyper-LGs (EPCAM, GNG4, SLFN13, USH1C, SPINT1, SLC39A4, LYZ, SPARCL1, DGKK, WNK2, DNAH9, STRA8, and ST8SIA3) were identified in NAFLD patients HCC with cirrhosis. Among them, EPCAM and SPARCL1 were identified as two hub genes and DGKK and HOXD9 were significantly correlated with prognosis. These genes may be involved in the development of NAFLD-related HCC with cirrhosis, which may be used in the clinical therapeutic targets. However, there are limitations to our study. Firstly, the sample size of the RT-PCR is small. Larger numbers of tissue and blood samples are further needed. Secondly, the potential molecular mechanism of identified genes is needed to explore in the animal models or cell experiment.
Availability of data and materials
The data supporting the conclusions of this article is included within the article and is available from the corresponding author on reasonable request.
Abbreviations
- NAFLD:
-
Nonalcoholic fatty liver disease
- NASH:
-
Nonalcoholic steatohepatitis
- HCC:
-
Hepatocellular carcinoma
- BCLB:
-
Bcl 2 like protein 10
- DEG:
-
Differentially expressed gene
- DMG:
-
Differentially methylated gene
- FC:
-
Fold change
- MDEG:
-
Methylated differentially expressed gene
- Hyper-LG:
-
Hypermethylated, lowly expressed gene
- HypoHG:
-
Hypomethylated, highly expressed gene
- PPI:
-
Proteinprotein interaction
- OS:
-
Overall survival
- EpCAM:
-
Epithelial cell adhesion molecule
References
Abdelgawad IA. Epithelial cell adhesion molecule mRNA can be a potential marker to predict metastasis in hepatocellular carcinoma patients. Asian Pac J Cancer Prev. 2020;21:861–6.
Anastasopoulos NT, Lianos GD, Tatsi V, Karampa A, Goussia A, Glantzounis GK. Clinical heterogeneity in patients with non-alcoholic fatty liver disease-associated hepatocellular carcinoma. Expert Rev Gastroenterol Hepatol. 2020;14:1025–33.
Baskiran A, Akbulut S. Effect of HBV-HDV co-infection on HBV-HCC co-recurrence in patients undergoing living donor liver transplantation. Hepatol Int. 2020;14:869–80.
Chen Y, Hao X, Sun R, Wei H, Tian Z. Natural killer cell-derived interferon-gamma promotes hepatocellular carcinoma through the epithelial cell adhesion molecule-epithelial-to-mesenchymal transition axis in hepatitis B virus transgenic mice. Hepatology. 2019;69:1735–50.
Corso-Díaz X, Jaeger C, Chaitankar V, Swaroop A. Epigenetic control of gene regulation during development and disease: a view from the retina. Prog Retin Eye Res. 2018;65:1–27.
Gao S, Gang J, Yu M, Xin G, Tan H. Computational analysis for identification of early diagnostic biomarkers and prognostic biomarkers of liver cancer based on GEO and TCGA databases and studies on pathways and biological functions affecting the survival time of liver cancer. BMC Cancer. 2021;21:791.
Hao XY, Li AQ, Shi H, Guo TK, Shen YF, Deng Y, et al. A novel DNA methylation-based model that effectively predicts prognosis in hepatocellular carcinoma. Biosci Rep. 2021;41:BSR20203945.
Hozyasz KK, Mostowska A, Kowal A, Mydlak D, Tsibulski A, Jagodzinski PP. Further evidence of the Association of the Diacylglycerol Kinase Kappa (DGKK) gene with hypospadias. Urol J. 2018;15:272–6.
Kanda T, Goto T, Hirotsu Y, Masuzaki R, Moriyama M, Omata M. Molecular mechanisms: connections between nonalcoholic fatty liver disease, Steatohepatitis and Hepatocellular Carcinoma. Int J Mol Sci. 2020;21:1525.
Kuo CC, Shih YL, Su HY, Yan MD, Hsieh CB, Liu CY, et al. Methylation of IRAK3 is a novel prognostic marker in hepatocellular carcinoma. World J Gastroenterol. 2015;21:3960–9.
Laugsand EA, Skorpen F, Kaasa S, Sabatowski R, Strasser F, Fayers P, et al. Genetic and non-genetic factors associated with constipation in cancer patients receiving opioids. Clin Transl Gastroenterol. 2015;6:e90.
Li Y, Farmer RW, Yang Y, Martin RC. Epithelial cell adhesion molecule in human hepatocellular carcinoma cell lines: a target of chemoresistence. BMC Cancer. 2016;16:228.
Lima LDP, Machado CJ, Rodrigues J, Vasconcellos LS, Junior EP, Vidigal PVT, et al. Immunohistochemical coexpression of epithelial cell adhesion molecule and alpha-fetoprotein in hepatocellular carcinoma. Can J Gastroenterol Hepatol. 2018;2018:5970852.
Liu B, Xiang L, Ji J, Liu W, Chen Y, Xia M, et al. Sparcl1 promotes nonalcoholic steatohepatitis progression in mice through upregulation of CCL2. J Clin Invest. 2021;131:e144801.
Liu JH, Chen PW, Asch SM, Busuttil RW, Ko CY. Surgery for hepatocellular carcinoma: does it improve survival? Ann Surg Oncol. 2004;11:298–303.
Liu X, Hu X, Kuang Y, Yan P, Li L, Li C, et al. BCLB, methylated in hepatocellular carcinoma, is a starvation stress sensor that induces apoptosis and autophagy through the AMPK-mTOR signaling cascade. Cancer Lett. 2017;395:63–71.
Long J, Zhang L, Wan X, Lin J, Bai Y, Xu W, et al. A four-gene-based prognostic model predicts overall survival in patients with hepatocellular carcinoma. J Cell Mol Med. 2018;22:5928–38.
Loomba R, Friedman SL, Shulman GI. Mechanisms and disease consequences of nonalcoholic fatty liver disease. Cell. 2021;184:2537–64.
Lv X, Li L, Lv L, Qu X, Jin S, Li K, et al. HOXD9 promotes epithelial-mesenchymal transition and cancer metastasis by ZEB1 regulation in hepatocellular carcinoma. J Exp Clin Cancer Res. 2015;34:133.
Mao XH, Ye Q, Zhang GB, Jiang JY, Zhao HY, Shao YF, et al. Identification of differentially methylated genes as diagnostic and prognostic biomarkers of breast cancer. World J Surg Oncol. 2021;19:29.
Marzese DM, Scolyer RA, Huynh JL, Huang SK, Hirose H, Chong KK, et al. Epigenome-wide DNA methylation landscape of melanoma progression to brain metastasis reveals aberrations on homeobox D cluster associated with prognosis. Hum Mol Genet. 2014;23:226–38.
McPherson S, McMullin MF, Mills K. Epigenetics in myeloproliferative neoplasms. J Cell Mol Med. 2017;21:1660–7.
Pal J, Patil V, Mondal B, Shukla S, Hegde AS, Arivazhagan A, et al. Epigenetically silenced GNG4 inhibits SDF1α/CXCR4 signaling in mesenchymal glioblastoma. Genes Cancer. 2016;7:136–47.
Plaz Torres MC, Bodini G, Furnari M, Marabotto E, Zentilin P, Strazzabosco M, et al. Surveillance for hepatocellular carcinoma in patients with non-alcoholic fatty liver disease: universal or selective? Cancers. 2020;12:1422.
Satilmis B, Sahin TT. Hepatocellular carcinoma tumor microenvironment and its implications in terms of anti-tumor immunity: future perspectives for new therapeutics. J Gastrointest Cancer. 2021;52:1198–205.
Shen J, Liu B, Sinclair A, Cunha G, Baskin LS, Choudhry S. Expression analysis of DGKK during external genitalia formation. J Urol. 2015;194:1728–36.
Shionoya T, Usuki T, Komenoi S, Isozaki T, Sakai H, Sakane F. Distinct expression and localization of the type II diacylglycerol kinase isozymes δ, η and κ in the mouse reproductive organs. BMC Dev Biol. 2015;15:6.
Simon TG, Roelstraete B, Sharma R, Khalili H, Hagström H, Ludvigsson JF. Cancer risk in patients with biopsy-confirmed nonalcoholic fatty liver disease: a population-based cohort study. Hepatology. 2021;74:2410–23.
Stravitz RT, Heuman DM, Chand N, Sterling RK, Shiffman ML, Luketic VA, et al. Surveillance for hepatocellular carcinoma in patients with cirrhosis improves outcome. Am J Med. 2008;121:119–26.
Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71:209–49.
Tabuse M, Ohta S, Ohashi Y, Fukaya R, Misawa A, Yoshida K, et al. Functional analysis of HOXD9 in human gliomas and glioma cancer stem cells. Mol Cancer. 2011;10:60.
Tanaka H, Kanda M, Miwa T, Umeda S, Sawaki K, Tanaka C, et al. G-protein subunit gamma-4 expression has potential for detection, prediction and therapeutic targeting in liver metastasis of gastric cancer. Br J Cancer. 2021;125:220–8.
Udali S, Guarini P, Ruzzenente A, Ferrarini A, Guglielmi A, Lotto V, et al. DNA methylation and gene expression profiles show novel regulatory pathways in hepatocellular carcinoma. Clin Epigenetics. 2015;7:43.
Wang H, Feng C, Lu M, Zhang B, Xu Y, Zeng Q, et al. Integrative single-cell transcriptome analysis reveals a subpopulation of fibroblasts associated with favorable prognosis of liver cancer patients. Transl Oncol. 2021;14:100981.
Wang J, Ding ZW, Chen K, Liu YZ, Li N, Hu MG. A predictive and prognostic model for hepatocellular carcinoma with microvascular invasion based TCGA database genomics. BMC Cancer. 2021;21:1337.
Wen D, Wang L, Tan S, Tang R, Xie W, Liu S, et al. HOXD9 aggravates the development of cervical cancer by transcriptionally activating HMCN1. Panminerva Med. 2020. https://doi.org/10.23736/S0031-0808.20.03911-7.
Yamashita T, Forgues M, Wang W, Kim JW, Ye Q, Jia H, et al. EpCAM and alpha-fetoprotein expression defines novel prognostic subtypes of hepatocellular carcinoma. Cancer Res. 2008;68:1451–61.
Younossi ZM, Otgonsuren M, Henry L, Venkatesan C, Mishra A, Erario M, et al. Association of nonalcoholic fatty liver disease (NAFLD) with hepatocellular carcinoma (HCC) in the United States from 2004 to 2009. Hepatology. 2015;62:1723–30.
Zhang T, Yuan K, Wang Y, Xu M, Cai S, Chen C, et al. Identification of candidate biomarkers and prognostic analysis in colorectal cancer liver metastases. Front Oncol. 2021;11:652354.
Zhang X, Liu T, Li Z, Feng Y, Corpe C, Liu S, et al. Hepatomas are exquisitely sensitive to pharmacologic ascorbate (P-AscH(-)). Theranostics. 2019;9:8109–26.
Zhao H, Sheng D, Qian Z, Ye S, Chen J, Tang Z. Identifying GNG4 might play an important role in colorectal cancer TMB. Cancer Biomark. 2021;32:435–50.
Zhou B, Zhu W, Yuan S, Wang Y, Zhang Q, Zheng H, et al. High GNG4 expression is associated with poor prognosis in patients with lung adenocarcinoma. Thorac Cancer. 2021;13:369–79.
Acknowledgements
Not applicable.
Funding
Not applicable.
Author information
Authors and Affiliations
Contributions
GZ and YH contributed to the study conception and design, material preparation, data collection, and analysis. The first draft of the manuscript was written by GZ and YH, and all authors commented on previous versions of the manuscript, read, and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
The study was approved by the ethics committee of the Second People’s Hospital of Liaocheng (2022-44). In addition, all individuals provided the informed consent of the patients and their families.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Zhang, G., Hou, Y. Screening for aberrantly methylated and differentially expressed genes in nonalcoholic fatty liver disease of hepatocellular carcinoma patients with cirrhosis. World J Surg Onc 20, 364 (2022). https://doi.org/10.1186/s12957-022-02828-3
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12957-022-02828-3