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.
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.
While cancer cells gain aggressiveness by mutations, abundant mutations release neoantigens, attracting anti-cancer immune cells. We hypothesized that in breast cancer (BC), where mutation is less ...common, tumors with high mutation rates demonstrate aggressive phenotypes and attract immune cells simultaneously. High mutation rates were defined as the top 10% of the mutation rate, utilizing TCGA and METABRIC transcriptomic data. Mutation rate did not impact survival although high mutation BCs were associated with aggressive clinical features, such as more frequent in ER-negative tumors (p < 0.01), in triple-negative subtype (p = 0.03), and increased MKI-67 mRNA expression (p < 0.01) in both cohorts. Tumors with high mutation rates were associated with APOBEC3B and homologous recombination deficiency, increasing neoantigen loads (all p < 0.01). Cell proliferation and immune activity pathways were enriched in BCs with high mutation rates. Furthermore, there were higher lymphocytes and M1 macrophage infiltration in high mutation BCs. Additionally, T-cell receptor diversity, cytolytic activity score (CYT), and T-cell exhaustion marker expression were significantly elevated in BCs with high mutation rates (all p < 0.01), indicating strong immunogenicity. In conclusion, enhanced immunity due to neoantigens can be one of possible forces to counterbalance aggressiveness of a high mutation rate, resulting in similar survival rates to low mutation BCs.
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
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 (
)-α and
-γ 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.
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.
Abstract
Cancer-associated adipocytes are known to cause inflammation; however, the role of adipogenesis, the formation of adipocytes, in breast cancer is unclear. We hypothesized that intra-tumoral ...adipogenesis reflects a different cancer biology than abundance of intra-tumoral adipocytes. The Molecular Signatures Database Hallmark adipogenesis gene set of gene set variant analysis was used to quantify adipogenesis. Total of 5,098 breast cancer patients in multiple cohorts (training; GSE96058 (
n
= 3273), validation; TCGA (
n
= 1069), treatment response; GSE25066 (
n
= 508) and GSE20194 (
n
= 248)) were analyzed. Adipogenesis did not correlate with abundance of adipocytes. Adipogenesis was significantly lower in triple negative breast cancer (TNBC). Elevated adipogenesis was significantly associated with worse survival in TNBC, but not in the other subtypes. High adipogenesis TNBC was significantly associated with low homologous recombination deficiency, but not with mutation load. High adipogenesis TNBC enriched metabolism-related gene sets, but neither of cell proliferation- nor inflammation-related gene sets, which were enriched to adipocytes. High adipogenesis TNBC was infiltrated with low CD8
+
T cells and high M2 macrophages. Although adipogenesis was not associated with neoadjuvant chemotherapy response, high adipogenesis TNBC was significantly associated with low expression of
PD-L1
and
PD-L2
genes, and immune checkpoint molecules index. In conclusion, adipogenesis in TNBC was associated with cancer metabolism and unfavorable tumor immune microenvironment, which is different from abundance of adipocytes.
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.
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.
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.