Cis-regulatory mutations associate with transcriptional and post-transcriptional deregulation of gene regulatory programs in cancers

Most cancer alterations occur in the noncoding portion of the human genome, where regulatory regions control gene expression. The discovery of noncoding mutations altering the cells' regulatory programs has been limited to few examples with high recurrence or high functional impact. Here, we show that transcription factor binding sites (TFBSs) have similar mutation loads to those in protein-coding exons. By combining cancer somatic mutations in TFBSs and expression data for protein-coding and miRNA genes, we evaluate the combined effects of transcriptional and post-transcriptional alterations on the regulatory programs in cancers. The analysis of seven TCGA cohorts culminates with the identification of protein-coding and miRNA genes linked to mutations at TFBSs that are associated with a cascading trans-effect deregulation on the cells' regulatory programs. Our analyses of cis-regulatory mutations associated with miRNAs recurrently predict 12 mature miRNAs (derived from 7 precursors) associated with the deregulation of their target gene networks. The predictions are enriched for cancer-associated protein-coding and miRNA genes and highlight cis-regulatory mutations associated with the dysregulation of key pathways associated with carcinogenesis. By combining transcriptional and post-transcriptional regulation of gene expression, our method predicts cis-regulatory mutations related to the dysregulation of key gene regulatory networks in cancer patients.

DNA copy number motifs are strong and independent predictors of survival in breast cancer

Somatic copy number alterations are a frequent sign of genome instability in cancer. A precise characterization of the genome architecture would reveal underlying instability mechanisms and provide an instrument for outcome prediction and treatment guidance. Here we show that the local spatial behavior of copy number profiles conveys important information about this architecture. Six filters were defined to characterize regional traits in copy number profiles, and the resulting Copy Aberration Regional Mapping Analysis (CARMA) algorithm was applied to tumors in four breast cancer cohorts (n = 2919). The derived motifs represent a layer of information that complements established molecular classifications of breast cancer. A score reflecting presence or absence of motifs provided a highly significant independent prognostic predictor. Results were consistent between cohorts. The nonsite-specific occurrence of the detected patterns suggests that CARMA captures underlying replication and repair defects and could have a future potential in treatment stratification.

Pan-cancer analysis of whole genomes

Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale1-3. Here we report the integrative analysis of 2,658 whole-cancer genomes and their matching normal tissues across 38 tumour types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). We describe the generation of the PCAWG resource, facilitated by international data sharing using compute clouds. On average, cancer genomes contained 4-5 driver mutations when combining coding and non-coding genomic elements; however, in around 5% of cases no drivers were identified, suggesting that cancer driver discovery is not yet complete. Chromothripsis, in which many clustered structural variants arise in a single catastrophic event, is frequently an early event in tumour evolution; in acral melanoma, for example, these events precede most somatic point mutations and affect several cancer-associated genes simultaneously. Cancers with abnormal telomere maintenance often originate from tissues with low replicative activity and show several mechanisms of preventing telomere attrition to critical levels. Common and rare germline variants affect patterns of somatic mutation, including point mutations, structural variants and somatic retrotransposition. A collection of papers from the PCAWG Consortium describes non-coding mutations that drive cancer beyond those in the TERT promoter4; identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation5,6; analyses timings and patterns of tumour evolution7; describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity8,9; and evaluates a range of more-specialized features of cancer genomes8,10-18.

Recommended Guidelines for Validation, Quality Control, and Reporting of TP53 Variants in Clinical Practice

Accurate assessment of TP53 gene status in sporadic tumors and in the germline of individuals at high risk of cancer due to Li-Fraumeni Syndrome (LFS) has important clinical implications for diagnosis, surveillance, and therapy. Genomic data from more than 20,000 cancer genomes provide a wealth of information on cancer gene alterations and have confirmed TP53 as the most commonly mutated gene in human cancer. Analysis of a database of 70,000 TP53 variants reveals that the two newly discovered exons of the gene, exons 9β and 9γ, generated by alternative splicing, are the targets of inactivating mutation events in breast, liver, and head and neck tumors. Furthermore, germline rearrange-ments in intron 1 of TP53 are associated with LFS and are frequently observed in sporadic osteosarcoma. In this context of constantly growing genomic data, we discuss how screening strategies must be improved when assessing TP53 status in clinical samples. Finally, we discuss how TP53 alterations should be described by using accurate nomenclature to avoid confusion in scientific and clinical reports. Cancer Res; 77(6); 1250-60. ©2017 AACR.

The Longitudinal Transcriptional Response to Neoadjuvant Chemotherapy with and without Bevacizumab in Breast Cancer

Purpose: Chemotherapy-induced alterations to gene expression are due to transcriptional reprogramming of tumor cells or subclonal adaptations to treatment. The effect on whole-transcriptome mRNA expression was investigated in a randomized phase II clinical trial to assess the effect of neoadjuvant chemotherapy with the addition of bevacizumab.Experimental Design: Tumor biopsies and whole-transcriptome mRNA profiles were obtained at three fixed time points with 66 patients in each arm. Altogether, 358 specimens from 132 patients were available, representing the transcriptional state before treatment start, at 12 weeks and after treatment (25 weeks). Pathologic complete response (pCR) in breast and axillary nodes was the primary endpoint.Results: pCR was observed in 15 patients (23%) receiving bevacizumab and chemotherapy and 8 patients (12%) receiving only chemotherapy. In the estrogen receptor-positive patients, 11 of 54 (20%) treated with bevacizumab and chemotherapy achieved pCR, while only 3 of 57 (5%) treated with chemotherapy reached pCR. In patients with estrogen receptor-positive tumors treated with combination therapy, an elevated immune activity was associated with good response. Proliferation was reduced after treatment in both treatment arms and most pronounced in the combination therapy arm, where the reduction in proliferation accelerated during treatment. Transcriptional alterations during therapy were subtype specific, and the effect of adding bevacizumab was most evident for luminal-B tumors.Conclusions: Clinical response and gene expression response differed between patients receiving combination therapy and chemotherapy alone. The results may guide identification of patients likely to benefit from antiangiogenic therapy. Clin Cancer Res; 23(16); 4662-70. ©2017 AACR.

TP53 Mutations in Breast and Ovarian Cancer

Breast and ovarian cancers are the second and fifth leading causes of cancer deaths among women. Both breast and ovarian cancers are highly heterogeneous and are presented with diverse morphology, natural history, and response to therapy. In recent years, international efforts have led to extensive molecular characterization of both breast and ovarian tumors and identified biologically and clinically relevant subtypes of the diseases based on these molecular features. The role of TP53 in tumor initiation and progression is context dependent, and abrogation of the TP53 pathway seems to be essential for the development of basal-like breast cancers and high-grade serous ovarian cancers. These subtypes of breast and ovarian cancer show several genomic similarities including high frequency of TP53 mutation, which seems to be an early, initiating, and driving alteration in these cancer subtypes.

Breast cancer genome and transcriptome integration implicates specific mutational signatures with immune cell infiltration

A recent comprehensive whole genome analysis of a large breast cancer cohort was used to link known and novel drivers and substitution signatures to the transcriptome of 266 cases. Here, we validate that subtype-specific aberrations show concordant expression changes for, for example, TP53, PIK3CA, PTEN, CCND1 and CDH1. We find that CCND3 expression levels do not correlate with amplification, while increased GATA3 expression in mutant GATA3 cancers suggests GATA3 is an oncogene. In luminal cases the total number of substitutions, irrespective of type, associates with cell cycle gene expression and adverse outcome, whereas the number of mutations of signatures 3 and 13 associates with immune-response specific gene expression, increased numbers of tumour-infiltrating lymphocytes and better outcome. Thus, while earlier reports imply that the sheer number of somatic aberrations could trigger an immune-response, our data suggests that substitutions of a particular type are more effective in doing so than others.

Direct Transcriptional Consequences of Somatic Mutation in Breast Cancer

Disordered transcriptomes of cancer encompass direct effects of somatic mutation on transcription, coordinated secondary pathway alterations, and increased transcriptional noise. To catalog the rules governing how somatic mutation exerts direct transcriptional effects, we developed an exhaustive pipeline for analyzing RNA sequencing data, which we integrated with whole genomes from 23 breast cancers. Using X-inactivation analyses, we found that cancer cells are more transcriptionally active than intermixed stromal cells. This is especially true in estrogen receptor (ER)-negative tumors. Overall, 59% of substitutions were expressed. Nonsense mutations showed lower expression levels than expected, with patterns characteristic of nonsense-mediated decay. 14% of 4,234 rearrangements caused transcriptional abnormalities, including exon skips, exon reusage, fusions, and premature polyadenylation. We found productive, stable transcription from sense-to-antisense gene fusions and gene-to-intergenic rearrangements, suggesting that these mutation classes drive more transcriptional disruption than previously suspected. Systematic integration of transcriptome with genome data reveals the rules by which transcriptional machinery interprets somatic mutation.

•

Greater transcriptional activity in cancer than stromal cells, particularly when ER-ve

•

Intron mutations only infrequently affect splicing, even at essential splice sites

•

Distinctive RNA effects of sense-to-antisense and gene-to-intergenic rearrangements

•

Exhaustive pipeline for identifying aberrant transcripts from RNA-sequencing data

Landscape of somatic mutations in 560 breast cancer whole-genome sequences

We analysed whole-genome sequences of 560 breast cancers to advance understanding of the driver mutations conferring clonal advantage and the mutational processes generating somatic mutations. We found that 93 protein-coding cancer genes carried probable driver mutations. Some non-coding regions exhibited high mutation frequencies, but most have distinctive structural features probably causing elevated mutation rates and do not contain driver mutations. Mutational signature analysis was extended to genome rearrangements and revealed twelve base substitution and six rearrangement signatures. Three rearrangement signatures, characterized by tandem duplications or deletions, appear associated with defective homologous-recombination-based DNA repair: one with deficient BRCA1 function, another with deficient BRCA1 or BRCA2 function, the cause of the third is unknown. This analysis of all classes of somatic mutation across exons, introns and intergenic regions highlights the repertoire of cancer genes and mutational processes operating, and progresses towards a comprehensive account of the somatic genetic basis of breast cancer.

Lymphocyte Invasion in IC10/Basal-Like Breast Tumors Is Associated with Wild-Type TP53

Lymphocytic infiltration is associated with better prognosis in several epithelial malignancies including breast cancer. The tumor suppressor TP53 is mutated in approximately 30% of breast adenocarcinomas, with varying frequency across molecular subtypes. In this study of 1,420 breast tumors, we tested for interaction between TP53 mutation status and tumor subtype determined by PAM50 and integrative cluster analysis. In integrative cluster 10 (IC10)/basal-like breast cancer, we identify an association between lymphocytic infiltration, determined by an expression score, and retention of wild-type TP53. The expression-derived score agreed with the degree of lymphocytic infiltration assessed by pathologic review, and application of the Nanodissect algorithm was suggestive of this infiltration being primarily of cytotoxic T lymphocytes (CTL). Elevated expression of this CTL signature was associated with longer survival in IC10/Basal-like tumors. These findings identify a new link between the TP53 pathway and the adaptive immune response in estrogen receptor (ER)-negative breast tumors, suggesting a connection between TP53 inactivation and failure of tumor immunosurveillance.

IMPLICATIONS

The association of lymphocytic invasion of ER-negative breast tumors with the retention of wild-type TP53 implies a novel protective connection between TP53 function and tumor immunosurveillance.

Abstract PR06: Lymphocyte invasion of basal breast tumors is associated with wild-type p53

The tumor suppressor p53 is mutated in approximately 30% of breast adenocarcinomas, but the clinical implications of p53 mutation depend on the subtype of breast tumor. Although p53 mutation is associated with worse outcome for patients with Luminal B and Her2 subtypes, the same is not true in Basal tumors, despite the fact that Basal tumors have the highest frequency of p53 mutation of any subtype. Lymphocytic infiltration is associated with better prognosis in several epithelial malignancies including triple-negative breast tumors. In this study we identify an association in Basal breast cancer between lymphocytic infiltration and the number of copies of wild-type p53 by combining gene expression, p53 sequencing, DNA copy number analysis, and semi-quantitative immune cell histology in 1,281 tumors. We use the Nanodissect algorithm to define signatures for T cell types and show that the infiltration consists primarily of cytotoxic T lymphocytes. Elevated expression of this gene signature is associated with better outcome in Basal tumors, and differential pathway activity analysis identifies significant activation of caveolin 1 signaling to TGF-beta in p53 wild-type Basal tumors. These findings identify a new link between the p53 pathway and beneficial adaptive immune response in breast tumors not responsive to the estrogen receptor, supporting a connection between a lack of p53 function and the failure of tumor immunosurveillance.

This abstract is also presented as Poster A35.

Citation Format: David A. Quigley, Charlie Vaske, Laxmi Silwa-Pandit, Ruth Dannenfelser, Anita Langerød, Hans Kristian M. Vollan, Olga Troyanskaya, Suet-Feung Chin, Carlos Caldas, Allan Balmain, Anne-Lise Børresen-Dale, Vessela Kristensen. Lymphocyte invasion of basal breast tumors is associated with wild-type p53. [abstract]. In: Abstracts: AACR Special Conference on Cellular Heterogeneity in the Tumor Microenvironment; 2014 Feb 26-Mar 1; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2015;75(1 Suppl):Abstract nr PR06. doi:10.1158/1538-7445.CHTME14-PR06

TP53 mutation spectrum in breast cancer is subtype specific and has distinct prognostic relevance

PURPOSE

In breast cancer, the TP53 gene is frequently mutated and the mutations have been associated with poor prognosis. The prognostic impact of the different types of TP53 mutations across the different molecular subtypes is still poorly understood. Here, we characterize the spectrum and prognostic significance of TP53 mutations with respect to the PAM50 subtypes and integrative clusters (IC).

EXPERIMENTAL DESIGN

TP53 mutation status was obtained for 1,420 tumor samples from the METABRIC cohort by sequencing all coding exons using the Sanger method.

RESULTS

TP53 mutations were found in 28.3% of the tumors, conferring a worse overall and breast cancer-specific survival [HR = 2.03; 95% confidence interval (CI), 1.65-2.48, P < 0.001], and were also found to be an independent marker of poor prognosis in estrogen receptor-positive cases (HR = 1.86; 95% CI, 1.39-2.49, P < 0.001). The mutation spectrum of TP53 varied between the breast cancer subtypes, and individual alterations showed subtype-specific association. TP53 mutations were associated with increased mortality in patients with luminal B, HER2-enriched, and normal-like tumors, but not in patients with luminal A and basal-like tumors. Similar observations were made in ICs, where mutation associated with poorer outcome in IC1, IC4, and IC5. The combined effect of TP53 mutation, TP53 LOH, and MDM2 amplification on mortality was additive.

CONCLUSION

This study reveals that TP53 mutations have different clinical relevance in molecular subtypes of breast cancer, and suggests diverse roles for TP53 in the biology underlying breast cancer development.

AZGP1 and SPDEF mRNA expression differentiates breast carcinoma from ovarian serous carcinoma

The ANPEP, AZGP1, and SPDEF genes were previously found to be overexpressed in breast compared to ovarian carcinoma effusions. The present study validated this finding in a larger cohort consisting of both primary and metastatic tumors. ANPEP, AZGP1, and SPDEF mRNA expression was investigated in 83 breast carcinomas (57 primary carcinomas and 26 effusions) and 40 ovarian carcinomas (20 primary carcinomas and 20 effusions) using qPCR. ANPEP protein expression was immunohistochemically analyzed in 53 breast carcinoma effusions and patient-matched primary carcinomas (n = 25) and lymph node metastases (n = 16). mRNA and protein levels were studied for association with tumor type and anatomic site, and for clinical role in breast carcinoma. AZGP1 and SPDEF mRNA was overexpressed in breast compared to ovarian carcinoma (both p < 0.001). AZGP1 mRNA was overexpressed in primary breast carcinoma compared to effusions (p < 0.001), with opposite findings for ANPEP (p = 0.044). AZGP1 mRNA expression correlated with positive ER status (p = 0.032) and grade 1 histology (p = 0.011), whereas SPDEF mRNA levels were associated with positive ER (p = 0.002) and PR (p = 0.013) status and tamoxifen treatment (p = 0.004). ANPEP protein expression was higher in breast carcinoma effusions compared to primary tumors and lymph node metastases (both p = 0.001). ANPEP, AZGP1, and SPDEF levels were unrelated to disease-free or overall survival. This is the first study documenting ANPEP, AZGP1, and SPDEF expression in breast carcinoma effusions. AZGP1 and SPDEF may be novel molecular markers for the differentiation of breast from ovarian carcinoma. ANPEP may be involved in breast carcinoma progression in view of its overexpression in effusions compared to solid specimens.

Mutational processes molding the genomes of 21 breast cancers

All cancers carry somatic mutations. The patterns of mutation in cancer genomes reflect the DNA damage and repair processes to which cancer cells and their precursors have been exposed. To explore these mechanisms further, we generated catalogs of somatic mutation from 21 breast cancers and applied mathematical methods to extract mutational signatures of the underlying processes. Multiple distinct single- and double-nucleotide substitution signatures were discernible. Cancers with BRCA1 or BRCA2 mutations exhibited a characteristic combination of substitution mutation signatures and a distinctive profile of deletions. Complex relationships between somatic mutation prevalence and transcription were detected. A remarkable phenomenon of localized hypermutation, termed “kataegis,” was observed. Regions of kataegis differed between cancers but usually colocalized with somatic rearrangements. Base substitutions in these regions were almost exclusively of cytosine at TpC dinucleotides. The mechanisms underlying most of these mutational signatures are unknown. However, a role for the APOBEC family of cytidine deaminases is proposed.

PaperClip

The life history of 21 breast cancers

Cancer evolves dynamically as clonal expansions supersede one another driven by shifting selective pressures, mutational processes, and disrupted cancer genes. These processes mark the genome, such that a cancer's life history is encrypted in the somatic mutations present. We developed algorithms to decipher this narrative and applied them to 21 breast cancers. Mutational processes evolve across a cancer's lifespan, with many emerging late but contributing extensive genetic variation. Subclonal diversification is prominent, and most mutations are found in just a fraction of tumor cells. Every tumor has a dominant subclonal lineage, representing more than 50% of tumor cells. Minimal expansion of these subclones occurs until many hundreds to thousands of mutations have accumulated, implying the existence of long-lived, quiescent cell lineages capable of substantial proliferation upon acquisition of enabling genomic changes. Expansion of the dominant subclone to an appreciable mass may therefore represent the final rate-limiting step in a breast cancer's development, triggering diagnosis.

PaperClip

The landscape of cancer genes and mutational processes in breast cancer

All cancers carry somatic mutations in their genomes. A subset, known as driver mutations, confer clonal selective advantage on cancer cells and are causally implicated in oncogenesis, and the remainder are passenger mutations. The driver mutations and mutational processes operative in breast cancer have not yet been comprehensively explored. Here we examine the genomes of 100 tumours for somatic copy number changes and mutations in the coding exons of protein-coding genes. The number of somatic mutations varied markedly between individual tumours. We found strong correlations between mutation number, age at which cancer was diagnosed and cancer histological grade, and observed multiple mutational signatures, including one present in about ten per cent of tumours characterized by numerous mutations of cytosine at TpC dinucleotides. Driver mutations were identified in several new cancer genes including AKT2, ARID1B, CASP8, CDKN1B, MAP3K1, MAP3K13, NCOR1, SMARCD1 and TBX3. Among the 100 tumours, we found driver mutations in at least 40 cancer genes and 73 different combinations of mutated cancer genes. The results highlight the substantial genetic diversity underlying this common disease.

The genomic and transcriptomic architecture of 2,000 breast tumours reveals novel subgroups

The elucidation of breast cancer subgroups and their molecular drivers requires integrated views of the genome and transcriptome from representative numbers of patients. We present an integrated analysis of copy number and gene expression in a discovery and validation set of 997 and 995 primary breast tumours, respectively, with long-term clinical follow-up. Inherited variants (copy number variants and single nucleotide polymorphisms) and acquired somatic copy number aberrations (CNAs) were associated with expression in ~40% of genes, with the landscape dominated by cis- and trans-acting CNAs. By delineating expression outlier genes driven in cis by CNAs, we identified putative cancer genes, including deletions in PPP2R2A, MTAP and MAP2K4. Unsupervised analysis of paired DNA–RNA profiles revealed novel subgroups with distinct clinical outcomes, which reproduced in the validation cohort. These include a high-risk, oestrogen-receptor-positive 11q13/14 cis-acting subgroup and a favourable prognosis subgroup devoid of CNAs. Trans-acting aberration hotspots were found to modulate subgroup-specific gene networks, including a TCR deletion-mediated adaptive immune response in the ‘CNA-devoid’ subgroup and a basal-specific chromosome 5 deletion-associated mitotic network. Our results provide a novel molecular stratification of the breast cancer population, derived from the impact of somatic CNAs on the transcriptome.

Functional characterization of the 19q12 amplicon in grade III breast cancers

INTRODUCTION

The 19q12 locus is amplified in a subgroup of oestrogen receptor (ER)-negative grade III breast cancers. This amplicon comprises nine genes, including cyclin E1 (CCNE1), which has been proposed as its 'driver'. The aim of this study was to identify the genes within the 19q12 amplicon whose expression is required for the survival of cancer cells harbouring their amplification.

METHODS

We investigated the presence of 19q12 amplification in a series of 313 frozen primary breast cancers and 56 breast cancer cell lines using microarray comparative genomic hybridisation (aCGH). The nine genes mapping to the smallest region of amplification on 19q12 were silenced using RNA interference in phenotypically matched breast cancer cell lines with (MDA-MB-157 and HCC1569) and without (Hs578T, MCF7, MDA-MB-231, ZR75.1, JIMT1 and BT474) amplification of this locus. Genes whose silencing was selectively lethal in amplified cells were taken forward for further validation. The effects of cyclin-dependent kinase 2 (CDK2) silencing and chemical inhibition were tested in cancer cells with and without CCNE1 amplification.

RESULTS

19q12 amplification was identified in 7.8% of ER-negative grade III breast cancer. Of the nine genes mapping to this amplicon, UQCRFS1, POP4, PLEKHF1, C19ORF12, CCNE1 and C19ORF2 were significantly over-expressed when amplified in primary breast cancers and/or breast cancer cell lines. Silencing of POP4, PLEKHF1, CCNE1 and TSZH3 selectively reduced cell viability in cancer cells harbouring their amplification. Cancer cells with CCNE1 amplification were shown to be dependent on CDK2 expression and kinase activity for their survival.

CONCLUSIONS

The 19q12 amplicon may harbour more than a single 'driver', given that expression of POP4, PLEKHF1, CCNE1 and TSZH3 is required for the survival of cancer cells displaying their amplification. The observation that cancer cells harbouring CCNE1 gene amplification are sensitive to CDK2 inhibitors provides a rationale for the testing of these chemical inhibitors in a subgroup of patients with ER-negative grade III breast cancers.

Mutant p53 disrupts mammary tissue architecture via the mevalonate pathway

p53 is a frequent target for mutation in human tumors, and mutant p53 proteins can actively contribute to tumorigenesis. We employed a three-dimensional culture model in which nonmalignant breast epithelial cells form spheroids reminiscent of acinar structures found in vivo, whereas breast cancer cells display highly disorganized morphology. We found that mutant p53 depletion is sufficient to phenotypically revert breast cancer cells to a more acinar-like morphology. Genome-wide expression analysis identified the mevalonate pathway as significantly upregulated by mutant p53. Statins and sterol biosynthesis intermediates reveal that this pathway is both necessary and sufficient for the phenotypic effects of mutant p53 on breast tissue architecture. Mutant p53 associates with sterol gene promoters at least partly via SREBP transcription factors. Finally, p53 mutation correlates with highly expressed sterol biosynthesis genes in human breast tumors. These findings implicate the mevalonate pathway as a therapeutic target for tumors bearing mutations in p53.

Association of the germline TP53 R72P and MDM2 SNP309 variants with breast cancer survival in specific breast tumor subgroups

The tumor suppressor gene TP53 and its regulator MDM2 are both important players in the DNA-damage repair "TP53 response pathway". Common germline polymorphisms in these genes may affect outcome in patients with tumors characterized by additional somatic changes in the same or a related pathway. To evaluate this hypothesis, we determined the effect of the common germline TP53 R72P and MDM2 SNP309 polymorphisms on breast cancer survival in a consecutive cohort of breast cancer patients (age at diagnosis <53 years, n = 295) with gene expression data available. Patients were classified in subgroups according to their tumor TP53 mutation status and three gene expression profiles; a TP53 mutation status expression signature, a PTEN/PI3K pathway signature and the 70-gene prognosis profile. Survival analyses were performed using Cox regression models adjusting for clinico-pathological characteristics and treatment. An increase in breast cancer-specific mortality was observed for carriers of the germline MDM2 SNP309 rare GG-genotype (range hazard ratios: 2-3) or TP53 R72P heterozygous GC-genotype (range hazard ratios: 1-2) compared to those having the common genotypes within subgroups of tumors displaying a "more aggressive phenotype" gene expression profile. There was no evidence of such an effect on survival within the TP53-mutated tumor group for TP53 R72P carriers but a suggestion of an effect for MDM2 SNP309 carriers (GG vs. TT-genotype HR 2.99, P = 0.06). These results indicate that common polymorphisms in specific pathways may add to the worse prognosis of patients with tumors in which these pathways are affected by somatic alterations.

Genomic architecture characterizes tumor progression paths and fate in breast cancer patients

Distinct molecular subtypes of breast carcinomas have been identified, but translation into clinical use has been limited. We have developed two platform-independent algorithms to explore genomic architectural distortion using array comparative genomic hybridization data to measure (i) whole-arm gains and losses [whole-arm aberration index (WAAI)] and (ii) complex rearrangements [complex arm aberration index (CAAI)]. By applying CAAI and WAAI to data from 595 breast cancer patients, we were able to separate the cases into eight subgroups with different distributions of genomic distortion. Within each subgroup data from expression analyses, sequencing and ploidy indicated that progression occurs along separate paths into more complex genotypes. Histological grade had prognostic impact only in the luminal-related groups, whereas the complexity identified by CAAI had an overall independent prognostic power. This study emphasizes the relation among structural genomic alterations, molecular subtype, and clinical behavior and shows that objective score of genomic complexity (CAAI) is an independent prognostic marker in breast cancer.

Combining gene signatures improves prediction of breast cancer survival

Background

Several gene sets for prediction of breast cancer survival have been derived from whole-genome mRNA expression profiles. Here, we develop a statistical framework to explore whether combination of the information from such sets may improve prediction of recurrence and breast cancer specific death in early-stage breast cancers. Microarray data from two clinically similar cohorts of breast cancer patients are used as training (n = 123) and test set (n = 81), respectively. Gene sets from eleven previously published gene signatures are included in the study.

Principal Findings

To investigate the relationship between breast cancer survival and gene expression on a particular gene set, a Cox proportional hazards model is applied using partial likelihood regression with an L2 penalty to avoid overfitting and using cross-validation to determine the penalty weight. The fitted models are applied to an independent test set to obtain a predicted risk for each individual and each gene set. Hierarchical clustering of the test individuals on the basis of the vector of predicted risks results in two clusters with distinct clinical characteristics in terms of the distribution of molecular subtypes, ER, PR status, TP53 mutation status and histological grade category, and associated with significantly different survival probabilities (recurrence: p = 0.005; breast cancer death: p = 0.014). Finally, principal components analysis of the gene signatures is used to derive combined predictors used to fit a new Cox model. This model classifies test individuals into two risk groups with distinct survival characteristics (recurrence: p = 0.003; breast cancer death: p = 0.001). The latter classifier outperforms all the individual gene signatures, as well as Cox models based on traditional clinical parameters and the Adjuvant! Online for survival prediction.

Conclusion

Combining the predictive strength of multiple gene signatures improves prediction of breast cancer survival. The presented methodology is broadly applicable to breast cancer risk assessment using any new identified gene set.

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

BACKGROUND

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 BRCA1-mutation carriers. The genomic instability of BRCA1-mutated tumors can be effectively targeted with DNA-damaging agents and poly-(ADP-ribose) polymerase 1 (PARP1) inhibitors. Molecular similarities between BLBCs and BRCA1-mutated tumors may therefore provide predictive markers for therapeutic response of BLBCs.

METHODS

There are several known molecular features characteristic for BRCA1-mutated breast tumors: 1) increased numbers of genomic aberrations, 2) a distinct pattern of genomic aberrations, 3) a high frequency of TP53 mutations and 4) a high incidence of complex, protein-truncating TP53 mutations. We compared the frequency of TP53 mutations and the pattern and amount of genomic aberrations between BRCA1-mutated breast tumors, BLBCs and luminal breast tumors by TP53 gene sequencing and array-based comparative genomics hybridization (aCGH) analysis.

RESULTS

We found that the high incidence of protein truncating TP53 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.

CONCLUSIONS

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.

p53-Dependent repression of focal adhesion kinase in response to estradiol in breast cancer cell-lines

Mutations in the TP53 suppressor gene are frequent in breast cancers. These mutations are associated with poor prognosis, thought to be due to proliferative advantage and poor response to chemotherapy associated with loss of p53 function. The focal adhesion kinase (FAK/PTK2), a tyrosine kinase, is over-expressed in a variety of human tumors including breast cancers. FAK is a critical regulator of adhesion and motility and its over-expression is associated with increased metastatic potential. Recently, FAK promoter has been shown to contain p53 responsive elements and to be down-regulated by DNA-damage in a p53-dependent manner. Here, we have used five estrogen-dependent breast cancer cells lines with different p53 status, including an isogenic model, to show that FAK expression was regulated in a p53-dependent manner in response to estradiol. FAK protein and mRNA expression were down-regulated by estradiol in wild-type but not mutant p53 cells. Moreover, silencing wild-type p53 increased FAK expression, while over expressing p53 repressed FAK expression. ChIP experiment showed that p53 bound to FAK promoter in the presence of estradiol in p53 wild-type but not in mutant p53 cells, suggesting a direct role of p53 in down regulating FAK mRNA expression. FAK mRNA expression was also found to correlate with TP53 mutation status in a series of breast tumors. Finally, loss of FAK down-regulation in p53 mutant cells was correlated with increased proliferation and invasion upon estradiol stimulation, while FAK silencing reduced invasion. These results suggest that p53 is an important down regulator of FAK and that loss of p53 function in breast cancer may contribute to the metastatic potential of estrogen-responsive tumors through uncontrolled FAK expression upon estrogens stimulation.

Epithelial-stromal interaction 1 (EPSTI1) substitutes for peritumoral fibroblasts in the tumor microenvironment

Tumor cells can activate stroma, yet the implication of this activation in terms of reciprocal induction of gene expression in tumor cells is poorly understood. Epithelial Stromal Interaction 1 (EPSTI1) is an interferon response gene originally isolated from heterotypic recombinant cultures of human breast cancer cells and activated breast myofibroblasts. Here we describe the first immunolocalization of EPSTI1 in normal and cancerous breast tissue, and we provide evidence for a role of this molecule in the regulation of tumor cell properties and epithelial-mesenchymal transition. In general, no EPSTI1 staining was observed in normal breast epithelial cells from reduction mammoplasties (n=25). However, in carcinomas, staining was positive in 22 of 40 biopsies and inversely correlated with the level of differentiation. To address the function of EPSTI1, we expressed EPSTI1 ectopically in one cell line and silenced endogenous EPSTI1 by RNA interference in another. Irrespective of the experimental approach, EPSTI1 expression led to an increase in tumorsphere formation-a property associated with breast stem/progenitor cells. Most remarkably, we show that EPSTI1, by conveying spread of tumor cells, can replace peritumoral activated fibroblasts in a tumor environment assay. These observations implicate EPSTI1 as a hitherto unappreciated regulator of tumor cell properties.

Focal amplification and oncogene dependency of GAB2 in breast cancer

DNA amplifications in breast cancer are frequent on chromosome 11q, in which multiple driver oncogenes likely reside in addition to cyclin D1 (CCND1). One such candidate, the scaffolding adapter protein, GRB2-associated binding protein 2 (GAB2), functions in ErbB signaling and was recently shown to enhance mammary epithelial cell proliferation, and metastasis of ERBB2 (HER2/neu)-driven murine breast cancer. However, the amplification status and function of GAB2 in the context of amplification remain undefined. In this study, by genomic profiling of 172 breast tumors, and fluorescence in situ hybridization validation in an independent set of 210 scorable cases, we observed focal amplification spanning GAB2 (11q14.1) independent of CCND1 (11q13.2) amplification, consistent with a driver role. Further, small interfering RNA (siRNA)-mediated knockdown of GAB2 in breast cancer lines (SUM52, SUM44PE and MDA468) with GAB2 amplification revealed a dependency on GAB2 for cell proliferation, cell-cycle progression, survival and invasion, likely mediated through altered phosphatidylinositol 3-kinase (PI3K) and mitogen-activated protein kinase (MAPK) signaling. GAB2 knockdown also reduced proliferation and survival in a cell line (BT474) with ERBB2 amplification, consistent with the possibility that GAB2 can function downstream of ERBB2. Our studies implicate focal amplification of GAB2 in breast carcinogenesis, and underscore an oncogenic role of scaffolding adapter proteins, and a potential new point of therapeutic intervention.

Somatic mutation databases as tools for molecular epidemiology and molecular pathology of cancer: proposed guidelines for improving data collection, distribution, and integration

There are currently less than 40 locus-specific databases (LSDBs) and one large general database that curate data on somatic mutations in human cancer genes. These databases have different scope and use different annotation standards and database systems, resulting in duplicated efforts in data curation, and making it difficult for users to find clear and consistent information. As data related to somatic mutations are generated at an increasing pace it is urgent to create a framework for improving the collecting of this information and making it more accessible to clinicians, scientists, and epidemiologists to facilitate research on biomarkers. Here we propose a data flow for improving the connectivity between existing databases and we provide practical guidelines for data reporting, database contents, and annotation standards. These proposals are based on common standards recommended by the Human Genome Variation Society (HGVS) with additions related to specific requirements of somatic mutations in cancer. Indeed, somatic mutations may be used in molecular pathology and clinical studies to characterize tumor types, help treatment choice, predict response to treatment and patient outcome, or in epidemiological studies as markers for tumor etiology or exposure assessment. Thus, specific annotations are required to cover these diverse research topics. This initiative is meant to promote collaboration and discussion on these issues and the development of adequate resources that would avoid the loss of extremely valuable information generated by years of basic and clinical research.

Full sequencing of TP53 identifies identical mutations within in situ and invasive components in breast cancer suggesting clonal evolution

In breast cancer, previous studies have suggested that somatic TP53 mutations are likely to be an early event. However, there are controversies regarding the cellular origin and linear course of breast cancer. The purpose of this study was to investigate tumor evolution in breast cancer by analyzing TP53 mutation status in tumors from various stages of the disease. The entire coding sequence of TP53 was sequenced in a cohort of pure ductal carcinoma in situ (DCIS), pure invasive cancer (≤15mm) and mixed lesions (i.e. invasive cancer with an in situ component). Of 118 tumor samples, 19 were found to harbor a TP53 mutation; 5 (15.6%) of the pure DCIS, 4 (10.5%) of the pure invasive cancers and 10 (20.8%) of the mixed lesions. In the mixed lesions, both the invasive and the DCIS components showed the same mutation in all 5 cases where the two components were successfully microdissected. Presence of the same mutation in both DCIS and invasive components from the same tumor indicates same cellular origin. The role of mutant TP53 in the progression of breast cancer is less clear and may vary between subtypes.

TP53 mutation is an early event in breast cancer progression

Abstract #1047

Background: In breast cancer, previous studies have shown that somatic TP53 mutations cause a more aggressive disease with poor clinical outcome and may impact treatment response. Although TP53 mutation is considered to be an early event in breast cancer, the timing of TP53 mutations is not known, and there are controversies regarding the cellular origin and linear model of breast cancer. The purpose of this study was to investigate whether TP53 mutations are early events in breast cancer progression.&#x2028; Methods: From a population-based cohort of women diagnosed between 1986 and 2004 either with a pure ductal carcinoma in situ (DCIS), a pure invasive cancer (&lt;15mm) or a mixed lesion (i.e. invasive cancer with a DCIS component), we included 118 women with stored frozen tissue. Mixed lesions were microdissected using LCM (laser capture microdissection) on a PALM slide to separate in situ and invasive tumor cells. DNA was isolated using phenol-chloroform extraction. The entire coding sequence of TP53 was analyzed for mutations by direct sequencing on a 3730 DNA analyzer.&#x2028; Results: Of 118 tumor samples, 19 were detected with a TP53 mutation; five 5 of 32 (15.6%) pure DCIS, 4 of 38 (10.5%) pure invasive cancers and 10 of 48 (20.8%) mixed lesions. In the mixed lesions, both the invasive and the DCIS component showed the same mutation in all 5 cases where we successfully could microdissect the two components separately. Pure DCIS demonstrated missense mutations (4/5, 80%) more frequently than pure invasive cancers (2/4, 50%) and mixed lesions (4/10, 40%), although this difference was not statistically significant (p=0.3). Also, the frequency of missense mutations in the DNA binding domain was not statistically different between the three groups.&#x2028; Conclusion: TP53 mutation is likely an early event in breast cancer, occurring previous to or in the in situ stage. Presence of the same mutation in both DCIS and invasive component from the same tumor indicates same cellular origin. The role of mutant TP53 in the progression into invasive cancer is less clear and may vary between subtypes of breast cancer.

Citation Information: Cancer Res 2009;69(2 Suppl):Abstract nr 1047.

EZH2 and BMI1 inversely correlate with prognosis and TP53 mutation in breast cancer

Introduction

PolycombGroup (PcG) proteins maintain gene repression through histone modifications and have been implicated in stem cell regulation and cancer. EZH2 is part of Polycomb Repressive Complex 2 (PRC2) and trimethylates H3K27. This histone mark recruits the BMI1-containing PRC1 that silences the genes marked by PRC2. Based on their role in stem cells, EZH2 and BMI1 have been predicted to contribute to a poor outcome for cancer patients.

Methods

We have analysed the expression of EZH2 and BMI1 in a well-characterised dataset of 295 human breast cancer samples.

Results

Interestingly, although EZH2 overexpression correlates with a poor prognosis in breast cancer, BMI1 overexpression correlates with a good outcome. Although this may reflect transformation of different cell types, we also observed a functional difference. The PcG-target genes INK4A and ARF are not expressed in tumours with high BMI1, but they are expressed in tumours with EZH2 overexpression. ARF expression results in tumour protein P53 (TP53) activation, and we found a significantly higher proportion of TP53 mutations in tumours with high EZH2. This may explain why tumours with high EZH2 respond poorly to therapy, in contrast to tumours with high BMI1.

Conclusions

Overall, our data highlight that whereas EZH2 and BMI1 may function in a 'linear' pathway in normal development, their overexpression has different functional consequences for breast tumourigenesis.

Human papillomavirus DNA and e6/e7 mRNA status in relation to survival of patients treated for cervical squamous cell carcinoma

The prognostic significance of human papillomavirus (HPV) DNA and E6/E7 mRNA, the presence of specific types, and the physical state of HPV DNA, were studied in 202 cervical squamous cell carcinomas. Absence or non-detectable levels of high-risk (types 16, 18, 31, 33, 35, 45, 52 and 58) E6/E7 mRNA, using the real-time nucleic acid sequence based amplification (NASBA) assay, and absence of HPV high-risk/HPV 6, 26, 66, 69, 73 (all methods collectively) were associated with poor overall survival in univariate analysis (P = 0.04 and P = 0.03, respectively) and had independent prognostic value in multivariate analysis (P = 0.01 and P = 0.03, respectively) including FIGO stage and age. Based on the individual results of type-specific PCR and in situ hybridization (ISH), the presence of HPV DNA was not found to be a prognostic factor. Likewise, concerning the presence of specific HPV types and the HPV integration status (determined by ISH), no prognostic significance was found. Mutation analyses of the TP53 gene revealed mutations in 3 of the 6 HPV negative samples (50%).

CAMK1D amplification implicated in epithelial-mesenchymal transition in basal-like breast cancer

Breast cancer exhibits clinical and molecular heterogeneity, where expression profiling studies have identified five major molecular subtypes. The basal-like subtype, expressing basal epithelial markers and negative for estrogen receptor (ER), progesterone receptor (PR) and HER2, is associated with higher overall levels of DNA copy number alteration (CNA), specific CNAs (like gain on chromosome 10p), and poor prognosis. Discovering the molecular genetic basis of tumor subtypes may provide new opportunities for therapy. To identify the driver oncogene on 10p associated with basal-like tumors, we analyzed genomic profiles of 172 breast carcinomas. The smallest shared region of gain spanned just seven genes at 10p13, including calcium/calmodulin-dependent protein kinase ID (CAMK1D), functioning in intracellular signaling but not previously linked to cancer. By microarray, CAMK1D was overexpressed when amplified, and by immunohistochemistry exhibited elevated expression in invasive carcinomas compared to carcinoma in situ. Engineered overexpression of CAMK1D in non-tumorigenic breast epithelial cells led to increased cell proliferation, and molecular and phenotypic alterations indicative of epithelial-mesenchymal transition (EMT), including loss of cell-cell adhesions and increased cell migration and invasion. Our findings identify CAMK1D as a novel amplified oncogene linked to EMT in breast cancer, and as a potential therapeutic target with particular relevance to clinically unfavorable basal-like tumors.

Common genetic variation in GATA-binding protein 3 and differential susceptibility to breast cancer by estrogen receptor alpha tumor status

GATA-binding protein 3 (GATA3) is a transcription factor and a putative tumor suppressor that is highly expressed in normal breast luminal epithelium and estrogen receptor alpha (ER)-positive breast tumors. We hypothesized that common genetic variation in GATA3 could influence breast carcinogenesis. Four tag single-nucleotide polymorphisms (SNP) in GATA3 and its 3' flanking gene FLJ4598 were genotyped in two case control studies in Norway and Poland (2,726 cases and 3,420 controls). Analyses of pooled data suggested a reduced risk of breast cancer associated with two intronic variants in GATA3 in linkage disequilibrium (rs3802604 in intron 3 and rs570613 in intron 4). Odds ratio (95% confidence interval) for rs570613 heterozygous and rare homozygous versus common homozygous were 0.85 (0.75-1.95) and 0.82 (0.62-0.96), respectively (P(trend)=0.004). Stronger associations were observed for subjects with ER-negative, than ER-positive, tumors (P(heterogeneity)=0.01 for rs3802604; P(heterogeneity)=0.09 for rs570613). Although no individual SNPs were associated with ER-positive tumors, two haplotypes (GGTC in 2% of controls and AATT in 7% of controls) showed significant and consistent associations with increased risk for these tumors when compared with the common haplotype (GATT in 46% of controls): 1.71 (1.27-2.32) and 1.26 (1.03-1.54), respectively. In summary, data from two independent study populations showed two intronic variants in GATA3 associated with overall decreases in breast cancer risk and suggested heterogeneity of these associations by ER status. These differential associations are consistent with markedly different levels of GATA3 protein by ER status. Additional epidemiologic studies are needed to clarify these intriguing relationships.

Common genetic variation in TP53 and its flanking genes, WDR79 and ATP1B2, and susceptibility to breast cancer

Germline mutations in the tumor suppressor gene TP53 are associated with high incidence of early-onset malignancies, and somatic mutations occur in 20-40% of all breast cancer cases. We investigated the association of common genetic variation in TP53 and its flanking genes, WDR79 and ATP1B2, with risk for breast cancer. Single nucleotide polymorphisms (SNPs) identified in a re-sequence analysis were genotyped in 2 large case-control studies including 731 cases and 1,124 controls from Norway, and 1,995 cases and 2,296 controls from Poland. Analyses of the pooled data showed no SNPs in TP53 to be significantly associated with risk for breast cancer. However, we found a significant and consistent association with risk for a SNP in exon 1 (R68G) of the 5' neighboring gene WDR79 (rs2287499, OR (95% CI) = 1.08 (0.95-1.23) for CG vs. CC and 1.60 (1.04-2.47) for GG vs. CC, p-trend = 0.01). Stratification by ER and PR status, showed these increases in risk to be limited to ER negative tumors (OR (95% CI) per variant allele: 1.42 (1.18-1.71) p-trend = 0.00009). In addition, 2 TP53 SNPs (rs17887200 3'of STP and rs12951053 in intron 7) showing weak and non-significant overall increases in risk, were also associated with ER negative tumors (1.48 (1.11-1.93) p-trend = 0.01 and 1.29 (1.06-1.58) p-trend = 0.009, respectively). In conclusion, this comprehensive evaluation of common genetic variation in TP53 and its flanking genes found no significant overall associations between SNPs in TP53 and breast cancer risk. However, data suggested that common variation in TP53 or WDR79 could be associated with ER negative breast cancers.

TP53 mutation status and gene expression profiles are powerful prognostic markers of breast cancer

Introduction

Gene expression profiling of breast carcinomas has increased our understanding of the heterogeneous biology of this disease and promises to impact clinical care. The aim of this study was to evaluate the prognostic value of gene expression-based classification along with established prognostic markers and mutation status of the TP53 gene (tumour protein p53) in a group of breast cancer patients with long-term (12 to 16 years) follow-up.

Methods

The clinical and histopathological parameters of 200 breast cancer patients were studied for their effects on clinical outcome using univariate/multivariate Cox regression. The prognostic impact of mutations in the TP53 gene, identified using temporal temperature gradient gel electrophoresis and sequencing, was also evaluated. Eighty of the samples were analyzed for gene expression using 42 K cDNA microarrays and the patients were assigned to five previously defined molecular expression groups. The strength of the gene expression based classification versus standard markers was evaluated by adding this variable to the Cox regression model used to analyze all samples.

Results

Both univariate and multivariate analysis showed that TP53 mutation status, tumor size and lymph node status were the strongest predictors of breast cancer survival for the whole group of patients. Analyses of the patients with gene expression data showed that TP53 mutation status, gene expression based classification, tumor size and lymph node status were significant predictors of survival. Breast cancer cases in the 'basal-like' and 'ERBB2⁺' gene expression subgroups had a very high mortality the first two years, while the 'highly proliferating luminal' cases developed the disease more slowly, showing highest mortality after 5 to 8 years. The TP53 mutation status showed strong association with the 'basal-like' and 'ERBB2⁺' subgroups, and tumors with mutation had a characteristic gene expression pattern.

Conclusion

TP53 mutation status and gene-expression based groups are important survival markers of breast cancer, and these molecular markers may provide prognostic information that complements clinical variables. The study adds experience and knowledge to an ongoing characterization and classification of the disease.

Somatic sequence alterations in twenty-one genes selected by expression profile analysis of breast carcinomas

INTRODUCTION

Genomic alterations have been observed in breast carcinomas that affect the capacity of cells to regulate proliferation, signaling, and metastasis. Re-sequence studies have investigated candidate genes based on prior genetic observations (changes in copy number or regions of genetic instability) or other laboratory observations and have defined critical somatic mutations in genes such as TP53 and PIK3CA.

METHODS

We have extended the paradigm and analyzed 21 genes primarily identified by expression profiling studies, which are useful for breast cancer subtyping and prognosis. This study conducted a bidirectional re-sequence analysis of all exons and 5', 3', and evolutionarily conserved regions (spanning more than 16 megabases) in 91 breast tumor samples.

RESULTS

Eighty-seven unique somatic alterations were identified in 16 genes. Seventy-eight were single base pair alterations, of which 23 were missense mutations; 55 were distributed across conserved intronic regions or the 5' and 3' regions. There were nine insertion/deletions. Because there is no a priori way to predict whether any one of the identified synonymous and noncoding somatic alterations disrupt function, analysis unique to each gene will be required to establish whether it is a tumor suppressor gene or whether there is no effect. In five genes, no somatic alterations were observed.

CONCLUSION

The study confirms the value of re-sequence analysis in cancer gene discovery and underscores the importance of characterizing somatic alterations across genes that are related not only by function, or functional pathways, but also based upon expression patterns.

Expression of full-length p53 and its isoform Deltap53 in breast carcinomas in relation to mutation status and clinical parameters

Background

The tumor suppressor gene p53 (TP53) controls numerous signaling pathways and is frequently mutated in human cancers. Novel p53 isoforms suggest alternative splicing as a regulatory feature of p53 activity.

Results

In this study we have analyzed mRNA expression of both wild-type and mutated p53 and its respective Δp53 isoform in 88 tumor samples from breast cancer in relation to clinical parameters and molecular subgroups. Three-dimensional structure differences for the novel internally deleted p53 isoform Δp53 have been predicted. We confirmed the expression of Δp53 mRNA in tumors using quantitative real-time PCR technique. The mRNA expression levels of the two isoforms were strongly correlated in both wild-type and p53-mutated tumors, with the level of the Δp53 isoform being approximately 1/3 of that of the full-length p53 mRNA. Patients expressing mutated full-length p53 and non-mutated (wild-type) Δp53, "mutational hybrids", showed a slightly higher frequency of patients with distant metastasis at time of diagnosis compared to other patients with p53 mutations, but otherwise did not differ significantly in any other clinical parameter. Interestingly, the p53 wild-type tumors showed a wide range of mRNA expression of both p53 isoforms. Tumors with mRNA expression levels in the upper or lower quartile were significantly associated with grade and molecular subtypes. In tumors with missense or in frame mutations the mRNA expression levels of both isoforms were significantly elevated, and in tumors with nonsense, frame shift or splice mutations the mRNA levels were significantly reduced compared to those expressing wild-type p53.

Conclusion

Expression of p53 is accompanied by the functionally different isoform Δp53 at the mRNA level in cell lines and human breast tumors. Investigations of "mutational hybrid" patients highlighted that wild-type Δp53 does not compensates for mutated p53, but rather may be associated with a worse prognosis. In tumors, both isoforms show strong correlations in different mutation-dependent mRNA expression patterns.

Discovery and validation of breast cancer subtypes

Background

Previous studies demonstrated breast cancer tumor tissue samples could be classified into different subtypes based upon DNA microarray profiles. The most recent study presented evidence for the existence of five different subtypes: normal breast-like, basal, luminal A, luminal B, and ERBB2⁺.

Results

Based upon the analysis of 599 microarrays (five separate cDNA microarray datasets) using a novel approach, we present evidence in support of the most consistently identifiable subtypes of breast cancer tumor tissue microarrays being: ESR1⁺/ERBB2^-, ESR1^-/ERBB2^-, and ERBB2⁺ (collectively called the ESR1/ERBB2 subtypes). We validate all three subtypes statistically and show the subtype to which a sample belongs is a significant predictor of overall survival and distant-metastasis free probability.

Conclusion

As a consequence of the statistical validation procedure we have a set of centroids which can be applied to any microarray (indexed by UniGene Cluster ID) to classify it to one of the ESR1/ERBB2 subtypes. Moreover, the method used to define the ESR1/ERBB2 subtypes is not specific to the disease. The method can be used to identify subtypes in any disease for which there are at least two independent microarray datasets of disease samples.

The clinical value of somatic TP53 gene mutations in 1,794 patients with breast cancer

To investigate the clinical value of somatic TP53 mutations in breast cancer, we assembled clinical and molecular data on 1,794 women with primary breast cancer with long-term follow-up and whose tumor has been screened for mutation in exons 5 to 8 of TP53 by gene sequencing. TP53 mutations were more frequent in tumors of ductal and medullar types, aggressive phenotype (high grade, large size, node positive cases, and low hormone receptor content) and in women <60 years old. TP53 mutations within exons 5 to 8 conferred an elevated risk of breast cancer-specific death of 2.27 (relative risk >10 years; P < 0.0001) compared with patients with no such mutation. The prognostic value of TP53 mutation was independent of tumor size, node status, and hormone receptor content, confirming and reconciling previous findings in smaller series. Moreover, an interaction between TP53 mutation and progesterone receptor (PR) status was revealed, TP53 mutation combined with the absence of progesterone receptor being associated with the worst prognosis. Whereas previous studies have emphasized the fact that missense mutations in the DNA-binding motifs have a worse prognosis than missense mutations outside these motifs, we show that non-missense mutations have prognostic value similar to missense mutations in DNA-binding motifs. Nonetheless, specific missense mutants (codon 179 and R248W) seem to be associated with an even worse prognosis. These results, obtained on the largest series analyzed thus far, show that TP53 mutations identified by gene sequencing have an independent prognostic value in breast cancer and could have potential uses in clinical practice.

Gene expression profiling of breast cancer in relation to estrogen receptor status and estrogen-metabolizing enzymes: clinical implications

Interactions between luminal epithelial cells and their surrounding microenvironment govern the normal development and function of the mammary gland. Estradiol plays a key role in abnormal intracellular signaling, which contributes to the development and progression of breast tumors. The present article summarizes the results from a microarray whole genome gene expression analysis as well as a quantitative analysis of the mRNA expression of members of the estradiol metabolic and signaling pathways in the tumors of postmenopausal breast cancer patients. The analysis of the variation in whole genome gene expression resulted in a tumor classification comprising several distinct groups with distinct expression of the estrogen receptor (ER). The parallel study on the expression of only nine mRNA transcripts of members of the estradiol pathways resulted in two main clusters, representing ER- and ER tumors. The mRNA expression of the estradiol-metabolizing enzymes did not follow the expression of the ER in all cases, leading to the recognition of several further subclasses of tumors. When the tumor classes obtained by whole genome gene expression analysis were compared with those obtained by independent quantitation of the estradiol-metabolizing enzymes, a statistically significant association between both classification groups was observed. These findings point to a possible association between development of a tumor with a particular expression profile and its capacity to synthesize estradiol as measured by the expression of the transcripts for the necessary key enzymes. Further, whole genome expression patterns were studied in 12 patients treated with anastrozole. Using significance analysis of microarrays, we identified 298 genes significantly differently expressed between partial response and progressive disease groups.