Hoxb13 is required for normal differentiation and secretory function of the ventral prostate

Effects of HOX family regulator-mediated modification patterns and immunity characteristics on tumor-associated cell type in endometrial cancer

Endometrial cancer (UCEC) is one of three major malignant tumors in women. The HOX gene regulates tumor development. However, the potential roles of HOX in the expression mechanism of multiple cell types and in the development and progression of tumor microenvironment (TME) cell infiltration in UCEC remain unknown. In this study, we utilized both the The Cancer Genome Atlas (TCGA) database and International Cancer Genome Consortium (ICGC) database to analyze transcriptome data of 529 patients with UCEC based on 39 HOX genes, combing clinical information, we discovered HOX gene were a pivotal factor in the development and progression of UCEC and in the formation of TME diversity and complexity. Here, a new scoring system was developed to quantify individual HOX patterns in UCEC. Our study found that patients in the low HOX score group had abundant anti-tumor immune cell infiltration, good tumor differentiation, and better prognoses. In contrast, a high HOX score was associated with blockade of immune checkpoints, which enhances the response to immunotherapy. The Real-Time quantitative PCR (RT-qPCR) and Immunohistochemistry (IHC) exhibited a higher expression of the HOX gene in the tumor patients. We revealed that the significant upregulation of the HOX gene in the epithelial cells can activate signaling pathway associated with tumour invasion and metastasis through single-cell RNA sequencing (scRNA-seq), such as nucleotide metabolic proce and so on. Finally, a risk prognostic model established by the positive relationship between HOX scores and cancer-associated fibroblasts (CAFs) can predict the prognosis of individual patients by scRNA-seq and transcriptome data sets. In sum, HOX gene may serve as a potential biomarker for the diagnosis and prediction of UCEC and to develop more effective therapeutic strategies.

Amino Terminal Acetylation of HOXB13 Regulates the DNA Damage Response in Prostate Cancer

Advanced localized prostate cancers (PC) recur despite chemotherapy, radiotherapy and/or androgen deprivation therapy. We recently reported HOXB13 lysine (K)13 acetylation as a gain-of-function modification that regulates interaction with the SWI/SNF chromatin remodeling complex and is critical for anti-androgen resistance. However, whether acetylated HOXB13 promotes PC cell survival following treatment with genotoxic agents is not known. Herein, we show that K13-acetylated HOXB13 is induced rapidly in PC cells in response to DNA damage induced by irradiation (IR). It colocalizes with the histone variant γH2AX at sites of double strand breaks (DSBs). Treatment of PCs with the Androgen Receptor (AR) antagonist Enzalutamide (ENZ) did not suppress DNA-damage-induced HOXB13 acetylation. In contrast, HOXB13 depletion or loss of acetylation overcame resistance of PC cells to ENZ and synergized with IR. HOXB13K13A mutants show diminished replication fork progression, impaired G2/M arrest with significant cell death following DNA damage. Mechanistically, we found that amino terminus regulates HOXB13 nuclear puncta formation that is essential for proper DNA damage response. Therefore, targeting HOXB13 acetylation with CBP/p300 inhibitors in combination with DNA damaging therapy may be an effective strategy to overcome anti-androgen resistance of PCs.

The Impact of Circulating Tumor Cell HOXB13 RNA Detection in Men with Metastatic Castration-Resistant Prostate Cancer (mCRPC) Treated with Abiraterone or Enzalutamide

PURPOSE

HOXB13 is an androgen receptor (AR) coregulator specifically expressed in cells of prostatic lineage. We sought to associate circulating tumor cell (CTC) HOXB13 expression with outcomes in men with mCRPC treated with abiraterone or enzalutamide.

EXPERIMENTAL DESIGN

We conducted a retrospective analysis of the multicenter prospective PROPHECY trial of mCRPC men (NCT02269982, n = 118) treated with abiraterone/enzalutamide. CTC detection and HOXB13 complementary DNA (cDNA) expression was measured using a modified Adnatest, grouping patients into 3 categories: CTC 0 (undetectable); CTC+ HOXB13 CTC low (<4 copies); or CTC+ HOXB13 CTC high. The HOXB13 threshold was determined by maximally selected rank statistics for prognostic associations with overall survival (OS) and progression-free survival (PFS).

RESULTS

We included 102 men with sufficient CTC HOXB13 cDNA, identifying 25%, 31%, and 44% of patients who were CTC 0, CTC+ HOXB13 low, and CTC+ HOXB13 high, respectively. Median OS were 25.7, 27.8, and 12.1 months whereas the median PFS were 9.0, 7.7, and 3.8 months, respectively. In subgroup analysis among men with CellSearch CTCs ≥5 copies/mL and adjusting for prior abi/enza treatment and Halabi clinical risk score, the multivariate HR for HOXB13 CTC detection was 2.39 (95% CI, 1.06-5.40) for OS and 2.78 (95% CI, 1.38-5.59) for PFS, respectively. Low HOXB13 CTC detection was associated with lower CTC PSA, PSMA, AR-FL, and AR-V7 detection, and more liver/lung metastases (41% vs. 25%).

CONCLUSIONS

Higher CTC HOXB13 expression is associated with AR-dependent biomarkers in CTCs and is adversely prognostic in the context of potent AR inhibition in men with mCRPC.

Characterization of HOXB13 expression patterns in localized and metastatic castration-resistant prostate cancer

HOXB13 is a key lineage homeobox transcription factor that plays a critical role in the differentiation of the prostate gland. Several studies have suggested that HOXB13 alterations may be involved in prostate cancer development and progression. Despite its potential biological relevance, little is known about the expression of HOXB13 across the disease spectrum of prostate cancer. To this end, we validated a HOXB13 antibody using genetic controls and investigated HOXB13 protein expression in murine and human developing prostates, localized prostate cancers, and metastatic castration-resistant prostate cancers. We observed that HOXB13 expression increases during later stages of murine prostate development. All localized prostate cancers showed HOXB13 protein expression. Interestingly, lower HOXB13 expression levels were observed in higher-grade tumors, although no significant association between HOXB13 expression and recurrence or disease-specific survival was found. In advanced metastatic prostate cancers, HOXB13 expression was retained in the majority of tumors. While we observed lower levels of HOXB13 protein and mRNA levels in tumors with evidence of lineage plasticity, 84% of androgen receptor-negative castration-resistant prostate cancers and neuroendocrine prostate cancers (NEPCs) retained detectable levels of HOXB13. Notably, the reduced expression observed in NEPCs was associated with a gain of HOXB13 gene body CpG methylation. In comparison to the commonly used prostate lineage marker NKX3.1, HOXB13 showed greater sensitivity in detecting advanced metastatic prostate cancers. Additionally, in a cohort of 837 patients, 383 with prostatic and 454 with non-prostatic tumors, we found that HOXB13 immunohistochemistry had a 97% sensitivity and 99% specificity for prostatic origin. Taken together, our studies provide valuable insight into the expression pattern of HOXB13 during prostate development and cancer progression. Furthermore, our findings support the utility of HOXB13 as a diagnostic biomarker for prostate cancer, particularly to confirm the prostatic origin of advanced metastatic castration-resistant tumors. © 2023 The Pathological Society of Great Britain and Ireland.

Hierarchical Phosphorylation of HOXB13 by mTOR Dictates Its Activity and Oncogenic Function in Prostate Cancer

Dysregulation of mTOR signaling plays a critical role in promoting prostate cancer growth. HOXB13, a homeodomain transcription factor, is known to influence the androgen response and prostate cancer development. Recently, HOXB13 was found to complex with mTOR on chromatin. However, the functional crosstalk between HOXB13 and mTOR remains elusive. We now report that mTOR directly interacts with and hierarchically phosphorylates HOXB13 at threonine 8 and 41 then serine 31 to promote its interaction with the E3 ligase SKP2 while enhancing its oncogenic properties. Expression of HOXB13 harboring phosphomimetic mutations at the mTOR-targeted sites stimulates prostate cancer cellular growth both in vitro and in murine xenografts. Transcriptional profiling studies revealed a phospho-HOXB13-dependent gene signature capable of robustly discriminating between normal prostate tissues, primary and metastatic prostate cancer samples. This work uncovers a previously unanticipated molecular cascade by which mTOR directly phosphorylates HOXB13 to dictate a specific gene program with oncogenic implications in prostate cancer.

IMPLICATIONS

Control of HOXB13 transcriptional activity via its direct phosphorylation by the mTOR kinase is a potential therapeutic avenue for the management of advanced prostate cancer.

YAP is required for prostate development, regeneration, and prostate stem cell function

Prostate development and regeneration depend on prostate stem cell function, the delicate balance of stem cell self-renewal and differentiation. However, mechanisms modulating prostate stem cell function remain poorly identified. Here, we explored the roles of Yes-associated protein 1 (YAP) in prostate stem cells, prostate development and regeneration. Using YAPfl/fl, CD133-CreER mice, we found that stem cell-specific YAP-deficient mice had compromised branching morphogenesis and epithelial differentiation, resulting in damaged prostate development. YAP inhibition also significantly affected the regeneration process of mice prostate, leading to impaired regenerated prostate. Furthermore, YAP ablation in prostate stem cells significantly reduced its self-renewal activity in vitro, and attenuated prostate regeneration of prostate grafts in vivo. Further analysis revealed a decrease in Notch and Hedgehog pathways expression in YAP inhibition cells, and treatment with exogenous Shh partially restored the self-renewal ability of prostate sphere cells. Taken together, our results revealed the roles of YAP in prostate stem cell function and prostate development and regeneration through regulation of the Notch and Hedgehog signaling pathways.

The prostate in women: an updated histological and immunohistochemical profile of the female periurethral glands and their relationship to an implanted midurethral sling

BACKGROUND

There is evidence of glandular tissue in the region of the anterior vaginal wall-female periurethral tissue (AVW-FPT) that has similar morphology and immunohistochemistry to the prostate in men. Surgical injury to this tissue has been suggested as a potential cause of sexual dysfunction following midurethral sling (MUS) procedures. However, the anatomy and embryology of these glands have not been fully resolved. This has led to difficulties in classifying this tissue as a prostate and defining its clinical significance related to MUS procedures.

AIM

To describe the histological and immunohistochemical characteristics of the female periurethral glands using markers of prostate tissue and innervation and to examine their anatomical relationships to an implanted MUS.

METHODS

Using gross and fine dissection, the AVW-FPT was dissected from 9 cadavers. Prior to dissection, 2 cadavers underwent simulation of the MUS procedure by a urogynecologist. Samples were paraffin embedded and serially sectioned. Immunohistochemistry was performed using markers of prostate tissue and innervation.

OUTCOMES

Redundant immunohistochemical localization of markers for prostatic tissue and innervation of the glandular tissue of the AVW-FPT, including the region of MUS implantation.

RESULTS

Female periurethral glands were immunoreactive for markers of male prostatic tissue, including prostate-specific antigen, androgen receptor, HOXB13, and NKX3.1. Markers of innervation (protein gene product 9.5, choline acetyl transferase, and vasoactive intestinal polypeptide) also localized to certain regions of the glandular tissue and associated blood supply. Surgical simulation of the MUS procedure demonstrated that some periurethral glands are located in close proximity to an implanted sling.

CLINICAL TRANSLATION

The AVW-FPT contains glandular tissue in the surgical field of MUS implantation. Iatrogenic damage to the female periurethral glands and the associated innervation during surgery could explain the negative impacts on sexual dysfunction reported following MUS procedures.

STRENGTHS AND LIMITATIONS

This is the first study to characterize the female periurethral glands using markers of prostatic tissue in concert with markers of general and autonomic innervation and characterize their anatomical relationships within the surgical field of MUS implantation. The small sample size is a limitation of this study.

CONCLUSION

We provide further evidence that the AVW-FPT contains innervated glands that are phenotypically similar to the male prostate and may share a common embryonic origin. The microscopic and immunohistochemical features of the periurethral glands may be indicative of their functional capacity in sexual responses. The location of these glands in the surgical field of MUS procedures underscores the clinical significance of this tissue.

HOXB13 facilitates hepatocellular carcinoma progression by activating AKT/mTOR signaling pathway

INTRODUCTION AND OBJECTIVES

Hepatocellular carcinoma (HCC) is one of the sixth most common malignancies worldwide and is accompanied by high mortality. Homeobox B13 (HOXB13) has been shown to be involved in the development of various cancers. This study aimed to investigate the role of HOXB13 in HCC progression.

MATERIALS AND METHODS

The expression of HOXB13 in HCC tumor tissues was analyzed using qRT-PCR and immunohistochemical staining . After overexpression or downregulation of HOXB13 in HCC cell lines, cell proliferation was detected by CCK8 assay and Ki67 staining and cell invasion ability were tested by transwell assay. Western blot assay was applied to analyze the effect of HOXB13 on related signaling pathways. In addition, the role of HOXB13 on HCC in vivo was explored using a HCC mouse model. IF and WB were performed to detect cell proliferation, apoptosis and related protein expression in mice tumor tissues.

RESULTS

The results showed that the expression of HOXB13 was significantly increased in HCC tissues compared with adjacent tissues and positively correlated with the tumor stage and survival of HCC patients. Overexpression of HOXB13 promoted the proliferation and invasion of HCC cells and up-regulated the protein expression of AKT, mTOR and MMP2. In contrast, the downregulation of HOXB13 resulted in the opposite results. In vivo experiments, HOXB13 significantly promoted tumor growth in mice bearing HCC by promoting cell proliferation and inhibiting cell apoptosis.

CONCLUSIONS

This study suggested that HOXB13 can facilitate HCC progression by activation of the AKT/mTOR signaling pathway. HOXB13 may be a novel target for HCC therapy.

A transcriptional network of cell cycle dysregulation in noninvasive papillary urothelial carcinoma

Human cancers display a restricted set of expression profiles, despite diverse mutational drivers. This has led to the hypothesis that select sets of transcription factors act on similar target genes as an integrated network, buffering a tumor's transcriptional state. Noninvasive papillary urothelial carcinoma (NIPUC) with higher cell cycle activity has higher risk of recurrence and progression. In this paper, we describe a transcriptional network of cell cycle dysregulation in NIPUC, which was delineated using the ARACNe algorithm applied to expression data from a new cohort (n = 81, RNA sequencing), and two previously published cohorts. The transcriptional network comprised 121 transcription factors, including the pluripotency factors SOX2 and SALL4, the sex hormone binding receptors ESR1 and PGR, and multiple homeobox factors. Of these 121 transcription factors, 65 and 56 were more active in tumors with greater and less cell cycle activity, respectively. When clustered by activity of these transcription factors, tumors divided into High Cell Cycle versus Low Cell Cycle groups. Tumors in the High Cell Cycle group demonstrated greater mutational burden and copy number instability. A putative mutational driver of cell cycle dysregulation, such as homozygous loss of CDKN2A, was found in only 50% of High Cell Cycle NIPUC, suggesting a prominent role of transcription factor activity in driving cell cycle dysregulation. Activity of the 121 transcription factors strongly associated with expression of EZH2 and other members of the PRC2 complex, suggesting regulation by this complex influences expression of the transcription factors in this network. Activity of transcription factors in this network also associated with signatures of pluripotency and epithelial-to-mesenchymal transition (EMT), suggesting they play a role in driving evolution to invasive carcinoma. Consistent with this, these transcription factors differed in activity between NIPUC and invasive urothelial carcinoma.

Genetic determinants of chromatin reveal prostate cancer risk mediated by context-dependent gene regulation

Many genetic variants affect disease risk by altering context-dependent gene regulation. Such variants are difficult to study mechanistically using current methods that link genetic variation to steady-state gene expression levels, such as expression quantitative trait loci (eQTLs). To address this challenge, we developed the cistrome-wide association study (CWAS), a framework for identifying genotypic and allele-specific effects on chromatin that are also associated with disease. In prostate cancer, CWAS identified regulatory elements and androgen receptor-binding sites that explained the association at 52 of 98 known prostate cancer risk loci and discovered 17 additional risk loci. CWAS implicated key developmental transcription factors in prostate cancer risk that are overlooked by eQTL-based approaches due to context-dependent gene regulation. We experimentally validated associations and demonstrated the extensibility of CWAS to additional epigenomic datasets and phenotypes, including response to prostate cancer treatment. CWAS is a powerful and biologically interpretable paradigm for studying variants that influence traits by affecting transcriptional regulation.

HOXB13 suppresses de novo lipogenesis through HDAC3-mediated epigenetic reprogramming in prostate cancer

HOXB13, a homeodomain transcription factor, critically regulates androgen receptor (AR) activities and androgen-dependent prostate cancer (PCa) growth. However, its functions in AR-independent contexts remain elusive. Here we report HOXB13 interaction with histone deacetylase HDAC3, which is disrupted by the HOXB13 G84E mutation that has been associated with early-onset PCa. Independently of AR, HOXB13 recruits HDAC3 to lipogenic enhancers to catalyze histone deacetylation and suppress lipogenic regulators such as fatty acid synthase. Analysis of human tissues reveals that the HOXB13 gene is hypermethylated and downregulated in approximately 30% of metastatic castration-resistant PCa. HOXB13 loss or G84E mutation leads to lipid accumulation in PCa cells, thereby promoting cell motility and xenograft tumor metastasis, which is mitigated by pharmaceutical inhibition of fatty acid synthase. In summary, we present evidence that HOXB13 recruits HDAC3 to suppress de novo lipogenesis and inhibit tumor metastasis and that lipogenic pathway inhibitors may be useful to treat HOXB13-low PCa.

A Review of Prostate Organogenesis and a Role for iPSC-Derived Prostate Organoids to Study Prostate Development and Disease

The prostate is vulnerable to two major age-associated diseases, cancer and benign enlargement, which account for significant morbidity and mortality for men across the globe. Prostate cancer is the most common cancer reported in men, with over 1.2 million new cases diagnosed and 350,000 deaths recorded annually worldwide. Benign prostatic hyperplasia (BPH), characterised by the continuous enlargement of the adult prostate, symptomatically afflicts around 50% of men worldwide. A better understanding of the biological processes underpinning these diseases is needed to generate new treatment approaches. Developmental studies of the prostate have shed some light on the processes essential for prostate organogenesis, with many of these up- or downregulated genes expressions also observed in prostate cancer and/or BPH progression. These insights into human disease have been inferred through comparative biological studies relying primarily on rodent models. However, directly observing mechanisms of human prostate development has been more challenging due to limitations in accessing human foetal material. Induced pluripotent stem cells (iPSCs) could provide a suitable alternative as they can mimic embryonic cells, and iPSC-derived prostate organoids present a significant opportunity to study early human prostate developmental processes. In this review, we discuss the current understanding of prostate development and its relevance to prostate-associated diseases. Additionally, we detail the potential of iPSC-derived prostate organoids for studying human prostate development and disease.

Genetics of prostate cancer and its utility in treatment and screening

Prostate cancer heritability is attributed to a combination of rare, moderate to highly penetrant genetic variants as well as commonly occurring variants conferring modest risks [single nucleotide polymorphisms (SNPs)]. Some of the former type of variants (e.g., BRCA2 mutations) predispose particularly to aggressive prostate cancer and confer poorer prognoses compared to men who do not carry mutations. Molecularly targeted treatments such as PARP inhibitors have improved outcomes in men carrying somatic and/or germline DNA repair gene mutations. Ongoing clinical trials are exploring other molecular targeted approaches based on prostate cancer somatic alterations. Genome wide association studies have identified >250 loci that associate with prostate cancer risk. Multi-ancestry analyses have identified shared as well as population specific risk SNPs. Prostate cancer risk SNPs can be used to estimate a polygenic risk score (PRS) to determine an individual's genetic risk of prostate cancer. The odds ratio of prostate cancer development in men whose PRS lies in the top 1% of the risk profile ranges from 9 to 11. Ongoing studies are investigating the utility of a prostate cancer PRS to target population screening to those at highest risk. With the advent of personalized medicine and development of DNA sequencing technologies, access to clinical genetic testing is increasing, and oncology guidelines from bodies such as NCCN and ESMO have been updated to provide criteria for germline testing of "at risk" healthy men as well as those with prostate cancer. Both germline and somatic prostate cancer research have significantly evolved in the past decade and will lead to further development of precision medicine approaches to prostate cancer treatment as well as potentially developing precision population screening models.

The HOXB13 variant X285K is associated with clinical significance and early age at diagnosis in African American prostate cancer patients

BACKGROUND

Recently, a novel HOXB13 variant (X285K) was observed in men of African descent with prostate cancer (PCa) in Martinique. Little is known about this or other variants in HOXB13 which may play a role in PCa susceptibility in African-American (AA) men.

METHODS

We sequenced HOXB13 in an AA population of 1048 men undergoing surgical treatment for PCa at Johns Hopkins Hospital.

RESULTS

Seven non-synonymous germline variants were observed in the patient population. While six of these variants were seen only once, X285K was found in eight patients. In a case-case analysis, we find that carriers of this latter variant are at increased risk of clinically significant PCa (1.2% carrier rate in Gleason Score ≥7 PCa vs. 0% in Gleason Score <7 PCa, odds ratio, OR = inf; 95% Confidence Interval, 95%CI:1.05-inf, P = 0.028), as well as PCa with early age at diagnosis (2.4% carrier rate in patients <50 year vs. 0.5% carrier rate in patients ≥50 year, OR = 5.25, 95% CI:1.00-28.52, P = 0.03).

CONCLUSIONS

While this variant is rare in the AA population (~0.2% MAF), its ancestry-specific occurrence and apparent preferential association with risk for the more aggressive disease at an early age emphasizes its translational potential as an important, novel PCa susceptibility marker in the high-risk AA population.

Stromal androgen and hedgehog signaling regulates stem cell niches in pubertal prostate development

Stromal androgen-receptor (AR) action is essential for prostate development, morphogenesis and regeneration. However, mechanisms underlying how stromal AR maintains the cell niche in support of pubertal prostatic epithelial growth are unknown. Here, using advanced mouse genetic tools, we demonstrate that selective deletion of stromal AR expression in prepubescent Shh-responsive Gli1-expressing cells significantly impedes pubertal prostate epithelial growth and development. Single-cell transcriptomic analyses showed that AR loss in these prepubescent Gli1-expressing cells dysregulates androgen signaling-initiated stromal-epithelial paracrine interactions, leading to growth retardation of pubertal prostate epithelia and significant development defects. Specifically, AR loss elevates Shh-signaling activation in both prostatic stromal and adjacent epithelial cells, directly inhibiting prostatic epithelial growth. Single-cell trajectory analyses further identified aberrant differentiation fates of prostatic epithelial cells directly altered by stromal AR deletion. In vivo recombination of AR-deficient stromal Gli1-lineage cells with wild-type prostatic epithelial cells failed to develop normal prostatic epithelia. These data demonstrate previously unidentified mechanisms underlying how stromal AR-signaling facilitates Shh-mediated cell niches in pubertal prostatic epithelial growth and development.

Developmental expression patterns of toolkit genes in male accessory gland of Drosophila parallels those of mammalian prostate

Conservation of genetic toolkits in disparate phyla may help reveal commonalities in organ designs transcending their extreme anatomical disparities. A male accessory sexual organ in mammals, the prostate, for instance, is anatomically disparate from its analogous, phylogenetically distant counterpart – the male accessory gland (MAG) – in insects like Drosophila. It has not been ascertained if the anatomically disparate Drosophila MAG shares developmental parallels with those of the mammalian prostate. Here we show that the development of Drosophila mesoderm-derived MAG entails recruitment of similar genetic toolkits of tubular organs like that seen in endoderm-derived mammalian prostate. For instance, like mammalian prostate, Drosophila MAG morphogenesis is marked by recruitment of fibroblast growth factor receptor (FGFR) – a signalling pathway often seen recruited for tubulogenesis – starting early during its adepithelial genesis. A specialisation of the individual domains of the developing MAG tube, on the other hand, is marked by the expression of a posterior Hox gene transcription factor, Abd-B, while Hh-Dpp signalling marks its growth. Drosophila MAG, therefore, reveals the developmental design of a unitary bud-derived tube that appears to have been co-opted for the development of male accessory sexual organs across distant phylogeny and embryonic lineages.

This article has an associated First Person interview with the first author of the paper.

Summary: We show genetic toolkit conservation between Drosophila MAG and mammalian prostate may suggest a common modular developmental design.

Single-cell analysis of mouse and human prostate reveals novel fibroblasts with specialized distribution and microenvironment interactions

Stromal-epithelial interactions are critical to the morphogenesis, differentiation, and homeostasis of the prostate, but the molecular identity and anatomy of discrete stromal cell types is poorly understood. Using single-cell RNA sequencing, we identified and validated the in situ localization of three smooth muscle subtypes (prostate smooth muscle, pericytes, and vascular smooth muscle) and two novel fibroblast subtypes in human prostate. Peri-epithelial fibroblasts (APOD+) wrap around epithelial structures, whereas interstitial fibroblasts (C7+) are interspersed in extracellular matrix. In contrast, the mouse displayed three fibroblast subtypes with distinct proximal-distal and lobe-specific distribution patterns. Statistical analysis of mouse and human fibroblasts showed transcriptional correlation between mouse prostate (C3+) and urethral (Lgr5+) fibroblasts and the human interstitial fibroblast subtype. Both urethral fibroblasts (Lgr5+) and ductal fibroblasts (Wnt2+) in the mouse contribute to a proximal Wnt/Tgfb signaling niche that is absent in human prostate. Instead, human peri-epithelial fibroblasts express secreted WNT inhibitors SFRPs and DKK1, which could serve as a buffer against stromal WNT ligands by creating a localized signaling niche around individual prostate glands. We also identified proximal-distal fibroblast density differences in human prostate that could amplify stromal signaling around proximal prostate ducts. In human benign prostatic hyperplasia, fibroblast subtypes upregulate critical immunoregulatory pathways and show distinct distributions in stromal and glandular phenotypes. A detailed taxonomy of leukocytes in benign prostatic hyperplasia reveals an influx of myeloid dendritic cells, T cells and B cells, resembling a mucosal inflammatory disorder. A receptor-ligand interaction analysis of all cell types revealed a central role for fibroblasts in growth factor, morphogen, and chemokine signaling to endothelia, epithelia, and leukocytes. These data are foundational to the development of new therapeutic targets in benign prostatic hyperplasia. © 2021 The Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Of Necks, Trunks and Tails: Axial Skeletal Diversity among Vertebrates

The axial skeleton of all vertebrates is composed of individual units known as vertebrae. Each vertebra has individual anatomical attributes, yet they can be classified in five different groups, namely cervical, thoracic, lumbar, sacral and caudal, according to shared characteristics and their association with specific body areas. Variations in vertebral number, size, morphological features and their distribution amongst the different regions of the vertebral column are a major source of the anatomical diversity observed among vertebrates. In this review I will discuss the impact of those variations on the anatomy of different vertebrate species and provide insights into the genetic origin of some remarkable morphological traits that often serve to classify phylogenetic branches or individual species, like the long trunks of snakes or the long necks of giraffes.

Identification of HOXA1 as a Novel Biomarker in Prognosis of Head and Neck Squamous Cell Carcinoma

Hox genes, a highly conserved homolog in most animals, play vital functions in cell development and organ formation. In recent years, researchers have discovered that it can act as a tumor regulator, and its members can participate in tumorigenesis by regulating receptor signaling, cell differentiation, apoptosis, migration, EMT, and angiogenesis. Hox genes and which major members play a vital role in the progress of head and neck squamous cell carcinoma (HNSCC) is still unclear. After analyzing the expression differences and prognostic value of all Hox genes through the TCGA-HNSC database, we use histochemistry stains in 52 pairs of HNSCC slices to verify the expression level of the key member-HOXA1. In correlation analysis, we found that high HOXA1 expression is related to poor pathological grade (p = 0.0077), advanced T stage (p = 0.021) and perineural invasion (PNI) (p = 0.0019). Furthermore, we used Cox univariate and multivariate regression analysis to confirm the independent predictive power of HOXA1 expression. To explore the underlying mechanisms behind HOXA1, we ran GSVA and GSEA and found fourteen mutual signaling pathways, including neuroprotein secretion and transport, tumor-associated signaling pathways, cell adhere junction and metabolic reprogramming. Finally, we found that the high expression of HOXA1 is significantly related to the decrease of CD8+ T cell infiltration and the decline of DNA methylation level. Our findings demonstrated that HOXA1, as a notable member of the HOX family, maybe an independent prognostic indicator in HNSCC.

HOXD13 suppresses prostate cancer metastasis and BMP4-induced epithelial-mesenchymal transition by inhibiting SMAD1

The HOX genes are a group of highly conserved Homeobox-containing genes that control the body plan organization during development. However, their contributions to tumorigenesis and tumor progression remain uncertain and controversial. Here we provided evidence of tumor-suppressive activity of HOXD13 in prostate cancer. HOXD13 depletion contributes to more aggressiveness of prostate cancer cells in vitro and in vivo. These effects were corroborated in a metastatic mice model, where we observed more bone metastatic lesions formed by prostate cancer cells with HOXD13 ablation. Mechanistically, HOXD13 prevents BMP4-induced epithelial-mesenchymal transition (EMT) by inhibiting mothers against decapentaplegic homolog 1 (SMAD1) transcription. Both bioinformation and our tissue microarray cohort data show that HOXD13 expression inversely correlated in advanced prostate cancer patient specimens. Our findings establish HOXD13 as a negative regulator of prostate cancer progression and metastasis by preventing BMP4/SMAD1 signaling, and potentially suggest new strategies for targeting metastatic prostate cancer.

Role of Runx2 in prostate development and stem cell function

BACKGROUND

RUNX2, a critical transcription factor in bone development, is also expressed in prostate and breast where it has been linked to cancer progression and cancer stem cells. However, its role in normal prostate biology has not been previously examined.

METHODS

Selective growth of murine prostate epithelium under non-adherent conditions was used to enrich for stem cells. Expression of runt domain transcription factors, stem cell and prostate marker messenger RNAs (mRNAs) was determined by quantitative reverse transcription polymerase chain reaction. Effects of Runx2 loss and gain-of-function on prostate epithelial cells were assessed using cells isolated from Runx2loxp/loxp mice transduced with Adeno-Cre or by Adeno-Runx2 transduction of wild type cells. Cellular distribution of RUNX2 and prostate-associated proteins was assessed using immunofluorescence microscopy. In vivo Runx2 knock out was achieved by tamoxifen treatment of Nkx3.1CreERT; Runx2loxp/loxp mice.

RESULTS

Prostate epithelium-derived spheroids, which are enriched in stem cells, were shown to contain elevated levels of Runx2 mRNA. Spheroid formation required Runx2 since adenovirus-Cre mediated knockout of Runx2 in prostatic epithelial cells from Runx2loxp/loxp mice severely reduced spheroid formation and stem cell markers while Runx2 overexpression was stimulatory. In vivo, Runx2 was detected during early prostate development (E16.5) and in adult mice where it was present in basal and luminal cells of ventral and anterior lobes. Prostate-selective deletion of Runx2 in tamoxifen-treated Nkx3.1CreERT; Runx2loxp/loxp mice severely inhibited growth and maturation of tubules in the anterior prostate and reduced expression of stem cell markers and prostate-associated genes.

CONCLUSION

This study demonstrates an important role for Runx2 in prostate development that may be explained by actions in prostate epithelial stem cells.

Neuroendocrine prostate cancer has distinctive, non-prostatic HOX code that is represented by the loss of HOXB13 expression

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.

Developmental estrogenization: Prostate gland reprogramming leads to increased disease risk with aging

While estrogens are involved in normal prostate morphogenesis and function, inappropriate early-life estrogenic exposures, either in type, dose or timing, can reprogram the prostate gland and lead to increased disease risk with aging. This process is referred to as estrogen imprinting or developmental estrogenization of the prostate gland. The present review discusses published and new evidence for prostatic developmental estrogenization that includes extensive research in rodent models combined with epidemiology findings that together have helped to uncover the architectural and molecular underpinnings that promote this phenotype. Complex interactions between steroid receptors, developmental morphoregulatory factors, epigenetic machinery and stem-progenitor cell targets coalesce to hard wire structural, cellular and epigenomic reorganization of the tissue which retains a life-long memory of early-life estrogens, ultimately predisposing the gland to prostatitis, hyperplasia and carcinogenesis with aging.

ACK1-AR and AR-HOXB13 signaling axes: epigenetic regulation of lethal prostate cancers

The androgen receptor (AR) is a critical transcription factor in prostate cancer (PC) pathogenesis. Its activity in malignant cells is dependent on interactions with a diverse set of co-regulators. These interactions fluctuate depending on androgen availability. For example, the androgen depletion increases the dependence of castration-resistant PCs (CRPCs) on the ACK1 and HOXB13 cell survival pathways. Activated ACK1, an oncogenic tyrosine kinase, phosphorylates cytosolic and nuclear proteins, thereby avoiding the inhibitory growth consequences of androgen depletion. Notably, ACK1-mediated phosphorylation of histone H4, which leads to epigenetic upregulation of AR expression, has emerged as a critical mechanism of CRPC resistance to anti-androgens. This resistance can be targeted using the ACK1-selective small-molecule kinase inhibitor (R)- 9b. CRPCs also deploy the bromodomain and extra-terminal domain protein BRD4 to epigenetically increase HOXB13 gene expression, which in turn activates the MYC target genes AURKA/AURKB. HOXB13 also facilitates ligand-independent recruitment of the AR splice variant AR-V7 to chromatin, compensating for the loss of the chromatin remodeling protein, CHD1, and restricting expression of the mitosis control gene HSPB8. These studies highlight the crosstalk between AR-ACK1 and AR-HOXB13 pathways as key mediators of CRPC recurrence.

Prostate cancer reactivates developmental epigenomic programs during metastatic progression

Epigenetic processes govern prostate cancer (PCa) biology, as evidenced by the dependency of PCa cells on the androgen receptor (AR), a prostate master transcription factor. We generated 268 epigenomic datasets spanning two state transitions-from normal prostate epithelium to localized PCa to metastases-in specimens derived from human tissue. We discovered that reprogrammed AR sites in metastatic PCa are not created de novo; rather, they are prepopulated by the transcription factors FOXA1 and HOXB13 in normal prostate epithelium. Reprogrammed regulatory elements commissioned in metastatic disease hijack latent developmental programs, accessing sites that are implicated in prostate organogenesis. Analysis of reactivated regulatory elements enabled the identification and functional validation of previously unknown metastasis-specific enhancers at HOXB13, FOXA1 and NKX3-1. Finally, we observed that prostate lineage-specific regulatory elements were strongly associated with PCa risk heritability and somatic mutation density. Examining prostate biology through an epigenomic lens is fundamental for understanding the mechanisms underlying tumor progression.

MEIS-mediated suppression of human prostate cancer growth and metastasis through HOXB13-dependent regulation of proteoglycans

The molecular roles of HOX transcriptional activity in human prostate epithelial cells remain unclear, impeding the implementation of new treatment strategies for cancer prevention and therapy. MEIS proteins are transcription factors that bind and direct HOX protein activity. MEIS proteins are putative tumor suppressors that are frequently silenced in aggressive forms of prostate cancer. Here we show that MEIS1 expression is sufficient to decrease proliferation and metastasis of prostate cancer cells in vitro and in vivo murine xenograft models. HOXB13 deletion demonstrates that the tumor-suppressive activity of MEIS1 is dependent on HOXB13. Integration of ChIP-seq and RNA-seq data revealed direct and HOXB13-dependent regulation of proteoglycans including decorin (DCN) as a mechanism of MEIS1-driven tumor suppression. These results define and underscore the importance of MEIS1-HOXB13 transcriptional regulation in suppressing prostate cancer progression and provide a mechanistic framework for the investigation of HOXB13 mutants and oncogenic cofactors when MEIS1/2 are silenced.

Shaping Chromatin States in Prostate Cancer by Pioneer Transcription Factors

The androgen receptor (AR) is a critical therapeutic target in prostate cancer that responds to antagonists in primary disease, but inevitably becomes reactivated, signaling onset of the lethal castration-resistant prostate cancer (CRPC) stage. Epigenomic investigation of the chromatin environment and interacting partners required for AR transcriptional activity has uncovered three pioneer factors that open up chromatin and facilitate AR-driven transcriptional programs. FOXA1, HOXB13, and GATA2 are required for normal AR transcription in prostate epithelial development and for oncogenic AR transcription during prostate carcinogenesis. AR signaling is dependent upon these three pioneer factors both before and after the clinical transition from treatable androgen-dependent disease to untreatable CRPC. Agents targeting their respective DNA binding or downstream chromatin-remodeling events have shown promise in preclinical studies of CRPC. AR-independent functions of FOXA1, HOXB13, and GATA2 are emerging as well. While all three pioneer factors exert effects that promote carcinogenesis, some of their functions may inhibit certain stages of prostate cancer progression. In all, these pioneer factors represent some of the most promising potential therapeutic targets to emerge thus far from the study of the prostate cancer epigenome.

HOXB13 expression is correlated with hepatic inflammatory activity of patients with hepatic fibrosis

Liver fibrosis is a common pathological process of chronic hepatic injury, preceded by the chronic inflammation. The homeobox B13 (HOXB13) gene, a member of HOX family, plays diverse biological roles in embryonic development, carcinogenesis, and many inflammatory diseases. However, the expression of HOXB13 in chronic liver diseases including hepatic fibrosis remains to be defined. In present study, 55 patients with hepatic fibrosis, 15 patients of hepatocellular carcinoma, and 17 healthy controls were enrolled in this study. Pathological specimens were collected through liver biopsy or surgical resection. The degree of hepatic inflammation (G0-G4) and fibrosis (S0-S4) of hepatic fibrosis was scored based on the modified histology activity index. Intrahepatic HOXB13 expression was analyzed using immunohistochemistry analysis. Compared with healthy subjects, both patients with hepatic fibrosis and patients with hepatocellular carcinoma exhibited significant accumulations of HOXB13+ cells in the liver (p < 0.05). Additionally, the number of HOXB13+ cell was significantly elevated along with the increment of hepatic inflammatory activities, but not fibrosis stages, among these liver fibrosis samples (p < 0.01). Furthermore, the quantity of HOXB13+ cells were also positively correlated with hepatic enzymes, alanine transaminase (r = 0.299, p = 0.041) and aspartate aminotransferase (r = 0.317, p = 0.013) in our cohort of hepatic fibrosis. In conclusion, our study identified a strong hepatic expression of HOXB13 among patients with hepatic fibrosis, which strongly associated with the degree of hepatic inflammatory activity for patients with hepatic fibrosis, suggesting an important role of HOXB13 during the pathogenesis of liver fibrogenesis.

The vertebrate tail: a gene playground for evolution

The tail of all vertebrates, regardless of size and anatomical detail, derive from a post-anal extension of the embryo known as the tail bud. Formation, growth and differentiation of this structure are closely associated with the activity of a group of cells that derive from the axial progenitors that build the spinal cord and the muscle-skeletal case of the trunk. Gdf11 activity switches the development of these progenitors from a trunk to a tail bud mode by changing the regulatory network that controls their growth and differentiation potential. Recent work in the mouse indicates that the tail bud regulatory network relies on the interconnected activities of the Lin28/let-7 axis and the Hox13 genes. As this network is likely to be conserved in other mammals, it is possible that the final length and anatomical composition of the adult tail result from the balance between the progenitor-promoting and -repressing activities provided by those genes. This balance might also determine the functional characteristics of the adult tail. Particularly relevant is its regeneration potential, intimately linked to the spinal cord. In mammals, known for their complete inability to regenerate the tail, the spinal cord is removed from the embryonic tail at late stages of development through a Hox13-dependent mechanism. In contrast, the tail of salamanders and lizards keep a functional spinal cord that actively guides the tail's regeneration process. I will argue that the distinct molecular networks controlling tail bud development provided a collection of readily accessible gene networks that were co-opted and combined during evolution either to end the active life of those progenitors or to make them generate the wide diversity of tail shapes and sizes observed among vertebrates.

DNA-PK, Nuclear mTOR, and the Androgen Pathway in Prostate Cancer

Androgen and its receptor (AR) are major drivers of prostate cancer (PCa), a leading cause of mortality in aging men. Thus, understanding the numerous mechanisms by which AR can promote the growth and proliferation of PCa cells and enable their escape from hormone-dependent therapies, eventually leading to metastasis and death of the patient, is essential to discover alternative therapeutic approaches. Recently, two structurally related members of the phosphatidylinositol 3-kinase-like protein kinase (PIKK) family, DNA-dependent protein kinase (DNA-PK) and mammalian target of rapamycin (mTOR), were shown to have a direct role in modulating AR activity on chromatin of PCa cells. In this review, the common features of DNA-PK and mTOR and the similarities in their noncanonical roles as transcription coregulators of the AR are highlighted. An outlook on how these findings could be translated into new approaches to manage and treat PCa is provided.

Genome-wide analysis of HOXC4 and HOXC6 regulated genes and binding sites in prostate cancer cells

Aberrant expression of HOXC6 and HOXC4 is commonly detected in prostate cancer. The high expression of these transcription factors is associated with aggressive prostate cancer and can predict cancer recurrence after treatment. Thus, HOXC4 and HOXC6 are clinically relevant biomarkers of aggressive prostate cancer. However, the molecular mechanisms by which these HOXC genes contribute to prostate cancer is not yet understood. To begin to address the role of HOXC4 and HOXC6 in prostate cancer, we performed RNA-seq analyses before and after siRNA-mediated knockdown of HOXC4 and/or HOXC6 and also performed ChIP-seq to identify genomic binding sites for both of these transcription factors. Our studies demonstrate that HOXC4 and HOXC6 co-localize with HOXB13, FOXA1 and AR, three transcription factors previously shown to contribute to the development of prostate cancer. We suggest that the aberrantly upregulated HOXC4 and HOXC6 proteins may compete with HOXB13 for binding sites, thus altering the prostate transcriptome. This competition model may be applicable to many different human cancers that display increased expression of a HOX transcription factor.

Cistrome Partitioning Reveals Convergence of Somatic Mutations and Risk Variants on Master Transcription Regulators in Primary Prostate Tumors

Thousands of noncoding somatic single-nucleotide variants (SNVs) of unknown function are reported in tumors. Partitioning the genome according to cistromes reveals the enrichment of somatic SNVs in prostate tumors as opposed to adjacent normal tissue cistromes of master transcription regulators, including AR, FOXA1, and HOXB13. This parallels enrichment of prostate cancer genetic predispositions over these transcription regulators' tumor cistromes, exemplified at the 8q24 locus harboring both risk variants and somatic SNVs in cis-regulatory elements upregulating MYC expression. However, Massively Parallel Reporter Assays reveal that few SNVs can alter the transactivation potential of individual cis-regulatory elements. Instead, similar to inherited risk variants, SNVs accumulate in cistromes of master transcription regulators required for prostate cancer development.

The Homeobox gene, HOXB13, Regulates a Mitotic Protein-Kinase Interaction Network in Metastatic Prostate Cancers

HOXB13, a homeodomain transcription factor, is linked to recurrence following radical prostatectomy. While HOXB13 regulates Androgen Receptor (AR) functions in a context dependent manner, its critical effectors in prostate cancer (PC) metastasis remain largely unknown. To identify HOXB13 transcriptional targets in metastatic PCs, we performed integrative bioinformatics analysis of differentially expressed genes (DEGs) in the proximity of the human prostate tumor-specific AR binding sites. Unsupervised Principal Component Analysis (PCA) led to a focused core HOXB13 target gene-set referred to as HOTPAM9 (HOXB13 Targets separating Primary And Metastatic PCs). HOTPAM9 comprised 7 mitotic kinase genes overexpressed in metastatic PCs, TRPM8, and the heat shock protein HSPB8, whose levels were significantly lower in metastatic PCs compared to the primary disease. The expression of a two-gene set, CIT and HSPB8 with an overall balanced accuracy of 98.8% and a threshold value of 0.2347, was sufficient to classify metastasis. HSPB8 mRNA expression was significantly increased following HOXB13 depletion in multiple metastatic CRPC models. Increased expression of HSPB8 by the microtubule inhibitor Colchicine or by exogenous means suppressed migration of mCRPC cells. Collectively, our results indicate that HOXB13 promotes metastasis of PCs by coordinated regulation of mitotic kinases and blockade of a putative tumor suppressor gene.

Paralogous HOX13 Genes in Human Cancers

Hox genes (HOX in humans), an evolutionary preserved gene family, are key determinants of embryonic development and cell memory gene program. Hox genes are organized in four clusters on four chromosomal loci aligned in 13 paralogous groups based on sequence homology (Hox gene network). During development Hox genes are transcribed, according to the rule of “spatio-temporal collinearity”, with early regulators of anterior body regions located at the 3’ end of each Hox cluster and the later regulators of posterior body regions placed at the distal 5’ end. The onset of 3’ Hox gene activation is determined by Wingless-type MMTV integration site family (Wnt) signaling, whereas 5’ Hox activation is due to paralogous group 13 genes, which act as posterior-inhibitors of more anterior Hox proteins (posterior prevalence). Deregulation of HOX genes is associated with developmental abnormalities and different human diseases. Paralogous HOX13 genes (HOX A13, HOX B13, HOX C13 and HOX D13) also play a relevant role in tumor development and progression. In this review, we will discuss the role of paralogous HOX13 genes regarding their regulatory mechanisms during carcinogenesis and tumor progression and their use as biomarkers for cancer diagnosis and treatment.

Homeobox B13 G84E Mutation and Prostate Cancer Risk

BACKGROUND

The homeobox B13 (HOXB13) G84E mutation has been recommended for use in genetic counselling for prostate cancer (PCa), but the magnitude of PCa risk conferred by this mutation is uncertain.

OBJECTIVE

To obtain precise risk estimates for mutation carriers and information on how these vary by family history and other factors.

DESIGN, SETTING, AND PARTICIPANTS

Two-fold: a systematic review and meta-analysis of published risk estimates, and a kin-cohort study comprising pedigree data on 11983 PCa patients enrolled during 1993-2014 from 189 UK hospitals and who had been genotyped for HOXB13 G84E.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS

Relative and absolute PCa risks. Complex segregation analysis with ascertainment adjustment to derive age-specific risks applicable to the population, and to investigate how these vary by family history and birth cohort.

RESULTS AND LIMITATIONS

A meta-analysis of case-control studies revealed significant heterogeneity between reported relative risks (RRs; range: 0.95-33.0, p<0.001) and differences by case selection (p=0.007). Based on case-control studies unselected for PCa family history, the pooled RR estimate was 3.43 (95% confidence interval [CI] 2.78-4.23). In the kin-cohort study, PCa risk for mutation carriers varied by family history (p<0.001). There was a suggestion that RRs decrease with age, but this was not significant (p=0.068). We found higher RR estimates for men from more recent birth cohorts (p=0.004): 3.09 (95% CI 2.03-4.71) for men born in 1929 or earlier and 5.96 (95% CI 4.01-8.88) for men born in 1930 or later. The absolute PCa risk by age 85 for a male HOXB13 G84E carrier varied from 60% for those with no PCa family history to 98% for those with two relatives diagnosed at young ages, compared with an average risk of 15% for noncarriers. Limitations include the reliance on self-reported cancer family history.

CONCLUSIONS

PCa risks for HOXB13 G84E mutation carriers are heterogeneous. Counselling should not be based on average risk estimates but on age-specific absolute risk estimates tailored to individual mutation carriers' family history and birth cohort.

PATIENT SUMMARY

Men who carry a hereditary mutation in the homeobox B13 (HOXB13) gene have a higher than average risk for developing prostate cancer. In our study, we examined a large number of families of men with prostate cancer recruited across UK hospitals, to assess what other factors may contribute to this risk and to assess whether we could create a precise model to help in predicting a man's prostate cancer risk. We found that the risk of developing prostate cancer in men who carry this genetic mutation is also affected by a family history of prostate cancer and their year of birth. This information can be used to assess more personalised prostate cancer risks to men who carry HOXB13 mutations and hence better counsel them on more personalised risk management options, such as tailoring prostate cancer screening frequency.

HOXB13 interaction with MEIS1 modifies proliferation and gene expression in prostate cancer

BACKGROUND

The recurrent p.Gly84Glu germline mutation (G84E) in HOXB13 is consistently associated with prostate cancer (PCa), although the mechanisms underlying such linkage remain elusive. The majority of the PCa-associated HOXB13 mutations identified are localized to two conserved domains in HOXB13 that have been shown to mediate the interaction with MEIS cofactors belonging to the TALE family of homeodomain transcription factors. In this study, we sought to interrogate the biochemical and functional interactions between HOXB13 and MEIS in prostatic cells with a goal of defining how the HOXB13-MEIS complex impacts PCa pathobiology and define the extent to which the oncogenic activity of G84E is related to its effect on HOXB13-MEIS interaction/function.

METHODS

HOXB13 and MEIS paralog expression in prostate epithelial cells and PCa cell lines was characterized by qPCR and immunoblot analyses. HOXB13 and MEIS1 co-expression in human prostate tissue was confirmed by IHC, followed by co-IP mapping of HOXB13-MEIS1 interactions. Proliferation of the PCa cell line LAPC4 following shRNA-mediated knockdown of each gene or both genes was assessed using DNA- and metabolic-based assays. Transcriptional targets of HOXB13 and MEIS1 were identified by gene expression profiling and qPCR. Finally, protein stability of HOXB13 in the context of MEIS1 was determined using pulse-chase assays.

RESULTS

HOXB13 and MEIS1 are co-expressed and interact in prostate cells. Both of the putative MEIS interacting domains (MID) within HOXB13 were shown to be capable of mediating the interaction between HOXB13 and MEIS1 independently and such interactions were not influenced by the G84E mutation. The inhibitory effect of either HOXB13 or MEIS1 knockdown on cellular proliferation was augmented by knockdown of both genes, and MEIS1 knockdown abolished HOXB13-driven regulation of BCHE and TNFSF10 mRNA expression. Notably, we demonstrated that MEIS1 stabilized the HOXB13 protein in LAPC4 cells.

CONCLUSIONS

Our study provides evidence for functional HOXB13-MEIS1 interactions in PCa. MEIS1 may contribute to the cancer-promoting actions of HOXB13 in cellular proliferation and gene regulation by prolonging HOXB13 half-life. Our data demonstrates that G84E is not a loss-of-function mutation that interferes with HOXB13 stability or ability to interact with MEIS1.

Roles of the HOXA10 gene during castrate-resistant prostate cancer progression

Homeobox A10 (HOXA10) is an important transcription factor that regulates the development of the prostate gland. However, it remains unknown whether it modulates prostate cancer (PCa) progression into castrate-resistant stages. In this study, we have applied RNA in situ hybridization assays to demonstrate that downregulation of HOXA10 expression is associated with castrate-resistant PCa. These findings are supported by public RNA-seq data showing that reduced HOXA10 expression is correlated with poor patient survival. We show that HOXA10 suppresses PCa cell proliferation, anchorage colony formation and xenograft growth independent to androgens. Using AmpliSeq transcriptome sequencing, we have found that gene groups associated with lipid metabolism and androgen receptor (AR) signaling are enriched in the HOXA10 transcriptome. Furthermore, we demonstrate that HOXA10 suppresses the transcription of the fatty acid synthase (FASN) gene by forming a protein complex with AR and prevents AR recruitment to the FASN gene promoter. These results lead us to conclude that downregulation of HOXA10 gene expression may enhance lipogenesis to promote PCa cell growth and tumor progression to castrate-resistant stage.

Targeting the BRD4-HOXB13 Coregulated Transcriptional Networks with Bromodomain-Kinase Inhibitors to Suppress Metastatic Castration-Resistant Prostate Cancer

Resistance to androgen receptor (AR) antagonists is a significant problem in the treatment of castration-resistant prostate cancers (CRPC). Identification of the mechanisms by which CRPCs evade androgen deprivation therapies (ADT) is critical to develop novel therapeutics. We uncovered that CRPCs rely on BRD4-HOXB13 epigenetic reprogramming for androgen-independent cell proliferation. Mechanistically, BRD4, a member of the BET bromodomain family, epigenetically promotes HOXB13 expression. Consistently, genetic disruption of HOXB13 or pharmacological suppression of its mRNA and protein expression by the novel dual-activity BET bromodomain-kinase inhibitors directly correlates with rapid induction of apoptosis, potent inhibition of tumor cell proliferation and cell migration, and suppression of CRPC growth. Integrative analysis revealed that the BRD4-HOXB13 transcriptome comprises a proliferative gene network implicated in cell-cycle progression, nucleotide metabolism, and chromatin assembly. Notably, although the core HOXB13 target genes responsive to BET inhibitors (HOTBIN10) are overexpressed in metastatic cases, in ADT-treated CRPC cell lines and patient-derived circulating tumor cells (CTC) they are insensitive to AR depletion or blockade. Among the HOTBIN10 genes, AURKB and MELK expression correlates with HOXB13 expression in CTCs of mCRPC patients who did not respond to abiraterone (ABR), suggesting that AURKB inhibitors could be used additionally against high-risk HOXB13-positive metastatic prostate cancers. Combined, our study demonstrates that BRD4-HOXB13-HOTBIN10 regulatory circuit maintains the malignant state of CRPCs and identifies a core proproliferative network driving ADT resistance that is targetable with potent dual-activity bromodomain-kinase inhibitors.

The lncRNA male-specific abdominal plays a critical role in Drosophila accessory gland development and male fertility

Although thousands of long non-coding RNAs (lncRNA) have been identified in the genomes of higher eukaryotes, the precise function of most of them is still unclear. Here, we show that a >65 kb, male-specific, lncRNA, called male-specific abdominal (msa) is required for the development of the secondary cells of the Drosophila male accessory gland (AG). msa is transcribed from within the Drosophila bithorax complex and shares much of its sequence with another lncRNA, the iab-8 lncRNA, which is involved in the development of the central nervous system (CNS). Both lncRNAs perform much of their functions via a shared miRNA embedded within their sequences. Loss of msa, or of the miRNA it contains, causes defects in secondary cell morphology and reduces male fertility. Although both lncRNAs express the same miRNA, the phenotype in the secondary cells and the CNS seem to reflect misregulation of different targets in the two tissues.

In many animals, the male seminal fluid induces physiology changes in the mated female that increase a male’s reproductive success. These changes are often referred to as the post-mating response (PMR). In Drosophila, the seminal fluid proteins responsible for generating the PMR are made in a specialized gland, analogous to the mammalian seminal vesicle and prostate, called the accessory gland (AG). In this work, we show that a male-specific, long, non-coding RNA (lncRNA), called msa, plays a critical role in the development and function of this gland, primarily through a microRNA (miRNA) encoded within its sequence. This same miRNA had previously been shown to be expressed in the central nervous system (CNS) via an alternative promoter, where its ability to repress homeotic genes is required for both male and female fertility. Here, we present evidence that the targets of this miRNA in the AG are likely different from those found in the CNS. Thus, the same miRNA seems to have been selected to affect Drosophila fertility through two different mechanisms. Although many non-coding RNAs have now been identified, very few can be shown to have function. Our work highlights a lncRNA that has multiple biological functions, affecting cellular morphology and fertility.

Prostate Organogenesis

The prostate is a male exocrine gland that secretes components of the seminal fluid. In men, prostate tumors are one of the most prevalent cancers. Studies on the development of the prostate have given a better understanding of the processes and genes that are important in the formation of this organ and have provided insights into the mechanisms of prostate tumorigenesis. These developmental studies have provided evidence that some of the genes and signaling pathways involved in development are reactivated or deregulated during prostate cancer. The prostate goes through a number of different stages during organogenesis, which include organ specification, epithelial budding, branching morphogenesis, canalization, and cytodifferentiation. During development, these processes are tightly regulated, many of which are controlled by the male hormone androgens. The majority of prostate tumors remain hormone regulated, and antiandrogen therapy is a first-line therapy, highlighting the important link between prostate organogenesis and cancer. In this review, we describe some of the data on genes that have important roles during prostate development that also have strong evidence linking them to prostate cancer.

The androgen receptor malignancy shift in prostate cancer

BACKGROUND

Androgens and the androgen receptor (AR) are necessary for the development, function, and homeostatic growth regulation of the prostate gland. However, once prostate cells are transformed, the AR is necessary for the proliferation and survival of the malignant cells. This change in AR function appears to occur in nearly every prostate cancer. We have termed this the AR malignancy shift.

METHODS

In this review, we summarize the current knowledge of the AR malignancy shift, including the DNA-binding patterns that define the shift, the transcriptome changes associated with the shift, the putative drivers of the shift, and its clinical implications.

RESULTS

In benign prostate epithelial cells, the AR primarily binds consensus AR binding sites. In carcinoma cells, the AR cistrome is dramatically altered, as the AR associates with FOXA1 and HOXB13 motifs, among others. This shift leads to the transcription of genes associated with a malignant phenotype. In model systems, some mutations commonly found in localized prostate cancer can alter the AR cistrome, consistent with the AR malignancy shift. Current evidence suggests that the AR malignancy shift is necessary but not sufficient for transformation of prostate epithelial cells.

CONCLUSIONS

Reinterpretation of prostate cancer genomic classification systems in light of the AR malignancy shift may improve our ability to predict clinical outcomes and treat patients appropriately. Identifying and targeting the molecular factors that contribute to the AR malignancy shift is not trivial but by doing so, we may be able to develop new strategies for the treatment or prevention of prostate cancer.

Two distinct DNA sequences recognized by transcription factors represent enthalpy and entropy optima

Most transcription factors (TFs) can bind to a population of sequences closely related to a single optimal site. However, some TFs can bind to two distinct sequences that represent two local optima in the Gibbs free energy of binding (ΔG). To determine the molecular mechanism behind this effect, we solved the structures of human HOXB13 and CDX2 bound to their two optimal DNA sequences, CAATAAA and TCGTAAA. Thermodynamic analyses by isothermal titration calorimetry revealed that both sites were bound with similar ΔG. However, the interaction with the CAA sequence was driven by change in enthalpy (ΔH), whereas the TCG site was bound with similar affinity due to smaller loss of entropy (ΔS). This thermodynamic mechanism that leads to at least two local optima likely affects many macromolecular interactions, as ΔG depends on two partially independent variables ΔH and ΔS according to the central equation of thermodynamics, ΔG = ΔH - TΔS.

Genes are sections of DNA that carry the instructions needed to build other molecules including all the proteins that the cell needs to fulfill its role. The information in the DNA is stored as a code consisting of four chemical bases, often referred to simply as “A”, “C”, “G” and “T”. The order or sequence of these bases determines the role of a protein. Many organisms – including humans – are built of many different types of cells that perform unique roles. Almost all cells carry the same genetic information, but proteins called transcription factors can regulate the activity of genes so that only a relevant subset of genes is switched on at a particular time.

Transcription factors glide along DNA and bind to short DNA sequences by attaching to the DNA bases directly or through bridges made up of water molecules. Two physical concepts known as enthalpy and entropy determine the strength of the connection. Enthalpy relates to how strong the chemical bonds that form between the transcription factors and the DNA bases are, compared to a situation where the transcription factor and DNA do not form a complex and bind to water molecules around them. Entropy measures the disorder of the system – the more disordered the solvent and protein-DNA complex are compared to solvent-containing free DNA and protein, the stronger the binding. A water molecule that bridges a DNA base with an amino-acid of a protein contributes to enthalpy, but results in loss of entropy, because the system becomes more ordered since the water molecule can no longer move freely.

Most transcription factors can only bind to DNA sequences that are very similar to each other, but some transcription factors can recognize several different kinds of sequences, and until now it was not clear how they could do this. Morgunova et al. studied four different human transcription factors that can each bind to two distinct DNA sequences. The results showed that the transcription factors bound to both DNA sequences with similar strength, but via different mechanisms. For one DNA sequence, an enthalpy-based mechanism essentially ‘froze’ the transcription factor to the DNA through rigid water bridges. The other DNA sequence was bound equally strongly but through moving water molecules, because this increased the entropy of the system. It is possible that these mechanisms could also apply to many other molecules that interact with each other through water-molecule bridges.

A better knowledge of the chemical bonds between transcription factors and DNA bases may in future help efforts to develop new treatments that depend on molecules being able to bind to other molecules. In addition, these findings may one day help scientists to predict how strongly two molecules will interact simply by knowing the structures of the molecules involved.

Morphoregulatory pathways in prostate ductal development

The mouse prostate is a male sex-accessory gland comprised of a branched ductal network arranged into three separate bilateral lobes: the anterior, dorsolateral, and ventral lobes. Prostate ductal development is the primary morphogenetic event in prostate development and requires a complex regulation of spatiotemporal factors. This review provides an overview of prostate development and the major genetic regulators and signaling pathways involved. To identify new areas for further study, we briefly highlight the likely important, but relatively understudied, role of the extracellular matrix (ECM). Finally, we point out the potential importance of the ECM in influencing the behavior and prognosis of prostate cancer. Developmental Dynamics 246:89-99, 2017. © 2016 Wiley Periodicals, Inc.

Impact of the G84E variant on HOXB13 gene and protein expression in formalin-fixed, paraffin-embedded prostate tumours

The HOXB13 G84E variant is associated with risk of prostate cancer (PCa), however the role this variant plays in PCa development is unknown. This study examined 751 cases, 450 relatives and 355 controls to determine the contribution of this variant to PCa risk in Tasmania and investigated HOXB13 gene and protein expression in tumours from nine G84E heterozygote variant and 13 wild-type carriers. Quantitative PCR and immunohistochemistry showed that HOXB13 gene and protein expression did not differ between tumour samples from variant and wild-type carriers. Allele-specific transcription revealed that two of seven G84E carriers transcribed both the variant and wild-type allele, while five carriers transcribed the wild-type allele. Methylation of surrounding CpG sites was lower in the variant compared to the wild-type allele, however overall methylation across the region was very low. Notably, tumour characteristics were less aggressive in the two variant carriers that transcribed the variant allele compared to the five that did not. This study has shown that HOXB13 expression does not differ between tumour tissue of G84E variant carriers and non-carriers. Intriguingly, the G84E variant allele was rarely transcribed in carriers, suggesting that HOXB13 expression may be driven by the wild-type allele in the majority of carriers.

Measuring quantitative effects of methylation on transcription factor-DNA binding affinity

Methylation of CpG (cytosine-phosphate-guanine) dinucleotides is a common epigenetic mark that influences gene expression. The effects of methylation on transcription factor (TF) binding are unknown for most TFs and, even when known, such knowledge is often only qualitative. In reality, methylation sensitivity is a quantitative effect, just as changes to the DNA sequence have quantitative effects on TF binding affinity. We describe Methyl-Spec-seq, an easy-to-use method that measures the effects of CpG methylation (mCPG) on binding affinity for hundreds to thousands of variants in parallel, allowing one to quantitatively assess the effects at every position in a binding site. We demonstrate its use on several important DNA binding proteins. We calibrate the accuracy of Methyl-Spec-seq using a novel two-color competitive fluorescence anisotropy method that can accurately determine the relative affinities of two sequences in solution. We also present software that extends standard methods for representing, visualizing, and searching for matches to binding site motifs to include the effects of methylation. These tools facilitate the study of the consequences for gene regulation of epigenetic marks on DNA.

Somatic molecular subtyping of prostate tumors from HOXB13 G84E carriers

A recurrent germline mutation (G84E) in the HOXB13 gene is associated with early onset and family history-positive prostate cancer in patients of European descent, occurring in up to 5% of prostate cancer families. To date, the molecular features of prostate tumors occurring in HOXB13 G84E carriers have not been studied in a large cohort of patients. We identified 101 heterozygous carriers of G84E who underwent radical prostatectomy for prostate cancer between 1985 and 2011 and matched these men by race, age and tumor grade to 99 HOXB13 wild-type controls. Immunostaining for HOXB13, PTEN, ERG, p53 and SPINK1 as well as RNA in situ hybridization for ETV1/4/5 were performed using genetically validated assays. Tumors from G84E carriers generally expressed HOXB13 protein at a level comparable to benign and wild-type glands. ETS gene expression (either ERG or ETV1/4/5) was seen in 36% (36/101) of tumors from G84E carriers compared to 68% (65/96) of the controls (p < 0.0001). PTEN was lost in 11% (11/101) of G84E carriers compared to 25% (25/99) of the controls (p = 0.014). PTEN loss was enriched among ERG-positive compared to ERG-negative tumors in both groups of patients. Nuclear accumulation of the p53 protein, indicative of underlying TP53 missense mutations, was uncommon in both groups, occurring in 1% (1/101) of the G84E carriers versus 2% (2/92) of the controls (p = NS). Taken together, these data suggest that genes other than ERG and PTEN may drive carcinogenesis/progression in the majority of men with germline HOXB13 mutations.

Prostate organogenesis: tissue induction, hormonal regulation and cell type specification

Prostate organogenesis is a complex process that is primarily mediated by the presence of androgens and subsequent mesenchyme-epithelial interactions. The investigation of prostate development is partly driven by its potential relevance to prostate cancer, in particular the apparent re-awakening of key developmental programs that occur during tumorigenesis. However, our current knowledge of the mechanisms that drive prostate organogenesis is far from complete. Here, we provide a comprehensive overview of prostate development, focusing on recent findings regarding sexual dimorphism, bud induction, branching morphogenesis and cellular differentiation.