Pooling breast cancer datasets has a synergetic effect on classification performance and improves signature stability

Background: Michiels et al. (Lancet 2005; 365: 488-92) employed a resampling strategy to show that the genes identified as predictors of prognosis from resamplings of a single gene expression dataset are highly variable. The genes most frequently identified in the separate resamplings were put forward as a 'gold standard'. On a higher level, breast cancer datasets collected by different institutions can be considered as resamplings from the underlying breast cancer population. The limited overlap between published prognostic signatures confirms the trend of signature instability identified by the resampling strategy. Six breast cancer datasets, totaling 947 samples, all measured on the Affymetrix platform, are currently available. This provides a unique opportunity to employ a substantial dataset to investigate the effects of pooling datasets on classifier accuracy, signature stability and enrichment of functional categories. Results: We show that the resampling strategy produces a suboptimal ranking of genes, which can not be considered to be a 'gold standard'. When pooling breast cancer datasets, we observed a synergetic effect on the classification performance in 73% of the cases. We also observe a significant positive correlation between the number of datasets that is pooled, the validation performance, the number of genes selected, and the enrichment of specific functional categories. In addition, we have evaluated the support for five explanations that have been postulated for the limited overlap of signatures. Conclusion: The limited overlap of current signature genes can be attributed to small sample size. Pooling datasets results in more accurate classification and a convergence of signature genes. We therefore advocate the analysis of new data within the context of a compendium, rather than analysis in isolatio

SIRAC: Supervised Identification of Regions of Aberration in aCGH datasets

<p>Abstract</p> <p>Background</p> <p>Array comparative genome hybridization (aCGH) provides information about genomic aberrations. Alterations in the DNA copy number may cause the cell to malfunction, leading to cancer. Therefore, the identification of DNA amplifications or deletions across tumors may reveal key genes involved in cancer and improve our understanding of the underlying biological processes associated with the disease.</p> <p>Results</p> <p>We propose a supervised algorithm for the analysis of aCGH data and the identification of regions of chromosomal alteration (SIRAC). We first determine the DNA-probes that are important to distinguish the classes of interest, and then evaluate in a systematic and robust scheme if these relevant DNA-probes are closely located, i.e. form a region of amplification/deletion. SIRAC does not need any preprocessing of the aCGH datasets, and requires only few, intuitive parameters.</p> <p>Conclusion</p> <p>We illustrate the features of the algorithm with the use of a simple artificial dataset. The results on two breast cancer datasets show promising outcomes that are in agreement with previous findings, but SIRAC better pinpoints the dissimilarities between the classes of interest.</p

Integrative molecular and functional profiling of ERBB2-amplified breast cancers identifies new genetic dependencies.

Overexpression of the receptor tyrosine kinase ERBB2 (also known as HER2) occurs in around 15% of breast cancers and is driven by amplification of the ERBB2 gene. ERBB2 amplification is a marker of poor prognosis, and although anti-ERBB2-targeted therapies have shown significant clinical benefit, de novo and acquired resistance remains an important problem. Genomic profiling has demonstrated that ERBB2+ve breast cancers are distinguished from ER+ve and 'triple-negative' breast cancers by harbouring not only the ERBB2 amplification on 17q12, but also a number of co-amplified genes on 17q12 and amplification events on other chromosomes. Some of these genes may have important roles in influencing clinical outcome, and could represent genetic dependencies in ERBB2+ve cancers and therefore potential therapeutic targets. Here, we describe an integrated genomic, gene expression and functional analysis to determine whether the genes present within amplicons are critical for the survival of ERBB2+ve breast tumour cells. We show that only a fraction of the ERBB2-amplified breast tumour lines are truly addicted to the ERBB2 oncogene at the mRNA level and display a heterogeneous set of additional genetic dependencies. These include an addiction to the transcription factor gene TFAP2C when it is amplified and overexpressed, suggesting that TFAP2C represents a genetic dependency in some ERBB2+ve breast cancer cell

A comprehensive analysis of prognostic signatures reveals the high predictive capacity of the proliferation, immune response and RNA splicing modules in breast cancer.

INTRODUCTION: Several gene expression signatures have been proposed and demonstrated to be predictive of outcome in breast cancer. In the present article we address the following issues: Do these signatures perform similarly? Are there (common) molecular processes reported by these signatures? Can better prognostic predictors be constructed based on these identified molecular processes? METHODS: We performed a comprehensive analysis of the performance of nine gene expression signatures on seven different breast cancer datasets. To better characterize the functional processes associated with these signatures, we enlarged each signature by including all probes with a significant correlation to at least one of the genes in the original signature. The enrichment of functional groups was assessed using four ontology databases. RESULTS: The classification performance of the nine gene expression signatures is very similar in terms of assigning a sample to either a poor outcome group or a good outcome group. Nevertheless the concordance in classification at the sample level is low, with only 50% of the breast cancer samples classified in the same outcome group by all classifiers. The predictive accuracy decreases with the number of poor outcome assignments given to a sample. The best classification performance was obtained for the group of patients with only good outcome assignments. Enrichment analysis of the enlarged signatures revealed 11 functional modules with prognostic ability. The combination of the RNA-splicing and immune modules resulted in a classifier with high prognostic performance on an independent validation set. CONCLUSIONS: The study revealed that the nine signatures perform similarly but exhibit a large degree of discordance in prognostic group assignment. Functional analyses indicate that proliferation is a common cellular process, but that other functional categories are also enriched and show independent prognostic ability. We provide new evidence of the potentially promising prognostic impact of immunity and RNA-splicing processes in breast cancer.RIGHTS : This article is licensed under the BioMed Central licence at http://www.biomedcentral.com/about/license which is similar to the 'Creative Commons Attribution Licence'. In brief you may : copy, distribute, and display the work; make derivative works; or make commercial use of the work - under the following conditions: the original author must be given credit; for any reuse or distribution, it must be made clear to others what the license terms of this work are

Age-Associated Salivary MicroRNA Biomarkers for Oculopharyngeal Muscular Dystrophy

Small non-coding microRNAs (miRNAs) are involved in the regulation of mRNA stability. Their features, including high stability and secretion to biofluids, make them attractive as potential biomarkers for diverse pathologies. This is the first study reporting miRNA as potential biomarkers for oculopharyngeal muscular dystrophy (OPMD), an adult-onset myopathy. We hypothesized that miRNA that is differentially expressed in affected muscles from OPMD patients is secreted to biofluids and those miRNAs could be used as biomarkers for OPMD. We first identified candidate miRNAs from OPMD-affected muscles and from muscles from an OPMD mouse model using RNA sequencing. We then compared the OPMD-deregulated miRNAs to the literature and, subsequently, we selected a few candidates for expression studies in serum and saliva biofluids using qRT-PCR. We identified 126 miRNAs OPMD-deregulated in human muscles, but 36 deregulated miRNAs in mice only (pFDR < 0.05). Only 15 OPMD-deregulated miRNAs overlapped between the in humans and mouse studies. The majority of the OPMD-deregulated miRNAs showed opposite deregulation direction compared with known muscular dystrophies miRNAs (myoMirs), which are associated. In contrast, similar dysregulation direction was found for 13 miRNAs that are common between OPMD and aging muscles. A significant age-association (p < 0.05) was found for 17 OPMD-deregulated miRNAs (13.4%), whereas in controls, only six miRNAs (1.4%) showed a significant age-association, suggesting that miRNA expression in OPMD is highly age-associated. miRNA expression in biofluids revealed that OPMD-associated deregulation in saliva was similar to that in muscles, but not in serum. The same as in muscle, miRNA expression levels in saliva were also found to be associated with age (p < 0.05). Moreover, the majority of OPMD-miRNAs were found to be associated with dysphagia as an initial symptom. We suggest that levels of specific miRNAs in saliva can mark muscle degeneration in general and dysphagia in OPMDFrench Muscular Dystrophy Association (AFM-Téléthon). Research grant to the eOPMD (European OPMD consortium, V.R. and B.G.M.v.E.)

IHC-based Ki67 as response biomarker to tamoxifen in breast cancer window trials enrolling premenopausal women

Window studies are gaining traction to assess (molecular) changes in short timeframes. Decreased tumor cell positivity for the proliferation marker Ki67 is often used as a proxy for treatment response. Immunohistochemistry (IHC)-based Ki67 on tissue from neo-adjuvant trials was previously reported to be predictive for long-term response to endocrine therapy for breast cancer in postmenopausal women, but none of these trials enrolled premenopausal women. Nonetheless, the marker is being used on this subpopulation. We compared pathologist assessed IHC-based Ki67 in samples from pre- and postmenopausal women in a neo-adjuvant, endocrine therapy focused trial (NCT00738777), randomized between tamoxifen, anastrozole, or fulvestrant. These results were compared with (1) IHC-based Ki67 scoring by AI, (2) mitotic figures, (3) mRNA-based Ki67, (4) five independent gene expression signatures capturing proliferation, and (5) blood levels for tamoxifen and its metabolites as well as estradiol. Upon tamoxifen, IHC-based Ki67 levels were decreased in both pre- and postmenopausal breast cancer patients, which was confirmed using mRNA-based cell proliferation markers. The magnitude of decrease of Ki67 IHC was smaller in pre- versus postmenopausal women. We found a direct relationship between post-treatment estradiol levels and the magnitude of the Ki67 decrease in tumors. These data suggest IHC-based Ki67 may be an appropriate biomarker for tamoxifen response in premenopausal breast cancer patients, but anti-proliferative effect size depends on estradiol levels.</p

Spatial interplay of lymphocytes and fibroblasts in estrogen receptor-positive HER2-negative breast cancer

In estrogen-receptor-positive, HER2-negative (ER+HER2−) breast cancer, higher levels of tumor infiltrating lymphocytes (TILs) are often associated with a poor prognosis and this phenomenon is still poorly understood. Fibroblasts represent one of the most frequent cells in breast cancer and harbor immunomodulatory capabilities. Here, we evaluate the molecular and clinical impact of the spatial patterns of TILs and fibroblast in ER+HER2− breast cancer. We used a deep neural network to locate and identify tumor, TILs, and fibroblasts on hematoxylin and eosin-stained slides from 179 ER+HER2− breast tumors (ICGC cohort) together with a new density estimation analysis to measure the spatial patterns. We clustered tumors based on their spatial patterns and gene set enrichment analysis was performed to study their molecular characteristics. We independently assessed the spatial patterns in a second cohort of ER+HER2− breast cancer (N = 630, METABRIC) and studied their prognostic value. The spatial integration of fibroblasts, TILs, and tumor cells leads to a new reproducible spatial classification of ER+HER2− breast cancer and is linked to inflammation, fibroblast meddling, or immunosuppression. ER+HER2− patients with high TIL did not have a significant improved overall survival (HR = 0.76, P = 0.212), except when they had received chemotherapy (HR = 0.447). A poorer survival was observed for patients with high fibroblasts that did not show a high level of TILs (HR = 1.661, P = 0.0303). Especially spatial mixing of fibroblasts and TILs was associated with a good prognosis (HR = 0.464, P = 0.013). Our findings demonstrate a reproducible pipeline for the spatial profiling of TILs and fibroblasts in ER+HER2− breast cancer and suggest that this spatial interplay holds a decisive role in their cancer-immune interactions

Adult-type granulosa cell tumor of the ovary : a FOXL2-centric disease

Adult-type granulosa cell tumors (aGCTs) account for 90% of malignant ovarian sex cord-stromal tumors and 2-5% of all ovarian cancers. These tumors are usually diagnosed at an early stage and are treated with surgery. However, one-third of patients relapse between 4 and 8 years after initial diagnosis, and there are currently no effective treatments other than surgery for these relapsed patients. As the majority of aGCTs (>95%) harbor a somatic mutation in FOXL2 (c.C402G; p.C134W), the aim of this study was to identify genetic mutations besides FOXL2 C402G in aGCTs that could explain the clinical diversity of this disease. Whole-genome sequencing of 10 aGCTs and their matched normal blood was performed to identify somatic mutations. From this analysis, a custom amplicon-based panel was designed to sequence 39 genes of interest in a validation cohort of 83 aGCTs collected internationally. KMT2D inactivating mutations were present in 10 of 93 aGCTs (10.8%), and the frequency of these mutations was similar between primary and recurrent aGCTs. Inactivating mutations, including a splice site mutation in candidate tumor suppressor WNK2 and nonsense mutations in PIK3R1 and NLRC5, were identified at a low frequency in our cohort. Missense mutations were identified in cell cycle-related genes TP53, CDKN2D, and CDK1. From these data, we conclude that aGCTs are comparatively a homogeneous group of tumors that arise from a limited set of genetic events and are characterized by the FOXL2 C402G mutation. Secondary mutations occur in a subset of patients but do not explain the diverse clinical behavior of this disease. As the FOXL2 C402G mutation remains the main driver of this disease, progress in the development of therapeutics for aGCT would likely come from understanding the functional consequences of the FOXL2 C402G mutation.Peer reviewe

BRCA1-mutated and basal-like breast cancers have similar aCGH profiles and a high incidence of protein truncating TP53 mutations

<p>Abstract</p> <p>Background</p> <p>Basal-like breast cancers (BLBC) are aggressive breast cancers for which, so far, no targeted therapy is available because they typically lack expression of hormone receptors and HER2. Phenotypic features of BLBCs, such as clinical presentation and early age of onset, resemble those of breast tumors from <it>BRCA1</it>-mutation carriers. The genomic instability of <it>BRCA1</it>-mutated tumors can be effectively targeted with DNA-damaging agents and poly-(ADP-ribose) polymerase 1 (PARP1) inhibitors. Molecular similarities between BLBCs and <it>BRCA1</it>-mutated tumors may therefore provide predictive markers for therapeutic response of BLBCs.</p> <p>Methods</p> <p>There are several known molecular features characteristic for <it>BRCA1</it>-mutated breast tumors: 1) increased numbers of genomic aberrations, 2) a distinct pattern of genomic aberrations, 3) a high frequency of <it>TP53 </it>mutations and 4) a high incidence of complex, protein-truncating <it>TP53 </it>mutations. We compared the frequency of <it>TP53 </it>mutations and the pattern and amount of genomic aberrations between <it>BRCA1</it>-mutated breast tumors, BLBCs and luminal breast tumors by <it>TP53 </it>gene sequencing and array-based comparative genomics hybridization (aCGH) analysis.</p> <p>Results</p> <p>We found that the high incidence of protein truncating <it>TP53 </it>mutations and the pattern and amount of genomic aberrations specific for BRCA1-mutated breast tumors are also characteristic for BLBCs and different from luminal breast tumors.</p> <p>Conclusions</p> <p>Complex, protein truncating TP53 mutations in BRCA1-mutated tumors may be a direct consequence of genomic instability caused by BRCA1 loss, therefore, the presence of these types of TP53 mutations in sporadic BLBCs might be a hallmark of BRCAness and a potential biomarker for sensitivity to PARP inhibition. Also, our data suggest that a small subset of genomic regions may be used to identify BRCA1-like BLBCs. BLBCs share molecular features that were previously found to be specific for BRCA1-mutated breast tumors. These features might be useful for the identification of tumors with increased sensitivity to (high-dose or dose-dense) alkylating agents and PARP inhibitors.</p

Carcinogenic Effects in a Phenylketonuria Mouse Model

Phenylketonuria (PKU) is a metabolic disorder caused by impaired phenylalanine hydroxylase (PAH). This condition results in hyperphenylalaninemia and elevated levels of abnormal phenylalanine metabolites, among which is phenylacetic acid/phenylacetate (PA). In recent years, PA and its analogs were found to have anticancer activity against a variety of malignancies suggesting the possibility that PKU may offer protection against cancer through chronically elevated levels of PA. We tested this hypothesis in a genetic mouse model of PKU (PAHenu2) which has a biochemical profile that closely resembles that of human PKU. Plasma levels of phenylalanine in homozygous (HMZ) PAHenu2 mice were >12-fold those of heterozygous (HTZ) littermates while tyrosine levels were reduced. Phenylketones, including PA, were also markedly elevated to the range seen in the human disease. Mice were subjected to 7,12 dimethylbenz[a]anthracene (DMBA) carcinogenesis, a model which is sensitive to the anticancer effects of the PA derivative 4-chlorophenylacetate (4-CPA). Tumor induction by DMBA was not significantly different between the HTZ and HMZ mice, either in total tumor development or in the type of cancers that arose. HMZ mice were then treated with 4-CPA as positive controls for the anticancer effects of PA and to evaluate its possible effects on phenylalanine metabolism in PKU mice. 4-CPA had no effect on the plasma concentrations of phenylalanine, phenylketones, or tyrosine. Surprisingly, the HMZ mice treated with 4-CPA developed an unexplained neuromuscular syndrome which precluded its use in these animals as an anticancer agent. Together, these studies support the use of PAHenu2 mice as a model for studying human PKU. Chronically elevated levels of PA in the PAHenu2 mice were not protective against cancer