The functional consequences of somatic non-coding mutations in ovarian cancer (OC) are unknown. To identify regulatory elements (RE) and genes perturbed by acquired non-coding variants, here we establish epigenomic and transcriptomic landscapes of primary OCs using H3K27ac ChIP-seq and RNA-seq, and then integrate these with whole genome sequencing data from 232 OCs. We identify 25 frequently mutated regulatory elements, including an enhancer at 6p22.1 which associates with differential expression of ZSCAN16 (P = 6.6 × 10-4) and ZSCAN12 (P = 0.02). CRISPR/Cas9 knockout of this enhancer induces downregulation of both genes. Globally, there is an enrichment of single nucleotide variants in active binding sites for TEAD4 (P = 6 × 10-11) and its binding partner PAX8 (P = 2×10-10), a known lineagespecific transcription factor in OC. In addition, the collection of cis REs associated with PAX8 comprise the most frequently mutated set of enhancers in OC (P = 0.003). These data indicate that non-coding somatic mutations disrupt the PAX8 transcriptional network during OC development.

The CaCTS algorithm nominates cancer cell master transcription factors and guides a model of ovarian cancer regulatory circuitry.

No abstract

Purpose: To identify molecular subclasses of clear cell ovarian carcinoma (CCOC) and assess their impact on clinical presentation and outcomes. Experimental Design: We profiled 421 primary CCOCs that passed quality control using a targeted deep sequencing panel of 163 putative CCOC driver genes and whole transcriptome sequencing of 211 of these tumors. Molecularly-defined subgroups were identified and tested for association with clinical characteristics and overall survival. Results: We detected a putative somatic driver mutation in at least one candidate gene in 95% (401 out of 421) of CCOC tumors including: ARID1A (in 49% of tumors), PIK3CA (49%), TERT (20%) and TP53 (16%). Clustering of cancer driver mutations and RNA expression converged upon two distinct subclasses of CCOC. The first was dominated by ARID1A-mutated tumors with enriched expression of canonical CCOC genes and markers of platinum resistance; the second was largely comprised of tumors with TP53-mutations and enriched for the expression of genes involved in extracellular matrix organization and mesenchymal differentiation. Compared to the ARID1A-mutated group, women with TP53-mutated tumors were more likely to have advanced stage disease, no antecedent history of endometriosis, and poorer survival, driven by their advanced stage at presentation. In women with ARID1A-mutated tumors, there was a trend towards lower response rate to first-line platinum-based therapy. Conclusions: Our study suggests that CCOC consists of two distinct molecular subclasses with distinct clinical presentation and outcomes, with potential relevance to both traditional and experimental therapy responsiveness.

The function of critical developmental regulators can be subverted by cancer cells to control expression of oncogenic transcriptional programs. These "master transcription factors" (MTFs) are often essential for cancer cell survival and represent vulnerabilities that can be exploited therapeutically. The current approaches to identify candidate MTFs examine super-enhancer associated transcription factor-encoding genes with high connectivity in network models. This relies on chromatin immunoprecipitation-sequencing (ChIP-seq) data, which is technically challenging to obtain from primary tumors, and is currently unavailable for many cancer types and clinically relevant subtypes. In contrast, gene expression data are more widely available, especially for rare tumors and subtypes where MTFs have yet to be discovered. We have developed a predictive algorithm called CaCTS (Cancer Core Transcription factor Specificity) to identify candidate MTFs using pancancer RNA-sequencing data from The Cancer Genome Atlas. The algorithm identified 273 candidate MTFs across 34 tumor types and recovered known tumor MTFs. We also made novel predictions, including for cancer types and subtypes for which MTFs have not yet been characterized. Clustering based on MTF predictions reproduced anatomic groupings of tumors that share 1-2 lineage-specific candidates, but also dictated functional groupings, such as a squamous group that comprised five tumor subtypes sharing 3 common MTFs. PAX8, SOX17, and MECOM were candidate factors in high-grade serous ovarian cancer (HGSOC), an aggressive tumor type where the core regulatory circuit is currently uncharacterized. PAX8, SOX17, and MECOM are required for cell viability and lie proximal to super-enhancers in HGSOC cells. ChIPseq revealed that these factors co-occupy HGSOC regulatory elements globally and co-bind at critical gene loci including MUC16 (CA-125). Addiction to these factors was confirmed in studies using THZ1 to inhibit transcription in HGSOC cells, suggesting early down-regulation of these genes may be responsible for cytotoxic effects of THZ1 on HGSOC models. Identification of MTFs across 34 tumor types and 140 subtypes, especially for those with limited understanding of transcriptional drivers paves the way to therapeutic targeting of MTFs in a broad spectrum of cancers.

RNA molecules function as messenger RNAs (mRNAs) that encode proteins and noncoding transcripts that serve as adaptor molecules, structural components, and regulators of genome organization and gene expression. Their function and regulation are largely mediated by RNA binding proteins (RBPs). Here we present RNA proximity labelling (RPL), an RNA-centric method comprising the endonuclease-deficient Type VI CRISPR-Cas protein dCas13b fused to engineered ascorbate peroxidase APEX2. RPL discovers target RNA proximal proteins in vivo via proximity-based biotinylation. RPL applied to U1 identified proteins involved in both U1 canonical and noncanonical functions. Profiling of poly(A) tail proximal proteins uncovered expected categories of RBPs and provided additional evidence for 5ʹ-3ʹ proximity and unexplored subcellular localizations of poly(A) + RNA. Our results suggest that RPL allows rapid identification of target RNA binding proteins in native cellular contexts, and is expected to pave the way for discovery of novel RNA-protein interactions important for health and disease.

PAX8 is a master transcription factor that is essential during embryogenesis and promotes neoplastic growth. It is expressed by the secretory cells lining the female reproductive tract, and its deletion during development results in atresia of reproductive tract organs. Nearly all ovarian carcinomas express PAX8, and its knockdown results in apoptosis of ovarian cancer cells. To explore the role of PAX8 in these tissues, we purified the PAX8 protein complex from nonmalignant fallopian tube cells and high-grade serous ovarian carcinoma cell lines. We found that PAX8 was a member of a large chromatin remodeling complex and preferentially interacted with SOX17, another developmental transcription factor. Depleting either PAX8 or SOX17 from cancer cells altered the expression of factors involved in angiogenesis and functionally disrupted tubule and capillary formation in cell culture and mouse models. PAX8 and SOX17 in ovarian cancer cells promoted the secretion of angiogenic factors by suppressing the expression of SERPINE1 , which encodes a proteinase inhibitor with anti a ngiogenic effects. The findings reveal a non–cell-autonomous function of these transcription factors in regulating angiogenesis in ovarian cancer.

Transcriptional regulation is highly disease and cell-type specific. We performed H3K27ac chromatin immunoprecipitation and transcriptomic sequencing in primary tumors for the four different subtypes of invasive epithelial ovarian cancer (OC). Histotype-specific regulatory elements (REs) were enriched in enhancers (P<0.001). In silico prediction of putative target genes for histotype-specific REs identified genes ( WFDC2 , P=5.5x10 -5 ) and pathways (PI3K-Akt signaling, P<0.002) known to be involved in OC development. Some genes (e.g. PAX8 and CA125 ) are associated with super-enhancers (SEs) in all OCs, while others are histotype-specific, including PPP1R3B which is associated with SEs specific to clear cell OC. Integrated analysis of active chromatin landscapes with somatic single nucleotide variants (SNVs) from whole genome sequencing (WGS) of 232 primary OCs identified frequently mutated REs, including the KLF6 promoter (P=8.2x10 -8 ) and a putative enhancer at chromosome 6p22.1 (P<0.05). In high-grade serous OCs, somatic SNVs clustered in binding sites for the PAX8 binding partner TEAD4 (P=6x10 -11 ), while the collection of cis regulatory elements associated with PAX8 was the most frequently mutated set of enhancers in OC (P=0.003).Functional analyses supported our findings: Knockdown of PPP1R3B in clear cell OC cells significantly reduced intracellular glycogen content, a signature feature of this histotype; and stable knockout of a 635 bp region in the 6p22.1 enhancer induced downregulation of two predicted target genes, ZSCAN16 and ZSCAN12 (P=6.6 x 10 -4 and P=0.02). In summary, we have characterized histotype-specific epigenomic and transcriptomic landscapes in OC and defined likely functional REs based on somatic mutation analysis of ovarian tumors.