DNA hypomethylation silences anti-tumor immune genes in early prostate cancer and CTCs.

Cancer is characterized by hypomethylation-associated silencing of large chromatin domains, whose contribution to tumorigenesis is uncertain. Through high-resolution genome-wide single-cell DNA methylation sequencing, we identify 40 core domains that are uniformly hypomethylated from the earliest detectable stages of prostate malignancy through metastatic circulating tumor cells (CTCs). Nested among these repressive domains are smaller loci with preserved methylation that escape silencing and are enriched for cell proliferation genes. Transcriptionally silenced genes within the core hypomethylated domains are enriched for immune-related genes; prominent among these is a single gene cluster harboring all five CD1 genes that present lipid antigens to NKT cells and four IFI16-related interferon-inducible genes implicated in innate immunity. The re-expression of CD1 or IFI16 murine orthologs in immuno-competent mice abrogates tumorigenesis, accompanied by the activation of anti-tumor immunity. Thus, early epigenetic changes may shape tumorigenesis, targeting co-located genes within defined chromosomal loci. Hypomethylation domains are detectable in blood specimens enriched for CTCs.

Widespread and Functional RNA Circularization in Localized Prostate Cancer.

The cancer transcriptome is remarkably complex, including low-abundance transcripts, many not polyadenylated. To fully characterize the transcriptome of localized prostate cancer, we performed ultra-deep total RNA-seq on 144 tumors with rich clinical annotation. This revealed a linear transcriptomic subtype associated with the aggressive intraductal carcinoma sub-histology and a fusion profile that differentiates localized from metastatic disease. Analysis of back-splicing events showed widespread RNA circularization, with the average tumor expressing 7,232 circular RNAs (circRNAs). The degree of circRNA production was correlated to disease progression in multiple patient cohorts. Loss-of-function screening identified 11.3% of highly abundant circRNAs as essential for cell proliferation; for ∼90% of these, their parental linear transcripts were not essential. Individual circRNAs can have distinct functions, with circCSNK1G3 promoting cell growth by interacting with miR-181. These data advocate for adoption of ultra-deep RNA-seq without poly-A selection to interrogate both linear and circular transcriptomes.

Differences in Tumor Microenvironment Dictate T Helper Lineage Polarization and Response to Immune Checkpoint Therapy.

Immune checkpoint therapy (ICT) shows encouraging results in a subset of patients with metastatic castration-resistant prostate cancer (mCRPC) but still elicits a sub-optimal response among those with bone metastases. Analysis of patients' bone marrow samples revealed increased Th17 instead of Th1 subsets after ICT. To further evaluate the different tumor microenvironments, we injected mice with prostate tumor cells either subcutaneously or intraosseously. ICT in the subcutaneous CRPC model significantly increases intra-tumoral Th1 subsets and improves survival. However, ICT fails to elicit an anti-tumor response in the bone CRPC model despite an increase in the intra-tumoral CD4 T cells, which are polarized to Th17 rather than Th1 lineage. Mechanistically, tumors in the bone promote osteoclast-mediated bone resorption that releases TGF-ß, which restrains Th1 lineage development. Blocking TGF-ß along with ICT increases Th1 subsets and promotes clonal expansion of CD8 T cells and subsequent regression of bone CRPC and improves survival.

The Evolutionary Landscape of Localized Prostate Cancers Drives Clinical Aggression.

The majority of newly diagnosed prostate cancers are slow growing, with a long natural life history. Yet a subset can metastasize with lethal consequences. We reconstructed the phylogenies of 293 localized prostate tumors linked to clinical outcome data. Multiple subclones were detected in 59% of patients, and specific subclonal architectures associate with adverse clinicopathological features. Early tumor development is characterized by point mutations and deletions followed by later subclonal amplifications and changes in trinucleotide mutational signatures. Specific genes are selectively mutated prior to or following subclonal diversification, including MTOR, NKX3-1, and RB1. Patients with low-risk monoclonal tumors rarely relapse after primary therapy (7%), while those with high-risk polyclonal tumors frequently do (61%). The presence of multiple subclones in an index biopsy may be necessary, but not sufficient, for relapse of localized prostate cancer, suggesting that evolution-aware biomarkers should be studied in prospective studies of low-risk tumors suitable for active surveillance.

A Somatically Acquired Enhancer of the Androgen Receptor Is a Noncoding Driver in Advanced Prostate Cancer.

Increased androgen receptor (AR) activity drives therapeutic resistance in advanced prostate cancer. The most common resistance mechanism is amplification of this locus presumably targeting the AR gene. Here, we identify and characterize a somatically acquired AR enhancer located 650 kb centromeric to the AR. Systematic perturbation of this enhancer using genome editing decreased proliferation by suppressing AR levels. Insertion of an additional copy of this region sufficed to increase proliferation under low androgen conditions and to decrease sensitivity to enzalutamide. Epigenetic data generated in localized prostate tumors and benign specimens support the notion that this region is a developmental enhancer. Collectively, these observations underscore the importance of epigenomic profiling in primary specimens and the value of deploying genome editing to functionally characterize noncoding elements. More broadly, this work identifies a therapeutic vulnerability for targeting the AR and emphasizes the importance of regulatory elements as highly recurrent oncogenic drivers.

Inactivation of CDK12 Delineates a Distinct Immunogenic Class of Advanced Prostate Cancer.

Using integrative genomic analysis of 360 metastatic castration-resistant prostate cancer (mCRPC) samples, we identified a novel subtype of prostate cancer typified by biallelic loss of CDK12 that is mutually exclusive with tumors driven by DNA repair deficiency, ETS fusions, and SPOP mutations. CDK12 loss is enriched in mCRPC relative to clinically localized disease and characterized by focal tandem duplications (FTDs) that lead to increased gene fusions and marked differential gene expression. FTDs associated with CDK12 loss result in highly recurrent gains at loci of genes involved in the cell cycle and DNA replication. CDK12 mutant cases are baseline diploid and do not exhibit DNA mutational signatures linked to defects in homologous recombination. CDK12 mutant cases are associated with elevated neoantigen burden ensuing from fusion-induced chimeric open reading frames and increased tumor T cell infiltration/clonal expansion. CDK12 inactivation thereby defines a distinct class of mCRPC that may benefit from immune checkpoint immunotherapy.

Risk SNP-Mediated Promoter-Enhancer Switching Drives Prostate Cancer through lncRNA PCAT19.

The prostate cancer (PCa) risk-associated SNP rs11672691 is positively associated with aggressive disease at diagnosis. We showed that rs11672691 maps to the promoter of a short isoform of long noncoding RNA PCAT19 (PCAT19-short), which is in the third intron of the long isoform (PCAT19-long). The risk variant is associated with decreased and increased levels of PCAT19-short and PCAT19-long, respectively. Mechanistically, the risk SNP region is bifunctional with both promoter and enhancer activity. The risk variants of rs11672691 and its LD SNP rs887391 decrease binding of transcription factors NKX3.1 and YY1 to the promoter of PCAT19-short, resulting in weaker promoter but stronger enhancer activity that subsequently activates PCAT19-long. PCAT19-long interacts with HNRNPAB to activate a subset of cell-cycle genes associated with PCa progression, thereby promoting PCa tumor growth and metastasis. Taken together, these findings reveal a risk SNP-mediated promoter-enhancer switching mechanism underlying both initiation and progression of aggressive PCa.

Genomic Hallmarks and Structural Variation in Metastatic Prostate Cancer.

While mutations affecting protein-coding regions have been examined across many cancers, structural variants at the genome-wide level are still poorly defined. Through integrative deep whole-genome and -transcriptome analysis of 101 castration-resistant prostate cancer metastases (109X tumor/38X normal coverage), we identified structural variants altering critical regulators of tumorigenesis and progression not detectable by exome approaches. Notably, we observed amplification of an intergenic enhancer region 624 kb upstream of the androgen receptor (AR) in 81% of patients, correlating with increased AR expression. Tandem duplication hotspots also occur near MYC, in lncRNAs associated with post-translational MYC regulation. Classes of structural variations were linked to distinct DNA repair deficiencies, suggesting their etiology, including associations of CDK12 mutation with tandem duplications, TP53 inactivation with inverted rearrangements and chromothripsis, and BRCA2 inactivation with deletions. Together, these observations provide a comprehensive view of how structural variations affect critical regulators in metastatic prostate cancer.

Nuclear Pores Promote Lethal Prostate Cancer by Increasing POM121-Driven E2F1, MYC, and AR Nuclear Import.

Nuclear pore complexes (NPCs) regulate nuclear-cytoplasmic transport, transcription, and genome integrity in eukaryotic cells. However, their functional roles in cancer remain poorly understood. We interrogated the evolutionary transcriptomic landscape of NPC components, nucleoporins (Nups), from primary to advanced metastatic human prostate cancer (PC). Focused loss-of-function genetic screen of top-upregulated Nups in aggressive PC models identified POM121 as a key contributor to PC aggressiveness. Mechanistically, POM121 promoted PC progression by enhancing importin-dependent nuclear transport of key oncogenic (E2F1, MYC) and PC-specific (AR-GATA2) transcription factors, uncovering a pharmacologically targetable axis that, when inhibited, decreased tumor growth, restored standard therapy efficacy, and improved survival in patient-derived pre-clinical models. Our studies molecularly establish a role of NPCs in PC progression and give a rationale for NPC-regulated nuclear import targeting as a therapeutic strategy for lethal PC. These findings may have implications for understanding how NPC deregulation contributes to the pathogenesis of other tumor types.

Biology and Clinical Implications of the 19q13 Aggressive Prostate Cancer Susceptibility Locus.

Genome-wide association studies (GWAS) have identified rs11672691 at 19q13 associated with aggressive prostate cancer (PCa). Here, we independently confirmed the finding in a cohort of 2,738 PCa patients and discovered the biological mechanism underlying this association. We found an association of the aggressive PCa-associated allele G of rs11672691 with elevated transcript levels of two biologically plausible candidate genes, PCAT19 and CEACAM21, implicated in PCa cell growth and tumor progression. Mechanistically, rs11672691 resides in an enhancer element and alters the binding site of HOXA2, a novel oncogenic transcription factor with prognostic potential in PCa. Remarkably, CRISPR/Cas9-mediated single-nucleotide editing showed the direct effect of rs11672691 on PCAT19 and CEACAM21 expression and PCa cellular aggressive phenotype. Clinical data demonstrated synergistic effects of rs11672691 genotype and PCAT19/CEACAM21 gene expression on PCa prognosis. These results provide a plausible mechanism for rs11672691 associated with aggressive PCa and thus lay the ground work for translating this finding to the clinic.

Structural Alterations Driving Castration-Resistant Prostate Cancer Revealed by Linked-Read Genome Sequencing.

Nearly all prostate cancer deaths are from metastatic castration-resistant prostate cancer (mCRPC), but there have been few whole-genome sequencing (WGS) studies of this disease state. We performed linked-read WGS on 23 mCRPC biopsy specimens and analyzed cell-free DNA sequencing data from 86 patients with mCRPC. In addition to frequent rearrangements affecting known prostate cancer genes, we observed complex rearrangements of the AR locus in most cases. Unexpectedly, these rearrangements include highly recurrent tandem duplications involving an upstream enhancer of AR in 70%-87% of cases compared with <2% of primary prostate cancers. A subset of cases displayed AR or MYC enhancer duplication in the context of a genome-wide tandem duplicator phenotype associated with CDK12 inactivation. Our findings highlight the complex genomic structure of mCRPC, nominate alterations that may inform prostate cancer treatment, and suggest that additional recurrent events in the non-coding mCRPC genome remain to be discovered.

A Unifying Theory of Branching Morphogenesis.

The morphogenesis of branched organs remains a subject of abiding interest. Although much is known about the underlying signaling pathways, it remains unclear how macroscopic features of branched organs, including their size, network topology, and spatial patterning, are encoded. Here, we show that, in mouse mammary gland, kidney, and human prostate, these features can be explained quantitatively within a single unifying framework of branching and annihilating random walks. Based on quantitative analyses of large-scale organ reconstructions and proliferation kinetics measurements, we propose that morphogenesis follows from the proliferative activity of equipotent tips that stochastically branch and randomly explore their environment but compete neutrally for space, becoming proliferatively inactive when in proximity with neighboring ducts. These results show that complex branched epithelial structures develop as a self-organized process, reliant upon a strikingly simple but generic rule, without recourse to a rigid and deterministic sequence of genetically programmed events.

2'-O-methylation at internal sites on mRNA promotes mRNA stability.

2'-O-methylation (Nm) is a prominent RNA modification well known in noncoding RNAs and more recently also found at many mRNA internal sites. However, their function and base-resolution stoichiometry remain underexplored. Here, we investigate the transcriptome-wide effect of internal site Nm on mRNA stability. Combining nanopore sequencing with our developed machine learning method, NanoNm, we identify thousands of Nm sites on mRNAs with a single-base resolution. We observe a positive effect of FBL-mediated Nm modification on mRNA stability and expression level. Elevated FBL expression in cancer cells is associated with increased expression levels for 2'-O-methylated mRNAs of cancer pathways, implying the role of FBL in post-transcriptional regulation. Lastly, we find that FBL-mediated 2'-O-methylation connects to widespread 3' UTR shortening, a mechanism that globally increases RNA stability. Collectively, we demonstrate that FBL-mediated Nm modifications at mRNA internal sites regulate gene expression by enhancing mRNA stability.

Hyd/UBR5 defines a tumor suppressor pathway that links Polycomb repressive complex to regulated protein degradation in tissue growth control and tumorigenesis.

Tumor suppressor genes play critical roles in normal tissue homeostasis, and their dysregulation underlies human diseases including cancer. Besides human genetics, model organisms such as Drosophila have been instrumental in discovering tumor suppressor pathways that were subsequently shown to be highly relevant in human cancer. Here we show that hyperplastic disc (Hyd), one of the first tumor suppressors isolated genetically in Drosophila and encoding an E3 ubiquitin ligase with hitherto unknown substrates, and Lines (Lin), best known for its role in embryonic segmentation, define an obligatory tumor suppressor protein complex (Hyd-Lin) that targets the zinc finger-containing oncoprotein Bowl for ubiquitin-mediated degradation, with Lin functioning as a substrate adaptor to recruit Bowl to Hyd for ubiquitination. Interestingly, the activity of the Hyd-Lin complex is directly inhibited by a micropeptide encoded by another zinc finger gene, drumstick (drm), which functions as a pseudosubstrate by displacing Bowl from the Hyd-Lin complex, thus stabilizing Bowl. We further identify the epigenetic regulator Polycomb repressive complex1 (PRC1) as a critical upstream regulator of the Hyd-Lin-Bowl pathway by directly repressing the transcription of the micropeptide drm Consistent with these molecular studies, we show that genetic inactivation of Hyd, Lin, or PRC1 resulted in Bowl-dependent hyperplastic tissue overgrowth in vivo. We also provide evidence that the mammalian homologs of Hyd (UBR5, known to be recurrently dysregulated in various human cancers), Lin (LINS1), and Bowl (OSR1/2) constitute an analogous protein degradation pathway in human cells, and that OSR2 promotes prostate cancer tumorigenesis. Altogether, these findings define a previously unrecognized tumor suppressor pathway that links epigenetic program to regulated protein degradation in tissue growth control and tumorigenesis.

Minor intron splicing is critical for survival of lethal prostate cancer.

The evolutionarily conserved minor spliceosome (MiS) is required for protein expression of ∼714 minor intron-containing genes (MIGs) crucial for cell-cycle regulation, DNA repair, and MAP-kinase signaling. We explored the role of MIGs and MiS in cancer, taking prostate cancer (PCa) as an exemplar. Both androgen receptor signaling and elevated levels of U6atac, a MiS small nuclear RNA, regulate MiS activity, which is highest in advanced metastatic PCa. siU6atac-mediated MiS inhibition in PCa in vitro model systems resulted in aberrant minor intron splicing leading to cell-cycle G1 arrest. Small interfering RNA knocking down U6atac was ∼50% more efficient in lowering tumor burden in models of advanced therapy-resistant PCa compared with standard antiandrogen therapy. In lethal PCa, siU6atac disrupted the splicing of a crucial lineage dependency factor, the RE1-silencing factor (REST). Taken together, we have nominated MiS as a vulnerability for lethal PCa and potentially other cancers.

PLEKHS1 drives PI3Ks and remodels pathway homeostasis in PTEN-null prostate.

The PIP3/PI3K network is a central regulator of metabolism and is frequently activated in cancer, commonly by loss of the PIP3/PI(3,4)P2 phosphatase, PTEN. Despite huge research investment, the drivers of the PI3K network in normal tissues and how they adapt to overactivation are unclear. We find that in healthy mouse prostate PI3K activity is driven by RTK/IRS signaling and constrained by pathway feedback. In the absence of PTEN, the network is dramatically remodeled. A poorly understood YXXM- and PIP3/PI(3,4)P2-binding PH domain-containing adaptor, PLEKHS1, became the dominant activator and was required to sustain PIP3, AKT phosphorylation, and growth in PTEN-null prostate. This was because PLEKHS1 evaded pathway-feedback and experienced enhanced PI3K- and Src-family kinase-dependent phosphorylation of Y258XXM, eliciting PI3K activation. hPLEKHS1 mRNA and activating Y419 phosphorylation of hSrc correlated with PI3K pathway activity in human prostate cancers. We propose that in PTEN-null cells receptor-independent, Src-dependent tyrosine phosphorylation of PLEKHS1 creates positive feedback that escapes homeostasis, drives PIP3 signaling, and supports tumor progression.

Higher-order SPOP assembly reveals a basis for cancer mutant dysregulation.

The speckle-type POZ protein (SPOP) functions in the Cullin3-RING ubiquitin ligase (CRL3) as a receptor for the recognition of substrates involved in cell growth, survival, and signaling. SPOP mutations have been attributed to the development of many types of cancers, including prostate and endometrial cancers. Prostate cancer mutations localize in the substrate-binding site of the substrate recognition (MATH) domain and reduce or prevent binding. However, most endometrial cancer mutations are dispersed in seemingly inconspicuous solvent-exposed regions of SPOP, offering no clear basis for their cancer-causing and peculiar gain-of-function properties. Herein, we present the first structure of SPOP in its oligomeric form, uncovering several new interfaces important for SPOP self-assembly and normal function. Given that many previously unaccounted-for cancer mutations are localized in these newly identified interfaces, we uncover molecular mechanisms underlying dysregulation of SPOP function, with effects ranging from gross structural changes to enhanced self-association, and heightened stability and activity.

A lncRNA from the FTO locus acts as a suppressor of the m6A writer complex and p53 tumor suppression signaling.

N6-methyladenosine (m6A) of mRNAs modulated by the METTL3-METTL14-WTAP-RBM15 methyltransferase complex and m6A demethylases such as FTO play important roles in regulating mRNA stability, splicing, and translation. Here, we demonstrate that FTO-IT1 long noncoding RNA (lncRNA) was upregulated and positively correlated with poor survival of patients with wild-type p53-expressing prostate cancer (PCa). m6A RIP-seq analysis revealed that FTO-IT1 knockout increased mRNA m6A methylation of a subset of p53 transcriptional target genes (e.g., FAS, TP53INP1, and SESN2) and induced PCa cell cycle arrest and apoptosis. We further showed that FTO-IT1 directly binds RBM15 and inhibits RBM15 binding, m6A methylation, and stability of p53 target mRNAs. Therapeutic depletion of FTO-IT1 restored mRNA m6A level and expression of p53 target genes and inhibited PCa growth in mice. Our study identifies FTO-IT1 lncRNA as a bona fide suppressor of the m6A methyltransferase complex and p53 tumor suppression signaling and nominates FTO-IT1 as a potential therapeutic target of cancer.

Allosteric interactions prime androgen receptor dimerization and activation.

The androgen receptor (AR) is a nuclear receptor that governs gene expression programs required for prostate development and male phenotype maintenance. Advanced prostate cancers display AR hyperactivation and transcriptome expansion, in part, through AR amplification and interaction with oncoprotein cofactors. Despite its biological importance, how AR domains and cofactors cooperate to bind DNA has remained elusive. Using single-particle cryo-electron microscopy, we isolated three conformations of AR bound to DNA, showing that AR forms a non-obligate dimer, with the buried dimer interface utilized by ancestral steroid receptors repurposed to facilitate cooperative DNA binding. We identify novel allosteric surfaces which are compromised in androgen insensitivity syndrome and reinforced by AR's oncoprotein cofactor, ERG, and by DNA-binding motifs. Finally, we present evidence that this plastic dimer interface may have been adopted for transactivation at the expense of DNA binding. Our work highlights how fine-tuning AR's cooperative interactions translate to consequences in development and disease.

PALI1 promotes tumor growth through competitive recruitment of PRC2 to G9A-target chromatin for dual epigenetic silencing.

PALI1 is a newly identified accessory protein of the Polycomb repressive complex 2 (PRC2) that catalyzes H3K27 methylation. However, the roles of PALI1 in cancer are yet to be defined. Here, we report that PALI1 is upregulated in advanced prostate cancer (PCa) and competes with JARID2 for binding to the PRC2 core subunit SUZ12. PALI1 further interacts with the H3K9 methyltransferase G9A, bridging the formation of a unique G9A-PALI1-PRC2 super-complex that occupies a subset of G9A-target genes to mediate dual H3K9/K27 methylation and gene repression. Many of these genes are developmental regulators required for cell differentiation, and their loss in PCa predicts poor prognosis. Accordingly, PALI1 and G9A drive PCa cell proliferation and invasion in vitro and xenograft tumor growth in vivo. Collectively, our study shows that PALI1 harnesses two central epigenetic mechanisms to suppress cellular differentiation and promote tumorigenesis, which can be targeted by dual EZH2 and G9A inhibition.