It has been shown that ovarian carcinoma subtypes are distinct pathologic entities with differing prognostic and therapeutic implications. Histotyping by pathologists has good reproducibility, but occasional cases are challenging and require immunohistochemistry and subspecialty consultation. Motivated by the need for more accurate and reproducible diagnoses and to facilitate pathologists' workflow, we propose an automatic framework for ovarian carcinoma classification.

We report a case of early-onset pancreatic ductal adenocarcinoma in a patient harboring biallelic MUTYH germline mutations, whose tumor featured somatic mutational signatures consistent with defective MUTYH-mediated base excision repair and the associated driver KRAS transversion mutation p.Gly12Cys. Analysis of an additional 730 advanced cancer cases (N = 731) was undertaken to determine whether the mutational signatures were also present in tumors from germline MUTYH heterozygote carriers or if instead the signatures were only seen in those with biallelic loss of function. We identified two patients with breast cancer each carrying a pathogenic germline MUTYH variant with a somatic MUTYH copy loss leading to the germline variant being homozygous in the tumor and demonstrating the same somatic signatures. Our results suggest that monoallelic inactivation of MUTYH is not sufficient for C:G>A:T transversion signatures previously linked to MUTYH deficiency to arise (N = 9), but that biallelic complete loss of MUTYH function can cause such signatures to arise even in tumors not classically seen in MUTYH-associated polyposis (N = 3). Although defective MUTYH is not the only determinant of these signatures, MUTYH germline variants may be present in a subset of patients with tumors demonstrating elevated somatic signatures possibly suggestive of MUTYH deficiency (e.g., COSMIC Signature 18, SigProfiler SBS18/SBS36, SignatureAnalyzer SBS18/SBS36).

Genomic analysis of cancer tissues is an essential aspect of personalized oncology treatment. Though it has been suggested that formalin fixation of patient tissues may be suboptimal for molecular studies, this tissue processing approach remains the industry standard. Therefore clinical molecular laboratories must be able to work with formalin fixed, paraffin embedded (FFPE) material. This study examines the effects of pre-analytic variables introduced by routine pathology processing on specimens used for clinical reports produced by next-generation sequencing technology. Tissue resected from three colorectal cancer patients was subjected to 2, 15, 24, and 48 hour fixation times in neutral buffered formalin. DNA was extracted from all tissues twice, once with uracil-N-glycosylase (UNG) treatment to counter deamination effects, and once without. Of note, deamination events at methylated cytosine, as found at CpG sites, remains unaffected by UNG. After extraction a two-step PCR targeted sequencing method was performed using the Illumina MiSeq and the data was analyzed via a custom-built bioinformatics pipeline, including filtration of reads with mapping quality <30. A larger baseline group of samples (n = 20) was examined to establish if there was a sample performance difference between the two DNA extraction methods, with/without UNG treatment. There was no statistical difference between sequencing performance of the two extraction methods when comparing read counts (raw, mapped, and filtered) and read quality (% mapped, % filtered). Analyzing mutation type, there was no significant difference between mutation calls until the 48 hour fixation treatment. At 48 hours there is a significant increase in C/G->T/A mutations that is not represented in DNA treated with UNG. This suggests these errors may be due to deamination events triggered by a longer fixation time. However the allelic frequency of these events remained below the limit of detection for reportable mutations in this assay (<2%). We do however recommend that suspected intratumoral heterogeneity events be verified by re-sequencing the same FFPE block.

Transitions between epithelial and mesenchymal cellular states (EMT/MET) contribute to cancer progression. We hypothesize that EMT followed by MET promotes cell population heterogeneity, favouring tumour growth. We developed an EMT model by on and off exposure of epithelial EpH4 cells (E-cells) to TGFβ1 that mimics phenotypic EMT (M-cells) and MET. We aimed at understanding whether phenotypic MET is accompanied by molecular and functional reversion back to epithelia by using RNA sequencing, immunofluorescence (IF), proliferation, wound healing, focus formation and mamosphere formation assays as well as cell xenografts in nude mice. Phenotypic reverted epithelial cells (RE-cells) obtained after MET induction presented epithelial morphologies and proliferation rates resembling E cells. However, the RE transcriptomic profile and IF staining of epithelial and mesenchymal markers revealed a uniquely heterogeneous mixture of cell subpopulations with a high self-renewal ability. RE cell heterogeneity was stably maintained for long periods after TGFβ1 removal both in vitro and in large tumours derived from the nude mice. Overall, we show that phenotypic reverted epithelial cells (RE cells) do not return to the molecular and functional epithelial state and present mesenchymal features related to aggressiveness and cellular heterogeneity that favour tumour growth in vivo. This work strengthens epithelial cell reprogramming and cellular heterogeneity fostered by inflammatory cues as a tumour growth-promoting factor in vivo.

HNF1B is overexpressed in clear cell epithelial ovarian cancer, and we observed epigenetic silencing in serous epithelial ovarian cancer, leading us to hypothesize that variation in this gene differentially associates with epithelial ovarian cancer risk according to histological subtype. Here we comprehensively map variation in HNF1B with respect to epithelial ovarian cancer risk and analyse DNA methylation and expression profiles across histological subtypes. Different single-nucleotide polymorphisms associate with invasive serous (rs7405776 odds ratio (OR)=1.13, P=3.1 × 10(-10)) and clear cell (rs11651755 OR=0.77, P=1.6 × 10(-8)) epithelial ovarian cancer. Risk alleles for the serous subtype associate with higher HNF1B-promoter methylation in these tumours. Unmethylated, expressed HNF1B, primarily present in clear cell tumours, coincides with a CpG island methylator phenotype affecting numerous other promoters throughout the genome. Different variants in HNF1B associate with risk of serous and clear cell epithelial ovarian cancer; DNA methylation and expression patterns are also notably distinct between these subtypes. These findings underscore distinct mechanisms driving different epithelial ovarian cancer histological subtypes.

Germline CDH1 point or small frameshift mutations can be identified in 30-50% of hereditary diffuse gastric cancer (HDGC) families. We hypothesized that CDH1 genomic rearrangements would be found in HDGC and identified 160 families with either two gastric cancers in first-degree relatives and with at least one diffuse gastric cancer (DGC) diagnosed before age 50, or three or more DGC in close relatives diagnosed at any age. Sixty-seven carried germline CDH1 point or small frameshift mutations. We screened germline DNA from the 93 mutation negative probands for large genomic rearrangements by Multiplex Ligation-Dependent Probe Amplification. Potential deletions were validated by RT-PCR and breakpoints cloned using a combination of oligo-CGH-arrays and long-range-PCR. In-silico analysis of the CDH1 locus was used to determine a potential mechanism for these rearrangements. Six of 93 (6.5%) previously described mutation negative HDGC probands, from low GC incidence populations (UK and North America), carried genomic deletions (UK and North America). Two families carried an identical deletion spanning 193 593 bp, encompassing the full CDH3 sequence and CDH1 exons 1 and 2. Other deletions affecting exons 1, 2, 15 and/or 16 were identified. The statistically significant over-representation of Alus around breakpoints indicates it as a likely mechanism for these deletions. When all mutations and deletions are considered, the overall frequency of CDH1 alterations in HDGC is approximately 46% (73/160). CDH1 large deletions occur in 4% of HDGC families by mechanisms involving mainly non-allelic homologous recombination in Alu repeat sequences. As the finding of pathogenic CDH1 mutations is useful for management of HDGC families, screening for deletions should be offered to at-risk families.

Gastric cancer is among the leading causes of cancer-related deaths worldwide. While heritable forms of gastric cancer are relatively rare, identifying the genes responsible for such cases can inform diagnosis and treatment for both hereditary and sporadic cases of gastric cancer. Mutations in the E-cadherin gene, CDH1, account for 40% of the most common form of familial gastric cancer (FGC), hereditary diffuse gastric cancer (HDGC). The genes responsible for the remaining forms of FGC are currently unknown. Here we examined a large family from Maritime Canada with FGC without CDH1 mutations, and identified a germline coding variant (p.P946L) in mitogen-activated protein kinase kinase kinase 6 (MAP3K6). Based on conservation, predicted pathogenicity and a known role of the gene in cancer predisposition, MAP3K6 was considered a strong candidate and was investigated further. Screening of an additional 115 unrelated individuals with non-CDH1 FGC identified the p.P946L MAP3K6 variant, as well as four additional coding variants in MAP3K6 (p.F849Sfs*142, p.P958T, p.D200Y and p.V207G). A somatic second-hit variant (p.H506Y) was present in DNA obtained from one of the tumor specimens, and evidence of DNA hypermethylation within the MAP3K6 gene was observed in DNA from the tumor of another affected individual. These findings, together with previous evidence from mouse models that MAP3K6 acts as a tumor suppressor, and studies showing the presence of somatic mutations in MAP3K6 in non-hereditary gastric cancers and gastric cancer cell lines, point towards MAP3K6 variants as a predisposing factor for FGC.

Dies1/VISTA induces embryonic stem-cell differentiation, via BMP-pathway, but also acts as inflammation regulator and immune-response modulator. Dies1 inhibition in a melanoma-mouse model led to increased tumour-infiltrating T-cells and decreased tumour growth, emphasizing Dies1 relevance in tumour-microenvironment. Dies1 is involved in cell de/differentiation, inflammation and cancer processes, which mimic those associated with Epithelial-to-Mesenchymal-Transition (EMT). Despite this axis linking Dies1 with EMT and cancer, its expression, modulation and relevance in these contexts is unknown. To address this, we analysed Dies1 expression, its regulation by promoter-methylation and miR-125a-5p overexpression, and its association with BMP-pathway downstream-effectors, in a TGFβ1-induced EMT-model, cancer cell-lines and primary samples. We detected promoter-methylation as a mechanism controlling Dies1 expression in our EMT-model and in several cancer cell-lines. We showed that the relationship between Dies1 expression and BMP-pathway effectors observed in the EMT-model, was not present in all cell-lines, suggesting that Dies1 has other cell-specific effectors, beyond the BMP-pathway. We further demonstrated that: Dies1 expression loss is a recurrent event in GC, caused by promoter methylation and/or miR-125a-5p overexpression and; GC-microenvironment myofibroblasts overexpress Dies1. Our findings highlight Dies1 as a novel player in GC, with distinct roles within tumour cells and in the tumour-microenvironment.

Chordoma is a malignant, often incurable bone tumour showing notochordal differentiation. Here, we defined the somatic driver landscape of 104 cases of sporadic chordoma. We reveal somatic duplications of the notochordal transcription factor brachyury (T) in up to 27% of cases. These variants recapitulate the rearrangement architecture of the pathogenic germline duplications of T that underlie familial chordoma. In addition, we find potentially clinically actionable PI3K signalling mutations in 16% of cases. Intriguingly, one of the most frequently altered genes, mutated exclusively by inactivating mutation, was LYST (10%), which may represent a novel cancer gene in chordoma.Chordoma is a rare often incurable malignant bone tumour. Here, the authors investigate driver mutations of sporadic chordoma in 104 cases, revealing duplications in notochordal transcription factor brachyury (T), PI3K signalling mutations, and mutations in LYST, a potential novel cancer gene in chordoma.

Epithelial ovarian cancer (EOC) is a heterogeneous group of malignancies, including high-grade serous ovarian cancer (HGSC). HGSC is often diagnosed at advanced stages and is linked to TP53 variants. While BRCA variants elevate risk, most HGSC cases occur in individuals without known genetic variants, necessitating prevention strategies for people without known high-risk genetic variants. Effective prevention programs are also needed due to the lack of traditional screening options. An emerging primary prevention strategy is opportunistic salpingectomy, which involves removing fallopian tubes during another planned pelvic surgery. Opportunistic salpingectomy offers a safe and cost-effective preventative option that is gaining global adoption. With the publication of the first cohort study of patients who underwent salpingectomy, specifically for cancer prevention, attention has turned to broadening opportunities for salpingectomy in addition to more targeted approaches. Prevention opportunities are promising with increasing adoption of salpingectomy and the increased understanding of the etiology of the distinct histotypes of ovarian cancer. Yet, further research on targeted risk-reducing salpingectomy with thoughtful consideration of equity is necessary to reduce death and suffering from ovarian cancer.

We present a novel hierarchical Bayes statistical model, xseq, to systematically quantify the impact of somatic mutations on expression profiles. We establish the theoretical framework and robust inference characteristics of the method using computational benchmarking. We then use xseq to analyse thousands of tumour data sets available through The Cancer Genome Atlas, to systematically quantify somatic mutations impacting expression profiles. We identify 30 novel cis-effect tumour suppressor gene candidates, enriched in loss-of-function mutations and biallelic inactivation. Analysis of trans-effects of mutations and copy number alterations with xseq identifies mutations in 150 genes impacting expression networks, with 89 novel predictions. We reveal two important novel characteristics of mutation impact on expression: (1) patients harbouring known driver mutations exhibit different downstream gene expression consequences; (2) expression patterns for some mutations are stable across tumour types. These results have critical implications for identification and interpretation of mutations with consequent impact on transcription in cancer.

We previously identified associations with ovarian cancer outcome at five genetic loci. To identify putatively causal genetic variants and target genes, we prioritized two ovarian outcome loci (1q22 and 19p12) for further study. Bioinformatic and functional genetic analyses indicated that MEF2D and ZNF100 are targets of candidate outcome variants at 1q22 and 19p12, respectively. At 19p12, the chromatin interaction of a putative regulatory element with the ZNF100 promoter region correlated with candidate outcome variants. At 1q22, putative regulatory elements enhanced MEF2D promoter activity and haplotypes containing candidate outcome variants modulated these effects. In a public dataset, MEF2D and ZNF100 expression were both associated with ovarian cancer progression-free or overall survival time. In an extended set of 6,162 epithelial ovarian cancer patients, we found that functional candidates at the 1q22 and 19p12 loci, as well as other regional variants, were nominally associated with patient outcome; however, no associations reached our threshold for statistical significance (p<1×10-5). Larger patient numbers will be needed to convincingly identify any true associations at these loci.

The aim of this study was to determine the prevalence of somatic and germline pathogenic variants (PVs) in high-grade serous cancer (HGSC) and to demonstrate the technical feasibility and effectiveness of a large-scale, population-based tumor testing program. It involved a retrospective review of genetic test results in 600 consecutive HGSC tumor samples and a subsequent comparison of germline and tumor results in a subset of 200 individuals. Tumor testing was successful in 95% of samples (570/600) with at least one BRCA1/2 PV identified in 16% (93/570) of cases. Among the 200 paired cases, BRCA1/2 PVs were detected in 38 tumors (19%); 58% were somatic (22/38); and 42% were germline (16/38). There was 100% concordance between germline and tumor test results. This is the largest series of BRCA1/2 testing in HGSC (tumor-only and paired cohorts), reported to date, and our data show that an effectively designed and validated population-based tumor testing program can be used to determine both treatment eligibility and hereditary cancer risk.

Ovarian low-grade serous carcinoma (LGSC) has fewer mutations than ovarian high-grade serous carcinoma (HGSC) and a less aggressive clinical course. However, an overwhelming majority of LGSC patients do not respond to conventional chemotherapy resulting in a poor long-term prognosis comparable to women diagnosed with HGSC. KRAS and BRAF mutations are common in LGSC, leading to clinical trials targeting the MAPK pathway. We assessed the stability of targetable somatic mutations over space and/or time in LGSC, with a view to inform stratified treatment strategies and clinical trial design.

Gene fusions created by somatic genomic rearrangements are known to play an important role in the onset and development of some cancers, such as lymphomas and sarcomas. RNA-Seq (whole transcriptome shotgun sequencing) is proving to be a useful tool for the discovery of novel gene fusions in cancer transcriptomes. However, algorithmic methods for the discovery of gene fusions using RNA-Seq data remain underdeveloped. We have developed deFuse, a novel computational method for fusion discovery in tumor RNA-Seq data. Unlike existing methods that use only unique best-hit alignments and consider only fusion boundaries at the ends of known exons, deFuse considers all alignments and all possible locations for fusion boundaries. As a result, deFuse is able to identify fusion sequences with demonstrably better sensitivity than previous approaches. To increase the specificity of our approach, we curated a list of 60 true positive and 61 true negative fusion sequences (as confirmed by RT-PCR), and have trained an adaboost classifier on 11 novel features of the sequence data. The resulting classifier has an estimated value of 0.91 for the area under the ROC curve. We have used deFuse to discover gene fusions in 40 ovarian tumor samples, one ovarian cancer cell line, and three sarcoma samples. We report herein the first gene fusions discovered in ovarian cancer. We conclude that gene fusions are not infrequent events in ovarian cancer and that these events have the potential to substantially alter the expression patterns of the genes involved; gene fusions should therefore be considered in efforts to comprehensively characterize the mutational profiles of ovarian cancer transcriptomes.

Phosphorylated Akt (P-Akt) is an attractive molecular target because it contributes to the development of breast cancer and confers resistance to conventional therapies. Akt also serves as a signalling intermediate for receptors such as human epidermal growth factor receptor (HER)-2, which is overexpressed in 30% of breast cancers; therefore, inhibitors to this pathway are being sought. New celecoxib analogues reportedly inhibit P-Akt in prostate cancer cells. We therefore examined the potential of these compounds in the treatment of breast cancer. The analogues were characterized in MDA-MB-453 cells because they overexpress HER-2 and have very high levels of P-Akt.

The International Endocervical Adenocarcinoma Criteria and Classification (IECC) separated endocervical adenocarcinomas into human papillomavirus (HPV) associated (HPVA) and non-HPV-associated (NHPVA) categories by morphology alone. Our primary objective was to assess the accuracy of HPV prediction by the IECC system compared to p16 immunohistochemistry and HPV RNA in-situ hybridization (RISH). Our secondary goal was to directly compare p16 and HPV RISH concordance.

MicroRNAs (miRNA) are 20 approximately 25 nucleotide non-coding RNAs that inhibit the translation of targeted mRNA, and they have been implicated in the development of human malignancies. High grade serous ovarian carcinomas, the most common and lethal subtype of ovarian cancer, can occur sporadically or in the setting of BRCA1/2 syndromes. Little is known regarding the miRNA expression profiles of high grade serous carcinoma in relation to BRCA1/2 status, and compared to normal tubal epithelium, the putative tissue of origin for high grade serous carcinomas.

Breast cancer accounts for approximately 15% of all cancer deaths. Currently, axillary nodal status is the most reliable prognostic indicator for breast cancer. Tumor size and histological grade are used to stage breast cancer. Estrogen receptor/progesterone receptor (ER/PR) and HER-2/neu status are useful in predicting patient survival and relapse. Ki67, an indicator of proliferative activity, also correlates well with prognosis. Connexin proteins form gap junction channels, permitting intercellular exchange of ions and small molecules. Reduced connexin protein levels and impaired gap junctional intercellular communication are associated with tumor phenotypes. This study investigated the prognostic value of connexin proteins as breast cancer markers. Tissue microarrays, containing 438 cases of invasive breast carcinoma, were stained with Cx26, Cx32, and Cx43 antibodies. The degree of connexin immunoreactivity was determined and then correlated with patient outcome, tumor grade, tumor size, lymph node status, and immunohistochemical markers, such as p53, ER/PR status, Ki67 and c-erbB-2 expression. Cx26, Cx32, or Cx43 did not correlate well with tumor grade, tumor size, p53 or c-erbB-2 status. There was an inverse correlation between Cx32 and lymph node status (P <0.05) and a positive correlation between Cx43 and PR status (P <0.01). Cx32 and Cx43 correlated positively with ER status (P <0.01). Cx43 correlated negatively with Ki67 expression (P <0.01). Cx26, Cx32, and Cx43 did not correlate with patient outcome. Based on our observations in this study, connexin proteins do not appear to be reliable indicators of breast cancer prognosis.

To identify molecular subclasses of clear cell ovarian carcinoma (CCOC) and assess their impact on clinical presentation and outcomes.