Survival analysis in breast cancer: evaluating ensemble learning techniques for prediction.

Breast cancer is most commonly faced with form of cancer amongst women worldwide. In spite of the fact that the breast cancer research and awareness have gained considerable momentum, there is still no one treatment due to disease heterogeneity. Survival data may be of specific interest in breast cancer studies to understand its dynamic and complex trajectories. This study copes with the most important covariates affecting the disease progression. The study utilizes the German Breast Cancer Study Group 2 (GBSG2) and the Molecular Taxonomy of Breast Cancer International Consortium dataset (METABRIC) datasets. In both datasets, interests lie in relapse of the disease and the time when the relapse happens. The three models, namely the Cox proportional hazards (PH) model, random survival forest (RSF) and conditional inference forest (Cforest) were employed to analyse the breast cancer datasets. The goal of this study is to apply these methods in prediction of breast cancer progression and compare their performances based on two different estimation methods: the bootstrap estimation and the bootstrap .632 estimation. The model performance was evaluated in concordance index (C-index) and prediction error curves (pec) for discrimination. The Cox PH model has a lower C-index and bigger prediction error compared to the RSF and the Cforest approach for both datasets. The analysis results of GBSG2 and METABRIC datasets reveal that the RSF and the Cforest algorithms provide non-parametric alternatives to Cox PH model for estimation of the survival probability of breast cancer patients.

Classification of HER2-negative breast cancers by ERBB2 copy number alteration status reveals molecular differences associated with chromosome 17 gene aberrations.

Background: Recently, HER2-negative breast cancers have been reclassified by protein expression into 'HER2-low' and 'HER2-zero' subgroups, but the consideration of HER2-low breast cancer as a distinct biological subtype with differing prognoses remains controversial. By contrast, non-neutral ERBB2 copy number alteration (CNA) status is associated with inferior survival outcomes compared to ERBB2 CNA-neutral breast cancer, providing an alternative approach to classification. Methods: Here, we investigated the molecular landscape of non-metastatic HER2-negative BCs in relation to ERBB2 CNA status to elucidate biological differences. Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) and The Cancer Genome Atlas (TCGA) TCGA-BRCA datasets (n = 1875) were analyzed. Results: Nearly two-fifths of the cohort harbored ERBB2 CNAs (39.4%), which were significantly enriched within hormone receptor-negative (56.1%) than within hormone receptor-positive BCs (35.5%; p < 0.0001). Globally, CNAs across the genome were significantly higher in ERBB2 non-neutral compared to neutral cohorts (p < 0.0001). Notably, genetic aberrations on chromosome 17 - BRCA1, NF1, TP53, MAP2K4, and NCOR1 - were widespread in the ERBB2 non-neutral cases. While chromosome 17q arm-level alterations were largely in tandem with ERBB2 CNA status, arm-level loss in chromosome 17p was prevalent regardless of ERBB2 gain, amplification, or loss. Differential gene expression analysis demonstrated that pathways involved in the cell cycle, proteasome, and DNA replication were upregulated in ERBB2 non-neutral cases. Conclusion: Classification of HER2-negative BCs according to ERBB2 CNA status reveals differences in the genomic landscape. The implications of concurrent aberrations in other genes on chromosome 17 merit further research in ERBB2 non-neutral BCs.

Glutaminase (GLS1) gene expression in primary breast cancer.

BACKGROUND: Tumor growth is mediated in part by glutamine, and glutaminase is an enzyme necessary for glutamine catabolism. We studied glutaminase (GLS1) gene expression in primary breast cancer to determine correlations with clinical and tumor characteristics, and gene associations in publicly available databases. A better understanding of glutaminase gene expression may help guide further exploration of glutaminase inhibitors in breast cancer. METHODS: GLS1 mRNA levels were evaluated in The Cancer Genome Atlas (n = 817) and METABRIC (n = 1992) datasets. Associations between GLS1 and tumor subtype (ANOVA followed by post-hoc Tukey test for pairwise comparisons) and selected genes involved in the pathogenesis of breast cancer (Pearson's correlations) were determined in both datasets. In METABRIC, associations with overall survival (Cox proportional hazard model) were determined. For all analyses, p < 0.05 was the threshold for statistical significance. RESULTS: GLS1 expression was significantly higher in triple negative breast cancer (TNBC) than hormone receptor (HR) +/HER2- and HER2+ breast cancer (p < 0.001) and basal versus luminal A, luminal B, and HER2 enriched breast cancer (p < 0.001) in both datasets. In METABRIC, higher GLS1 expression was associated with improved overall survival (HR 0.91, 95% CI: 0.85-0.97, p = 0.005) and this association remained significant in the TNBC subset (HR 0.83, 95% CI: 0.71-0.98, p = 0.032). GLS1 had significant positive gene correlations with immune, proliferative, and basal genes, and inverse correlations with luminal genes and genes involved in metabolism. CONCLUSION: GLS1 expression is highest in TNBC and basal breast cancer, supporting ongoing clinical investigation of GLS1 inhibition in TNBC. GLS1 may have prognostic implications but further research is needed to validate this finding. GLS1 had significant positive gene correlations with immune genes, which may have implications for potential combinations of glutaminase inhibition and immunotherapy.

Establishing the Link between X-Chromosome Aberrations and TP53 Status, with Breast Cancer Patient Outcomes.

Ubiquitous to normal female human somatic cells, X-chromosome inactivation (XCI) tightly regulates the transcriptional silencing of a single X chromosome from each pair. Some genes escape XCI, including crucial tumour suppressors. Cancer susceptibility can be influenced by the variability in the genes that escape XCI. The mechanisms of XCI dysregulation remain poorly understood in complex diseases, including cancer. Using publicly available breast cancer next-generation sequencing data, we show that the status of the major tumour suppressor TP53 from Chromosome 17 is highly associated with the genomic integrity of the inactive X (Xi) and the active X (Xa) chromosomes. Our quantification of XCI and XCI escape demonstrates that aberrant XCI is linked to poor survival. We derived prognostic gene expression signatures associated with either large deletions of Xi; large amplifications of Xa; or abnormal X-methylation. Our findings expose a novel insight into female cancer risks, beyond those associated with the standard molecular subtypes.

NRF2 signaling pathway: A comprehensive prognostic and gene expression profile analysis in breast cancer.

Breast cancer is the most frequently diagnosed malignant tumor in women and a major public health concern. NRF2 axis is a cellular protector signaling pathway protecting both normal and cancer cells from oxidative damage. NRF2 is a transcription factor that binds to the gene promoters containing antioxidant response element-like sequences. In this report, differential expression of NRF2 signaling pathway elements, as well as the correlation of NRF2 pathway mRNAs with various clinicopathologic characteristics, including molecular subtypes, tumor grade, tumor stage, and methylation status, has been investigated in breast cancer using METABRIC and TCGA datasets. In the current report, our findings revealed the deregulation of several NRF2 signaling elements in breast cancer patients. Moreover, there were negative correlations between the methylation of NRF2 genes and mRNA expression. The expression of NRF2 genes significantly varied between different breast cancer subtypes. In conclusion, substantial deregulation of NRF2 signaling components suggests an important role of these genes in breast cancer. Because of the clear associations between mRNA expression and methylation status, DNA methylation could be one of the mechanisms that regulate the NRF2 pathway in breast cancer. Differential expression of Hippo genes among various breast cancer molecular subtypes suggests that NRF2 signaling may function differently in different subtypes of breast cancer. Our data also highlights an interesting link between NRF2 components' transcription and tumor grade/stage in breast cancer.

Loss of genes in chromosome arms 5q and 16q in breast cancer.

INTRODUCTION: Losses of genetic material from chromosomes 5q and 16q commonly occur in sub-sets of breast cancer. Their significance from a pathophysiologic point of view is not well-defined. METHODS: This study uses publicly available genomic data from extensive breast cancer datasets to define the landscape of losses in chromosomal arms 5q and 16q in the two sub-types of breast cancer they most commonly occur, basal-like and luminal A cancers, respectively. RESULTS: It is shown that dozens of genes from these chromosomal arms are putatively hemi-deleted in few samples each. No individual gene from either 5q or 16q shows an incidence of deep deletion above 10% in the cohorts with basal-like and luminal A cancers or in the whole cohorts. A few tumor suppressor genes are deleted in a small number of samples, less than 5% in each cohort. Losses of 5q or 16q confer no survival advantage in either the basal-like or the luminal A cohorts from TCGA that harbor them. CONCLUSION: Results suggest that there are no individual genes in chromosomes 5q and 16q whose loss can be implicated in a dominant pathophysiologic sequence of events in breast cancer or its sub-sets.

Hippo signaling pathway: A comprehensive gene expression profile analysis in breast cancer.

Breast cancer (BC) is the most frequently diagnosed malignancy in women and a major public health concern. The Hippo pathway is an evolutionarily conserved signaling pathway that serves as a key regulator for a wide variety of biological processes. Hippo signaling has been shown to have both oncogenic and tumor-suppressive functions in various cancers. Core components of the Hippo pathway consist of various kinases and downstream effectors such as YAP/TAZ. In the current report, differential expression of Hippo pathway elements as well as the correlation of Hippo pathway mRNAs with various clinicopathologic characteristics, including molecular subtypes, receptor status, and methylation status, has been investigated in BC using METABRIC and TCGA datasets. In this review, we note deregulation of several Hippo signaling elements in BC patients. Moreover, we see examples of negative correlations between methylation of Hippo genes and mRNA expression. The expression of Hippo genes significantly varies between different receptor subgroups. Because of the clear associations between mRNA expression and methylation status, DNA methylation may be one of the mechanisms that regulate the Hippo pathway in BC cells. Differential expression of Hippo genes among various BC molecular subtypes suggests that Hippo signaling may function differently in different subtypes of BC. Our data also highlights an interesting link between Hippo components' transcription and ER negativity in BC. In conclusion, substantial deregulation of Hippo signaling components suggests an important role of these genes in breast cancer.

NR4A Family Genes: A Review of Comprehensive Prognostic and Gene Expression Profile Analysis in Breast Cancer.

This report analyzes nuclear receptor (NR) subfamily 4A's potential role in treating those diagnosed with breast cancer. Here we reviewed the current literature on NR4 family members. We also examined the relative gene expression of the NR4A receptor subfamily in the basal, HER2 (human epidermal growth factor receptor 2) positive, luminal A, and luminal B subtypes using data from tumor samples in The Cancer Genome Atlas (TCGA) and Molecular Taxonomy of Breast Cancer International Consortium (METABRIC). These data showed a positive link between NR4A1-NR4A3 expression and increased overall survival and relapse-free survival in breast cancer patients. In addition, we observed that high expression of NR4A1, NR4A2, and NR4A3 led to better survival. Furthermore, NR4A family genes seem to play an essential regulatory role in glycolysis and oxidative phosphorylation in breast cancer. The novel prognostic role of the NR4A1-NR4A3 receptors implicates these receptors as important mediators controlling breast cancer metabolic reprograming and its progression. The review establishes a strong clinical basis for the investigation of the cellular, molecular, and physiological roles of NR4A genes in breast cancer.

HER2 expression, copy number variation and survival outcomes in HER2-low non-metastatic breast cancer: an international multicentre cohort study and TCGA-METABRIC analysis.

BACKGROUND: HER2-low breast cancer (BC) is currently an area of active interest. This study evaluated the impact of low expression of HER2 on survival outcomes in HER2-negative non-metastatic breast cancer (BC). METHODS: Patients with HER2-negative non-metastatic BC from 6 centres within the Asian Breast Cancer Cooperative Group (ABCCG) (n = 28,280) were analysed. HER2-low was defined as immunohistochemistry (IHC) 1+ or 2+ and in situ hybridization non-amplified (ISH-) and HER2-zero as IHC 0. Relapse-free survival (RFS) and overall survival (OS) by hormone receptor status and HER2 IHC 0, 1+ and 2+ ISH- status were the main outcomes. A combined TCGA-BRCA and METABRIC cohort (n = 1967) was also analysed to explore the association between HER2 expression, ERBB2 copy number variation (CNV) status and RFS. RESULTS: ABCCG cohort median follow-up was 6.6 years; there were 12,260 (43.4%) HER2-low BC and 16,020 (56.6%) HER2-zero BC. The outcomes were better in HER2-low BC than in HER2-zero BC (RFS: centre-adjusted hazard ratio (HR) 0.88, 95% CI 0.82-0.93, P < 0.001; OS: centre-adjusted HR 0.82, 95% CI 0.76-0.89, P < 0.001). On multivariable analysis, HER2-low status was prognostic (RFS: HR 0.90, 95% CI 0.85-0.96, P = 0.002; OS: HR 0.86, 95% CI 0.79-0.93, P < 0.001). These differences remained significant in hormone receptor-positive tumours and for OS in hormone receptor-negative tumours. Superior outcomes were observed for HER2 IHC1+ BC versus HER2-zero BC (RFS: HR 0.89, 95% CI 0.83-0.96, P = 0.001; OS: HR 0.85, 95% CI 0.78-0.93, P = 0.001). No significant differences were seen between HER2 IHC2+ ISH- and HER2-zero BCs. In the TCGA-BRCA and METABRIC cohorts, ERBB2 CNV status was an independent RFS prognostic factor (neutral versus non-neutral HR 0.71, 95% CI 0.59-0.86, P < 0.001); no differences in RFS by ERBB2 mRNA expression levels were found. CONCLUSIONS: HER2-low BC had a superior prognosis compared to HER2-zero BC in the non-metastatic setting, though absolute differences were modest and driven by HER2 IHC 1+ BC. ERBB2 CNV merits further investigation in HER2-negative BC.

Comprehensive Analysis of the Implication of PGRMC1 in Triple-Negative Breast Cancer.

TNBC represents the most malignant subtype of breast cancer with heterogenicity and poor prognosis. PGRMC1 has been reported to predict worse prognosis and correlate with MHT mediated signal transduction in breast cancer, whereas its involvement in TNBC remains poorly explored. The purpose of the study was to explore the roles of PGRMC1 in TNBC. Bioinformatic approaches were performed to analyzed the expression of PGRMC1 among different subtypes of breast cancers using RNA-seq data from the TCGA, METABRIC and GEO databases. PGRMC1 mRNA expression and survival in breast cancer were analyzed. Furthermore, we analyzed the expression of PGRMC1 in TNBC by single cell RNA-seq data and immunohistochemistry. The expression of PGRMC1 in TNBC group was significantly higher compared with that of Luminal subtypes, especially in the epithelia cells, which was further proved by IHC at protein level. Better overall survival (p = 0.027) was observed in the patients with lower expression of PGRMC1. Different states of hormone and Her2 receptors contributed to the distinct functions of PGRMC1. In TNBC, PGRMC1 might play an important role in mitochondrial functions. In summary, this study revealed the correlation between PGRMC1 expression and its clinical significance in TNBC, probably through mitochondria-associated pathway, which may provide new ideas for prognosis and therapy of TNBC.

Low expression of miR-195 is associated with cell proliferation, glycolysis and poor survival in estrogen receptor (ER)-positive but not in triple negative breast cancer.

MiR-195 is a tumor suppressive microRNA in breast cancer. Its clinical relevance remains debatable as it has only been studied via in vitro experiments or small cohort studies. We analyzed a total of 2,038 patients in the TCGA and METABRIC cohorts to assess whether low miR-195 expressing tumors are associated with aggressive cancer characteristics and poor prognostic outcomes. The median cutoff of miR-195 expression was used to split the groups into miR-195 high and low groups. Low miR-19 expressing tumors demonstrated high cell proliferating features by enriching the gene sets associated with cell proliferation, MKI67 expression and pathological grade. One-third of the top target miR-195 genes were related to cell proliferation. Low miR-195 expressing tumors were associated with both pro-cancerous and anti-cancerous immune cells. Low miR-195 expressing tumors were associated with enhanced glycolysis and poor survival in ER-positive tumors, but not other subtypes of breast cancer. In conclusion, low expression of miR-195 in ER-positive breast cancer was associated with enhanced cancer cell proliferation, glycolysis, and worse overall survival.

Enhanced Thermogenesis in Triple-Negative Breast Cancer Is Associated with Pro-Tumor Immune Microenvironment.

Mild cold stress induced by housing mice with a 4T1 triple-negative breast cancer (TNBC) cell implantation model at 22 °C increases tumor growth rate with a pro-tumorigenic immune microenvironment (lower CD8 +T cells, higher myeloid-derived suppressor cells (MDSCs) and regulatory T-cells (Tregs)). Since cold stress also activates thermogenesis, we hypothesized that enhanced thermogenesis is associated with more aggressive cancer biology and unfavorable tumor microenvironment (TME) in TNBC patients. A total of 6479 breast cancer patients from METABRIC, TCGA, GSE96058, GSE20194, and GSE25066 cohorts were analyzed using Kyoto Encyclopedia of Genes and Genomes (KEGG) thermogenesis score. High-thermogenesis TNBC was associated with a trend towards worse survival and with angiogenesis, adipogenesis, and fatty acid metabolism pathways. On the other hand, low-thermogenesis TNBC enriched most of the hallmark cell-proliferation-related gene sets (i.e., mitotic spindle, E2F targets, G2M checkpoint, MYC targets), as well as immune-related gene sets (i.e., IFN-α and IFN-γ response). Favorable cytotoxic T-cell-attracting chemokines CCL5, CXCL9, CXCL10, and CXCL11 were lower; while the MDSC- and Treg-attracting chemokine CXCL12 was higher. There were higher M2 but lower M1 macrophages and Tregs. In conclusion, high-thermogenesis TNBC is associated with pro-tumor immune microenvironment and may serve as biomarker for testing strategies to overcome this immunosuppression.

MIEAP, a p53-downstream gene, is associated with suppression of breast cancer cell proliferation and better survival.

Mitochondria-eating protein (MIEAP; also known as SPATA18), a p53-downstream gene, is involved in mitochondrial quality control (MQC). Enforced MIEAP expression induces caspase-dependent cell death in vitro, and impairment of the p53/MIEAP-regulated MQC pathway is frequently observed in breast cancer (BC), resulting in poor disease-free survival (DFS). To investigate the clinical significance of MIEAP in BC, we identified 2,980 patients from two global, large-scale primary BC cohorts: the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC; n=1,904) and the Cancer Genome Atlas (TCGA; n=1,076). We divided patients in each cohort into high and low groups based on median gene expression levels and analyzed the association between MIEAP expression and clinical outcomes. Compared with normal tumors, MIEAP expression was significantly downregulated in all patients with p53-mutant BC regardless of subtype. MIEAP expression was negatively correlated with KI67 expression. Gene set enrichment analysis demonstrated that cell cycle- and proliferation-associated gene sets were significantly enriched in MIEAP-low tumors compared to MIEAP-high tumors. Patients with MIEAP-high luminal subtype were associated with significantly longer DFS than those with MIEAP-low luminal tumors in both cohorts, whereas significantly longer overall survival was observed only in the METABRIC cohort, which has roughly double the number of samples. These results indicated that the mechanistic role of MIEAP is clinically relevant in the two independent cohorts. This is the first study to use large cohorts to demonstrate the association between MIEAP expression and survival in patients with luminal subtype BC.

Intratumoral Adipocyte-High Breast Cancer Enrich for Metastatic and Inflammation-Related Pathways but Associated with Less Cancer Cell Proliferation.

Cancer-associated adipocytes are known to cause inflammation, leading to cancer progression and metastasis. The clinicopathological and transcriptomic data from 2256 patients with breast cancer were obtained based on three cohorts: The Cancer Genome Atlas (TCGA), GSE25066, and a study by Yau et al. For the current study, we defined the adipocyte, which is calculated by utilizing a computational algorithm, xCell, as "intratumoral adipocyte". These intratumoral adipocytes appropriately reflected mature adipocytes in a bulk tumor. The amount of intratumoral adipocytes demonstrated no relationship with survival. Intratumoral adipocyte-high tumors significantly enriched for metastasis and inflammation-related gene sets and are associated with a favorable tumor immune microenvironment, especially in the ER+/HER2- subtype. On the other hand, intratumoral adipocyte-low tumors significantly enriched for cell cycle and cell proliferation-related gene sets. Correspondingly, intratumoral adipocyte-low tumors are associated with advanced pathological grades and inversely correlated with MKI67 expression. In conclusion, a high amount of intratumoral adipocytes in breast cancer was associated with inflammation, metastatic pathways, cancer stemness, and favorable tumor immune microenvironment. However, a low amount of adipocytes was associated with a highly proliferative tumor in ER-positive breast cancer. This cancer biology may explain the reason why patient survival did not differ by the amount of adipocytes.

Contribution of Immune Cells to Glucocorticoid Receptor Expression in Breast Cancer.

Breast cancer (BC) patients experience increased stress with elevated cortisol levels, increasing risk of cancer recurrence. Cortisol binds to a cytoplasmic receptor, glucocorticoid receptor (GR) encoded by GR gene (NR3C1). We hypothesized that not only cancer cells, but even immune cells in the tumor microenvironment (TME) may contribute to GR expression in bulk tumor and influence prognosis. To test this, mRNA expression data was accessed from METABRIC and TCGA. "High" and "low" expression was based on highest and lowest quartiles of NR3C1 gene expression, respectively. Single-cell sequencing data were obtained from GSE75688 and GSE114725 cohorts. Computer algorithms CIBERSORT, Gene Set Enrichment Analysis and TIMER were used. GR-high BC has better median disease-free and disease-specific survival. Single cell sequencing data showed higher GR expression on immune cells compared to cancer and stromal cells. Positive correlation between GR-high BC and CD8+ T-cells was noted. In GR-high tumors, higher cytolytic activity (CYT) with decreased T-regulatory and T-follicular helper cells was observed. High GR expression was associated with lower proliferation index Ki67, enriched in IL-2_STAT5, apoptosis, KRAS, TGF-ß signaling, and epithelial-to-mesenchymal transition. Immune cells significantly contribute to GR expression of bulk BC. GR-high BC has a favorable TME with higher CYT with favorable outcomes.

Transcriptional profiling reveals a subset of human breast tumors that retain wt TP53 but display mutant p53-associated features.

TP53 gene mutations are very common in human cancer. While such mutations abrogate the tumor suppressive activities of the wild-type (wt) p53 protein, some of them also endow the mutant (mut) protein with oncogenic gain of function (GOF), facilitating cancer progression. Yet, p53 may acquire altered functionality even without being mutated; in particular, experiments with cultured cells revealed that wtp53 can be rewired to adopt mut-like features in response to growth factors or cancer-mimicking genetic manipulations. To assess whether such rewiring also occurs in human tumors, we interrogated gene expression profiles and pathway deregulation patterns in the METABRIC breast cancer (BC) dataset as a function of TP53 gene mutation status. Harnessing the power of machine learning, we optimized a gene expression classifier for ER+Her2- patients that distinguishes tumors carrying TP53 mutations from those retaining wt TP53. Interestingly, a small subset of wt TP53 tumors displayed gene expression and pathway deregulation patterns markedly similar to those of TP53-mutated tumors. Moreover, similar to TP53-mutated tumors, these 'pseudomutant' cases displayed a signature for enhanced proliferation and had worse prognosis than typical wtp53 tumors. Notably, these tumors revealed upregulation of genes which, in BC cell lines, were reported to be positively regulated by p53 GOF mutants. Thus, such tumors may benefit from mut p53-associated activities without having to accrue TP53 mutations.

High Expression of microRNA-143 is Associated with Favorable Tumor Immune Microenvironment and Better Survival in Estrogen Receptor Positive Breast Cancer.

microRNA-143(miR-143) is a well-known tumor suppressive microRNA that exhibits anti-tumoral function by targeting KRAS signaling pathways in various malignancies. We hypothesized that miR-143 suppresses breast cancer progression by targeting KRAS and its effector molecules. We further hypothesized that high expression of miR-143 is associated with a favorable tumor immune microenvironment of estrogen receptor (ER)-positive breast cancer patients which result in improved survival. Two major publicly available breast cancer cohorts; The Cancer Genome Atlas (TCGA) and Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) were used. The miR-143 high expression group was associated with increased infiltration of anti-cancer immune cells and decreased pro-cancer immune cells, as well as enrichment of the genes relating to T helper (Th1) cells resulting in improved overall survival (OS) in ER-positive breast cancer patients. To the best of our knowledge, this is the first study to demonstrate that high expression of miR-143 in cancer cells associates with a favorable tumor immune microenvironment, upregulation of anti-cancer immune cells, and suppression of the pro-cancer immune cells, associating with better survival of the breast cancer patients.

High Expression of miR-34a Associated with Less Aggressive Cancer Biology but Not with Survival in Breast Cancer.

Most breast cancer (BC) patients succumb to metastatic disease. MiR-34a is a well-known tumor suppressive microRNA which exerts its anti-cancer functions by playing a role in p53, apoptosis induction, and epithelial-mesenchymal transition (EMT) suppression. Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) and The Cancer Genome Atlas (TCGA) cohorts were used to test our hypothesis that miR-34a high BCs translate to less aggressive cancer biology and better survival in large cohorts. There was no association between miR-34a expression levels and clinicopathological features of BC patients except for HER2 positivity. MiR-34a high expressing tumors were associated with lower Nottingham pathological grades and lower MKI67 expression. In agreement, high miR-34a tumors demonstrated lower GSVA scores of cell cycle and cell proliferation-related gene sets. High miR-34a tumors enriched the p53 pathway and apoptosis gene sets. Unexpectedly, high miR-34a tumors also associated with elevated EMT pathway score and ZEB1 and two expressions. MiR-34a expression did not associate with any distant metastasis. Further, high miR-34a tumors did not associate with better survival compared with miR-34a low tumors. In conclusion, the clinical relevance of miR-34a high expressing tumors was associated with suppressed cell proliferation, enhanced p53 pathway and apoptosis, but enhanced EMT and these findings did not reflect better survival outcomes in large BC patient cohorts.

KRAS signaling enriched triple negative breast cancer is associated with favorable tumor immune microenvironment and better survival.

KRAS signaling is associated with cancer progression in several cancers. Upregulation of KRAS signaling is often seen in cancers that harbor high KRAS mutation rate, such as pancreatic cancer and non-small cell lung cancer (NSCLC). Less than 2% of breast cancers have KRAS mutation, however, the alteration of the effector signaling such as PI3K/AKT and MAPK pathways are well known. Mutated KRAS is known to function as immune suppressor in other cancers, but the role of KRAS signaling on tumor immune microenvironment (TIME) in breast cancer is not known. We hypothesize that the enrichment of KRAS signaling is associated with reduced patient survival as well as TIME in triple negative breast cancer (TNBC). Patient cohorts from Molecular Taxonomy of Breast Cancer International Consortium (METABRIC; n = 1903) and The Cancer Genome Atlas (TCGA; n = 982) were used. Higher expression of KRAS in breast cancer cell-lines (MCF7, BT474, and MDA-MB231) compared to MCF10A, which is a model of benign mammary cells was found. Both MEK and PI3K inhibitors suppressed MB231 cell proliferation in dose dependent manner. Gene Set Variant Analysis (GSVA) of the patient cohorts demonstrated two peaks by KRAS_SIGNALING_UP gene sets which were divided into KRAS-high and -low groups using median cutoff. There was no difference in KRAS mutation between KRAS-high and low. Despite its cell proliferation promoting role, KRAS-high patients demonstrated significantly better Disease-Free Survival and Overall Survival in triple negative breast cancer (TNBC). KRAS-high TNBC was associated with favorable tumor immune microenvironment with elevated B cells and CD8 T cells, monocytes, or M1 macrophage. It was associated with decreased CD4 central memory T-cells, but not Regulatory T-cells, or M2 macrophage detected by xCell. To elucidate the mechanism of this association, Gene Set Enrichment Analysis was performed. Inflammatory response, IL6/JAK-STAT3 signaling, and Interferon gamma response gene sets were enriched in KRAS-high TNBC patients in both METABRIC and TCGA cohorts. In agreement, cytolytic activity score, interferon gamma response score, and lymphocyte infiltrating signature score, were all significantly elevated in KRAS-high TNBC. In conclusion, we found that patients with enrichment of KRAS signaling gene sets were associated with inflammation and favorable tumor immune microenvironment as well as improved survival in TNBC.

The G Protein-Coupled Estrogen Receptor (GPER) Expression Correlates with Pro-Metastatic Pathways in ER-Negative Breast Cancer: A Bioinformatics Analysis.

The G protein-coupled estrogen receptor (GPER, formerly known as GPR30) is a seven-transmembrane receptor that mediates estrogen signals in both normal and malignant cells. In particular, GPER has been involved in the activation of diverse signaling pathways toward transcriptional and biological responses that characterize the progression of breast cancer (BC). In this context, a correlation between GPER expression and worse clinical-pathological features of BC has been suggested, although controversial data have also been reported. In order to better assess the biological significance of GPER in the aggressive estrogen receptor (ER)-negative BC, we performed a bioinformatics analysis using the information provided by The Invasive Breast Cancer Cohort of The Cancer Genome Atlas (TCGA) project and Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) datasets. Gene expression correlation and the statistical analysis were carried out with R studio base functions and the tidyverse package. Pathway enrichment analysis was evaluated with Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway on the Database for Annotation, Visualization and Integrated Discovery (DAVID) website, whereas gene set enrichment analysis (GSEA) was performed with the R package phenoTest. The survival analysis was determined with the R package survivALL. Analyzing the expression data of more than 2500 primary BC, we ascertained that GPER levels are associated with pro-migratory and metastatic genes belonging to cell adhesion molecules (CAMs), extracellular matrix (ECM)-receptor interaction, and focal adhesion (FA) signaling pathways. Thereafter, evaluating the disease-free interval (DFI) in ER-negative BC patients, we found that the subjects expressing high GPER levels exhibited a shorter DFI in respect to those exhibiting low GPER levels. Overall, our results may pave the way to further dissect the network triggered by GPER in the breast malignancies lacking ER toward a better assessment of its prognostic significance and the action elicited in mediating the aggressive features of the aforementioned BC subtype.