Non-coding Somatic Mutations Converge on the PAX8 Pathway in Epithelial Ovarian Cancer
ABSTRACTTranscriptional 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 silicoprediction of putative target genes for histotype-specific REs identified genes (WFDC2, P=5.5×10-5) and pathways (PI3K-Akt signaling, P<0.002) known to be involved in OC development. Some genes (e.g.PAX8andCA125) are associated with super-enhancers (SEs) in all OCs, while others are histotype-specific, includingPPP1R3Bwhich 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 theKLF6promoter (P=8.2×10-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=6×10-11), while the collection ofcisregulatory elements associated withPAX8was the most frequently mutated set of enhancers in OC (P=0.003). Functional analyses supported our findings: Knockdown ofPPP1R3Bin 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,ZSCAN16andZSCAN12(P=6.6 x 10-4and 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.