Promoter Demethylation Upregulates STEAP1 Gene Expression in Human Prostate Cancer: In Vitro and In Silico Analysis

The Six Transmembrane Epithelial Antigen of the Prostate (STEAP1) is an oncogene overexpressed in several human tumors, particularly in prostate cancer (PCa). However, the mechanisms involved in its overexpression remain unknown. It is well known that epigenetic modifications may result in abnormal gene expression patterns, contributing to tumor initiation and progression. Therefore, this study aimed to analyze the methylation pattern of the STEAP1 gene in PCa versus non-neoplastic cells. Bisulfite amplicon sequencing of the CpG island at the STEAP1 gene promoter showed a higher methylation level in non-neoplastic PNT1A prostate cells than in human PCa samples. Bioinformatic analysis of the GEO datasets also showed the STEAP1 gene promoter as being demethylated in human PCa, and a negative association with STEAP1 mRNA expression was observed. These results are supported by the treatment of non-neoplastic PNT1A cells with DNMT and HDAC inhibitors, which induced a significant increase in STEAP1 mRNA expression. In addition, the involvement of HDAC in the regulation of STEAP1 mRNA expression was corroborated by a negative association between STEAP1 mRNA expression and HDAC4,5,7 and 9 in human PCa. In conclusion, our work indicates that STEAP1 overexpression in PCa can be driven by the hypomethylation of STEAP1 gene promoter.


Introduction
Prostate cancer (PCa) is the second most common cancer worldwide, and it is the sixth leading cause of cancer death globally [1]. Prostatic carcinogenesis arises from precursor preneoplastic lesions that have the distinct molecular characteristics of normal prostate cells, and may give rise to localized cancer and eventually acquire the potential of invasion and metastasis [2]. The molecular alterations underpinning these cellular events include mutations, gene deletions and amplifications, chromosomal rearrangements and epigenetic Life 2021, 11, 1251 2 of 12 modifications in key genes controlling cell fate [3][4][5]. The main epigenetic mechanisms are DNA methylation and histone modifications, the main function of which is to ensure the proper regulation of gene expression by changing the chromatin structure [6]. DNA methylation is catalyzed by the family of enzymes known as DNA methyltransferases (DNMTs) and by histone modifications, including mainly histone acetylation, which is regulated by two groups of enzymes exerting opposite effects, histone acetyltransferases (HATs) and histone deacetylases (HDACs) [5,6].
In tumor cells, hypermethylation is observed in promoters of specific genes, particularly tumor suppressor genes, and a global hypomethylation contributes to genomic instability and the activation of oncogenes [7]. The disruption of epigenetic mechanisms may conduct the deregulation of gene expression, leading to tumor development and progression [6].
The Six Transmembrane Epithelial Antigen of the Prostate 1 (STEAP1) gene was identified as overexpressed in PCa compared to non-malignant tissues [8]. The STEAP1 protein is mainly located in the plasma membrane of epithelial cells, particularly at cell-cell junctions where it may act as an ion channel or transporter protein [8]. In fact, it was reported that STEAP1 may allow the transport of small molecules between adjacent cells, indicating that STEAP1 may be involved in intercellular communication [9,10]. Several studies have demonstrated the role of STEAP1 in cancer. In Ewing tumors, STEAP1 protein seems to promote cell growth and invasiveness by increasing intracellular reactive oxygen species levels. The oxidative stress that results from STEAP1 overexpression may enhance tumor aggressiveness through the activation of genes involved in cell proliferation and invasion [11]. Additionally, in gastric tumors, it was demonstrated that the upregulation of STEAP1 increased cell proliferation, migration and invasion [12]. In human ovarian and lung cancers, the STEAP1 gene is highly expressed, and it is associated with metastasis and epithelial-mesenchymal transition [12][13][14]. It was also demonstrated that the knockdown of STEAP1 expression on prostate tumor cells is associated with antitumor effects, such as enhanced apoptosis, and reduced proliferation, migration and invasion [12][13][14][15]. Previously, it was shown that post-transcriptional and post-translational modifications may contribute to STEAP1 overexpression, as STEAP1 mRNA and protein stability are higher in neoplastic LNCaP cells than in non-neoplastic PNT1A cells [16]. However, these alterations do not justify the overexpression of STEAP1 in tumor cells, suggesting that other mechanisms may be involved. Thus, we hypothesized that epigenetic alterations of the STEAP1 gene result in its overexpression in PCa. The present study aimed to analyze the methylation pattern of the STEAP1 gene in PCa cells. Therefore, we analyzed the methylation levels of the STEAP1 gene promoter in prostate cell lines and human samples of PCa. In addition, we analyzed the association between methylation levels of STEAP1 and gene expression using publicly available datasets. Additionally, non-neoplastic PNT1A cells were used to demonstrate that demethylation of the STEAP1 gene, as well as a synergistic effect between DNA demethylation and inhibition of class I and II HDACs, induces STEAP1 overexpression.

Prostate Cell Lines and Treatments
The human LNCaP PCa cell line and the immortalized non-neoplastic PNT1A prostate epithelial cell line were purchased from the European Collection of Cell Cultures (ECACC, Salisbury, UK). LNCaP and PNT1A cell lines were cultured at 37 • C in a 5% CO 2 atmosphere with RPMI 1640 phenol-red medium (Sigma-Aldrich, St. Louis, MO, USA) supplemented with 10% fetal bovine serum (FBS, Sigma-Aldrich, USA) and 1% penicillin/streptomycin (Sigma-Aldrich, USA). For the treatment with the demethylation and the histone deacetylation drug (epidrugs), approximately 3 × 10 5 PNT1A cells were seeded in six-well plates until reaching about 60% confluence. After that, PNT1A cells were exposed to one treatment with 5 µM 5-aza-dC (Sigma-Aldrich, USA) for 72 h, and the other with 5 µM 5-aza-dC for 48 h followed by 24 h with 1 µM TSA (Sigma-Aldrich, USA). In the control group, the medium was replaced by RPMI medium with DMSO for 72 h.

Patients and Tissue Sample Collection
Prostate tissue samples from five patients diagnosed with clinically localized PCa and primary treatment with radical prostatectomy, at the Portuguese Oncology Institute of Porto (IPO-Porto), were used in this study. Informed consent was obtained from all participants, according to institutional regulations. This study was approved by the institutional review board (Comissão de Ética para a Saúde-(IRB-CES-IPOFG-EPE 019/08)) of IPO-Porto.

DNA Extraction and Bisulfite Conversion
DNA extraction from PNT1A and LNCaP cells and clinical samples was carried out using the Gentra Puregene Cell Kit (Qiagen, Hilden, Germany) according to the manufacturer's instructions. To evaluate the methylation pattern, 1 µg of genomic DNA was modified using the EZ DNA Methylation-Gold kit (ZYMO RESEARCH, Irvine, CA, USA) according to the manufacturer's instructions. The modified DNA was stored at −80 • C.

Polymerase Chain Reaction (PCR) Amplification and Cloning Products
PCR reactions were performed using 200 ng of bisulfite-modified DNA in 25 µL reaction containing 1 U of TrueStart Hot Start Taq DNA Polymerase (Thermo Scientific, Waltham, MA, USA), 2.5 mM of MgCl 2 , 10 mM dNTPs and 300 nM of each primer (−338 fw/+74 rv). The primer sequences and characteristics are described in Table 1. The PCR products were purified using NucleoSpin Gel and the PCR Clean-up kit (Macherey-Nagel, Germany) according to the manufacturer's instructions, cloned into a pNZY28 vector and transformed in NZYStar Competent Cells. After heat shock, the cells were plated onto LB agar plates containing 100 µg/mL ampicillin, 80 µg/mL X-gal and 0.5 mM IPTG and incubated at 37 • C overnight. All components used in cloning and transformation were purchased from Nzytech, Portugal.

DNA Sequencing
Colony screening was performed by PCR, with standard vector primers (T7 and M13), to confirm and amplify the DNA insert. Thus, white colonies were selected and incubated in 10 µL TE buffer at 100 • C for 2 min. Afterwards, 1 µL was used for PCR reaction with Speedy NZYTaq 2x Green Master Mix (Nzytech, Portugal). Sequencing of PCR products was carried out using the CEQ Dye Terminator Cycle Sequencing Quick Start Kit (Beckman Coulter, Fullerton, USA) according to the manufacturer's instructions. For each sample, the DNA sequencing reaction was performed in both strands with universal sequencing primers (T7 and M13). The sequencing products were separated on an automated capillary DNA sequencer (GenomeLabTM GeXP, Genetic Analysis System; Beckman Coulter, Fullerton, CA, USA). The sequencing data analysis was performed using the Clustal Omega software to align the PCR product sequences with the sequence of STEAP1 gene modified DNA.

Datasets and Bioinformatic Analysis
Three PCa datasets (GSE52955, GSE76938 and GSE38240) were downloaded from the public repository NCBI Gene Expression Omnibus (GEO) databases (https://www.ncbi. nlm.nih.gov/geo/, assessed on 12 November 2021). All these datasets were based on the GPL13534 platform (Illumina HumanMethylation450 BeadChip). For each dataset, only the samples associated with prostate were selected. Methylation status was determined through the interactive web tool GEO2R, which allows the comparison of the two groups defined, pathologic condition and normal. The methylation status varies between 0 and 1, where 0 means a low degree of methylation and 1 indicates a high degree of methylation. The results were exported, and graphs were constructed with GraphPad Prism 8.0.1. The details of each dataset used in the present study are described in Table 2. The correlation of DNA methylation status with STEAP1 mRNA expression was assessed using another public repository, the Prostate Adenocarcinoma (TCGA, Cell 2015) [17] dataset, available from cBioPortal for Cancer Genomic (https://www.cbioportal.org/, assessed on 12 November 2021). This dataset comprises data from 333 primary prostate carcinomas. * Obtained from patients submitted to cystoprostatectomy due to bladder cancer. # Obtained from organ donor with no evidence of prostate cancer.

Total RNA Extraction, cDNA Synthesis and Quantitative Real-Time PCR (qPCR)
Total RNA extraction from PNT1A cells treated with epidrugs (5-aza-dC and TSA) was carried out using TRI reagent (Sigma-Aldrich, USA) according to the manufacturer's instructions. cDNA synthesis was performed using the NZY First-Strand cDNA Synthesis KIT (Nzytech, Portugal) according to the manufacturer's instructions, and qPCR was carried out to evaluate the expression of STEAP1 mRNA (hSTEAP1) in PNT1A cells treated with 5-aza-dC alone and 5-aza-dC plus TSA. To normalize the expression of the STEAP1 gene, human GAPDH (hGAPDH) and human beta-2-microglobulin (hβ2M) primers were used as internal controls. qPCR reactions were carried out using 1 µL of cDNA synthesized in a 20 µL reaction containing 10 µL of Maxima SYBR Green/Fluorescein qPCR Master Mix (Thermo Scientific) and 300 nM of primer for each gene. Fold differences were calculated following the mathematical model proposed by Pfaffl using the formula: 2-(∆∆Ct) [21]. The primer sequence for each gene and respective amplicon sizes used in qPCR are described in Table 1.

Statistical Analysis
Data analysis was performed using GraphPad Prism version 8.0.1. for Windows (GraphPad Software, San Diego, CA, USA). The statistical significance of differences in STEAP1 mRNA expression for the treatment with 5-aza-dC and TSA in non-neoplastic PNT1A cells was assessed by student's t-test. Significant differences were considered when p < 0.05 compared to control values. All experimental data are shown as mean ± SEM.

Methylation Analysis of STEAP1 Gene in Neoplastic Tissue/Cells Compared with Non-Neoplastic Cells
To evaluate if DNA methylation plays a role in STEAP1 gene regulation and if there are alterations in PCa, the methylation pattern of STEAP1 was determined in PCa tissue samples, LNCaP and PNT1A cells. For this purpose, cytosine-rich regions of STEAP1 gene promoter and primer design were performed using the Methyl Primer Express Software v1.0 (Applied Biosystems). This analysis indicated that part of the promoter region and the first exon of the STEAP1 gene contain a large CpG island containing 24 CpG dinucleotides (Figure 1), which could provide a large number of sites for the methylation modification of this gene. are alterations in PCa, the methylation pattern of STEAP1 was determined in PCa tissue samples, LNCaP and PNT1A cells. For this purpose, cytosine-rich regions of STEAP1 gene promoter and primer design were performed using the Methyl Primer Express Software v1.0 (Applied Biosystems). This analysis indicated that part of the promoter region and the first exon of the STEAP1 gene contain a large CpG island containing 24 CpG dinucleotides (Figure 1), which could provide a large number of sites for the methylation modification of this gene. The methylation pattern of the STEAP1 was analyzed through the bisulfite sequencing PCR method from position −338 (promoter region) to +74 (exon 1). The analysis of the methylation pattern of STEAP1 revealed some differences between neoplastic and nonneoplastic samples. In non-neoplastic PNT1A cells, the results show that some of the CpG dinucleotides located in the promoter region are methylated, but in PCa tissue or in neoplastic LNCaP cells the CpG dinucleotides are completely demethylated (Figure 2). The methylation pattern of the STEAP1 was analyzed through the bisulfite sequencing PCR method from position −338 (promoter region) to +74 (exon 1). The analysis of the methylation pattern of STEAP1 revealed some differences between neoplastic and nonneoplastic samples. In non-neoplastic PNT1A cells, the results show that some of the CpG dinucleotides located in the promoter region are methylated, but in PCa tissue or in neoplastic LNCaP cells the CpG dinucleotides are completely demethylated (Figure 2).

Analysis of STEAP1 Promoter Methylation Levels in PCa and Normal Prostate Tissues from the GEO Database
In order to support and validate the results above, the methylation pattern of STEAP1 gene promoter methylation was evaluated in datasets from public PCa databases. The GSE52955, GSE76938 and GSE38240 datasets were downloaded and analyzed by GEO2R online software. Four CpG probes on the STEAP1 gene promoter were selected: cg15089950, cg19317433, cg19532731 and cg24286372 located on −314 to −193 bp, −285 to −163 bp, −301 to −180 bp and −250 to −129 bp upstream of the transcription start, respectively. As shown in Figure 3, there were significant differences in the levels of CpG methylation between prostate tumor and normal tissue. In three datasets analyzed, the results reveal lower methylation levels of the STEAP1 gene promoter in neoplastic tissue than in

Analysis of STEAP1 Promoter Methylation Levels in PCa and Normal Prostate Tissues from the GEO Database
In order to support and validate the results above, the methylation pattern of STEAP1 gene promoter methylation was evaluated in datasets from public PCa databases. The GSE52955, GSE76938 and GSE38240 datasets were downloaded and analyzed by GEO2R online software. Four CpG probes on the STEAP1 gene promoter were selected: cg15089950, Life 2021, 11, 1251 6 of 12 cg19317433, cg19532731 and cg24286372 located on −314 to −193 bp, −285 to −163 bp, −301 to −180 bp and −250 to −129 bp upstream of the transcription start, respectively. As shown in Figure 3, there were significant differences in the levels of CpG methylation between prostate tumor and normal tissue. In three datasets analyzed, the results reveal lower methylation levels of the STEAP1 gene promoter in neoplastic tissue than in normal tissue.

Correlation between STEAP1 Gene Promoter Demethylation and Its Expression in PCa Tissue from the TCGA Database
In an attempt to better understand whether DNA methylation status may have an impact on STEAP1 gene expression, the Prostate Adenocarcinoma dataset was analyzed from The Cancer Genome Atlas (TCGA, Cell 2015) [17], accessed through the cBioPortal. In this dataset, differential methylation levels between the unaltered (n = 286) and altered (n = 43) expression of the STEAP1 gene were observed (Figure 4a). Additionally, a negative correlation (Spearman coefficient of −0.42, and Pearson coefficient of −0.44) was observed between the methylation levels and the STEAP1 gene expression (Figure 4b).

Correlation between STEAP1 Gene Promoter Demethylation and Its Expression in PCa Tissue from the TCGA Database
In an attempt to better understand whether DNA methylation status may have an impact on STEAP1 gene expression, the Prostate Adenocarcinoma dataset was analyzed from The Cancer Genome Atlas (TCGA, Cell 2015) [17], accessed through the cBioPorta In this dataset, differential methylation levels between the unaltered (n = 286) and altered (n = 43) expression of the STEAP1 gene were observed (Figure 4a). Additionally, a nega tive correlation (Spearman coefficient of −0.42, and Pearson coefficient of −0.44) was ob served between the methylation levels and the STEAP1 gene expression (Figure 4b).

Effect of Epigenetic-Modulating Drugs in STEAP1 Gene Expression in Non-Neoplasti PNT1A Cells
In order to support that epigenetic modifications contribute to the regulation o STEAP1 expression, non-neoplastic PNT1A cells were used to evaluate the effect o DNMT and HDAC inhibitors (5-aza-dC and TSA, respectively) on STEAP1 mRNA ex pression by qPCR. As shown in Figure 5, treatment with the demethylation agent 5-aza dC induced a three-fold increase in STEAP1 mRNA levels when compared to the contro group (p < 0.01). Moreover, treatment with both epidrugs (5-aza-dC + TSA), which con tributes to demethylation and histone hyperacetylation, induced a strong increase (15-fold variation relative to control, p < 0.001) in STEAP1 mRNA levels.

Effect of Epigenetic-Modulating Drugs in STEAP1 Gene Expression in Non-Neoplastic PNT1A Cells
In order to support that epigenetic modifications contribute to the regulation of STEAP1 expression, non-neoplastic PNT1A cells were used to evaluate the effect of DNMT and HDAC inhibitors (5-aza-dC and TSA, respectively) on STEAP1 mRNA expression by qPCR. As shown in Figure 5, treatment with the demethylation agent 5-aza-dC induced a three-fold increase in STEAP1 mRNA levels when compared to the control group (p < 0.01). Moreover, treatment with both epidrugs (5-aza-dC + TSA), which contributes to demethylation and histone hyperacetylation, induced a strong increase (15-fold variation relative to control, p < 0.001) in STEAP1 mRNA levels.

Analysis of Co-Expression between STEAP1 mRNA Expression and HDACs
The TSA drug is a potent and specific inhibitor of HDAC classes I and II, which include ten isoforms (HDAC1-HDAC10) [22]. Considering that the results above suggest a negative association between histone deacetylation and STEAP1 mRNA expression, we intended to analyze the association between STEAP1 mRNA expression and HDAC mRNA expression, using the Prostate Adenocarcinoma (TCGA, Cell 2015) [17] dataset. Of the ten isoforms analyzed, STEAP1 mRNA expression showed a positive association with HDAC8, as highlighted in Table 3. On the other hand, STEAP1 mRNA expression was negatively associated with HDAC4,5,7 and 9 (Table 3). . Effect of treatment with 5-aza-dC and TSA (DNMT and HDAC inhibitors, respectively) on STEAP1 mRNA expression in PNT1A cells. Relative STEAP1 mRNA expression was determined by qPCR analysis after normalization with the GAPDH and β2M housekeeping genes. Results are expressed as fold-variation relative to the control group. Error bars indicate mean ± SEM (n = 6). ** p < 0.01 and *** p < 0.001 relative to control. $$$ p < 0.001 relative to 5-aza-dC.

Analysis of Co-Expression between STEAP1 mRNA Expression and HDACs
The TSA drug is a potent and specific inhibitor of HDAC classes I and II, which include ten isoforms (HDAC1-HDAC10) [22]. Considering that the results above suggest a negative association between histone deacetylation and STEAP1 mRNA expression, we intended to analyze the association between STEAP1 mRNA expression and HDAC mRNA expression, using the Prostate Adenocarcinoma (TCGA, Cell 2015) [17] dataset. Of the ten isoforms analyzed, STEAP1 mRNA expression showed a positive association with HDAC8, as highlighted in Table 3. On the other hand, STEAP1 mRNA expression was negatively associated with HDAC4,5,7 and 9 (Table 3). −0.064 0.245 Figure 5. Effect of treatment with 5-aza-dC and TSA (DNMT and HDAC inhibitors, respectively) on STEAP1 mRNA expression in PNT1A cells. Relative STEAP1 mRNA expression was determined by qPCR analysis after normalization with the GAPDH and β2M housekeeping genes. Results are expressed as fold-variation relative to the control group. Error bars indicate mean ± SEM (n = 6). ** p < 0.01 and *** p < 0.001 relative to control. $$$ p < 0.001 relative to 5-aza-dC.

Discussion
The methylation pattern observed in normal tissues undergoes relevant modifications in cancer, leading to changes in the regulation of the transcription of numerous genes [23]. Recent studies have shown that STEAP1 acts as an oncogene, showing that its overexpression in several human cancers contributes to tumor progression and aggressiveness through the inhibition of apoptosis and stimulation of cell proliferation, invasion and epithelial-mesenchymal transition [11][12][13][14][15]. Although the stability of the STEAP1 gene and protein is higher in LNCaP PCa cells than in PNT1A cells, contributing to STEAP1 overexpression, other mechanisms underlying its overexpression in cancer must be involved. As epigenetics has been pointed out as a major hallmark in cancer, affecting genes involved in all cellular pathways [24,25], our main goal was to assess whether epigenetic mechanisms are involved in the regulation of the STEAP1 gene expression in PCa, and if there are changes between normal and PCa cells. As a first approach, two different cell lines were chosen, neoplastic LNCaP and non-neoplastic PNT1A, which express high and low levels of STEAP1, respectively [16], to analyze the methylation status of the STEAP1 gene. Furthermore, five human PCa samples were also analyzed. The major CpG island located from position −338 of the promoter region to position +74 of the first exon of STEAP1 revealed differences in the methylation status. While in non-neoplastic PNT1A cells, the CpG dinucleotides near the transcription start site are methylated, in neoplastic LNCaP cells and PCa samples these are demethylated. This result suggests that demethylation of the STEAP1 gene promoter may lead to its overexpression in PCa. Hypomethylation/demethylation-dependent overexpression of several oncogenes has already been described in several cancer types. In PCa cells, Wingless-related MMTV integration site 5A, S100 calcium-binding protein P and cysteine-rich protein 1 were found to be hypomethylated [26]. hsa-miR-191 was also hypomethylated in 63% of hepatocellular carcinoma tissue samples, associated with its increased expression [27], and demethylation of the miR-128a promoter region drives the upregulation of miR-128a expression in the human T lymphocyte Jurkat cell line [28].
The possibility that demethylation of the STEAP1 gene promoter may be involved in its overexpression in PCa was also corroborated by a bioinformatic analysis, using public datasets from the GEO database and cBioPortal platform. This analysis showed that CpG dinucleotides in the STEAP1 gene of PCa samples have low levels of methylation, negatively correlated with STEAP1 mRNA expression. These results suggest, once more, that demethylation of the STEAP1 gene promoter may contribute to its overexpression in PCa. Our results are in line with other studies that also analyzed genes of the STEAP family, which showed that epigenetic alterations may be responsible for changes in gene expression. Using combined analysis of GEO and TCGA datasets, it was shown that STEAP3 is hypomethylated and consequently upregulated in glioblastoma, and may be used as a potential methylation-based prognostic biomarker. In addition, the authors suggested that STEAP3 is a potential target for glioblastoma treatment [29]. Another gene that was also reported to be deregulated due to epigenetic changes is the STEAP4 gene. Tamura et al. showed no CpG methylation in the STEAP4 promoter region in LNCaP cells, suggesting that demethylation may activate the expression of the STEAP4 gene in PCa cells [30]. On the other hand, a more recent study showed that STEAP4 was hypermethylated and downregulated in the hepatocellular carcinoma when compared to the non-tumor liver tissues [31]. This study also demonstrated that expression of STEAP4 was restored with a DNA methyltransferase (DNMTs) inhibitor and a histone deacetylase (HDACs) inhibitor, suggesting that aberrant DNA methylation suppressed the expression of the STEAP4 gene [31].
It is well established that there is a crosstalk between DNA methylation and histone modifications in the control of gene expression [32,33]. Aberrant CpG-island methylation by recruiting DNMTs and HDACs might be directly targeted by consequence of gene expression alterations [34]. Therefore, it was hypothesized that treatment with both DNMT and HDAC inhibitors might enhance STEAP1 expression in cells with reduced levels of STEAP1. Thus, the non-neoplastic PNT1A cell line was used to evaluate the effects of 5-aza-dC (DNMT inhibitor) plus TSA (HDAC inhibitor) on the expression of STEAP1.
The results indicate an increase in STEAP1 gene expression in response to treatment with DNMT and HDAC inhibitors, suggesting a synergistic effect of combined hypomethylation and histone hyperacetylation.
HDACs are enzymes that remove acetyl groups from the tails of histones, resulting in a more closed chromatin structure and repression of gene expression [22,35,36]. Class I of the HDAC family includes HDAC1, 2, 3 and 8; class II includes HDAC4, 5, 6, 7, 9 and 10 [36]. To the best of our knowledge, the role played by both HDAC class I and II enzymes in the regulation of STEAP1 gene expression is still unknown. So, we aimed to explore the association between HDAC and STEAP1 gene expression. Our results reveal a positive association of HDAC8 and a negative association of HDAC4,5,7 and 9 with STEAP1 mRNA expression. In fact, the HDAC family modulates several genes involved in cancer Life 2021, 11, 1251 10 of 12 development/progression via angiogenesis, cell adhesion, migration and invasion [22]. Some studies support our results, showing that increased expression of autotaxin in PCa cell lines is mediated by the downregulation of HDAC7 and HDAC3 [37]; additionally, HDAC5 is downregulated in human cancer, namely PCa [38], and decreased expression of HDAC4 increases the growth of PCa cells [39]. On the other hand, there are also studies showing an overexpression of HDAC in several types of cancer, suggesting the use of HDAC inhibitors as a promising class of compounds for cancer treatment [40][41][42][43]. Thus, the effects that can be triggered by these inhibitors on oncogenes should not be ignored and more studies are required to clarify their effectiveness in cancer treatment.
To summarize, this study showed that the STEAP1 gene is hypomethylated in PCa cells when compared to non-neoplastic cells, contributing to the overexpression of STEAP1 in PCa. Furthermore, our results suggest a putative involvement of HDCA4,5,7 and 9 in the regulation of STEAP1 in PCa. Considering the complexity of the mechanisms associated with HDAC, more studies are required to clarify their role in STEAP1 regulation, as well as to elucidate this association with PCa development and progression.