Abstract
HOX gene-encoded homeobox proteins control body patterning during embryonic development; the specific expression pattern of HOX genes may correspond to tissue identity. In this study, using RNAseq data of 1019 human cancer cell lines that originated from 24 different anatomic sites, we established HOX codes for various types of tissues. We applied these HOX codes to the transcriptomic profiles of prostate cancer (PCa) samples and found that the majority of prostate adenocarcinoma (AdPCa) samples sustained a prostate-specific HOX code whereas the majority of neuroendocrine prostate cancer (NEPCa) samples did not, which reflects the anaplastic nature of NEPCa. Also, our analysis showed that the NEPCa samples did not correlate well with the HOX codes of any other tissue types, indicating that NEPCa tumors lose their prostate identities but do not gain new tissue identities. Additionally, using immunohistochemical staining, we evaluated the prostatic expression of HOXB13, the most prominently changed HOX gene in NEPCa. We found that HOXB13 was expressed in both benign prostatic tissues and AdPCa but its expression was reduced or lost in NEPCa. Furthermore, we treated PCa cells with all trans retinoic acid (ATRA) and found that the reduced HOXB13 expression can be reverted. This suggests that ATRA is a potential therapeutic agent for the treatment of NEPCa tumors by reversing them to a more treatable AdPCa.