- Primary research
- Open Access
DNA methylation biomarkers for hepatocellular carcinoma
Cancer Cell International volume 18, Article number: 140 (2018)
Aberrant methylation of DNA is a key driver of hepatocellular carcinoma (HCC). In this study, we sought to integrate four cohorts profile datasets to identify such abnormally methylated genes and pathways associated with HCC.
To this end, we downloaded microarray datasets examining gene expression (GSE84402, GSE46408) and gene methylation (GSE73003, GSE57956) from the GEO database. Abnormally methylated differentially expressed genes (DEGs) were sorted and pathways were analyzed. The String database was then used to perform enrichment and functional analysis of identified pathways and genes. Cytoscape software was used to create a protein–protein interaction network, and MCODE was used for module analysis. Finally, overall survival analysis of hub genes was performed by the OncoLnc online tool.
In total, we identified 19 hypomethylated highly expressed genes and 14 hypermethylated lowly expressed genes at the screening step, and finally found six mostly changed hub genes including MAD2L1, CDC20, CCNB1, CCND1, AR and ESR1. Pathway analysis showed that aberrantly methylated-DEGs mainly associated with the cell cycle process, p53 signaling, and MAPK signaling in HCC. After validation in TCGA database, the methylation and expression status of hub genes was significantly altered and same with our results. Patients with high expression of MAD2L1, CDC20 and CCNB1 and low expression of CCND1, AR, and ESR1 was associated with shorter overall survival.
Taken together, we have identified novel aberrantly methylated genes and pathways linked to HCC, potentially offering novel insights into the molecular mechanisms governing HCC progression and serving as novel biomarkers for precision diagnosis and disease treatment.
Hepatocellular carcinoma (HCC), an inflammation-driven disease, is the third deadliest cancer worldwide, and HCC prevalence is predicted to continue to rise in coming years, serving as a major economic burden [1, 2]. Most cases of HCC occur in developing countries, such as China, and the leading cause of HCC is infection with hepatitis B virus (HBV); in contrast, the main cause in developed countries, such as the USA, is infection with hepatitis C virus (HCV) [3, 4]. Other risk factors for developing HCC include exposure to aflatoxin, alcohol intake, smoking, and diabetes . The best curative treatments to date in early stage HCC patients involve surgical resection, tumor ablation, and potentially liver transplantation [6, 7]. However, the prognosis after curative therapy for HCC remains unsatisfactory because of a high postoperative recurrence rate. An improved understanding of the basic biology of HCC is needed to enhance prognostic predictions and to enhance therapeutic efficacy against this deadly disease.
The term epigenetics refers to heritable gene expression alterations no associated with DNA sequence changes . The DNA methylation is closely related to embryonic development , regulation of gene expression , X-chromosome inactivation , genomic imprinting , and genomic stability . Altered DNA methylation such as tumor suppressor gene hypermethylation or oncogene hypomethylation is thought to promote tumorigenesis. Genes including P15, P16, Ras association domain family 1 isoform A (RASSF1A), and Retinoblastoma 1 are inactivated in HCC due to promoter hypermethylation of these genes [14,15,16,17]. Given that methylation is potentially reversible, detection of such aberrant DNA methylation of tumor suppressors and oncogenes in HCC could be useful as a therapeutic target.
While altered methylation of many genes has been demonstrated to date in the context of HCC, a complete interaction network documenting the relationship between said genes remains to be produced. The comprehensive analysis of multiple datasets offers the power needed to properly identify and assess pertinent pathways and genes mediating the biological processes associated with HCC. To this end, we used datasets from microarrays examining gene expression (GSE84402, GSE46408) and gene methylation (GSE73003, GSE57956) to assess HCC gene and epigenetic signatures, allowing for identification of genes and pathways that were both abnormally methylated and differentially expressed. Using a protein–protein interaction network we were also able to identify key so-called “hub” genes central to these signaling events. Through this analysis, we believe it is possible to identify novel differentially methylated genes associated with HCC, offering key insights into the molecular mechanisms governing HCC development and progression.
Gene expression profiling datasets GSE84402 and GSE46408, and gene methylation profiling datasets GSE73003 and GSE57956, were downloaded from the gene expression omnibus (GEO, https://www.ncbi.nlm.nih.gov/geo/) [18, 19]. 14 HCC and 14 normal specimens were obtained in GSE84402, while 6 HCC and 6 normal samples were obtained in GSE46408. GSE57956 included a total of 59 primary HCC tumor samples and 59 adjacent normal tissue samples, while GSE73003 included 40 paired normal and HCC samples from 20 patients.
Data acquisition and processing
Raw gene expression profiling datasets GSE84402 and GSE46408 were downloaded from GEO public repositories. Data processing was performed using robust multi-array average (RMA) in GeneSpring GX 11.5 (Agilent Technologies Pty Ltd), including background adjustment, normalization and log transformation of the values. The threshold set for up and down-regulated genes was a | log FC | ≥ 2 and P ≤ 0.05. For gene methylation profiling dataset (GSE73003 and GSE57956), the GEO2R software was used to analyze the raw data and to identify differentially methylated genes (DMGs). GEO2R is an interactive web tool which allows users to compare different groups of samples in a GEO series to screen genes that are differentially expressed in experimental conditions. The adjusted P-values (adj. P) and Benjamini and Hochberg false discovery rate were applied to provide a balance between discovery of statistically significant genes and limitations of false-positives. The DMGs cutoff criteria were P < 0.05 and |t| > 2. Additionally, one-sided tests were used to categorize the upregulated or downregulated DMGs. Finally, a Venn diagram was used to identify hypomethylated highly-expressed genes and hypermethylated lowly-expressed genes.
Functional and pathway enrichment analysis
Gene ontology (GO) analyses and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were conducted on identified genes with altered methylation/expression using the Search Tool for the Retrieval of Interacting Genes (STRING, https://string-db.org/). This allowed for protein–protein interaction (PPI) network generation and functional annotation of genes of interest. P < 0.05 was the threshold for statistical significance.
Generation and analysis of a protein–protein interaction (PPI) network
In order to interpret the drivers of carcinogenesis in a meaningful manner, a functional PPI analysis was necessary. The STRING database allowed us to generate PPI networks for both hypomethylated/and highly-expressed genes as well as for hypermethylated and lowly-expressed genes. A 0.4 interaction score served as the cutoff prior to visualization. Modules were then screened using Molecular Complex Detection (MCODE) in the Cytoscape software package, with MCODE score > 4 and number of nodes > 5. The top 3 Hub genes were selected by CytoHubba app in Cytoscape software. To determine the expression pattern of six hub genes in HCC, we used the datasets in the Oncomine (https://www.oncomine.org) and UALCAN (http://ualcan.path.uab.edu/) database. Oncomine and UALCAN are an online database consisting of previously published and open-access microarray data. The analysis enables multiple comparisons of gene expression between different studies; the significance of the gene expression across the available studies was also taken into account [20, 21]. The results were filtered by selecting hepatocellular carcinoma vs. normal tissue.
Hub gene validation
The Cancer Genome Atlas (TCGA) database has generated comprehensive, multi-dimensional maps of the key genomic changes in various types of cancers. MEXPRESS (http://mexpress.be/) is a data visualization tool designed for the easy visualization of TCGA expression, DNA methylation and clinical data, as well as the relationships between them. To confirm our results, we used the MEXPRESS to validate hypermethylation/low-expression hub genes and hypomethylation/high-expression hub genes in TCGA database. The probability of survival and significance were calculated using the OncoLnc database. OncoLnc (http://www.oncolnc.org/) is an online tool for interactively exploring survival correlations, containing survival data for 8647 patients from 21 cancer studies.
Identification of abnormally methylated and differentially expressed genes in HCC
A study design flowchart is shown in Fig. 1. In gene expression microarrays, 127 genes were up-regulated in both datasets (207 in GSE84402, 851 in GSE46408) and 169 were down-regulated in both datasets (307 in GSE84402, 503 in GSE46408). For the gene methylation microarrays, 2139 overlapping hypermethylated genes (2468 in GSE73003, 3157 in GSE57956) and 4091 overlapping hypomethylated genes (4527 in GSE73003, 5334 in GSE57956) were identified. By comparing these two sets of genes, a total of 19 hypomethylated, highly-expressed genes and 14 hypermethylated, low-expression genes were identified (Fig. 2). A heat map of these genes in GSE84402 is shown in Fig. 3.
Gene ontology and pathway functional enrichment analysis
The GO annotation and pathway enrichment analyses of identified aberrantly methylated and differentially expressed genes were implemented using the online tool STRING. Genes that were hypomethylated and highly expressed were enriched for cell cycle regulation, while the hypermethylated low-expression genes were primarily linked to cell proliferation, gene expression, and signal transduction. Cell component enrichment analysis indicated that hypomethylated highly-expressed genes were correlated with cytoskeletal part, whereas hypermethylation/low-expression genes were predominant at intracellular organelle lumen. As for molecular function, hypomethylation/high-expression genes were enriched mainly in roundabout binding, histone kinase activity, and protein binding, while hypermethylation/low-expression genes were mostly enriched in cyclin binding, receptor binding, and signal transducer activity (Table 1). The pathway analysis revealed that hypomethylated highly-expressed genes were linked to oocyte meiosis, cell cycle, and ubiquitin mediated proteolysis, while hypermethylation/low-expression genes significantly enriched in pathways in cancer, p53 signaling pathway, MAPK signaling pathway, and proteoglycans in cancer (Table 2).
Construction and analysis of PPI networks
The STRING database was used for PPI network construction, with MCODE used for module analysis. Hub genes were selected using the cytoHubba Cytoscape package. The PPI network for genes that were hypomethylated and highly expressed is shown in Fig. 4a, with corresponding modules shown in Fig. 5a. The most significantly enriched functional modules were those linked to progesterone-mediated oocyte maturation, the cell cycle, HTLV-I infection, oocyte meiosis, and ubiquitin mediated proteolysis (Fig. 5b). Top three hub genes were MAD2L1, CDC20 and CCNB1 (Fig. 4b). The PPI network for genes that were hypermethylated and expressed at low levels is shown in Fig. 4c, with corresponding modules shown in Fig. 5c, e. Significant vital modules showed functions including p53 signaling pathway, cell cycle, viral carcinogenesis, transcriptional misregulation in cancer, neurotrophin signaling pathway, MAPK signaling pathway, apoptosis and inflammatory mediator regulation of TRP channels (Fig. 5d, f). Top three hub genes were CCND1, AR and ESR1 (Fig. 4d). Furthermore, we use Oncomine and UALCAN database to confirm the expression of hub genes in HCC (Fig. 6, Additional file 1: Figure S1). The data are in agreement with our results.
Hub gene validation
Hub genes for the Hypermethylated/low-expression and Hypomethylation/high-expression datasets were then validated using the TCGA database in order to affirm the validity of these findings. The MEXPRESS web tool allowed for easy and rapid assessment of the relationships between methylation and gene expression within TCGA at a single-gene level. By ordering samples based on gene expression levels, MEXPRESS makes it clear that hub gene expression was negatively correlated with methylation, and a Pearson correlation confirmed this finding. For the hypermethylation/low-expression hub genes, gene expression tended to be higher in normal samples than in tumor samples (Additional file 2: Figure S2). However, tumor samples tended to have higher expression than normal samples for hypomethylation/high-expression hub genes (Additional file 3: Figure S3). The outcomes are summarized in Table 3. Gene expression and methylation were altered and consistent with our findings, indicating that these results are reliable and reproducible.
The relationship between hub genes and the survival in hepatocellular carcinoma
The prognostic value of six hub genes was assessed by OncoLnc database. The threshold was adjusted to cox P value < 0.05. Patients with high expression of MAD2L1, CDC20 and CCNB1 and low expression of CCND1, AR, and ESR1 was associated with shorter overall survival (Fig. 7). The survival data showed that MAD2L1, CDC20 and CCNB1 played an oncogenic role, while CCND1, AR, and ESR1 genes were associated with better survival in HCC.
Exploring the mechanisms underlying HCC development and progression not only has prognostic implications, but may also be helpful in monitoring treatment response, surveillance for tumor recurrence, and guidance of clinical decisions. Modern advances in sequencing technologies and microarray development have provided ample high-throughput opportunities to study disease-related biology, allowing for simultaneous assessment of gene methylation and expression for thousands of genes in the human genome. In the present study, we identified 19 hypomethylated, highly-expressed genes and 14 hypermethylated, low-expression genes, using bioinformatic analysis. Functional enrichment of these genes revealed that aberrant methylation indeed affects certain pathways and hub genes. These results can provide novel insight into the explanation of HCC pathogenesis.
The GO enrichment analysis revealed that the primary biological processes of the hypomethylated/highly-expressed genes were the regulation of cell cycle processes, chromosome segregation, and mitotic nuclear division while the hypermethylated/low-expression genes were involved mainly in controlling cell proliferation, gene expression, and the mitotic G1/S transition. This is expected given that the chromosome segregation process occurs during mitosis, which is a part of the cell cycle. The G1/S cell cycle transition is tightly controlled. Deregulation of this key checkpoint can allow cells to undergo transformation, thereby permitting tumorigenesis. This finding is consistent with the fundamental role played by cell cycle regulators in cell proliferation, invasion, and metastasis in HCC. KEGG pathway analysis of hypomethylation/high-expression genes revealed that they were linked to the cell cycle, oocyte meiosis, and ubiquitin-mediated proteolysis. The cell cycle and oocyte meiosis are vital for cell proliferation in tumor cells, and the ubiquitin proteasome pathway functions to regulate cell cycle control and the DNA damage response in tumor genesis . The KEGG pathway analysis of hypermethylation/low-expression genes suggested that methylation may be linked to HCC development through the p53 and MAPK signaling pathways. p53 is a tumor suppressor to conserve genome stability by preventing mutations caused by cellular stress or DNA damage. Together, these results suggest that hypermethylation and hypomethylation are key mediators of cancer development and progression.
The PPI network of hypomethylated/highly-expressed genes provides insight into the functional associations between them, and from this, the top three hub genes were selected: MAD2L1, CDC20, and CCNB1. Mitotic Arrest Deficient 2 Like 1 (MAD2L1) and Cell Division Cycle 20 (CDC20) are two key mitotic checkpoint genes. Both MAD2L1 and CDC20 were more highly expressed in higher grade tumors than in low-grade tumors. High MAD2L1 or CDC20 levels may allow for the development of multi-organ aggressive tumors, including those affecting the breasts, lungs, liver, and stomach [23,24,25]. Collectively, these findings suggest that MAD2L1 and CDC20 may be key regulators of tumor severity, ultimately dictating patient survival. Cyclin B1 (CCNB1), complexing with CDC2, is a G2/M-phase checkpoint regulator that is vital for regulation of proliferation and DNA synthesis. CCNB1 overexpression has been found to occur in HCC  and many other cancer, often being linked to progression, recurrence, and to poor prognoses . Therefore, MAD2L1, CDC20, and CCNB1 may all be abnormally methylated genes that modulate the cell cycle and proliferation in HCC. With regard to the hypermethylated/low-expression genes, the most prominent hub genes were CCND1, AR, and ESR1. Cyclin D1 (CCND1) is a proto-oncogene regulating G1 to S phase progression; it participates in the Wnt/β-catenin pathway [28, 29]. Androgen receptor (AR) is a steroid hormone receptor superfamily member that is involved in human hepatocarcinogenesis . It alters the AR-dependent transcriptome and stimulates oncogenic growth. Estrogen receptor 1 (ESR1) functions as a transcription factor, regulating cell cycle, cell proliferation, apoptosis, and inflammation-associated gene expression . Research has shown that estrogen-depleted postmenopausal women undergo more rapid progression from HCV-infection to HCC development . Aberrant expression of ESR subtypes may contribute to the progression of HCC. These three genes are related to prognosis, tumorigenesis, and metastasis of HCC. Furthermore, survival analysis of hub genes revealed that MAD2L1, CDC20 and CCNB1 play an oncogenic role, while CCND1, AR, and ESR1 genes were associated with favorable patient survival in HCC.
The core PPI network module for hypomethylated/highly-expressed genes was linked to the cell cycle, oocyte meiosis, and ubiquitin-mediated proteolysis, indicating that these pathways are key targets of hypermethylation. The top two modules of the hypermethylated/low-expression gene PPI network were those linked to the p53 signaling pathways, viral carcinogenesis and neurotrophin signaling pathways. p53 signaling conserves the stability of the genome. The leading cause of HCC is infection with HBV or HCV. It is reasonable that viral carcinogenesis is involved in the development of HCC. Neurotrophins and neurotrophin receptors are found on tumor and stromal cells, and are linked to many kinds of tumor development. Brain-derived neurotrophic factor and nerve growth factor are both neurotrophins linked to tumor development, promoting proliferation, angiogenesis, and invasion.
While previously groups have assessed arrays cataloging gene expression or methylation, the two have not been examined simultaneously. Furthermore, single studies generally lack the power needed to identify critical regulatory genes and signaling pathways. Our research used a bioinformatics workflow to jointly analyze extant gene expression and gene methylation profiling microarrays, allowing for more powerful and precise insights into these screening results. However, we only validated candidate abnormally methylated genes that were differentially expressed using the TCGA database. Further experiments will be necessary in order to confirm that these genes and pathways are linked to HCC.
In summary, this study provides a comprehensive bioinformatics analysis of aberrantly methylated DEGs that may be involved in the progression and development of HCC. In addition, six mostly changed hub genes were identified, including MAD2L1, CDC20, CCNB1, CCND1, AR, and ESR1. These novel findings may contribute to the unraveling of the pathogenesis of HCC, and these candidate genes may be optimal abnormal methylation-based biomarkers that can be used to accurately diagnose and treat HCC.
differentially expressed genes
differentially methylated genes
Ras association domain family 1 isoform A
gene expression omnibus
Kyoto Encyclopedia of Genes and Genomes
Search Tool for the Retrieval of Interacting Genes
Molecular Complex Detection
The Cancer Genome Atlas
Mitotic Arrest Deficient 2 Like 1
Cell Division Cycle 20
estrogen receptor 1
nerve growth factor
brain-derived neurotrophic factor
Addeo R, Caraglia M, Del Prete S. Highlights of regional meeting of Italian Southern Oncological Group (GOIM): focus on hepatocellular carcinoma: biological and clinical background, therapeutic guide-lines and perspectives. 7 November 2008, Naples, Italy. Expert Opin Investig Drugs. 2009;18(3):373–8.
Ringelhan M, Pfister D, O’Connor T, Pikarsky E, Heikenwalder M. The immunology of hepatocellular carcinoma. Nat Immunol. 2018;19(3):222–32.
Jaquet A, Tchounga B, Tanon A, Bagny A, Ekouevi DK, Traore HA, Sasco AJ, Maiga M, Dabis F. Etiology of hepatocellular carcinoma in West Africa, a case-control study. Int J Cancer. 2018;143(4):869–77.
Ligat G, Schuster C, Baumert TF. HBV core variants, liver fibrosis and hepatocellular carcinoma. Hepatology. 2018. https://doi.org/10.1002/hep.30231.
Jemal A, Bray F, Center MM, Ferlay J, Ward E, Forman D. Global cancer statistics. CA Cancer J Clin. 2011;61(2):69–90.
Venook AP, Papandreou C, Furuse J, de Guevara LL. The incidence and epidemiology of hepatocellular carcinoma: a global and regional perspective. Oncologist. 2010;15(Suppl 4):5–13.
Younossi Z, Anstee QM, Marietti M, Hardy T, Henry L, Eslam M, George J, Bugianesi E. Global burden of NAFLD and NASH: trends, predictions, risk factors and prevention. Nat Rev Gastroenterol Hepatol. 2018;15(1):11–20.
Wong CC, Qian Y, Yu J. Interplay between epigenetics and metabolism in oncogenesis: mechanisms and therapeutic approaches. Oncogene. 2017;36(24):3359–74.
Yin LJ, Zhang Y, Lv PP, He WH, Wu YT, Liu AX, Ding GL, Dong MY, Qu F, Xu CM, et al. Insufficient maintenance DNA methylation is associated with abnormal embryonic development. BMC Med. 2012;10:26.
Shi B, Thomas AJ, Benninghoff AD, Sessions BR, Meng Q, Parasar P, Rutigliano HM, White KL, Davies CJ. Genetic and epigenetic regulation of major histocompatibility complex class I gene expression in bovine trophoblast cells. Am J Reprod Immunol. 2018;79(1):e12779.
Pasque V, Karnik R, Chronis C, Petrella P, Langerman J, Bonora G, Song J, Vanheer L, Sadhu Dimashkie A, Meissner A, et al. X chromosome dosage influences DNA methylation dynamics during reprogramming to mouse iPSCs. Stem Cell Rep. 2018;10(5):1537–50.
Lee CC, Peng SH, Shen L, Lee CF, Du TH, Kang ML, Xu GL, Upadhyay AK, Cheng X, Yan YT, et al. The role of N-alpha-acetyltransferase 10 protein in DNA methylation and genomic imprinting. Mol Cell. 2017;68(1):89–103.
Hernandez-Saavedra D, Strakovsky RS, Ostrosky-Wegman P, Pan YX. Epigenetic regulation of centromere chromatin stability by dietary and environmental factors. Adv Nutr. 2017;8(6):889–904.
Xu T, Cui T, Peng L, Kong S, Zou J, Tian X. The anti-hepatocellular carcinoma activity of Mel-P15 is mediated by natural killer cells. Oncol Lett. 2017;14(6):6901–6.
Sun G, Zhang C, Feng M, Liu W, Xie H, Qin Q, Zhao E, Wan L. Methylation analysis of p16, SLIT2, SCARA5, and Runx3 genes in hepatocellular carcinoma. Medicine. 2017;96(41):e8279.
Zhang X, Guo C, Wu X, Li AX, Liu L, Tsark W, Dammann R, Shen H, Vonderfecht SL, Pfeifer GP. Analysis of liver tumor-prone mouse models of the hippo kinase scaffold proteins RASSF1A and SAV1. Cancer Res. 2016;76(9):2824–35.
Ishak CA, Marshall AE, Passos DT, White CR, Kim SJ, Cecchini MJ, Ferwati S, MacDonald WA, Howlett CJ, Welch ID, et al. An RB-EZH2 complex mediates silencing of repetitive DNA sequences. Mol Cell. 2016;64(6):1074–87.
Wang H, Huo X, Yang XR, He J, Cheng L, Wang N, Deng X, Jin H, Wang C, Zhao F, et al. STAT3-mediated upregulation of lncRNA HOXD-AS1 as a ceRNA facilitates liver cancer metastasis by regulating SOX4. Mol Cancer. 2017;16(1):136.
Mah WC, Thurnherr T, Chow PK, Chung AY, Ooi LL, Toh HC, Teh BT, Saunthararajah Y, Lee CG. Methylation profiles reveal distinct subgroup of hepatocellular carcinoma patients with poor prognosis. PLoS ONE. 2014;9(8):e104158.
Roessler S, Jia HL, Budhu A, Forgues M, Ye QH, Lee JS, Thorgeirsson SS, Sun Z, Tang ZY, Qin LX, et al. A unique metastasis gene signature enables prediction of tumor relapse in early-stage hepatocellular carcinoma patients. Cancer Res. 2010;70(24):10202–12.
Chen X, Cheung ST, So S, Fan ST, Barry C, Higgins J, Lai KM, Ji J, Dudoit S, Ng IO, et al. Gene expression patterns in human liver cancers. Mol Biol Cell. 2002;13(6):1929–39.
Tu Y, Chen C, Pan J, Xu J, Zhou ZG, Wang CY. The Ubiquitin Proteasome Pathway (UPP) in the regulation of cell cycle control and DNA damage repair and its implication in tumorigenesis. Int J Clin Exp Pathol. 2012;5(8):726–38.
Liu X, Zhou W, Zhang X, Ding Y, Du Q, Hu R. 1-L-MT, an IDO inhibitor, prevented colitis-associated cancer by inducing CDC20 inhibition-mediated mitotic death of colon cancer cells. Int J Cancer. 2018;143(6):1516–29.
Wang Z, Katsaros D, Shen Y, Fu Y, Canuto EM, Benedetto C, Lu L, Chu WM, Risch HA, Yu H. Biological and clinical significance of MAD2L1 and BUB1, genes frequently appearing in expression signatures for breast cancer prognosis. PLoS ONE. 2015;10(8):e0136246.
Schvartzman JM, Sotillo R, Benezra R. Mitotic chromosomal instability and cancer: mouse modelling of the human disease. Nat Rev Cancer. 2010;10(2):102–15.
Covini G, Chan EK, Nishioka M, Morshed SA, Reed SI, Tan EM. Immune response to cyclin B1 in hepatocellular carcinoma. Hepatology. 1997;25(1):75–80.
Egloff AM, Vella LA, Finn OJ. Cyclin B1 and other cyclins as tumor antigens in immunosurveillance and immunotherapy of cancer. Cancer Res. 2006;66(1):6–9.
Karami F, Mohammadi-Yeganeh S, Abedi N, Koochaki A, Kia V, Paryan M. Bioinformatics prediction and in vitro analysis revealed that miR-17 targets cyclin D1 mRNA in triple negative breast cancer cells. Chem Biol Drug Des. 2016;87(3):317–20.
Karimkhanloo H, Mohammadi-Yeganeh S, Ahsani Z, Paryan M. Bioinformatics prediction and experimental validation of microRNA-20a targeting Cyclin D1 in hepatocellular carcinoma. Tumour Biol. 2017;39(4):1010428317698361.
Zhang H, Li XX, Yang Y, Zhang Y, Wang HY, Zheng XFS. Significance and mechanism of androgen receptor overexpression and androgen receptor/mechanistic target of rapamycin cross-talk in hepatocellular carcinoma. Hepatology. 2018;67(6):2271–86.
Heldring N, Pike A, Andersson S, Matthews J, Cheng G, Hartman J, Tujague M, Strom A, Treuter E, Warner M, et al. Estrogen receptors: how do they signal and what are their targets. Physiol Rev. 2007;87(3):905–31.
Iyer JK, Kalra M, Kaul A, Payton ME, Kaul R. Estrogen receptor expression in chronic hepatitis C and hepatocellular carcinoma pathogenesis. World J Gastroenterol. 2017;23(37):6802–16.
Conceived and designed the study strategy: XC, CDW; Acquisition of data: statistical analysis and interpretation of data GRF, YQT; Drafting or revision of the manuscript: YQT, GRF; Reference collection and data management: CC, HYS; Wrote the manuscript: YQT; Prepared the tables and figures: GRF; Study supervision: XC, CDW; All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets generated and/or analysed during the current study are available in the GEO repository, https://www.ncbi.nlm.nih.gov/geo/.
Consent for publication
Ethics approval and consent to participate
This work was supported by the National Natural Science Foundation of China (No. 81500796) and Research fund of Wuhan Union hospital (No. 02.03.2017-321).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.