Distinct gene expression profiles between primary breast cancers and brain metastases from pair-matched samples

Spatial Profiling Reveals Equivalence-Derived Molecular Signatures of Brain Mimicry and Adaptation in Breast Cancer Brain Metastases

Brain metastases (BrMets), common for advanced-stage breast cancer patients, are associated with poor median survival and accompanied by severe neurologic decline. Halting the progression of breast cancer brain metastases (BCBMs) may require modulation of the tumor microenvironment (TME), yet little is known about the impact of the primary breast TME on brain tropism, or how, once there, metastatic breast cancer cells coexist with brain-resident cells (e.g., neurons and glia). Traditionally, studies in this space have focused on differential expression analysis, overlooking potential insights gained from investigating genes with equivalent expression between groups. This is particularly crucial in distant metastasis, where tumor cells may co-opt the transcriptional programs of the host organ (e.g., brain) to facilitate successful seeding and outgrowth. Prior to our work, no computational framework existed to determine biologically-relevant equivalent gene expression. To resolve molecular mechanisms of BCBM enabled by metastatic cancer cells and/or resident brain cells, we leveraged Nanostring GeoMx to perform spatially-resolved transcriptomic profiling on 235 patient-derived tissue cores from BCBM (including adjacent normal brain), primary invasive breast cancers, and normal (non-cancer) brain; analyzing 18,677 RNAs in 450 areas of interest (AOIs). We introduce the "Equivalent Expression Index" a highly specific and accurate algorithm that identifies statistically significant "Equivalently-Expressed Genes". This method facilitated the identification of molecular remodeling and mimicry genes within tissue-specific TMEs. By integrating differential expression analysis with the Equivalent Expression Index, we discovered multiple novel gene signatures associated with BCBM and primary tumor brain-metastatic potential. We demonstrate that the Equivalent Expression Index is a powerful tool to uncover shared gene expression programs representing the adaptation of metastatic cells and brain-resident cells to the BCBM microenvironment.

Distant Metastases of Breast Cancer Resemble Primary Tumors in Cancer Cell Composition but Differ in Immune Cell Phenotypes

Breast cancer is the most commonly diagnosed cancer in women, with distant metastasis being the main cause of breast cancer-related deaths. Elucidating the changes in the tumor and immune ecosystems that are associated with metastatic disease is essential to improve understanding and ultimately treatment of metastasis. Here, we developed an in-depth, spatially resolved single-cell atlas of the phenotypic diversity of tumor and immune cells in primary human breast tumors and matched distant metastases, using imaging mass cytometry to analyze a total of 75 unique antibody targets. Although the same tumor cell phenotypes were typically present in primary tumors and metastatic sites, suggesting a strong founder effect of the primary tumor, their proportions varied between matched samples. Notably, the metastatic site did not influence tumor phenotype composition, except for the brain. Metastatic sites exhibited a lower number of immune cells overall but had a higher proportion of myeloid cells as well as exhausted and cytotoxic T cells. Myeloid cells showed distinct tissue-specific compositional signatures and increased presence of potentially matrix remodeling phenotypes in metastatic sites. This analysis of tumor and immune cell phenotypic composition of metastatic breast cancer highlights the heterogeneity of the disease within patients and across distant metastatic sites, indicating myeloid cells as the predominant immune modulators that could potentially be targeted at these sites. Significance: Multiplex imaging analysis of matched primary and metastatic breast tumors provides a phenotypic and spatial map of tumor microenvironments, revealing similar compositions of cancer cells and divergent immunologic features between matched samples.

Prognostic value of the 21-Gene Breast Recurrence Score® assay for hormone receptor-positive/human epidermal growth factor 2-negative advanced breast cancer: subanalysis from Japan Breast Cancer Research Group-M07 (FUTURE trial)

PURPOSE

This study aimed to determine whether the 21-Gene Breast Recurrence Score® assay from primary breast tissue predicts the prognosis of patients with hormone receptor-positive and human epidermal growth factor 2-negative advanced breast cancers (ABCs) treated with fulvestrant monotherapy (Group A) and the addition of palbociclib combined with fulvestrant (Group B), which included those who had progression in Group A from the Japan Breast Cancer Research Group-M07 (FUTURE trial).

METHODS

Progression-free survival (PFS) and overall survival (OS) were compared using the log-rank test and Cox regression analysis based on original recurrence score (RS) categories (Low: 0-17, Intermediate: 18-30, High: 31-100) by treatment groups (A and B) and types of ABCs (recurrence and de novo stage IV).

RESULTS

In total, 102 patients [Low: n = 44 (43.1%), Intermediate: n = 38 (37.5%), High: n = 20 (19.6%)] in Group A, and 45 in Group B, who had progression in Group A were analyzed. The median follow-up time was 23.8 months for Group A and 8.9 months for Group B. Multivariate analysis in Group A showed that low-risk [hazard ratio (HR) 0.15, 95% confidence interval (CI) 0.04-0.53, P = 0.003] and intermediate-risk (HR 0.22, 95% CI 0.06-0.78) with de novo stage IV breast cancer were significantly associated with better prognosis compared to high-risk. However, no significant difference was observed among patients with recurrence. No prognostic significance was observed in Group B.

CONCLUSION

We found a distinct prognostic value of the 21-Gene Breast Recurrence Score® assay by the types of ABCs and a poor prognostic value of the high RS for patients with de novo stage IV BC treated with fulvestrant monotherapy. Further validations of these findings are required.

An essential gene signature of breast cancer metastasis reveals targetable pathways

BACKGROUND

The differential gene expression profile of metastatic versus primary breast tumors represents an avenue for discovering new or underappreciated pathways underscoring processes of metastasis. However, as tumor biopsy samples are a mixture of cancer and non-cancer cells, most differentially expressed genes in metastases would represent confounders involving sample biopsy site rather than cancer cell biology.

METHODS

By paired analysis, we defined a top set of differentially expressed genes in breast cancer metastasis versus primary tumors using an RNA-sequencing dataset of 152 patients from The Breast International Group Aiming to Understand the Molecular Aberrations dataset (BIG-AURORA). To filter the genes higher in metastasis for genes essential for breast cancer proliferation, we incorporated CRISPR-based data from breast cancer cell lines.

RESULTS

A significant fraction of genes with higher expression in metastasis versus paired primary were essential by CRISPR. These 264 genes represented an essential signature of breast cancer metastasis. In contrast, nonessential metastasis genes largely involved tumor biopsy site. The essential signature predicted breast cancer patient outcome based on primary tumor expression patterns. Pathways underlying the essential signature included proteasome degradation, the electron transport chain, oxidative phosphorylation, and cancer metabolic reprogramming. Transcription factors MYC, MAX, HDAC3, and HCFC1 each bound significant fractions of essential genes.

CONCLUSIONS

Associations involving the essential gene signature of breast cancer metastasis indicate true biological changes intrinsic to cancer cells, with important implications for applying existing therapies or developing alternate therapeutic approaches.

Unlocking molecular mechanisms and identifying druggable targets in matched-paired brain metastasis of breast and lung cancers

Introduction

The incidence of brain metastases in cancer patients is increasing, with lung and breast cancer being the most common sources. Despite advancements in targeted therapies, the prognosis remains poor, highlighting the importance to investigate the underlying mechanisms in brain metastases. The aim of this study was to investigate the differences in the molecular mechanisms involved in brain metastasis of breast and lung cancers. In addition, we aimed to identify cancer lineage-specific druggable targets in the brain metastasis.

Methods

To that aim, a cohort of 44 FFPE tissue samples, including 22 breast cancer and 22 lung adenocarcinoma (LUAD) and their matched-paired brain metastases were collected. Targeted gene expression profiles of primary tumors were compared to their matched-paired brain metastases samples using nCounter PanCancer IO 360™ Panel of NanoString technologies. Pathway analysis was performed using gene set analysis (GSA) and gene set enrichment analysis (GSEA). The validation was performed by using Immunohistochemistry (IHC) to confirm the expression of immune checkpoint inhibitors.

Results

Our results revealed the significant upregulation of cancer-related genes in primary tumors compared to their matched-paired brain metastases (adj. p ≤ 0.05). We found that upregulated differentially expressed genes in breast cancer brain metastasis (BM-BC) and brain metastasis from lung adenocarcinoma (BM-LUAD) were associated with the metabolic stress pathway, particularly related to the glycolysis. Additionally, we found that the upregulated genes in BM-BC and BM-LUAD played roles in immune response regulation, tumor growth, and proliferation. Importantly, we identified high expression of the immune checkpoint VTCN1 in BM-BC, and VISTA, IDO1, NT5E, and HDAC3 in BM-LUAD. Validation using immunohistochemistry further supported these findings.

Conclusion

In conclusion, the findings highlight the significance of using matched-paired samples to identify cancer lineage-specific therapies that may improve brain metastasis patients outcomes.

The chromatin architectural regulator SND1 mediates metastasis in triple-negative breast cancer by promoting CDH1 gene methylation

BACKGROUND

SND1 participates in tumorigenesis, tumour invasion and metastasis in different cancers. Previous studies have shown that SND1 can promote the invasion and migration of breast cancer cells. Triple-negative breast cancer (TNBC) is a specific breast cancer subtype with high metastatic potential and poor prognosis. However, the specific roles and mechanisms of SND1 in TNBC metastasis remain unaddressed.

METHODS

Immunostaining was used to detect the SND1 expression in tissue samples of 58 TNBC and 10 glioblastomas (GBM) as positive control. The correlation between SND1 expression and patient prognosis was assessed using the Kaplan-Meier estimator. The gene expression was evaluated by qRT-PCR, Western blot and immunofluorescence analyses. Gene Ontology analysis, ChIP, a dual-luciferase reporter assay, EMSA, and 3C analysis were applied to identify SND1-activated target genes. Bisulfite sequencing PCR and MeDIP were used to detect DNA methylation. We also used wound healing, Transwell and orthotopic implantation assays to investigate the function of SND1 in TNBC cell migration and invasion.

RESULTS

The data of immunohistochemistry manifested that SND1 is the overexpression in metastasized TNBC and an independent factor for TNBC prognosis. SND1 knockdown inhibited the migration and invasion of TNBC cells. We found that SND1 promotes the metastatic phenotype of TNBC cells by epigenetically altering chromatin conformational interactions, which in turn activates DNMT3A transcription. Then, DNMT3A attenuates CCND1 expression by inducing CCND1 gene methylation, leading to TNBC metastasis.

CONCLUSION

SND1 can promote the invasion and migration of TNBC cells by promoting DNMT3A expression and suppressing CDH1 activity. SND1 is a potential biomarker and a promising therapeutic target for TNBC.

The Expression of Trace Amine-Associated Receptors (TAARs) in Breast Cancer Is Coincident with the Expression of Neuroactive Ligand-Receptor Systems and Depends on Tumor Intrinsic Subtype

Currently, the contribution of trace amine-associated receptors (TAARs) to breast cancer (BC) is recognized, but their associations with various pathological characteristics are not yet understood. There is accumulated transcriptomic data for BC tumors, which are represented in publicly accessible databases. We estimated TAARs' (including TAAR1, TAAR2, TAAR5, TAAR6, TAAR8, and TAAR9) associations with BC stage, grade, and molecular subtypes in these data and identified that the expression of all TAARs was associated with more unfavorable cancer subtypes, including basal-like and HER2-positive tumors. Also, the significant upregulation of all TAARs was demonstrated in circulating tumor cells compared to the metastatic lesions. Considering that co-expressed genes are more likely to be involved in the same biologic processes, we analyzed genes that are co-expressed with TAARs in BC. These gene sets were enriched with the genes of the olfactory transduction pathway and neuroactive ligand-receptor interaction participants. TAARs are co-expressed with G-protein-coupled receptors of monoamine neurotransmitters including dopamine, norepinephrine, and serotonin as well as with other neuroactive ligand-specific receptors. Since TAAR1 is able to modulate the activity of monoamine receptors that are involved in the regulation of BC growth, TAAR1 and potentially other TAARs may be regarded as prospective therapeutic targets for breast cancer.

Whole-Exome Sequencing Reveals High Mutational Concordance between Primary and Matched Recurrent Triple-Negative Breast Cancers

PURPOSE

Triple-negative breast cancer (TNBC) is a molecularly complex and heterogeneous breast cancer subtype with distinct biological features and clinical behavior. Although TNBC is associated with an increased risk of metastasis and recurrence, the molecular mechanisms underlying TNBC metastasis remain unclear. We performed whole-exome sequencing (WES) analysis of primary TNBC and paired recurrent tumors to investigate the genetic profile of TNBC.

METHODS

Genomic DNA extracted from 35 formalin-fixed paraffin-embedded tissue samples from 26 TNBC patients was subjected to WES. Of these, 15 were primary tumors that did not have recurrence, and 11 were primary tumors that had recurrence (nine paired primary and recurrent tumors). Tumors were analyzed for single-nucleotide variants and insertions/deletions.

RESULTS

The tumor mutational burden (TMB) was 7.6 variants/megabase in primary tumors that recurred (n = 9); 8.2 variants/megabase in corresponding recurrent tumors (n = 9); and 7.3 variants/megabase in primary tumors that did not recur (n = 15). MUC3A was the most frequently mutated gene in all groups. Mutations in MAP3K1 and MUC16 were more common in our dataset. No alterations in PI3KCA were detected in our dataset.

CONCLUSIONS

We found similar mutational profiles between primary and paired recurrent tumors, suggesting that genomic features may be retained during local recurrence.

Type I interferon response in astrocytes promotes brain metastasis by enhancing monocytic myeloid cell recruitment

Cancer metastasis to the brain is a significant clinical problem. Metastasis is the consequence of favorable interactions between invaded cancer cells and the microenvironment. Here, we demonstrate that cancer-activated astrocytes create a sustained low-level activated type I interferon (IFN) microenvironment in brain metastatic lesions. We further confirm that the IFN response in astrocytes facilitates brain metastasis. Mechanistically, IFN signaling in astrocytes activates C-C Motif Chemokine Ligand 2 (CCL2) production, which further increases the recruitment of monocytic myeloid cells. The correlation between CCL2 and monocytic myeloid cells is confirmed in clinical brain metastasis samples. Lastly, genetically or pharmacologically inhibiting C-C Motif Chemokine Receptor 2 (CCR2) reduces brain metastases. Our study clarifies a pro-metastatic effect of type I IFN in the brain even though IFN response has been considered to have anti-tumor effects. Moreover, this work expands our understandings on the interactions between cancer-activated astrocytes and immune cells in brain metastasis.

Oral IRAK-4 Inhibitor CA-4948 Is Blood-Brain Barrier Penetrant and Has Single-Agent Activity against CNS Lymphoma and Melanoma Brain Metastases

PURPOSE

An ongoing challenge in cancer is the management of primary and metastatic brain malignancies. This is partly due to restrictions of the blood-brain barrier and their unique microenvironment. These challenges are most evident in cancers such as lymphoma and melanoma, which are typically responsive to treatment in systemic locations but resistant when established in the brain. We propose interleukin-1 receptor-associated kinase-4 (IRAK-4) as a potential target across these diseases and describe the activity and mechanism of oral IRAK-4 inhibitor CA-4948.

EXPERIMENTAL DESIGN

Human primary central nervous system lymphoma (PCNSL) and melanoma brain metastases (MBM) samples were analyzed for expression of IRAK-4 and downstream transcription pathways. We next determined the central nervous system (CNS) applicability of CA-4948 in naïve and tumor-bearing mice using models of PCNSL and MBM. The mechanistic effect on tumors and the tumor microenvironment was then analyzed.

RESULTS

Human PCNSL and MBM have high expression of IRAK-4, IRAK-1, and nuclear factor kappa B (NF-κB). This increase in inflammation results in reflexive inhibitory signaling. Similar profiles are observed in immunocompetent murine models. Treatment of tumor-bearing animals with CA-4948 results in the downregulation of mitogen-activated protein kinase (MAPK) signaling in addition to decreased NF-κB. These intracellular changes are associated with a survival advantage.

CONCLUSIONS

IRAK-4 is an attractive target in PCNSL and MBM. The inhibition of IRAK-4 with CA-4948 downregulates the expression of important transcription factors involved in tumor growth and proliferation. CA-4948 is currently being investigated in clinical trials for relapsed and refractory lymphoma and warrants further translation into PCNSL and MBM.

Pan-cancer molecular subtypes of metastasis reveal distinct and evolving transcriptional programs

Molecular mechanisms underlying cancer metastasis span diverse tissues of origin. Here, we synthesize and collate the transcriptomes of patient-derived xenografts and patient tumor metastases, and these data collectively represent 38 studies and over 3,000 patients and 4,000 tumors. We identify four expression-based subtypes of metastasis transcending tumor lineage. The first subtype has extensive copy alterations, higher expression of MYC transcriptional targets and DNA repair genes, and bromodomain inhibitor response association. The second subtype has higher expression of genes involving metabolism and prostaglandin synthesis and regulation. The third subtype has evidence of neuronal differentiation, higher expression of DNA and histone methylation genes and EZH2 transcriptional targets, and BCL2 inhibitor response association. The fourth subtype has higher expression of immune checkpoint and Notch pathway genes. The metastasis subtypes reflect expression differences from paired primaries, with subtype switching being common. These subtypes facilitate understanding of the molecular underpinnings of metastases beyond tissue-oriented domains, with therapeutic implications.

An integrated investigation of structural and pathway alteration caused by PIK3CA and TP53 mutations identified in cfDNA of metastatic breast cancer

In peripheral blood, cell-free DNA (cfDNA) contains circulating tumor DNA (ctDNA), which indicates molecular abnormalities in metastatic breast tumor tissue. The sequencing of cfDNA of Metastatic Breast Cancer (MBC) patients allows assessment of therapy response and noninvasive treatment. In the proposed study, clinically significant alterations in PIK3CA and TP53 genes associated with MBC resulting in a missense substitution of His1047Arg and Arg282Trp from an next-generation sequencing-based multi-gene panel were reported in a cfDNA of a patient with MBC. To investigate the impact of the reported mutation, we used molecular docking, molecular dynamics simulation, network analysis, and pathway analysis. Molecular Docking analysis determined the distinct binding pattern revealing H1047R-ATP complex has a higher number of Hydrogen bonds (H-bonds) and binding affinity with a slight difference compared to the PIK3CA-ATP complex. Following, molecular dynamics simulation for 200 ns, of which H1047R-ATP complex resulted in the instability of PIK3CA. Similarly, for TP53 mutant R282W, the zinc-free state (apo) and zinc-bounded (holo) complexes were investigated for conformational change between apo and holo complexes, of which the holo complex mutant R282W was unstable. To validate the conformational change of PIK3CA and TP53, 80% mutation of H1047R in the kinase domain of p110α expressed ubiquitously in PIK3CA protein that alters PI3K pathway, while R282W mutation in DNA binding helix (H2) region of P53 protein inhibits the transcription factor in P53 pathway causing MBC. According to our findings, the extrinsic (hypoxia, oxidative stress, and acidosis); intrinsic factors (MYC amplification) in PIK3CA and TP53 mutations will provide potential insights for developing novel therapeutic methods for MBC therapy.

Distinct tumor architectures for metastatic colonization of the brain

Brain metastasis is a dismal cancer complication, hinging on the initial survival and outgrowth of disseminated cancer cells. To understand these crucial early stages of colonization, we investigated two prevalent sources of cerebral relapse, triple-negative (TNBC) and HER2+ breast cancer (HER2BC). We show that these tumor types colonize the brain aggressively, yet with distinct tumor architectures, stromal interfaces, and autocrine growth programs. TNBC forms perivascular sheaths with diffusive contact with astrocytes and microglia. In contrast, HER2BC forms compact spheroids prompted by autonomous extracellular matrix components and segregating stromal cells to their periphery. Single-cell transcriptomic dissection reveals canonical Alzheimer's disease-associated microglia (DAM) responses. Differential engagement of tumor-DAM signaling through the receptor AXL suggests specific pro-metastatic functions of the tumor architecture in both TNBC perivascular and HER2BC spheroidal colonies. The distinct spatial features of these two highly efficient modes of brain colonization have relevance for leveraging the stroma to treat brain metastasis.

Imputation of single-cell transcriptome data enables the reconstruction of networks predictive of breast cancer metastasis

Single-cell transcriptome data provide a unique opportunity to explore the gene networks of a particular cell type. However, insufficient capture rate and high dimensionality of single-cell RNA sequencing (scRNA-seq) data challenge cell-type-specific gene network (CGN) reconstruction. Here, we demonstrated that the imputation of scRNA-seq data enables reconstruction of CGNs by effective retrieval of gene functional associations. We reconstructed CGNs for seven primary and nine metastatic breast cancer cell lines using scRNA-seq data with imputation. Key genes for primary or metastatic cell lines were prioritized based on network centrality measures and CGN hub genes that were presumed to be the major determinant of cell type characteristics. To identify novel genes in breast cancer metastasis, we used the average rank difference of centrality between the primary and metastatic cell lines. Genes predicted using CGN centrality analysis were more enriched for known breast cancer metastatic genes than those predicted using differential expression. The molecular chaperone CCT2 was identified as a novel gene for breast metastasis during knockdown assays of several candidate genes. Overall, our study demonstrated an effective CGN reconstruction technique with imputation of scRNA-seq data and the feasibility of identifying key genes for particular cell subsets using single-cell network analysis.

Regulation of micropatterned curvature-dependent FA heterogeneity on cytoskeleton tension and nuclear DNA synthesis of malignant breast cancer cells

Breast cancer is considered as a worldwide disease due to its high incidence and malignant metastasis. Although numerous techniques have been developed well to conduct breast cancer therapy, the influence of micropattern-induced interfacial heterogeneity on the molecular mechanism and nuclear signalling transduction of carcinogenesis is rarely announced. In this study, PDMS stencil-assisted micropatterns were fabricated on tissue culture plates to manage cell clustering colony by adjusting initial cell seeding density and the size of microholes. The curvature of each microholes was controlled to construct the interfacial heterogeneity of MDA-MB231 cancer cells at the periphery of micropatterned colony. The distinguished focal adhesion (FA) and cytoskeleton distribution at the central and peripheral regions of the cell colony were regulated by heterogeneous properties. The interfacial heterogeneity of FA and cytoskeleton would induce the biased tension force to encourage more ezrin expression at the periphery and further promote DNA synthesis, therefore disclosing a stem-like phenotype in heterogeneous cells. This study will provide a value source of information for the development of micropattern-induced heterogeneity and the interpretation of metastatic mechanism in malignant breast cancer cells.

Re-purposing the pro-senescence properties of doxorubicin to introduce immunotherapy in breast cancer brain metastasis

An increasing number of breast cancer patients develop brain metastases (BM). Standard-of-care treatments are largely inefficient, and breast cancer brain metastasis (BCBM) patients are considered untreatable. Immunotherapies are not successfully employed in BCBM, in part because breast cancer is a "cold" tumor and also because the brain tissue has a unique immune landscape. Here, we generate and characterize immunocompetent models of BCBM derived from PyMT and Neu mammary tumors to test how harnessing the pro-senescence properties of doxorubicin can be used to prime the specific immune BCBM microenvironment. We reveal that BCBM senescent cells, induced by doxorubicin, trigger the recruitment of PD1-expressing T cells to the brain. Importantly, we demonstrate that induction of senescence with doxorubicin improves the efficacy of immunotherapy with anti-PD1 in BCBM in a CD8 T cell-dependent manner, thereby providing an optimized strategy to introduce immune-based treatments in this lethal disease. In addition, our BCBM models can be used for pre-clinical testing of other therapeutic strategies in the future.

Decellularized organ biomatrices facilitate quantifiable in vitro 3D cancer metastasis models

Metastatic cancers are chemoresistant, involving complex interplay between disseminated cancer cell aggregates and the distant organ microenvironment (extracellular matrix and stromal cells). Conventional metastasis surrogates (scratch/wound healing, Transwell migration assays) lack 3D architecture and ECM presence. Metastasis studies can therefore significantly benefit from biomimetic 3D in vitro models recapitulating the complex cascade of distant organ invasion and colonization by collective clusters of cells. We aimed to engineer reproducible and quantifiable 3D models of highly therapy-resistant cancer processes: (i) colorectal cancer liver metastasis; and (ii) breast cancer lung metastasis. Metastatic seeds are engineered using 3D tumor spheroids to recapitulate the 3D aggregation of cancer cells both in the tumor and in circulation throughout the metastatic cascade of many cancers. Metastatic soil was engineered by decellularizing porcine livers and lungs to generate biomatrix scaffolds, followed by extensive materials characterization. HCT116 colorectal and MDA-MB-231 breast cancer spheroids were generated on hanging drop arrays to initiate clustered metastatic seeding into liver and lung biomatrix scaffolds, respectively. Between days 3-7, biomatrix cellular colonization was apparent with increased metabolic activity and the presence of cellular nests evaluated via multiphoton microscopy. HCT116 and MDA-MB-231 cells colonized liver and lung biomatrices, and at least 15% of the cells invaded more than 20 μm from the surface. Engineered metastases also expressed increased signatures of genes associated with the metastatic epithelial to mesenchymal transition (EMT). Importantly, inhibition of matrix metalloproteinase-9 inhibited metastatic invasion into the biomatrix. Furthermore, metastatic nests were significantly more chemoresistant (>3 times) to the anti-cancer drug oxaliplatin, compared to 3D spheroids. Together, our data indicated that HCT116 and MDA-MB-231 spheroids invade, colonize, and proliferate in livers and lungs establishing metastatic nests in 3D settings in vitro. The metastatic nature of these cells was confirmed with functional readouts regarding EMT and chemoresistance. Modeling the dynamic metastatic cascade in vitro has potential to identify therapeutic targets to treat or prevent metastatic progression in chemoresistant metastatic cancers.

CILP, a Putative Gene Associated With Immune Infiltration in Breast Cancer Brain Metastases

Breast cancer (BC) is the second leading cause of brain metastases (BM), with high morbidity and mortality. The aim of our study was to explore the effect of the cartilage intermediate layer protein (CILP) on breast cancer brain metastases (BCBM). Using a weighted gene coexpression network analysis (WGCNA) in GSE100534 and GSE125989 datasets, we found that the yellow module was closely related to the occurrence of BCBM, and CILP was a hub gene in the yellow module. Low CILP expression was associated with a poor prognosis, and it was an independent prognostic factor for stage III-IV BC determined using Cox regression analysis. A nomogram model including CILP expression was established to predict the 5-, 7-, and 10-year overall survival (OS) probabilities of stage III-IV BC patients. We found that CILP mRNA expression was downregulated in BCBM through GSE100534, GSE125989, and GSE43837 datasets. In addition, we found that CILP mRNA expression was negatively correlated with vascular endothelial growth factor A (VEGFA), which is involved in regulating the development of BM. UALCAN analysis showed that CILP expression was downregulated in HER2-positive (HER2+) and triple-negative breast cancer (TNBC), which are more prone to BM. The vitro experiments demonstrated that CILP significantly inhibited BC cell proliferation and metastasis. Western blot (WB) results further showed that the mesenchymal protein marker vimentin was significantly downregulated following CILP overexpression, suggesting that CILP could participate in migration through epithelial-mesenchymal transition (EMT). A comparison of CILP expression using immunohistochemistry in BC and BCBM showed that CILP was significantly downregulated in BCBM. In addition, gene set variation analysis (GSVA) revealed that CILP was associated with the T-cell receptor signaling pathway in BCBM and BC, indicating that CILP may be involved in BCBM through immune effects. BCBM showed lower immune infiltration than BC. Moreover, CILP expression was positively correlated with HLA-II, T helper cells (CD4⁺ T cells), and Type II IFN Response in BCBM. Collectively, our study indicates that CILP is associated with immune infiltration and may be a putative gene involved in BCBM. CILP offers new insights into the pathogenesis of BCBM, which will facilitate the development of novel targets for BCBM patients.

DNMT3A and DNMT3B in Breast Tumorigenesis and Potential Therapy

Breast cancer has become a leading cause of cancer-related deaths in women worldwide. DNA methylation has been revealed to play an enormously important role in the development and progression of breast cancer. DNA methylation is regulated by DNA methyltransferases (DNMTs), including DNMT1, DNMT2, and DNMT3. DNMT3 family has three members: DNMT3A, DNMT3B, and DNMT3L. The roles and functions of DNMT1 in breast cancer have been well reviewed. In this article, the roles of DNMT3A and DNMT3B in breast tumorigenesis and development are reviewed. We also discuss the SNP and mutations of DNMT3A and DNMT3B in breast cancer. In addition, we summarize how DNMT3A and DNMT3B are regulated by non-coding RNAs and signaling pathways in breast cancer, and targeting the expression levels of DNMT3A and DNMT3B may be a promising therapeutic approach for breast cancer. This review will provide reference for further studies on the biological functions and molecular mechanisms of DNMT3A and DNMT3B in breast cancer.

Distinct Gene Expression Profiles of Matched Primary and Metastatic Triple-Negative Breast Cancers

Background: Although triple-negative breast cancer (TNBC) is associated with an increased risk of recurrence and metastasis, the molecular mechanisms underlying metastasis in TNBC remain unknown. To identify transcriptional changes and genes regulating metastatic progression in TNBC, we compared the transcriptomic profiles of primary and matched metastatic tumors using massively parallel RNA sequencing. Methods: We performed gene expression profiling using formalin-fixed paraffin-embedded (FFPE) TNBC tissues of patients from two cohorts: the Zurich cohort (n = 31) and the Stavanger cohort (n = 5). Among the 31 patients in the Zurich cohort, 18 had primary TNBC tumors that did not metastasize, and 13 had primary tumors that metastasized (11 paired primary and locoregional recurrences). The Stavanger cohort included five matched primary and metastatic TNBC tumors. Significantly differentially expressed genes (DEGs; absolute fold change ≥2, p < 0.05) were identified and subjected to functional analyses. We investigated if there was any overlap between DEGs from both the cohorts with epithelial-to-mesenchymal-to-amoeboid transition (EMAT) gene signature. xCell was used to estimate relative fractions of 64 immune and stromal cell types in each RNA-seq sample. Results: In the Zurich cohort, we identified 1624 DEGs between primary TNBC tumors and matched metastatic lesions. xCell analysis revealed a significantly higher immune scores for metastatic lesions compared to paired primary tumors in the Zurich cohort. We also found significant upregulation of three MammaPrint signature genes (HRASLS, TGFB3 and RASSF7) in primary tumors that metastasized compared to primary tumors that remained metastasis-free. In the Stavanger cohort, we identified 818 DEGs between primary tumors and matched metastatic lesions. No significant differences in xCell immune scores were observed. We found that 21 and 14 DEGs from Zurich and Stavanger cohort, respectively, overlapped with the EMAT gene signature. In both cohorts, genes belonging to the MMP, FGF, and PDGFR families were upregulated in primary tumors compared to matched metastatic lesions. Conclusions: Our results suggest that distinct gene expression patterns exist between primary TNBCs and matched metastatic tumors. Further studies are warranted to explore whether these discrete expression profiles underlie or result from disease status.

The double agents in liquid biopsy: promoter and informant biomarkers of early metastases in breast cancer

Breast cancer continues to be a major global problem with significant mortality associated with advanced stage and metastases at clinical presentation. However, several findings suggest that metastasis is indeed an early occurrence. The standard diagnostic techniques such as invasive core needle biopsy, serological protein marker assays, and non-invasive radiological imaging do not provide information about the presence and molecular profile of small fractions of early metastatic tumor cells which are prematurely dispersed in the circulatory system. These circulating tumor cells (CTCs) diverge from the primary tumors as clusters with a defined secretome comprised of circulating cell-free nucleic acids and small microRNAs (miRNAs). These circulatory biomarkers provide a blueprint of the mutational profile of the tumor burden and tumor associated alterations in the molecular signaling pathways involved in oncogenesis. Amidst the multitude of circulatory biomarkers, miRNAs serve as relatively stable and precise biomarkers in the blood for the early detection of CTCs, and promote step-wise disease progression by executing paracrine signaling that transforms the microenvironment to guide the metastatic CTCs to anchor at a conducive new organ. Random sampling of easily accessible patient blood or its serum/plasma derivatives and other bodily fluids collectively known as liquid biopsy (LB), forms an efficient alternative to tissue biopsies. In this review, we discuss in detail the divergence of early metastases as CTCs and the involvement of miRNAs as detectable blood-based diagnostic biomarkers that warrant a timely screening of cancer, serial monitoring of therapeutic response, and the dynamic molecular adaptations induced by miRNAs on CTCs in guiding primary and second-line systemic therapy.

Altered propionate metabolism contributes to tumour progression and aggressiveness

The alteration of metabolic pathways is a critical strategy for cancer cells to attain the traits necessary for metastasis in disease progression. Here, we find that dysregulation of propionate metabolism produces a pro-aggressive signature in breast and lung cancer cells, increasing their metastatic potential. This occurs through the downregulation of methylmalonyl coenzyme A epimerase (MCEE), mediated by an extracellular signal-regulated kinase 2-driven transcription factor Sp1/early growth response protein 1 transcriptional switch driven by metastatic signalling at its promoter level. The loss of MCEE results in reduced propionate-driven anaplerotic flux and intracellular and intratumoral accumulation of methylmalonic acid, a by-product of propionate metabolism that promotes cancer cell invasiveness. Altogether, we present a previously uncharacterized dysregulation of propionate metabolism as an important contributor to cancer and a valuable potential target in the therapeutic treatment of metastatic carcinomas.

Intracranial Metastatic Disease: Present Challenges, Future Opportunities

Intracranial metastatic disease (IMD) is a prevalent complication of cancer that significantly limits patient survival and quality of life. Over the past half-century, our understanding of the epidemiology and pathogenesis of IMD has improved and enabled the development of surveillance and treatment algorithms based on prognostic factors and tumor biomolecular characteristics. In addition to advances in surgical resection and radiation therapy, the treatment of IMD has evolved to include monoclonal antibodies and small molecule antagonists of tumor-promoting proteins or endogenous immune checkpoint inhibitors. Moreover, improvements in the sensitivity and specificity of imaging as well as the development of new serological assays to detect brain metastases promise to revolutionize IMD diagnosis. In this review, we will explore current treatment principles in patients with IMD, including the emerging role of targeted and immunotherapy in select primary cancers, and discuss potential areas for further investigation.

NOTCH induced MDSC recruitment after oHSV virotherapy in CNS cancer models modulates anti-tumor immunotherapy

PURPOSE

Oncolytic herpes simplex virus-1 (oHSV) infection of brain tumors activates NOTCH, however the consequences of NOTCH on oHSV-induced immunotherapy is largely unknown. Here we evaluated the impact of NOTCH blockade on virus-induced immunotherapy.

EXPERIMENTAL DESIGN

RNA sequencing (RNA-seq), TCGA data analysis, flow cytometry, Luminex- and ELISA-based assays, brain tumor animal models, and serum analysis of patients with recurrent glioblastoma (GBM) treated with oHSV was used to evaluate the effect of NOTCH signaling on virus-induced immunotherapy.

RESULTS

TCGA data analysis of patients with grade IV glioma and oHSV treatment of experimental brain tumors in mice showed that NOTCH signaling significantly correlated with a higher myeloid cell infiltration. Immunofluorescence staining and RNA-seq uncovered a significant induction of Jag1 (NOTCH ligand) expression in infiltrating myeloid cells upon oHSV infection. Jag1-expressing macrophages further spread NOTCH activation in the tumor microenvironment (TME). NOTCH-activated macrophages increased the secretion of CCL2, which further amplified myeloid-derived suppressor cells. CCL2 and IL10 induction was also observed in serum of patients with recurrent GBM treated with oHSV (rQnestin34.5; NCT03152318). Pharmacologic blockade of NOTCH signaling rescued the oHSV-induced immunosuppressive TME and activated a CD8-dependent antitumor memory response, resulting in a therapeutic benefit.

CONCLUSIONS

NOTCH-induced immunosuppressive myeloid cell recruitment limited antitumor immunity. Translationally, these findings support the use of NOTCH inhibition in conjunction with oHSV therapy.

Discordance of PIK3CA and TP53 mutations between breast cancer brain metastases and matched primary tumors

There is limited knowledge of the biology of breast cancer (BC) brain metastasis (BM). We primarily aimed to determine the mutations in BCBM and to compare the mutational pattern with the matched primary breast cancer (BC). Secondary aims were to determine mutations in each subgroup (Luminal A-/B-like, HER2+ and TNBC) of BCBM, and to determine survival according to specific mutations. We investigated 57 BCBMs, including 46 cases with matched primary tumors (PT) by targeted Next Generation Sequencing (NGS) using the Cancer Hotspot Panel v2 (ThermoFisher Scientific) covering 207 targeted regions in 50 cancer related genes. Subtype according to immunohistochemistry was re-evaluated. NGS results fulfilling sequencing quality criteria were obtained from 52 BM and 41 PT, out of which 37 were matched pairs. Pathogenic mutations were detected in 66% of PTs (27/41), and 62% of BMs (32/52). TP53 mutations were most frequent; 49% (20/41) of PTs and 48% (25/52) in BMs, followed by PIK3CA mutations; 22% (9/42) in PTs and 25% (13/52) in BMs. Mutations in CDH1, EGFR, HRAS, RB1 CDKN2A and PTEN were detected in single pairs or single samples. Mutational pattern was discordant in 24% of matched pairs. We show a discordance of PIK3CA and TP53 mutations of roughly 25% indicating the need to develop methods to assess mutational status in brain metastasis where analysis of cell-free DNA from cerebrospinal fluid (CSF) has shown promising results.

Identification of potential genes related to breast cancer brain metastasis in breast cancer patients

Brain metastases (BMs) usually develop in breast cancer (BC) patients. Thus, the molecular mechanisms of breast cancer brain metastasis (BCBM) are of great importance in designing therapeutic strategies to treat or prevent BCBM. The present study attempted to identify novel diagnostic and prognostic biomarkers of BCBM. Two datasets (GSE125989 and GSE100534) were obtained from the Gene Expression Omnibus (GEO) database to find differentially expressed genes (DEGs) in cases of BC with and without brain metastasis (BM). A total of 146 overlapping DEGs, including 103 up-regulated and 43 down-regulated genes, were identified. Functional enrichment analysis showed that these DEGs were mainly enriched for functions including extracellular matrix (ECM) organization and collagen catabolic fibril organization. Using protein-protein interaction (PPI) and principal component analysis (PCA) analysis, we identified ten key genes, including LAMA4, COL1A1, COL5A2, COL3A1, COL4A1, COL5A1, COL5A3, COL6A3, COL6A2, and COL6A1. Additionally, COL5A1, COL4A1, COL1A1, COL6A1, COL6A2, and COL6A3 were significantly associated with the overall survival of BC patients. Furthermore, COL6A3, COL5A1, and COL4A1 were potentially correlated with BCBM in human epidermal growth factor 2 (HER2) expression. Additionally, the miR-29 family might participate in the process of metastasis by modulating the cancer microenvironment. Based on datasets in the GEO database, several DEGs have been identified as playing potentially important roles in BCBM in BC patients.

A Novel Brain-Permeant Chemotherapeutic Agent for the Treatment of Brain Metastasis in Triple-Negative Breast Cancer

Development of metastases to central nervous system (CNS) is an increasing clinical issue following the diagnosis of advanced breast cancer. The propensity to metastasize to CNS varies by breast cancer subtype. Of the four breast cancer subtypes, triple-negative breast cancers (TNBC) have the highest rates of both parenchymal brain metastasis and leptomeningeal metastasis (LM). LM is rapidly fatal due to poor detection and limited therapeutic options. Therapy of TNBC brain metastasis and LM is challenged by multifocal brain metastasis and diffuse spread of LM, and must balance brain penetration, tumor cytotoxicity, and the avoidance of neurotoxicity. Thus, there is an urgent need for novel therapeutic options in TNBCs CNS metastasis. QBS10072S is a novel chemotherapeutic that leverages TNBC-specific defects in DNA repair and LAT1 (L-amino acid transporter type 1)-dependent transport into the brain. In our study, activity of QBS10072S was investigated in vitro with various cell lines including the human TNBC cell line MDA-MB-231 and its brain-tropic derivative MDA-MB-231-BR3. QBS10072S was preferentially toxic to TNBC cells. The efficacy of QBS10072S against brain metastasis and LM was tested using a model of brain metastasis based on the internal carotid injection of luciferase-expressing tumor cells into NuNu mice. The compound was well tolerated, delayed tumor growth and reduced leptomeningeal dissemination, resulting in significant extension of survival. Given that current treatments for LM are palliative with only few studies reporting a survival benefit, QBS10072S is planned to be investigated in clinical trials as a therapeutic for TNBC LM. SIGNIFICANCE: TNBC brain metastasis often involves dissemination into leptomeninges. Treatment options for TNBC leptomeningeal metastasis are limited and are mostly palliative. Our study demonstrates significant efficacy of the brain-penetrating agent QBS10072S against TNBC brain metastasis and leptomeningeal spread.

CancerSiamese: one-shot learning for predicting primary and metastatic tumor types unseen during model training

BACKGROUND

The state-of-the-art deep learning based cancer type prediction can only predict cancer types whose samples are available during the training where the sample size is commonly large. In this paper, we consider how to utilize the existing training samples to predict cancer types unseen during the training. We hypothesize the existence of a set of type-agnostic expression representations that define the similarity/dissimilarity between samples of the same/different types and propose a novel one-shot learning model called CancerSiamese to learn this common representation. CancerSiamese accepts a pair of query and support samples (gene expression profiles) and learns the representation of similar or dissimilar cancer types through two parallel convolutional neural networks joined by a similarity function.

RESULTS

We trained CancerSiamese for cancer type prediction for primary and metastatic tumors using samples from the Cancer Genome Atlas (TCGA) and MET500. Network transfer learning was utilized to facilitate the training of the CancerSiamese models. CancerSiamese was tested for different N-way predictions and yielded an average accuracy improvement of 8% and 4% over the benchmark 1-Nearest Neighbor (1-NN) classifier for primary and metastatic tumors, respectively. Moreover, we applied the guided gradient saliency map and feature selection to CancerSiamese to examine 100 and 200 top marker-gene candidates for the prediction of primary and metastatic cancers, respectively. Functional analysis of these marker genes revealed several cancer related functions between primary and metastatic tumors.

CONCLUSION

This work demonstrated, for the first time, the feasibility of predicting unseen cancer types whose samples are limited. Thus, it could inspire new and ingenious applications of one-shot and few-shot learning solutions for improving cancer diagnosis, prognostic, and our understanding of cancer.

Diseasome and comorbidities complexities of SARS-CoV-2 infection with common malignant diseases

With the increasing number of immunoinflammatory complexities, cancer patients have a higher risk of serious disease outcomes and mortality with SARS-CoV-2 infection which is still not clear. In this study, we aimed to identify infectome, diseasome and comorbidities between COVID-19 and cancer via comprehensive bioinformatics analysis to identify the synergistic severity of the cancer patient for SARS-CoV-2 infection. We utilized transcriptomic datasets of SARS-CoV-2 and different cancers from Gene Expression Omnibus and Array Express Database to develop a bioinformatics pipeline and software tools to analyze a large set of transcriptomic data and identify the pathobiological relationships between the disease conditions. Our bioinformatics approach revealed commonly dysregulated genes (MARCO, VCAN, ACTB, LGALS1, HMOX1, TIMP1, OAS2, GAPDH, MSH3, FN1, NPC2, JUND, CHI3L1, GPNMB, SYTL2, CASP1, S100A8, MYO10, IGFBP3, APCDD1, COL6A3, FABP5, PRDX3, CLEC1B, DDIT4, CXCL10 and CXCL8), common gene ontology (GO), molecular pathways between SARS-CoV-2 infections and cancers. This work also shows the synergistic complexities of SARS-CoV-2 infections for cancer patients through the gene set enrichment and semantic similarity. These results highlighted the immune systems, cell activation and cytokine production GO pathways that were observed in SARS-CoV-2 infections as well as breast, lungs, colon, kidney and thyroid cancers. This work also revealed ribosome biogenesis, wnt signaling pathway, ribosome, chemokine and cytokine pathways that are commonly deregulated in cancers and COVID-19. Thus, our bioinformatics approach and tools revealed interconnections in terms of significant genes, GO, pathways between SARS-CoV-2 infections and malignant tumors.

An immune-related gene expression signature predicts brain metastasis in lung adenocarcinoma patients after surgery: gene expression profile and immunohistochemical analyses

Background

Lung adenocarcinoma (LUAD) with brain metastasis (BM) occurs frequently and has a poor prognosis. In this study, we aimed to assess the correlation between gene expression signatures and the development of BM after surgical resection of LUAD.

Methods

We analyzed the immune-related gene expression profiles of 72 LUADs with and without BM after surgery and verified them using NanoString method and immunohistochemistry (IHC). We matched the Tumor, Node, Metastasis (TNM) stage in the groups with and without BM to minimize the effect of TNM stage. Pathway enrichment studies were also performed.

Results

In the NanoString results, we identified 11 upregulated immune-related gene signature that correlated specifically with BM in the discovery and validation sets [area under the curve (AUC) =0.750 and 0.787, respectively]. The discovery set achieved 94% sensitivity and 62% specificity and the validation set displayed 100% sensitivity and 50% specificity. Eight out of the 11 genes were verified by IHC and had profiles similar to the gene expression profile results (AUC =0.844 for the discovery set and AUC =0.795 for the validation set). Subgroup analysis revealed that 11 immune-related gene signature enabled prediction of BM at all TNM stages. There were no differences in the 11 immune-related gene expression signatures between the primary LUAD samples and the matched brain samples. Pathway enrichment analysis revealed that the cytokine-cytokine receptor interaction pathway was closely correlated with BM.

Conclusions

The 11 identified immune-related gene expression signatures may be potentially clinically useful predictors for BM and can provide patient-specific treatment options.

Metastasis-Specific Gene Expression in Autochthonous and Allograft Mouse Mammary Tumor Models: Stratification and Identification of Targetable Signatures

Breast cancer metastasis is a leading cause of cancer-related death of women in the United States. A hurdle in advancing metastasis-targeted intervention is the phenotypic heterogeneity between primary and secondary lesions. To identify metastasis-specific gene expression profiles we performed RNA-sequencing of breast cancer mouse models; analyzing metastases from models of various drivers and routes. We contrasted the models and identified common, targetable signatures. Allograft models exhibited more mesenchymal-like gene expression than genetically engineered mouse models (GEMM), and primary culturing of GEMM-derived metastatic tissue induced mesenchymal-like gene expression. In addition, metastasis-specific transcriptomes differed between tail vein and orthotopic injection of the same cell line. Gene expression common to models of spontaneous metastasis included sildenafil response and nicotine degradation pathways. Strikingly, in vivo sildenafil treatment significantly reduced metastasis by 54%, while nicotine significantly increased metastasis by 46%. These data suggest that (i) actionable metastasis-specific pathways can be readily identified, (ii) already available drugs may have great potential to alleviate metastatic incidence, and (iii) metastasis may be influenced greatly by lifestyle choices such as the choice to consume nicotine products. In summary, while mouse models of breast cancer metastasis vary in ways that must not be ignored, there are shared features that can be identified and potentially targeted therapeutically. IMPLICATIONS: The data we present here exposes critical variances between preclinical models of metastatic breast cancer and identifies targetable pathways integral to metastatic spread. VISUAL OVERVIEW: http://mcr.aacrjournals.org/content/molcanres/18/9/1278/F1.large.jpg.

Genomic landscape of the immune microenvironments of brain metastases in breast cancer

Background

This study was intended to investigate the genomic landscape of the immune microenvironments of brain metastases in breast cancer.

Methods

Three gene expression profile datasets (GSE76714, GSE125989 and GSE43837) of breast cancer with brain metastases were downloaded from Gene Expression Omnibus (GEO) database. After differential expression analysis, the tumor immune microenvironment and immune cell infiltration were analyzed. Then immune-related genes were identified, followed by function analysis, transcription factor (TF)-miRNA–mRNA co-regulatory network analysis, and survival analysis of metastatic recurrence.

Results

The present results showed that the tumor immune microenvironment in brain metastases was immunosuppressed compared with primary caner. Compared with primary cancer samples, the infiltration ratio of plasma cells in brain metastases samples was significantly higher, while the infiltration ratio of macrophages M2 cells in brain metastases samples was significantly lower. Total 42 immune-related genes were identified, such as THY1 and NEU2. CD1B, THY1 and DOCK2 were found to be implicated in the metastatic recurrence of breast cancer.

Conclusions

Targeting macrophages or plasma cells may be new strategies for immunotherapy of breast cancer with brain metastases. THY1 and NEU2 may be potential therapeutic targets for breast cancer with brain metastases, and THY1, CD1B and DOCK2 may serve as potential prognostic markers for improvement of brain metastases survival.

Tumorigenic and Metastatic Role of CD44-/low/CD24-/low Cells in Luminal Breast Cancer

Cells with high CD44 but low CD24 expression (CD44^high/CD24^-/low) and high aldehyde dehydrogenase activity (ALDH^br) are widely considered to be drivers of metastasis, therapy resistance and tumor recurrence in breast cancer. However, the role of the CD44^high/CD24^-/low and ALDH^br phenotypes in identifying tumorigenic cells in breast cancer remains controversial due to the discrepancy in their distribution and tumorigenic potential in intrinsic breast cancer subtypes. In this study, we analyzed the cells expressing these markers in six different breast cancer cell lines representing major breast cancer subtypes (T47D, MCF-7, BT-474, AU-565, Hs578T and MDA-MB-231). CD44^high/CD24^-/low, ALDH^br and CD44^-/low/CD24^-/low cell populations were isolated by flow cytometry and analyzed for hallmark stem cell characteristics of differentiation, migration, invasiveness and metastasis using in vitro and in vivo techniques. Our results demonstrate that the CD44^-/low/CD24^-/low cell population, which is enriched in luminal cell lines (T47D, MCF-7 and BT-474), possesses metastatic and tumorigenic properties. We also show that, contrary to previous claims, the expression of the ALDH1 isoform ALDH1A1 does not affect the tumorigenic potential of cell lines with high ALDH activity (BT-474 and AU-565). Further transcriptomic and clinical studies are needed to determine the potential of these markers as early diagnostic tools and treatment targets.

Potent Anticancer Effect of the Natural Steroidal Saponin Gracillin Is Produced by Inhibiting Glycolysis and Oxidative Phosphorylation-Mediated Bioenergetics

Metabolic rewiring to utilize aerobic glycolysis is a hallmark of cancer. However, recent findings suggest the role of mitochondria in energy generation in cancer cells and the metabolic switch to oxidative phosphorylation (OXPHOS) in response to the blockade of glycolysis. We previously demonstrated that the antitumor effect of gracillin occurs through the inhibition of mitochondrial complex II-mediated energy production. Here, we investigated the potential of gracillin as an anticancer agent targeting both glycolysis and OXPHOS in breast and lung cancer cells. Along with the reduction in adenosine triphosphate (ATP) production, gracillin markedly suppresses the production of several glycolysis-associated metabolites. A docking analysis and enzyme assay suggested phosphoglycerate kinase 1 (PGK1) is a potential target for the antiglycolytic effect of gracillin. Gracillin reduced the viability and colony formation ability of breast cancer cells by inducing apoptosis. Gracillin displayed efficacious antitumor effects in mice bearing breast cancer cell line or breast cancer patient-derived tumor xenografts with no overt changes in body weight. An analysis of publicly available datasets further suggested that PGK1 expression is associated with metastasis status and poor prognosis in patients with breast cancer. These results suggest that gracillin is a natural anticancer agent that inhibits both glycolysis and mitochondria-mediated bioenergetics.