Infiltration of Common Myeloid Progenitor (CMP) Cells is Associated With Less Aggressive Tumor Biology, Lower Risk of Brain Metastasis, Better Response to Immunotherapy and Higher Patient Survival in Breast Cancer.

OBJECTIVE: To investigate the clinical relevance of common myeloid progenitor (CMP) cells in breast tumor microenvironment (TME). BACKGROUND: The role of rare cells in TME is less studied. In Silico transcriptomic analyses of real-world data enable us to detect and quantify rare cells, including CMP cells. METHODS: Total of 5,176 breast cancer (BC) patients from SCAN-B, METABRIC, and 5 single-cell sequence cohorts were analyzed using xCell algorithm. High group was defined as more than two thirds of CMP score in each cohort. RESULTS: CMP cells consist of 0.07-0.25% of bulk breast tumor cells, more in Estrogen Receptor-positive (ER+) compared with triple-negative (TN) subtype (0.1-0.75%, 0.18-0.33% of immune cells, respectively). CMP cells did not correlate with any of myeloid lineage nor stem cells in TME. CMP infiltration was higher in smaller tumors, with lower Nottingham grade, and in ER+/HER2- than in TNBC consistently in both SCAN-B and METABRIC cohorts. High CMP was significantly associated with lower risk of brain metastasis and with better survival, particularly in ER+/HER2- . High CMP enriched epithelial-to-mesenchymal transition and angiogenesis pathways, and less cell proliferation and DNA repair gene sets. High CMP ER+/HER2- was associated with less immune cell infiltration, and cytolytic activity (P<0.001). CMP infiltration correlated with neoadjuvant chemoimmunotherapy response for both ER+/HER2- and TNBC in the ISPY-2 cohort (AUC=0.69 and 0.74, respectively). CONCLUSIONS: CMP in BC is inversely associated with cell-proliferation and brain metastasis, better response to immunotherapy and survival. This is the first to report the clinical relevance of CMP infiltration in BC.

Eribulin promotes proliferation of CD8+ T cells and potentiates T cell-mediated anti-tumor activity against triple-negative breast cancer cells.

PURPOSE: Chemotherapeutic agents exert immunomodulatory effects on triple-negative breast cancer (TNBC) cells and immune cells. Eribulin favorably affects the immunological status of patients with breast cancer. However, the effects of eribulin on the immune cells remain unexplored. The aim of this study was to investigate the effects of eribulin on immune cells. METHODS: Peripheral blood mononuclear cells (PBMCs) from healthy donors and mouse splenocytes were stimulated with anti-CD3 and anti-CD28 antibodies. The effects of eribulin and paclitaxel on cell proliferation and differentiation status were analyzed using flow cytometry. RNA sequencing was performed to assess alterations in gene expression in CD8+ T cells following eribulin and paclitaxel treatment. Using TNBC cell lines (MDA-MB-231, Hs578T, and MDA-MB-157), the anti-tumor activity of CD3/CD28-stimulated T cells combined with eribulin or paclitaxel was evaluated. RESULTS: Eribulin did not affect CD3/CD28-stimulated PBMCs proliferation. However, eribulin significantly decreased the CD4/CD8 ratio in T cells, indicating that eribulin facilitates CD8+ T cell proliferation. Furthermore, eribulin significantly increased the frequency of less differentiated CD45RA+, CCR7+, and TCF1+ subsets of CD8+ T cells. RNA sequencing revealed that eribulin enhanced the expression of gene sets related to cell proliferation and immune responses. Moreover, eribulin augmented the anti-tumor effects of CD3/CD28-stimulated T cells against TNBC cells. These results were not observed in experiments using paclitaxel. CONCLUSIONS: Eribulin promoted CD8+ T cell proliferation, repressed effector T cell differentiation, and harnessed T cell-mediated anti-tumor effects. These mechanisms may be one of the cues that eribulin can improve the immunological status of tumor-bearing hosts.

CD133 expression is associated with less DNA repair, better response to chemotherapy and survival in ER-positive/HER2-negative breast cancer.

PURPOSE: CD133, a cancer stem cells (CSC) marker, has been reported to be associated with treatment resistance and worse survival in triple-negative breast cancer (BC). However, the clinical relevance of CD133 expression in ER-positive/HER2-negative (ER + /HER2-) BC, the most abundant subtype, remains unknown. METHODS: The BC cohorts from the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC, n = 1904) and The Cancer Genome Atlas (TCGA, n = 1065) were used to obtain biological variables and gene expression data. RESULTS: Epithelial cells were the exclusive source of CD133 gene expression in a bulk BC. CD133-high ER + /HER2- BC was associated with CD24, NOTCH1, DLL1, and ALDH1A1 gene expressions, as well as with WNT/ß-Catenin, Hedgehog, and Notch signaling pathways, all characteristic for CSC. Consistent with a CSC phenotype, CD133-low BC was enriched with gene sets related to cell proliferation, such as G2M Checkpoint, MYC Targets V1, E2F Targets, and Ki67 gene expression. CD133-low BC was also linked with enrichment of genes related to DNA repair, such as BRCA1, E2F1, E2F4, CDK1/2. On the other hand, CD133-high tumors had proinflammatory microenvironment, higher activity of immune cells, and higher expression of genes related to inflammation and immune response. Finally, CD133-high tumors had better pathological complete response after neoadjuvant chemotherapy in GSE25066 cohort and better disease-free survival and overall survival in both TCGA and METABRIC cohorts. CONCLUSION: CD133-high ER + /HER2- BC was associated with CSC phenotype such as less cell proliferation and DNA repair, but also with enhanced inflammation, better response to neoadjuvant chemotherapy and better prognosis.

Enhanced cancer cell proliferation and aggressive phenotype counterbalance in breast cancer with high BRCA1 gene expression.

PURPOSE: While comprehensive research exists on the mutation of the DNA repair gene BRCA1, limited information is available regarding the clinical significance of BRCA1 gene expression. Given that cancer cell proliferation is aggrevated by DNA repair, we hypothesized that high BRCA1 gene expression breast cancer (BC) might be linked with aggressive tumor biology and poor clinical outcomes. METHODS: The cohorts: The Cancer Genome Atlas (TCGA, n = 1069), METABRIC (n = 1903), and SCAN-B (n = 3273) were utilzed to obtain data of 6245 BC patients. RESULTS: BC patients without BRCA1 mutation exhibited higher BRCA1 expression, which was associated with DNA repair functionality. However, no such correlation was observed with BRCA2 expression. The association of high BRCA1 expression with cancer cell proliferation was evidenced by significant enrichment of cell proliferation-related gene sets, higher histological grade, and proliferation score. Furthermore, increased levels of homologous recombination deficiency, intratumoral heterogeneity, and altered fractions were associated with high BRCA1 expression. Moreover, BC with high BRCA1 expression exhibited reduced infiltration of dendritic cells and CD8 T-cells, while showing increased infiltration of Th1 cells. Surprisingly, BRCA1 expression was not associated with the survival of BC irrespective of the subtypes. Conversely, BC with low BRCA1 expression enriched cancer aggravating pathway gene sets, such as Cancer Stem Cell-related signaling (NOTCH and HEDGEHOG), Angiogenesis, Epithelial-Mesenchymal Transition, Inflammatory Response, and TGF-beta signaling. CONCLUSION: Despite being linked to heightened proliferation of cancer cells and unassertive phenotype, BRCA1 expression did not show any association with survival in BC.

SLC7A11 and the glutathione pathway as novel prognostic markers in resectable pancreatic ductal adenocarcinoma: A metabolomics study of clinical specimens.

BACKGROUND/OBJECTIVES: Despite the poor prognosis associated with pancreatic ductal adenocarcinoma (PDAC), there remains a lack of clarity regarding the metabolic pathways and their significant impact on its phenotype. Therefore, we aimed to utilize metabolomics to capture changes in clinical PDAC tissues and elucidate the significant metabolic pathways close to its phenotypes. METHODS: This basic research was retrospectively validated using database research, immunohistochemistry, and protein analysis based on the findings obtained from metabolomics using clinical tissues collected from prospectively registered patients with PDAC. mRNA expression analysis using a database and protein analysis using archived clinical specimens was performed to validate the candidate pathways identified using metabolomics. Between-group comparisons were analyzed using paired t-tests and log-rank test, and Kaplan-Meier curves illustrated survival times. RESULTS: Patients subjected to metabolomics revealed a significant increase in glutathione disulfide levels in PDAC tissues when compared to normal pancreatic tissues. The Cancer Genome Atlas database analysis revealed significant changes in glutathione pathway-related mRNAs in PDAC compared to that in the normal pancreas. Protein analysis of previously resected specimens demonstrated a significant increase in SLC7A11 expression in PDAC tissues. The abundance ratio of SLC7A11 isoforms was associated with the post-operative prognosis in resectable PDAC. CONCLUSION: Glutathione disulfide levels were significantly increased in clinical PDAC metabolomics. Additionally, increased mRNA and protein expression in SLC7A11 was observed in PDAC. Furthermore, the SLC7A11 isoform abundance ratio may be a valuable prognostic marker in patients with resectable PDAC.

Emerging measurements for tumor-infiltrating lymphocytes in breast cancer.

Tumor-infiltrating lymphocytes are a general term for lymphocytes or immune cells infiltrating the tumor microenvironment. Numerous studies have demonstrated tumor-infiltrating lymphocytes to be robust prognostic and predictive biomarkers in breast cancer. Recently, immune checkpoint inhibitors, which directly target tumor-infiltrating lymphocytes, have become part of standard of care treatment for triple-negative breast cancer. Surprisingly, tumor-infiltrating lymphocytes quantified by conventional methods do not predict response to immune checkpoint inhibitors, which highlights the heterogeneity of tumor-infiltrating lymphocytes and the complexity of the immune network in the tumor microenvironment. Tumor-infiltrating lymphocytes are composed of diverse immune cell populations, including cytotoxic CD8-positive T lymphocytes, B cells and myeloid cells. Traditionally, tumor-infiltrating lymphocytes in tumor stroma have been evaluated by histology. However, the standardization of this approach is limited, necessitating the use of various novel technologies to elucidate the heterogeneity in the tumor microenvironment. This review outlines the evaluation methods for tumor-infiltrating lymphocytes from conventional pathological approaches that evaluate intratumoral and stromal tumor-infiltrating lymphocytes such as immunohistochemistry, to the more recent advancements in computer tissue imaging using artificial intelligence, flow cytometry sorting and multi-omics analyses using high-throughput assays to estimate tumor-infiltrating lymphocytes from bulk tumor using immune signatures or deconvolution tools. We also discuss higher resolution technologies that enable the analysis of tumor-infiltrating lymphocytes heterogeneity such as single-cell analysis and spatial transcriptomics. As we approach the era of personalized medicine, it is important for clinicians to understand these technologies.

Intratumoral Tumor Infiltrating Lymphocytes (TILs) are Associated With Cell Proliferation and Better Survival But Not Always With Chemotherapy Response in Breast Cancer.

OBJECTIVE: To investigate the clinical relevance of intratumoral tumor infiltrating lymphocytes (TILs) in breast cancer as measured by computational deconvolution of bulk tumor transcriptomes. SUMMARY BACKGROUND DATA: Commonly assessed TILs, located in tumor stroma without direct contact with cancer cells (stromal TILs), correlate with breast cancer treatment response and survival. The clinical relevance of intratumoral TILs has been less studied partly due to their rarity; however, they may have nonnegligible effects given their direct contact with cancer cells. METHODS: In all, 5870 breast cancer patients from TCGA, METABRIC, GSE96058, GSE25066, GSE163882, GSE123845, and GSE20271 cohorts were analyzed and validated. RESULTS: The intratumoral TIL score was established by the sum of all types of lymphocytes using the xCell algorithm. This score was the highest in triple-negative breast cancer (TNBC) and the lowest in the ER-positive/HER2-negative subtype. It correlated with cytolytic activity and infiltrations of dendritic cells, macrophages, and monocytes, and uniformly enriched immune-related gene sets regardless of subtype. Intratumoral TIL-high tumors correlated with higher mutation rates and significant cell proliferation on biological, pathological, and molecular analyses only in the ER-positive/HER2-negative subtype. It was significantly associated with pathological complete response after anthracycline- and taxane-based neoadjuvant chemotherapy in about half of the cohorts, regardless of the subtype. Intratumoral TIL-high tumors correlated with better overall survival in HER2-positive and TNBC subtypes consistently in 3 cohorts. CONCLUSIONS: Intratumoral TILs estimated by transcriptome computation were associated with increased immune response and cell proliferation in ER-positive/HER2-negative and better survival in HER2-positive and TNBC subtypes, but not always with pathological complete response after neoadjuvant chemotherapy.

Intratumoral PDGFB gene predominantly expressed in endothelial cells is associated with angiogenesis and lymphangiogenesis, but not with metastasis in breast cancer.

PURPOSE: Platelet-derived growth factor B (PDGFB) is known to play essential roles in angiogenesis and lymphangiogenesis during development, and tumor growth and vessel stabilization in experimental models. However, whether these findings could be translated to breast cancer patients remains unclear. We hypothesized that PDGFB gene expression is associated with angiogenesis, cell proliferation, and clinical outcomes in breast cancer patients. METHODS: A total of 7635 primary breast cancer patients with full transcriptome and clinical data available from 13 independent cohorts were analyzed using in silico approach. The median value was used to divide each cohort into high and low PDGFB expression groups. RESULTS: High PDGFB gene expression was associated with increased expression of angiogenesis-related genes, higher amount of vascular cell infiltrations, and with enrichment of angiogenesis gene set, lymphangiogenesis-related gene expressions, lymphangiogenesis-related cell infiltrations, and enrichmentof lymphangiogenesis gene set in GSE96058 and validated by TCGA cohorts; however, not with lymphatic metastasis. PDGFB expression was neither associated with cell proliferation as assessed by Ki67 expression nor with Nottingham histological grade, or response to neoadjuvant chemotherapy. We found that PDGFB was most extensively expressed by endothelial and perivascular-like cells in the tumor microenvironment, and minimally by cancer cells consistently in two single-cell sequence cohorts. High PDGFB expression enriched TGFß, epithelial-mesenchymal transition, hypoxia, and cancer stem cell-associated pathways. However, no association with distant metastasis was observed. Disease-specific and disease-free survival were worse in the high PDGFB expression group consistently in TCGA and METABRIC cohorts. CONCLUSION: PDGFB is predominantly expressed in endothelial cells and is associated with angiogenesis and lymphangiogenesis, but not with cellular proliferation or metastasis in breast cancer.

High RAD51 gene expression is associated with aggressive biology and with poor survival in breast cancer.

PURPOSE: Although the DNA repair mechanism is important in preventing carcinogenesis, its activation in established cancer cells may support their proliferation and aggravate cancer progression. RAD51 cooperates with BRCA2 and is essential in the homologous recombination of DNA repair. To this end, we hypothesized that RAD51 gene expression is associated with cancer cell proliferation and poor prognosis of breast cancer (BC) patients. METHODS: A total of 8515 primary BC patients with transcriptome and clinical data from 17 independent cohorts were analyzed. The median value was used to divide each cohort into high and low RAD51 expression groups. RESULTS: High RAD51 expression enriched the DNA repair gene set and was correlated with DNA repair-related genes. Nottingham histological grade, Ki67 expression and cell proliferation-related gene sets (E2F Targets, G2M Checkpoint and Myc Targets) were all significantly associated with the high RAD51 BC group. RAD51 expression was positively correlated with Homologous Recombination Deficiency, as well as both mutational burden and neoantigens that accompanied a higher infiltration of immune cells. Primary BC with lymph node metastases was associated with high expression of RAD51 in two cohorts. There was no strong correlation between RAD51 expression and drug sensitivity in cell lines, and RAD51 expression was lower after the neoadjuvant chemotherapy compared to before the treatment. High RAD51 BC was associated with poor prognosis consistently in three independent cohorts. CONCLUSION: RAD51 gene expression is associated with aggressive cancer biology, cancer cell proliferation, and poor survival in breast cancer.

Abundance of reactive oxygen species (ROS) is associated with tumor aggressiveness, immune response, and worse survival in breast cancer.

PURPOSE: Reactive oxygen species (ROS) are oxygen-containing molecules that have high reactivity and play roles in protection or harm the cancer cells. We aimed to clarify the clinical relevance of ROS in breast cancer (BC) tumor microenvironment (TME). We hypothesized that it is associated with worse BC patient outcomes. METHODS: ROS score was generated by Gene Set Variation Analysis of Hallmark ROS pathway gene set and a total of 6245 BC patients were analyzed. RESULTS: High ROS BC significantly enriched cell proliferation-related gene sets (MYC targets v1 and v2, G2M checkpoint, E2F targets), pro-cancer-related gene sets (DNA repair, unfolded protein response, MTORC1 signaling, PI3K/AKT/MTOR signaling, glycolysis, and oxidative phosphorylation), immune-related gene sets (inflammatory response, allograft rejection, interferon-α and γ responses, complement, and IL6/JAK/STAT3 signaling), and infiltrated immune cells (CD4+ memory and CD8+ T cells, Th1 and Th2, dendritic cells, Tregs, M1 and M2 macrophages) and B cells, as well as elevated cytolytic activity consistently in both METABRIC and GSE96058 cohorts. Cancer cells were the major source of ROS in BC TME of single-cell sequence (GSE75688) cohort. High ROS was associated with intratumor heterogeneity, homologous recombination defects, mutation rates, and neoantigens, and with clinical aggressiveness in AJCC stage, Nottingham grade and Ki67 expression, as well as worse overall survival in both GSE96058 and METABRIC, and with worse disease-specific survival in METABRIC. CONCLUSION: Abundant ROS in BC patients is associated with abundant mutations, aggressive cancer biology, immune response, and worse survival.

Low RUFY3 expression level is associated with lymph node metastasis in older women with invasive breast cancer.

PURPOSE: Sentinel lymph node biopsy is omitted in older women (≥ 70 years old) with clinical lymph node (LN)-negative hormone receptor-positive breast cancer as it does not influence adjuvant treatment decision-making. However, older women are heterogeneous in frailty while the chance of recurrence increase with improving longevity. Therefore, a biomarker that identifies LN metastasis may facilitate treatment decision-making. RUFY3 is associated with cancer progression. We evaluated RUFY3 expression level as a biomarker for LN-positive breast cancer in older women. METHODS: Clinical and transcriptomic data of breast cancer patients were obtained from the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC, n = 1903) and The Cancer Genome Atlas (TCGA, n = 1046) Pan-cancer study cohorts. RESULTS: A total of 510 (METABRIC) and 211 (TCGA) older women were identified. LN-positive breast cancer, which represented 51.4% (METABRIC) and 48.4% (TCGA), demonstrated worse disease-free, disease-specific, and overall survival. RUFY3 levels were significantly lower in LN-positive tumors regardless of age. The area under the curve for the receiver operator characteristic (AUC-ROC) curves showed RUFY3-predicted LN metastasis. Low RUFY3 enriched oxidative phosphorylation, DNA repair, MYC targets, unfolded protein response, and mtorc1 signaling gene sets, was associated with T helper type 1 cell infiltration, and with intratumor heterogeneity and fraction altered. Low RUFY3 expression was associated with LN-positive breast cancer and with worse disease-specific survival among older women. CONCLUSION: Older women with breast cancers who had low expression level of RUFY3 were more frequently diagnosed with LN-positive tumors, which translated into worse prognosis.

High BRCA2 Gene Expression is Associated with Aggressive and Highly Proliferative Breast Cancer.

BACKGROUND: Mutations of BRCA genes are the most studied in breast cancer, but the clinical relevance of BRCA2 gene expression has been less well studied. Given that BRCA2 is a DNA repair gene, we hypothesized that high BRCA2 expression is associated with highly proliferative and aggressive biology in breast cancer. MATERIALS AND METHODS: A total of 4342 breast cancer patients were analyzed from The Cancer Genome Atlas (TCGA, n = 1069) as the testing cohort and Gene Expression Omnibus (GEO) dataset GSE96058 (n = 3273) as a validation cohort. RESULTS: There was no relationship between BRCA2 mutation and BRCA2 gene expression. BRCA2 high expression breast cancer was associated with higher Nottingham grade (p < 0.001), with high proliferation (MKI-67, p < 0.001), and with higher intratumor heterogeneity, homologous recombination deficiency, mutation rate, fraction altered, and neoantigens (all p < 0.001). BRCA2 high expression was associated with E2F1, RB1, PALB2, and PARP, as well as cell proliferation-related gene sets, E2F targets, G2M checkpoints, and mitotic spindle, and with less ESR1 and ER response early and late gene sets. They were associated with significantly reduced number of stroma cells and with high infiltration of both favorable and unfavorable immune cells. BRCA2 high expression significantly correlated with response to olaparib, a PARP inhibitor, and inversely with cyclophosphamide in ER-positive/HER2-negative tumors, but not in TNBC. CONCLUSIONS: BRCA2 high gene expression was associated with highly proliferative and aggressive breast cancer that was highly immunogenic with better response to PARP inhibitors in ER-positive patients. BRCA2 gene expression may become a candidate marker for aggressive biology and to tailor specific treatment strategies in the future.

Orthotopic Implantation Achieves Better Engraftment and Faster Growth Than Subcutaneous Implantation in Breast Cancer Patient-Derived Xenografts.

Patient-Derived Xenograft (PDX) is now accepted as a murine model that better mimics human cancer when compared to a conventional cancer cell-line inoculation model. Some claim the advantage of orthotopic site implantation of patient tumor (OS) over ectopic implantation into the subcutaneous space (SQ); however, there has been no study that describes a head-to-head comparison of oncological differences between these two models to date. We hypothesize that OS tumors re-transplant and grow better than SQ tumors and are therefore a better model to evaluate tumor aggressiveness. Breast cancer PDXs were generated using the tumors derived from 11 patients into NOD scid gamma (NSG) mice. We used six ER(+)HER2(-) tumors and five triple negative (TN) tumors for a total of 11 tumors. Five PDX lines grew for an overall engraftment rate of 45%. We present our OS implantation method in detail. The re-transplantation rate of TN tumors in each transplant site was significantly higher in OS when compared to SQ tumors (70.1% vs. 32.1%, p < 0.01). OS tumors grow significantly faster than SQ tumors. Similarly, OS tumors demonstrated significantly more mitotic figures and Ki-67 positive cells than SQ tumors. The tumor re-transplantation rate significantly increased by the second and third generations with the OS method. The time from implantation to development of a palpable tumor dramatically decreased after the first passage. PDX of ER(+) tumors demonstrated significantly lower engraftment rates and slower tumor growth than TN tumors, which remarkably improved by the first passage. Orthotopically implanted PDX tumors showed better re-transplantation rates, greater tumor size, and more significant growth compared to the subcutaneously implanted model.

Molecular Biological Features of Nottingham Histological Grade 3 Breast Cancers.

INTRODUCTION: Cancer biology dominates the behavior and prognosis of a tumor. Although Nottingham histological grade is a subjective pathological determination, it has been accepted as a surrogate model for cancer biology. As such, histologic grade was incorporated into the latest 8th edition of the American Joint Committee on Cancer breast cancer staging system. In this study, we hypothesized that grade 3 breast cancers demonstrate aggressive molecular biological profiles, reflecting worse biology and possible underlying immunogenicity. METHODS: Transcriptomic and clinical data were obtained from the Molecular Taxonomy of Breast Cancer International Consortium, and the findings were validated by The Cancer Genome Atlas breast cancer cohort and GSE25066. RESULTS: Overall, 2876 patients were analyzed in this study. Grade 3 tumors were more common in estrogen receptor (ER)-negative, advanced-stage patients, and were associated with human epidermal growth factor receptor 2 and basal subtypes by the PAM50 classifier, as well as with increased MKI67 expression (all p <0.001). Disease-free survival was significantly worse in grade 3 tumors (all cohorts). Gene set enrichment analysis demonstrated that grade 3 tumors were significantly enriched with not only cell proliferation and cell cycle-related gene sets but also immune activity-related gene sets. CIBERSORT confirmed that grade 3 tumors were infiltrated with macrophage M1, follicular helper T cells, and activated natural killer cells (all p <0.001). Furthermore, grade 3 tumors were associated with more diverse T cell receptors (p =0.001) and increased cytolytic activity (p <0.001). Lastly, major T-cell exhaustion markers were significantly elevated in grade 3 breast cancers (p <0.001). CONCLUSION: Grade 3 breast cancers demonstrated aggressive transcriptomic features with enhanced immunogenicity and elevated T-cell exhaustion markers.

MYC Targets Scores Are Associated with Cancer Aggressiveness and Poor Survival in ER-Positive Primary and Metastatic Breast Cancer.

MYC is one of the most studied oncogenes that is known to promote cell proliferation. We utilized MYC targets v1 and MYC targets v2 scores of gene set variation analysis and hypothesized that these scores correlate with tumor aggressiveness and survival outcomes. We examined a total of 3109 breast cancer patients from TCGA, METABRIC, and GSE124647 cohorts. In each cohort, the patients were divided into high- and low-score groups using the upper third value as the cut off. As expected, higher scores were related to increased cell proliferation and worse clinical and pathologic features. High MYC targets scores were associated with worse survival, specifically in primary ER-positive breast cancer, consistently in both TCGA and METABRIC cohorts. In ER-positive breast cancer, high MYC targets v1, but not v2 score, was associated with high mutation load, and high MYC targets v1 and v2 scores were both associated with increased infiltration of pro- and anti-cancerous immune cells. We found that high MYC scores were associated with worse survival in metastatic breast cancer. Our findings show that the MYC targets v1 and v2 scores are associated with tumor aggressiveness and poor prognosis in ER-positive primary tumors, as well as in metastatic breast cancer.

CD8 T Cell Score as a Prognostic Biomarker for Triple Negative Breast Cancer.

CD8 T cell is an essential component of tumor-infiltrating lymphocytes (TIL) and tumor immune microenvironment (TIME). Using the xCell CD8 T cell score of whole tumor gene expression data, we estimated these cells in total of 3837 breast cancer patients from TCGA, METABRIC and various GEO cohorts. The CD8 score correlated strongly with expression of CD8 genes. The score was highest for triple-negative breast cancer (TNBC), and a high score was associated with high tumor immune cytolytic activity and better survival in TNBC but not other breast cancer subtypes. In TNBC, tumors with a high CD8 score had enriched expression of interferon (IFN)-α and IFN-γ response and allograft rejection gene sets, and greater infiltration of anti-cancerous immune cells. The score strongly correlated with CD4 memory T cells in TNBC, and tumors with both a high CD8 score and high CD4 memory T cell abundance had significantly better survival. Finally, a high CD8 score was significantly associated with high expression of multiple immune checkpoint molecules. In conclusion, a high CD8 T cell score is associated with better survival in TNBC, particularly when tumor CD4 memory T cells were elevated. Our findings also suggest a possible use of the score as a predictive biomarker for response to immune checkpoint therapy.

Intra-Tumoral Angiogenesis Is Associated with Inflammation, Immune Reaction and Metastatic Recurrence in Breast Cancer.

Angiogenesis is one of the hallmarks of cancer. We hypothesized that intra-tumoral angiogenesis correlates with inflammation and metastasis in breast cancer patients. To test this hypothesis, we generated an angiogenesis pathway score using gene set variation analysis and analyzed the tumor transcriptome of 3999 breast cancer patients from The Cancer Genome Atlas Breast Cancer (TCGA-BRCA), Molecular Taxonomy of Breast Cancer International Consortium (METABRIC), GSE20194, GSE25066, GSE32646, and GSE2034 cohorts. We found that the score correlated with expression of various angiogenesis-, vascular stability-, and sphingosine-1-phosphate (S1P)-related genes. Surprisingly, the angiogenesis score was not associated with breast cancer subtype, Nottingham pathological grade, clinical stage, response to neoadjuvant chemotherapy, or patient survival. However, a high score was associated with a low fraction of both favorable and unfavorable immune cell infiltrations except for dendritic cell and M2 macrophage, and with Leukocyte Fraction, Tumor Infiltrating Lymphocyte Regional Fraction and Lymphocyte Infiltration Signature scores. High-score tumors had significant enrichment for unfavorable inflammation-related gene sets (interleukin (IL)6, and tumor necrosis factor (TNF)α- and TGFß-signaling), as well as metastasis-related gene sets (epithelial mesenchymal transition, and Hedgehog-, Notch-, and WNT-signaling). High score was significantly associated with metastatic recurrence particularly to brain and bone. In conclusion, using the angiogenesis pathway score, we found that intra-tumoral angiogenesis is associated with immune reaction, inflammation and metastasis-related pathways, and metastatic recurrence in breast cancer.

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.

G2M Cell Cycle Pathway Score as a Prognostic Biomarker of Metastasis in Estrogen Receptor (ER)-Positive Breast Cancer.

The vast majority of breast cancer death is a result of metastasis. Thus, accurate identification of patients who are likely to have metastasis is expected to improve survival. The G2M checkpoint plays a critical role in cell cycle. We hypothesized that breast cancer tumors with high activity of G2M pathway genes are more aggressive and likely to metastasize. To test this, we used the single-sample gene set variation analysis method to calculate the score for the Hallmark G2M checkpoint pathway using gene expression data of a total of 4626 samples from 12 human breast cancer cohorts. As expected, a high G2M pathway score correlated with enriched tumor expression of other cell proliferation-related gene sets. The score was significantly associated with clinical aggressive features of tumors and patient survival in estrogen receptor (ER)-positive/human epidermal growth factor receptor 2 (HER2)-negative breast cancer. Interestingly, a high G2M score of metastasis tumors was also significantly associated with worse survival. In primary as well as metastasis tumors with high scores, the infiltration of both pro- and anti-cancerous immune cells increased. Tumor G2M score was also associated with treatment response to systemic chemotherapy in ER-positive/HER2-negative cancer, and was predictive of response to cyclin-dependent kinase inhibition therapy.

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.