[BACK]
BIOCELL
DOI:10.32604/biocell.2022.020198
images
Article

Immune-related DNA methylation signature associated with APLN expression predicts prognostic of hepatocellular carcinoma

FEIFEI TIAN1, HUAN HU1, DI WANG1, HUAN DING1, QINGJIA CHI2, HUAPING LIANG3 and WENLI ZENG4,*

1School of Life Science and Engineering, Southwest Jiaotong University, Chengdu, 610000, China
2Department of Engineering Structure and Mechanics, School of Science, Wuhan University of Technology, Wuhan, 430070, China
3State Key Laboratory of Trauma, Burns and Combined Injury, Research Institute of Surgery, Daping Hospital, The Army Medical University, Chongqing, 400042, China
4Renmin Hospital of Wuhan University Nursing Department, Wuhan, 430070, China
*Address correspondence to: Wenli Zeng, Zengwenli600@163.com
Received: 18 November 2021; Accepted: 11 February 2022

Abstract: This study used transcriptome and epigenetic data to predict the prognosis of immune-related genes (IRGs) Apelin (APLN) in patients with hepatocellular carcinoma (HCC). The TCGA database has gene expression and clinical data for HCC. And DNA methylation 450 k data for HCC was download from the University of California Santa Cruz (UCSC) Xena browser. Performing clinical and prognostic analysis of APLN expression, results show that APLN is highly expressed in tumor samples. And it has an increasing trend with the development of clinical stage and T stage. To explore the prognostic role of APLN, the Immune-related DNA methylation (DNAm) sites associated with APLN analyzed by bioinformatics. Univariate COX screened the methylation sites that are related to both APLN and survival. The risk score related to methylation site signature was determined according to their least absolute shrinkage and selection operator (LASSO) coefficients. Then the patients were divided into high-risk groups and low-risk groups. Significant differences in overall survival (OS) were found in the training cohort. Nomogram shows that APLN or methylation signature can effectively predict the prognosis of HCC patients. In summary, APLN may be a diagnostic and prognostic marker for HCC.

Keywords: TIME; Immunity; Prognosis; Bioinformatics

Abbreviations

HCC: hepatocellular carcinoma
K-M: Kaplan–Meier
LASSO: least absolute shrinkage and selection operator
TCGA: the cancer genome atlas
WGCNA: weighted gene co-expression network analysis
DNAm: DNA methylation
UCSC: University of California, Santa Cruz

Introduction

Hepatocellular carcinoma is one of the common primary tumors, causing more than 700,000 deaths each year. Tobacco, heredity, epigenetic factors will affect its development. Although many detection and treatment methods of computer tomography (CT) and emerging carbon ion radiation have made good progress (Malouff et al., 2020). HCC patients are prone to recurrence and metastasis of hemangioma after surgery or ablation, causing the poor prognosis (Korean Liver Cancer Association and National Cancer Center, 2019; Brown et al., 2019; Muto et al., 2014). Currently, the study’s diagnostic markers are limited and only for some populations (Rawat et al., 2018). Therefore, it is of great significance to study specific diagnostic biomarkers to predict the survival and diagnosis of HCC.

Epigenetic regulation does not change DNA sequence but changes the transcription mechanism and affects the tumor’s progress (Sharma et al., 2020; Dzobo et al., 2020; Baylin and Jones, 2011). Among them, DNA methylation (5-methylcytosine DNAm, m5C) and RNA methylation (6-methyladenine RNAm, m6A) are two crucial nucleic acid modifications, which are epigenetic important regulation mechanism (Liu and He, 2020). DNAm modification affects gene transcription (alternative splicing or translation), leading to abnormal gene expression and causing diseases. DNA methylation sites were served as prognostic markers in colorectal cancer and kidney cancer with high sensitivities, e.g., cg10673833 and cg04448376 (Luo et al., 2020; Wang and Zhang, 2020; Ibrahim et al., 2011; Ashokan et al., 2021). Compared with a single omics biomaker to explore HCC prognosis, the strategy combining DNA methylation data and transcriptomics’ data is believed to be more promising.

The APLN gene locates on chromosome Xq25-q26.1 and can encode 77 amino acid pro-peptide Apelin (Xu et al., 2017). APLN may be a potential marker in ovarian cancer, colorectal cancer, kidney cancer and other cancers (Chen et al., 2021; Ran et al., 2021; Bai et al., 2020; Zuurbier et al., 2017). Not only can APLN gene activate the pathway of APLN/APLNR to promote the HCC development, but also it activates different signal pathways or plays anti-apoptotic and anti-inflammatory functions by promoting tumor angiogenesis (Cabiati et al., 2021; Wang et al., 2020b). Many factors can affect the development of HCC. Immature and unstable tumor vasculature is associated with abnormal tumor microenvironment, which is liable to pose drug resistance during the anti-cancer treatment. The new peptide apelin and its receptor can enhance immune efficacy by inducing the morphological and functional maturation of tumor blood vessels (Kidoya et al., 2012). APLN/APLNR will be up-regulated during tumor angiogenesis. And APLN stimulates the proliferation, migration and invasion of colon adenocarcinomas (Picault et al., 2014; Podgorska et al., 2018). In addition, increased expression of APLN has been observed in muscle-invasive bladder cancer and is associated with poor clinical outcomes (Yang et al., 2019). Studies found that epigenetic modification and Expression of APLN can affect the progression of lung cancer (Miller et al., 2018). The degree of CpG island DNAm in APLN is related to the infiltration of pneumonia cells. And the DNA methylation sites in important regions can regulate gene function (Mishra et al., 2015). APLN can inhibit tumor progression by regulating Tregs immune cells’ expression and VCAM-1 (Amoozgar et al., 2019). We studied the expression changes of APLN in the development of HCC and its ability to predict prognosis. The methylation sites signatures in epigenetics are used to predict the viability of HCC, which is closely related to APLN.

Our team determined the molecular interaction mechanisms and biomarkers of certain diseases through bioinformatics methods (Tian et al., 2015). This study is based on the transcriptome, genetic data of the TCGA dataset and the UCSC dataset to study the prognosis of HCC. First, we used GEPIA and TIMER to screen out suitable genes. Then we screened out APLN gene through PubMed. CIBERSORT’s deconvolution can estimate the abundance of 22 kinds of immune cells and differential immune cells. Immune cells and checkpoints are used as clinical traits for WGCNA. And the univariate Cox, Lasso Cox were utilized to construct the risk score. The Nomogram was used to identify the prognostic ability of DNAm signature. The flow chart of this research is shown in the Fig. 1.

images

Figure 1: Workflow chart for prognostic signature. Based on the UCSC DNAm dataset and TCGA RNA-seq datasets, Univariate COX was used for identifying methylation sites. Differential immune cells were screened out through the results of CIBERSORT. Afterwards, the modules with the highest correlation with immune cells and immune checkpoints were screened through WGCN. And LASSO COX was implemented to select robust survival methylation site characteristics. Finally, the prognostic ability of APLN and APLN-related immune methylation sites is determined.

Materials and Methods

HCC data and preprocess

Gene expression data (N = 422) and clinical data (N = 377) of HCC were download from the TCGA database. The Illumina Infinium HumanDNAm450 data in the USUC database (http://xena.ucsc.edu/) was used for the screening of DNA methylation sites. The Beta value matrix data is composed of genomic Matrix (ROWs: site identifier and COLUMNs: sample). Immune genes were download from the Immport database. We deleted samples with survival time less than or equal to 1 day and DNAm Beta Value missing greater than 70%. GEPIA and ENCORI screened out the genes based on survival and expression data of HCC.

APLN and its expression and survival

In this study, 365 overlapping genes in the TCGA database and the Immport database were used for subsequent analysis. GEPIA (http://gepia2.cancer-pku.cn/) and ENCORI (http://starbase.sysu.edu.cn/panGeneSurvivalExp.php) websites used the “limma” method to analyze overlapping genes. DEGs between the normal group and the cancer group were retained, and genes with a significant difference in survival in the K-M curve (P < 0.05) were retained. Then we conducted literature research through PubMed, and deleted genes not related to immunity or HCC.

Analysis of tumor immune microenvironment

A certain proportion of immune cells in cancer samples were analyzed for immune infiltration. According to the median expression of APLN, tumor samples were divided into high expression group and low expression group. The tumor microenvironment immune cell infiltration was analyzed by CIBERSORT (https://cibersort.stanford.edu/) that based on the deconvolution of gene expression and linear support vector regression. It evaluates the correlation between cell types’ abundance and survival rate in 22 immune cell molecular subtypes. Then we screened out the types of immune cells with significant differences.

Key immune methylation sites

(1) To obtain the methylation sites related to APLN gene, the univariate COX regression model was used. (2) The overall survival (OS) and survival status were jointly-used as dependent variables. And the Univariate COX regression model with the methylation sites were performed again. Then we obtained the methylation sites that are significantly related to prognosis. (3) Venn analysis was implemented to yield the methylation sites related to both the APLN and survival. (4) According to the weighted gene co-expression network analysis, i.e., WGCNA, the essential modules and the immune-related methylation sites were acquired. On the previous step, we used WGCNA to determine the relationship with 3 immune cells and programmed cell death ligand 1 (PD-1), PD-L1, cell Toxic T lymphocyte-associated protein 4 (CTLA-4), the costimulatory molecule B7-1 immune checkpoint had a significant correlation with the key immune methylation sites. Finally, the highest positive/negative correlation module was selected, deemed the most relevant module for immunity.

Methylation sites signature and risk score

We used LAASO COX to analyze the key immune methylation sites. Then we obtained the signature of the methylation site and the LASSO risk coefficient. Kaplan-Meier (K-M) curve and ROC curve verify the accuracy and predictive power of the risk scoring model. 1, 3, and 5 years of AUC value was judged whether the methylation sites signature has diagnostic value. The LASSO COX retains all feature variables through regularization and reduces all parameters in the loss function. Thereby reducing the loss of feature information and preventing data overfitting. As shown in formula (1):

12m[i=1mgθ(A(i)B(i))2+γj=1K|Wj|] (1)

Among them is the L1 regularization term; m is the number of samples; k is the number of parameters; γ is the target of the balanced fitting training and the regularization parameter that keeps the parameter value small.

Prognostic value of APLN and risk score

In the univariate regression analysis, “forestplot” R package analyzed the impact of APLN and risk score on survival. The nomogram analyzed the prognosis of APLN and methylation site signatures. The 1, 3, and 5-year predicted survival status results are compared with the actual 1, 3, and 5-year survival status to determine the prognostic value of APLN and risk score.

Statistical analysis

All data statistical processing uses R version 3.6.1 software. Kaplan-Meier curves were generated uses the “survival” R software package. The R package “timeROC” performs the ROC analysis to assess the predictive performance of APLN gene and the DNA methylation sites signature. The R package “glmnet” implements the Lasso analysis that performs 10-fold cross-validation to obtain the best Lambda value (minimum Lambda value). The ward.D algorithm obtains the relationship between 22 immune cells. For all calculated statistical P values, it is statistically significant to set P < 0.05.

Results

Identify the APLN gene

The screening process of APLN gene is shown in the Fig. 2. TCGA database and Immport database had 365 overlapping genes. The expression of 9 genes in HCC had significant impact on survival. We Deleted genes unrelated to immunity and genes unrelated to HCC.

images

Figure 2: APLN gene screening process. The 365 overlapping genes in the TCGA database and the Immport database are analyzed in the ENCORI and GEPIA2 databases. There are 9 genes with significant differences in gene expression and survival. The above genes were analyzed by TIMER and pubmed. Fially, APLN gene was screened.

The expression of APLN in HCC is related to clinical features and prognosis

The expression of APLN was significantly up-regulated in 371 tumor samples (Fig. 3A; P < 0.05). APLN tends to increase with the extent of the primary tumor and the grade of cancer tissue (Figs. 3B and 3C, P < 0.05). The Kaplan-Meier curve in high expression group shown a lower survival rate than low expression one (Fig. 3D). And the effective area under the ROC curve (AUC1 = 0.64, AUC3 = 0.61, AUC5 = 0.55) shows that APLN predict patient survival greatly (Fig. 3E). The Univariate Cox shown that 66365 DNA methylation sites relate to the expression of APLN and 46208 DNA methylation sites related to survival. A total of 7,229 important DNA methylation sites overlapped in the above two results (Fig. 3F).

images

Figure 3: The expression of APLN in HCC patients and its relationship with clinical and prognosis. (A) Based on the transcription level of APLN in tumor and normal tissues in the TCGA-LIHC data set (371 tumor samples and 50 normal samples); (B) Different disease stages in TCGA; (C) The expression of APLN in the range of the primary tumor; (D) The impact of APLN gene expression on survival in TCGA; (E) ROC analysis of OS prediction for 1, 3, and 5 years. OS, overall survival; (F) Overlapping methylation sites in TCGA related to survival and APLN expression.

The relationship between APLN and infiltrating immune cells

The median expression of APLN divides the samples into high expression group and low expression one. Immune infiltration in the immune microenvironment describes the pattern of immune cells. The proportion of tumor immune cells is significantly different between the high expression and low expression. The highest proportions are T cell CD4 memory resting (23.37%), Macrophages M2 (18.57%), and Macrophages M0 (11.75%) (Fig. 4A). APLN affects expression of CD4 T cell memory resting (P < 0.01), T cell follicular helper (P = 0.001), Tregs (P = 0.002). Three tumor immune cells have a significant impact (Fig. 4D). The relative content of different tumor immune cells in the sample is different (Fig. 4B). Spearman method analysis of 22 tumor immune cells reveal the correlation between 22 tumor immune cells (Fig. 4C).

images

Figure 4: Tumor Infiltrating Immune Cell Analysis by CIBERSORT. (A) Infiltration ratio of 22 immune cells in tumor tissue; (B) Infiltration ratio of 22 immune cells in each tumor sample; (C) Correlation analysis for 22 immune cells; (D) Different immune cells based on APLN expression.

Key methylation sites associated with immune infiltrating cells

According to univariate COX regression, 19754 methylation sites related to APLN and 56,106 methylation sites related to survival were obtained. The WGCNA analyzed overlapping 6797 DNA methylation sites in 365 samples. HCC is closely related to inflammation, immune cells and immune checkpoints may have important therapeutic value. PD-L1 is the main checkpoint in the tumor immune microenvironment (Han et al., 2020). And the PD-1/PD-L1 axis is responsible for cancer immune escape and has a huge impact on cancer treatment (Gao et al., 2009). When the threshold (power) is 26, it is closer to the non-network scale distribution (Figs. 5A and 5B). There are 7 highly correlated modules, of which turquoise and brown modules include 106 and 27 important DNA methylation sites (Fig. 5C).

images

Figure 5: Weighted gene co-expression network analysis of methylation sites (WGCNA). (A) Fitting index β of the law of scale-free network; (B) Module clustering and module-to-module correlation matrix; (C) Correlation analysis between modules and methylation sites, among which turquoise modules and brown modules are most related to immune infiltration.

Construction and verification of risk score

LASSO COX regression analysis of the above 133 methylation sites finally identified 10 important methylation site signatures (cg02316066, cg06194738, cg07311615, cg07570723, cg15629460, cg16987524, cg18437792, cg21860560, cg22580629). (Figs. 6A and 6B). And the best LASSO coefficient (Table 1). The score is determined based on the signature of 10 methylation sites, risk score = cg02316066 × (–2.19) + cg06194738 × 0.07 + cg07311615 × 0.24 + cg07570723 × (0.99) + cg15629460 × 1.99 + cg16987524 × 0.10 + cg18437792 × (0.04) + cg21860560 × 0.98 + cg22580629 × (−0.09). The results of the K-M curve show that the proportion of high risk in the cancer group is high (P = 0.00014 (Fig. 6C)) and good predictive effect (AUC1 = 0.61, AUC3 = 0.64, AUC5 = 0.63) (Fig. 6D). The 365 samples were randomly divided into a verification set and a training set at a ratio of 1:2. The validation set results are consistent with those of 365 samples, and the risk score can predict the prognosis great (Figs. 7A and 7B).

images

Figure 6: LASSO model and verification. (A) LASSO coefficient distribution map of 133 methylation sites; (B) Penalty graph of 133 prognostic methylation sites of the LASSO model, with error bars representing standard errors; (C) Kaplan-Meier curve of OS survival status between the high-risk score group and the low-risk score group; (D) Predict the time-dependent ROC curve analysis of OS for 1, 3, and 5 years based on the risk score.

images

images

Figure 7: Validation set. (A) Kaplan-Meier curve of OS survival status between the high-risk score group and the low-risk score group; (B) The time-dependent ROC curve analysis of OS for 1, 3, and 5 years based on the risk score.

DNA methylation sites signature and the prognostic value of APLN

Univariate analysis was used to screen out the variables related to patient survival. The results showed that APLN expression, risk score, and tumor stage T all had significant effects on patient survival. (Figs. 8A and 8B). Factors with significant differences are used as Nomogram predictors. The results showed that the survival prediction is in good agreement with the actual observations (Fig. 8C).

images

Figure 8: Prognostic analysis of APLN and risk score. (A) One-way COX prognostic regression model based on TCGA and UCSCD data sets; (B) APLN, risk score, stage, T stage combined with nomogram survival prediction; (C) The prediction and actual 1, 3, 5-year survival rate.

Discussion

Conventional treatment for HCC is highly invasive. And early prediction and personalized therapy of targeted immune checkpoints have research significance (Yin et al., 2019; Wang et al., 2021). APLN and DNAm are abnormally expressed in a variety of cancers. But the role of APLN in the immune microenvironment development of HCC is still unknown. APLN is up-regulated in HCC and HCC38/KMUH cells and may be used as a therapeutic target (Lin et al., 2012; Chen et al., 2019). APLN is a carcinogen that is highly expressed and its expression trend is increasing with HCC development. APLN, as a transcriptional target of the Wnt/β-catenin pathway, promotes HCC by activating PI3K/Akt signaling. The pharmacological targeting of APLN strongly inhibits the proliferation of HCC cells and tumor growth, which means that APLN may be a potential drug target for HCC. The gene has significant differences in the expression and clinical characteristics between the tumor and normal tissues. It has the potential prognostic value. The Cox regression model predicts the risk prediction characteristics significantly related to OS. Finally, the prognostic value of the risk prediction was verified by the nomogram.

The study determined that the proportion of total T cells in HCC cells is greater than in healthy tissues. And T cell subsets mainly include T cell CD4 memory resting, follicular helper T cells, etc., which impact the prognosis of HCC (Jia et al., 2015a; Rohr-Udilova et al., 2018; Garnelo et al., 2017). APLN affects three types of immune cells, including T cell CD4 memory resting, follicular helper T cells, and Tregs, highlighting the poor clinical correlation between T cell-mediated immunity and tumors (Li et al., 2021). Both PD-1 and CTLA-4 are related to some subsets of T cells, and mainly affect the development of tumors through PD-1/PD-L1 acting on T cells (Kim and Chen, 2016). PD-1 and PD-L1 interaction can resulting in inactivation of the immune response (Poureau and Metges, 2021). PD-1 is the main checkpoint in TME. Follicular helper T cells can inhibit the expression of Treg cells by affecting the expression of PD-1. And CTLA-4 can also enhance the immunosuppressive activity of Treg (a key mediator of CTLA-4) by regulating follicular helper T cells (Jia et al., 2015b; Hervouet, 1970; Yang et al., 2014; Topalian et al., 2016). The expression of APLN and specific immune cells play an essential role in the tumor progression of HCC (Wang et al., 2020a). APLN can regulate or induce the morphological and functional maturation of tumor blood vessels. And the recruitment of dendritic immune cells helps to enhance the effect of immunotherapy (Kidoya et al., 2012). Closely related to the HCC prognosis, the immune microenvironment is complex and usually heterogeneous with several genomic alterations, therefore making meaningful the integrated investigation of immune microenvironment with the DNA methylation sites of APLN in studies of HCC prognosis.

Not only does tumor-related macrophages promote thrombosis and endocytosis, angiogenesis, invasion and metastasis, but also it inhibits anti-tumor immune responses (Lin et al., 2019). Tumor-associated macrophages are mainly polarized into Macrophages M1 and Macrophages M2. The former plays an anti-tumor effect, while the latter has a tumor-promoting effect. Macrophages M2 mainly induces anti-inflammatory responses by producing IL-10 and (TGF)-β, meanwhile participating in the polarization of Th2 responses to promote tumor growth (Wildes et al., 2020; Biswas and Mantovani, 2010; Wolf-Dennen et al., 2020). Taking a higher percentage in HBV, non-polarized Macrophages M0 also affect the HCC prognosis (Ding et al., 2021; Huo et al., 2020). APLNs induce tumor vascular maturation and promote the infiltration of natural killer T cells into the central region of the tumor. This process enhances the efficacy of cancer dendritic cell-based immunotherapies (Mastrella et al., 2019). APLN inhibits and reduces neutrophil recruitment, thereby reducing the inflammatory and fibrotic response of pancreatitis. It also inhibits macrophage infiltration and monocytes, thereby enhancing the stable phenotype of atherosclerotic plaques (Ioannidis et al., 2010).

The signatures that DNA methylation sites analyzed by LASSO COX prevent data from overfitting. Studies found that the combination of selected central genes and DNAm conditions has an advantage in predicting the survival of HNSCC (Mou et al., 2020; Cao et al., 2020). The human m6A near the mRNA stop codon and the 3’ untranslated region will affect the transcription and function of genes through the promoter DNAm (Ondo et al., 2021). Usually, DNAm inhibits gene expression. APLN and the probe cg05688478 can affect gene expression in human pancreatic islets (Hall et al., 2014). The interaction of genome (gene expression) and epigenome (DNAm level) to control Apelin endothelium is one of the most critical ways for non-triple-negative breast cancer (Wu et al., 2020). DNA methylation silences the APLN expression, leading to increased insulin secretion in women (Hall et al., 2014). In APLN, rs3761581 G and rs2235312 T are relative to the degree of CpG island. Low levels of apelin-13 will cause greater methylation of CpG islands in the 5’UTR, which contributes to the pathophysiology of HA lung water. Abnormal methylation may prevent apelin from binding to hypoxia response elements, thereby inhibiting APLN transcription (Mishra et al., 2015). In downstream of DNA damage, the increased methylation of the APLN promoter may lead to reductions in protective signaling of the apeliner energy system. Then it lead to the pulmonary edema observed after exposure to oxidative air pollution (Miller et al., 2018).

The interaction between the genome and the epigenome will also change the evolution and control the development of different cancers. Epigenetic controlled immunotherapy has potential effects such as restoring anti-tumor signals (Gomez et al., 2020). Many sites in the DNAm region and epigenetic changes play an essential role in the early stage of liver cancer (Seal et al., 2020). The DNAm mechanism would regulate the immune cell response of cancer. Different DNAm patterns will also affect the characteristics of immune infiltration (Wu et al., 2020; Du et al., 2021). Epigenetics and gene expression influence each other (Mishra et al., 2015; Wang et al., 2020c). Studies have found that the CpG island of APLN has a more significant impact on lung tissue. The combined study of APLN and epigenetic DNA methylation sites are significantly related to survival. Research using epigenetics, gene transcriptome, and immune data, predicting prognostic markers of HCC has great potential.

The study still has limitations. The samples we studied are entirely retrospective. The data bias may affect survival and DNAm signatures. Although we determined the selected DNA methylation sites and APLN as predictive features, more experimental studies are still needed.

Conclusions

In summary, our research combines transcription, epigenetics, and immune data to predict HCC prognosis. The novel biomarkers immune-related APLN may be reliable prognostic markers for HCC.

Author Contribution: The authors confirm contribution to the paper as follows: Feifei Tian and Huan Hu wrote the manuscript and performed the statistically analyses; Di Wang and Huan Ding took part in data collection. Qingjia Chi and Huaping Liang prepared some figures and contributed to revise the manuscript; Wenli Zeng designed the work. All the authors read and approved the final manuscript.

Availability of Data and Materials: The datasets generated during the current study are available from the corresponding author on reasonable request.

Funding Statement: This work was supported by the Fund of Biosecurity Specialized Project of PLA (No. 19SWAQ18).

Conflicts of Interest: The authors declare that they have no conflicts of interest to report regarding the present study.

References

Amoozgar Z, Jain RK, Duda DG (2019). Role of apelin in glioblastoma vascularization and invasion after anti-VEGF therapy: What is the impact on the immune system? Cancer Research 79: 2104–2106. DOI 10.1158/0008-5472.CAN-19-0749. [Google Scholar] [CrossRef]

Ashokan A, Harisankar HS, Kameswaran M, Aradhyam GK (2021). Critical APJ receptor residues in extracellular domains that influence effector selectivity. FEBS Journal 288: 6543–6562. DOI 10.1111/febs.16048. [Google Scholar] [CrossRef]

Bai S, Wu Y, Yan Y, Shao S, Zhang J et al. (2020). Construct a circRNA/miRNA/mRNA regulatory network to explore potential pathogenesis and therapy options of clear cell renal cell carcinoma. Scientific Reports 10: 13659. DOI 10.1038/s41598-020-70484-2. [Google Scholar] [CrossRef]

Baylin SB, Jones PA (2011). A decade of exploring the cancer epigenome-biological and translational implications. Nature Reviews Cancer 11: 726–734. DOI 10.1038/nrc3130. [Google Scholar] [CrossRef]

Biswas SK, Mantovani A (2010). Macrophage plasticity and interaction with lymphocyte subsets: Cancer as a paradigm. Nature Immunology 11: 889–896. DOI 10.1038/ni.1937. [Google Scholar] [CrossRef]

Brown ZJ, Greten TF, Heinrich B (2019). Adjuvant treatment of hepatocellular carcinoma: Prospect of immunotherapy. Hepatology 70: 1437–1442. DOI 10.1002/hep.30633. [Google Scholar] [CrossRef]

Cabiati M, Gaggini M, de Simone P, Del RS (2021). Evaluation of Apelin/APJ system expression in hepatocellular carcinoma as a function of clinical severity. Clinical and Experimental Medicine 21: 269–275. DOI 10.1007/s10238-020-00672-x. [Google Scholar] [CrossRef]

Cao Z, Ao Y, Guo Y, Zhou S (2020). Comprehensive analysis of mRNA expression profiles in head and neck cancer by using robust rank aggregation and weighted gene coexpression network analysis. BioMed Research International 2020: 4908427. DOI 10.1155/2020/4908427. [Google Scholar] [CrossRef]

Chen H, Wong CC, Liu D, Go M, Wu B et al. (2019). APLN promotes hepatocellular carcinoma through activating PI3K/Akt pathway and is a druggable target. Theranostics 9: 5246–5260. DOI 10.7150/thno.34713. [Google Scholar] [CrossRef]

Chen Y, Xu T, Xie F, Wang L, Liang Z et al. (2021). Evaluating the biological functions of the prognostic genes identified by the Pathology Atlas in bladder cancer. Oncology Reports 45: 191–201. DOI 10.3892/or.2020.7853. [Google Scholar] [CrossRef]

Ding Z, Dong Z, Chen Z, Hong J, Yan L et al. (2021). Viral status and efficacy of immunotherapy in hepatocellular carcinoma: A systematic review with meta-analysis. Frontiers in Immunology 12: 733530. DOI 10.3389/fimmu.2021.733530. [Google Scholar] [CrossRef]

Du J, Ji H, Ma S, Jin J, Mi S et al. (2021). m6A regulator-mediated methylation modification patterns and characteristics of immunity and stemness in low-grade glioma. Briefings in Bioinformatics 22: bbab013. DOI 10.1093/bib/bbab013. [Google Scholar] [CrossRef]

Dzobo K, Senthebane DA, Ganz C, Thomford NE, Wonkam A et al. (2020). Advances in therapeutic targeting of cancer stem cells within the tumor microenvironment: An updated review. Cells 9: 1896. DOI 10.3390/cells9081896. [Google Scholar] [CrossRef]

Gao Q, Wang XY, Qiu SJ, Yamato I, Sho M et al. (2009). Overexpression of PD-L1 significantly associates with tumor aggressiveness and postoperative recurrence in human hepatocellular carcinoma. Clinical Cancer Research 15: 971–979. DOI 10.1158/1078-0432.CCR-08-1608. [Google Scholar] [CrossRef]

Garnelo M, Tan A, Her Z, Yeong J, Lim CJ et al. (2017). Interaction between tumour-infiltrating B cells and T cells controls the progression of hepatocellular carcinoma. Gut 66: 342–351. DOI 10.1136/gutjnl-2015-310814. [Google Scholar] [CrossRef]

Gomez S, Tabernacki T, Kobyra J, Roberts P, Chiappinelli KB (2020). Combining epigenetic and immune therapy to overcome cancer resistance. Seminars in Cancer Biology 65: 99–113. DOI 10.1016/j.semcancer.2019.12.019. [Google Scholar] [CrossRef]

Hall E, Volkov P, Dayeh T, Esguerra JL, Salo S et al. (2014). Sex differences in the genome-wide DNA methylation pattern and impact on gene expression, microRNA levels and insulin secretion in human pancreatic islets. Genome Biology 15: 522. DOI 10.1186/s13059-014-0522-z. [Google Scholar] [CrossRef]

Han Y, Liu D, Li L (2020). PD-1/PD-L1 pathway: Current researches in cancer. American Journal of Cancer Research 10: 727–742. [Google Scholar]

Hervouet F (1970). Treatment of blepharospasm by resection of the musculus orbicularis. Bulletin des societes d’ophtalmologie de France 70: 846. [Google Scholar]

Huo J, Wu L, Zang Y (2020). Development and validation of a novel immune-gene pairs prognostic model associated with CTNNB1 alteration in hepatocellular carcinoma. Medical Science Monitor 26: e925494. DOI 10.12659/MSM.925494. [Google Scholar] [CrossRef]

Ibrahim AE, Arends MJ, Silva AL, Wyllie AH, Greger L et al. (2011). Sequential DNA methylation changes are associated with DNMT3B overexpression in colorectal neoplastic progression. Gut 60: 499–508. DOI 10.1136/gut.2010.223602. [Google Scholar] [CrossRef]

Ioannidis JP, Tatsioni A, Karassa FB (2010). A vision for the European journal of clinical investigation: Note from the new editors. European Journal of Clinical Investigation 40: 1–3. DOI 10.1111/j.1365-2362.2009.02229.x. [Google Scholar] [CrossRef]

Jia Y, Zeng Z, Li Y, Li Z, Jin L et al. (2015a). Impaired function of CD4+ T follicular helper (Tfh) cells associated with hepatocellular carcinoma progression. PLoS One 10: e0117458. DOI 10.1371/journal.pone.0117458. [Google Scholar] [CrossRef]

Jia Y, Zeng Z, Li Y, Li Z, Jin L et al. (2015b). Impaired function of CD4+ T follicular helper (Tfh) cells associated with hepatocellular carcinoma progression. PLoS One 10: e0117458. DOI 10.1371/journal.pone.0117458. [Google Scholar] [CrossRef]

Kidoya H, Kunii N, Naito H, Muramatsu F, Okamoto Y et al. (2012). The apelin/APJ system induces maturation of the tumor vasculature and improves the efficiency of immune therapy. Oncogene 31: 3254–3264. DOI 10.1038/onc.2011.489. [Google Scholar] [CrossRef]

Kim JM, Chen DS (2016). Immune escape to PD-L1/PD-1 blockade: Seven steps to success (or failure). Annals of Oncology 27: 1492–1504. DOI 10.1093/annonc/mdw217. [Google Scholar] [CrossRef]

Korean Liver Cancer Association (2019). National cancer center korea practice guidelines for the management of hepatocellular carcinoma. Korean Journal of Radiology 20: 1042–1113. DOI 10.3348/kjr.2019.0140. [Google Scholar] [CrossRef]

Li XT, Zhang JT, Yan HH, Su J, Cheng ML et al. (2021). Gene co-expression modules integrated with immunoscore predicts survival of non-small cell lung cancer. Cancer Treatment and Research Communications 26: 100297. DOI 10.1016/j.ctarc.2020.100297. [Google Scholar] [CrossRef]

Lin Y, Xu J, Lan H (2019). Tumor-associated macrophages in tumor metastasis: Biological roles and clinical therapeutic applications. Journal of Hematology & Oncology 12: 76. DOI 10.1186/s13045-019-0760-3. [Google Scholar] [CrossRef]

Lin ZY, Chuang YH, Chuang WL (2012). Cancer-associated fibroblasts up-regulate CCL2, CCL26, IL6 and LOXL2 genes related to promotion of cancer progression in hepatocellular carcinoma cells. Biomedicine & Pharmacotherapy 66: 525–529. DOI 10.1016/j.biopha.2012.02.001. [Google Scholar] [CrossRef]

Liu J, He Z (2020). Small DNA methylation, big player in plant abiotic stress responses and memory. Frontiers in Plant Science 11: 595603. DOI 10.3389/fpls.2020.595603. [Google Scholar] [CrossRef]

Luo H, Zhao Q, Wei W, Zheng L, Yi S et al. (2020). Circulating tumor DNA methylation profiles enable early diagnosis, prognosis prediction, and screening for colorectal cancer. Science Translational Medicine 12: 556. DOI 10.1126/scitranslmed.aax7533. [Google Scholar] [CrossRef]

Malouff TD, Vallow LA, Seneviratne D, Mahajan A, Foote RL et al. (2020). Estimating the number of patients eligible for carbon ion radiotherapy in the United States. International Journal of Particle Therapy 7: 31–41. DOI 10.14338/IJPT-19-00079.1. [Google Scholar] [CrossRef]

Mastrella G, Hou M, Li M, Stoecklein VM, Zdouc N et al. (2019). Targeting APLN/APLNR improves antiangiogenic efficiency and blunts proinvasive side effects of VEGFA/VEGFR2 blockade in glioblastoma. Cancer Research 79: 2298–2313. DOI 10.1158/0008-5472.CAN-18-0881. [Google Scholar] [CrossRef]

Miller CN, Dye JA, Schladweiler MC, Richards JH, Ledbetter AD et al. (2018). Acute inhalation of ozone induces DNA methylation of apelin in lungs of Long-Evans rats. Inhalation Toxicology 30: 178–186. DOI 10.1080/08958378.2018.1483984. [Google Scholar] [CrossRef]

Mishra A, Kohli S, Dua S, Thinlas T, Mohammad G et al. (2015). Genetic differences and aberrant methylation in the apelin system predict the risk of high-altitude pulmonary edema. Proceedings of the National Academy of Sciences of the United States of America 112: 6134–6139. DOI 10.1073/pnas.1422759112. [Google Scholar] [CrossRef]

Mou MA, Keya NA, Islam M, Hossain MJ, Al HM et al. (2020). Validation of CSN1S1 transcriptional expression, promoter methylation, and prognostic power in breast cancer using independent datasets. Biochemistry and Biophysics Reports 24: 100867. DOI 10.1016/j.bbrep.2020.100867. [Google Scholar] [CrossRef]

Muto J, Shirabe K, Yoshizumi T, Ikegami T, Aishima S et al. (2014). The apelin-APJ system induces tumor arteriogenesis in hepatocellular carcinoma. Anticancer Research 34: 5313–5320. [Google Scholar]

Ondo K, Isono M, Nakano M, Hashiba S, Fukami T et al. (2021). The N(6)-methyladenosine modification posttranscriptionally regulates hepatic UGT2B7 expression. Biochemical Pharmacology 189: 114402. DOI 10.1016/j.bcp.2020.114402. [Google Scholar] [CrossRef]

Picault FX, Chaves-Almagro C, Projetti F, Prats H, Masri B et al. (2014). Tumour co-expression of apelin and its receptor is the basis of an autocrine loop involved in the growth of colon adenocarcinomas. European Journal of Cancer 50: 663–674. DOI 10.1016/j.ejca.2013.11.017. [Google Scholar] [CrossRef]

PoureauP, MetgesJ(2021). Fundamentals of Digestive Cancers Immunology, Especially Gastric and Hepatocellular CarcinomasFondamentaux de l’immunologie des Cancers Digestifs (Gastriques et Hépatocellulaires). Oncologie 23: 47–59. DOI 10.32604/Oncologie.2021.15525. [Google Scholar] [CrossRef]

Podgorska M, Pietraszek-Gremplewicz K, Nowak D (2018). Apelin effects migration and invasion abilities of colon cancer cells. Cells 7: 113. DOI 10.3390/cells7080113. [Google Scholar] [CrossRef]

Ran J, Li Y, Liu L, Zhu Y, Ni Y et al. (2021). Apelin enhances biological functions in lung cancer A549 cells by downregulating exosomal miR-15a-5p. Carcinogenesis 42: 243–253. DOI 10.1093/carcin/bgaa089. [Google Scholar] [CrossRef]

Rawat D, Shrivastava S, Naik RA, Chhonker SK, Mehrotra A et al. (2018). An overview of natural plant products in the treatment of hepatocellular carcinoma. Anti-cancer Agents in Medicinal Chemistry 18: 1838–1859. DOI 10.2174/1871520618666180604085612. [Google Scholar] [CrossRef]

Rohr-Udilova N, Klinglmuller F, Schulte-Hermann R, Stift J, Herac M et al. (2018). Deviations of the immune cell landscape between healthy liver and hepatocellular carcinoma. Science Reports 8: 6220. DOI 10.1038/s41598-018-24437-5. [Google Scholar] [CrossRef]

Seal DB, Das V, Goswami S, De RK (2020). Estimating gene expression from DNA methylation and copy number variation: A deep learning regression model for multi-omics integration. Genomics 112: 2833–2841. DOI 10.1016/j.ygeno.2020.03.021. [Google Scholar] [CrossRef]

Sharma A, Sharma KL, Bansal C, Kumar A (2020). Updates on “Cancer Genomics and Epigenomics”. World Journal of Clinical Oncology 11: 890–897. DOI 10.5306/wjco.v11.i11.890. [Google Scholar] [CrossRef]

Tian F, Zhou P, Kang W, Luo L, Fan X et al. (2015). The small-molecule inhibitor selectivity between IKKalpha and IKKbeta kinases in NF-kappaB signaling pathway. Journal of Receptor and Signal Transduction 35: 307–318. DOI 10.3109/10799893.2014.980950. [Google Scholar] [CrossRef]

Topalian SL, Taube JM, Anders RA, Pardoll DM (2016). Mechanism-driven biomarkers to guide immune checkpoint blockade in cancer therapy. Nature Reviews Cancer 16: 275–287. DOI 10.1038/nrc.2016.36. [Google Scholar] [CrossRef]

WangZ, HeY, SongY, WangY, ChenF(2021). MicroRNA in HCC: Biomarkers and Therapeutic Targets. Oncologie23: 177–184. DOI 10.32604/Oncologie.2021.014773. [Google Scholar] [CrossRef]

Wang R, Zhao H, Zhang Y, Zhu H, Su Q et al. (2020a). Identification of MicroRNA-92a-3p as an essential regulator of tubular epithelial cell pyroptosis by targeting Nrf1 via HO-1. Frontiers in Genetics 11: 616947. DOI 10.3389/fgene.2020.616947. [Google Scholar] [CrossRef]

Wang X, Liu X, Song Z, Shen X, Lu S et al. (2020b). Emerging roles of APLN and APELA in the physiology and pathology of the female reproductive system. PeerJ 8: e10245. DOI 10.7717/peerj.10245. [Google Scholar] [CrossRef]

Wang Y, Zhang X, Xiao X, Zhang FR, Yan X et al. (2020c). Accurately estimating the length distributions of genomic micro-satellites by tumor purity deconvolution. BMC Bioinformatics 21: 82. DOI 10.1186/s12859-020-3349-5. [Google Scholar] [CrossRef]

Wang YL, Zhang YY (2020). cg04448376, cg24387542, cg08548498, and cg14621323 as a novel signature to predict prognosis in kidney renal papillary cell carcinoma. BioMed Research International 2020: 4854390. DOI 10.1155/2020/4854390. [Google Scholar] [CrossRef]

Wildes TJ, Dyson KA, Francis C, Wummer B, Yang C et al. (2020). Immune escape after adoptive T-cell therapy for malignant gliomas. Clinical Cancer Research 26: 5689–5700. DOI 10.1158/1078-0432.CCR-20-1065. [Google Scholar] [CrossRef]

Wolf-Dennen K, Gordon N, Kleinerman ES (2020). Exosomal communication by metastatic osteosarcoma cells modulates alveolar macrophages to an M2 tumor-promoting phenotype and inhibits tumoricidal functions. OncoImmunology 9: 1747677. DOI 10.1080/2162402X.2020.1747677. [Google Scholar] [CrossRef]

Wu J, Mamidi T, Zhang L, Hicks C (2020). Unraveling the genomic-epigenomic interaction landscape in triple negative and non-triple negative breast cancer. Cancers 12: 1559. DOI 10.3390/cancers12061559. [Google Scholar] [CrossRef]

Xu RH, Wei W, Krawczyk M, Wang W, Luo H et al. (2017). Circulating tumour DNA methylation markers for diagnosis and prognosis of hepatocellular carcinoma. Nature Materials 16: 1155–1161. DOI 10.1038/nmat4997. [Google Scholar] [CrossRef]

Yang L, Li YL, Li XQ, Zhang Z (2019). Corrigendum to “High apelin level indicates a poor prognostic factor in muscle-invasive bladder cancer”. Disease Markers 2019: 5974238. DOI 10.1155/2019/5974238. [Google Scholar] [CrossRef]

Yang X, Yang J, Chu Y, Xue Y, Xuan D et al. (2014). T follicular helper cells and regulatory B cells dynamics in systemic lupus erythematosus. PLoS One 9: e88441. DOI 10.1371/journal.pone.0088441. [Google Scholar] [CrossRef]

Yin L, He N, Chen C, Zhang N, Lin Y et al. (2019). Identification of novel blood-based HCC-specific diagnostic biomarkers for human hepatocellular carcinoma. Artificial Cells, Nanomedicine, and Biotechnology 47: 1908–1916. DOI 10.1080/21691401.2019.1613421. [Google Scholar] [CrossRef]

Zuurbier L, Rahman A, Cordes M, Scheick J, Wong TJ et al. (2017). Apelin: A putative novel predictive biomarker for bevacizumab response in colorectal cancer. Oncotarget 8: 42949–42961. DOI 10.18632/oncotarget.17306. [Google Scholar] [CrossRef]

images This work is licensed under a Creative Commons Attribution 4.0 International License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.