FOXM1 predicts overall and disease specific survival in muscle-invasive urothelial carcinoma and presents a differential expression between bladder cancer subtypes

Cyclin A2 Expression as Predictive Biomarker in Muscle-Invasive Upper Tract Urothelial Carcinoma

INTRODUCTION

The aim was to evaluate the prognostic value of altered Cyclin A2 (CCNA2) gene expression in upper tract urothelial carcinoma (UTUC) and to assess its predictive potential as a prognostic factor for overall survival (OS) and disease-free survival.

METHODS

62 patients who underwent surgical treatment for UTUC were included. Gene expression of CCNA2, MKI67, and p53 was analyzed by quantitative reverse transcriptase polymerase chain reaction. Survival analyses were performed using the Kaplan-Meier method and the log-rank test. For Cox regression analyses, uni- and multivariable hazard ratios were calculated. Spearman correlation was used to analyze correlation of CCNA2 expression with MKI67 and p53.

RESULTS

The median age of the cohort was 73 years, and it consisted of 48 males (77.4%) and 14 females (22.6%). Patients with high CCNA2 expression levels showed longer OS (HR 0.33; 95% CI: 0.15-0.74; p = 0.0073). Multivariable Cox regression analyses identified CCNA2 overexpression (HR 0.37; 95% CI: 0.16-0.85; p = 0.0189) and grading G2 (vs. G3) (HR 0.39; 95% CI: 0.17-0.87; p = 0.0168) to be independent predictors for longer OS. CCNA2 expression correlated positively with MKI67 expression (Rho = 0.4376, p = 0.0005).

CONCLUSION

Low CCNA2 expression is significantly associated with worse OS. Thus, CCNA2 might serve as a potential biomarker in muscle-invasive UTUC and may be used to characterize a subset of patients having an unfavorable outcome and for future risk assessment scores.

Classification of bladder cancer cell lines according to regulon activity

Bladder cancer is one of the most frequent cancers and causes more than 150.000 deaths each year. During the last decade, several studies provided important aspects about genomic characterization, consensus subgroup definition, and transcriptional regulation of bladder cancer. Still, much more research needs to be done to characterize molecular signatures of this cancer in depth. At this point, the use of bladder cancer cell lines is quite useful for the identification and test of new signatures. In this study, we classified the bladder cancer cell lines according to the activities of regulons implicated in the regulation of primary bladder tumors. Our regulon gene expression-based classification revealed three groups, neuronal-basal (NB), luminal-papillary (LP), and basal-squamous (BS). These regulon gene expression-based classifications showed a quite good concordance with the consensus subgroups assigned by the primary bladder cancer classifier. Importantly, we identified FGFR1 regulon to be involved in the characterization of the NB group, where neuroendocrine signature genes were significantly upregulated, and further β-catenin was shown to have significantly higher nuclear localization. LP groups were mainly driven by the regulons ERBB2, FOXA1, GATA3, and PPARG, and they showed upregulation of the genes involved in epithelial differentiation and urogenital development, while the activity of EGFR, FOXM1, STAT3, and HIF1A was implicated for the regulation of BS group. Collectively, our results and classifications may serve as an important guide for the selection and use of bladder cancer cell lines for experimental strategies, which aim to manipulate regulons critical for bladder cancer development.

High IL-22RA1 gene expression is associated with poor outcome in muscle invasive bladder cancer

BACKGROUND

The cell surface interleukin 22 (IL-22) receptor complex is mainly expressed in epithelial and tissue cells like pancreatitis cells. Recent studies described that IL-22R was overexpressed in malignant diseases and was associated with a poor overall survival (OS). The role of IL-22RA1 gene expression in muscle invasive bladder cancer (MIBC) has not been investigated, yet.

OBJECTIVES

The aim of this study was to analyze the role of IL-22RA1 gene expression in patients with MIBC.

METHODS

In a cohort of 114 patients with MIBC who underwent radical cystectomy, IL-22RA1 gene expression was analyzed with qRT-PCR and correlated with clinical parameters. Furthermore, Kaplan-Meier and Cox regression analysis were performed. For validation, an in silico dataset (TCGA 2017, n=407) was reanalyzed.

RESULTS

IL-22RA1 gene expression was independent of clinicopathological parameters like age (P=0.2681), T stage (P=0.2130), nodal status (P=0.3238) and lymph vascular invasion (LVI, P=0.5860) in patients with MIBC. A high expression of IL-22RA1 was associated with a shorter OS (P=0.0040) and disease-specific survival (P=0.0385). Furthermore, a shorter disease-free survival (DFS) was also associated with a high expression of IL-22RA1 (P=0.0102). In the multivariable analysis, IL-22RA1 expression was an independent prognostic predictors regarding OS (P=0.0096, HR=0.48). In the TCGA cohort, IL-22RA1 expression was independent regarding to OS and DFS.

CONCLUSION

A high IL-22RA1 gene expression was associated with worse outcome. Furthermore, IL-22RA1 represented an independent predictor regarding OS in our cohort and therefore might be used for risk stratification in patients with MIBC.

RNA Expression of DNA Damage Response Genes in Muscle-Invasive Bladder Cancer: Influence on Outcome and Response to Adjuvant Cisplatin-Based Chemotherapy

BACKGROUND

Perioperative cisplatin-based chemotherapy (CBC) can improve the outcome of patients with muscle-invasive bladder cancer (MIBC), but it is still to be defined which patients benefit. Mutations in DNA damage response genes (DDRG) can predict the response to CBC. The value of DDRG expression as a marker of CBC treatment effect remains unclear.

MATERIAL AND METHODS

RNA expression of the nine key DDRG (BCL2, BRCA1, BRCA2, ERCC2, ERCC6, FOXM1, RAD50, RAD51, and RAD52) was assessed by qRT-PCR in a cohort of 61 MICB patients (median age 66 y, 48 males, 13 females) who underwent radical cystectomy in a tertiary care center. The results were validated in the The Cancer Genome Atlas (TCGA) cohort of MIBC (n = 383). Gene expression was correlated with disease-free survival (DFS) and overall survival (OS). Subgroup analyses were performed in patients who received adjuvant cisplatin-based chemotherapy (ACBC) (Mannheim n = 20 and TCGA n = 75).

RESULTS

Low expression of RAD52 was associated with low DFS in both the Mannheim and the TCGA cohorts (Mannheim: p = 0.039; TCGA: p = 0.017). This was especially apparent in subgroups treated with ACBC (Mannheim: p = 0.0059; TCGA: p = 0.012). Several other genes showed an influence on DFS in the Mannheim cohort (BRCA2, ERCC2, FOXM1) where low expression was associated with poor DFS (p < 0.05 for all). This finding was not fully supported by the data in the TCGA cohort, where high expression of FOXM1 and BRCA2 correlated with poor DFS.

CONCLUSION

Low expression of RAD52 correlated with decreased DFS in the Mannheim and the TCGA cohort. This effect was especially pronounced in the subset of patients who received ACBC, making it a promising indicator for response to ACBC on the level of gene expression.

Neuropilin-2 and Its Transcript Variants Correlate with Clinical Outcome in Bladder Cancer

Urothelial bladder cancer ranks among the 10 most frequently diagnosed cancers worldwide. In our previous study, the transmembrane protein neuropilin-2 (NRP2) emerged as a predictive marker in patients with bladder cancer. NRP2 consists of several splice variants; the most abundant of these, NRP2a and NRP2b, are reported to have different biological functions in lung cancer progression. For other cancer types, there are no published data on the role of these transcript variants in cancer progression and the clinical outcome. Here, we correlate NRP2 and its two most abundant transcript variants, NRP2A and NRP2B, with the clinical outcome using available genomic data with subsequent validation in our own cohort of patients with muscle-invasive bladder cancer. In addition to NRP2, NRP1 and the NRP ligands PDGFC and PDGFD were studied. Only NRP2A emerged as an independent prognostic marker for shorter cancer-specific survival in muscle-invasive bladder cancer in our cohort of 102 patients who underwent radical cystectomy between 2008 and 2014 with a median follow-up time of 82 months. Additionally, we demonstrate that high messenger expression of NRP2, NRP1, PDGFC and PDGFD associates with a more aggressive disease (i.e., a high T stage, positive lymph node status and reduced survival).

Immune Cell-Associated Protein Expression Helps to Predict Survival in Muscle-Invasive Urothelial Bladder Cancer Patients after Radical Cystectomy and Optional Adjuvant Chemotherapy

Simple Summary

Adjuvant chemotherapy following radical cystectomy is a common therapy for muscle invasive bladder cancer (MIBC) patients. No applicable biomarkers exist to predict which patients will benefit from chemotherapy. Three immune cell markers, the chemokine CC motif ligand 2 (CCL2), the pan macrophage marker CD68, and the M2 macrophage marker CD163, were examined using immunohistochemistry to determine their predictive value for chemotherapy responses in different nodal stage and tumor stage subgroups. The presence of tumor-infiltrating immune cells, characterized by the markers CD68, CD163, and CCL2, was associated with a superior prognosis, and chemotherapy may not add an advantage for prognosis. However, a depleted immune microenvironment, here represented as a reduction or loss of macrophages, helped to predict the benefit of chemotherapy in N1 + 2 stage patients. Altogether, it is meaningful to consider the abundance of immune cells, such as macrophages, to better predict the response to chemotherapy for bladder cancer (BCa) patients after radical treatment.

Abstract

Bladder cancer (BCa) is the tenth most commonly diagnosed malignant cancer worldwide. Although adjuvant chemotherapy following radical cystectomy is a common therapy for muscle invasive bladder cancer patients, no applicable biomarkers exist to predict which patients will benefit from chemotherapy. In this study, we examined three immune cell markers, the chemokine CC motif ligand 2 (CCL2), the pan macrophage marker cluster of differentiation 68 (CD68) and the M2 macrophage marker cluster of differentiation 163 (CD163), using immunohistochemistry to determine their predictive value for the chemotherapy response in different nodal stage (pN0 vs. pN1 + 2) and tumor stage subgroups (pT2 vs. pT3 + 4). The prognosis was studied in terms of the overall survival (OS), disease-specific survival (DSS), and recurrence-free-survival (RFS) in 168 muscle invasive BCa patients. Chemotherapy was associated with a poorer prognosis in patients with a higher expression of the immune markers CCL2 (RFS), CD68 (DSS and RFS), and CD163 (DSS and RFS) in the N0 group and with poorer survival in patients with a higher expression of the immune markers CCL2 (OS, DSS, and RFS), CD68 (OS, DSS, and RFS), and CD163 (OS, DSS, and RFS) in the pT2 group when compared with treatments without chemotherapy. In contrast, chemotherapy was associated with a better prognosis in patients with a low expression of the immune markers CCL2 (DSS and RFS), CD68 (OS, DSS, and RFS), and CD163 (OS) in the N1 + 2 group. In addition, chemotherapy was associated with improved survival in patients with a low expression of the immune marker CD68 (OS and DSS) and there was a trend for a better prognosis in patients with a low expression of CD163 (OS) in the pT3 + 4 group compared to patients not treated with chemotherapy. Interestingly, CD68 appeared to be the most applicable immune marker to stratify patients by the outcome of chemotherapy in the nodal stage and tumor stage groups. Overall, we suggest that, in addition to the clinical factors of tumor stage and nodal stage, it is also meaningful to consider the abundance of immune cells, such as macrophages, to better predict the response to chemotherapy for BCa patients after radical treatment.

Experimental and bioinformatics considerations in cancer application of single cell genomics

Single cell genomics offers an unprecedented resolution to interrogate genetic heterogeneity in a patient's tumour at the intercellular level. However, the DNA yield per cell is insufficient for today's sequencing library preparation protocols. This necessitates DNA amplification which is a key source of experimental noise. We provide an evaluation of two protocols using micro-fluidics based amplification for whole exome sequencing, which is an experimental scenario commonly used in single cell genomics. The results highlight their respective biases and relative strengths in identification of single nucleotide variations. Towards this end, we introduce a workflow SoVaTSiC, which allows for quality evaluation and somatic variant identification of single cell data. As proof of concept, the framework was applied to study a lung adenocarcinoma tumour. The analysis provides insights into tumour phylogeny by identifying key mutational events in lung adenocarcinoma evolution. The consequence of this inference is supported by the histology of the tumour and demonstrates usefulness of the approach.

ANLN and TLE2 in Muscle Invasive Bladder Cancer: A Functional and Clinical Evaluation Based on In Silico and In Vitro Data

Anilin actin binding protein (ANLN) and transducing-like enhancer protein 2 (TLE2) are associated with cancer patient survival and progression. The impact of their gene expression on progression-free survival (PFS) of patients with muscle invasive bladder cancer (MIBC) treated with radical cystectomy (RC) and subtype association has not yet been investigated. qRT-PCR was used to measure the transcript levels of ANLN and TLE2 in the Mannheim cohort, and validated in silico by The Cancer Genome Atlas (TCGA) cohort. Uni- and multivariate Cox regression analyses identified predictors for disease-specific survival (DSS) and overall survival (OS). In the Mannheim cohort, tumors with high ANLN expression were associated with lower OS and DSS, while high TLE2 expression was associated with a favorable OS. The TCGA cohort confirmed that high ANLN and low TLE2 expression was associated with shorter OS and disease-free survival (DFS). In both cohorts, multivariate analyses showed ANLN and TLE2 expression as independent outcome predictors. Furthermore, ANLN was more highly expressed in cell lines and patients with the basal subtype, while TLE2 expression was higher in cell lines and patients with the luminal subtype. ANLN and TLE2 are promising biomarkers for individualized bladder cancer therapy including cancer subclassification and informed MIBC prognosis.

[Molecular tumor board-urothelial cancer]

BACKGROUND

Molecular tumor boards (MTB) are becoming more common. There are several molecular alterations in urothelial cancer a molecular tumor board can potentially rely on.

OBJECTIVES

The aim is to specify molecular alterations and their correlations with different clinical endpoints and to highlight potential questions addressed to a MTB for urothelial cancer.

MATERIALS AND METHODS

Descriptive review of the literature based on PubMed.

RESULTS

The landscape of molecular alterations in urothelial cancer is heterogeneous. Thus, recent biomarker research has been focusing on biomarker panels and classifiers instead of single biomarkers. Recently, molecular subtypes of urothelial cancer have been identified and correlated with different clinical endpoints. Furthermore, circulating tumor cells and tumor DNA are under investigation as potential biomarkers. In addition to treatment response and prognosis, molecular markers are also needed to improve clinical staging prior to radical cystectomy or for proper patient selection for neoadjuvant chemotherapy. Erdafitinib is the first targeted therapy (fibroblast growth factor receptor [FGFR] alteration) in urothelial cancer that was recently approved (in the USA).

CONCLUSIONS

Due to the lack of external validation, none of the identified biomarkers is currently established in clinical routine. In addition, there is no single driver mutation in urothelial cancer that facilitates the development of biomarkers and targeted therapies.

Analysis of the prognostic relevance of sex-steroid hormonal receptor mRNA expression in muscle-invasive urothelial carcinoma of the urinary bladder

Muscle-invasive urothelial carcinoma of the urinary bladder (UCB) often recurs following radical cystectomy (RC). An altered expression of sex-steroid hormone receptors has been associated with oncological outcomes of UCB and may represent therapeutic targets. Here the expression of different hormone receptors was measured on mRNA levels in patients treated by RC and associated with outcomes. Androgen receptor (AR), estrogen receptor 1 (ESR1), and progesterone receptor (PGR) mRNA expression was assessed by quantitative reverse transcription polymerase chain reaction (RT-qPCR) in RC samples of 87 patients with a median age of 66 (39-88) years. Univariate and multivariate analyses were performed to test associations with pathological and clinical characteristics as well as recurrence-free (RFS) and disease-specific survival (DSS). AR mRNA expression was lower in comparison with ESR1 and PGR expression (p < 0.0001). In univariate analysis, high expression levels of AR were associated with reduced RFS (HR 2.8, p = 0.015) and DSS (HR 2.8, p = 0.010). High AR mRNA expression and a positive lymph node status were independent predictors for reduced RFS (HR 2.5, p = 0.0049) and DSS (HR 3.4, p = 0.009). In patients with low AR mRNA expression, an increased ESR1 and PGR mRNA expression were associated with reduced RFS and DSS. High expression levels of AR are significantly associated with adverse outcome in patients with muscle-invasive UCB following RC. ESR1 and PGR expression status can further stratify patients with low AR expression into subgroups with significantly reduced RFS and DSS. Therapeutic targeting of AR may influence outcomes in patients with UCB.

Emerging Biomarkers in Bladder Cancer Identified by Network Analysis of Transcriptomic Data

Bladder cancer is a very common malignancy. Although new treatment strategies have been developed, the identification of new therapeutic targets and reliable diagnostic/prognostic biomarkers for bladder cancer remains a priority. Generally, they are found among differentially expressed genes between patients and healthy subjects or among patients with different tumor stages. However, the classical approach includes processing these data taking into consideration only the expression of each single gene regardless of the expression of other genes. These complex gene interaction networks can be revealed by a recently developed systems biology approach called Weighted Gene Co-expression Network Analysis (WGCNA). It takes into account the expression of all genes assessed in an experiment in order to reveal the clusters of co-expressed genes (modules) that, very probably, are also co-regulated. If some genes are co-expressed in controls but not in pathological samples, it can be hypothesized that a regulatory mechanism was altered and that it could be the cause or the effect of the disease. Therefore, genes within these modules could play a role in cancer and thus be considered as potential therapeutic targets or diagnostic/prognostic biomarkers. Here, we have reviewed all the studies where WGCNA has been applied to gene expression data from bladder cancer patients. We have shown the importance of this new approach in identifying candidate biomarkers and therapeutic targets. They include both genes and miRNAs and some of them have already been identified in the literature to have a role in bladder cancer initiation, progression, metastasis, and patient survival.

FOXM1 overexpression is associated with adverse outcome and predicts response to intravesical instillation therapy in stage pT1 non-muscle-invasive bladder cancer

OBJECTIVE

To investigate the role of forkhead box protein M1 (FOXM1) mRNA expression and its prognostic value in stage pT1 non-muscle-invasive bladder cancer (NMIBC).

PATIENTS AND METHODS

Clinical data and formalin-fixed paraffin-embedded tissues from transurethral resection of the bladder from patients with stage pT1 NMIBC, treated with an organ-preserving approach, were analysed retrospectively. Total RNA was isolated using commercial RNA extraction kits, and mRNA expression of FOXM1, MKI67, KRT20 and KRT5 was measured by single-step quantitative RT-PCR using RNA-specific TaqMan Assays. Statistical analysis was performed using Spearman's Rho, Wilcoxon or Kruskal-Wallis tests, Kaplan-Meier estimates of recurrence-free (RFS), progression-free (PFS) and cancer-specific survival (CSS) and Cox regression analysis.

RESULTS

Data from 296 patients (79.4% men, median age 72 years) were available for the final evaluation. Spearman correlation analysis showed that mRNA expression of FOXM1 was significantly correlated with MKI67 (ρ: 0.6530, P < 0.001) and with the luminal subtype, reflected by the positive correlation with KRT20 (ρ: 0.2113, P < 0.001). Furthermore, there was also a strong correlation of FOXM1 expression with adverse clinical and pathological variables, such as concomitant carcinoma in situ (P = 0.05), multifocal tumours (P = 0.005) and World Health Organization 1973 grade 3 disease (P < 0.001). Kaplan-Meier analysis showed overexpression of FOMX1 to be associated with worse PFS (P = 0.028) and worse CSS (P = 0.015). FOXM1 overexpression was also shown to be a predictive risk factor for CSS (hazard ratio 1.61 [1.13-2.34], L-R chi-squared: 7.19, P = 0.007). FOXM1 overexpression identified a subgroup of patients within the luminal subtype with worse RFS (P = 0.017), PFS (P < 0.001) and CSS (P = 0.015). Patients with low FOXM1 expression had better outcomes, irrespective of instillation therapy, whereas patients with high FOXM1 expression benefitted from intravesical chemotherapy with mitomycin C.

CONCLUSION

High FOXM1 expression was associated with adverse clinical and pathological features and worse outcomes, and predicted response to intravesical instillation therapy in patients with stage pT1 NMIBC.

FOXM1 predicts disease progression in non-muscle invasive bladder cancer

Purpose

The proto-oncogene forkhead box M1 (FOXM1) is associated with poor survival in many cancers. The impact of FOXM1 expression on progression-free survival (PFS) of non-muscle invasive bladder cancer (NMIBC) has not yet been investigated. The differential expression of FOXM1 between the different molecular NMIBC subtypes has further been assessed.

Methods

Transcript levels of FOXM1 and MKI67 were determined in 460 NMIBC patients (UROMOL cohort) by RNA-Seq and validated in silico by the Chungbuk and Lund cohort (n = 277). FOXM1 and MKI67 cutoffs were identified by the minimal p value method. Variables were evaluated by multivariable Cox regression analyses in order to identify independent predictors.

Results

FOXM1 is an independent predictor for PFS superior to current histological, clinical and molecular staging methods. Patients with high FOXM1 expression have a 6- to 8-fold higher risk of progression in multivariable analysis (p < 0.03). Highest transcript levels were found in the Class 2 and genomically unstable molecular NMIBC subtype (p < 0.03). The proto-oncogene further positively correlated with tumor grade and stage. NMIBCs with high FOXM1 expression showed a PFS advantage when treated with intravesical BCG instillation.

Conclusion

FOXM1 is a highly prognostic marker for disease progression of NMIBC superior to current histological, clinical and molecular staging methods and MKI67. It is mainly expressed in the Class 2 and genomically unstable molecular bladder cancer subtypes. Its role in drug resistance development makes FOXM1 valuable biomarker for NMIBC risk stratification.

Functional transcriptomic annotation and protein-protein interaction network analysis identify NEK2, BIRC5, and TOP2A as potential targets in obese patients with luminal A breast cancer

PURPOSE

Although obesity is a risk factor for breast cancer, little effort has been made in the identification of druggable molecular alterations in obese-breast cancer patients. Tumors are controlled by their surrounding microenvironment, in which the adipose tissue is a main component. In this work, we intended to describe molecular alterations at a transcriptomic and protein-protein interaction (PPI) level between obese and non-obese patients.

METHODS AND RESULTS

Gene expression data of 269 primary breast tumors were compared between normal-weight (BMI < 25, n = 130) and obese (IMC > 30, n = 139) patients. No significant differences were found for the global breast cancer population. However, within the luminal A subtype, upregulation of 81 genes was observed in the obese group (FC ≥ 1.4). Next, we explored the association of these genes with patient outcome, observing that 39 were linked with detrimental outcome. Their PPI map formed highly compact cluster and functional annotation analyses showed that cell cycle, cell proliferation, cell differentiation, and cellular response to extracellular stimuli were the more altered functions. Combined analyses of genes within the described functions are correlated with poor outcome. PPI network analyses for each function were to search for druggable opportunities. We identified 16 potentially druggable candidates. Among them, NEK2, BIRC5, and TOP2A were also found to be amplified in breast cancer, suggesting that they could act as strategic players in the obese-deregulated transcriptome.

CONCLUSION

In summary, our in silico analysis describes molecular alterations of luminal A tumors and proposes a druggable PPI network in obese patients with potential for translation to the clinical practice.

Downregulation of FoxM1 inhibits cell growth and migration and invasion in bladder cancer cells

The FoxM1 (Forkhead Box M1) transcription factor plays a key role in regulation of cell growth, cell cycle, and transformation. Higher expression of FoxM1 has been observed in various types of human cancers including bladder cancer. However, the exact function of FoxM1 in bladder cancer has not been elucidated. To investigate the cellular and molecular function of FoxM1 in bladder cancer, we measured the consequences of downregulation and upregulation of FoxM1 in bladder cancer cells using MTT assay, wound healing assay, and invasion assay. We found that downregulation of FoxM1 inhibited cell growth, but induced apoptosis in bladder cancer cells. Moreover, we found that inhibition of FoxM1 retarded cell migration and invasion. In line with this, upregulation of FoxM1 led to cell growth promotion and inhibited cell apoptosis in bladder cancer cells. Consistently, upregulation of FoxM1 led to increased cell migration and invasion. Our Western blotting results identified that downregulation of FoxM1 increased p27 level and inhibited VEGF, while overexpression of FoxM1 reduced p27 level and increased VEGF. Our findings suggest that FoxM1 could be a useful target for the treatment of bladder cancer.

FoxM1 is a promising candidate target in the treatment of breast cancer

Forkhead box protein M1(FoxM1) is a member of forkhead superfamily transcription factors. Emerging evidences have progressively contributed to our understanding on a central role of FoxM1 in human cancers. However, perspectives on the function of FoxM1 in breast cancer (BC) remain conflicting, and mostly were from basic research. Here, we explored the expression profile and prognostic values of FoxM1 based on analysis of pooled clinical datasets derived from online accessible databases, including ONCOMINE, Breast Cancer Gene-Expression Miner v4.0, and Kaplan-Meier plotter. It was found that, FoxM1 mRNA expression was significantly higher in breast tumor versus normal control. FoxM1expression profile presented a distinct pattern in different molecular subtypes of BC patients. Higher expression of FoxM1 was correlated to low mRNA expression of estrogen receptor 1 (ESR1), erb-B2 receptor tyrosine kinase 2 (ERBB2), and was inversely associated with the expression of classical luminal regulators forkhead box protein A1 (FoxA1) and GATA binding protein 3 (GATA3). Elevated FoxM1 expression predicted shorter distance metastasis free survival (DMFS) in BC patients, particularly with estrogen receptor (ER) positive and Luminal A, Luminal B subtypes of BC. More interestingly, elevated FoxM1 expression predicted poor survival in breast cancer patients, especially in the ER (+), progesterone receptor (PR) (+) subgroups and BC patients received adjuvant chemotherapy only or treated with tamoxifen only. These results implied that FoxM1 is an essential prognostic factor and promising candidate target in the treatment of breast cancer.