Amplified loci on chromosomes 8 and 17 predict early relapse in ER-positive breast cancers

TONSL Is an Immortalizing Oncogene and a Therapeutic Target in Breast Cancer

Study of genomic aberrations leading to immortalization of epithelial cells has been technically challenging due to the lack of isogenic models. To address this, we used healthy primary breast luminal epithelial cells of different genetic ancestry and their hTERT-immortalized counterparts to identify transcriptomic changes associated with immortalization. Elevated expression of TONSL (Tonsoku-like, DNA repair protein) was identified as one of the earliest events during immortalization. TONSL, which is located on chromosome 8q24.3, was found to be amplified in approximately 20% of breast cancers. TONSL alone immortalized primary breast epithelial cells and increased telomerase activity, but overexpression was insufficient for neoplastic transformation. However, TONSL-immortalized primary cells overexpressing defined oncogenes generated estrogen receptor-positive adenocarcinomas in mice. Analysis of a breast tumor microarray with approximately 600 tumors revealed poor overall and progression-free survival of patients with TONSL-overexpressing tumors. TONSL increased chromatin accessibility to pro-oncogenic transcription factors, including NF-κB and limited access to the tumor-suppressor p53. TONSL overexpression resulted in significant changes in the expression of genes associated with DNA repair hubs, including upregulation of several genes in the homologous recombination (HR) and Fanconi anemia pathways. Consistent with these results, TONSL-overexpressing primary cells exhibited upregulated DNA repair via HR. Moreover, TONSL was essential for growth of TONSL-amplified breast cancer cell lines in vivo, and these cells were sensitive to TONSL-FACT complex inhibitor CBL0137. Together, these findings identify TONSL as a regulator of epithelial cell immortalization to facilitate cancer initiation and as a target for breast cancer therapy.

SIGNIFICANCE

The chr.8q24.3 amplicon-resident gene TONSL is upregulated during the initial steps of tumorigenesis to support neoplastic transformation by increasing DNA repair and represents a potential therapeutic target for treating breast cancer.

Aberrant Cyclin D1 splicing in cancer: from molecular mechanism to therapeutic modulation

Cyclin D1 (CCND1), a crucial mediator of cell cycle progression, possesses many mutation types with different mutation frequencies in human cancers. The G870A mutation is the most common mutation in CCND1, which produces two isoforms: full-length CCND1a and divergent C-terminal CCND1b. The dysregulation of the CCND1 isoforms is associated with multiple human cancers. Exploring the molecular mechanism of CCND1 isoforms has offer new insight for cancer treatment. On this basis, the alterations of CCND1 gene are described, including amplification, overexpression, and mutation, especially the G870A mutation. Subsequently, we review the characteristics of CCND1 isoforms caused by G870A mutation. Additionally, we summarize cis-regulatory elements, trans-acting factors, and the splice mutation involved in splicing regulation of CCND1. Furthermore, we highlight the function of CCND1 isoforms in cell cycle, invasion, and metastasis in cancers. Importantly, the clinical role of CCND1 isoforms is also discussed, particularly concerning prognosis, chemotherapy, and radiotherapy. Last, emphasis is given to the corrective strategies that modulate the cancerous CCND1 isoforms. Thus, it is highlighting significance of aberrant isoforms of CCND1 as targets for cancer therapy.

The Revelation of Continuously Organized, Co-Overexpressed Protein-Coding Genes with Roles in Cellular Communications in Breast Cancer

Many human cancers, including breast cancer, are polygenic and involve the co-dysregulation of multiple regulatory molecules and pathways. Though the overexpression of genes and amplified chromosomal regions have been closely linked in breast cancer, the notion of the co-upregulation of genes at a single locus remains poorly described. Here, we describe the co-overexpression of 34 continuously organized protein-coding genes with diverse functions at 8q.24.3(143437655-144326919) in breast and other cancer types, the CanCord34 genes. In total, 10 out of 34 genes have not been reported to be overexpressed in breast cancer. Interestingly, the overexpression of CanCord34 genes is not necessarily associated with genomic amplification and is independent of hormonal or HER2 status in breast cancer. CanCord34 genes exhibit diverse known and predicted functions, including enzymatic activities, cell viability, multipotency, cancer stem cells, and secretory activities, including extracellular vesicles. The co-overexpression of 33 of the CanCord34 genes in a multivariant analysis was correlated with poor survival among patients with breast cancer. The analysis of the genome-wide RNAi functional screening, cell dependency fitness, and breast cancer stem cell databases indicated that three diverse overexpressed CanCord34 genes, including a component of spliceosome PUF60, a component of exosome complex EXOSC4, and a ribosomal biogenesis factor BOP1, shared roles in cell viability, cell fitness, and stem cell phenotypes. In addition, 17 of the CanCord34 genes were found in the microvesicles (MVs) secreted from the mesenchymal stem cells that were primed with MDA-MB-231 breast cancer cells. Since these MVs were important in the chemoresistance and dedifferentiation of breast cancer cells into cancer stem cells, these findings highlight the significance of the CanCord34 genes in cellular communications. In brief, the persistent co-overexpression of CanCord34 genes with diverse functions can lead to the dysregulation of complementary functions in breast cancer. In brief, the present study provides new insights into the polygenic nature of breast cancer and opens new research avenues for basic, preclinical, and therapeutic studies in human cancer.

Association of 17q22 Amplicon Via Cell-Free DNA With Platinum Chemotherapy Response in Metastatic Triple-Negative Breast Cancer

PURPOSE

To determine whether specific somatic copy-number alterations detectable in circulating tumor DNA (ctDNA) from patients with metastatic triple-negative breast cancer (mTNBC) are associated with sensitivity to platinum chemotherapy.

MATERIALS AND METHODS

In this secondary analysis of a large cohort of patients with mTNBC whose ctDNA underwent ultralow-pass whole-genome sequencing, tumor fraction and somatic copy-number alterations were derived with the ichorCNA algorithm. Seventy-two patients were identified who had received a platinum-based chemotherapy regimen in the metastatic setting. Gene-level copy-number analyses were performed with GISTIC2.0. Cytobands were associated with progression-free survival (PFS) to platinum chemotherapy using Cox proportional hazards models. The Cancer Genome Atlas and Molecular Taxonomy of Breast Cancer International Consortium data sets were interrogated for frequency of significant cytobands in primary triple-negative breast cancer (pTNBC) tumors.

RESULTS

Among 71 evaluable patients, 17q21 and 17q22 amplifications were most strongly associated with improved PFS with platinum chemotherapy. There were no significant differences in clinicopathologic features or (neo)adjuvant chemotherapy among patients with 17q22 amplification. Patients with 17q22 amplification (n = 17) had longer median PFS with platinum (7.0 v 3.8 months; log-rank P = .015) than patients without 17q22 amplification (n = 54), an effect that remained significant in multivariable analyses (PFS hazard ratio 0.37; 95% CI, 0.16 to 0.84; P = .02). Among 39 patients who received the nonplatinum chemotherapy agent capecitabine, there was no association between 17q22 amplification and capecitabine PFS (log-rank P = .69). In The Cancer Genome Atlas and Molecular Taxonomy of Breast Cancer International Consortium, 17q22 amplification occurred in more than 20% of both pTNBC and mTNBC tumors, whereas 17q21 was more frequently amplified in mTNBC relative to pTNBC (16% v 8.1%, P = .015).

CONCLUSION

The 17q22 amplicon, detected by ctDNA, is associated with improved PFS with platinum chemotherapy in patients with mTNBC and warrants further investigation.

Rearrangement-mediated cis-regulatory alterations in advanced patient tumors reveal interactions with therapy

The global impact of somatic structural variants (SVs) on gene regulation in advanced tumors with complex treatment histories has been mostly uncharacterized. Here, using whole-genome and RNA sequencing from 570 recurrent or metastatic tumors, we report the altered expression of hundreds of genes in association with nearby SV breakpoints, including oncogenes and G-protein-coupled receptor-related genes such as PLEKHG2. A significant fraction of genes with SV-expression associations correlate with worse patient survival in primary and advanced cancers, including SRD5A1. In many instances, SV-expression associations involve retrotransposons being translocated near genes. High overall SV burden is associated with treatment with DNA alkylating agents or taxanes and altered expression of metabolism-associated genes. SV-expression associations within tumors from topoisomerase I inhibitor-treated patients include chromatin-related genes. Within anthracycline-treated tumors, SV breakpoints near chromosome 1p genes include PDE4B. Patient treatment and history can help understand the widespread SV-mediated cis-regulatory alterations found in cancer.

RRM2B Is Frequently Amplified Across Multiple Tumor Types: Implications for DNA Repair, Cellular Survival, and Cancer Therapy

RRM2B plays a crucial role in DNA replication, repair and oxidative stress. While germline RRM2B mutations have been implicated in mitochondrial disorders, its relevance to cancer has not been established. Here, using TCGA studies, we investigated RRM2B alterations in cancer. We found that RRM2B is highly amplified in multiple tumor types, particularly in MYC-amplified tumors, and is associated with increased RRM2B mRNA expression. We also observed that the chromosomal region 8q22.3–8q24, is amplified in multiple tumors, and includes RRM2B, MYC along with several other cancer-associated genes. An analysis of genes within this 8q-amplicon showed that cancers that have both RRM2B-amplified along with MYC have a distinct pattern of amplification compared to cancers that are unaltered or those that have amplifications in RRM2B or MYC only. Investigation of curated biological interactions revealed that gene products of the amplified 8q22.3–8q24 region have important roles in DNA repair, DNA damage response, oxygen sensing, and apoptosis pathways and interact functionally. Notably, RRM2B-amplified cancers are characterized by mutation signatures of defective DNA repair and oxidative stress, and at least RRM2B-amplified breast cancers are associated with poor clinical outcome. These data suggest alterations in RR2MB and possibly the interacting 8q-proteins could have a profound effect on regulatory pathways such as DNA repair and cellular survival, highlighting therapeutic opportunities in these cancers.

High mRNA expression of LY6 gene family is associated with overall survival outcome in pancreatic ductal adenocarcinoma

Pancreatic cancer ranks one of the worst in overall survival outcome with a 5 year survival rate being less than 10%. Pancreatic cancer faces unique challenges in its diagnosis and treatment, such as the lack of clinically validated biomarkers and the immensely immunosuppressive tumor microenvironment. Recently, the LY6 gene family has received increasing attention for its multi-faceted roles in cancer development, stem cell maintenance, immunomodulation, and association with more aggressive and hard-to-treat cancers. A detailed study of mRNA expression of LY6 gene family and its association with overall survival (OS) outcome in pancreatic cancers is lacking. We used publicly available clinical datasets to analyze the mRNA expression of a set of LY6 genes and its effect on OS outcome in the context of the tumor microenvironment and immunomodulation. We used web-based tools Kaplan-Meier Plotter, cBioPortal, Oncomine and R-programming to analyze copy number alterations, mRNA expression and its association with OS outcome in pancreatic cancer. These analyses demonstrated that high expression of LY6 genes is associated with OS and disease free survival (DFS) outcome. High expression of LY6 genes and their association with OS outcome is dependent on the composition of tumor microenvironment. Considering that LY6 proteins are anchored to the outer cell membrane or secreted, making them readily accessible, these findings highlight the potential of LY6 family members in the future of pancreatic cancer diagnosis and treatment.

The genomic landscape of metastasis in treatment-naïve breast cancer models

Metastasis remains the principle cause of mortality for breast cancer and presents a critical challenge because secondary lesions are often refractory to conventional treatments. While specific genetic alterations are tightly linked to primary tumor development and progression, the role of genetic alteration in the metastatic process is not well-understood. The theory of tumor evolution postulated by Peter Nowell in 1976 has yet to be proven in the context of metastasis. Therefore, in order to investigate how somatic evolution contributes to breast cancer metastasis, we performed exome, whole genome, and RNA sequencing of matched metastatic and primary tumors from pre-clinical mouse models of breast cancer. Here we show that in a treatment-naïve setting, recurrent single nucleotide variants and copy number variation, but not gene fusion events, play key metastasis-driving roles in breast cancer. For instance, we identified recurrent mutations in Kras, a known driver of colorectal and lung tumorigenesis that has not been previously implicated in breast cancer metastasis. However, in a set of in vivo proof-of-concept experiments we show that the Kras G12D mutation is sufficient to significantly promote metastasis using three syngeneic allograft models. The work herein confirms the existence of metastasis-driving mutations and presents a novel framework to identify actionable metastasis-targeted therapies.

TuBA: Tunable biclustering algorithm reveals clinically relevant tumor transcriptional profiles in breast cancer

BACKGROUND

Traditional clustering approaches for gene expression data are not well adapted to address the complexity and heterogeneity of tumors, where small sets of genes may be aberrantly co-expressed in specific subsets of tumors. Biclustering algorithms that perform local clustering on subsets of genes and conditions help address this problem. We propose a graph-based Tunable Biclustering Algorithm (TuBA) based on a novel pairwise proximity measure, examining the relationship of samples at the extremes of genes' expression profiles to identify similarly altered signatures.

RESULTS

TuBA's predictions are consistent in 3,940 breast invasive carcinoma samples from 3 independent sources, using different technologies for measuring gene expression (RNA sequencing and Microarray). More than 60% of biclusters identified independently in each dataset had significant agreement in their gene sets, as well as similar clinical implications. Approximately 50% of biclusters were enriched in the estrogen receptor-negative/HER2-negative (or basal-like) subtype, while >50% were associated with transcriptionally active copy number changes. Biclusters representing gene co-expression patterns in stromal tissue were also identified in tumor specimens.

CONCLUSIONS

TuBA offers a simple biclustering method that can identify biologically relevant gene co-expression signatures not captured by traditional unsupervised clustering approaches. It complements biclustering approaches that are designed to identify constant or coherent submatrices in gene expression datasets, and outperforms them in identifying a multitude of altered transcriptional profiles that are associated with observed genomic heterogeneity of diseased states in breast cancer, both within and across tumor subtypes, a promising step in understanding disease heterogeneity, and a necessary first step in individualized therapy.

Eukaryotic initiation factor 4E-binding protein as an oncogene in breast cancer

BACKGROUND

Eukaryotic Initiation Factor 4E-Binding Protein (EIF4EBP1, 4EBP1) is overexpressed in many human cancers including breast cancer, yet the role of 4EBP1 in breast cancer remains understudied. Despite the known role of 4EBP1 as a negative regulator of cap-dependent protein translation, 4EBP1 is predicted to be an essential driving oncogene in many cancer cell lines in vitro, and can act as a driver of cancer cell proliferation. EIF4EBP1 is located within the 8p11-p12 genomic locus, which is frequently amplified in breast cancer and is known to predict poor prognosis and resistance to endocrine therapy.

METHODS

Here we evaluated the effect of 4EBP1 targeting using shRNA knock-down of expression of 4EBP1, as well as response to the mTORC targeted drug everolimus in cell lines representing different breast cancer subtypes, including breast cancer cells with the 8p11-p12 amplicon, to better define a context and mechanism for oncogenic 4EBP1.

RESULTS

Using a genome-scale shRNA screen on the SUM panel of breast cancer cell lines, we found 4EBP1 to be a strong hit in the 8p11 amplified SUM-44 cells, which have amplification and overexpression of 4EBP1. We then found that knock-down of 4EBP1 resulted in dramatic reductions in cell proliferation in 8p11 amplified breast cancer cells as well as in other luminal breast cancer cell lines, but had little or no effect on the proliferation of immortalized but non-tumorigenic human mammary epithelial cells. Kaplan-Meier analysis of EIF4EBP1 expression in breast cancer patients demonstrated that overexpression of this gene was associated with reduced relapse free patient survival across all breast tumor subtypes.

CONCLUSIONS

These results are consistent with an oncogenic role of 4EBP1 in luminal breast cancer and suggests a role for this protein in cell proliferation distinct from its more well-known role as a regulator of cap-dependent translation.

Clinical Significance of the Wild Type p53-Induced Phosphatase 1 Expression in Invasive Breast Cancer

BACKGROUND

Wild type p53-induced phosphatase 1 (Wip1), encoded by the protein phosphatase magnesium dependent 1 delta (PPM1D), inhibits p53. PPM1D amplification has been reported in breast cancer. Breast cancer can sometimes develop without a tumor protein 53 (TP53) mutation. In these cases, the p53 pathway might be disrupted by alternative mechanisms, and Wip1 is reported to be a key molecule involved.

MATERIALS AND METHODS

Primary invasive ductal carcinoma specimens were obtained from 201 cases, for which archival tissue samples for immunohistochemistry were available. We evaluated Wip1 and p21 protein expression (201 cases), Wip1 mRNA expression (63 cases), PPM1D DNA copy number (71 cases) and TP53 status (36 cases) using available samples among the 201 cases, and analyzed their relationships with clinicopathological factors and prognosis.

RESULTS

The nuclear expression of Wip1 protein was positive in 21 cases (10.4%). The PPM1D DNA copy number was significantly correlated with Wip1 protein expression. All cases with PPM1D amplification by single-nucleotide polymorphism comparative genomic hybridization array showed positive nuclear Wip1 expression. Wip1 protein expression was positively correlated with p21 expression. The tumors with positive Wip1 and negative p21 expression showed the poorest prognosis among all tumor types.

CONCLUSION

The protein expression of Wip1 might be regulated by PPM1D amplification, independent of TP53 status. Positive Wip1 and negative p21 expression was associated with the poorest prognosis and suggests the loss of p53 function.

MiSTIC, an integrated platform for the analysis of heterogeneity in large tumour transcriptome datasets

Genome-wide transcriptome profiling has enabled non-supervised classification of tumours, revealing different sub-groups characterized by specific gene expression features. However, the biological significance of these subtypes remains for the most part unclear. We describe herein an interactive platform, Minimum Spanning Trees Inferred Clustering (MiSTIC), that integrates the direct visualization and comparison of the gene correlation structure between datasets, the analysis of the molecular causes underlying co-variations in gene expression in cancer samples, and the clinical annotation of tumour sets defined by the combined expression of selected biomarkers. We have used MiSTIC to highlight the roles of specific transcription factors in breast cancer subtype specification, to compare the aspects of tumour heterogeneity targeted by different prognostic signatures, and to highlight biomarker interactions in AML. A version of MiSTIC preloaded with datasets described herein can be accessed through a public web server (http://mistic.iric.ca); in addition, the MiSTIC software package can be obtained (github.com/iric-soft/MiSTIC) for local use with personalized datasets.

Development of mammary hyperplasia, dysplasia, and invasive ductal carcinoma in transgenic mice expressing the 8p11 amplicon oncogene NSD3

PURPOSE

NSD3 has been implicated as a candidate driver oncogene from the 8p11-p12 locus, and we have previously published evidence for its amplification and overexpression in human breast cancer. This aim of this study was to further characterize the transforming function of NSD3 in vivo.

METHODS

We generated a transgenic mouse model in which NSD3 gene expression was driven by the MMTV promoter and expressed in mammary epithelium of FVB mice. Mammary glands were fixed and whole mounts were stained with carmine to visualize gland structure. Mammary tumors were formalin-fixed, and paraffin embedded (FFPE) tumors were stained with hematoxylin and eosin.

RESULTS

Pups born to transgenic females were significantly underdeveloped compared to pups born to WT females due to a lactation defect in transgenic female mice. Whole mount analysis of the mammary glands of transgenic female mice revealed a profound defect in functional differentiation of mammary gland alveoli that resulted in the lactation defect. We followed parous and virgin NSD3 transgenic and control mice to 50 weeks of age and observed that several NSD3 parous females developed mammary tumors. Whole mount analysis of the mammary glands of tumor-bearing mice revealed numerous areas of mammary hyperplasia and ductal dysplasia. Histological analysis showed that mammary tumors were high-grade ductal carcinomas, and lesions present in other mammary glands exhibited features of alveolar hyperplasia, ductal dysplasia, and carcinoma in situ.

CONCLUSIONS

Our results are consistent with our previous studies and demonstrate that NSD3 is a transforming breast cancer oncogene.

Long non-coding RNA CCAT1 is overexpressed in oral squamous cell carcinomas and predicts poor prognosis

Oral squamous cell carcinoma (SCC) is the most common malignant tumor in India with 5-year survival rates totaling <50%. Recently, dysregulation of non-coding RNA was reported as a potential hallmark of carcinogenesis. Colon Cancer Associated Transcript 1 (CCAT1), an lncRNA located in chromosome 8q24, close to the c-Myc gene, has been reported to be overexpressed in many human cancers. In the present study, the authors analyzed the expression of CCAT1, c-Myc and the miRNAs miR155-5p, let7b-5p, miR490-3p and miR218-5p sponged by CCAT1 in 60 oral tumor and 8 normal tissue samples by reverse transcription-quantitative polymerase chain reaction. CCAT1 was significantly overexpressed in 27% (16/60) of oral tumors. Interestingly, a high level of c-Myc expression was observed in all CCAT1 overexpressing cases (P=0.0473). Furthermore, CCAT1 overexpression significantly downregulated miR155-5p (P=0.03) and let7b-5p (P<0.0001). Oral cancer cases expressing high level of CCAT1 (P=0.01) presented poor therapeutic outcome. To the best of the authors' knowledge, this is the first study to report the overexpression of the CCAT1 in oral SCCs, and the results suggested that CCAT1 overexpression may sponge miR155-5p and let7b-5p, and may account for poor treatment response.

A systematic comparison of copy number alterations in four types of female cancer

Background

Detection and localization of genomic alterations and breakpoints are crucial in cancer research. The purpose of this study was to investigate, in a methodological and biological perspective, different female, hormone-dependent cancers to identify common and diverse DNA aberrations, genes, and pathways.

Methods

In this work, we analyzed tissue samples from patients with breast (n = 112), ovarian (n = 74), endometrial (n = 84), or cervical (n = 76) cancer. To identify genomic aberrations, the Circular Binary Segmentation (CBS) and Piecewise Constant Fitting (PCF) algorithms were used and segmentation thresholds optimized. The Genomic Identification of Significant Targets in Cancer (GISTIC) algorithm was applied to the segmented data to identify significantly altered regions and the associated genes were analyzed by Ingenuity Pathway Analysis (IPA) to detect over-represented pathways and functions within the identified gene sets.

Results and Discussion

Analyses of high-resolution copy number alterations in four different female cancer types are presented. For appropriately adjusted segmentation parameters the two segmentation algorithms CBS and PCF performed similarly. We identified one region at 8q24.3 with focal aberrations that was altered at significant frequency across all four cancer types. Considering both, broad regions and focal peaks, three additional regions with gains at significant frequency were revealed at 1p21.1, 8p22, and 13q21.33, respectively. Several of these events involve known cancer-related genes, like PPP2R2A, PSCA, PTP4A3, and PTK2. In the female reproductive system (ovarian, endometrial, and cervix [OEC]), we discovered three common events: copy number gains at 5p15.33 and 15q11.2, further a copy number loss at 8p21.2. Interestingly, as many as 75% of the aberrations (75% amplifications and 86% deletions) identified by GISTIC were specific for just one cancer type and represented distinct molecular pathways.

Conclusions

Our results disclose that some prominent copy number changes are shared in the four examined female, hormone-dependent cancer whereas others are definitive to specific cancer types.

Electronic supplementary material

The online version of this article (doi:10.1186/s12885-016-2899-4) contains supplementary material, which is available to authorized users.

Emerging Role of Genomic Rearrangements in Breast Cancer: Applying Knowledge from Other Cancers

Significant advances in our knowledge of cancer genomes are rapidly changing the way we think about tumor biology and the heterogeneity of cancer. Recent successes in genomically-guided treatment approaches accompanied by more sophisticated sequencing techniques have paved the way for deeper investigation into the landscape of genomic rearrangements in cancer. While considerable research on solid tumors has focused on point mutations that directly alter the coding sequence of key genes, far less is known about the role of somatic rearrangements. With many recurring alterations observed across tumor types, there is an obvious need for functional characterization of these genomic biomarkers in order to understand their relevance to tumor biology, therapy, and prognosis. As personalized therapy approaches are turning toward genomic alterations for answers, these biomarkers will become increasingly relevant to the practice of precision medicine. This review discusses the emerging role of genomic rearrangements in breast cancer, with a particular focus on fusion genes. In addition, it raises several key questions on the therapeutic value of such rearrangements and provides a framework to evaluate their significance as predictive and prognostic biomarkers.

Evaluation of CCND1 amplification and CyclinD1 expression: diffuse and strong staining of CyclinD1 could have same predictive roles as CCND1 amplification in ER positive breast cancers

CCND1 is amplified in around 10-20% of primary breast cancers and preferentially occurs in ER positive tumors. Though CCND1 amplification was reported predicting poor response of adjuvant tamoxifen treatment and poor prognosis in ER positive breast cancers, there were controversial data regarding the predicting value of CyclinD1 protein overexpression. In this study, we detected CyclinD1 expression using immunohistochemistry and CCND1 gene copy number using fluorescence in situ hybridization (FISH) in 355 invasive breast cancers with foci ductal carcinoma in situ (DCIS). CCND1 amplification was founded in 52 (14.6%) cases all of which showed moderate to strong CyclinD1 expression. However, majority of CCND1- tumors exhibited mild to moderate CyclinD1 staining. There were identical alterations in DCIS and the invasive lesions within the same tumor. CCND1 amplification was positively correlated with ER, PR and lymph node status (P<0.001) while negatively correlated with HER-2 amplification and p53 status (P<0.05). The majority of the CCND1 amplification/high CyclinD1 breast cancers were luminal B type while basal-like type often lost the expression of this protein. The ROC curve analysis showed that a cut-off point at which the immunostaining score of CyclinD1 is 6.5 could predict CCND1 gene amplification in breast cancer. This study indicated loss expression of CyclinD1 might be an important event in the tumorigenesis in basal-like breast cancers. Further, we confirmed an optimal cut-off point of immunostaining scores of CyclinD1 protein which could be used to predict the status of CCND1 gene and identify a subgroup of ER positive breast cancers with poor response to endocrine agents.

The Major Prognostic Features of Nuclear Receptor NR5A2 in Infiltrating Ductal Breast Carcinomas

Background. Gene expression profiles of 181 breast cancer samples were analyzed to identify prognostic features of nuclear receptors NR5A1 and NR5A2 based upon their associated transcriptional networks. Methods. A supervised network analysis approach was used to build the NR5A-mediated transcriptional regulatory network. Other bioinformatic tools and statistical methods were utilized to confirm and extend results from the network analysis methodology. Results. NR5A2 expression is a negative factor in breast cancer prognosis in both ER(-) and ER(-)/ER(+) mixed cohorts. The clinical and cohort significance of NR5A2-mediated transcriptional activities indicates that it may have a significant role in attenuating grade development and cancer related signal transduction pathways. NR5A2 signature that conditions poor prognosis was identified based upon results from 15 distinct probes. Alternatively, the expression of NR5A1 predicts favorable prognosis when concurrent NR5A2 expression is low. A favorable signature of eight transcription factors mediated by NR5A1 was also identified. Conclusions. Correlation of poor prognosis and NR5A2 activity is identified by NR5A2-mediated 15-gene signature. NR5A2 may be a potential drug target for treating a subset of breast cancer tumors across breast cancer subtypes, especially ER(-) breast tumors. The favorable prognostic feature of NR5A1 is predicted by NR5A1-mediated 8-gene signature.

Identification of Copy Number Aberrations in Breast Cancer Subtypes Using Persistence Topology

DNA copy number aberrations (CNAs) are of biological and medical interest because they help identify regulatory mechanisms underlying tumor initiation and evolution. Identification of tumor-driving CNAs (driver CNAs) however remains a challenging task, because they are frequently hidden by CNAs that are the product of random events that take place during tumor evolution. Experimental detection of CNAs is commonly accomplished through array comparative genomic hybridization (aCGH) assays followed by supervised and/or unsupervised statistical methods that combine the segmented profiles of all patients to identify driver CNAs. Here, we extend a previously-presented supervised algorithm for the identification of CNAs that is based on a topological representation of the data. Our method associates a two-dimensional (2D) point cloud with each aCGH profile and generates a sequence of simplicial complexes, mathematical objects that generalize the concept of a graph. This representation of the data permits segmenting the data at different resolutions and identifying CNAs by interrogating the topological properties of these simplicial complexes. We tested our approach on a published dataset with the goal of identifying specific breast cancer CNAs associated with specific molecular subtypes. Identification of CNAs associated with each subtype was performed by analyzing each subtype separately from the others and by taking the rest of the subtypes as the control. Our results found a new amplification in 11q at the location of the progesterone receptor in the Luminal A subtype. Aberrations in the Luminal B subtype were found only upon removal of the basal-like subtype from the control set. Under those conditions, all regions found in the original publication, except for 17q, were confirmed; all aberrations, except those in chromosome arms 8q and 12q were confirmed in the basal-like subtype. These two chromosome arms, however, were detected only upon removal of three patients with exceedingly large copy number values. More importantly, we detected 10 and 21 additional regions in the Luminal B and basal-like subtypes, respectively. Most of the additional regions were either validated on an independent dataset and/or using GISTIC. Furthermore, we found three new CNAs in the basal-like subtype: a combination of gains and losses in 1p, a gain in 2p and a loss in 14q. Based on these results, we suggest that topological approaches that incorporate multiresolution analyses and that interrogate topological properties of the data can help in the identification of copy number changes in cancer.

Wip1 phosphatase in breast cancer

Understanding the factors contributing to tumor initiation, progression and evolution is of paramount significance. Among them, wild-type p53-induced phosphatase 1 (Wip1) is emerging as an important oncogene by virtue of its negative control on several key tumor suppressor pathways. Originally discovered as a p53-regulated gene, Wip1 has been subsequently found amplified and more recently mutated in a significant fraction of human cancers including breast tumors. Recent development in the field further uncovered the utility of anti-Wip1-directed therapies in delaying tumor onset or in reducing the tumor burden. Furthermore, Wip1 could be an important factor that contributes to tumor heterogeneity, suggesting that its inhibition may decrease the rate of cancer evolution. These effects depend on several signaling pathways modulated by Wip1 phosphatase in a spatial and temporal manner. In this review we discuss the recent development in understanding how Wip1 contributes to tumorigenesis with its relevance to breast cancer.

Amplification of distant estrogen response elements deregulates target genes associated with tamoxifen resistance in breast cancer

A causal role of gene amplification in tumorigenesis is well known, whereas amplification of DNA regulatory elements as an oncogenic driver remains unclear. In this study, we integrated next-generation sequencing approaches to map distant estrogen response elements (DEREs) that remotely control the transcription of target genes through chromatin proximity. Two densely mapped DERE regions located on chromosomes 17q23 and 20q13 were frequently amplified in estrogen receptor-α-positive luminal breast cancer. These aberrantly amplified DEREs deregulated target gene expression potentially linked to cancer development and tamoxifen resistance. Progressive accumulation of DERE copies was observed in normal breast progenitor cells chronically exposed to estrogenic chemicals. These findings may extend to other DNA regulatory elements, the amplification of which can profoundly alter target transcriptome during tumorigenesis.

Biomolecular events in cancer revealed by attractor metagenes

Mining gene expression profiles has proven valuable for identifying signatures serving as surrogates of cancer phenotypes. However, the similarities of such signatures across different cancer types have not been strong enough to conclude that they represent a universal biological mechanism shared among multiple cancer types. Here we present a computational method for generating signatures using an iterative process that converges to one of several precise attractors defining signatures representing biomolecular events, such as cell transdifferentiation or the presence of an amplicon. By analyzing rich gene expression datasets from different cancer types, we identified several such biomolecular events, some of which are universally present in all tested cancer types in nearly identical form. Although the method is unsupervised, we show that it often leads to attractors with strong phenotypic associations. We present several such multi-cancer attractors, focusing on three that are prominent and sharply defined in all cases: a mesenchymal transition attractor strongly associated with tumor stage, a mitotic chromosomal instability attractor strongly associated with tumor grade, and a lymphocyte-specific attractor.

Cancer is known to be characterized by several unifying biological capabilities or “hallmarks.” However, attempts to computationally identify patterns, such as gene expression signatures, shared across many different cancer types have been largely unsuccessful. A typical approach has been to classify samples into mutually exclusive subtypes, each of which is characterized by a particular gene signature. Although occasional similarities of such signatures in different cancer types exist, these similarities have not been sufficiently strong to conclude that they reflect the same biological event. By contrast, we have developed a computational methodology that has identified some signatures of co-expressed genes exhibiting remarkable similarity across many different cancer types. These signatures appear as stable “attractors” of an iterative computational procedure that tends to collect mutually associated genes, so that its convergence can point to the core (“heart”) of the underlying biological co-expression mechanism. One of these “pan-cancer” attractors corresponds to a transdifferentiation of cancer cells empowering them with invasiveness and motility. Another represents a mitotic chromosomal instability of cancer cells. A third attractor is lymphocyte-specific.

Gene expression signatures that predict outcome of tamoxifen-treated estrogen receptor-positive, high-risk, primary breast cancer patients: a DBCG study

Background

Tamoxifen significantly improves outcome for estrogen receptor-positive (ER+) breast cancer, but the 15-year recurrence rate remains 30%. The aim of this study was to identify gene profiles that accurately predicted the outcome of ER+ breast cancer patients who received adjuvant Tamoxifen mono-therapy.

Methodology/Principal Findings

Post-menopausal breast cancer patients diagnosed no later than 2002, being ER+ as defined by >1% IHC staining and having a frozen tumor sample with >50% tumor content were included. Tumor samples from 108 patients treated with adjuvant Tamoxifen were analyzed for the expression of 59 genes using quantitative-PCR. End-point was clinically verified recurrence to distant organs or ipsilateral breast. Gene profiles were identified using a model building procedure based on conditional logistic regression and leave-one-out cross-validation, followed by a non-parametric bootstrap (1000x re-sampling). The optimal profiles were further examined in 5 previously-reported datasets containing similar patient populations that were either treated with Tamoxifen or left untreated (n = 623). Three gene signatures were identified, the strongest being a 2-gene combination of BCL2-CDKN1A, exhibiting an accuracy of 75% for prediction of outcome. Independent examination using 4 previously-reported microarray datasets of Tamoxifen-treated patient samples (n = 503) confirmed the potential of BCL2-CDKN1A. The predictive value was further determined by comparing the ability of the genes to predict recurrence in an additional, previously-published, cohort consisting of Tamoxifen-treated (n = 58, p = 0.015) and untreated patients (n = 62, p = 0.25).

Conclusions/Significance

A novel gene expression signature predictive of outcome of Tamoxifen-treated patients was identified. The validation suggests that BCL2-CDKN1A exhibit promising predictive potential.

Major chromosomal breakpoint intervals in breast cancer co-localize with differentially methylated regions

Solid tumors exhibit chromosomal rearrangements resulting in gain or loss of multiple chromosomal loci (copy number variation, or CNV), and translocations that occasionally result in the creation of novel chimeric genes. In the case of breast cancer, although most individual tumors each have unique CNV landscape, the breakpoints, as measured over large datasets, appear to be non-randomly distributed in the genome. Breakpoints show a significant regional concentration at genomic loci spanning perhaps several megabases. The proximal cause of these breakpoint concentrations is a subject of speculation, but is, as yet, largely unknown. To shed light on this issue, we have performed a bio-statistical analysis on our previously published data for a set of 119 breast tumors and normal controls (Wiedswang et al., 2003), where each sample has both high-resolution CNV and methylation data. The method examined the distribution of closeness of breakpoint regions with differentially methylated regions (DMR), coupled with additional genomic parameters, such as repeat elements and designated "fragile sites" in the reference genome. Through this analysis, we have identified a set of 93 regional loci called breakpoint enriched DMR (BEDMRs) characterized by altered DNA methylation in cancer compared to normal cells that are associated with frequent breakpoint concentrations within a distance of 1 Mb. BEDMR loci are further associated with local hypomethylation (66%), concentrations of the Alu SINE repeats within 3 Mb (35% of the cases), and tend to occur near a number of cancer related genes such as the protocadherins, AKT1, DUB3, GAB2. Furthermore, BEDMRs seem to deregulate members of the histone gene family and chromatin remodeling factors, e.g., JMJD1B, which might affect the chromatin structure and disrupt coordinate signaling and repair. From this analysis we propose that preference for chromosomal breakpoints is related to genome structure coupled with alterations in DNA methylation and hence, chromatin structure, associated with tumorigenesis.