- Primary research
- Open Access
Downexpression of HSD17B6 correlates with clinical prognosis and tumor immune infiltrates in hepatocellular carcinoma
Cancer Cell International volume 20, Article number: 210 (2020)
- The Correction to this article has been published in Cancer Cell International 2020 20:479
Hydroxysteroid 17-Beta Dehydrogenase 6 (HSD17B6), a key protein involved in synthetizing dihydrotestosterone, is abundant in the liver. Previous studies have suggested a role for dihydrotestosterone in modulating progress of various malignancies, and HSD17B6 dysfunction was associated with lung cancer and prostate cancer. However, little is known about the detailed role of HSD17B6 in hepatocellular carcinoma (HCC).
Clinical implication and survival data related to HSD17B6 expression in patients with HCC were obtained through TCGA, ICGC, ONCOMINE, GEO and HPA databases. Survival analysis plots were drawn with Kaplan–Meier Plotter. The ChIP-seq data were obtained from Cistrome DB. Protein–Protein Interaction and gene functional enrichment analyses were performed in STRING database. The correlations between HSD17B6 and tumor immune infiltrates was investigated via TIMER and xCell. The proliferation, migration and invasion of liver cancer cells transfected with HSD17B6 were evaluated by the CCK8 assay, wound healing test and transwell assay respectively. Expression of HSD17B6, TGFB1 and PD-L1 were assessed by quantitative RT-PCR.
HSD17B6 expression was lower in HCC compared to normal liver and correlated with tumor stage and grade. Lower expression of HSD17B6 was associated with worse OS, PFS, RFS and DSS in HCC patients. HNF4A bound to enhancer and promoter regions of HSD17B6 gene, activating its transcription, and DNA methylation of HSD17B6 promoter negatively controlled the expression. HSD17B6 and its interaction partners were involved in androgen metabolism and biosynthesis in liver. HSD17B6 inhibited tumor cell proliferation, migration and invasion in liver cancer cells and low expression of HSD17B6 correlated with high immune cells infiltration, relative reduction of immune responses and multiple immune checkpoint genes expression in HCC, probably by regulating the expression of TGFB1.
This study indicate that HSD17B6 could be a new biomarker for the prognosis of HCC and an important negative regulator of immune responses in HCC.
Liver cancer is the sixth most common cancer and the fourth leading cause of cancer-related deaths worldwide, with 700,000 annual deaths in recent years [1, 2]. It is also the second most lethal tumor with only 18% of 5-year survival . Hepatocellular carcinoma (HCC) accounts for 75–85% of primary liver cancers [2, 4], and various pathogenic factors can lead to its occurrence and development, such as chronic hepatitis B or C virus infections and alcohol abuse. Patients with HCC have sustained hepatic inflammation, fibrosis, which cause a series of genetic and epigenetic events. Despite tremendous efforts in the past decades, the incidence of HCC is increasing rapidly, and patients with advanced stage HCC have very poor outcomes . Advances in HCC treatment will depend on a better understanding of the precise molecular mechanism involved in HCC. Although HCC is a clinically and biologically heterogeneous malignancy, most HCCs share common features regarding genetic and epigenetic alterations. Identification of these common molecular alterations in HCC might provide a rational strategy for the development of effective molecular targeted therapy in HCC.
HSD17B6 is a gene located on chromosome 12q13.3 and encodes a protein with both oxidoreductase and epimerase activities, involved in steroid metabolism. The oxidoreductase activity can convert 3 alpha-adiol to dihydrotestosterone (DHT), while the epimerase activity can convert androsterone to epi-androsterone [5, 6]. DHT was closely related to the development of many tumors, including prostate cancer, breast cancer, endometrial cancer, lung cancer and colorectal cancer . Polymorphisms in the HSD17B6 gene are associated with polycystic ovary syndrome (PCOS) and key clinical phenotypes of the disorder [8, 9]. Expression of HSD17B6 in prostate cancer samples with bone metastasis were significantly lower than that in non-metastatic counterparts , indicating the dysfunction of HSD17B6 in cancer metastases. It has also been identified to be associated with non-small-cell lung cancer  and its SNP was significantly associated with liver fibrosis risk . However, little is known about the role of HSD17B6 in the pathogenesis of HCC.
In the present study, we evaluated the expression patterns, potential functions, and prognostic values of HSD17B6 in HCC in detail, and found that HSD17B6 was frequently down-regulated and predicted poor prognosis in HCC. Moreover, our study indicated that transcription of HSD17B6 was positively regulated by hepatocyte nuclear factor 4 alpha (HNF4A), a liver-enriched transcription factor. Besides, DNA hypermethylation of the HSD17B6 promoter also contributed a lot to the low expression of HSD17B6 in HCC. Moreover, downregulation of HSD17B6 affected the cell proliferation, migration, invasion, androgen metabolism and biosynthesis in HCC, and was associated with infiltration of tumor immune cells and expression of multiple immune checkpoint genes, indicating an important role in HCC development, progression and immune therapy.
Oncomine gene expression array datasets (https://www.oncomine.org, an online cancer microarray database) were used to analyze the transcription levels of gene expression in normal liver and liver disease tissues. The mRNA expressions of HSD17B6 in clinical cancer specimens were compared with those in normal controls, using a Student’s t test to generate a p value.
The Kaplan–Meier plotter
The prognostic value of HSD17B6 mRNA expression in TCGA-LIHC was evaluated using an online database, Kaplan–Meier Plotter (http://www.kmplot.com) , which contained gene expression data and survival information of liver cancer patients (http://kmplot.com/analysis/index.php?p=service&cancer=liver_rnaseq). To analyze the OS (Overall survival), PFS (Progression-free survival), RFS (Relapse-free survival) and DSS (Disease specific survival) of patients with HCC, patient samples were split into two groups by auto select best cutoff (high versus low expression) and assessed by a Kaplan–Meier survival plot, with the hazard ratio (HR) with 95% confidence intervals (CIs) and log rank p value.
The cancer genome atlas and GTEx data
The cancer genome atlas (TCGA) had both sequencing and pathological data on 30 different cancers . The liver hepatocellular carcinoma (The cancer genome atlas, TCGA-LIHC) dataset, including data from 438 cases with pathology reports, was selected for further analyses of HSD17B6 using Xena Functional Genomics Explorer, which is a visual exploration resource for both public and private omics data, supported through the web-based Xena Browser and multiple turn-key Xena Hubs. The Genotype Tissue Expression (GTEx) Project is a data resource and tissue bank established by the National Institutes of Health Common Fund, and total of 53 human normal tissues from nearly 1000 individuals have been studied by genomic and RNA sequencing. The GTEx Project was supported by the Common Fund of the Office of the Director of the National Institutes of Health, and by NCI, NHGRI, NHLBI, NIDA, NIMH, and NINDS. The data used for the analyses described in this manuscript were obtained from GTEx Portal V8 Release on 10/10/2019.
The basic expression of protein in liver and cancer tissues
We used the online tool “The human protein atlas” (HPA, https://www.proteinatlas.org) to analyze the basic protein expression level of HSD17B6 in human liver and liver cancer tissues [15,16,17]. The HPA is Swedish-based program initiated in 2003 with the aim to map all human proteins in cells, tissues, and organs using the integration of various omics technologies including antibody-based imaging. We worked with version 19 in this study, which was launched in September 2019 [15, 17]. With the purpose of analyzing the protein expression differences of HSD17B6 in normal liver and liver cancer tissues in live cancer patients, immunohistochemistry pictures were downloaded from the Tissue Atlas and Pathology Atlas in HPA.
Protein–protein interaction (PPI) network construction and gene function and pathway enrichment analysis
We used the search tool for the retrieval of interacting genes (STRING) database to analyze functional interactions between proteins . We retrieved interactions with confidence scores greater or equal to 0.9 and screened the resulting PPI network for analysis. And, we then performed a functional enrichment analysis of these genes using STRING database.
Computational deconvolution of infiltrating immune cells
In order to evaluate the correlation of the infiltrating immune cell subsets in HCC samples with the expression of HSD17B6, we performed the deconvolution analysis through two web-based tools, TIMER (http://timer.cistrome.org/) and xCell (http://xcell.ucsf.edu/), to infer the presence in TCGA LIHC and ICGC LIRI-JP datasets. TIMER could estimate the abundances of six immune cell types, including B cells, CD4+ T cells, CD8+ T cells, neutrophils, macrophages, and dendritic cells via “Estimation” modules. xCell could estimate the abundance scores of 64 cell types including 34 kinds of immune cells. The correlations between HSD17B6 expression and abundance scores of immune cells evaluated by Spearman’s correlation.
The HCCDB database (http://lifeome.net/database/hccdb/home.html) archived 15 public HCC gene expression datasets containing totally 3917 samples (13 GEO microarray datasets and two RNA-Seq datasets, TCGA-LIHC and ICGC LIRI-JP), serving as a one-stop online resource for exploring HCC gene expression with user-friendly interfaces . HCCDB was used to analyze the expression levels of HSD17B6 and over survival in different datasets.
Gene set enrichment analysis (GSEA)
Correlation coefficients between each gene and HSD17B6 expression generated using Spearman correlation analysis were enrolled to create the pre-ranked gene list. Then, a pre-ranked GSEA analysis was performed on GO-term database (c5.bp.v7.0.symbols.gmt) of Molecular Signatures Database (MSigDB) with the pre-ranked list. Gene sets with FDR < 0.25 and specifically enriched at the beginning and end of the ranked list were considered to be enrichment significance .
Cell culture and lentivirus transfection
The HCC cell line HepG2, which were obtained from ATCC (catalog number HB-8065), were cultured in DMEM medium plus 10% fetal bovine serum (cat. no. 10099‑141; Thermo Fisher Scientific, Inc.) at 37 °C with 5% CO2. Antibiotic or antifungal agents were not used, to avoid the potential effects of these agents on gene expression and cytotoxicity assays. The recombinant lentivirus expressing HSD17B6 and empty lentivector expressing GFP alone were purchased from HanBio (Shanghai, China). To establish stable transfecting cell lines, HepG2 cells were transduced with lentivirus in the presence of 5 μg/ml polybrene. Puromycin was added into medium at the concentration of 2 mg/ml for stable cell line selection. After antibiotic selection for 2 weeks, stable transfected HepG2 cell lines were established.
Cell proliferation assay
Cells in the logarithmic growth phase were seeded in 96-well plates at a cell density of 2 × 103/well (in triplicate) to allow adhesion. At 0, 24, 48, 72 and 96 h, cells were incubated with 10 µl CCK‑8 at 37 °C for an extra 2 h. The optical density was then measured with a microplate reader (Tecan Group Ltd.) at 450 nm. The cells were then cultured with fresh medium until the next round of measurements. The mean and standard deviation of the triplet measurements were calculated and plotted.
Confluent cells were serum‑starved in DMEM medium for 10–12 h and scratched using the tip of a 10-µl pipette. After being washed twice with PBS to remove non-adherent cells, the plates were added 500 µl DMEM medium plus 10% FBS. The wound area was photographed at 0, 24 and 48 h under an Olympus IX73 inverted microscope. A cell-free region was drawn and measured by CellSens Standard software (Olympus).
Transwell invasion assay
A BioCoat™ Matrigel invasion chamber (cat. no. 40480; BD Biosciences) was used according to the manufacturer’s protocol. Briefly, 4 × 104 cells were trypsinized, washed, suspended in 200 µl DMEM medium, and seeded in the upper portion of the invasion chamber. The lower portion of the chamber contained 500 µl of DMEM medium plus 10% FBS, which served as a chemo-attractant. After 48 h, the non-invasive cells were removed from the upper surface of the membrane with a cotton swab. The invasive cells on the lower surface of the membrane were stained with 0.1% crystal violet for 30 min at room temperature, and counted in four separate areas with an inverted microscope.
Reverse transcription‑quantitative PCR (RT‑qPCR)
Total RNA was isolated from the cells at the logarithmic growth phase using TRNzol‑A + reagent (cat. no. DP421; Tiangen Biotech Co., Ltd.), and converted to cDNA with TGFB1/PD-L1 primer and β-actin primers using the HiScript II 1st Strand cDNA Synthesis kit (cat. no. R211‑01; Vazyme). The relative expression level of TGFB1/PD-L1 was normalized to β-actin using the 2−ΔΔCt method. The sequences of the primers were listed as follows (5′ → 3′): TGFB1 forward, ACCAACTATTGCTTCAGCTC, reverse, TTATGCTGGTTGTACAGG; PD-L1 forward, GGCATTTGCTGAACGCAT, reverse, CAATTAGTGCAGCCAGGT; β-actin forward, GCCCATCTACGAGGGGTATG, reverse, GAGGTAGTCAGTCAGGTCCCG.
HepG2 cells were treated with 50 mM 5‑aza‑2′‑deoxycytidine (5‑aza‑dC; cat. no. A3656; Sigma‑Aldrich; Merck KGaA) for 72 h with a change of culture medium every 24 h as previously described .
Western blot analysis
Cells were lysed in a solution of 60 mM Tris–HCl, pH 6.8, 2.00% sodium dodecyl sulfate, 20.00% glycerol, 0.25% bromophenol blue, 1.25% 2-mercaptoethanol and heated at 100 °C for 10 min. Protein concentrations were determined using a BCA protein assay kit. After being separated by 10 or 12% SDS-PAGE, the protein (40-50 µg) was transferred to a PVDF membrane (cat. no. IPVH00010; EMD Millipore). The membranes were blocked in 1× PBS buffer containing 5% BSA (cat. no. A1933; Sigma‑Aldrich; Merck KGaA) and 0.05% Tween‑20 (cat. no. A100777; Sangon Biotech Co., Ltd.) for 1 h at room temperature, and then incubated with the following primary antibodies: HSD17B6 (rabbit anti‑human ployclonal; 1:1000; cat. no. 11855-1-AP; ProteinTech Group, Inc.); GAPDH (mouse anti‑human monoclonal; 1:2000; cat. no. 60004‑1‑Ig; ProteinTech Group, Inc.) overnight at 4 °C. The blots were washed with PBST three times for 10 min each and incubated with secondary antibodies anti‑rabbit IgG (1:3000; cat. no. SA00001‑2); and anti‑mouse IgG (1:3000, cat. no. SA00001‑1; both from ProteinTech Group, Inc.) for 1 h at room temperature. The target bands were revealed by SuperSignal West Pico PLUS chemiluminescence substrate (cat. no. 34580; Thermo Fisher Scientific, Inc.).
Data are presented as mean ± standard error of the mean (SEM), and the statistical significance of differences between groups was assessed by student’s t-test. The correlation of gene expression was evaluated by Spearman’s correlation and statistical significance. In all figures, values of p ≤ 0.05 were considered as statistically significant with *(p < 0.05), **(p < 0.01), ***(p < 0.001), ****(p < 0.0001). All the analyses were undertaken using PRISM version 8.02 (GraphPad software, version 8, San Diego, CA. USA).
HSD17B6 is down-regulated in patients with hepatocellular carcinoma
The male dominance in HCC occurrence suggested an important role of androgen and/or their receptors in the development of HCC . A previous report suggested that androgen signals might suppress the cell invasion during the later stages of HCC . HSD17B6 has both oxidoreductase and epimerase activities and is involved in androgen catabolism. It is mainly expressed in the liver (Fig. 1a) according to the GTEx, a data resource and tissue bank to study gene expression in multiple human normal tissues. To determine differences of HSD17B6 expression in tumor and normal tissues, the HSD17B6 mRNA levels in different tumors and normal tissues of multiple cancer types were analyzed in TCGA using the TIMER “DiffExp” module. HSD17B6 expression was significantly lower in most of tumor tissues, including LIHC (Liver hepatocellular carcinoma), compared to the normal tissues (p < 0.001, Fig. 1b). Higher expression of HSD17B6 was only observed in BRCA (Breast invasive carcinoma) and HNSC (Head and Neck squamous cell carcinoma) compared to the normal (Fig. 1b).
To further evaluate HSD17B6 expression in human liver cancers, the mRNA expression level of the HSD17B6 gene was examined using several other bioinformatics web resources. First, we analyzed four microarray datasets from the Oncomine Platform, and found that HSD17B6 transcript levels were significantly reduced in patients with HCC in all four datasets (p < 0.0001, Fig. 2a–d). Especially, in Wurmbach Liver and Chen Liver dataset, although HSD17B6 was downregulated in HCC versus normal liver tissue, its expression level in liver cancer precursors was not significantly reduced (Fig. 2c, d). HSD17B6 expression was also analyzed in multiple datasets with HCCDB tool, a one-stop online platform for exploring HCC gene expression . Totally, 9 of 10 datasets showed that HSD17B6 expression in HCC tissues was significantly decreased compared with in adjacent/healthy tissues (Table 1). This trend was further verified by pathologic analyses. The IHC staining results from the Human Protein Atlas project suggested that HSD17B6 was strongly expressed in normal liver specimens (negative rate of 0%) but weakly expressed or not detected in HCC tissues (negative rate of 42%) (Fig. 2e). These results collectively suggest that reduced expression of HSD17B6 may promote hepatocarcinogenesis, acting as a tumor suppressor.
HSD17B6 acts as a putative prognostic factor in hepatocellular carcinoma
To the best of our knowledge, there have been no clear reports of the relationship between HSD17B6 expression and the clinical prognosis of HCC. As shown in plots using the Kaplan–Meier Plotter platform, an intriguing TCGA analysis tool, down-regulation of HSD17B6 expression was closely correlated with shorter OS (p = 3.9e−5), PFS (p = 0.00024), RFS (p = 0.00022) and DSS (p = 0.00023) (Fig. 3a–d), indicating a poorer prognosis with lower HSD17B6 level in HCC. Shorter OS was also observed in patients with lower expression of HSD27B6 in other two HCC datasets (p = 0.0245 for GSE14520 and p = 0.0373 for ICGC-LIRI-JP, respectively), analyzed by HCCDB tool (Fig. 3e, f).
To investigate what alterations in HCC relate to HSD17B6, we analyzed the association between HSD17B6 mRNA expression levels and the clinical characteristics of HCC patients in TCGA and four liver disease microarray datasets from the Oncomine Platform. Decreased HSD17B6 mRNA expression was associated with increased tumor grade (Fig. 4a, data from TCGA and Fig. 5a, data from Wurmbach liver dataset) and pathological stage (Fig. 4b, data from TCGA; Fig. 5b, c, data from Wurmbach Liver and Jia liver datasets). It was also associated with worse Barcelona Clinic Liver Cancer Stage, which is the most common staging system for liver cancer, in Chiang Liver dataset (Fig. 5d). Lower expression of HSD17B6 was associated with worse T stage (Fig. 4c, data from TCGA, and Fig. 5e, data from Liao Liver dataset), which refers to the size and extent of the main tumor. In accordance with this, lower expression of HSD17B6 was associated with larger tumor size in Wurmbach liver dataset (Fig. 5f). Lower expression of HSD17B6 was also associated with worse N stage in TCGA (Fig. 4e), which refers to the degree of spread to regional lymph nodes. Accordantly, HCC samples with more vascular invasion/satellites had significantly lower expression of HSD17B6 (Fig. 5g, h, data from Wurmbach liver dataset). Moreover, its expression level was lower in metastatic tumor tissues than in primary HCC (Fig. 5i).
Collectively, these results provide insight into the functional role of HSD17B6 in HCC and suggest that the reduced expression of HSD17B6 in HCC patients might be a valuable prognostic factor.
Transcriptional regulation of HSD17B6 in hepatocellular carcinoma
In order to shed light on mechanism involved in the regulation of HSD17B6, we first tried to find out transcription factors regulating HSD17B6. It has been recently reported that HNF4A, a liver-enriched transcription factor, could upregulate gene expression in hepatocytes by binding the enhancers of target genes . There was a significant positive correlation between HSD17B6 and HNF4A transcript levels both in normal and HCC tissues from both TCGA and ICGC LIRI-JP datasets (Fig. 6a–d). Analysis of three liver disease microarray datasets from the Oncomine Platform confirmed this positive correlation (Additional file 1: Fig. S1). Additionally, analysis of HSD17B6 expression from GDS1916, a dataset form GEO database, showed that knockout of HNF4A in liver decreased the expression of HSD17B6 (Fig. 6e). Next, analyzing previously published, H3K4me1 , H3K4me3 , H3K27ac  and three HNF4A [27, 28] ChIP-seq datasets in HepG2 cells (a liver hepatocellular carcinoma cell line), we found the enrichment of modifications around HNF4A bound regions (HBRs) using CistromeDB [29, 30] (Fig. 6f). In particular, H3K4me1 and H3K4me3 mark enhancers and promoters respectively . And H3K27ac is associated with the active state of both elements. Analysis of HNF4A ChIP-seq data revealed its preferential binding at first intron and promoter of HSD17B6 gene (Fig. 6f). Examination of HBRs demonstrated a high enrichment of H3K4me1 and H3K27ac flanking these regions (Fig. 6f), suggestive of active enhancers. The active status of these enhancers was further confirmed by examination of a HepG2 DNase CHIP data  (Fig. 6f), which showed very high signals, demonstrative of open chromatin at HBRs (Fig. 6f).
DNA methylation is an epigenetic mechanism used by cells to control gene expression. We found that the methylation status of promoter of HSD17B6 was negatively correlated with expression of HSD17B6 in TCGA LIHC dataset (r = − 0.2543, Fig. 7a). Moreover, HSD17B6 expression was also negatively correlated with expression of DNA methyltransferases (DNMT1, DNMT3A and DNMT3B) in TCGA and ICGC LIRI-JP liver datasets (Fig. 7e–j). It has been reported that DNMT3B recruitment through E2F6 transcriptional repressor mediates DNA methylation and gene silencing . ChIP-seq results from six independent investigations, accessed using CistromeDB, showed E2F6 binding at the HSD17B6 promoter (Fig. 8c), which harbors a CpG island (Fig. 8c). HSD17B6 expression was also negatively correlated with expression of E2F6 in TCGA and ICGC LIRI-JP liver datasets (Fig. 8a, b). Together, these results suggest that E2F6 might repress HSD17B6 expression by recruiting DNMT3B to HSD17B6 promoter. In addition, treatment with 5-aza-2′-deoxycytidine, a DNA-hypomethylating agent, upregulated the mRNA expression of HSD17B6 in live cancer cell lines (HepG2, Hep3B and HuH-7) based on the data from GSE5230 and GSE35313 datasets (p < 0.05, Fig. 7b, c). Confirming this, treatment with 5-aza-2′-deoxycytidine increased the protein level of HSD17B6 in HepG2 cells (Fig. 7d).
Taken together, the expression of HSD17B6 in HCC is regulated by both HNF4A and DNA methylation.
HSD17B6 is Associated with the Androgen metabolic and biosynthetic processes
Previous studies have demonstrated that androgen (a steroid hormone) and/or their receptors play important roles in the development of HCC [22, 23]. The protein encoded by HSD17B6 has both oxidoreductase and epimerase activities and is involved in androgen catabolism. We then used the Search Tool for the Retrieval of Interacting Genes (STRING) database to search for functional interaction proteins with HSD17B6, and found that 10 most related interaction proteins, which were RD5A2, CYP11B1, SRD5A1, CYP11B2, CYP17A1, HSD3B2, HSD3B1, CYP19A1, SULT1E1 and AKR1D1 (Fig. 9a). Further GO (Biological Process) analysis showed that all 11 proteins were involved in hormone (steroid) metabolic process. Among these 11 proteins, 9 were involved in androgen metabolic process, and 10 in steroid biosynthetic process. Particularly, 6 proteins (HSD17B6, SRD5A1, SRD5A2, HSD3B1, HSD3B2, and CYP17A1) were involved in androgen biosynthetic process (a process consisting of 11 proteins) (Fig. 9b). GO (Molecular Function) analysis showed that 10 of 11 proteins had oxidoreductase activity, 5 had steroid dehydrogenase activity, and 4 had hydroxylase activity (Fig. 9c). All these results indicate that HSD17B6 plays a very important role in androgen biosynthesis and androgen signaling pathway in HCC. Moreover, GO (Cellular Component) analysis shown that 6 of 11 proteins were located in endoplasmic reticulum (Fig. 9d). As endoplasmic reticulum stress contributes to the development and progression of HCC , HSD17B6 and its interaction proteins are likely linked to this process.
HSD17B6 inhibits cell proliferation, migration, and invasion and correlates with immune response in HCC
To determine which biological processes (BP) are associated with HSD17B6, we searched for cellular signaling pathways altered in tumors with low HSD17B6 expression. We used Gene Set Enrichment Analysis (GSEA) to identify gene sets whose expression is enriched or depleted in low HSD17B6 expression tumors. This analysis was performed in TCGA LIHC and ICGC LIRI-JP datasets, and the common biological processes in these two datasets were used for further analysis.
Eight of the ten biological processes most negatively correlating with expression of HSD17B6 were those implicated in mitosis and cell cycle regulation, such as chromosome segregation, mitotic spindle organization, and microtubule cytoskeleton organization involved in mitosis (Fig. 10a, d, e and Additional file 2: Table S1b). The levels of cyclin B1 (CCNB1) and CDC20, two markers of mitosis, were also negatively correlated with HSD17B6 expression in TCGA LIHC and ICGC LIRI-JP (Additional file 3: Fig. S2). To further explore the biological role of HSD17B6 in liver cancer, HSD17B6 expression was stably up-regulated in HepG2 cell through transfection (Fig. 11a). Cell proliferation was evaluated by CCK-8 assay. As shown in Fig. 11b, cell proliferation was significantly inhibited by HSD17B6 overexpression (p < 0.05) (Fig. 11b). These results suggest that overexpression of HSD17B6 could inhibit the proliferation of liver cancer cells. Moreover, the effects of HSD17B6 on the migrated and invasive abilities of liver cancer cells were assessed by wound‑healing assay and transwell assay respectively. Overexpression of HSD17B6 in HepG2 cells caused a significant reduction in cell migration (Fig. 11c, d) and cell invasion (Fig. 11e, f).
Two distinct groups of biological processes positively correlating with HSD17B6 expression emerged from this analysis. First, the biological processes associated with metabolism was significantly decreased in low HSD17B6 expression tumors (Fig. 10a, b and Additional file 2: Table S1a). Seven of the ten biological processes most positively correlating with expression of HSD17B6 were those implicated in metabolism, such as fatty acid catabolic process, xenobiotic metabolic process and lipid oxidation (Fig. 10a, b). Second, four of the twenty biological processes most positively correlating with expression of HSD17B6 were those implicated in immunity response, including complement activation, regulation of humoral immune response, humoral immune response mediated by circulating immunoglobulin and acute phase response (Fig. 10a, c). These indicated that in HCC with low expression of HSD17B6, the immune response to tumors was reduced.
HSD17B6 is associated with tumor infiltrating immune cells and immune suppression in hepatocellular carcinoma
Tumor-infiltrating immune cells (TIICs) are a part of the complex microenvironment that regulates tumor development and progression . Since recent studies revealed that TIICs were preferentially enriched and potentially clonally expanded in HCC , and accumulation of TIICs correlated with poor prognosis for HCC , we quantified the TIICs from RNA sequencing data through two widely used bioinformatics tools: “TIMER” and “xCell” [38, 39] in TCGA LIHC and ICGC LIRI-JP datasets. “TIMER” could estimate the abundances of six immune cell types, including B-cells, CD4+ T-cells, CD8+ T-cells, dendritic cells, macrophages and neutrophils. The analysis showed that infiltration level of B-cells, CD4+ T-cells and dendritic cells were negatively correlated with HSD17B6 expression in both TCGA LIHC and ICGC LIRI-JP datasets (Fig. 12). Further analysis using “xCell” also showed that infiltration of 20 kinds of immune cells, including B-cells, CD4+ T-cells and dendritic cells, was negatively correlated with HSD17B6 expression (Fig. 13, Additional file 4: Fig. S3). However, infiltration of 3 kinds of immune cells, including conventional dendritic cells, macrophages, macrophages M2, was positively correlated with HSD17B6 expression (Fig. 13, Additional file 4: Fig. S3).
Tumors could evade immune responses by taking advantage of immune checkpoint genes, such as PD-1 and CTLA-4. We analyzed HSD17B6 and immune checkpoint genes expression in TCGA LIHC and ICGC LIRI-JP. HSD17B6 expression was negatively correlates with expression of PD-1(r = − 0.2861 in TCGA LIHC and r = − 0.3062 in ICGC), CTLA4(r = − 0.3207 in TCGA LIHC and r = − 0.4395 in ICGC), ICOS(r = − 0.1731 in TCGA LIHC and r = − 0.2704 in ICGC), B7-H3(r = − 0.3066 in TCGA LIHC and r = − 0.3307 in ICGC) and ENTPD1(r = − 0.3172 in TCGA LIHC and r = − 0.5825 in ICGC) (Fig. 14), suggesting that low HSD17B6 expression potentially plays an important role in mediating immune evasion.
Transforming growth factor-β (TGF-β) has long been shown to play an essential role in establishing immune suppression and evasion [40, 41]. It has been reported that treatment with DHT, a product of HSD17B6, inhibited the expression of TGFB1 in prostate cancer . Our analysis showed a significant negative correlation between HSD17B6 and TGFB1 transcript levels both in normal and HCC tissues from TCGA and ICGC LIRI-JP datasets (Fig. 15a–d). And both HSD17B6 overexpression and DHT treatment inhibited the expression of TGFB1 and PD-L1 in HepG2 cells (Fig. 15e, f).
This study was conducted to determine the feasibility of HSD17B6 as a potential biomarker in HCC patients. Our study is the first to analyze the differential expression of HSD17B6 in HCC tissues using different public datasets and examine its roles in cancer-related signaling pathways to identify its likely biological significance in HCC carcinogenesis. Using the TCGA, ICGC, Oncomine, GEO and HPA datasets, we have shown that HSD17B6 was significantly down-regulated in HCC tissues and that its expression level was negatively correlated with the OS, PFS, RFS and DSS of HCC patients. Using public CHIP datasets, we have proven that HSD17B6 was regulated by HNF4A and DNA methylation. Furthermore, we showed here that loss of HSD17B6 affected the androgen biosynthesis and androgen signaling pathway, promoted tumor cell proliferation, migration, invasion, immune cell infiltration and immune evasion partially through TGFB1.
HSD17B6 has both oxidoreductase and epimerase activities, involved in steroid metabolism. The oxidoreductase activity can convert 3 alpha-adiol to DHT, while the epimerase activity can convert androsterone to epi-androsterone [5, 6]. Previous studies have shown that polymorphisms of HSD17B6 gene are associated with polycystic ovary syndrome (PCOS) [8, 9]. For example, it was associated with reduced fasting glucose-insulin ratio and increased homeostasis model assessment and body mass index in the polycystic ovary syndrome group. The expression of HSD17B6 was significantly reduced in prostate cancer and it was undetectable in prostate cancers of Gleason grade higher than three . In men receiving androgen-deprived therapy (ADT), patients showing biochemical progression had a higher HSD17B6 score than those without progression, suggesting that lower expression of HSD17B6 correlated with worse prognosis . Its mRNA levels in castration resistant prostate cancer with bone metastatic samples were significantly lower than those without metastasis , indicating that HSD17B6 dysfunction involved in prostate cancer metastases. Downregulated HSD17B6 has also been identified to be associated with non-small-cell lung cancer , implying that it might be the key gene contributing to tumorigenesis or tumor progression in lung cancer . Moreover, SNP of HSD17B6 was significantly associated with liver fibrosis risk . In this study, using bioinformatics analyses of public datasets and experimentation, we have illustrated the key tumor suppressive roles of HSD17B6 in HCC samples (Fig. 15f).
DNA methylation status of CpGs, which is the best-characterized epigenetic mechanism, may play a pivotal role for downregulating HSD17B6 expression, as its expression was negatively correlated with the ratio of DNA methylation and expression of DNA methyltransferases (DNMT1, DNMT3A and DNMT3B) (Fig. 7e–j). And treatment with 5-aza, an inhibitor of DNA methylation, upregulated the HSD17B6 mRNA and protein expression of liver cancer cells (Fig. 7b–d). Actually, viral infection, chronic inflammation, and oxidative stress have been shown to affect DNMTs expression in HCC .Viral oncoproteins, such as the HCV core and the HBx protein of HBV, upregulate DNMT1 and DNMT3B . HBx can directly interact with DNMT3A to induce abnormal hypermethylation at tumor suppressor gene (TSG) promoters . Besides, HNF4A is also involved in regulation of HSD17B6, as their expression was positively correlated, and HNF4A binds the enhancer and promoter of HSD17B6. However, it is still unknown whether hypermethylation in the promoter region of HSD17B6 reduced the binding of HNF4A in HCC. Further studies are needed to more comprehensively explore the detailed molecular mechanisms of this altered biomarker in the progression of HCC.
In this study, biological pathway analyses had revealed clearly that the pathways related to androgen metabolic and biosynthetic processes in liver were the important pathways modulated by HSD17B6 (Fig. 9). Actually, cirrhotic patients with HCC had significantly lower plasma concentrations of several kinds of androgen including DHT, which was produced by HSD17B6, than patients with cirrhosis alone . Moreover, high levels of DHT could induced cell cycle arrest and apoptosis in human liver cells by activation of PKR/eIF2alpha signaling cascade , and DHT gives growth advantage to cells by inhibiting TGF-beta-induced apoptosis in rat hepatocellular carcinoma cells . However, another study showed that DHT could enhance HCC cell growth and apoptotic resistance . Some previous reports confirmed the key roles of androgen signaling in the etiology of human cancers, including liver cancer . Results from different studies indicated that androgen signaling might promote the initiation and development of HCC during the early stage , yet might suppress the cell invasion during the later stages of HCC . These results and our study suggested that HSD17B6 contributed to HCC development by affecting androgen signaling pathway. Further research is needed to elucidate the detailed mechanisms.
Recent researches have revealed that tumor-infiltrating immune cells (TIICs) could regulates tumor development and progression . Moreover, TIICs are preferentially enriched and potentially clonally expanded in HCC , and accumulation of TIICs correlates with poor prognosis for HCC . Accordingly, we found that HSD17B6 expression in liver cancer was significantly negatively correlated with multiple immune cell infiltration, such as dendritic cells, CD4+ T cells, and B cells (Figs. 12, 13). Interestingly, recent studies shown that androgen receptor expression inversely correlates with immune cell infiltration in breast cancer , and low level of DHT could promote the infiltration of CD8+ T cells in prostate of benign prostatic hyperplasia (BPH) tissues . What’s more, TGF-β promoted infiltration of immune cells and cancer-associated fibroblasts in the tumor microenvironments , and our present study showed that HSD17B6 could inhibit the expression of TGFB1 (Fig. 15). This may explain the negative correlation between the expression of HSD17B6 and the infiltration of immune cells.
The acquisition of genetic and epigenetic changes in cancer would lead to the expression of novel antigens, which are recognized by immune system, enabling the elimination of many nascent tumors . Thus, immune evasion is one of the hallmarks of tumors including hepatocellular carcinoma [56, 57]. The most studied and clinically relevant immunosuppressive factors implicated in HCC are CTLA-4, PD-1, PD-L1 and the immunosuppressive cytokine TGF-β [55, 58,59,60,61,62]. Immune checkpoint proteins, including CTLA-4, PD-1, and PD-L1, could inhibit the activities of multiple immune cells, such as CD8+ T cells, natural killer cells and dendritic cells [63, 64]. TGF-β is a crucial cytokine to inhibit the expansion and function of many components of the immune system. It is central to tumor immune evasion and poor responses to cancer immunotherapy [64,65,66]. For example, TGF-β could upregulate expression of some immune checkpoint genes, such as TIM-3, CTLA-4, B7-H3, ENTPD1(CD39) and PD-1/PD-L1 [67,68,69,70,71,72,73]. And TGF-β1 is the most relevant member of the TGF-β family for immune regulation . For example, TGF-β1 was shown to mediate tumor immune evasion from adaptive immunity . Our present study has shown that HSD17B6 overexpression or treatment with its product DHT could inhibit the expression of TGFB1 and PD-L1 (Fig. 15). Furthermore, HS17B6 expression was positively correlated with immune response pathways in HCC (Fig. 10a, c). These data indicated that low expression of HSD17B6 promoted immune evasion through upregulating the expression of TGFB1 in HCC tumor microenvironment.
Our results must be carefully explained with great caution. There are some limitations in the present study. First, we showed here that low expression of HSD17B6 led to high expression of TGFB1 through DHT, and then promoted tumor development in HCC. However, the protein encoded by HSD17B6 has both oxidoreductase and epimerase activities. The oxidoreductase activity can convert 3 alpha-adiol to DHT, while the epimerase activity can convert androsterone to epi-androsterone. Further experiments should be carried out to determine which genes other than TGFB1 are affected by DHT and epiandrosterone, and their role in hepatocarcinogenesis. Second, we used bioinformatics tools to calculate the proportion of various immune cells in liver cancer samples. Further experiments, such as immunohistochemistry (IHC), immune fluorescence (IF), and flow cytometry, are needed to confirm these calculations.
In the present study, we identified HSD17B6 as a promising molecular marker with prognostic value in HCC. Analysis of the public liver cancer data provides a meaningful method to screen key genes associated with the onset of human malignant diseases in the future. Moreover, HSD17B6, could be a novel biomarker for the biological behavior of HCC, serving as a putative tumor suppressor gene.
In conclusion, the data presented in this study show that HSD17B6 is frequently down-regulated in HCC. Moreover, our study indicated that transcription of HSD17B6 was positively regulated by hepatocyte nuclear factor 4 alpha (HNF4A), a liver-enriched transcription factor. Besides, DNA hypermethylation of the HSD17B6 promoter also contributed a lot to the low expression of HSD17B6 in HCC. HCC with low expression of HSD17B6 has worse tumor stage and prognosis. HSD17B6 may inhibit the expression of TGFB1 or other genes, and then inhibit tumor cell proliferation and invasion, immune cell infiltration and immune evasion (Fig. 15g). Taken together, this study demonstrates that HSD17B6 is an important player of liver carcinogenesis and helps better understanding of the molecular pathways involved in HCC development. The findings presented in this study add novel insight into the regulation of immune infiltration and suppression in HCC. We hope that it can help to make a better diagnosis and also help to predict the prognosis which might provide essential information regarding personalized treatment decisions for individual patients, especially for immune therapy.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Hydroxysteroid 17-beta dehydrogenase 6
The cancer genome atlas
International Cancer Genome Consortium
Gene expression omnibus
The human protein atlas
Polycystic ovary syndrome
Hepatocyte nuclear factor 4 alpha
The genotype tissue expression
Search tool for the retrieval of interacting genes
Tumor-infiltrating immune cells
Villanueva A. Hepatocellular carcinoma. N Engl J Med. 2019;380(15):1450–62.
Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. Cancer J Clin. 2018;68(6):394–424.
Jemal A, Ward EM, Johnson CJ, Cronin KA, Ma J, Ryerson B, Mariotto A, Lake AJ, Wilson R, Sherman RL, et al. Annual report to the nation on the status of cancer, 1975–2014, Featuring Survival. J Natl Cancer Inst. 2017;109:9.
Llovet JM, Zucman-Rossi J, Pikarsky E, Sangro B, Schwartz M, Sherman M, Gores G. Hepatocellular carcinoma. Nat Rev Dis Primers. 2016;2:16018.
Huang XF, Luu-The V. Gene structure, chromosomal localization and analysis of 3-ketosteroid reductase activity of the human 3(alpha– > beta)-hydroxysteroid epimerase. Biochem Biophys Acta. 2001;1520(2):124–30.
Huang XF, Luu-The V. Molecular characterization of a first human 3(alpha– > beta)-hydroxysteroid epimerase. J Biol Chem. 2000;275(38):29452–7.
Chan YX, Yeap BB. Dihydrotestosterone and cancer risk. Curr Opin Endocrinol Diabetes Obes. 2018;25(3):209–17.
Jones MR, Italiano L, Wilson SG, Mullin BH, Mead R, Dudbridge F, Watts GF, Stuckey BG. Polymorphism in HSD17B6 is associated with key features of polycystic ovary syndrome. Fertil Steril. 2006;86(5):1438–46.
Jones MR, Mathur R, Cui J, Guo X, Azziz R, Goodarzi MO. Independent confirmation of association between metabolic phenotypes of polycystic ovary syndrome and variation in the type 6 17beta-hydroxysteroid dehydrogenase gene. J Clin Endocrinol Metabol. 2009;94(12):5034–8.
Jernberg E, Thysell E, Bovinder Ylitalo E, Rudolfsson S, Crnalic S, Widmark A, Bergh A, Wikstrom P. Characterization of prostate cancer bone metastases according to expression levels of steroidogenic enzymes and androgen receptor splice variants. PLoS ONE. 2013;8(11):e77407.
Ma Q, Xu Y, Liao H, Cai Y, Xu L, Xiao D, Liu C, Pu W, Zhong X, Guo X. Identification and validation of key genes associated with non-small-cell lung cancer. J Cell Physiol. 2019;234(12):22742–52.
White DL, Liu Y, Garcia J, El-Serag HB, Jiao L, Tsavachidis S, Franco LM, Lee JS, Tavakoli-Tabasi S, Moore D, et al. Sex hormone pathway gene polymorphisms are associated with risk of advanced hepatitis C-related liver disease in males. Int J Mol Epidemiol Genet. 2014;5(3):164–76.
Nagy A, Lanczky A, Menyhart O, Gyorffy B. Validation of miRNA prognostic power in hepatocellular carcinoma using expression data of independent datasets. Sci Rep. 2018;8(1):9227.
Cancer Genome Atlas N. Comprehensive molecular portraits of human breast tumours. Nature. 2012;490(7418):61–70.
Uhlen M, Fagerberg L, Hallstrom BM, Lindskog C, Oksvold P, Mardinoglu A, Sivertsson A, Kampf C, Sjostedt E, Asplund A, et al. Proteomics. Tissue-based map of the human proteome. Science. 2015;347(6220):1260419.
Thul PJ, Akesson L, Wiking M, Mahdessian D, Geladaki A, Ait Blal H, Alm T, Asplund A, Bjork L, Breckels LM, et al. A subcellular map of the human proteome. Science. 2017;356:6340.
Uhlen M, Zhang C, Lee S, Sjostedt E, Fagerberg L, Bidkhori G, Benfeitas R, Arif M, Liu Z, Edfors F, et al. A pathology atlas of the human cancer transcriptome. Science. 2017;357:6352.
Szklarczyk D, Morris JH, Cook H, Kuhn M, Wyder S, Simonovic M, Santos A, Doncheva NT, Roth A, Bork P, et al. The STRING database in 2017: quality-controlled protein-protein association networks, made broadly accessible. Nucleic Acids Res. 2017;45(D1):D362–8.
Lian Q, Wang S, Zhang G, Wang D, Luo G, Tang J, Chen L, Gu J. HCCDB: a database of hepatocellular carcinoma expression atlas. Genomics Proteomics Bioinf. 2018;16(4):269–75.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005;102(43):15545–50.
Zhang JX, Chen ZH, Xu Y, Chen JW, Weng HW, Yun M, Zheng ZS, Chen C, Wu BL, Li EM, et al. Downregulation of MicroRNA-644a promotes esophageal squamous cell carcinoma aggressiveness and stem cell-like phenotype via dysregulation of PITX2. Clin Cancer Res. 2017;23(1):298–310.
Li Z, Tuteja G, Schug J, Kaestner KH. Foxa1 and Foxa2 are essential for sexual dimorphism in liver cancer. Cell. 2012;148(1–2):72–83.
Ma WL, Hsu CL, Yeh CC, Wu MH, Huang CK, Jeng LB, Hung YC, Lin TY, Yeh S, Chang C. Hepatic androgen receptor suppresses hepatocellular carcinoma metastasis through modulation of cell migration and anoikis. Hepatology. 2012;56(1):176–85.
Thakur A, Wong JCH, Wang EY, Lotto J, Kim D, Cheng JC, Mingay M, Cullum R, Moudgil V, Ahmed N, et al. Hepatocyte nuclear factor 4-alpha is essential for the active epigenetic state at enhancers in mouse liver. Hepatology. 2019;70(4):1360–76.
Ernst J, Kheradpour P, Mikkelsen TS, Shoresh N, Ward LD, Epstein CB, Zhang X, Wang L, Issner R, Coyne M, et al. Mapping and analysis of chromatin state dynamics in nine human cell types. Nature. 2011;473(7345):43–9.
Thurman RE, Rynes E, Humbert R, Vierstra J, Maurano MT, Haugen E, Sheffield NC, Stergachis AB, Wang H, Vernot B, et al. The accessible chromatin landscape of the human genome. Nature. 2012;489(7414):75–82.
Gertz J, Savic D, Varley KE, Partridge EC, Safi A, Jain P, Cooper GM, Reddy TE, Crawford GE, Myers RM. Distinct properties of cell-type-specific and shared transcription factor binding sites. Mol Cell. 2013;52(1):25–36.
Schmidt D, Schwalie PC, Ross-Innes CS, Hurtado A, Brown GD, Carroll JS, Flicek P, Odom DT. A CTCF-independent role for cohesin in tissue-specific transcription. Genome Res. 2010;20(5):578–88.
Mei S, Qin Q, Wu Q, Sun H, Zheng R, Zang C, Zhu M, Wu J, Shi X, Taing L, et al. Cistrome data browser: a data portal for ChIP-Seq and chromatin accessibility data in human and mouse. Nucleic Acids Res. 2017;45(D1):D658–62.
Zheng R, Wan C, Mei S, Qin Q, Wu Q, Sun H, Chen CH, Brown M, Zhang X, Meyer CA, et al. Cistrome data browser: expanded datasets and new tools for gene regulatory analysis. Nucleic Acids Res. 2019;47(D1):D729–35.
Robertson AG, Bilenky M, Tam A, Zhao Y, Zeng T, Thiessen N, Cezard T, Fejes AP, Wederell ED, Cullum R, et al. Genome-wide relationship between histone H3 lysine 4 mono- and tri-methylation and transcription factor binding. Genome Res. 2008;18(12):1906–17.
Mercer TR, Neph S, Dinger ME, Crawford J, Smith MA, Shearwood AM, Haugen E, Bracken CP, Rackham O, Stamatoyannopoulos JA, et al. The human mitochondrial transcriptome. Cell. 2011;146(4):645–58.
Velasco G, Hube F, Rollin J, Neuillet D, Philippe C, Bouzinba-Segard H, Galvani A, Viegas-Pequignot E, Francastel C. Dnmt3b recruitment through E2F6 transcriptional repressor mediates germ-line gene silencing in murine somatic tissues. Proc Natl Acad Sci USA. 2010;107(20):9281–6.
Maurel M, Samali A, Chevet E. Endoplasmic reticulum stress: at the crossroads of inflammation and metabolism in hepatocellular carcinoma development. Cancer Cell. 2014;26(3):301–3.
Domingues P, Gonzalez-Tablas M, Otero A, Pascual D, Miranda D, Ruiz L, Sousa P, Ciudad J, Goncalves JM, Lopes MC, et al. Tumor infiltrating immune cells in gliomas and meningiomas. Brain Behav Immun. 2016;53:1–15.
Zheng C, Zheng L, Yoo JK, Guo H, Zhang Y, Guo X, Kang B, Hu R, Huang JY, Zhang Q, et al. Landscape of infiltrating T cells in liver cancer revealed by single-cell sequencing. Cell. 2017;169(7):1342–56.
Sun H, Liu L, Huang Q, Liu H, Huang M, Wang J, Wen H, Lin R, Qu K, Li K, et al. Accumulation of tumor-infiltrating CD49a(+) NK cells correlates with poor prognosis for human hepatocellular carcinoma. Cancer Immunol Res. 2019;7(9):1535–46.
Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, Li B, Liu XS. TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res. 2017;77(21):e108–10.
Aran D, Hu Z, Butte AJ. xCell: digitally portraying the tissue cellular heterogeneity landscape. Genome Biol. 2017;18(1):220.
Sia D, Jiao Y, Martinez-Quetglas I, Kuchuk O, Villacorta-Martin C, Castro de Moura M, Putra J, Camprecios G, Bassaganyas L, Akers N, et al. Identification of an immune-specific class of hepatocellular carcinoma. Based on molecular features. Gastroenterology. 2017;153(3):812–26.
Batlle E, Massague J. Transforming growth factor-beta signaling in immunity and cancer. Immunity. 2019;50(4):924–40.
Chen WY, Tsai YC, Yeh HL, Suau F, Jiang KC, Shao AN, Huang J, Liu YN. Loss of SPDEF and gain of TGFBI activity after androgen deprivation therapy promote EMT and bone metastasis of prostate cancer. Sci Signal. 2017;10:492.
Muthusamy S, Andersson S, Kim HJ, Butler R, Waage L, Bergerheim U, Gustafsson JA. Estrogen receptor beta and 17beta-hydroxysteroid dehydrogenase type 6, a growth regulatory pathway that is lost in prostate cancer. Proc Natl Acad Sci USA. 2011;108(50):20090–4.
Ishizaki F, Nishiyama T, Kawasaki T, Miyashiro Y, Hara N, Takizawa I, Naito M, Takahashi K. Androgen deprivation promotes intratumoral synthesis of dihydrotestosterone from androgen metabolites in prostate cancer. Sci Rep. 2013;3:1528.
Hlady RA, Tiedemann RL, Puszyk W, Zendejas I, Roberts LR, Choi JH, Liu C, Robertson KD. Epigenetic signatures of alcohol abuse and hepatitis infection during human hepatocarcinogenesis. Oncotarget. 2014;5(19):9425–43.
Guechot J, Peigney N, Ballet F, Vaubourdolle M, Giboudeau J, Poupon R. Sex hormone imbalance in male alcoholic cirrhotic patients with and without hepatocellular carcinoma. Cancer. 1988;62(4):760–2.
Dai R, Yan D, Li J, Chen S, Liu Y, Chen R, Duan C, Wei M, Li H, He T. Activation of PKR/eIF2alpha signaling cascade is associated with dihydrotestosterone-induced cell cycle arrest and apoptosis in human liver cells. J Cell Biochem. 2012;113(5):1800–8.
Lim IK, Joo HJ, Choi KS, Sueoka E, Lee MS, Ryu MS, Fujiki H. Protection of 5alpha-dihydrotestosterone against TGF-beta-induced apoptosis in FaO cells and induction of mitosis in HepG2 cells. Int J Cancer. 1997;72(2):351–5.
Jie X, Lang C, Jian Q, Chaoqun L, Dehua Y, Yi S, Yanping J, Luokun X, Qiuping Z, Hui W, et al. Androgen activates PEG10 to promote carcinogenesis in hepatic cancer cells. Oncogene. 2007;26(39):5741–51.
Ma WL, Lai HC, Yeh S, Cai X, Chang C. Androgen receptor roles in hepatocellular carcinoma, fatty liver, cirrhosis and hepatitis. Endocr Relat Cancer. 2014;21(3):R165–82.
Ma WL, Hsu CL, Wu MH, Wu CT, Wu CC, Lai JJ, Jou YS, Chen CW, Yeh S, Chang C. Androgen receptor is a new potential therapeutic target for the treatment of hepatocellular carcinoma. Gastroenterology. 2008;135(3):947–55.
van Rooijen JM, Qiu SQ, Timmer-Bosscha H, van der Vegt B, Boers JE, Schroder CP, de Vries EGE. Androgen receptor expression inversely correlates with immune cell infiltration in human epidermal growth factor receptor 2-positive breast cancer. Eur J Cancer. 2018;103:52–60.
Fan Y, Hu S, Liu J, Xiao F, Li X, Yu W, Cui Y, Sun M, Lv T, He Q, et al. Low intraprostatic DHT promotes the infiltration of CD8+ T cells in BPH tissues via modulation of CCL5 secretion. Mediat Inflamm. 2014;2014:397815.
Yang L, Pang Y, Moses HL. TGF-beta and immune cells: an important regulatory axis in the tumor microenvironment and progression. Trends Immunol. 2010;31(6):220–7.
Pardoll DM. The blockade of immune checkpoints in cancer immunotherapy. Nat Rev Cancer. 2012;12(4):252–64.
Hanahan D, Weinberg RA. Hallmarks of cancer: the next generation. Cell. 2011;144(5):646–74.
Ihling C, Naughton B, Zhang Y, Rolfe PA, Frick-Krieger E, Terracciano LM, Dussault I. Observational study of PD-L1, TGF-beta, and immune cell infiltrates in hepatocellular carcinoma. Front Med (Lausanne). 2019;6:15.
Jung HI, Jeong D, Ji S, Ahn TS, Bae SH, Chin S, Chung JC, Kim HC, Lee MS, Baek MJ. Overexpression of PD-L1 and PD-L2 is associated with poor prognosis in patients with hepatocellular carcinoma. Cancer Res Treat. 2017;49(1):246–54.
Katz LH, Likhter M, Jogunoori W, Belkin M, Ohshiro K, Mishra L. TGF-beta signaling in liver and gastrointestinal cancers. Cancer Lett. 2016;379(2):166–72.
Wang J, Li M, Wang Y, Liu X. Integrating subpathway analysis to identify candidate agents for hepatocellular carcinoma. Onco Targets Ther. 2016;9:1221–30.
Agdashian D, ElGindi M, Xie C, Sandhu M, Pratt D, Kleiner DE, Figg WD, Rytlewski JA, Sanders C, Yusko EC, et al. The effect of anti-CTLA4 treatment on peripheral and intra-tumoral T cells in patients with hepatocellular carcinoma. Cancer Immunol Immunother. 2019;68(4):599–608.
Chen J, Gingold JA, Su X. Immunomodulatory TGF-beta signaling in hepatocellular carcinoma. Trends Mol Med. 2019;25(11):1010–23.
Tiegs G, Lohse AW. Immune tolerance: what is unique about the liver. J Autoimmun. 2010;34(1):1–6.
Lugade AA, Kalathil S, Miller A, Iyer R, Thanavala Y. High immunosuppressive burden in advanced hepatocellular carcinoma patients: can effector functions be restored? Oncoimmunology. 2013;2(7):e24679.
Chen CH, Seguin-Devaux C, Burke NA, Oriss TB, Watkins SC, Clipstone N, Ray A. Transforming growth factor beta blocks Tec kinase phosphorylation, Ca2+ influx, and NFATc translocation causing inhibition of T cell differentiation. J Exp Med. 2003;197(12):1689–99.
Gorelik L, Constant S, Flavell RA. Mechanism of transforming growth factor beta-induced inhibition of T helper type 1 differentiation. J Exp Med. 2002;195(11):1499–505.
Wiener Z, Kohalmi B, Pocza P, Jeager J, Tolgyesi G, Toth S, Gorbe E, Papp Z, Falus A. TIM-3 is expressed in melanoma cells and is upregulated in TGF-beta stimulated mast cells. J Invest Dermatol. 2007;127(4):906–14.
Zhou X, Mao Y, Zhu J, Meng F, Chen Q, Tao L, Li R, Fu F, Liu C, Hu Y, et al. TGF-beta1 promotes colorectal cancer immune escape by elevating B7-H3 and B7-H4 via the miR-155/miR-143 axis. Oncotarget. 2016;7(41):67196–211.
Peres RS, Donate PB, Talbot J, Cecilio NT, Lobo PR, Machado CC, Lima KWA, Oliveira RD, Carregaro V, Nakaya HI, et al. TGF-beta signalling defect is linked to low CD39 expression on regulatory T cells and methotrexate resistance in rheumatoid arthritis. J Autoimmun. 2018;90:49–58.
Chen W, Jin W, Hardegen N, Lei KJ, Li L, Marinos N, McGrady G, Wahl SM. Conversion of peripheral CD4+ CD25− naive T cells to CD4+ CD25+ regulatory T cells by TGF-beta induction of transcription factor Foxp3. J Exp Med. 2003;198(12):1875–86.
Park HB, Paik DJ, Jang E, Hong S, Youn J. Acquisition of anergic and suppressive activities in transforming growth factor-beta-costimulated CD4+ CD25− T cells. Int Immunol. 2004;16(8):1203–13.
Zheng SG, Gray JD, Ohtsuka K, Yamagiwa S, Horwitz DA. Generation ex vivo of TGF-beta-producing regulatory T cells from CD4+ CD25− precursors. J Immunol. 2002;169(8):4183–9.
Song S, Yuan P, Wu H, Chen J, Fu J, Li P, Lu J, Wei W. Dendritic cells with an increased PD-L1 by TGF-beta induce T cell anergy for the cytotoxicity of hepatocellular carcinoma cells. Int Immunopharmacol. 2014;20(1):117–23.
Donkor MK, Sarkar A, Savage PA, Franklin RA, Johnson LK, Jungbluth AA, Allison JP, Li MO. T cell surveillance of oncogene-induced prostate cancer is impeded by T cell-derived TGF-beta1 cytokine. Immunity. 2011;35(1):123–34.
This work was supported by the National Natural Science Foundation of China (81402327 to QYY, 81602230 to LL, 31600680 to YZ).
Ethics approval and consent to participate
The Ethics approval and consent to participate of the current study was approved and consented by the ethics committee of Anhui Medical University.
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Correlation of HSD17B6 expression with HNF4A expression in three Oncomine liver datasets.
List of top 20 common biological processes correlating with expression of HSD17B6 in TCGA LIHC and ICGC LIRI-JP datasets.
Correlation of HSD17B6 expression with CCNB1/CDC20 expression in TCGA LIHC and ICGC-LIRI-JP datasets.
Correlation of HSD17B6 expression with infiltration level of immune cells by xCell in ICGC LIRI-JP dataset.
About this article
Cite this article
Lv, L., Zhao, Y., Wei, Q. et al. Downexpression of HSD17B6 correlates with clinical prognosis and tumor immune infiltrates in hepatocellular carcinoma. Cancer Cell Int 20, 210 (2020). https://doi.org/10.1186/s12935-020-01298-5
- DNA methylation
- Hepatocellular carcinoma (HCC)
- Immune suppression