ROCK: a resource for integrative breast cancer data analysis

Tumor cell intrinsic and extrinsic features predict prognosis in estrogen receptor positive breast cancer

Although estrogen-receptor-positive (ER+) breast cancer is generally associated with favorable prognosis, clinical outcome varies substantially among patients. Genomic assays have been developed and applied to predict patient prognosis for personalized treatment. We hypothesize that the recurrence risk of ER+ breast cancer patients is determined by both genomic mutations intrinsic to tumor cells and extrinsic immunological features in the tumor microenvironment. Based on the Cancer Genome Atlas (TCGA) breast cancer data, we identified the 72 most common genomic aberrations (including gene mutations and indels) in ER+ breast cancer and defined sample-specific scores that systematically characterized the deregulated pathways intrinsic to tumor cells. To further consider tumor cell extrinsic features, we calculated immune infiltration scores for six major immune cell types. Many individual intrinsic features are predictive of patient prognosis in ER+ breast cancer, and some of them achieved comparable accuracy with the Oncotype DX assay. In addition, statistical learning models that integrated these features predicts the recurrence risk of patients with significantly better performance than the Oncotype DX assay (our optimized random forest model AUC = 0.841, Oncotype DX model AUC = 0.792, p = 0.04). As a proof-of-concept, our study indicates the great potential of genomic and immunological features in prognostic prediction for improving breast cancer precision medicine. The framework introduced in this work can be readily applied to other cancers.

A framework to predict the applicability of Oncotype DX, MammaPrint, and E2F4 gene signatures for improving breast cancer prognostic prediction

To improve cancer precision medicine, prognostic and predictive biomarkers are critically needed to aid physicians in deciding treatment strategies in a personalized fashion. Due to the heterogeneous nature of cancer, most biomarkers are expected to be valid only in a subset of patients. Furthermore, there is no current approach to determine the applicability of biomarkers. In this study, we propose a framework to improve the clinical application of biomarkers. As part of this framework, we develop a clinical outcome prediction model (CPM) and a predictability prediction model (PPM) for each biomarker and use these models to calculate a prognostic score (P-score) and a confidence score (C-score) for each patient. Each biomarker’s P-score indicates its association with patient clinical outcomes, while each C-score reflects the biomarker applicability of the biomarker’s CPM to a patient and therefore the confidence of the clinical prediction. We assessed the effectiveness of this framework by applying it to three biomarkers, Oncotype DX, MammaPrint, and an E2F4 signature, which have been used for predicting patient response, pathologic complete response versus residual disease to neoadjuvant chemotherapy (a classification problem), and recurrence-free survival (a Cox regression problem) in breast cancer, respectively. In both applications, our analyses indicated patients with higher C scores were more likely to be correctly predicted by the biomarkers, indicating the effectiveness of our framework. This framework provides a useful approach to develop and apply biomarkers in the context of cancer precision medicine.

CCR5 is a prognostic biomarker and an immune regulator for triple negative breast cancer

This study aims to explore the clinical implications and potential mechanistic functions of CCR5 in triple negative breast cancer. Briefly, we demonstrated that CCR5 is overexpressed in TNBC and is associated with better prognosis of TNBC. CCR5 expression is positively correlated with tumor immune cell infiltration and tumor immune response related pathways. Multi-omics data analyses identified CCR5 associated genomic and proteomic changes. CCR5 overexpression was associated with better overall survival in TNBC patients with TP53 mutation. We also summarized the latest findings on ICB efficacy related genes and explored the association between CCR5 and those genes. These results indicated that CCR5 is a potential tumor suppressor gene and individualized therapeutic strategy could be established based on multi-omics background and expression pattern of ICB related genes. In conclusion, CCR5 is associated with better survival of TNBC patients with TP53 mutation, which may exert its roles through tumor immune environment.

The androgen receptor is a tumor suppressor in estrogen receptor-positive breast cancer

The role of the androgen receptor (AR) in estrogen receptor (ER)-α-positive breast cancer is controversial, constraining implementation of AR-directed therapies. Using a diverse, clinically relevant panel of cell-line and patient-derived models, we demonstrate that AR activation, not suppression, exerts potent antitumor activity in multiple disease contexts, including resistance to standard-of-care ER and CDK4/6 inhibitors. Notably, AR agonists combined with standard-of-care agents enhanced therapeutic responses. Mechanistically, agonist activation of AR altered the genomic distribution of ER and essential co-activators (p300, SRC-3), resulting in repression of ER-regulated cell cycle genes and upregulation of AR target genes, including known tumor suppressors. A gene signature of AR activity positively predicted disease survival in multiple clinical ER-positive breast cancer cohorts. These findings provide unambiguous evidence that AR has a tumor suppressor role in ER-positive breast cancer and support AR agonism as the optimal AR-directed treatment strategy, revealing a rational therapeutic opportunity.

Clinical Implications of Chromosomal Instability (CIN) and Kinetochore Abnormalities in Breast Cancers

Genetic instability is a defining property of cancer cells and is the basis of various lesions including point mutations, copy number alterations and translocations. Chromosomal instability (CIN) is part of the genetic instability of cancer and consists of copy number alterations in whole or parts of cancer cell chromosomes. CIN is observed in differing degrees in most cancers. In breast cancer, CIN is commonly part of the genomic landscape of the disease and has a higher incidence in aggressive sub-types. Tumor suppressors that are commonly mutated or disabled in cancer, such as p53 and pRB, play roles in protection against CIN, and as a result, their dysfunction contributes to the establishment or tolerance of CIN. Several structural and regulatory proteins of the centromeres and kinetochore, the complex structure that is responsible for the correct distribution of genetic material in the daughter cells during mitosis, are direct or, mostly, indirect transcription targets of p53 and pRB. Thus, despite the absence of structural defects in genes encoding for centromere and kinetochore components, dysfunction of these tumor suppressors may have profound implications for the correct function of the mitotic apparatus contributing to CIN. CIN and its prognostic and therapeutic implications in breast cancer are discussed in this article.

Increased expression of the HDAC9 gene is associated with antiestrogen resistance of breast cancers

Estrogens play a pivotal role in breast cancer etiology, and endocrine therapy remains the main first line treatment for estrogen receptor‐alpha (ERα)‐positive breast cancer. ER are transcription factors whose activity is finely regulated by various regulatory complexes, including histone deacetylases (HDACs). Here, we investigated the role of HDAC9 in ERα signaling and response to antiestrogens in breast cancer cells. Various Michigan Cancer Foundation‐7 (MCF7) breast cancer cell lines that overexpress class IIa HDAC9 or that are resistant to the partial antiestrogen 4‐hydroxy‐tamoxifen (OHTam) were used to study phenotypic changes in response to ER ligands by using transcriptomic and gene set enrichment analyses. Kaplan–Meier survival analyses were performed using public transcriptomic datasets from human breast cancer biopsies. In MCF7 breast cancer cells, HDAC9 decreased ERα mRNA and protein expression and inhibited its transcriptional activity. Conversely, HDAC9 mRNA was strongly overexpressed in OHTam‐resistant MCF7 cells and in ERα‐negative breast tumor cell lines. Moreover, HDAC9‐overexpressing cells were less sensitive to OHTam antiproliferative effects compared with parental MCF7 cells. Several genes (including MUC1, SMC3 and S100P) were similarly deregulated in OHTam‐resistant and in HDAC9‐overexpressing MCF7 cells. Finally, HDAC9 expression was positively associated with genes upregulated in endocrine therapy‐resistant breast cancers and high HDAC9 levels were associated with worse prognosis in patients treated with OHTam. These results demonstrate the complex interactions of class IIa HDAC9 with ERα signaling in breast cancer cells and its effect on the response to hormone therapy.

Gasdermin B expression predicts poor clinical outcome in HER2-positive breast cancer

Around, 30–40% of HER2-positive breast cancers do not show substantial clinical benefit from the targeted therapy and, thus, the mechanisms underlying resistance remain partially unknown. Interestingly, ERBB2 is frequently co-amplified and co-expressed with neighbour genes that may play a relevant role in this cancer subtype. Here, using an in silico analysis of data from 2,096 breast tumours, we reveal a significant correlation between Gasdermin B (GSDMB) gene (located 175 kilo bases distal from ERBB2) expression and the pathological and clinical parameters of poor prognosis in HER2-positive breast cancer. Next, the analysis of three independent cohorts (totalizing 286 tumours) showed that approximately 65% of the HER2-positive cases have GSDMB gene amplification and protein over-expression. Moreover, GSDMB expression was also linked to poor therapeutic responses in terms of lower relapse free survival and pathologic complete response as well as positive lymph node status and the development of distant metastasis under neoadjuvant and adjuvant treatment settings, respectively. Importantly, GSDMB expression promotes survival to trastuzumab in different HER2-positive breast carcinoma cells, and is associated with trastuzumab resistance phenotype in vivo in Patient Derived Xenografts. In summary, our data identifies the ERBB2 co-amplified and co-expressed gene GSDMB as a critical determinant of poor prognosis and therapeutic response in HER2-positive breast cancer.

Computational Investigation of Homologous Recombination DNA Repair Deficiency in Sporadic Breast Cancer

BRCAness has important implications in the management and treatment of patients with breast and ovarian cancer. In this study, we propose a computational framework to measure the BRCAness of breast and ovarian tumor samples based on their gene expression profiles. We define a characteristic profile for BRCAness by comparing gene expression differences between BRCA1/2 mutant familial tumors and sporadic breast cancer tumors while adjusting for relevant clinical factors. With this BRCAness profile, our framework calculates sample-specific BRCA scores, which indicates homologous recombination (HR)-mediated DNA repair pathway activity of samples. We found that in sporadic breast cancer high BRCAness score is associated with aberrant copy number of HR genes rather than somatic mutation and other genomic features. Moreover, we observed significant correlations of BRCA score with genome instability and neoadjuvant chemotherapy. More importantly, BRCA score provides significant prognostic value in both breast and ovarian cancers after considering established clinical variables. In summary, the inferred BRCAness from our framework can be used as a robust biomarker for the prediction of prognosis and treatment response in breast and ovarian cancers.

Exceptionally high UBE2C expression is a unique phenomenon in basal-like type breast cancer and is regulated by BRCA1

Ubiquitin-conjugating enzyme 2C (UBE2C) is overexpressed in various types of cancer, leading to poor outcomes and drug resistance. UBE2C may also have a critical role in phenotypes associated with poor prognosis in breast cancer; however, the relationship between UBE2C expression and clinical outcome in breast cancer subtypes has not previously been investigated. We firstly analyzed breast cancer patient data and immunohistochemistry of breast cancer patient samples. We demonstrated that UBE2C was associated with poor prognosis in breast cancer, particularly basal-like breast cancer, a subtype with aggressive clinical features. Interestingly, we found that there was a close relationship between the expression of BRCA1 and UBE2C in the MCF-7 and MDA-MB-231 breast cancer cell lines. Upregulation of BRCA1 could inhibit the expression of UBE2C. In cells with BRCA1 silenced down, expression of UBE2C was obviously increased, with a concurrent decrease in cellular sensitivity to doxorubicin. Suppression of UBE2C expression by RNA interference led to decrease the mRNA expressions of BCRP, MRP1 and P-gp in doxorubicin-treated MDA-MB-231 cells. Moreover, treatment with 1μg/ml doxorubicin led to increased expression of UBE2C. The results show high expression of UBE2C is a potential prognostic factor of poor outcome in basal-like breast cancer. Moreover, loss of BRCA1 function results in an increase in UBE2C expression and chemical resistance to doxorubicin in breast cancer cells.

Basal-like breast cancer: molecular profiles, clinical features and survival outcomes

BACKGROUND

Basal-like constitutes an important molecular subtype of breast cancer characterised by an aggressive behaviour and a limited therapy response. The outcome of patients within this subtype is, however, divergent. Some individuals show an increased risk of dying in the first five years, and others a long-term survival of over ten years after the diagnosis. In this study, we aim at identifying markers associated with basal-like patients' survival and characterising subgroups with distinct disease outcome.

METHODS

We explored the genomic and transcriptomic profiles of 351 basal-like samples from the METABRIC and ROCK data sets. Two selection methods, labelled Differential and Survival filters, were employed to determine genes/probes that are differentially expressed in tumour and control samples, and are associated with overall survival. These probes were further used to define molecular subgroups, which vary at the microRNA level and in DNA copy number.

RESULTS

We identified the expression signature of 80 probes that distinguishes between two basal-like subgroups with distinct clinical features and survival outcomes. Genes included in this list have been mainly linked to cancer immune response, epithelial-mesenchymal transition and cell cycle. In particular, high levels of CXCR6, HCST, C3AR1 and FPR3 were found in Basal I; whereas HJURP, RRP12 and DNMT3B appeared over-expressed in Basal II. These genes exhibited the highest betweenness centrality and node degree values and play a key role in the basal-like breast cancer differentiation. Further molecular analysis revealed 17 miRNAs correlated to the subgroups, including hsa-miR-342-5p, -150, -155, -200c and -17. Additionally, increased percentages of gains/amplifications were detected on chromosomes 1q, 3q, 8q, 10p and 17q, and losses/deletions on 4q, 5q, 8p and X, associated with reduced survival.

CONCLUSIONS

The proposed signature supports the existence of at least two subgroups of basal-like breast cancers with distinct disease outcome. The identification of patients at a low risk may impact the clinical decisions-making by reducing the prescription of high-dose chemotherapy and, consequently, avoiding adverse effects. The recognition of other aggressive features within this subtype may be also critical for improving individual care and for delineating more effective therapies for patients at high risk.

Pentraxin-3 is a PI3K signaling target that promotes stem cell-like traits in basal-like breast cancers

Basal-like breast cancers (BLBCs) exhibit hyperactivation of the phosphoinositide 3-kinase (PI3K) signaling pathway because of the frequent mutational activation of the PIK3CA catalytic subunit and the genetic loss of its negative regulators PTEN (phosphatase and tensin homolog) and INPP4B (inositol polyphosphate-4-phosphatase type II). However, PI3K inhibitors have had limited clinical efficacy in BLBC management because of compensatory amplification of PI3K downstream signaling loops. Therefore, identification of critical PI3K mediators is paramount to the development of effective BLBC therapeutics. Using transcriptomic analysis of activated PIK3CA-expressing BLBC cells, we identified the gene encoding the humoral pattern recognition molecule pentraxin-3 (PTX3) as a critical target of oncogenic PI3K signaling. We found that PTX3 abundance is stimulated, in part, through AKT- and nuclear factor κB (NF-κB)-dependent pathways and that presence of PTX3 is necessary for PI3K-induced stem cell-like traits. We further showed that PTX3 expression is greater in tumor samples from patients with BLBC and that it is prognostic of poor patient survival. Our results thus reveal PTX3 as a newly identified PI3K-regulated biomarker and a potential therapeutic target in BLBC.

Contextual Refinement of Regulatory Targets Reveals Effects on Breast Cancer Prognosis of the Regulome

Gene expression regulators, such as transcription factors (TFs) and microRNAs (miRNAs), have varying regulatory targets based on the tissue and physiological state (context) within which they are expressed. While the emergence of regulator-characterizing experiments has inferred the target genes of many regulators across many contexts, methods for transferring regulator target genes across contexts are lacking. Further, regulator target gene lists frequently are not curated or have permissive inclusion criteria, impairing their use. Here, we present a method called iterative Contextual Transcriptional Activity Inference of Regulators (icTAIR) to resolve these issues. icTAIR takes a regulator’s previously-identified target gene list and combines it with gene expression data from a context, quantifying that regulator’s activity for that context. It then calculates the correlation between each listed target gene’s expression and the quantitative score of regulatory activity, removes the uncorrelated genes from the list, and iterates the process until it derives a stable list of refined target genes. To validate and demonstrate icTAIR’s power, we use it to refine the MSigDB c3 database of TF, miRNA and unclassified motif target gene lists for breast cancer. We then use its output for survival analysis with clinicopathological multivariable adjustment in 7 independent breast cancer datasets covering 3,430 patients. We uncover many novel prognostic regulators that were obscured prior to refinement, in particular NFY, and offer a detailed look at the composition and relationships among the breast cancer prognostic regulome. We anticipate icTAIR will be of general use in contextually refining regulator target genes for discoveries across many contexts. The icTAIR algorithm can be downloaded from https://github.com/icTAIR.

Gene expression regulators, such as transcription factors and microRNAs, are critical actors in cellular physiology and pathophysiology and act by modulating the expression levels of sets of target genes. Given their significance, numerous experiments have sought to characterize the specific target genes of specific regulators, which in turn has led to regulator target gene list databases. Unfortunately, these lists are plagued by poor curation and validation. Further, all lists suffer from the fundamental issue that regulator targets vary across tissue type and physiological state, or “context”, making them poor for conducting downstream, context-specific analyses. To address this issue, here we present a method called icTAIR that contextually-refines regulator target gene lists. To demonstrate its value, we use icTAIR to take the largest-available database of regulator target gene lists, refine it for the breast cancer context, and use both the pre-refined and refined lists for downstream survival analyses in over 3,400 tumors. We find that icTAIR improves the statistical power of the analyses by multiple orders of magnitude. This in turn lets us map the relational network of breast cancer regulators and identify regulators with prognostic effects even after clinicopathological adjustment. We anticipate icTAIR will be broadly useful in regulator studies.

Loss of CtIP disturbs homologous recombination repair and sensitizes breast cancer cells to PARP inhibitors

Breast cancer is one of the leading causes of death worldwide, and therefore, new and improved approaches for the treatment of breast cancer are desperately needed. CtIP (RBBP8) is a multifunctional protein that is involved in various cellular functions, including transcription, DNA replication, DNA repair and the G1 and G2 cell cycle checkpoints. CtIP plays an important role in homologous recombination repair by interacting with tumor suppressor protein BRCA1. Here, we analyzed the expression profile of CtIP by data mining using published microarray data sets. We found that CtIP expression is frequently decreased in breast cancer patients, and the patient group with low-expressing CtIP mRNA is associated with a significantly lower survival rate. The knockdown of CtIP in breast cancer MCF7 cells reduced Rad51 foci numbers and enhanced f H2AX foci formation after f-irradiation, suggesting that deficiency of CtIP decreases homologous recombination repair and delays DNA double strand break repair. To explore the effect of CtIP on PARP inhibitor therapy for breast cancer, CtIP-depleted MCF7 cells were treated with PARP inhibitor olaparib (AZD2281) or veliparib (ABT-888). As in BRCA mutated cells, PARP inhibitors showed cytotoxicity to CtIP-depleted cells by preventing cells from repairing DNA damage, leading to decreased cell viability. Further, a xenograft tumor model in mice with MCF7 cells demonstrated significantly increased sensitivity towards PARP inhibition under CtIP deficiency. In summary, this study shows that low level of CtIP expression is associated with poor prognosis in breast cancer, and provides a rationale for establishing CtIP expression as a biomarker of PARP inhibitor response, and consequently offers novel therapeutic options for a significant subset of patients.

Application of RNAi-induced gene expression profiles for prognostic prediction in breast cancer

Homologous recombination (HR) is the primary pathway for repairing double-strand DNA breaks implicating in the development of cancer. RNAi-based knockdowns of BRCA1 and RAD51 in this pathway have been performed to investigate the resulting transcriptomic profiles. Here we propose a computational framework to utilize these profiles to calculate a score, named RNA-Interference derived Proliferation Score (RIPS), which reflects cell proliferation ability in individual breast tumors. RIPS is predictive of breast cancer classes, prognosis, genome instability, and neoadjuvant chemosensitivity. This framework directly translates the readout of knockdown experiments into potential clinical applications and generates a robust biomarker in breast cancer.

Centromere and kinetochore gene misexpression predicts cancer patient survival and response to radiotherapy and chemotherapy

Chromosomal instability (CIN) is a hallmark of cancer that contributes to tumour heterogeneity and other malignant properties. Aberrant centromere and kinetochore function causes CIN through chromosome missegregation, leading to aneuploidy, rearrangements and micronucleus formation. Here we develop a Centromere and kinetochore gene Expression Score (CES) signature that quantifies the centromere and kinetochore gene misexpression in cancers. High CES values correlate with increased levels of genomic instability and several specific adverse tumour properties, and prognosticate poor patient survival for breast and lung cancers, especially early-stage tumours. They also signify high levels of genomic instability that sensitize cancer cells to additional genotoxicity. Thus, the CES signature forecasts patient response to adjuvant chemotherapy or radiotherapy. Our results demonstrate the prognostic and predictive power of the CES, suggest a role for centromere misregulation in cancer progression, and support the idea that tumours with extremely high CIN are less tolerant to specific genotoxic therapies.

Centromeres and kinetochores are important in maintaining chromosomal stability. Here, the authors show that overexpression of a subset of centromere and kinetochore genes is associated with chromosomal instability and mutation burden in cancer, and predict patient survival and response to genotoxic therapies.

Extensive Transcriptomic and Genomic Analysis Provides New Insights about Luminal Breast Cancers

Despite constituting approximately two thirds of all breast cancers, the luminal A and B tumours are poorly classified at both clinical and molecular levels. There are contradictory reports on the nature of these subtypes: some define them as intrinsic entities, others as a continuum. With the aim of addressing these uncertainties and identifying molecular signatures of patients at risk, we conducted a comprehensive transcriptomic and genomic analysis of 2,425 luminal breast cancer samples. Our results indicate that the separation between the molecular luminal A and B subtypes—per definition—is not associated with intrinsic characteristics evident in the differentiation between other subtypes. Moreover, t-SNE and MST-kNN clustering approaches based on 10,000 probes, associated with luminal tumour initiation and/or development, revealed the close connections between luminal A and B tumours, with no evidence of a clear boundary between them. Thus, we considered all luminal tumours as a single heterogeneous group for analysis purposes. We first stratified luminal tumours into two distinct groups by their HER2 gene cluster co-expression: HER2-amplified luminal and ordinary-luminal. The former group is associated with distinct transcriptomic and genomic profiles, and poor prognosis; it comprises approximately 8% of all luminal cases. For the remaining ordinary-luminal tumours we further identified the molecular signature correlated with disease outcomes, exhibiting an approximately continuous gene expression range from low to high risk. Thus, we employed four virtual quantiles to segregate the groups of patients. The clinico-pathological characteristics and ratios of genomic aberrations are concordant with the variations in gene expression profiles, hinting at a progressive staging. The comparison with the current separation into luminal A and B subtypes revealed a substantially improved survival stratification. Concluding, we suggest a review of the definition of luminal A and B subtypes. A proposition for a revisited delineation is provided in this study.

Integrative analysis of breast cancer reveals prognostic haematopoietic activity and patient-specific immune response profiles

Transcriptional programmes active in haematopoietic cells enable a variety of functions including dedifferentiation, innate immunity and adaptive immunity. Understanding how these programmes function in the context of cancer can provide valuable insights into host immune response, cancer severity and potential therapy response. Here we present a method that uses the transcriptomes of over 200 murine haematopoietic cells, to infer the lineage-specific haematopoietic activity present in human breast tumours. Correlating this activity with patient survival and tumour purity reveals that the transcriptional programmes of many cell types influence patient prognosis and are found in environments of high lymphocytic infiltration. Collectively, these results allow for a detailed and personalized assessment of the patient immune response to a tumour. When combined with routinely collected patient biopsy genomic data, this method can enable a richer understanding of the complex interplay between the host immune system and cancer.

The Discovery of Novel Biomarkers Improves Breast Cancer Intrinsic Subtype Prediction and Reconciles the Labels in the METABRIC Data Set

BACKGROUND

The prediction of breast cancer intrinsic subtypes has been introduced as a valuable strategy to determine patient diagnosis and prognosis, and therapy response. The PAM50 method, based on the expression levels of 50 genes, uses a single sample predictor model to assign subtype labels to samples. Intrinsic errors reported within this assay demonstrate the challenge of identifying and understanding the breast cancer groups. In this study, we aim to: a) identify novel biomarkers for subtype individuation by exploring the competence of a newly proposed method named CM1 score, and b) apply an ensemble learning, as opposed to the use of a single classifier, for sample subtype assignment. The overarching objective is to improve class prediction.

METHODS AND FINDINGS

The microarray transcriptome data sets used in this study are: the METABRIC breast cancer data recorded for over 2000 patients, and the public integrated source from ROCK database with 1570 samples. We first computed the CM1 score to identify the probes with highly discriminative patterns of expression across samples of each intrinsic subtype. We further assessed the ability of 42 selected probes on assigning correct subtype labels using 24 different classifiers from the Weka software suite. For comparison, the same method was applied on the list of 50 genes from the PAM50 method.

CONCLUSIONS

The CM1 score portrayed 30 novel biomarkers for predicting breast cancer subtypes, with the confirmation of the role of 12 well-established genes. Intrinsic subtypes assigned using the CM1 list and the ensemble of classifiers are more consistent and homogeneous than the original PAM50 labels. The new subtypes show accurate distributions of current clinical markers ER, PR and HER2, and survival curves in the METABRIC and ROCK data sets. Remarkably, the paradoxical attribution of the original labels reinforces the limitations of employing a single sample classifiers to predict breast cancer intrinsic subtypes.

IDEA: Integrated Drug Expression Analysis-Integration of Gene Expression and Clinical Data for the Identification of Therapeutic Candidates

Cancer drug discovery is an involved process spanning efforts from several fields of study and typically requires years of research and development. However, the advent of high-throughput genomic technologies has allowed for the use of in silico, genomics-based methods to screen drug libraries and accelerate drug discovery. Here we present a novel approach to computationally identify drug candidates for the treatment of breast cancer. In particular, we developed a Drug Regulatory Score similarity metric to evaluate gene expression profile similarity, in the context of drug treatment, and incorporated time-to-event patient survival information to develop an integrated analysis pipeline: Integrated Drug Expression Analysis (IDEA). We were able to predict drug candidates that have been known and those that have not been known in the literature to exhibit anticancer effects. Overall, our method enables quick preclinical screening of drug candidates for breast cancer and other diseases by using the most important indicator of drug efficacy: survival.

Integrative analysis of survival-associated gene sets in breast cancer

BACKGROUND

Patient gene expression information has recently become a clinical feature used to evaluate breast cancer prognosis. The emergence of prognostic gene sets that take advantage of these data has led to a rich library of information that can be used to characterize the molecular nature of a patient's cancer. Identifying robust gene sets that are consistently predictive of a patient's clinical outcome has become one of the main challenges in the field.

METHODS

We inputted our previously established BASE algorithm with patient gene expression data and gene sets from MSigDB to develop the gene set activity score (GSAS), a metric that quantitatively assesses a gene set's activity level in a given patient. We utilized this metric, along with patient time-to-event data, to perform survival analyses to identify the gene sets that were significantly correlated with patient survival. We then performed cross-dataset analyses to identify robust prognostic gene sets and to classify patients by metastasis status. Additionally, we created a gene set network based on component gene overlap to explore the relationship between gene sets derived from MSigDB. We developed a novel gene set based on this network's topology and applied the GSAS metric to characterize its role in patient survival.

RESULTS

Using the GSAS metric, we identified 120 gene sets that were significantly associated with patient survival in all datasets tested. The gene overlap network analysis yielded a novel gene set enriched in genes shared by the robustly predictive gene sets. This gene set was highly correlated to patient survival when used alone. Most interestingly, removal of the genes in this gene set from the gene pool on MSigDB resulted in a large reduction in the number of predictive gene sets, suggesting a prominent role for these genes in breast cancer progression.

CONCLUSIONS

The GSAS metric provided a useful medium by which we systematically investigated how gene sets from MSigDB relate to breast cancer patient survival. We used this metric to identify predictive gene sets and to construct a novel gene set containing genes heavily involved in cancer progression.

E2F4 regulatory program predicts patient survival prognosis in breast cancer

INTRODUCTION

Genetic and molecular signatures have been incorporated into cancer prognosis prediction and treatment decisions with good success over the past decade. Clinically, these signatures are usually used in early-stage cancers to evaluate whether they require adjuvant therapy following surgical resection. A molecular signature that is prognostic across more clinical contexts would be a useful addition to current signatures.

METHODS

We defined a signature for the ubiquitous tissue factor, E2F4, based on its shared target genes in multiple tissues. These target genes were identified by chromatin immunoprecipitation sequencing (ChIP-seq) experiments using a probabilistic method. We then computationally calculated the regulatory activity score (RAS) of E2F4 in cancer tissues, and examined how E2F4 RAS correlates with patient survival.

RESULTS

Genes in our E2F4 signature were 21-fold more likely to be correlated with breast cancer patient survival time compared to randomly selected genes. Using eight independent breast cancer datasets containing over 1,900 unique samples, we stratified patients into low and high E2F4 RAS groups. E2F4 activity stratification was highly predictive of patient outcome, and our results remained robust even when controlling for many factors including patient age, tumor size, grade, estrogen receptor (ER) status, lymph node (LN) status, whether the patient received adjuvant therapy, and the patient's other prognostic indices such as Adjuvant! and the Nottingham Prognostic Index scores. Furthermore, the fractions of samples with positive E2F4 RAS vary in different intrinsic breast cancer subtypes, consistent with the different survival profiles of these subtypes.

CONCLUSIONS

We defined a prognostic signature, the E2F4 regulatory activity score, and showed it to be significantly predictive of patient outcome in breast cancer regardless of treatment status and the states of many other clinicopathological variables. It can be used in conjunction with other breast cancer classification methods such as Oncotype DX to improve clinical outcome prediction.

USP9X downregulation renders breast cancer cells resistant to tamoxifen

Tamoxifen is one of the most widely used endocrine agents for the treatment of estrogen receptor α (ERα)-positive breast cancer. Although effective in most patients, resistance to tamoxifen is a clinically significant problem and the mechanisms responsible remain elusive. To address this problem, we performed a large scale loss-of-function genetic screen in ZR-75-1 luminal breast cancer cells to identify candidate resistance genes. In this manner, we found that loss of function in the deubiquitinase USP9X prevented proliferation arrest by tamoxifen, but not by the ER downregulator fulvestrant. RNAi-mediated attenuation of USP9X was sufficient to stabilize ERα on chromatin in the presence of tamoxifen, causing a global tamoxifen-driven activation of ERα-responsive genes. Using a gene signature defined by their differential expression after USP9X attenuation in the presence of tamoxifen, we were able to define patients with ERα-positive breast cancer experiencing a poor outcome after adjuvant treatment with tamoxifen. The signature was specific in its lack of correlation with survival in patients with breast cancer who did not receive endocrine therapy. Overall, our findings identify a gene signature as a candidate biomarker of response to tamoxifen in breast cancer.

Integrated analysis of expression profiling data identifies three genes in correlation with poor prognosis of triple-negative breast cancer

Triple-negative breast cancer (TNBC) shows more aggressive clinical behavior and poorer outcome than non-triple-negative breast cancer (NTNBC), and cannot be treated either via endocrine therapy or by Trastuzumab. For TNBC, chemotherapy is currently the mainstay of systemic medical treatment, the lack of more efficient options of treatment has been a problem in breast cancer prevention. In this study, we aimed to find genes related to prognosis in TNBC by bioinformatic analysis and to provide therapeutic candidates for TNBC treatment. We compared the differences in gene expression levels between cancer patients and healthy individuals across five breast cancer microarray databases to generate a gene cohort specifically upregulated in the NTNBC subtype, whose expression levels are ≥2-fold higher in TNBC compared to NTNBC and healthy individuals. Another two databases with clinical information were applied for following Kaplan-Meier analysis, and high expression of BIRC5, CENPA and FAM64A in this cohort were found to be related to poor survival (OS, DMFS, DFS and RFS). This correlation was also seen in patients at early stages and grades. On the other hand, the outcome of patients with synchronous upregulation of these three genes was the worst, while those with synchronous low gene level was the best. In conclusion, BIRC5, CENPA and FAM64A are specifically upregulated in TNBC, and the high expression of these three genes is associated with poor breast cancer prognosis, suggesting their clinical implication as therapeutic targets in TNBC.

Effect of aromatase inhibition on functional gene modules in estrogen receptor-positive breast cancer and their relationship with antiproliferative response

PURPOSE

To investigate potential associations between gene modules representing key biologic processes and response to aromatase inhibitors (AI) in estrogen receptor-positive (ER(+)) breast cancer.

PATIENTS AND METHODS

Paired gene expression and Ki67 protein expression were available from 69 postmenopausal women with ER(+) early breast cancer, at baseline and 2 weeks post-anastrozole treatment, in the presurgical setting. Functional gene modules (n = 26) were retrieved from published studies and their module scores were computed before and after elimination of proliferation-associated genes (PAG). Ki67 and module scores were assessed at baseline and 2 weeks post-anastrozole. Unsupervised clustering was used to assess associations between modules and Ki67.

RESULTS

Proliferation-based modules were highly correlated with Ki67 expression both pretreatment and on-treatment. At baseline with and without PAGs, Ki67 expression was significantly inversely correlated with ERG, ESR1.2, SET, and PIK3CA modules. Modules measuring estrogen signaling strongly predicted antiproliferative response to therapy with and without PAGs. Baseline expression of insulin-like growth factor-1 (IGF-I) module predicted a poor change in Ki67-implicating genes within the module as involved in de novo resistance to AIs. High expression of Immune.2.STAT1 module pretreatment predicted poor antiproliferative response to therapy. A significant association between estrogen-regulated genes modules (ESR1, ESR1-2, SET, and ERG) was evident post AI.

CONCLUSIONS

Multiple processes and pathways are affected by AI treatment in ER(+) breast cancer. Modules closely associated with ESR1 expression were predictive of good antiproliferative response to AIs, but modules representing immune activity and IGF-I/MAPK were predictive of poor Ki67 response, supporting their therapeutic targeting in combination with AIs.

CCNA2 is a prognostic biomarker for ER+ breast cancer and tamoxifen resistance

Identification of effective prognostic biomarkers and targets are of crucial importance to the management of estrogen receptor positive (ER+) breast cancer. CCNA2 (also known as CyclinA2) belongs to the highly conserved cyclin family and is significantly overexpressed in various cancer types. In this study, we demonstrated that CCNA2 had significant predictive power in distant metastasis free survival, disease free survival, recurrence free survival and overall survival of ER+ breast cancer patients. We also found that CCNA2 was closely associated with tamoxifen resistance. In addition, gene set enrichment analysis (GSEA) revealed that its expression was positively associated with genes overexpressed in endocrine therapy resistant samples. Finally, though CCNA2-Drug interaction network, we demonstrated the interactions between CCNA2 and several available cancer drugs. Overall, we suggest that CCNA2 is a biomarker for the prognosis of ER+ breast cancer and monitoring of tamoxifen efficacy. It's also a promising target for developing new strategies to prevent or even reverse tamoxifen resistance. Moreover, CCNA2 expression may help monitoring tamoxifen efficacy and directing personalized therapies. Nevertheless, in vivo and in vitro experiments and multi-center randomized controlled clinical trials are still needed before its application in clinical settings.