MYBL2 Drives Prostate Cancer Plasticity: Inhibiting Its Transcriptional Target CDK2 for RB1-Deficient Neuroendocrine Prostate Cancer.

Phenotypic plasticity is a recognized mechanism driving therapeutic resistance in patients with prostate cancer. Although underlying molecular causations driving phenotypic plasticity have been identified, therapeutic success is yet to be achieved. To identify putative master regulator transcription factors (MR-TF) driving phenotypic plasticity in prostate cancer, this work utilized a multiomic approach using genetically engineered mouse models of prostate cancer combined with patient data to identify MYB proto-oncogene like 2 (MYBL2) as a significantly enriched transcription factor in prostate cancer exhibiting phenotypic plasticity. Genetic inhibition of Mybl2 using independent murine prostate cancer cell lines representing phenotypic plasticity demonstrated Mybl2 loss significantly decreased in vivo growth as well as cell fitness and repressed gene expression signatures involved in pluripotency and stemness. Because MYBL2 is currently not druggable, a MYBL2 gene signature was employed to identify cyclin-dependent kinase-2 (CDK2) as a potential therapeutic target. CDK2 inhibition phenocopied genetic loss of Mybl2 and significantly decreased in vivo tumor growth associated with enrichment of DNA damage. Together, this work demonstrates MYBL2 as an important MR-TF driving phenotypic plasticity in prostate cancer. Furthermore, high MYBL2 activity identifies prostate cancer that would be responsive to CDK2 inhibition. SIGNIFICANCE: Prostate cancers that escape therapy targeting the androgen receptor signaling pathways via phenotypic plasticity are currently untreatable. Our study identifies MYBL2 as a MR-TF in phenotypic plastic prostate cancer and implicates CDK2 inhibition as a novel therapeutic target for this most lethal subtype of prostate cancer.

Sexually dimorphic DNA methylation and gene expression patterns in human first trimester placenta.

BACKGROUND: Fetal sex and placental development impact pregnancy outcomes and fetal-maternal health, but the critical timepoint of placenta establishment in first trimester is understudied in human pregnancies. METHODS: Pregnant subjects were recruited in late first trimester (weeks 10-14) at time of chorionic villus sampling, a prenatal diagnostic test. Leftover placenta tissue was collected and stored until birth outcomes were known, then DNA and RNA were isolated from singleton, normal karyotype pregnancies resulting in live births. DNA methylation was measured with the Illumina Infinium MethylationEPIC BeadChip array (n = 56). Differential methylation analysis compared 25 females versus 31 males using a generalized linear model on 743,461 autosomal probes. Gene expression sex differences were analyzed with RNA-sequencing (n = 74). An integrated analysis was performed using linear regression to correlate gene expression and DNA methylation in 51 overlapping placentas. RESULTS: Methylation analysis identified 151 differentially methylated probes (DMPs) significant at false discovery rate < 0.05, including 89 (59%) hypermethylated in females. Probe cg17612569 (GABPA, ATP5J) was the most significant CpG site, hypermethylated in males. There were 11 differentially methylated regions affected by fetal sex, with transcription factors ZNF300 and ZNF311 most significantly hypermethylated in males and females, respectively. RNA-sequencing identified 152 genes significantly sexually dimorphic at false discovery rate < 0.05. The 151 DMPs were associated with 18 genes with gene downregulation (P < 0.05) in the direction of hypermethylation, including 2 genes significant at false discovery rate < 0.05 (ZNF300 and CUB and Sushi multiple domains 1, CSMD1). Both genes, as well as Family With Sequence Similarity 228 Member A (FAM228A), showed significant correlation between DNA methylation and sexually dimorphic gene expression, though FAM228A DNA methylation was less sexually dimorphic. Comparison with other sex differences studies found that cg17612569 is male-hypermethylated across gestation in placenta and in human blood up to adulthood. CONCLUSIONS: Overall, sex dimorphic differential methylation with associated differential gene expression in the first trimester placenta is small, but there remain significant genes that may be regulated through methylation leading to differences in the first trimester placenta.

Fetal sex and placenta development affect pregnancy outcomes for both the fetus and mother throughout pregnancy, including risk of miscarriages, preterm birth, preeclampsia, and other outcomes. Epigenetics, the "overlay" of regulatory signals on DNA which affects how DNA is read, is not well understood in early pregnancy when critical placenta developments are happening that affect the rest of pregnancy. Here, we use leftover placenta biopsy samples (n = 56) donated by Cedars-Sinai patients with informed consent to learn about first trimester human placenta DNA methylation differences due to fetal sex. Out of the total 743,461 sites analyzed, we identified 151 sites significantly affected by fetal sex after correcting p-values to reduce false positives (false discovery rate < 0.05). We also performed an analysis to look at multiple sites and identified 11 regions across the genome with significant DNA methylation changes due to fetal sex. Furthermore, because DNA methylation is a regulatory mark on DNA which typically dampens gene expression, we also compared the DNA methylation sex differences to placental RNA-sequencing gene expression analysis using the same tissue from a mostly overlapping patient group (n = 74 total sequenced, n = 51 overlap). We identify 18 genes which show both significant DNA methylation differences and gene expression changes. The most significant gene was transcription factor ZNF300 with higher DNA methylation in males and reduced gene expression in males (and thus higher gene expression in females). This study identifies some sex differences that continue until later pregnancy and others that are unique to first trimester.

Rare germline genetic variation in PAX8 transcription factor binding sites and susceptibility to epithelial ovarian cancer.

Common genetic variation throughout the genome together with rare coding variants identified to date explain about a half of the inherited genetic component of epithelial ovarian cancer risk. It is likely that rare variation in the non-coding genome will explain some of the unexplained heritability, but identifying such variants is challenging. The primary problem is lack of statistical power to identifying individual risk variants by association as power is a function of sample size, effect size and allele frequency. Power can be increased by using burden tests which test for association of carriers of any variant in a specified genomic region. This has the effect of increasing the putative effect allele frequency. PAX8 is a transcription factor that plays a critical role in tumour progression, migration and invasion. Furthermore, regulatory elements proximal to target genes of PAX8 are enriched for common ovarian cancer risk variants. We hypothesised that rare variation in PAX8 binding sites are also associated with ovarian cancer risk, but unlikely to be associated with risk of breast, colorectal or endometrial cancer. We have used publicly available, whole-genome sequencing data from the UK 100,000 Genomes Project to evaluate the burden of rare variation in PAX8 binding sites across the genome. Data were available for 522 ovarian cancers, 2,984 breast cancers, 2,696 colorectal cancers, 836 endometrial cancers and 2253 non-cancer controls. Active binding sites were defined using data from multiple PAX8 and H3K27 ChIPseq experiments. We found no association between the burden of rare variation in PAX8 binding sites (defined in several ways) and risk of ovarian, breast or endometrial cancer. An apparent association with colorectal cancer was likely to be a technical artefact as a similar association was also detected for rare variation in random regions of the genome. Despite the null result this study provides a proof-of -principle for using burden testing to identify rare, non-coding germline genetic variation associated with disease. Larger sample sizes available from large-scale sequencing projects together with improved understanding of the function of the non-coding genome will increase the potential of similar studies in the future.

High-throughput mRNA sequencing of human placenta shows sex differences across gestation.

INTRODUCTION: Fetal sex affects fetal and maternal health outcomes in pregnancy, but this connection remains poorly understood. As the placenta is the route of fetomaternal communication and derives from the fetal genome, placental gene expression sex differences may explain these outcomes. OBJECTIVES: We utilized next generation sequencing to study the normal human placenta in both sexes in first and third trimester to generate a normative transcriptome based on sex and gestation. STUDY DESIGN: We analyzed 124 first trimester (T1, 59 female and 65 male) and 43 third trimester (T3, 18 female and 25 male) samples for sex differences within each trimester and sex-specific gestational differences. RESULTS: Placenta shows more significant sexual dimorphism in T1, with 94 T1 and 26 T3 differentially expressed genes (DEGs). The sex chromosomes contributed 60.6% of DEGs in T1 and 80.8% of DEGs in T3, excluding X/Y pseudoautosomal regions. There were 6 DEGs from the pseudoautosomal regions, only significant in T1 and all upregulated in males. The distribution of DEGs on the X chromosome suggests genes on Xp (the short arm) may be particularly important in placental sex differences. Dosage compensation analysis of X/Y homolog genes shows expression is primarily contributed by the X chromosome. In sex-specific analyses of first versus third trimester, there were 2815 DEGs common to both sexes upregulated in T1, and 3263 common DEGs upregulated in T3. There were 7 female-exclusive DEGs upregulated in T1, 15 female-exclusive DEGs upregulated in T3, 10 male-exclusive DEGs upregulated in T1, and 20 male-exclusive DEGs upregulated in T3. DISCUSSION: This is the largest cohort of placentas across gestation from healthy pregnancies defining the normative sex dimorphic gene expression and sex common, sex specific and sex exclusive gene expression across gestation. The first trimester has the most sexually dimorphic transcripts, and the majority were upregulated in females compared to males in both trimesters. The short arm of the X chromosome and the pseudoautosomal region is particularly critical in defining sex differences in the first trimester placenta. As pregnancy is a dynamic state, sex specific DEGs across gestation may contribute to sex dimorphic changes in overall outcomes.

MYBL2 drives prostate cancer plasticity and identifies CDK2 as a therapeutic vulnerability in RB1-loss and neuroendocrine prostate cancer.

Phenotypic plasticity is a recognized mechanism driving therapeutic resistance in prostate cancer (PCa) patients. While underlying molecular causations driving phenotypic plasticity have been identified, therapeutic success is yet to be achieved. To identify putative master regulator transcription factors (MR-TF) driving phenotypic plasticity in PCa, this work utilized a multiomic approach using genetically engineered mouse models of prostate cancer combined with patient data to identify MYBL2 as a significantly enriched transcription factor in PCa exhibiting phenotypic plasticity. Genetic inhibition of Mybl2 using independent murine PCa cell lines representing phenotypic plasticity demonstrated Mybl2 loss significantly decreased in vivo growth as well as cell fitness and repressed gene expression signatures involved in pluripotency and stemness. Because MYBL2 is currently not druggable, a MYBL2 gene signature was employed to identify cyclin-dependent kinase-2 (CDK2) as a potential therapeutic target. CDK2 inhibition phenocopied genetic loss of Mybl2 and significantly decreased in vivo tumor growth associated with enrichment of DNA damage. Together, this work demonstrates MYBL2 as an important MR-TF driving phenotypic plasticity in PCa. Further, high MYBL2 activity identifies PCa that would be responsive to CDK2 inhibition. Significance: PCa that escapes therapy targeting the androgen receptor signaling pathways via phenotypic plasticity are currently untreatable. Our study identifies MYBL2 as a MR-TF in phenotypic plastic PCa and implicates CDK2 inhibition as novel therapeutic target for this most lethal subtype of PCa.

Fallopian tube single cell analysis reveals myeloid cell alterations in high-grade serous ovarian cancer.

Most high-grade serous ovarian cancers (HGSCs) likely initiate from fallopian tube (FT) epithelia. While epithelial subtypes have been characterized using single-cell RNA-sequencing (scRNA-Seq), heterogeneity of other compartments and their involvement in tumor progression are poorly defined. Integrated analysis of human FT scRNA-Seq and HGSC-related tissues, including tumors, revealed greater immune and stromal transcriptional diversity than previously reported. We identified abundant monocytes in FTs across two independent cohorts. The ratio of macrophages to monocytes is similar between benign FTs, ovaries, and adjacent normal tissues but significantly greater in tumors. FT-defined monocyte and macrophage signatures, cell-cell communication, and gene set enrichment analyses identified monocyte- and macrophage-specific interactions and functional pathways in paired tumors and adjacent normal tissues. Further reanalysis of HGSC scRNA-Seq identified monocyte and macrophage subsets associated with neoadjuvant chemotherapy. Taken together, our work provides data that an altered FT myeloid cell composition could inform the discovery of early detection markers for HGSC.

High-throughput mRNA-seq atlas of human placenta shows vast transcriptome remodeling from first to third trimester†.

The placenta, composed of chorionic villi, changes dramatically across gestation. Understanding differences in ongoing pregnancies are essential to identify the role of chorionic villi at specific times in gestation and develop biomarkers and prognostic indicators of maternal-fetal health. The normative mRNA profile is established using next-generation sequencing of 124 first trimester and 43 third trimester human placentas from ongoing healthy pregnancies. Stably expressed genes (SEGs) not different between trimesters and with low variability are identified. Differential expression analysis of first versus third trimester adjusted for fetal sex is performed, followed by a subanalysis with 23 matched pregnancies to control for subject variability using the same genetic and environmental background. Placenta expresses 14,979 polyadenylated genes above sequencing noise (transcripts per million > 0.66), with 10.7% SEGs across gestation. Differentially expressed genes (DEGs) account for 86.7% of genes in the full cohort [false discovery rate (FDR) < 0.05]. Fold changes highly correlate between the full cohort and subanalysis (Pearson = 0.98). At stricter thresholds (FDR < 0.001, fold change > 1.5), there remains 50.1% DEGs (3353 upregulated in first and 4155 upregulated in third trimester). This is the largest mRNA atlas of healthy human placenta across gestation, controlling for genetic and environmental factors, demonstrating substantial changes from first to third trimester in chorionic villi. Specific differences and SEGs may be used to understand the specific role of the chorionic villi throughout gestation and develop first trimester biomarkers of placental health that transpire across gestation, which can be used for future development of biomarkers for maternal-fetal health.

α-Ketoglutarate-Dependent KDM6 Histone Demethylases and Interferon-Stimulated Gene Expression in Lupus.

OBJECTIVE: We aimed to investigate the hypothesis that interferon (IFN)-stimulated gene (ISG) expression in systemic lupus erythematosus (SLE) monocytes is linked to changes in metabolic reprogramming and epigenetic regulation of ISG expression. METHODS: Monocytes from healthy volunteers and patients with SLE at baseline or following IFNα treatment were analyzed by extracellular flux analysis, proteomics, metabolomics, chromatin immunoprecipitation, and gene expression. The histone demethylases KDM6A/B were inhibited using glycogen synthase kinase J4 (GSK-J4). GSK-J4 was tested in pristane and resiquimod (R848) models of IFN-driven SLE. RESULTS: SLE monocytes had enhanced rates of glycolysis and oxidative phosphorylation compared to healthy control monocytes, as well as increased levels of isocitrate dehydrogenase and its product, α-ketoglutarate (α-KG). Because α-KG is a required cofactor for histone demethylases KDM6A and KDM6B, we hypothesized that IFNα may be driving "trained immune" responses through altering histone methylation. IFNα priming (day 1) resulted in a sustained increase in the expression of ISGs in primed cells (day 5) and enhanced expression on restimulation with IFNα. Importantly, decreased H3K27 trimethylation was observed at the promoters of ISGs following IFNα priming. Finally, GSK-J4 (KDM6A/B inhibitor) resulted in decreased ISG expression in SLE patient monocytes, as well as reduced autoantibody production, ISG expression, and kidney pathology in R848-treated BALB/c mice. CONCLUSION: Our study suggests long-term IFNα exposure alters the epigenetic regulation of ISG expression in SLE monocytes via changes in immunometabolism, a mechanism reflecting trained immunity to type I IFN. Importantly, it opens the possibility that targeting histone-modifying enzymes, such as KDM6A/B, may reduce IFN responses in SLE.

The tumor multi-omic landscape of endometrial cancers developed on a germline genetic background of adiposity.

High body mass index (BMI) is a causal risk factor for endometrial cancer but the tumor molecular mechanisms affected by adiposity and their therapeutic relevance remain poorly understood. Here we characterize the tumor multi-omic landscape of endometrial cancers that have developed on a background of lifelong germline genetic exposure to elevated BMI. We built a polygenic score (PGS) for BMI in women using data on independent, genome-wide significant variants associated with adult BMI in 434,794 women. We performed germline (blood) genotype quality control and imputation on data from 354 endometrial cancer cases from The Cancer Genome Atlas (TCGA). We assigned each case in this TCGA cohort their genetically predicted life-course BMI based on the BMI PGS. Multivariable generalized linear models adjusted for age, stage, microsatellite status and genetic principal components were used to test for associations between the BMI germline PGS and endometrial cancer tumor genome-wide genomic, transcriptomic, proteomic, epigenomic and immune traits in TCGA. High BMI germline PGS was associated with (i) upregulated tumor gene expression in the IL6-JAK-STAT3 pathway (FDR=4.2×10-7); (ii) increased estimated intra-tumor activated mast cell infiltration (FDR=0.008); (iii) increased single base substitution (SBS) mutational signatures 1 (FDR=0.03) and 5 (FDR=0.09) and decreased SBS13 (FDR=0.09), implicating age-related and APOBEC mutagenesis, respectively; and (iv) decreased tumor EGFR protein expression (FDR=0.07). Alterations in IL6-JAK-STAT3 signaling gene and EGFR protein expression were, in turn, significantly associated with both overall survival and progression-free interval. Thus, we integrated germline and somatic data using a novel study design to identify associations between genetically predicted lifelong exposure to higher BMI and potentially actionable endometrial cancer tumor molecular features. These associations inform our understanding of how high BMI may influence the development and progression of this cancer, impacting endometrial tumor biology and clinical outcomes.

A biallelic multiple nucleotide length polymorphism explains functional causality at 5p15.33 prostate cancer risk locus.

To date, single-nucleotide polymorphisms (SNPs) have been the most intensively investigated class of polymorphisms in genome wide associations studies (GWAS), however, other classes such as insertion-deletion or multiple nucleotide length polymorphism (MNLPs) may also confer disease risk. Multiple reports have shown that the 5p15.33 prostate cancer risk region is a particularly strong expression quantitative trait locus (eQTL) for Iroquois Homeobox 4 (IRX4) transcripts. Here, we demonstrate using epigenome and genome editing that a biallelic (21 and 47 base pairs (bp)) MNLP is the causal variant regulating IRX4 transcript levels. In LNCaP prostate cancer cells (homozygous for the 21 bp short allele), a single copy knock-in of the 47 bp long allele potently alters the chromatin state, enabling de novo functional binding of the androgen receptor (AR) associated with increased chromatin accessibility, Histone 3 lysine 27 acetylation (H3K27ac), and ~3-fold upregulation of IRX4 expression. We further show that an MNLP is amongst the strongest candidate susceptibility variants at two additional prostate cancer risk loci. We estimated that at least 5% of prostate cancer risk loci could be explained by functional non-SNP causal variants, which may have broader implications for other cancers GWAS. More generally, our results underscore the importance of investigating other classes of inherited variation as causal mediators of human traits.

High-throughput mRNA-seq atlas of human placenta shows vast transcriptome remodeling from first to third trimester.

Background: The placenta, composed of chorionic villi, changes dramatically across gestation. Understanding differences in ongoing pregnancies are essential to identify the role of chorionic villi at specific times in gestation and develop biomarkers and prognostic indicators of maternal- fetal health. Methods: The normative mRNA profile is established using next-generation sequencing of 124 first trimester and 43 third trimester human placentas from ongoing healthy pregnancies. Stably expressed genes not different between trimesters and with low variability are identified. Differential expression analysis of first versus third trimester adjusted for fetal sex is performed, followed by a subanalysis with 23 matched pregnancies to control for subject variability using the same genetic and environmental background. Results: Placenta expresses 14,979 mRNAs above sequencing noise (TPM>0.66), with 1,545 stably expressed genes across gestation. Differentially expressed genes account for 86.7% of genes in the full cohort (FDR<0.05). Fold changes highly correlate between the full cohort and subanalysis (Pearson = 0.98). At stricter thresholds (FDR<0.001, fold change>1.5), there are 6,941 differentially expressed protein coding genes (3,206 upregulated in first and 3,735 upregulated in third trimester). Conclusion: This is the largest mRNA atlas of healthy human placenta across gestation, controlling for genetic and environmental factors, demonstrating substantial changes from first to third trimester in chorionic villi. Specific differences and stably expressed genes may be used to understand the specific role of the chorionic villi throughout gestation and develop first trimester biomarkers of placental health that transpire across gestation, which can be used for future development of biomarkers in maternal-fetal disease.

Rewiring of master transcription factor cistromes during high-grade serous ovarian cancer development.

The transcription factors MECOM, PAX8, SOX17 and WT1 are candidate master regulators of high-grade serous 'ovarian' cancer (HGSC), yet their cooperative role in the hypothesized tissue of origin, the fallopian tube secretory epithelium (FTSEC) is unknown. We generated 26 epigenome (CUT&TAG, CUT&RUN, ATAC-seq and HiC) data sets and 24 profiles of RNA-seq transcription factor knock-down followed by RNA sequencing in FTSEC and HGSC models to define binding sites and gene sets regulated by these factors in cis and trans. This revealed that MECOM, PAX8, SOX17 and WT1 are lineage-enriched, super-enhancer associated master regulators whose cooperative DNA-binding patterns and target genes are re-wired during tumor development. All four TFs were indispensable for HGSC clonogenicity and survival but only depletion of PAX8 and WT1 impaired FTSEC cell survival. These four TFs were pharmacologically inhibited by transcriptional inhibitors only in HGSCs but not in FTSECs. Collectively, our data highlights that tumor-specific epigenetic remodeling is tightly related to MECOM, PAX8, SOX17 and WT1 activity and these transcription factors are targetable in a tumor-specific manner through transcriptional inhibitors.

KRAS mutations and endometriosis burden of disease.

The clinical phenotype of somatic mutations in endometriosis is unknown. The objective was to determine whether somatic KRAS mutations were associated with greater disease burden in endometriosis (i.e. more severe subtypes and higher stage). This prospective longitudinal cohort study included 122 subjects undergoing endometriosis surgery at a tertiary referral center between 2013 and 2017, with 5-9 years of follow-up. Somatic activating KRAS codon 12 mutations were detected in endometriosis lesions using droplet digital PCR. KRAS mutation status for each subject was coded as present (KRAS mutation in at least one endometriosis sample in a subject) or absent. Standardized clinical phenotyping for each subject was carried out via linkage to a prospective registry. Primary outcome was anatomic disease burden, based on distribution of subtypes (deep infiltrating endometriosis, ovarian endometrioma, and superficial peritoneal endometriosis) and surgical staging (Stages I-IV). Secondary outcomes were markers of surgical difficulty, demographics, pain scores, and risk of re-operation. KRAS mutation presence was higher in subjects with deep infiltrating endometriosis or endometrioma lesions only (57.9%; 11/19) and subjects with mixed subtypes (60.6%; 40/66), compared with those with superficial endometriosis only (35.1%; 13/37) (p = 0.04). KRAS mutation was present in 27.6% (8/29) of Stage I cases, in comparison to 65.0% (13/20) of Stage II, 63.0% (17/27) of Stage III, and 58.1% (25/43) of Stage IV cases (p = 0.02). KRAS mutation was also associated with greater surgical difficulty (ureterolysis) (relative risk [RR] = 1.47, 95% CI: 1.02-2.11) and non-Caucasian ethnicity (RR = 0.64, 95% CI: 0.47-0.89). Pain severities did not differ based on KRAS mutation status, at either baseline or follow-up. Re-operation rates were low overall, occurring in 17.2% with KRAS mutation compared with 10.3% without (RR = 1.66, 95% CI: 0.66-4.21). In conclusion, KRAS mutations were associated with greater anatomic severity of endometriosis, resulting in increased surgical difficulty. Somatic cancer-driver mutations may inform a future molecular classification of endometriosis.

Single-cell transcriptomic analysis of endometriosis.

Endometriosis is a common condition in women that causes chronic pain and infertility and is associated with an elevated risk of ovarian cancer. We profiled transcriptomes of >370,000 individual cells from endometriomas (n = 8), endometriosis (n = 28), eutopic endometrium (n = 10), unaffected ovary (n = 4) and endometriosis-free peritoneum (n = 4), generating a cellular atlas of endometrial-type epithelial cells, stromal cells and microenvironmental cell populations across tissue sites. Cellular and molecular signatures of endometrial-type epithelium and stroma differed across tissue types, suggesting a role for cellular restructuring and transcriptional reprogramming in the disease. Epithelium, stroma and proximal mesothelial cells of endometriomas showed dysregulation of pro-inflammatory pathways and upregulation of complement proteins. Somatic ARID1A mutation in epithelial cells was associated with upregulation of pro-angiogenic and pro-lymphangiogenic factors and remodeling of the endothelial cell compartment, with enrichment of lymphatic endothelial cells. Finally, signatures of ciliated epithelial cells were enriched in ovarian cancers, reinforcing epidemiologic associations between these two diseases.

Epigenomic charting and functional annotation of risk loci in renal cell carcinoma.

While the mutational and transcriptional landscapes of renal cell carcinoma (RCC) are well-known, the epigenome is poorly understood. We characterize the epigenome of clear cell (ccRCC), papillary (pRCC), and chromophobe RCC (chRCC) by using ChIP-seq, ATAC-Seq, RNA-seq, and SNP arrays. We integrate 153 individual data sets from 42 patients and nominate 50 histology-specific master transcription factors (MTF) to define RCC histologic subtypes, including EPAS1 and ETS-1 in ccRCC, HNF1B in pRCC, and FOXI1 in chRCC. We confirm histology-specific MTFs via immunohistochemistry including a ccRCC-specific TF, BHLHE41. FOXI1 overexpression with knock-down of EPAS1 in the 786-O ccRCC cell line induces transcriptional upregulation of chRCC-specific genes, TFCP2L1, ATP6V0D2, KIT, and INSRR, implicating FOXI1 as a MTF for chRCC. Integrating RCC GWAS risk SNPs with H3K27ac ChIP-seq and ATAC-seq data reveals that risk-variants are significantly enriched in allelically-imbalanced peaks. This epigenomic atlas in primary human samples provides a resource for future investigation.

circCsnk1g3- and circAnkib1-regulated interferon responses in sarcoma promote tumorigenesis by shaping the immune microenvironment.

Exonic circular RNAs (circRNAs) produce predominantly non-coding RNA species that have been recently profiled in many tumors. However, their functional contribution to cancer progression is still poorly understood. Here, we identify the circRNAs expressed in soft tissue sarcoma cells and explore how the circRNAs regulate sarcoma growth in vivo. We show that circCsnk1g3 and circAnkib1 promote tumor growth by shaping a pro-tumorigenic microenvironment, possibly due to their capabilities to regulate tumor-promoting elements extrinsic to the tumor cells. Accordingly, circCsnk1g3 and circAnkib1 can control the expression of interferon-related genes and pro-inflammatory factors in the sarcoma cells, thus directing immune cell recruitment into the tumor mass, and hence their activation. Mechanistically, circRNAs may repress pro-inflammatory elements by buffering activation of the pathways mediated by RIG-I, the cytosolic viral RNA sensor. The current findings suggest that the targeting of specific circRNAs could augment the efficacy of tumor and immune response to mainstay therapies.

DNA methylation and transcriptomic features are preserved throughout disease recurrence and chemoresistance in high grade serous ovarian cancers.

BACKGROUND: Little is known about the role of global DNA methylation in recurrence and chemoresistance of high grade serous ovarian cancer (HGSOC). METHODS: We performed whole genome bisulfite sequencing and transcriptome sequencing in 62 primary and recurrent tumors from 28 patients with stage III/IV HGSOC, of which 11 patients carried germline, pathogenic BRCA1 and/or BRCA2 mutations. RESULTS: Landscapes of genome-wide methylation (on average 24.2 million CpGs per tumor) and transcriptomes in primary and recurrent tumors showed extensive heterogeneity between patients but were highly preserved in tumors from the same patient. We identified significant differences in the burden of differentially methylated regions (DMRs) in tumors from BRCA1/2 compared to non-BRCA1/2 carriers (mean 659 DMRs and 388 DMRs in paired comparisons respectively). We identified overexpression of immune pathways in BRCA1/2 carriers compared to non-carriers, implicating an increased immune response in improved survival (P = 0.006) in these BRCA1/2 carriers. CONCLUSION: These findings indicate methylome and gene expression programs established in the primary tumor are conserved throughout disease progression, even after extensive chemotherapy treatment, and that changes in methylation and gene expression are unlikely to serve as drivers for chemoresistance in HGSOC.

chromMAGMA: regulatory element-centric interrogation of risk variants.

Candidate causal risk variants from genome-wide association studies reside almost exclusively in noncoding regions of the genome and innovative approaches are necessary to understand their biological function. Multi-marker analysis of genomic annotation (MAGMA) is a widely used program that nominates candidate risk genes by mapping single-nucleotide polymorphism summary statistics from genome-wide association studies to gene bodies. We augmented MAGMA to create chromatin-MAGMA (chromMAGMA), a method to nominate candidate risk genes based on the presence of risk variants within noncoding regulatory elements (REs). We applied chromMAGMA to a genetic susceptibility dataset for epithelial ovarian cancer (EOC), a rare gynecologic malignancy characterized by high mortality. This identified 155 unique candidate EOC risk genes across five EOC histotypes; 83% (105/127) of high-grade serous ovarian cancer risk genes had not previously been implicated in this EOC histotype. Risk genes nominated by chromMAGMA converged on mRNA splicing and transcriptional dysregulation pathways. chromMAGMA is a pipeline that nominates candidate risk genes through a gene regulation-focused approach and helps interpret the biological mechanism of noncoding risk variants for complex diseases.

Molecular Subclasses of Clear Cell Ovarian Carcinoma and Their Impact on Disease Behavior and Outcomes.

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/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 with 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 toward a lower rate of response 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. See related commentary by Lheureux, p. 4838.