Aberrant Activation of a Gastrointestinal Transcriptional Circuit in Prostate Cancer Mediates Castration Resistance

Targeting FGFR4 Abrogates HNF1A-driven Metastasis in Pancreatic Ductal Adenocarcinoma

Purpose

We previously identified an oncogenic role for the transcription factor HNF1A in pancreatic ductal adenocarcinoma (PDAC). However, the role of HNF1A in the metastatic progression of PDAC remains unknown and targeting modalities for HNF1A -dependent phenotypes have yet to be identified.

Experimental Design

Transwell chambers were used to assess the effects of HNF1A and FGFR4 modulation on the migration and invasion of ATCC and patient-derived PDAC cells in vitro . An intrasplenic injection xenograft model was used to evaluate the impact of HNF1A knockdown and overexpression on metastatic tumor burden. Single-cell RNA sequencing, tissue microarray (TMA) data, and UMAP spatial profiling were used to identify FGFR4 as an HNF1A target gene upregulated in metastatic cells. RNAi and two FGFR4 inhibiting modalities (H3B-6527 and U3- 1784) were utilized to demonstrate the efficacy of FGFR4 inhibiting agents at reducing HNF1A- driven metastasis.

Results

Knockdown of HNF1A significantly decreases and HNF1A overexpression significantly increases PDAC cell migration and invasion. In vivo studies show that HNF1A knockdown significantly abrogates metastasis, while overexpression significantly promotes metastasis. Single-cell RNAseq shows that FGFR4 is upregulated in metastatic PDAC cells and staining for HNF1A and FGFR4 in a PDAC TMA reveals significant correlation between HNF1A and FGFR4 in PDAC patients. Further, knockdown and inhibition of FGFR4 significantly decreases HNF1A- mediated cell migration and invasion, and blocks HNF1A-driven metastasis in vivo .

Conclusions

These findings demonstrate that HNF1A drives PDAC metastasis via upregulation of FGFR4, and FGFR4 inhibition is a potential mechanism to target metastasis in PDAC patients.

Translational Relevance

Pancreatic ductal adenocarcinoma (PDAC) is one of the most lethal malignancies, made even more devastating when metastases overwhelm major organs. The vast majority of PDAC patients either present with metastases or will relapse with recurrent metastatic PDAC after primary tumor resection. Unfortunately, toxic and largely ineffective chemotherapies are currently the only approved treatment options for these patients and therefore there exists a critical and unmet clinical need for targeted therapies against pro-metastatic pathways in PDAC. In the current study, we identify HNF1A as an oncogenic transcription factor that drives metastasis in PDAC, and it does so through upregulation of the receptor tyrosine kinase FGFR4. Importantly, FGFR4 is a targetable vulnerability and treatment with an FGFR4 blocking antibody reduces HNF1A-driven metastasis. These findings suggest that FGFR4 inhibitors could be an efficacious treatment for PDAC patients for the prevention or delay of metastatic tumor development.

Loss of Kmt2c or Kmt2d primes urothelium for tumorigenesis and redistributes KMT2A-menin to bivalent promoters

Members of the KMT2C/D-KDM6A complex are recurrently mutated in urothelial carcinoma and in histologically normal urothelium. Here, using genetically engineered mouse models, we demonstrate that Kmt2c/d knockout in the urothelium led to impaired differentiation, augmented responses to growth and inflammatory stimuli and sensitization to oncogenic transformation by carcinogen and oncogenes. Mechanistically, KMT2D localized to active enhancers and CpG-poor promoters that preferentially regulate the urothelial lineage program and Kmt2c/d knockout led to diminished H3K4me1, H3K27ac and nascent RNA transcription at these sites, which leads to impaired differentiation. Kmt2c/d knockout further led to KMT2A-menin redistribution from KMT2D localized enhancers to CpG-high and bivalent promoters, resulting in derepression of signal-induced immediate early genes. Therapeutically, Kmt2c/d knockout upregulated epidermal growth factor receptor signaling and conferred vulnerability to epidermal growth factor receptor inhibitors. Together, our data posit that functional loss of Kmt2c/d licenses a molecular 'field effect' priming histologically normal urothelium for oncogenic transformation and presents therapeutic vulnerabilities.

Castrate-resistant prostate cancer response to taxane is determined by an HNF1-dependent apoptosis resistance circuit

Metastatic castrate-resistant prostate cancer (mCRPC) is a genetically and phenotypically heterogeneous cancer where advancements are needed in biomarker discovery and targeted therapy. A critical and often effective component of treatment includes taxanes. We perform a high-throughput screen across a cohort of 30 diverse patient-derived castrate-resistant prostate cancer (CRPC) organoids to a library of 78 drugs. Combining quantitative response measures with transcriptomic analyses demonstrates that HNF1 homeobox A (HNF1A) drives a transcriptional program of taxane resistance, commonly dependent upon cellular inhibitor of apoptosis protein 2 (cIAP2). Monotherapy with cIAP2 inhibitor LCL161 is sufficient to treat HNF1A+ models of mCRPC previously resistant to docetaxel. These data may be useful in future clinical trial designs.

LBFextract: Unveiling transcription factor dynamics from liquid biopsy data

Motivation

The analysis of circulating cell-free DNA (cfDNA) holds immense promise as a non-invasive diagnostic tool across various human conditions. However, extracting biological insights from cfDNA fragments entails navigating complex and diverse bioinformatics methods, encompassing not only DNA sequence variation, but also epigenetic characteristics like nucleosome footprints, fragment length, and methylation patterns.

Results

We introduce Liquid Biopsy Feature extract (LBFextract), a comprehensive package designed to streamline feature extraction from cfDNA sequencing data, with the aim of enhancing the reproducibility and comparability of liquid biopsy studies. LBFextract facilitates the integration of preprocessing and postprocessing steps through alignment fragment tags and a hook mechanism. It incorporates various methods, including coverage-based and fragment length-based approaches, alongside two novel feature extraction methods: an entropy-based method to infer TF activity from fragmentomics data and a technique to amplify signals from nucleosome dyads. Additionally, it implements a method to extract condition-specific differentially active TFs based on these features for biomarker discovery. We demonstrate the use of LBFextract for the subtype classification of advanced prostate cancer patients using coverage signals at transcription factor binding sites from cfDNA. We show that LBFextract can generate robust and interpretable features that can discriminate between different clinical groups. LBFextract is a versatile and user-friendly package that can facilitate the analysis and interpretation of liquid biopsy data.

Data and Code Availability and Implementation

LBFextract is freely accessible at https://github.com/Isy89/LBF. It is implemented in Python and compatible with Linux and Mac operating systems. Code and data to reproduce these analyses have been uploaded to 10.5281/zenodo.10964406.

Single-cell analysis of treatment-resistant prostate cancer: Implications of cell state changes for cell surface antigen-targeted therapies

Targeting cell surface molecules using radioligand and antibody-based therapies has yielded considerable success across cancers. However, it remains unclear how the expression of putative lineage markers, particularly cell surface molecules, varies in the process of lineage plasticity, wherein tumor cells alter their identity and acquire new oncogenic properties. A notable example of lineage plasticity is the transformation of prostate adenocarcinoma (PRAD) to neuroendocrine prostate cancer (NEPC)-a growing resistance mechanism that results in the loss of responsiveness to androgen blockade and portends dismal patient survival. To understand how lineage markers vary across the evolution of lineage plasticity in prostate cancer, we applied single-cell analyses to 21 human prostate tumor biopsies and two genetically engineered mouse models, together with tissue microarray analysis on 131 tumor samples. Not only did we observe a higher degree of phenotypic heterogeneity in castrate-resistant PRAD and NEPC than previously anticipated but also found that the expression of molecules targeted therapeutically, namely PSMA, STEAP1, STEAP2, TROP2, CEACAM5, and DLL3, varied within a subset of gene-regulatory networks (GRNs). We also noted that NEPC and small cell lung cancer subtypes shared a set of GRNs, indicative of conserved biologic pathways that may be exploited therapeutically across tumor types. While this extreme level of transcriptional heterogeneity, particularly in cell surface marker expression, may mitigate the durability of clinical responses to current and future antigen-directed therapies, its delineation may yield signatures for patient selection in clinical trials, potentially across distinct cancer types.

FOXA1 regulates ribosomal RNA transcription in prostate cancer

BACKGROUND

Ribosome biogenesis is excessively activated in tumor cells, yet it is little known whether oncogenic transcription factors (TFs) are involved in the ribosomal RNA (rRNA) transactivation.

METHODS

Nucleolar proteomics data and large-scale immunofluorescence were re-analyzed to jointly identify the proteins localized at nucleolus. RNA-Seq data of five prostate cancer (PCa) cohorts were combined and integrated with multi-dimensional data to define the upregulated nucleolar TFs in PCa tissues. Then, ChIP-Seq data of PCa cell lines and two PCa clinical cohorts were re-analyzed to reveal the TF binding patterns at ribosomal DNA (rDNA) repeats. The TF binding at rDNA was validated by ChIP-qPCR. The effect of the TF on rRNA transcription was determined by rDNA luciferase reporter, nascent RNA synthesis, and global protein translation assays.

RESULTS

In this study, we reveal the role of oncogenic TF FOXA1 in regulating rRNA transcription within nucleolar organization regions. By analyzing human TFs in prostate cancer clinical datasets and nucleolar proteomics data, we identified that FOXA1 is partially localized in the nucleolus and correlated with global protein translation. Our extensive FOXA1 ChIP-Seq analysis provides robust evidence of FOXA1 binding across rDNA repeats in prostate cancer cell lines, primary tumors, and castration-resistant variants. Notably, FOXA1 occupancy at rDNA repeats correlates with histone modifications associated with active transcription, namely H3K27ac and H3K4me3. Reducing FOXA1 expression results in decreased transactivation at rDNA, subsequently diminishing global protein synthesis.

CONCLUSIONS

Our results suggest FOXA1 regulates aberrant ribosome biogenesis downstream of oncogenic signaling in prostate cancer.

Molecular subtypes of neuroendocrine carcinomas: A cross-tissue classification framework based on five transcriptional regulators

Neuroendocrine carcinomas (NECs) are extremely lethal malignancies that can arise at almost any anatomic site. Characterization of NECs is hindered by their rarity and significant inter- and intra-tissue heterogeneity. Herein, through an integrative analysis of over 1,000 NECs originating from 31 various tissues, we reveal their tissue-independent convergence and further unveil molecular divergence driven by distinct transcriptional regulators. Pan-tissue NECs are therefore categorized into five intrinsic subtypes defined by ASCL1, NEUROD1, HNF4A, POU2F3, and YAP1. A comprehensive portrait of these subtypes is depicted, highlighting subtype-specific transcriptional programs, genomic alterations, evolution trajectories, therapeutic vulnerabilities, and clinicopathological presentations. Notably, the newly discovered HNF4A-dominated subtype-H exhibits a gastrointestinal-like signature, wild-type RB1, unique neuroendocrine differentiation, poor chemotherapeutic response, and prevalent large-cell morphology. The proposal of uniform classification paradigm illuminates transcriptional basis of NEC heterogeneity and bridges the gap across different lineages and cytomorphological variants, in which context-dependent prevalence of subtypes underlies their phenotypic disparities.

A prostate cancer gastrointestinal transcriptional phenotype may be associated with diminished response to AR-targeted therapy

Background

Prostate cancer is a heterogenous disease, but once it becomes metastatic it eventually becomes treatment resistant. One mechanism of resistance to AR-targeting therapy is lineage plasticity, where the tumor undergoes a transformation to an AR-indifferent phenotype, most studied in the context of neuroendocrine prostate cancer (NEPC). However, activation of additional de- or trans-differentiation programs, including a gastrointestinal (GI) gene expression program, has been suggested as an alternative method of resistance. In this study, we explored the previously identified GI prostate cancer phenotype (PCa-GI) in a large cohort of metastatic castration-resistant prostate cancer (mCRPC) patient biopsy samples.

Methods

We analyzed a dataset of 634 mCRPC samples with batch effect corrected gene expression data from the West Coast Dream Team (WCDT), the East Coast Dream Team (ECDT), the Fred Hutchinson Cancer Research Center (FHCRC) and the Weill Cornell Medical center (WCM). Survival data was available from the WCDT and ECDT cohorts. We calculated a gene expression GI score using the sum of z-scores of genes from a published set of PCa-GI-defining genes (N=38). Survival analysis was performed using the Kaplan-Meier method and Cox proportional hazards regression with endpoint overall survival from time of biopsy to death of any cause.

Results

We found that the PCa-GI score had a bimodal distribution, identifying a distinct set of tumors with an activated GI expression pattern. Approximately 35% of samples were classified as PCa-GI high, which was concordant with prior reports. Liver metastases had the highest median score but after excluding liver samples, 29% of the remaining samples were still classified as PCa-GI high, suggesting a distinct phenotype not exclusive to liver metastases. No correlation was observed between GI score and proliferation, AR signaling, or NEPC scores. Furthermore, the PCa-GI score was not associated with genomic alterations in AR, FOXA1, RB1, TP53 or PTEN. However, tumors with MYC amplifications showed significantly higher GI scores (p=0.0001). Patients with PCa-GI tumors had a shorter survival (HR=1.5 [1.1-2.1], p=0.02), but this result was not significant after adjusting for the liver as metastatic site (HR=1.2 [0.82-1.7], p=0.35). Patients with PCa-GI low samples had a better outcome after androgen receptor signaling inhibitors (ASI, abiraterone or enzalutamide) than other therapies (HR=0.37 [0.22-0.61], p=0.0001) while the benefit of ASI was smaller and non-significant for PCa-GI high samples (HR=0.55 [0.29-1.1], p=0.07). A differential pathway analysis identified FOXA2 signaling to be upregulated PCa-GI high tumors (FDR = 3.7 × 10-13).

Conclusions

The PCa-GI phenotype is prevalent in clinical mCRPC samples and may represent a distinct biological entity. PCa-GI tumors may respond less to ASI and could offer a strategy to study novel therapeutic targets.

Understanding the molecular regulators of neuroendocrine prostate cancer

Worldwide, prostate cancer (PCa) remains a leading cause of death in men. Histologically, the majority of PCa cases are classified as adenocarcinomas, which are mainly composed of androgen receptor-positive luminal cells. PCa is initially driven by the androgen receptor axis, where androgen-mediated activation of the receptor is one of the primary culprits for disease progression. Therefore, in advanced stage PCa, patients are generally treated with androgen deprivation therapies alone or in combination with androgen receptor pathway inhibitors. However, after an initial decrease, the cancer recurs for majority patients. At this stage, cancer is known as castration-resistant prostate cancer (CRPC). Majority of CRPC tumors still depend on androgen receptor axis for its progression to metastasis. However, in around 20-30% of cases, CRPC progresses via an androgen receptor-independent pathway and is often presented as neuroendocrine cancer (NE). This NE phenotype is highly aggressive with poor overall survival as compared to CRPC adenocarcinoma. NE cancers are resistant to standard taxane chemotherapies, which are often used to treat metastatic disease. Pathologically and morphologically, NE cancers are highly diverse and often co-exist with adenocarcinoma. Due to the lack of proper biomarkers, it is often difficult to make an early diagnosis of this lethal disease. Moreover, increased tumor heterogeneity and admixtures of adeno and NE subtypes in the same tumor make early detection of NE tumors very difficult. With the advancement of our knowledge and sequencing technology, we are now able to better understand the molecular mediators of this transformation pathway. This current study will give an update on how various molecular regulators are involved in these lineage transformation processes and what challenges we are still facing to detect and treat this cancer.

Alternative splicing in prostate cancer progression and therapeutic resistance

Prostate cancer (CaP) remains the second leading cause of cancer deaths in western men. CaP mortality results from diverse molecular mechanisms that mediate resistance to the standard of care treatments for metastatic disease. Recently, alternative splicing has been recognized as a hallmark of CaP aggressiveness. Alternative splicing events cause treatment resistance and aggressive CaP behavior and are determinants of the emergence of the two major types of late-stage treatment-resistant CaP, namely castration-resistant CaP (CRPC) and neuroendocrine CaP (NEPC). Here, we review recent multi-omics data that are uncovering the complicated landscape of alternative splicing events during CaP progression and the impact that different gene transcript isoforms can have on CaP cell biology and behavior. We discuss renewed insights in the molecular machinery by which alternative splicing occurs and contributes to the failure of systemic CaP therapies. The potential for alternative splicing events to serve as diagnostic markers and/or therapeutic targets is explored. We conclude by considering current challenges and promises associated with splicing-modulating therapies, and their potential for clinical translation into CaP patient care.

Single Cell Analysis of Treatment-Resistant Prostate Cancer: Implications of Cell State Changes for Cell Surface Antigen Targeted Therapies

Targeting cell surface molecules using radioligand and antibody-based therapies has yielded considerable success across cancers. However, it remains unclear how the expression of putative lineage markers, particularly cell surface molecules, varies in the process of lineage plasticity, wherein tumor cells alter their identity and acquire new oncogenic properties. A notable example of lineage plasticity is the transformation of prostate adenocarcinoma (PRAD) to neuroendocrine prostate cancer (NEPC)--a growing resistance mechanism that results in the loss of responsiveness to androgen blockade and portends dismal patient survival. To understand how lineage markers vary across the evolution of lineage plasticity in prostate cancer, we applied single cell analyses to 21 human prostate tumor biopsies and two genetically engineered mouse models, together with tissue microarray analysis (TMA) on 131 tumor samples. Not only did we observe a higher degree of phenotypic heterogeneity in castrate-resistant PRAD and NEPC than previously anticipated, but also found that the expression of molecules targeted therapeutically, namely PSMA, STEAP1, STEAP2, TROP2, CEACAM5, and DLL3, varied within a subset of gene-regulatory networks (GRNs). We also noted that NEPC and small cell lung cancer (SCLC) subtypes shared a set of GRNs, indicative of conserved biologic pathways that may be exploited therapeutically across tumor types. While this extreme level of transcriptional heterogeneity, particularly in cell surface marker expression, may mitigate the durability of clinical responses to novel antigen-directed therapies, its delineation may yield signatures for patient selection in clinical trials, potentially across distinct cancer types.

An integrative proteomics approach identifies tyrosine kinase KIT as a therapeutic target for SPINK1-positive prostate cancer

Elevated serine peptidase inhibitor, Kazal type 1 (SPINK1) levels in ∼10%-25% of prostate cancer (PCa) patients associate with aggressive phenotype, for which there are limited treatment choices and dismal clinical outcomes. Using an integrative proteomics approach involving label-free phosphoproteome and proteome profiling, we delineated the downstream signaling pathways involved in SPINK1-mediated tumorigenesis and identified tyrosine kinase KIT as highly enriched. Furthermore, high to moderate levels of KIT expression were detected in ∼85% of SPINK1-positive PCa specimens. We show KIT signaling orchestrates SPINK1-mediated oncogenesis, and treatment with KIT inhibitor reduces tumor growth and metastases in preclinical mice models. Mechanistically, KIT signaling modulates WNT/β-catenin pathway and confers stemness-related features in PCa. Notably, inhibiting KIT signaling led to restoration of AR/REST levels, forming a feedback loop enabling SPINK1 repression. Overall, we uncover the role of KIT signaling downstream of SPINK1 in maintaining lineage plasticity and provide distinct treatment modalities for advanced-stage SPINK1-positive patients.

BET inhibitors as a therapeutic intervention in gastrointestinal gene signature-positive castration-resistant prostate cancer

A subgroup of castration-resistant prostate cancer (CRPC) aberrantly expresses a gastrointestinal (GI) transcriptome governed by two GI-lineage-restricted transcription factors, HNF1A and HNF4G. In this study, we found that expression of GI transcriptome in CRPC correlates with adverse clinical outcomes to androgen receptor signaling inhibitor treatment and shorter overall survival. Bromo- and extra-terminal domain inhibitors (BETi) downregulated HNF1A, HNF4G, and the GI transcriptome in multiple CRPC models, including cell lines, patient-derived organoids, and patient-derived xenografts, while AR and the androgen-dependent transcriptome were largely spared. Accordingly, BETi selectively inhibited growth of GI transcriptome-positive preclinical models of prostate cancer. Mechanistically, BETi inhibited BRD4 binding at enhancers globally, including both AR and HNF4G bound enhancers while gene expression was selectively perturbed. Restoration of HNF4G expression in the presence of BETi rescued target gene expression without rescuing BRD4 binding. This suggests that inhibition of master transcription factors expression underlies the selective transcriptional effects of BETi.

SIGNIFICANCE

GI transcriptome expression in CRPC is regulated by the HNF1A-HNF4G-BRD4 axis and correlates with worse clinical outcomes. Accordingly, BET inhibitors significantly reduce tumor cell growth in multiple GI-transcriptome-positive preclinical models of CRPC. Our studies point that expression of GI transcriptome could serve as a predictive biomarker to BETi therapy response.

Nucleosome Patterns in Circulating Tumor DNA Reveal Transcriptional Regulation of Advanced Prostate Cancer Phenotypes

Advanced prostate cancers comprise distinct phenotypes, but tumor classification remains clinically challenging. Here, we harnessed circulating tumor DNA (ctDNA) to study tumor phenotypes by ascertaining nucleosome positioning patterns associated with transcription regulation. We sequenced plasma ctDNA whole genomes from patient-derived xenografts representing a spectrum of androgen receptor active (ARPC) and neuroendocrine (NEPC) prostate cancers. Nucleosome patterns associated with transcriptional activity were reflected in ctDNA at regions of genes, promoters, histone modifications, transcription factor binding, and accessible chromatin. We identified the activity of key phenotype-defining transcriptional regulators from ctDNA, including AR, ASCL1, HOXB13, HNF4G, and GATA2. To distinguish NEPC and ARPC in patient plasma samples, we developed prediction models that achieved accuracies of 97% for dominant phenotypes and 87% for mixed clinical phenotypes. Although phenotype classification is typically assessed by IHC or transcriptome profiling from tumor biopsies, we demonstrate that ctDNA provides comparable results with diagnostic advantages for precision oncology.

SIGNIFICANCE

This study provides insights into the dynamics of nucleosome positioning and gene regulation associated with cancer phenotypes that can be ascertained from ctDNA. New methods for classification in phenotype mixtures extend the utility of ctDNA beyond assessments of somatic DNA alterations with important implications for molecular classification and precision oncology. This article is highlighted in the In This Issue feature, p. 517.

Enhancer profiling identifies epigenetic markers of endocrine resistance and reveals therapeutic options for metastatic castration-resistant prostate cancer patients

Androgen Receptor (AR) signaling inhibitors, including enzalutamide, are treatment options for patients with metastatic castration-resistant prostate cancer (mCRPC), but resistance inevitably develops. Using metastatic samples from a prospective phase II clinical trial, we epigenetically profiled enhancer/promoter activities with H3K27ac chromatin immunoprecipitation followed by sequencing, before and after AR-targeted therapy. We identified a distinct subset of H3K27ac-differentially marked regions that associated with treatment responsiveness. These data were successfully validated in mCRPC patient-derived xenograft models (PDX). In silico analyses revealed HDAC3 as a critical factor that can drive resistance to hormonal interventions, which we validated in vitro . Using cell lines and mCRPC PDX tumors in vitro , we identified drug-drug synergy between enzalutamide and the pan-HDAC inhibitor vorinostat, providing therapeutic proof-of-concept. These findings demonstrate rationale for new therapeutic strategies using a combination of AR and HDAC inhibitors to improve patient outcome in advanced stages of mCRPC.

HNF4G accelerates glioma progression by facilitating NRP1 transcription

Hepatocyte nuclear factor 4γ (HNF4G) is considered to be a transcription factor and functions as an oncogene in certain types of human cancer. However, the precise functions and the potential molecular mechanisms of HNF4G in glioma remain unclear. Therefore, the present study aimed to elucidate the role of HNF4G in glioma and the underlying mechanism. Western blotting and reverse transcription-quantitative PCR (RT-qPCR) demonstrated that HNF4G was highly expressed in glioma tissues and cell lines. The overexpression of HNF4G in LN229 and U251 glioma cells promoted cell proliferation and cell cycle progression, and inhibited apoptosis, while the knockdown of HNF4G suppressed cell proliferation, cell cycle progression and tumor growth, and induced apoptosis. A significant positive association was detected between HNF4G and neuropilin-1 (NRP1) mRNA expression in glioma tissues. Bioinformatics analysis, chromatin immunoprecipitation-RT-qPCR and promoter reporter assays confirmed that HNF4G promoted NRP1 transcription in glioma by binding to its promoter. NRP1 overexpression facilitated glioma cell proliferation and cell cycle progression, and suppressed apoptosis in vitro, while the knockdown of NRP1 inhibited cell proliferation and cell cycle progression, and facilitated apoptosis. NRP1 overexpression reversed the effects induced by HNF4G knockdown on glioma cell proliferation, cell cycle progression and apoptosis. In summary, the present study demonstrated that HNF4G promotes glioma cell proliferation and suppresses apoptosis by activating NRP1 transcription. These findings indicate that HNF4G acts as an oncogene in glioma and may thus be an effective therapeutic target for glioma.

HNF4G stimulates the development of pancreatic cancer by promoting IGF2BP2 transcription

OBJECTIVE

Pancreatic cancer is a devastating and lethal malignancy. Our study investigated the effective mechanism of HNF4G on pancreatic cancer cell functions through the IGF2BP2 transcription.

METHODS

HNF4G and IGF2BP2 expressions in pancreatic cancer were examined. The relationship between HNF4G expression and pancreatic cancer patients' clinicopathological characteristics was evaluated. After interfering with HNF4G expression in pancreatic cancer cells, the cell proliferative, migratory, and invasive capabilities were evaluated. Also, the expression of proliferation-related gene PCNA and migration and invasion-related gene MMP2 was determined. The binding relation between HNF4G and HNF4G promoter was forecasted and testified. A tumorigenesis assay in nude mice was performed to detect the HNF4G interference's effect on the subcutaneous tumorigenic capacity of pancreatic cancer cells.

RESULTS

HNF4G and IGF2BP2 expressions were up-regulated in pancreatic cancer. Specifically, interfering with HNF4G inhibited PANC-1 cell proliferative, invasive and migratory behaviors, and decreased PCNA and MMP2 expression. Mechanistically, HNF4G as a transcription factor could specifically bind to IGF2BP2 and promote its expression. Rescue assay findings showed that IGF2BP2 overexpression could reverse the inhibiting effect of HNF4G interference on pancreatic cancer cells. For the in vivo finding, interfering HNF4G expression retarded the subcutaneous tumorigenic ability of pancreatic cancer cells.

CONCLUSION

We summarize that HNF4G as a transcription factor regulates IGF2BP2 expression to promote pancreatic cancer cell proliferation and migration capacities.

Intrinsic Molecular Subtypes of Metastatic Castration-Resistant Prostate Cancer

PURPOSE

Although numerous biology-driven subtypes have been described previously in metastatic castration-resistant prostate cancer (mCRPC), unsupervised molecular subtyping based on gene expression has been less studied, especially using large cohorts. Thus, we sought to identify the intrinsic molecular subtypes of mCRPC and assess molecular and clinical correlates in the largest combined cohort of mCRPC samples with gene expression data available to date.

EXPERIMENTAL DESIGN

We combined and batch-effect corrected gene expression data from four mCRPC cohorts from the Fred Hutchinson Cancer Research Center (N = 157), a small-cell neuroendocrine (NE) prostate cancer (SCNC)-enriched cohort from Weill Cornell Medicine (N = 49), and cohorts from the Stand Up 2 Cancer/Prostate Cancer Foundation East Coast Dream Team (N = 266) and the West Coast Dream Team (N = 162).

RESULTS

Hierarchical clustering of RNA-sequencing data from these 634 mCRPC samples identified two distinct adenocarcinoma subtypes, one of which (adeno-immune) was characterized by higher gene expression of immune pathways, higher CIBERSORTx immune scores, diminished ASI benefit, and non-lymph node metastasis tropism compared with an adeno-classic subtype. We also identified two distinct subtypes with enrichment for an NE phenotype, including an NE-liver subgroup characterized by liver metastasis tropism, PTEN loss, and APC and SPOP mutations compared with an NE-classic subgroup.

CONCLUSIONS

Our results emphasize the heterogeneity of mCRPC beyond currently accepted molecular phenotypes, and suggest that future studies should consider incorporating transcriptome-wide profiling to better understand how these differences impact treatment responses and outcomes.

The roles of mutated SPINK1 gene in prostate cancer cells

SPINK1-positive prostate cancer (PCa) has been identified as an aggressive PCa subtype. However, there is a lack of definite studies to elucidate the underlying mechanism of the loss of SPINK1 expression in most PCa cells except 22Rv1 cells, which are derived from a human prostatic carcinoma xenograft, CWR22R. The aim of this study was to investigate the mechanisms of SPINK1 protein positive/negative expression and its biological roles in PCa cell lines. SPINK1 mRNA was highly expressed in 22Rv1 cells compared with LNCaP, C4-2B, DU145, and PC-3 cells, and the protein was only detected in 22Rv1 cells. Among these cell lines, the wild-type SPINK1 coding sequence was only found in 22Rv1 cells, and two mutation sites, the c.194G>A missense mutation and the c.210T>C synonymous mutation, were found in other cell lines. Our further research showed that the mutations were associated with a reduction in SPINK1 mRNA and protein levels. Functional experiments indicated that SPINK1 promoted PC-3 cell proliferation, migration, and invasion, while knockdown of SPINK1 attenuated 22Rv1 cell proliferation, migration, and invasion. The wild-type SPINK1 gene can promote the malignant behaviors of cells more than the mutated ones. Cell cycle analysis by flow cytometry showed that SPINK1 decreased the percentage of cells in the G0/G1 phase and increased the percentage of S phase cells. We demonstrated that the c.194G>A and c.210T>C mutations in the SPINK1 gene decreased the mRNA and protein levels. The wild-type SPINK1 gene is related to aggressive biological behaviors of PCa cells and may be a potential therapeutic target for PCa.

3D-Epigenomic Regulation of Gene Transcription in Hepatocellular Carcinoma

The fundamental cause of transcription dysregulation in hepatocellular carcinoma (HCC) remains elusive. To investigate the underlying mechanisms, comprehensive 3D-epigenomic analyses are performed in cellular models of THLE2 (a normal hepatocytes cell line) and HepG2 (a hepatocellular carcinoma cell line) using integrative approaches for chromatin topology, genomic and epigenomic variation, and transcriptional output. Comparing the 3D-epigenomes in THLE2 and HepG2 reveal that most HCC-associated genes are organized in complex chromatin interactions mediated by RNA polymerase II (RNAPII). Incorporation of genome-wide association studies (GWAS) data enables the identification of non-coding genetic variants that are enriched in distal enhancers connecting to the promoters of HCC-associated genes via long-range chromatin interactions, highlighting their functional roles. Interestingly, CTCF binding and looping proximal to HCC-associated genes appear to form chromatin architectures that overarch RNAPII-mediated chromatin interactions. It is further demonstrated that epigenetic variants by DNA hypomethylation at a subset of CTCF motifs proximal to HCC-associated genes can modify chromatin topological configuration, which in turn alter RNAPII-mediated chromatin interactions and lead to dysregulation of transcription. Together, the 3D-epigenomic analyses provide novel insights of multifaceted interplays involving genetics, epigenetics, and chromatin topology in HCC cells.

The 5-Hydroxymethylcytosine Landscape of Prostate Cancer

Analysis of DNA methylation is a valuable tool to understand disease progression and is increasingly being used to create diagnostic and prognostic clinical biomarkers. While conversion of cytosine to 5-methylcytosine (5mC) commonly results in transcriptional repression, further conversion to 5-hydroxymethylcytosine (5hmC) is associated with transcriptional activation. Here we perform the first study integrating whole-genome 5hmC with DNA, 5mC, and transcriptome sequencing in clinical samples of benign, localized, and advanced prostate cancer. 5hmC is shown to mark activation of cancer drivers and downstream targets. Furthermore, 5hmC sequencing revealed profoundly altered cell states throughout the disease course, characterized by increased proliferation, oncogenic signaling, dedifferentiation, and lineage plasticity to neuroendocrine and gastrointestinal lineages. Finally, 5hmC sequencing of cell-free DNA from patients with metastatic disease proved useful as a prognostic biomarker able to identify an aggressive subtype of prostate cancer using the genes TOP2A and EZH2, previously only detectable by transcriptomic analysis of solid tumor biopsies. Overall, these findings reveal that 5hmC marks epigenomic activation in prostate cancer and identify hallmarks of prostate cancer progression with potential as biomarkers of aggressive disease.

SIGNIFICANCE

In prostate cancer, 5-hydroxymethylcytosine delineates oncogene activation and stage-specific cell states and can be analyzed in liquid biopsies to detect cancer phenotypes. See related article by Wu and Attard, p. 3880.

Wntless expression promotes lineage plasticity and is associated with neuroendocrine prostate cancer

Resistance to androgen receptor (AR) targeted therapies remains as the main reason for most prostate cancer related deaths. Lineage plasticity resulting in altered, treatment insensitive prostate tumor cell phenotypes such neuroendocrine differentiated prostate cancer is a common manifestation within resistant tumors upon AR-targeted therapies. The mechanisms responsible for lineage plasticity in prostate cancer remain incompletely understood. Here we demonstrate that the enzalutamide resistant MDVR cell line possesses lineage plastic characteristics associated with overexpression of the Wnt transporter Wntless (WLS). Furthermore, we present evidence that overexpression of WLS is common in varying cell line models of lineage plastic prostate cancer, is higher in neuroendocrine patient samples, and positively correlates with the neuroendocrine marker SYP in clinical data. Targeting WLS in lineage plastic cellular models reduces viability and represses lineage plasticity associated gene expression. Our study provides insight into the importance of WLS to the development of lethal resistant prostate cancer and provides a potential target for the treatment of advanced disease.

Fatty Acid Signaling Impacts Prostate Cancer Lineage Plasticity in an Autocrine and Paracrine Manner

Prostate cancer (PCa) affects an estimated 250,000 men every year and causes 34,000 deaths annually. A high-fat diet and obesity are associated with PCa progression and mortality. This study's premise was the novel observation of crosstalk between PCa epithelia and cancer-associated fibroblasts (CAF) in response to palmitate-mediated lineage plasticity. We found that cholesterol activated canonical Hedgehog (Hh) signaling by increasing cilium Gli activity in PCa cells, while palmitate activated Hh independent of Gli. Exogenous palmitate activated SOX2, a known mediator of lineage plasticity, in PCa cells cocultured with CAF. Stroma-derived Wnt5a was upregulated in CAF while cocultured with PCa cells and treated with palmitate. Wnt5a knockdown in CAF inhibited Hh and SOX2 expression in PCa cells from cocultures. These findings supported our proposed mechanism of a high-fat diet promoting Hh signaling-mediated transformation within the tumor microenvironment. SOX2 and Wnt5a expression were limited by the CD36 neutralizing antibody. Mice xenografted with PCa epithelia and CAF tumors were fed a high-fat diet, leading to elevated SOX2 expression and lineage plasticity reprogramming compared to mice fed an isocaloric rodent diet. CD36 inhibition with enzalutamide elevated apoptosis by TUNEL, but limited proliferation and SOX2 expression compared to enzalutamide alone. This study revealed a mechanism for a high-fat diet to affect prostate cancer progression. We found that saturated fat induced lineage plasticity reprogramming of PCa by interaction with CAF through Wnt5a and Hh signaling.

Chromatin profiles classify castration-resistant prostate cancers suggesting therapeutic targets

In castration-resistant prostate cancer (CRPC), the loss of androgen receptor (AR) dependence leads to clinically aggressive tumors with few therapeutic options. We used ATAC-seq (assay for transposase-accessible chromatin sequencing), RNA-seq, and DNA sequencing to investigate 22 organoids, six patient-derived xenografts, and 12 cell lines. We identified the well-characterized AR-dependent and neuroendocrine subtypes, as well as two AR-negative/low groups: a Wnt-dependent subtype, and a stem cell-like (SCL) subtype driven by activator protein-1 (AP-1) transcription factors. We used transcriptomic signatures to classify 366 patients, which showed that SCL is the second most common subtype of CRPC after AR-dependent. Our data suggest that AP-1 interacts with the YAP/TAZ and TEAD proteins to maintain subtype-specific chromatin accessibility and transcriptomic landscapes in this group. Together, this molecular classification reveals drug targets and can potentially guide therapeutic decisions.

Response to supraphysiological testosterone is predicted by a distinct androgen receptor cistrome

The androgen receptor (AR) is a master transcription factor that regulates prostate cancer (PC) development and progression. Inhibition of AR signaling by androgen deprivation is the first-line therapy with initial efficacy for advanced and recurrent PC. Paradoxically, supraphysiological levels of testosterone (SPT) also inhibit PC progression. However, as with any therapy, not all patients show a therapeutic benefit, and responses differ widely in magnitude and duration. In this study, we evaluated whether differences in the AR cistrome before treatment can distinguish between SPT-responding (R) and -nonresponding (NR) tumors. We provide the first preclinical evidence to our knowledge that SPT-R tumors exhibit a distinct AR cistrome when compared with SPT-NR tumors, indicating a differential biological role of the AR. We applied an integrated analysis of ChIP-Seq and RNA-Seq to the pretreatment tumors and identified an SPT-R signature that distinguishes R and NR tumors. Because transcriptomes of SPT-treated clinical specimens are not available, we interrogated available castration-resistant PC (CRPC) transcriptomes and showed that the SPT-R signature is associated with improved survival and has the potential to identify patients who would respond to SPT. These findings provide an opportunity to identify the subset of patients with CRPC who would benefit from SPT therapy.

PGC1 alpha coactivates ERG fusion to drive antioxidant target genes under metabolic stress

The presence of ERG gene fusion; from developing prostatic intraepithelial neoplasia (PIN) lesions to hormone resistant high grade prostate cancer (PCa) dictates disease progression, altered androgen metabolism, proliferation and metastasis^1–3. ERG driven transcriptional landscape may provide pro-tumorigenic cues in overcoming various strains like hypoxia, nutrient deprivation, inflammation and oxidative stress. However, insights on the androgen independent regulation and function of ERG during stress are limited. Here, we identify PGC1α as a coactivator of ERG fusion under various metabolic stress. Deacetylase SIRT1 is necessary for PGC1α-ERG interaction and function. We reveal that ERG drives the expression of antioxidant genes; SOD1 and TXN, benefitting PCa growth. We observe increased expression of these antioxidant genes in patients with high ERG expression correlates with poor survival. Inhibition of PGC1α-ERG axis driven transcriptional program results in apoptosis and reduction in PCa xenografts. Here we report a function of ERG under metabolic stress which warrants further studies as a therapeutic target for ERG fusion positive PCa.

PGC1α acts as a co-activator of the ERG transcription factor during metabolic stress resulting in antioxidant functionsand inhibition of the PGC1α-ERG driven transcriptional program reduces prostate cancer growth by inducing ROS mediated apoptosis.

Past, Current, and Future Strategies to Target ERG Fusion-Positive Prostate Cancer

Simple Summary

In addition to its role in development and in the vascular and hematopoietic systems, ERG plays a central role in prostate cancer. Approximately 40–50% of prostate cancer cases are characterized by ERG gene fusions, which lead to ERG overexpression. Importantly, inhibition of ERG activity in prostate cancer cells decreases their viability. Therefore, inhibiting ERG might represent an important step to improve treatment efficacy for patients with ERG-positive prostate tumors. Here, we summarize the attempts made over the past years to repress ERG activity, the current use of ERG fusion detection and the strategies that might be utilized in the future to treat ERG fusion-positive tumors.

Abstract

The ETS family member ERG is a transcription factor with physiological roles during development and in the vascular and hematopoietic systems. ERG oncogenic activity characterizes several malignancies, including Ewing’s sarcoma, leukemia and prostate cancer (PCa). In PCa, ERG rearrangements with androgen-regulated genes—mostly TMPRSS2—characterize a large subset of patients across disease progression and result in androgen receptor (AR)-mediated overexpression of ERG in the prostate cells. Importantly, PCa cells overexpressing ERG are dependent on ERG activity for survival, further highlighting its therapeutic potential. Here, we review the current understanding of the role of ERG and its partners in PCa. We discuss the strategies developed in recent years to inhibit ERG activity, the current therapeutic utility of ERG fusion detection in PCa patients, and the possible future approaches to target ERG fusion-positive tumors.

Hepatocyte nuclear factor 4 gamma (HNF4G) is correlated with poor prognosis and promotes tumor cell growth by inhibiting caspase-dependent intrinsic apoptosis in colorectal cancer

The hepatocyte nuclear factor 4 gamma (HNF4G), a member of orphan nuclear receptors, is up-regulated and functions as an oncoprotein in a variety of tumors. Recent advances in understanding the biologic function and action mechanism of HNF4G in colorectal cancer (CRC) have not been fully elucidated. In the present study, we observed that HNF4G expression levels were significantly increased in CRC tissues compared with adjacent normal tissues, and HNF4G overexpression correlated with worse prognosis in colorectal cancer. Transfection with a small interference RNA (siRNA) targeting HNF4G in HCT116 and SW480 CRC cell lines significantly inhibited cell proliferation and promoted apoptosis in vitro. In contrast, overexpression of HNF4G increased cell proliferation and decreased the percentage of apoptotic cells. Moreover, we discovered that HNF4G was involved in CRC cell apoptosis via the caspase-dependent intrinsic pathway. Finally, knockdown of HNF4G expression led to attenuated colorectal cancer growth and promoted apoptosis in a xenograft mouse model. Collectively, these results indicate that HNF4G exerts as an oncogenic role in colorectal cancer and provides a potential therapeutic target.

Pseudogene RPL32P3 regulates the blood-tumor barrier permeability via the YBX2/HNF4G axis

The existence of the blood–tumor barrier (BTB) severely hinders the transport of anti-tumor drugs to brain tumor tissues. Selectively opening BTB is of great significance to improve the chemotherapy effect of glioma. Pseudogenes have been recognized as important regulators in various biologic processes. In this study, we identified that ribosomal protein L32 pseudogene 3 (RPL32P3) was highly expressed in glioma-exposed endothelial cells (GECs). Knockdown of RPL32P3 decreased the expression of tight junction-related proteins (TJPs) and increased BTB permeability. Subsequent analysis of the underlying mechanism indicated that RPL32P3 recruited lysine methyltransferase 2 A (KMT2A) to the Y-box binding protein 2 (YBX2) promoter region and mediated H3K4me3 to promote YBX2 transcription. Highly expressed YBX2 bound and stabilized hepatocyte nuclear factor 4 gamma (HNF4G) mRNA. Highly expressed HNF4G directly bound to the promoters of TJPs ZO-1, occludin and claudin-5 to promote their transcriptional activities and regulated BTB permeability. The simultaneous knockdown of RPL32P3, YBX2, and HNF4G combined with doxorubicin (DOX) increased the apoptosis of glioma cells. In conclusion, the current study indicated that RPL32P3 knockdown increased BTB permeability through the YBX2/HNF4G pathway. These findings may provide new targets for the comprehensive treatment of glioma.

Opposing transcriptional programs of KLF5 and AR emerge during therapy for advanced prostate cancer

Endocrine therapies for prostate cancer inhibit the androgen receptor (AR) transcription factor. In most cases, AR activity resumes during therapy and drives progression to castration-resistant prostate cancer (CRPC). However, therapy can also promote lineage plasticity and select for AR-independent phenotypes that are uniformly lethal. Here, we demonstrate the stem cell transcription factor Krüppel-like factor 5 (KLF5) is low or absent in prostate cancers prior to endocrine therapy, but induced in a subset of CRPC, including CRPC displaying lineage plasticity. KLF5 and AR physically interact on chromatin and drive opposing transcriptional programs, with KLF5 promoting cellular migration, anchorage-independent growth, and basal epithelial cell phenotypes. We identify ERBB2 as a point of transcriptional convergence displaying activation by KLF5 and repression by AR. ERBB2 inhibitors preferentially block KLF5-driven oncogenic phenotypes. These findings implicate KLF5 as an oncogene that can be upregulated in CRPC to oppose AR activities and promote lineage plasticity.

Allele-specific epigenetic activity in prostate cancer and normal prostate tissue implicates prostate cancer risk mechanisms

Genome-wide association studies (GWASs) of prostate cancer have identified >250 significant risk loci, but the causal variants and mechanisms for these loci remain largely unknown. Here, we sought to identify and characterize risk-harboring regulatory elements by integrating epigenomes from primary prostate tumor and normal tissues of 27 individuals across the H3K27ac, H3K4me3, and H3K4me2 histone marks and FOXA1 and HOXB13 transcription factors. We identified 7,371 peaks with significant allele specificity (allele-specific quantitative trait locus [asQTL] peaks). Showcasing their relevance to prostate cancer risk, H3K27ac T-asQTL peaks were the single annotation most enriched for prostate cancer GWAS heritability (40×), significantly higher than corresponding non-asQTL H3K27ac peaks (14×) or coding regions (14×). Surprisingly, fine-mapped GWAS risk variants were most significantly enriched for asQTL peaks observed in tumors, including asQTL peaks that were differentially imbalanced with respect to tumor-normal states. These data pinpointed putative causal regulatory elements at 20 GWAS loci, of which 11 were detected only in the tumor samples. More broadly, tumor-specific asQTLs were enriched for expression QTLs in benign tissues as well as accessible regions found in stem cells, supporting a hypothesis where some germline variants become reactivated during or after transformation and can be captured by epigenomic profiling of the tumor. Our study demonstrates the power of allele specificity in chromatin signals to uncover GWAS mechanisms, highlights the relevance of tumor-specific regulation in the context of cancer risk, and prioritizes multiple loci for experimental follow-up.

Transcriptional network involving ERG and AR orchestrates Distal-less homeobox-1 mediated prostate cancer progression

Distal-less homeobox-1 (DLX1) is a well-established non-invasive biomarker for prostate cancer (PCa) diagnosis, however, its mechanistic underpinnings in disease pathobiology are not known. Here, we reveal the oncogenic role of DLX1 and show that abrogating its function leads to reduced tumorigenesis and metastases. We observed that ~60% of advanced-stage and metastatic patients display higher DLX1 levels. Moreover, ~96% of TMPRSS2-ERG fusion-positive and ~70% of androgen receptor (AR)-positive patients show elevated DLX1, associated with aggressive disease and poor survival. Mechanistically, ERG coordinates with enhancer-bound AR and FOXA1 to drive transcriptional upregulation of DLX1 in ERG-positive background. However, in ERG-negative context, AR/AR-V7 and FOXA1 suffice to upregulate DLX1. Notably, inhibiting ERG/AR-mediated DLX1 transcription using BET inhibitor (BETi) or/and anti-androgen drugs reduce its expression and downstream oncogenic effects. Conclusively, this study establishes DLX1 as a direct-target of ERG/AR with an oncogenic role and demonstrates the clinical significance of BETi and anti-androgens for DLX1-positive patients.

Mapping of m6A and Its Regulatory Targets in Prostate Cancer Reveals a METTL3-Low Induction of Therapy Resistance

Recent evidence has highlighted the role of N 6-methyladenosine (m6A) in the regulation of mRNA expression, stability, and translation, supporting a potential role for posttranscriptional regulation mediated by m6A in cancer. Here, we explore prostate cancer as an exemplar and demonstrate that low levels of N 6-adenosine-methyltransferase (METTL3) is associated with advanced metastatic disease. To investigate this relationship, we generated the first prostate m6A maps, and further examined how METTL3 regulates expression at the level of transcription, translation, and protein. Significantly, transcripts encoding extracellular matrix proteins are consistently upregulated with METTL3 knockdown. We also examined the relationship between METTL3 and androgen signaling and discovered the upregulation of a hepatocyte nuclear factor-driven gene signature that is associated with therapy resistance in prostate cancer. Significantly, METTL3 knockdown rendered the cells resistant to androgen receptor antagonists via an androgen receptor-independent mechanism driven by the upregulation of nuclear receptor NR5A2/LRH-1. IMPLICATIONS: These findings implicate changes in m6A as a mechanism for therapy resistance in metastatic prostate cancer.

Elevated HNF1A expression promotes radiation-resistance via driving PI3K/AKT signaling pathway in esophageal squamous cell carcinoma cells

Purpose: Radiotherapy is a major modality for treatment of local advanced esophageal squamous cell carcinoma (ESCC). Hepatocyte nuclear factor 1-alpha (HNF1A) is involved in regulation of tumor cell proliferation, apoptosis, cycle distribution, invasion metastasis and chemical resistance. The aim of this study was to investigate the effect of HNF1A on radiosensitivity of ESCC cells. Methods: In our study, HNF1A expression was verified from GEPIA in multiple types of cancer. The prognostic value of HNF1A in ESCC was obtained by TCGA database. In addition, the expression of HNF1A in ESCC cell lines was verified by western blot. Subsequently, lentiviruses were used to construct HNF1A overexpressed cell lines TE1 and KYSE150.Then, the roles of HNF1A on cell proliferation, invasion, apoptosis, cell cycle distribution and radiosensitivity were verified. Furthermore, the relationship between HNF1A and γH2AX were determined by western blot and immunofluorescence. We also detected the expression changes of key factors in PI3K/AKT pathway after overexpression of HNF1A. Results: The results showed that the overexpression of HNF1A promoted cell proliferation and invasion with or without irradiation (IR), and potently radiation-resistance ESCC cells with a sensitization enhancement ratio (SER) of 0.76 and 0.87. In addition, HNF1A regulated Cyclin D1 and CDK4 proteins to promote the transition from radiation-induced G0/G1 phase arrest to S phase, and coordinated BAX and BCL2 proteins to reduce the occurrence of radiation-induced apoptosis. It was worth noting that HNF1A might be involved in radiation-induced DNA damage repair by regulating γH2AX though PI3K/AKT signal pathway. Conclusion: Our study preliminarily suggested that HNF1A was associated with the progression and radiosensitivity of ESCC cells, and it might reduce the radiosensitivity of ESCC cells by promoting cell proliferation, releasing G0/G1 phase arrest, reducing apoptosis, and regulating the expression of γH2AX protein though driving PI3K/AKT signal pathway.

microRNA-320b suppresses HNF4G and IGF2BP2 expression to inhibit angiogenesis and tumor growth of lung cancer

We examined the effect of microRNA-320b (miR-320b) on tumor growth and angiogenesis in lung cancer and also determined its downstream molecular mechanisms. Lung cancer tissues and adjacent non-cancerous tissues were collected from 66 patients with lung cancer. miR-320b expression was experimentally determined to be expressed at low level in cancer tissues. The results of gain-of-function experiments suggested that miR-320b overexpression suppressed cancer cell invasion, tube formation, tumor volume and angiogenesis in xenografted nude mice. Hepatocyte nuclear factor 4 gamma (HNF4G) was identified as a target of miR-320b based on in silico analysis. Dual-luciferase reporter gene assays further identified the binding relationship between HNF4G and miR-320b. Lung cancer tissues exhibited increased expression of HNF4G and insulin-like growth factor 2 mRNA-binding protein 2 (IGF2BP2). Meanwhile, HNF4G knockdown suppressed IGF2BP2 expression, thereby repressing cancer cell invasion and tube formation. Furthermore, IGF2BP2 modified m6A to increase the expression of thymidine kinase 1 (TK1), thus promoting angiogenesis. In nude mice, restoration of TK1 reversed the suppressive effect of miR-320b overexpression on tumor growth rate and CD31 expression. In conclusion, miR-320b suppresses lung cancer growth and angiogenesis by inhibiting HNF4G, IGF2BP2 and TK1.

Metformin inhibits pancreatic cancer metastasis caused by SMAD4 deficiency and consequent HNF4G upregulation

Pancreatic ductal adenocarcinoma (PDAC) has poor prognosis due to limited therapeutic options. This study examines the roles of genome-wide association study identified PDAC-associated genes as therapeutic targets. We have identified HNF4G gene whose silencing most effectively repressed PDAC cell invasiveness. HNF4G overexpression is induced by the deficiency of transcriptional factor and tumor suppressor SMAD4. Increased HNF4G are correlated with SMAD4 deficiency in PDAC tumor samples and associated with metastasis and poor survival time in xenograft animal model and in patients with PDAC (log-rank P = 0.036; HR = 1.60, 95% CI = 1.03-2.47). We have found that Metformin suppresses HNF4G activity via AMPK-mediated phosphorylation-coupled ubiquitination degradation and inhibits in vitro invasion and in vivo metastasis of PDAC cells with SMAD4 deficiency. Furthermore, Metformin treatment significantly improve clinical outcomes and survival in patients with SMAD4-deficient PDAC (log-rank P = 0.022; HR = 0.31, 95% CI = 0.14-0.68) but not in patients with SMAD4-normal PDAC. Pathway analysis shows that HNF4G may act in PDAC through the cell-cell junction pathway. These results indicate that SMAD4 deficiency-induced overexpression of HNF4G plays a critical oncogenic role in PDAC progression and metastasis but may form a druggable target for Metformin treatment.

Transcriptional programmes underlying cellular identity and microbial responsiveness in the intestinal epithelium

The intestinal epithelium serves the unique and critical function of harvesting dietary nutrients, while simultaneously acting as a cellular barrier separating tissues from the luminal environment and gut microbial ecosystem. Two salient features of the intestinal epithelium enable it to perform these complex functions. First, cells within the intestinal epithelium achieve a wide range of specialized identities, including different cell types and distinct anterior-posterior patterning along the intestine. Second, intestinal epithelial cells are sensitive and responsive to the dynamic milieu of dietary nutrients, xenobiotics and microorganisms encountered in the intestinal luminal environment. These diverse identities and responsiveness of intestinal epithelial cells are achieved in part through the differential transcription of genes encoded in their shared genome. Here, we review insights from mice and other vertebrate models into the transcriptional regulatory mechanisms underlying intestinal epithelial identity and microbial responsiveness, including DNA methylation, chromatin accessibility, histone modifications and transcription factors. These studies are revealing that most transcription factors involved in intestinal epithelial identity also respond to changes in the microbiota, raising both opportunities and challenges to discern the underlying integrative transcriptional regulatory networks.

Machine Learning Identifies Robust Matrisome Markers and Regulatory Mechanisms in Cancer

The expression and regulation of matrisome genes-the ensemble of extracellular matrix, ECM, ECM-associated proteins and regulators as well as cytokines, chemokines and growth factors-is of paramount importance for many biological processes and signals within the tumor microenvironment. The availability of large and diverse multi-omics data enables mapping and understanding of the regulatory circuitry governing the tumor matrisome to an unprecedented level, though such a volume of information requires robust approaches to data analysis and integration. In this study, we show that combining Pan-Cancer expression data from The Cancer Genome Atlas (TCGA) with genomics, epigenomics and microenvironmental features from TCGA and other sources enables the identification of "landmark" matrisome genes and machine learning-based reconstruction of their regulatory networks in 74 clinical and molecular subtypes of human cancers and approx. 6700 patients. These results, enriched for prognostic genes and cross-validated markers at the protein level, unravel the role of genetic and epigenetic programs in governing the tumor matrisome and allow the prioritization of tumor-specific matrisome genes (and their regulators) for the development of novel therapeutic approaches.

Control of Cell Identity by the Nuclear Receptor HNF4 in Organ Pathophysiology

Hepatocyte Nuclear Factor 4 (HNF4) is a transcription factor (TF) belonging to the nuclear receptor family whose expression and activities are restricted to a limited number of organs including the liver and gastrointestinal tract. In this review, we present robust evidence pointing to HNF4 as a master regulator of cellular differentiation during development and a safekeeper of acquired cell identity in adult organs. Importantly, we discuss that transient loss of HNF4 may represent a protective mechanism upon acute organ injury, while prolonged impairment of HNF4 activities could contribute to organ dysfunction. In this context, we describe in detail mechanisms involved in the pathophysiological control of cell identity by HNF4, including how HNF4 works as part of cell-specific TF networks and how its expression/activities are disrupted in injured organs.

Patient-specific multi-omics models and the application in personalized combination therapy

The rapid advancement of high-throughput technologies and sharp decrease in cost have opened up the possibility to generate large amount of multi-omics data on an individual basis. The development of high-throughput -omics, including genomics, epigenomics, transcriptomics, proteomics, metabolomics and microbiomics, enables the application of multi-omics technologies in the clinical settings. Combination therapy, defined as disease treatment with two or more drugs to achieve efficacy with lower doses or lower drug toxicity, is the basis for the care of diseases like cancer. Patient-specific multi-omics data integration can help the identification and development of combination therapies. In this review, we provide an overview of different -omics platforms, and discuss the methods for multi-omics, high-throughput, data integration, personalized combination therapy.

dbInDel: a database of enhancer-associated insertion and deletion variants by analysis of H3K27ac ChIP-Seq

SUMMARY

Cancer hallmarks rely on its specific transcriptional programs, which are dysregulated by multiple mechanisms, including genomic aberrations in the DNA regulatory regions. Genome-wide association studies have shown many variants are found within putative enhancer elements. To provide insights into the regulatory role of enhancer-associated non-coding variants in cancer epigenome, and to facilitate the identification of functional non-coding mutations, we present dbInDel, a database where we have comprehensively analyzed enhancer-associated insertion and deletion variants for both human and murine samples using ChIP-Seq data. Moreover, we provide the identification and visualization of upstream TF binding motifs in InDel-containing enhancers. Downstream target genes are also predicted and analyzed in the context of cancer biology. The dbInDel database promotes the investigation of functional contributions of non-coding variants in cancer epigenome.

AVAILABILITY AND IMPLEMENTATION

The database, dbInDel, can be accessed from http://enhancer-indel.cam-su.org/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.

Androgen deprivation upregulates SPINK1 expression and potentiates cellular plasticity in prostate cancer

Emergence of an aggressive androgen receptor (AR)-independent neuroendocrine prostate cancer (NEPC) after androgen-deprivation therapy (ADT) is well-known. Nevertheless, the majority of advanced-stage prostate cancer patients, including those with SPINK1-positive subtype, are treated with AR-antagonists. Here, we show AR and its corepressor, REST, function as transcriptional-repressors of SPINK1, and AR-antagonists alleviate this repression leading to SPINK1 upregulation. Increased SOX2 expression during NE-transdifferentiation transactivates SPINK1, a critical-player for maintenance of NE-phenotype. SPINK1 elicits epithelial-mesenchymal-transition, stemness and cellular-plasticity. Conversely, pharmacological Casein Kinase-1 inhibition stabilizes REST, which in cooperation with AR causes SPINK1 transcriptional-repression and impedes SPINK1-mediated oncogenesis. Elevated levels of SPINK1 and NEPC markers are observed in the tumors of AR-antagonists treated mice, and in a subset of NEPC patients, implicating a plausible role of SPINK1 in treatment-related NEPC. Collectively, our findings provide an explanation for the paradoxical clinical-outcomes after ADT, possibly due to SPINK1 upregulation, and offers a strategy for adjuvant therapies.

Gene set based systematic analysis of prostate cancer and its subtypes

Aim: A gene set based systematic analysis strategy is used to investigate prostate tumors and its subclusters with focuses on similarities and differences of biological functions. Results: Dysregulation of methylation status, as well as RAS/RAF/ERK and PI3K-ATK signaling pathways, were found to be the most dramatic changes during prostate cancer tumorigenesis. Besides, neural and inflammation microenvironment is also significantly divergent between tumor and adjacent tissues. Insights of subclasses within prostate tumor cohorts revealed four different clusters with distinct gene expression patterns. We found that samples are mainly clustered by immune environments and proliferation traits. Conclusion: The findings of this article may help to advance the progress of identifying better diagnosis biomarkers and therapeutic targets.

The Role of Lineage Plasticity in Prostate Cancer Therapy Resistance

Lineage plasticity has emerged as an important mechanism of treatment resistance in prostate cancer. Treatment-refractory prostate cancers are increasingly associated with loss of luminal prostate markers, and in many cases induction of developmental programs, stem cell-like phenotypes, and neuroendocrine/neuronal features. Clinically, lineage plasticity may manifest as low PSA progression, resistance to androgen receptor (AR) pathway inhibitors, and sometimes small cell/neuroendocrine pathologic features observed on metastatic biopsy. This mechanism is not restricted to prostate cancer as other malignancies also demonstrate lineage plasticity during resistance to targeted therapies. At present, there is no established therapeutic approach for patients with advanced prostate cancer developing lineage plasticity or small cell neuroendocrine prostate cancer (NEPC) due to knowledge gaps in the underlying biology. Few clinical trials address questions in this space, and the outlook for patients remains poor. To move forward, urgently needed are: (i) a fundamental understanding of how lineage plasticity occurs and how it can best be defined; (ii) the temporal contribution and cooperation of emerging drivers; (iii) preclinical models that recapitulate biology of the disease and the recognized phenotypes; (iv) identification of therapeutic targets; and (v) novel trial designs dedicated to the entity as it is defined. This Perspective represents a consensus arising from the NCI Workshop on Lineage Plasticity and Androgen Receptor-Independent Prostate Cancer. We focus on the critical questions underlying lineage plasticity and AR-independent prostate cancer, outline knowledge and resource gaps, and identify strategies to facilitate future collaborative clinical translational and basic studies in this space.

Cellular and Molecular Mechanisms Underlying Prostate Cancer Development: Therapeutic Implications

Prostate cancer is the most frequent nonskin cancer and second most common cause of cancer-related deaths in man. Prostate cancer is a clinically heterogeneous disease with many patients exhibiting an aggressive disease with progression, metastasis, and other patients showing an indolent disease with low tendency to progression. Three stages of development of human prostate tumors have been identified: intraepithelial neoplasia, adenocarcinoma androgen-dependent, and adenocarcinoma androgen-independent or castration-resistant. Advances in molecular technologies have provided a very rapid progress in our understanding of the genomic events responsible for the initial development and progression of prostate cancer. These studies have shown that prostate cancer genome displays a relatively low mutation rate compared with other cancers and few chromosomal loss or gains. The ensemble of these molecular studies has led to suggest the existence of two main molecular groups of prostate cancers: one characterized by the presence of ERG rearrangements (~50% of prostate cancers harbor recurrent gene fusions involving ETS transcription factors, fusing the 5' untranslated region of the androgen-regulated gene TMPRSS2 to nearly the coding sequence of the ETS family transcription factor ERG) and features of chemoplexy (complex gene rearrangements developing from a coordinated and simultaneous molecular event), and a second one characterized by the absence of ERG rearrangements and by the frequent mutations in the E3 ubiquitin ligase adapter SPOP and/or deletion of CDH1, a chromatin remodeling factor, and interchromosomal rearrangements and SPOP mutations are early events during prostate cancer development. During disease progression, genomic and epigenomic abnormalities accrued and converged on prostate cancer pathways, leading to a highly heterogeneous transcriptomic landscape, characterized by a hyperactive androgen receptor signaling axis.

Molecular Underpinnings Governing Genetic Complexity of ETS-Fusion-Negative Prostate Cancer

Inter- and intra-patient molecular heterogeneity of primary and metastatic prostate cancer (PCa) confers variable clinical outcome and poses a formidable challenge in disease management. High-throughput integrative genomics and functional approaches have untangled the complexity involved in this disease and revealed a spectrum of diverse aberrations prevalent in various molecular subtypes, including ETS fusion negative. Emerging evidence indicates that SPINK1 upregulation, mutations in epigenetic regulators or chromatin modifiers, and SPOP are associated with the ETS-fusion negative subtype. Additionally, patients with defects in a DNA-repair pathway respond to poly-(ADP-ribose)-polymerase (PARP) inhibition therapies. Furthermore, a new class of immunogenic subtype defined by CDK12 biallelic loss has also been identified in ETS-fusion-negative cases. This review focuses on the emerging molecular underpinnings driving key oncogenic aberrations and advancements in therapeutic strategies of this disease.

Epigenetic Silencing of miRNA-338-5p and miRNA-421 Drives SPINK1-Positive Prostate Cancer

PURPOSE

Serine peptidase inhibitor, Kazal type-1 (SPINK1) overexpression defines the second most recurrent and aggressive prostate cancer subtype. However, the underlying molecular mechanism and pathobiology of SPINK1 in prostate cancer remains largely unknown.

EXPERIMENTAL DESIGN

miRNA prediction tools were employed to examine the SPINK1-3'UTR for miRNA binding. Luciferase reporter assays were performed to confirm the SPINK1-3'UTR binding of shortlisted miR-338-5p/miR-421. Furthermore, miR-338-5p/-421-overexpressing cancer cells (SPINK1-positive) were evaluated for oncogenic properties using cell-based functional assays and a mouse xenograft model. Global gene expression profiling was performed to unravel the biological pathways altered by miR-338-5p/-421. IHC and RNA in situ hybridization were carried out on prostate cancer patients' tissue microarray for SPINK1 and EZH2 expression, respectively. Chromatin immunoprecipitation assay was performed to examine EZH2 occupancy on the miR-338-5p/-421-regulatory regions. Bisulfite sequencing and methylated DNA immunoprecipitation were performed on prostate cancer cell lines and patients' specimens.

RESULTS

We established a critical role of miRNA-338-5p/-421 in posttranscriptional regulation of SPINK1. Ectopic expression of miRNA-338-5p/-421 in SPINK1-positive cells abrogates oncogenic properties including cell-cycle progression, stemness, and drug resistance, and shows reduced tumor burden and distant metastases in a mouse model. Importantly, we show that patients with SPINK1-positive prostate cancer exhibit increased EZH2 expression, suggesting its role in epigenetic silencing of miRNA-338-5p/-421. Furthermore, presence of CpG dinucleotide DNA methylation marks on the regulatory regions of miR-338-5p/-421 in SPINK1-positive prostate cancer cells and patients' specimens confirms epigenetic silencing.

CONCLUSIONS

Our findings revealed that miRNA-338-5p/-421 are epigenetically silenced in SPINK1-positive prostate cancer, although restoring the expression of these miRNAs using epigenetic drugs or synthetic mimics could abrogate SPINK1-mediated oncogenesis.See related commentary by Bjartell, p. 2679.