We report a novel computational method, RegNetDriver, to identify tumorigenic drivers using the combined effects of coding and non-coding single nucleotide variants, structural variants, and DNA methylation changes in the DNase I hypersensitivity based regulatory network. Integration of multi-omics data from 521 prostate tumor samples indicated a stronger regulatory impact of structural variants, as they affect more transcription factor hubs in the tissue-specific network. Moreover, crosstalk between transcription factor hub expression modulated by structural variants and methylation levels likely leads to the differential expression of target genes. We report known prostate tumor regulatory drivers and nominate novel transcription factors (ERF, CREB3L1, and POU2F2), which are supported by functional validation.

Mutations in the transcription factor FOXA1 define a unique subset of prostate cancers but the functional consequences of these mutations and whether they confer gain or loss of function is unknown^1-9. Here, by annotating the landscape of FOXA1 mutations from 3,086 human prostate cancers, we define two hotspots in the forkhead domain: Wing2 (around 50% of all mutations) and the highly conserved DNA-contact residue R219 (around 5% of all mutations). Wing2 mutations are detected in adenocarcinomas at all stages, whereas R219 mutations are enriched in metastatic tumours with neuroendocrine histology. Interrogation of the biological properties of wild-type FOXA1 and fourteen FOXA1 mutants reveals gain of function in mouse prostate organoid proliferation assays. Twelve of these mutants, as well as wild-type FOXA1, promoted an exaggerated pro-luminal differentiation program, whereas two different R219 mutants blocked luminal differentiation and activated a mesenchymal and neuroendocrine transcriptional program. Assay for transposase-accessible chromatin using sequencing (ATAC-seq) of wild-type FOXA1 and representative Wing2 and R219 mutants revealed marked, mutant-specific changes in open chromatin at thousands of genomic loci and exposed sites of FOXA1 binding and associated increases in gene expression. Of note, ATAC-seq peaks in cells expressing R219 mutants lacked the canonical core FOXA1-binding motifs (GTAAAC/T) but were enriched for a related, non-canonical motif (GTAAAG/A), which was preferentially activated by R219-mutant FOXA1 in reporter assays. Thus, FOXA1 mutations alter its pioneering function and perturb normal luminal epithelial differentiation programs, providing further support for the role of lineage plasticity in cancer progression.

Background

Cancer cell metabolism requires sustained pools of intracellular nicotinamide adenine dinucleotide (NAD⁺) which is maintained by a balance of NAD⁺ hydrolase activity and NAD⁺ salvage activity. We recently reported that human prostate cancer can be initiated following oncogene expression in progenitor-like luminal cells marked by low expression of the NAD⁺-consuming enzyme CD38. CD38 expression is reduced in prostate cancer compared to benign prostate, suggesting that tumor cells may reduce CD38 expression in order to enhance pools of NAD⁺. However, little is known about how CD38 expression is repressed in advanced prostate cancer and whether CD38 plays a role in regulating NAD⁺ levels in prostate epithelial cells.

Methods

CD38 expression, its association with recurrence after prostatectomy for clinically localized prostate cancer, and DNA methylation of the CD38 promoter were evaluated in human prostate tissues representing various stages of disease progression. CD38 was inducibly over-expressed in benign and malignant human prostate cell lines in order to determine the effects on cell proliferation and levels of NAD⁺ and NADH. NAD⁺ and NADH were also measured in urogenital tissues from wild-type and CD38 knockout mice.

Results

CD38 mRNA expression was reduced in metastatic castration-resistant prostate cancer compared to localized prostate cancer. In a large cohort of men undergoing radical prostatectomy, CD38 protein expression was inversely correlated with recurrence. We identified methylation of the CD38 promoter in primary and metastatic prostate cancer. Over-expression of wild-type CD38, but not an NAD⁺ hydrolase-deficient mutant, depleted extracellular NAD⁺ levels in benign and malignant prostate cell lines. However, expression of CD38 did not significantly alter intracellular NAD⁺ levels in human prostate cell lines grown in vitro and in urogenital tissues isolated from wild-type and CD38 knockout mice.

Conclusions

CD38 protein expression in prostate cancer is associated with risk of recurrence. Methylation results suggest that CD38 is epigenetically regulated in localized and metastatic prostate cancer tissues. Our study provides support for CD38 as a regulator of extracellular, but not intracellular, NAD⁺ in epithelial cells. These findings suggest that repression of CD38 by methylation may serve to increase the availability of extracellular NAD⁺ in prostate cancer tissues.

To broaden our understanding of human neurodevelopment, we profiled transcriptomic and epigenomic landscapes across brain regions and/or cell types for the entire span of prenatal and postnatal development. Integrative analysis revealed temporal, regional, sex, and cell type-specific dynamics. We observed a global transcriptomic cup-shaped pattern, characterized by a late fetal transition associated with sharply decreased regional differences and changes in cellular composition and maturation, followed by a reversal in childhood-adolescence, and accompanied by epigenomic reorganizations. Analysis of gene coexpression modules revealed relationships with epigenomic regulation and neurodevelopmental processes. Genes with genetic associations to brain-based traits and neuropsychiatric disorders (including MEF2C, SATB2, SOX5, TCF4, and TSHZ3) converged in a small number of modules and distinct cell types, revealing insights into neurodevelopment and the genomic basis of neuropsychiatric risks.

There is substantial heterogeneity among primary prostate cancers, evident in the spectrum of molecular abnormalities and its variable clinical course. As part of The Cancer Genome Atlas (TCGA), we present a comprehensive molecular analysis of 333 primary prostate carcinomas. Our results revealed a molecular taxonomy in which 74% of these tumors fell into one of seven subtypes defined by specific gene fusions (ERG, ETV1/4, and FLI1) or mutations (SPOP, FOXA1, and IDH1). Epigenetic profiles showed substantial heterogeneity, including an IDH1 mutant subset with a methylator phenotype. Androgen receptor (AR) activity varied widely and in a subtype-specific manner, with SPOP and FOXA1 mutant tumors having the highest levels of AR-induced transcripts. 25% of the prostate cancers had a presumed actionable lesion in the PI3K or MAPK signaling pathways, and DNA repair genes were inactivated in 19%. Our analysis reveals molecular heterogeneity among primary prostate cancers, as well as potentially actionable molecular defects.