High MYC mRNA Expression Is More Clinically Relevant than MYC DNA Amplification in Triple-Negative Breast Cancer

C-terminally phosphorylated p27 activates self-renewal driver genes to program cancer stem cell expansion, mammary hyperplasia and cancer

In many cancers, a stem-like cell subpopulation mediates tumor initiation, dissemination and drug resistance. Here, we report that cancer stem cell (CSC) abundance is transcriptionally regulated by C-terminally phosphorylated p27 (p27pT157pT198). Mechanistically, this arises through p27 co-recruitment with STAT3/CBP to gene regulators of CSC self-renewal including MYC, the Notch ligand JAG1, and ANGPTL4. p27pTpT/STAT3 also recruits a SIN3A/HDAC1 complex to co-repress the Pyk2 inhibitor, PTPN12. Pyk2, in turn, activates STAT3, creating a feed-forward loop increasing stem-like properties in vitro and tumor-initiating stem cells in vivo. The p27-activated gene profile is over-represented in STAT3 activated human breast cancers. Furthermore, mammary transgenic expression of phosphomimetic, cyclin-CDK-binding defective p27 (p27CK-DD) increases mammary duct branching morphogenesis, yielding hyperplasia and microinvasive cancers that can metastasize to liver, further supporting a role for p27pTpT in CSC expansion. Thus, p27pTpT interacts with STAT3, driving transcriptional programs governing stem cell expansion or maintenance in normal and cancer tissues.

Dual activity of Minnelide chemosensitize basal/triple negative breast cancer stem cells and reprograms immunosuppressive tumor microenvironment

Triple negative breast cancer (TNBC) subtype is characterized with higher EMT/stemness properties and immune suppressive tumor microenvironment (TME). Women with advanced TNBC exhibit aggressive disease and have limited treatment options. Although immune suppressive TME is implicated in driving aggressive properties of basal/TNBC subtype and therapy resistance, effectively targeting it remains a challenge. Minnelide, a prodrug of triptolide currently being tested in clinical trials, has shown anti-tumorigenic activity in multiple malignancies via targeting super enhancers, Myc and anti-apoptotic pathways such as HSP70. Distinct super-enhancer landscape drives cancer stem cells (CSC) in TNBC subtype while inducing immune suppressive TME. We show that Minnelide selectively targets CSCs in human and murine TNBC cell lines compared to cell lines of luminal subtype by targeting Myc and HSP70. Minnelide in combination with cyclophosphamide significantly reduces the tumor growth and eliminates metastasis by reprogramming the tumor microenvironment and enhancing cytotoxic T cell infiltration in 4T1 tumor-bearing mice. Resection of residual tumors following the combination treatment leads to complete eradication of disseminated tumor cells as all mice are free of local and distant recurrences. All control mice showed recurrences within 3 weeks of post-resection while single Minnelide treatment delayed recurrence and one mouse was free of tumor. We provide evidence that Minnelide targets tumor intrinsic pathways and reprograms the immune suppressive microenvironment. Our studies also suggest that Minnelide in combination with cyclophosphamide may lead to durable responses in patients with basal/TNBC subtype warranting its clinical investigation.

Single-cell profile of tumor and immune cells in primary breast cancer, sentinel lymph node, and metastatic lymph node

PURPOSE

Little is known about the host-tumor interaction in the lymph-node basin at a single cell level. This study examines single cell sequences in breast cancer nodal metastases of a patient with triple-negative breast cancer.

METHODS

The primary breast tumor, sentinel lymph node, an adjacent lymph node with metastatic involvement and a clinically normal-appearing lymph node were collected during surgery. Single-cell sequencing was performed on all four specimens.

RESULTS

14,016 cells were clustered into 6 cell subpopulations. Cancer cells demonstrated the molecular characteristics of TNBC basal B subtype and highly expressed genes in the MAPK signaling cascade. Tumor-associated macrophages regulated antigen processing and presentation and other immune-related pathways to promote tumor invasion. CD8 + and CD4 + T lymphocytes concentrated more in sentinel lymph node and mainly stratified into two transcriptional states. The immune-cell amount variation among primary tumor, sentinel and normal lymph nodes showed a similar tendency between the sc-RNA-seq profile of TNBC samples and a previous reported bulk RNA-seq profile of a breast cancer cohort, including all four breast cancer subtype samples.

DISCUSSION

Single-cell sequencing analysis suggested that the sentinel lymph node was the initial meeting site of tumor infiltration and immune response, where partial T lymphocytes perform anti-tumor activity, while other T cells exhibit an exhausted state. We proposed a molecular explanation to the well-established clinical principle that the 5-year and 10-year survival outcomes were noninferior between SLND and ALND.

Inhibition of MYC suppresses programmed cell death ligand-1 expression and enhances immunotherapy in triple-negative breast cancer

BACKGROUND

Cancer immunotherapy has emerged as a promising strategy against triple-negative breast cancer (TNBC). One of the immunosuppressive pathways involves programmed cell death-1 (PD-1) and programmed cell death ligand-1 (PD-L1), but many patients derived little benefit from PD-1/PD-L1 checkpoint blockades treatment. Prior research has shown that MYC, a master transcription amplifier highly expressed in TNBC cells, can regulate the tumor immune microenvironment and constrain the efficacy of immunotherapy. This study aims to investigate the regulatory relationship between MYC and PD-L1, and whether a cyclin-dependent kinase (CDK) inhibitor that inhibits MYC expression in combination with anti-PD-L1 antibodies can enhance the response to immunotherapy.

METHODS

Public databases and TNBC tissue microarrays were used to study the correlation between MYC and PD-L1. The expression of MYC and PD-L1 in TNBCs was examined by quantitative real-time polymerase chain reaction and Western blotting. A patient-derived tumor xenograft (PDTX) model was used to evaluate the influence of a CDK7 inhibitor THZ1 on PD-L1 expression. Cell proliferation and migration were detected by 5-ethynyl-2'-deoxyuridine (EdU) cell proliferation and cell migration assays. Tumor xenograft models were established for in vivo verification.

RESULTS

A high MYC expression level was associated with a poor prognosis and could alter the proportion of tumor-infiltrating immune cells (TIICs). The positive correlation between MYC and PD-L1 was confirmed by immunostaining samples from 165 TNBC patients. Suppression of MYC in TNBC caused a reduction in the levels of both PD-L1 messenger RNA and protein. In addition, antitumor immune response was enhanced in the TNBC cancer xenograft mouse model with suppression of MYC by CDK7 inhibitor THZ1.

CONCLUSIONS

The combined therapy of CDK7 inhibitor THZ1 and anti-PD-L1 antibody appeared to have a synergistic effect, which might offer new insight for enhancing immunotherapy in TNBC.

Gastric cancer with enhanced apical junction pathway has increased metastatic potential and worse clinical outcomes

Excessive intercellular connection at confluency may be limiting further cell growth or a sign of aggressive biology in the cell culture. As apical junction complex is a main component of cell-to-cell connection, we aimed to investigate gastric cancer biology using Apical Junction Pathway score that we generated using Gene set variant analysis (GSVA) of the "Hallmark Apical Junction" gene set. 1,239 gastric cancer patients from the Cancer Genome Atlas (TCGA) and two GSE cohorts were included in this study. The cohorts were dichotomized using the median of the score. Apical Junction Pathway score high gastric cancer was not consistently associated with increased cell proliferation or immune cell infiltration. On the other hand, Apical Junction Pathway score high gastric cancer was associated with significantly higher infiltration of stromal cells, such as endothelial cells; hence, increased neovascularization and angiogenesis in the tumor microenvironment (TME) were speculated. Gene set enrichment analysis (GSEA) confirmed increased expression of epithelial mesenchymal transition (EMT) and angiogenesis in the high Apical Junction Pathway score group (false discovery rate (FDR) <0.25). Lastly, the high Apical Junction Pathway score group was associated with more aggressive clinicopathological characteristics, such as significantly higher American Joint Committee on Cancer (AJCC) T-category and higher pathological stage, leading to worse disease-specific survival and overall survival (P<0.05, respectively). In conclusion, enhanced Apical Junction Pathway score gastric cancer was associated with aggressive clinical characteristics leading to shorter survival likely due to increased metastatic potential from EMT and angiogenesis.

Inhibition of the CDK9-cyclin T1 protein-protein interaction as a new approach against triple-negative breast cancer

Cyclin-dependent kinase 9 (CDK9) activity is correlated with worse outcomes of triple-negative breast cancer (TNBC) patients. The heterodimer between CDK9 with cyclin T1 is essential for maintaining the active state of the kinase and targeting this protein-protein interaction (PPI) may offer promising avenues for selective CDK9 inhibition. Herein, we designed and generated a library of metal complexes bearing the 7-chloro-2-phenylquinoline CˆN ligand and tested their activity against the CDK9-cyclin T1 PPI. Complex 1 bound to CDK9 via an enthalpically-driven binding mode, leading to disruption of the CDK9-cyclin T1 interaction in vitro and in cellulo. Importantly, complex 1 showed promising anti-metastatic activity against TNBC allografts in mice and was comparably active compared to cisplatin. To our knowledge, 1 is the first CDK9-cyclin T1 PPI inhibitor with anti-metastatic activity against TNBC. Complex 1 could serve as a new platform for the future design of more efficacious kinase inhibitors against cancer, including TNBC.

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.

Combined mRNAs and clinical factors model on predicting prognosis in patients with triple-negative breast cancer

OBJECTIVE

Triple-negative breast cancer (TNBC) is aggressive cancer usually diagnosed in young women with no effective prognosis prediction model to use. The present study was performed to develop a useful prognostic model for predicting overall survival (OS) for TNBC patients.

METHODS

The Cancer Genome Atlas (TCGA) and Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) databases were used as training and validation data sets, respectively, in which the gene expression levels and clinical prognostic information of TNBC were collected. Differentially expressed genes (DEGs) between TNBC and non-TNBC (NTNBC) were identified with the thresholds of false discovery rate < 0.05 and |log2 Fold Change| > 1. DEGs in AmiGO2 and the Kyoto Encyclopedia of Genes and Genomes (KEGG) databases were retained for further study. Univariate, multivariate Cox, and logistic regression analysis were conducted for detecting DEG signature with the threshold of log-rank P < 0.05. The prognosis models of mRNA signature, clinical factors were constructed and compared.

RESULTS

One five-DEG signature, including CHST4, COCH, CST9, SOX11, and TDGF1 was identified in DEG prognosis model. Stratified analysis showed that the patients aged over 60, with higher pathologic stage (III-IV) and recurrence induced a significantly lower survival rate than those aged below 60, lower pathologic stage and without recurrence. Compared with patients with low-risk scores, those presented high-risk scores demonstrated significantly lower survival rate in the subgroup aged over 60 [HR = 3.780 (1.801-7.933), P < 0.0001]. For patients who obtained a higher pathologic stage and recurrence, high-risk scores were correlated with a significantly lower survival rate than patients with low-risk scores. The five-mRNA signature combined with clinical model (AUC = 0.950) predicted better than single clinical model (AUC = 0.795) or five-mRNA signature model (AUC = 0.823).

CONCLUSION

Our present study identified a prognostic prediction model (combined with five-mRNA signature and clinical factors) for TNBC patients receiving immunotherapy, which will benefit future research and clinical therapies.

Clinical relevance of stem cell surface markers CD133, CD24, and CD44 in colorectal cancer

Colon cancer stem cells (CSC) identified by cell surface markers CD133, CD24, and CD44, have been shown to be involved with tumor formation, chemotherapy resistance, and the progression of metastatic disease. Using an in silico translational approach, we hypothesize that a combination of these CSC markers has prognostic value in a large cohort of patients with colorectal cancer. Clinicopathologic and RNA expression data from a total of 594 colorectal cancer (CRC) patients from TCGA were analyzed. The expression of CD133, CD24, and CD44 was individually defined as "high" or "low" based on the median expression. Disease specific survival (DSS) and overall survival (OS) were not associated with tumors that are CD133-high or CD44-high alone. Patients with CD24-high tumors have significantly better DSS (P<0.001) and OS (P = 0.043). CD24-high, CD44-high and CD133-high tumors were associated with significantly greater EGFR, KRAS and Ki67 expression (all P<0.001). CD133, CD24 and CD44-high tumors were independently enriched for conventional stemness-related signaling pathways such as Wnt/β-catenin and Hedgehog signaling pathways. There was no survival difference linked to CD133-high/CD44-low patients, but CD44-high/CD24-low patients have worse DSS (P = 0.005) compared with CD44-low/CD24-high patients. CD133-high/CD24-low tumors show significant negative enrichment of MYC targets, E2F targets, G2M checkpoint and mitotic spindle gene sets, suggesting less cell proliferation in these tumors. Patients with CD133-high/CD24-low tumors have worse DSS (P = 0.004) and OS (P = 0.044), and are more likely to have early and late recurrences. In conclusion, we demonstrated that CD133-high/CD24-low tumors may predict colorectal cancer prognosis.

A prognostic score based on long-term survivor unique transcriptomic signatures predicts patient survival in pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) is known for its poor prognosis with few long-term survivors. This study aimed to establish a prognostic score using unique transcriptomic profiles of long-term survivors to be used as a patient selection tool for meaningful clinical intervention in PDAC. In TCGA PDAC cohort, 16 genes were significantly upregulated in the long-term survivor tumors. A prognostic score was established using these 16 genes by LASSO Cox regression, and PHKG1, HOXA4, ISL2, DMRT3 and TRA2A gene expressions were included in the score. The prognostic value was confirmed in both testing and validation cohorts. The characteristics of the high score tumor was investigated by bioinformatical approach. The high score tumor was associated with TP53 mutation but not with other commonly enhanced signaling pathways in PDAC. The high score tumor was associated with higher tumor mutational burden and unfavorable tumor microenvironment (TME), such as lower infiltration of CD8-positive T cells and dendritic cells, and less cell composition of mature blood vessels and fibroblasts. The high score tumor was also associated with enhanced cell proliferation and margin positivity after surgery. The impact of score component genes on the cell proliferation was investigated by in vitro experiments. Silencing of the score component genes promoted cell proliferation. In conclusion, the prognostic score predicted PDAC patient survival and was associated with cancer aggressiveness such as unfavorable TME and enhanced cell proliferation.

Increased apoptosis is associated with robust immune cell infiltration and cytolytic activity in breast cancer

Tumor infiltrating immune cells plays a critical role in cancer progression. Apoptosis is an autonomous cell death that counteracts tumor growth. To this end, we hypothesized that increased apoptosis in breast cancer is associated with immune cell killing. Apoptosis score of MSigDB Hallmark collection was used to analyze METABRIC cohort (n=1904) and TCGA (n=1069) as validation cohort. High apoptosis tumors enriched cancer promoting signaling pathways; hypoxia, KRAS, TGF-β, PI3K signaling, and was associated with low MKI67 expression and less cell proliferation gene sets, less homologous recombination defects, and less altered fraction. High apoptosis tumors also enriched angiogenesis and high infiltration of vascular endothelial cells, pericytes and stromal cells and significantly enriched inflammation and immune response-related gene sets and high infiltration of CD8, CD4 memory, dendritic cells, M1 and M2 macrophages and significant elevation of cytolytic activity and immune checkpoint molecules, consistently in both cohorts. In conclusion, breast cancer patients with high apoptosis are associated with angiogenesis, immune response, high immune cell infiltration and cytolytic activity. To the best of our knowledge, this is the first study to utilize in silico translational approach to demonstrate the clinical relevance of apoptosis in breast cancer patients in large cohorts.

Cytoprotective Effect of Vitamin D on Doxorubicin-Induced Cardiac Toxicity in Triple Negative Breast Cancer

Background: Doxorubicin (Dox) is a first-line treatment for triple negative breast cancer (TNBC), but its use may be limited by its cardiotoxicity mediated by the production of reactive oxygen species. We evaluated whether vitamin D may prevent Dox-induced cardiotoxicity in a mouse TNBC model. Methods: Female Balb/c mice received rodent chow with vitamin D₃ (1500 IU/kg; vehicle) or chow supplemented with additional vitamin D₃ (total, 11,500 IU/kg). the mice were inoculated with TNBC tumors and treated with intraperitoneal Dox (6 or 10 mg/kg). Cardiac function was evaluated with transthoracic echocardiography. The cardiac tissue was evaluated with immunohistochemistry and immunoblot for levels of 4-hydroxynonenal, NAD(P)H quinone oxidoreductase (NQO1), C-MYC, and dynamin-related protein 1 (DRP1) phosphorylation. Results: At 15 to 18 days, the mean ejection fraction, stroke volume, and fractional shortening were similar between the mice treated with vitamin D + Dox (10 mg/kg) vs. vehicle but significantly greater in mice treated with vitamin D + Dox (10 mg/kg) vs. Dox (10 mg/kg). Dox (10 mg/kg) increased the cardiac tissue levels of 4-hydroxynonenal, NQO1, C-MYC, and DRP1 phosphorylation at serine 616, but these increases were not observed with vitamin D + Dox (10 mg/kg). A decreased tumor volume was observed with Dox (10 mg/kg) and vitamin D + Dox (10 mg/kg). Conclusions: Vitamin D supplementation decreased Dox-induced cardiotoxicity by decreasing the reactive oxygen species and mitochondrial damage, and did not decrease the anticancer efficacy of Dox against TNBC.

Low expression of miR-29a is associated with aggressive biology and worse survival in gastric cancer

Advanced gastric cancer (GC) is one of the most lethal cancer types, thus a better understanding of its biology in patients is urgently needed. MicroRNA (miR)-29a is a known tumor suppressive miR that is related to metastasis, but its clinical relevance in GC remains ambiguous. Here, using a large GC patient cohort we hypothesized that low expression of miR-29a in GC is associated with aggressive cancer biology and worse survival. We demonstrated that low miR-29a GC enriched cell proliferation, apoptosis, metastasis, and angiogenesis related gene sets, as well as the higher expression of related genes. Low miR-29a GC was associated with less anti-cancer immune cell infiltration as well as immune related scoring. Low miR-29a GC demonstrated a worse overall survival (OS) as well as disease specific survival (DSS) compared with high expressing miR-29a GC. Notably, low miR-29a expression was the only factor, other than residual tumor status, to be an independent prognostic biomarker of worse OS and DSS. In conclusion, low miR-29a GC was associated with aggressive cancer biology and worse OS as well as DSS. Additionally, low expression of miR-29a was an independent prognostic biomarker of OS and DSS in gastric cancer patients.

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.

Cytotoxic T-lymphocyte infiltration and chemokine predict long-term patient survival independently of tumor mutational burden in triple-negative breast cancer

Background:

Cytotoxic T-lymphocyte (CTL) infiltration into tumor is a positive prognostic factor in breast cancer. High tumor mutational burden (TMB) is also considered as a predictor of tumor immunogenicity and response to immunotherapy. However, it is unclear whether the infiltration of functional CTL simply reflects the TMB or represents an independent prognostic value.

Methods:

Utilizing The Cancer Genome Atlas (TCGA) breast cancer cohort, we established the Functional Hotness Score (FHS). The associations of FHS and breast cancer patient prognosis as well as distinct immunity markers were analyzed in a total of 3011 breast cancer patients using TCGA, METABRIC and metastatic breast cancer (MBC) cohort GSE110590.

Results:

We established FHS, based on CD8A, GZMB and CXCL10 gene expression levels of bulk tumors, which delivered the best prognostic value among some gene combinations. Breast cancer patients with the high-FHS tumors showed significantly better survival. FHS was lower in the MBCs. Triple-negative breast cancer (TNBC) showed the highest FHS among subtypes. FHS predicted patient survival in hormone receptor (HR)-negative, especially in TNBC, but not in HR-positive breast cancer. FHS predicted patient prognosis independently in TNBC. The high-FHS TNBCs showed not only higher CD8+ T cell infiltration, but also enhanced broader type-1 anti-cancer immunity. The patients with the high-FHS tumors showed better prognosis not only in high-TMB tumors but also in low-TMB TNBCs. The combination of high-TMB with high-FHS identified a unique subset of patients who do not recur over time in TNBC.

Conclusion:

TNBCs with high FHS based on the expression levels of CD8A, GZMB and CXCL10 showed improved prognosis with enhanced anti-cancer immunity regardless of TMB. FHS constitutes an independent prognostic marker of survival, particularly robustly when combined with TMB in TNBC.

Th2 cell infiltrations predict neoadjuvant chemotherapy response of estrogen receptor-positive breast cancer

Background

High infiltration of Th2 is linked to breast cancer progression and metastasis through the induction of cytokine release and T-cell anergy. The estrogen receptor (ER)-positive subtype, which accounts for 70% of breast cancer, is known to respond less to neoadjuvant chemotherapy (NAC) due to its low potential for proliferation. We hypothesized that Th2 high tumors are highly proliferative, and thus more likely to respond to NAC in ER-positive breast cancer.

Methods

We obtained clinicopathological data and overall survival information on 1,069 breast cancer patients from The Cancer Genome Atlas (TCGA). Computational algorithms and CIBERSORT were used to estimate immune cell infiltration. Additionally, xCell was used for validation.

Results

Th2 high tumors did not consistently associate with an unfavorable immune cell composition and tumor immune microenvironment but were found to be significantly elevated in the cancer stage. Th2 high tumors also correlated with high Nottingham pathological grade, as well as with Ki-67 and proliferation score in ER-positive subtypes. High Th2 tumors achieved a pathological complete response (pCR) significantly higher in ER-positive breast cancer.

Conclusions

In conclusion, high levels of Th2 are associated with aggressive features of breast cancer. Th2 levels may be a biomarker in patient selection for NAC in ER-positive breast cancer.

Inflammation Is Associated with Worse Outcome in the Whole Cohort but with Better Outcome in Triple-Negative Subtype of Breast Cancer Patients

Inflammation has been linked with cancer, but whether it is part of the problem or part of the solution remains to be a matter of debate in breast cancer. Our group and others have demonstrated that inflammation aggravates cancer progression; however, some claim that inflammation may support immune cell infiltration and suppress cancer. We defined the gene set variation analysis of the Molecular Signatures Database Hallmark inflammatory response gene set as the inflammatory pathway score and analyzed 3632 tumors in total from 4 breast cancer cohorts (METABRIC, TCGA, GSE25066, and GSE21094). In the whole breast cancer cohort, high-score tumors were associated with aggressive clinical characteristics, such as worse disease specific survival, higher Nottingham histological grade, and younger age. Inflammatory score was significantly higher in triple-negative (TNBC) as well as basal and normal subtypes compared with the other subtypes, which suggest that the detrimental effect of high level of inflammation may be because it includes a more aggressive subtype. On the contrary, high score within TNBC was significantly associated with better survival. TNBC with high score enriched not only IFN-α, IFN-γ response, IL-2/STAT5 signaling, Allograft rejection, Complement, p53 pathway, Reactive Oxygen, and Apoptosis but also TNF-α signaling, IL6-JAK-STAT signaling, TGF-β signaling, Coagulation, Angiogenesis, EMT, KRAS signaling, and PI3K-AKT-MTOR signaling gene sets. High score was associated with mainly favorable anticancerous immune cell infiltration as well as Leukocyte fraction, TIL regional fraction, Lymphocyte infiltration, IFN-γ response, TGF-β response, and cytolytic activity scores. Although the inflammatory pathway score was not associated with neoadjuvant treatment response, it associated with expressions of immune checkpoint molecules. In conclusion, inflammation was associated with worse outcome in the whole breast cancer cohort, but with better outcome in TNBC, which was associated with favorable anticancerous immune response and immune cell infiltrations.

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.

Abundance of Regulatory T Cell (Treg) as a Predictive Biomarker for Neoadjuvant Chemotherapy in Triple-Negative Breast Cancer

Regulatory CD4+ T cell (Treg), a subset of tumor-infiltrating lymphocytes (TILs), are known to suppress anticancer immunity but its clinical relevance in human breast cancer remains unclear. In this study, we estimated the relative abundance of Tregs in breast cancer of multiple patient cohorts by using the xCell algorithm on bulk tumor gene expression data. In total, 5177 breast cancer patients from five independent cohorts (TCGA-BRCA, GSE96058, GSE25066, GSE20194, and GSE110590) were analyzed. Treg abundance was not associated with cancer aggressiveness, patient survival, or immune activity markers, but it was lower in metastatic tumors when compared to matched primary tumors. Treg was associated with a high mutation rate of TP53 genes and copy number mutations as well as with increased tumor infiltration of M2 macrophages and decreased infiltration of T helper type 1 (Th1) cells. Pathological complete response (pCR) after neoadjuvant chemotherapy (NAC) was significantly associated with low Treg abundance in triple negative breast cancer (TNBC) but not in ER-positive/Her2-negative subtype. High Treg abundance was significantly associated with high tumor expression of multiple immune checkpoint inhibitor genes. In conclusion, Treg abundance may have potential as a predictive biomarker of pCR after NAC in TNBC.

M1 Macrophage and M1/M2 ratio defined by transcriptomic signatures resemble only part of their conventional clinical characteristics in breast cancer

Tumor associated macrophages (TAMs) play a critical role in biology of various cancers, including breast cancer. In the current study, we defined "M1" macrophage and "M1"/"M2" ratio by transcriptomic signatures using xCell. We investigated the association between high level of "M1" macrophage or "M1"/"M2" ratio and the tumor immune microenvironment by analyzing the transcriptome of publicly available cohorts, TCGA and METABRIC. We found that "M1" high tumors were not associated with prolonged survival compared with "M1" low tumors, or with the response to neoadjuvant chemotherapy. "M1" high tumors were associated with clinically aggressive features and "M1" high tumors enriched the cell proliferation and cell cycle related gene sets in GSEA. At the same time, "M1" high tumors were associated with high immune activity and favorable tumor immune microenvironment, as well as high expression of immune check point molecules. Strikingly, all these results were mirrored in "M1"/"M2" ratio high tumors. In conclusion, transcriptomically defined "M1" or "M1"/"M2" high tumors were associated with aggressive cancer biology and favorable tumor immune microenvironment but not with survival benefit, which resembled only part of their conventional clinical characteristics.

High G2M Pathway Score Pancreatic Cancer is Associated with Worse Survival, Particularly after Margin-Positive (R1 or R2) Resection

Pancreatic cancer is highly mortal due to uncontrolled cell proliferation. The G2M checkpoint pathway is an essential part of the cell cycle. We hypothesized that a high G2M pathway score is associated with cell proliferation and worse survival in pancreatic cancer patients. Gene set variation analysis using the Hallmark G2M checkpoint gene set was used as a score to analyze a total of 390 human pancreatic cancer patients from 3 cohorts (TCGA, GSE62452, GSE57495). High G2M score tumors enriched other cell proliferation genes sets as well as MKI67 expression, pathological grade, and proliferation score. Independent of other prognostic factors, G2M score was predictive of disease-specific survival in pancreatic cancer. High G2M tumor was associated with high mutation rate of KRAS and TP53 and significantly enriched these pathway gene sets, as well as high infiltration of Th2 cells. High G2M score consistently associated with worse overall survival in 3 cohorts, particularly in R1/2 resection, but not in R0. High G2M tumor in R1/2 highly enriched metabolic and cellular components' gene sets compared to R0. To our knowledge, this is the first study to use gene set variation analysis as a score to examine the clinical relevancy of the G2M pathway in pancreatic cancer.

ITPKC as a Prognostic and Predictive Biomarker of Neoadjuvant Chemotherapy for Triple Negative Breast Cancer

Triple negative breast cancer (TNBC) is the most aggressive subtype of breast cancer with higher mortality than the others. Pathological complete response (pCR) to neoadjuvant chemotherapy (NAC) is considered as a surrogate to predict survival. Inositol 1,4,5-trisphosphate 3-kinase C (ITPKC) is a negative regulator of T cell activation, and reduction in ITPKC function is known to promote Kawasaki disease. Given the role of tumor infiltrating lymphocytes in NAC and since TNBC has the most abundant immune cell infiltration in breast cancer, we hypothesized that the ITPKC expression level is associated with NAC response and prognosis in TNBC. The ITPKC gene was expressed in the mammary gland, but its expression was highest in breast cancer cells among other stromal cells in a bulk tumor. ITPKC expression was highest in TNBC, associated with its survival, and was its independent prognostic factor. Although high ITPKC was not associated with immune function nor with any immune cell fraction, low ITPKC significantly enriched cell proliferation-related gene sets in TNBC. TNBC with low ITPKC achieved a significantly higher pCR rate after NAC. To the best of our knowledge, this is the first report to demonstrate that ITPKC gene expression may be useful as a prognostic and predictive biomarker in TNBC.

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.

The E2F Pathway Score as a Predictive Biomarker of Response to Neoadjuvant Therapy in ER+/HER2- Breast Cancer

E2F transcription factors play critical roles in the cell cycle. Therefore, their activity is expected to reflect tumor aggressiveness and responsiveness to therapy. We scored 3905 tumors of nine breast cancer cohorts for this activity based on their gene expression for the Hallmark E2F targets gene set. As expected, tumors with a high score had an increased expression of cell proliferation-related genes. A high score was significantly associated with shorter patient survival, greater MKI67 expression, histological grade, stage, and genomic aberrations. Furthermore, metastatic tumors had higher E2F scores than the primary tumors from which they arose. Although tumors with a high score had greater infiltration by both pro- and anti-cancerous immune cells, they had an increased expression of immune checkpoint genes. Estrogen receptor (ER)-positive/human epidermal growth factor receptor 2 (HER2)-negative cancer with a high E2F score achieved a significantly higher pathological complete response (pCR) rate to neoadjuvant chemotherapy. The E2F score was significantly associated with the expression of cyclin-dependent kinase (CDK)-related genes and strongly correlated with sensitivity to CDK inhibition in cell lines. In conclusion, the E2F score is a marker of breast cancer aggressiveness and predicts the responsiveness of ER-positive/HER2-negative patients to neoadjuvant chemotherapy and possibly to CDK and immune checkpoint inhibitors.

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.

Clinical relevance of tumor microenvironment: immune cells, vessels, and mouse models

The tumor microenvironment (TME) plays a crucial role in tumor progression, therapeutic response, and patient outcomes. TME includes immune cells, blood and lymphatic vessels, and so on. There are anti-cancer and pro-cancer immune cells. In general, infiltration of anti-cancer immune cells, such as cytotoxic T cells (CTLs), is associated with a favorable patient prognosis. In contrast, infiltration of pro-cancer immune cells, such as regulatory T cells (Tregs), tumor-associated macrophages (TAMs) and myeloid-derived suppressor cells (MDSCs), is associated with a worse prognosis. However, some immune cells, which play an ambivalent role in cancer immunity, have demonstrated contradictory impacts on patient prognosis. Blood and lymphatic vessels play crucial roles in TME not only as delivery and draining systems of fluid and molecules, but also allowing cancer cells access to systematic circulation to metastasize. Angiogenesis promotes cancer aggressiveness and is associated with a worse prognosis. Its targeted therapy shows a benefit in some cancers, however, because the target can vary by caner type, a benefit of anti-angiogenesis therapy is limited in the current standard of care. Lymphangiogenesis plays a role in lymph node metastasis, thus, it is associated with a poor prognosis in some cancers. To study TME, the mouse model is one of the most commonly used tools. The choice of appropriate mouse model depends on the hypothesis being tested and the scientific question being asked. Here, we review recent studies that investigated the clinical relevance of TME components and introduce mouse models to study TME.

A Novel 4-Gene Score to Predict Survival, Distant Metastasis and Response to Neoadjuvant Therapy in Breast Cancer

We generated a 4-gene score with genes upregulated in LM2-4, a metastatic variant of MDA-MB-231 (DOK 4, HCCS, PGF, and SHCBP1) that was strongly associated with disease-free survival (DFS) in TCGA cohort (hazard ratio [HR]>1.2, p < 0.02). The 4-gene score correlated with overall survival of TCGA (HR = 1.44, p < 0.001), which was validated with DFS and disease-specific survival of METABRIC cohort. The 4-gene score was able to predict worse survival or clinically aggressive tumors, such as high Nottingham pathological grade and advanced cancer staging. High score was associated with worse survival in the hormonal receptor (HR)-positive/Her2-negative subtype. High score enriched cell proliferation-related gene sets in GSEA. The score was high in primary tumors that originated, in and metastasized to, brain and lung, and it predicted worse progression-free survival for metastatic tumors. Good tumor response to neoadjuvant chemotherapy or hormonal therapy was accompanied by score reduction. High scores were also predictive of response to neoadjuvant chemotherapy for HR-positive/Her2-negative subtype. High score tumors had increased expression of T cell exhaustion marker genes, suggesting that the score may also be a biomarker for immunotherapy response. Our novel 4-gene score with both prognostic and predictive values may, therefore, be clinically useful particularly in HR-positive breast cancer.

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.

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.

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.