Background

Most pelvic high-grade serous (HGS) carcinomas have been proposed to arise from tubal primaries that progress rapidly to advanced disease. However, the temporal sequence of ovarian and peritoneal metastases is not well characterized.

Methods

To establish the sequence of metastases, phylogenetic relationships among the ovarian and peritoneal carcinomas were determined from single-nucleotide variations (SNVs) in nine tumor regions from each patient with pelvic HGS carcinomas. Somatic SNVs from each tumor sample were used to reconstruct phylogenies of samples from each patient. Variant allele frequencies were used to reconstruct subclone phylogenies in each tumor sample.

Results

We show that pelvic HGS carcinomas are highly heterogeneous, only sharing less than 4% of somatic SNVs among all nine carcinoma implants in one patient. TP53 mutations are found in all nine carcinoma implants in each patient. The phylogenetic analyses reveal that peritoneal metastases arose from early branching events that preceded branching events for ovarian carcinomas in some patients. Finally, subclone phylogenies indicate the presence of multiple subclones at each tumor implant and early tumor clones in peritoneal implants.

Conclusion

The genetic evidence that peritoneal implants arose before or concurrently with ovarian implants is consistent with the emerging concept of the extra-ovarian origin of pelvic HGS cancer. Our results challenge the concept of stepwise spatial progression from the fallopian primary to ovarian carcinomas to peritoneal dissemination and suggest an alternative progression model where peritoneal spreading of early clones occurs before or in parallel with ovarian metastases.

Background

Both genetic and epigenetic factors influence the development and progression of epithelial ovarian cancer (EOC). However, there is an incomplete understanding of the interrelationship between these factors and the extent to which they interact to impact disease risk. In the present study, we aimed to gain insight into this relationship by identifying DNA methylation marks that are candidate mediators of ovarian cancer genetic risk.

Methods

We used 214 cases and 214 age-matched controls from the Mayo Clinic Ovarian Cancer Study. Pretreatment, blood-derived DNA was profiled for genome-wide methylation (Illumina Infinium HumanMethylation27 BeadArray) and single nucleotide polymorphisms (SNPs, Illumina Infinium HD Human610-Quad BeadArray). The Causal Inference Test (CIT) was implemented to distinguish CpG sites that mediate genetic risk, from those that are consequential or independently acted on by genotype.

Results

Controlling for the estimated distribution of immune cells and other key covariates, our initial epigenome-wide association analysis revealed 1,993 significantly differentially methylated CpGs that between cases and controls (FDR, q < 0.05). The relationship between methylation and case-control status for these 1,993 CpGs was found to be highly consistent with the results of previously published, independent study that consisted of peripheral blood DNA methylation signatures in 131 pretreatment cases and 274 controls. Implementation of the CIT test revealed 17 CpG/SNP pairs, comprising 13 unique CpGs and 17 unique SNPs, which represent potential methylation-mediated relationships between genotype and EOC risk. Of these 13 CpGs, several are associated with immune related genes and genes that have been previously shown to exhibit altered expression in the context of cancer.

Conclusions

These findings provide additional insight into EOC etiology and may serve as novel biomarkers for EOC susceptibility.

Background

The mechanisms of recurrence have been under-studied in rare histologies of invasive epithelial ovarian cancer (EOC) (endometrioid, clear cell, mucinous, and low-grade serous). We hypothesised the existence of an expression signature predictive of outcome in the rarer histologies.

Methods

In split discovery and validation analysis of 131 Mayo Clinic EOC cases, we used clustering to determine clinically relevant transcriptome classes using microarray gene expression measurements. The signature was validated in 967 EOC tumours (91 rare histological subtypes) with recurrence information.

Results

We found two validated transcriptome classes associated with progression-free survival (PFS) in the Mayo Clinic EOC cases (P=8.24 × 10(-3)). This signature was further validated in the public expression data sets involving the rare EOC histologies, where these two classes were also predictive of PFS (P=1.43 × 10(-3)). In contrast, the signatures were not predictive of PFS in the high-grade serous EOC cases. Moreover, genes upregulated in Class-1 (with better outcome) were showed enrichment in steroid hormone biosynthesis (false discovery rate, FDR=0.005%) and WNT signalling pathway (FDR=1.46%); genes upregulated in Class-2 were enriched in cell cycle (FDR=0.86%) and toll-like receptor pathways (FDR=2.37%).

Conclusions

These findings provide important biological insights into the rarer EOC histologies that may aid in the development of targeted treatment options for the rarer histologies.

Background

It is unclear whether the transcriptional subtypes of high grade serous ovarian cancer (HGSOC) apply to high grade clear cell (HGCCOC) or high grade endometrioid ovarian cancer (HGEOC). We aim to delineate transcriptional profiles of HGCCOCs and HGEOCs.

Methods

We used Agilent microarrays to determine gene expression profiles of 276 well annotated ovarian cancers (OCs) including 37 HGCCOCs and 66 HGEOCs. We excluded low grade OCs as these are known to be distinct molecular entities. We applied the prespecified TCGA and CLOVAR gene signatures using consensus non-negative matrix factorization (NMF).

Results

We confirm the presence of four TCGA transcriptional subtypes and their significant prognostic relevance (p<0.001) across all three histological subtypes (HGSOC, HGCCOC and HGEOCs). However, we also demonstrate that 22/37 (59%) HGCCOCs and 30/67 (45%) HGEOCs form 2 additional separate clusters with distinct gene signatures. Importantly, of the HGCCOC and HGEOCs that clustered separately 62% and 65% were early stage (FIGO I/II), respectively. These finding were confirmed using the reduced CLOVAR gene set for classification where most early stage HGCCOCs and HGEOCs formed a distinct cluster of their own. When restricting the analysis to the four TCGA signatures (ssGSEA or NMF with CLOVAR genes) most early stage HGCCOCs and HGEOC were assigned to the differentiated subtype.

Conclusions

Using transcriptional profiling the current study suggests that HGCCOCs and HGEOCs of advanced stage group together with HGSOCs. However, HGCCOCs and HGEOCs of early disease stages may have distinct transcriptional signatures similar to those seen in their low grade counterparts.

Purpose

Ovarian cancer has a high recurrence and mortality rate. A barrier to improved outcomes includes a lack of accurate models for preclinical testing of novel therapeutics.

Experimental design

Clinically relevant, patient-derived tumorgraft models were generated from sequential patients and the first 168 engrafted models are described. Fresh ovarian, primary peritoneal, and fallopian tube carcinomas were collected at the time of debulking surgery and injected intraperitoneally into severe combined immunodeficient mice.

Results

Tumorgrafts demonstrated a 74% engraftment rate with microscopic fidelity of primary tumor characteristics. Low-passage tumorgrafts also showed comparable genomic aberrations with the corresponding primary tumor and exhibit gene set enrichment of multiple ovarian cancer molecular subtypes, similar to patient tumors. Importantly, each of these tumorgraft models is annotated with clinical data and for those that have been tested, response to platinum chemotherapy correlates with the source patient.

Conclusions

Presented herein is the largest known living tumor bank of patient-derived, ovarian tumorgraft models that can be applied to the development of personalized cancer treatment.

To reveal biologic mechanisms underlying clinical outcome of high-grade serous (HGS) epithelial ovarian carcinomas (EOC), we evaluated the association between tumor epigenetic changes and time to recurrence (TTR). We assessed methylation at approximately 450,000 genome-wide CpGs in tumors of 337 Mayo Clinic (Rochester, MN) patients. Semi-supervised clustering of discovery (n=168) and validation (n=169) sets was used to determine clinically relevant methylation classes. Clustering identified two methylation classes based on 60 informative CpGs, which differed in TTR in the validation set [R vs. L class, P=2.9×10(-3), HR=0.52; 95% confidence interval (CI), 0.34-0.80]. Follow-up analyses considered genome-wide tumor mRNA expression (n=104) and CD8 T-cell infiltration (n=89) in patient subsets. Hypomethylation of CpGs located in 6p21.3 in the R class associated with cis upregulation of genes enriched in immune response processes (TAP1, PSMB8, PSMB9, HLA-DQB1, HLA-DQB2, HLA-DMA, and HLA-DOA), increased CD8 T-cell tumor infiltration (P=7.6×10(-5)), and trans-regulation of genes in immune-related pathways (P=1.6×10(-32)). This is the most comprehensive assessment of clinical outcomes with regard to epithelial ovarian carcinoma tumor methylation to date. Collectively, these results suggest that an epigenetically mediated immune response is a predictor of recurrence and, possibly, treatment response for HGS EOC.

Molecular classification of high-grade serous ovarian cancer (HGSOC) using transcriptional profiling has proven to be complex and difficult to validate across studies. We determined gene expression profiles of 174 well-annotated HGSOCs and demonstrate prognostic significance of the prespecified TCGA Network gene signatures. Furthermore, we confirm the presence of four HGSOC transcriptional subtypes using a de novo classification. Survival differed statistically significantly between de novo subtypes (log rank, P = .006) and was the best for the immunoreactive-like subtype, but statistically significantly worse for the proliferative- or mesenchymal-like subtypes (adjusted hazard ratio = 1.89, 95% confidence interval = 1.18 to 3.02, P = .008, and adjusted hazard ratio = 2.45, 95% confidence interval = 1.43 to 4.18, P = .001, respectively). More prognostic information was provided by the de novo than the TCGA classification (Likelihood Ratio tests, P = .003 and P = .04, respectively). All statistical tests were two-sided. These findings were replicated in an external data set of 185 HGSOCs and confirm the presence of four prognostically relevant molecular subtypes that have the potential to guide therapy decisions.

To determine early somatic changes in high-grade serous ovarian cancer (HGSOC), we performed whole genome sequencing on a rare collection of 16 low stage HGSOCs. The majority showed extensive structural alterations (one had an ultramutated profile), exhibited high levels of p53 immunoreactivity, and harboured a TP53 mutation, deletion or inactivation. BRCA1 and BRCA2 mutations were observed in two tumors, with nine showing evidence of a homologous recombination (HR) defect. Combined Analysis with The Cancer Genome Atlas (TCGA) indicated that low and late stage HGSOCs have similar mutation and copy number profiles. We also found evidence that deleterious TP53 mutations are the earliest events, followed by deletions or loss of heterozygosity (LOH) of chromosomes carrying TP53, BRCA1 or BRCA2. Inactivation of HR appears to be an early event, as 62.5% of tumours showed a LOH pattern suggestive of HR defects. Three tumours with the highest ploidy had little genome-wide LOH, yet one of these had a homozygous somatic frame-shift BRCA2 mutation, suggesting that some carcinomas begin as tetraploid then descend into diploidy accompanied by genome-wide LOH. Lastly, we found evidence that structural variants (SV) cluster in HGSOC, but are absent in one ultramutated tumor, providing insights into the pathogenesis of low stage HGSOC.

Ovarian cancer is a clinically and molecularly heterogeneous disease. The driving forces behind this variability are unknown. Here, we report wide variation in the expression of the DNA cytosine deaminase APOBEC3B, with elevated expression in the majority of ovarian cancer cell lines (three SDs above the mean of normal ovarian surface epithelial cells) and high-grade primary ovarian cancers. APOBEC3B is active in the nucleus of several ovarian cancer cell lines and elicits a biochemical preference for deamination of cytosines in 5'-TC dinucleotides. Importantly, examination of whole-genome sequence from 16 ovarian cancers reveals that APOBEC3B expression correlates with total mutation load as well as elevated levels of transversion mutations. In particular, high APOBEC3B expression correlates with C-to-A and C-to-G transversion mutations within 5'-TC dinucleotide motifs in early-stage high-grade serous ovarian cancer genomes, suggesting that APOBEC3B-catalyzed genomic uracil lesions are further processed by downstream DNA "repair" enzymes including error-prone translesion polymerases. These data identify a potential role for APOBEC3B in serous ovarian cancer genomic instability.

Although ovarian cancer is the most lethal of gynecologic malignancies, wide variation in outcome following conventional therapy continues to exist. The presence of tumor-infiltrating regulatory T cells (Tregs) has a role in outcome of this disease, and a growing body of data supports the existence of inherited prognostic factors. However, the role of inherited variants in genes encoding Treg-related immune molecules has not been fully explored. We analyzed expression quantitative trait loci (eQTL) and sequence-based tagging single nucleotide polymorphisms (tagSNPs) for 54 genes associated with Tregs in 3,662 invasive ovarian cancer cases. With adjustment for known prognostic factors, suggestive results were observed among rarer histological subtypes; poorer survival was associated with minor alleles at SNPs in RGS1 (clear cell, rs10921202, p=2.7×10(-5)), LRRC32 and TNFRSF18/TNFRSF4 (mucinous, rs3781699, p=4.5×10(-4), and rs3753348, p=9.0×10(-4), respectively), and CD80 (endometrioid, rs13071247, p=8.0×10(-4)). Fo0r the latter, correlative data support a CD80 rs13071247 genotype association with CD80 tumor RNA expression (p=0.006). An additional eQTL SNP in CD80 was associated with shorter survival (rs7804190, p=8.1×10(-4)) among all cases combined. As the products of these genes are known to affect induction, trafficking, or immunosuppressive function of Tregs, these results suggest the need for follow-up phenotypic studies.