An Atlas of Accessible Chromatin in Advanced Prostate Cancer Reveals the Epigenetic Evolution during Tumor Progression.

Metastatic castration-resistant prostate cancer (mCRPC) is a lethal disease that resists therapy targeting androgen signaling, the primary driver of prostate cancer. mCRPC resists androgen receptor (AR) inhibitors by amplifying AR signaling or by evolving into therapy-resistant subtypes that do not depend on AR. Elucidation of the epigenetic underpinnings of these subtypes could provide important insights into the drivers of therapy resistance. In this study, we produced chromatin accessibility maps linked to the binding of lineage-specific transcription factors (TF) by performing ATAC sequencing on 70 mCRPC tissue biopsies integrated with transcriptome and whole genome sequencing. mCRPC had a distinct global chromatin accessibility profile linked to AR function. Analysis of TF occupancy across accessible chromatin revealed 203 TFs associated with mCRPC subtypes. Notably, ZNF263 was identified as a putative prostate cancer TF with a significant impact on gene activity in the double-negative (AR- neuroendocrine-) subtype, potentially activating MYC targets. Overall, this analysis of chromatin accessibility in mCRPC provides valuable insights into epigenetic changes that occur during progression to mCRPC.

HER2 overexpression in urothelial carcinoma with GATA3 and PPARG copy number gains.

HER2, encoded by the ERBB2 gene, is an important druggable driver of human cancer gaining increasing importance as a therapeutic target in urothelial carcinoma (UC). The genomic underpinnings of HER2 overexpression in ERBB2 nonamplified UC are poorly defined. To address this knowledge gap, we investigated 172 UC tumors from patients treated at the University of California San Francisco, using immunohistochemistry and next-generation sequencing. We found that GATA3 and PPARG copy number gains individually predicted HER2 protein expression independently of ERBB2 amplification. To validate these findings, we interrogated the Memorial Sloan Kettering/The Cancer Genome Atlas (MSK/TCGA) dataset and found that GATA3 and PPARG copy number gains individually predicted ERBB2 mRNA expression independently of ERBB2 amplification. Our findings reveal a potential link between the luminal marker HER2 and the key transcription factors GATA3 and PPARG in UC and highlight the utility of examining GATA3 and PPARG copy number states to identify UC tumors that overexpress HER2 in the absence of ERBB2 amplification. In summary, we found that an increase in copy number of GATA3 and PPARG was independently associated with higher ERBB2 expression in patient samples of UC. This finding provides a potential explanation for HER2 overexpression in UC tumors without ERBB2 amplification and a way to identify these tumors for HER2-targeted therapies.

Integrative analysis of ultra-deep RNA-seq reveals alternative promoter usage as a mechanism of activating oncogenic programmes during prostate cancer progression.

Transcription factor (TF) proteins regulate gene activity by binding to regulatory regions, most importantly at gene promoters. Many genes have alternative promoters (APs) bound by distinct TFs. The role of differential TF activity at APs during tumour development is poorly understood. Here we show, using deep RNA sequencing in 274 biopsies of benign prostate tissue, localized prostate tumours and metastatic castration-resistant prostate cancer, that AP usage increases as tumours progress and APs are responsible for a disproportionate amount of tumour transcriptional activity. Expression of the androgen receptor (AR), the key driver of prostate tumour activity, is correlated with elevated AP usage. We identified AR, FOXA1 and MYC as potential drivers of AP activation. DNA methylation is a likely mechanism for AP activation during tumour progression and lineage plasticity. Our data suggest that prostate tumours activate APs to magnify the transcriptional impact of tumour drivers, including AR and MYC.

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.

Single cell-transcriptomic analysis informs the lncRNA landscape in metastatic castration resistant prostate cancer.

Metastatic castration-resistant prostate cancer (mCRPC) is a lethal form of prostate cancer. Although long-noncoding RNAs (lncRNAs) have been implicated in mCRPC, past studies have relied on bulk sequencing methods with low depth and lack of single-cell resolution. Hence, we performed a lncRNA-focused analysis of single-cell RNA-sequencing data (n = 14) from mCRPC biopsies followed by integration with bulk multi-omic datasets. This yielded 389 cell-enriched lncRNAs in prostate cancer cells and the tumor microenvironment (TME). These lncRNAs demonstrated enrichment with regulatory elements and exhibited alterations during prostate cancer progression. Prostate-lncRNAs were correlated with AR mutational status and response to treatment with enzalutamide, while TME-lncRNAs were associated with RB1 deletions and poor prognosis. Finally, lncRNAs identified between prostate adenocarcinomas and neuroendocrine tumors exhibited distinct expression and methylation profiles. Our findings demonstrate the ability of single-cell analysis to refine our understanding of lncRNAs in mCRPC and serve as a resource for future mechanistic studies.

Transcriptomic Heterogeneity of Expansile Cribriform and Other Gleason Pattern 4 Prostate Cancer Subtypes.

BACKGROUND: Prostate cancers featuring an expansile cribriform (EC) pattern are associated with worse clinical outcomes following radical prostatectomy (RP). However, studies of the genomic characteristics of Gleason pattern 4 subtypes are limited. OBJECTIVE: To explore transcriptomic characteristics and heterogeneity within Gleason pattern 4 subtypes (fused/poorly formed, glomeruloid, small cribriform, EC/intraductal carcinoma [IDC]) and the association with biochemical recurrence (BCR)-free survival. DESIGN, SETTING, AND PARTICIPANTS: This was a retrospective cohort study including 165 men with grade group 2-4 prostate cancer who underwent RP at a single academic institution (2016-2020) and Decipher testing of the RP specimen. Patients with Gleason pattern 5 were excluded. IDC and EC patterns were grouped. Median follow-up was 2.5 yr after RP for patients without BCR. OUTCOMES MEASUREMENTS AND STATISTICAL ANALYSIS: Prompted by heterogeneity within pattern 4 subtypes identified via exploratory analyses, we investigated transcriptomic consensus clusters using partitioning around medoids and hallmark gene set scores. The primary clinical outcome was BCR, defined as two consecutive prostate-specific antigen measurements >0.2 ng/ml at least 8 wk after RP, or any additional treatment. Multivariable Cox proportional-hazards models were used to determine factors associated with BCR-free survival. RESULTS AND LIMITATIONS: In this cohort, 99/165 patients (60%) had EC and 67 experienced BCR. Exploratory analyses and clustering demonstrated transcriptomic heterogeneity within each Gleason pattern 4 subtype. In the multivariable model controlled for pattern 4 subtype, margin status, Cancer of the Prostate Risk Assessment Post-Surgical score, and Decipher score, a newly identified steroid hormone-driven cluster (hazard ratio 2.35 95% confidence interval 1.01-5.47) was associated with worse BCR-free survival. The study is limited by intermediate follow-up, no validation cohort, and lack of accounting for intratumoral and intraprostatic heterogeneity. CONCLUSIONS: Transcriptomic heterogeneity was present within and across each Gleason pattern 4 subtype, demonstrating there is additional biologic diversity not captured by histologic subtypes. This heterogeneity can be used to develop novel signatures and to classify transcriptomic subtypes, which may help in refining risk stratification following RP to further guide decision-making on adjuvant and salvage treatments. PATIENT SUMMARY: We studied prostatectomy specimens and found that tumors with similar microscopic appearance can have genetic differences that may help to predict outcomes after prostatectomy for prostate cancer. Our results demonstrate that further gene expression analysis of prostate cancer subtypes may improve risk stratification after prostatectomy. Future studies are needed to develop novel gene expression signatures and validate these findings in independent sets of patients.

AARDVARK: an automated reversion detector for variants affecting resistance kinetics.

SUMMARY: Resistance to two classes of FDA-approved therapies that target DNA repair-deficient tumors is caused by mutations that restore the tumor cell's DNA repair function. Identifying these "reversion" mutations currently requires manual annotation of patient tumor sequence data. Here we present AARDVARK, an R package that automatically identifies reversion mutations from DNA sequence data. AVAILABILITY AND IMPLEMENTATION: AARDVARK is implemented in R (≥3.5). It is available on GitHub at https://github.com/davidquigley/aardvark. It is licensed under the MIT license.

The Genomic and Epigenomic Landscape of Double-Negative Metastatic Prostate Cancer.

Systemic targeted therapy in prostate cancer is primarily focused on ablating androgen signaling. Androgen deprivation therapy and second-generation androgen receptor (AR)-targeted therapy selectively favor the development of treatment-resistant subtypes of metastatic castration-resistant prostate cancer (mCRPC), defined by AR and neuroendocrine (NE) markers. Molecular drivers of double-negative (AR-/NE-) mCRPC are poorly defined. In this study, we comprehensively characterized treatment-emergent mCRPC by integrating matched RNA sequencing, whole-genome sequencing, and whole-genome bisulfite sequencing from 210 tumors. AR-/NE- tumors were clinically and molecularly distinct from other mCRPC subtypes, with the shortest survival, amplification of the chromatin remodeler CHD7, and PTEN loss. Methylation changes in CHD7 candidate enhancers were linked to elevated CHD7 expression in AR-/NE+ tumors. Genome-wide methylation analysis nominated Krüppel-like factor 5 (KLF5) as a driver of the AR-/NE- phenotype, and KLF5 activity was linked to RB1 loss. These observations reveal the aggressiveness of AR-/NE- mCRPC and could facilitate the identification of therapeutic targets in this highly aggressive disease. SIGNIFICANCE: Comprehensive characterization of the five subtypes of metastatic castration-resistant prostate cancer identified transcription factors that drive each subtype and showed that the double-negative subtype has the worst prognosis.

Reversible epigenetic alterations mediate PSMA expression heterogeneity in advanced metastatic prostate cancer.

Prostate-specific membrane antigen (PSMA) is an important cell surface target in prostate cancer. There are limited data on the heterogeneity of PSMA tissue expression in metastatic castration-resistant prostate cancer (mCRPC). Furthermore, the mechanisms regulating PSMA expression (encoded by the FOLH1 gene) are not well understood. Here, we demonstrate that PSMA expression is heterogeneous across different metastatic sites and molecular subtypes of mCRPC. In a rapid autopsy cohort in which multiple metastatic sites per patient were sampled, we found that 13 of 52 (25%) cases had no detectable PSMA and 23 of 52 (44%) cases showed heterogeneous PSMA expression across individual metastases, with 33 (63%) cases harboring at least 1 PSMA-negative site. PSMA-negative tumors displayed distinct transcriptional profiles with expression of druggable targets such as MUC1. Loss of PSMA was associated with epigenetic changes of the FOLH1 locus, including gain of CpG methylation and loss of histone 3 lysine 27 (H3K27) acetylation. Treatment with histone deacetylase (HDAC) inhibitors reversed this epigenetic repression and restored PSMA expression in vitro and in vivo. Collectively, these data provide insights into the expression patterns and regulation of PSMA in mCRPC and suggest that epigenetic therapies - in particular, HDAC inhibitors - can be used to augment PSMA levels.

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.

Resistance to ATR Inhibitors Is Mediated by Loss of the Nonsense-Mediated Decay Factor UPF2.

Over one million cases of gastric cancer are diagnosed each year globally, and the metastatic disease continues to have a poor prognosis. A significant proportion of gastric tumors have defects in the DNA damage response pathway, creating therapeutic opportunities through synthetic lethal approaches. Several small-molecule inhibitors of ATR, a key regulator of the DNA damage response, are now in clinical development as targeted agents for gastric cancer. Here, we performed a large-scale CRISPR interference screen to discover genetic determinants of response and resistance to ATR inhibitors (ATRi) in gastric cancer cells. Among the top hits identified as mediators of ATRi response were UPF2 and other components of the nonsense-mediated decay (NMD) pathway. Loss of UPF2 caused ATRi resistance across multiple gastric cancer cell lines. Global proteomic, phosphoproteomic, and transcriptional profiling experiments revealed that cell-cycle progression and DNA damage responses were altered in UPF2-mutant cells. Further studies demonstrated that UPF2-depleted cells failed to accumulate in G1 following treatment with ATRi. UPF2 loss also reduced transcription-replication collisions, which has previously been associated with ATRi response, thereby suggesting a possible mechanism of resistance. Our results uncover a novel role for NMD factors in modulating response to ATRi in gastric cancer, highlighting a previously unknown mechanism of resistance that may inform the clinical use of these drugs. SIGNIFICANCE: Loss of NMD proteins promotes resistance to ATR inhibitors in gastric cancer cells, which may provide a combination of therapeutic targets and biomarkers to improve the clinical utility of these drugs.

Transcriptional profiling of matched patient biopsies clarifies molecular determinants of enzalutamide-induced lineage plasticity.

The androgen receptor (AR) signaling inhibitor enzalutamide (enza) is one of the principal treatments for metastatic castration-resistant prostate cancer (CRPC). Several emergent enza clinical resistance mechanisms have been described, including lineage plasticity in which the tumors manifest reduced dependency on the AR. To improve our understanding of enza resistance, herein we analyze the transcriptomes of matched biopsies from men with metastatic CRPC obtained prior to treatment and at progression (n = 21). RNA-sequencing analysis demonstrates that enza does not induce marked, sustained changes in the tumor transcriptome in most patients. However, three patients' progression biopsies show evidence of lineage plasticity. The transcription factor E2F1 and pathways linked to tumor stemness are highly activated in baseline biopsies from patients whose tumors undergo lineage plasticity. We find a gene signature enriched in these baseline biopsies that is strongly associated with poor survival in independent patient cohorts and with risk of castration-induced lineage plasticity in patient-derived xenograft models, suggesting that tumors harboring this gene expression program may be at particular risk for resistance mediated by lineage plasticity and poor outcomes.

Why ARNT Prostate Tumors Responding to Enzalutamide?

SUMMARY: Prostate tumors can develop resistance to androgen receptor (AR)-targeted therapies through treatment-induced changes in transcription factor activity that promote transcriptional and morphologic features of a neuroendocrine lineage. This study identifies an unexpected role for the circadian protein ARNTL in resistance to enzalutamide, a second-generation AR-targeted therapy. See related article by Linder et al., p. 2074 (4).

Deep whole-genome ctDNA chronology of treatment-resistant prostate cancer.

Circulating tumour DNA (ctDNA) in blood plasma is an emerging tool for clinical cancer genotyping and longitudinal disease monitoring1. However, owing to past emphasis on targeted and low-resolution profiling approaches, our understanding of the distinct populations that comprise bulk ctDNA is incomplete2-12. Here we perform deep whole-genome sequencing of serial plasma and synchronous metastases in patients with aggressive prostate cancer. We comprehensively assess all classes of genomic alterations and show that ctDNA contains multiple dominant populations, the evolutionary histories of which frequently indicate whole-genome doubling and shifts in mutational processes. Although tissue and ctDNA showed concordant clonally expanded cancer driver alterations, most individual metastases contributed only a minor share of total ctDNA. By comparing serial ctDNA before and after clinical progression on potent inhibitors of the androgen receptor (AR) pathway, we reveal population restructuring converging solely on AR augmentation as the dominant genomic driver of acquired treatment resistance. Finally, we leverage nucleosome footprints in ctDNA to infer mRNA expression in synchronously biopsied metastases, including treatment-induced changes in AR transcription factor signalling activity. Our results provide insights into cancer biology and show that liquid biopsy can be used as a tool for comprehensive multi-omic discovery.

Predicting cancer drug TARGETS - TreAtment Response Generalized Elastic-neT Signatures.

We are now in an era of molecular medicine, where specific DNA alterations can be used to identify patients who will respond to specific drugs. However, there are only a handful of clinically used predictive biomarkers in oncology. Herein, we describe an approach utilizing in vitro DNA and RNA sequencing and drug response data to create TreAtment Response Generalized Elastic-neT Signatures (TARGETS). We trained TARGETS drug response models using Elastic-Net regression in the publicly available Genomics of Drug Sensitivity in Cancer (GDSC) database. Models were then validated on additional in-vitro data from the Cancer Cell Line Encyclopedia (CCLE), and on clinical samples from The Cancer Genome Atlas (TCGA) and Stand Up to Cancer/Prostate Cancer Foundation West Coast Prostate Cancer Dream Team (WCDT). First, we demonstrated that all TARGETS models successfully predicted treatment response in the separate in-vitro CCLE treatment response dataset. Next, we evaluated all FDA-approved biomarker-based cancer drug indications in TCGA and demonstrated that TARGETS predictions were concordant with established clinical indications. Finally, we performed independent clinical validation in the WCDT and found that the TARGETS AR signaling inhibitors (ARSI) signature successfully predicted clinical treatment response in metastatic castration-resistant prostate cancer with a statistically significant interaction between the TARGETS score and PSA response (p = 0.0252). TARGETS represents a pan-cancer, platform-independent approach to predict response to oncologic therapies and could be used as a tool to better select patients for existing therapies as well as identify new indications for testing in prospective clinical trials.

Prognosis Associated With Luminal and Basal Subtypes of Metastatic Prostate Cancer.

IMPORTANCE: Luminal and basal subtypes of primary prostate cancer have been shown to be molecularly distinct and clinically important in predicting response to therapy. These subtypes have not been described in metastatic prostate cancer. OBJECTIVES: To identify clinical and molecular correlates of luminal and basal subtypes in metastatic castration-resistant prostate cancer (mCRPC) and investigate differences in survival, particularly after treatment with androgen-signaling inhibitors (ASIs). DESIGN, SETTING, AND PARTICIPANTS: In this cohort study, a retrospective analysis was conducted of 4 cohorts with mCRPC (N = 634) across multiple academic centers. Treatment was at the physicians' discretion. Details of the study cohorts have been published elsewhere between 2016 and 2019. Data were analyzed from March 2018 to February 2021. MAIN OUTCOMES AND MEASURES: The primary clinical end point was overall survival from the date of tissue biopsy/molecular profiling. Luminal and basal subtypes were also stratified by postbiopsy ASI treatment. The primary molecular analyses included associations with small cell/neuroendocrine prostate cancer (SCNC), molecular pathways, and DNA alterations. RESULTS: In the 634 patients, 288 (45%) had tumors classified as luminal, and 346 (55%) had tumors classified as basal. However, 53 of 59 (90%) SCNC tumors were basal (P < .001). Similar to primary prostate cancer, luminal tumors exhibited overexpression of AR pathway genes. In basal tumors, a significantly higher rate of RB1 loss (23% basal vs 4% luminal; P < .001), FOXA1 alterations (36% basal vs 27% luminal; P = .03) and MYC alterations (73% basal vs 56% luminal; P < .001) were identified. Patients with basal tumors had worse overall survival compared with those with luminal tumors only in patients treated with an ASI postbiopsy (East Coast Dream Team: hazard ratio [HR], 0.39; 95% CI, 0.20-0.74; P = .004; West Coast Dream Team: HR, 0.57; 95% CI, 0.33-0.97; P = .04). Among patients with luminal tumors, those treated with an ASI had significantly better survival (HR, 0.27; 95% CI, 0.14-0.53; P < .001), whereas patients with basal tumors did not (HR, 0.62; 95% CI, 0.36-1.04, P = .07). The interaction term between subtype and ASI treatment was statistically significant (HR, 0.42; 95% CI, 0.20-0.89; P = .02). CONCLUSIONS AND RELEVANCE: These findings represent the largest integrated clinical, transcriptomic, and genomic analysis of mCRPC samples to date, and suggest that mCRPC can be classified as luminal and basal tumors. Analogous to primary prostate cancer, these data suggest that the benefit of ASI treatment is more pronounced in luminal tumors and support the use of ASIs in this population. In the basal tumors, a chemotherapeutic approach could be considered in some patients given the similarity to SCNC and the diminished benefit of ASI therapy. Further validation in prospective clinical trials is warranted.

Circadian PERformance in breast cancer: a germline and somatic genetic study of PER3VNTR polymorphisms and gene co-expression.

Polymorphisms in the PER3 gene have been associated with several human disease phenotypes, including sleep disorders and cancer. In particular, the long allele of a variable number of tandem repeat (VNTR) polymorphism has been previously linked to an increased risk of breast cancer. Here we carried out a combined germline and somatic genetic analysis of the role of the PER3VNRT polymorphism in breast cancer. The combined data from 8284 individuals showed a non-significant trend towards increased breast cancer risk in the 5-repeat allele homozygous carriers (OR = 1.17, 95% CI: 0.97-1.42). We observed allelic imbalance at the PER3 locus in matched blood and tumor DNA samples, showing a significant retention of the long variant (risk) allele in tumor samples, and a preferential loss of the short repetition allele (p = 0.0005). Gene co-expression analysis in healthy and tumoral breast tissue samples uncovered significant associations between PER3 expression levels with those from genes which belong to several cancer-associated pathways. Finally, relapse-free survival (RFS) analysis showed that low expression levels of PER3 were linked to a significant lower RSF in luminal A (p = 3 × 10-12) but not in the rest of breast cancer subtypes.

An integrated functional and clinical genomics approach reveals genes driving aggressive metastatic prostate cancer.

Genomic sequencing of thousands of tumors has revealed many genes associated with specific types of cancer. Similarly, large scale CRISPR functional genomics efforts have mapped genes required for cancer cell proliferation or survival in hundreds of cell lines. Despite this, for specific disease subtypes, such as metastatic prostate cancer, there are likely a number of undiscovered tumor specific driver genes that may represent potential drug targets. To identify such genetic dependencies, we performed genome-scale CRISPRi screens in metastatic prostate cancer models. We then created a pipeline in which we integrated pan-cancer functional genomics data with our metastatic prostate cancer functional and clinical genomics data to identify genes that can drive aggressive prostate cancer phenotypes. Our integrative analysis of these data reveals known prostate cancer specific driver genes, such as AR and HOXB13, as well as a number of top hits that are poorly characterized. In this study we highlight the strength of an integrated clinical and functional genomics pipeline and focus on two top hit genes, KIF4A and WDR62. We demonstrate that both KIF4A and WDR62 drive aggressive prostate cancer phenotypes in vitro and in vivo in multiple models, irrespective of AR-status, and are also associated with poor patient outcome.

CRISPRi screens reveal a DNA methylation-mediated 3D genome dependent causal mechanism in prostate cancer.

Prostate cancer (PCa) risk-associated SNPs are enriched in noncoding cis-regulatory elements (rCREs), yet their modi operandi and clinical impact remain elusive. Here, we perform CRISPRi screens of 260 rCREs in PCa cell lines. We find that rCREs harboring high risk SNPs are more essential for cell proliferation and H3K27ac occupancy is a strong indicator of essentiality. We also show that cell-line-specific essential rCREs are enriched in the 8q24.21 region, with the rs11986220-containing rCRE regulating MYC and PVT1 expression, cell proliferation and tumorigenesis in a cell-line-specific manner, depending on DNA methylation-orchestrated occupancy of a CTCF binding site in between this rCRE and the MYC promoter. We demonstrate that CTCF deposition at this site as measured by DNA methylation level is highly variable in prostate specimens, and observe the MYC eQTL in the 8q24.21 locus in individuals with low CTCF binding. Together our findings highlight a causal mechanism synergistically driven by a risk SNP and DNA methylation-mediated 3D genome architecture, advocating for the integration of genetics and epigenetics in assessing risks conferred by genetic predispositions.