High-resolution clonal mapping of multi-organ metastasis in triple negative breast cancer

Precise nano-system-based drug delivery and synergistic therapy against androgen receptor-positive triple-negative breast cancer

Targeting androgen receptor (AR) has shown great therapeutic potential in triple-negative breast cancer (TNBC), yet its efficacy remains unsatisfactory. Here, we aimed to identify promising targeted agents that synergize with enzalutamide, a second-generation AR inhibitor, in TNBC. By using a strategy for screening drug combinations based on the Sensitivity Index (SI), we found that MK-8776, a selective checkpoint kinase1 (CHK1) inhibitor, showed favorable synergism with enzalutamide in AR-positive TNBC. The combination of enzalutamide and MK-8776 was found to exert more significant anti-tumor effects in TNBC than the single application of enzalutamide or MK-8776, respectively. Furthermore, a nanoparticle-based on hyaluronic acid (HA)-modified hollow-manganese dioxide (HMnO2), named HMnE&M@H, was established to encapsulate and deliver enzalutamide and MK-8776. This HA-modified nanosystem managed targeted activation via pH/glutathione responsiveness. HMnE&M@H repressed tumor growth more obviously than the simple addition of enzalutamide and MK-8776 without a carrier. Collectively, our study elucidated the synergy of enzalutamide and MK-8776 in TNBC and developed a novel tumor-targeted nano drug delivery system HMnE&M@H, providing a potential therapeutic approach for the treatment of TNBC.

Targeting neddylation and sumoylation in chemoresistant triple negative breast cancer

Triple negative breast cancer (TNBC) accounts for 15-20% of breast cancer cases in the United States. Systemic neoadjuvant chemotherapy (NACT), with or without immunotherapy, is the current standard of care for patients with early-stage TNBC. However, up to 70% of TNBC patients have significant residual disease once NACT is completed, which is associated with a high risk of developing recurrence within two to three years of surgical resection. To identify targetable vulnerabilities in chemoresistant TNBC, we generated longitudinal patient-derived xenograft (PDX) models from TNBC tumors before and after patients received NACT. We then compiled transcriptomes and drug response profiles for all models. Transcriptomic analysis identified the enrichment of aberrant protein homeostasis pathways in models from post-NACT tumors relative to pre-NACT tumors. This observation correlated with increased sensitivity in vitro to inhibitors targeting the proteasome, heat shock proteins, and neddylation pathways. Pevonedistat, a drug annotated as a NEDD8-activating enzyme (NAE) inhibitor, was prioritized for validation in vivo and demonstrated efficacy as a single agent in multiple PDX models of TNBC. Pharmacotranscriptomic analysis identified a pathway-level correlation between pevonedistat activity and post-translational modification (PTM) machinery, particularly involving neddylation and sumoylation targets. Elevated levels of both NEDD8 and SUMO1 were observed in models exhibiting a favorable response to pevonedistat compared to those with a less favorable response in vivo. Moreover, a correlation emerged between the expression of neddylation-regulated pathways and tumor response to pevonedistat, indicating that targeting these PTM pathways may prove effective in combating chemoresistant TNBC.

Spatial multi-omics at subcellular resolution via high-throughput in situ pairwise sequencing

Technology for spatial multi-omics aids the discovery of new insights into cellular functions and disease mechanisms. Here we report the development and applicability of multi-omics in situ pairwise sequencing (MiP-seq), a method for the simultaneous detection of DNAs, RNAs, proteins and biomolecules at subcellular resolution. Compared with other in situ sequencing methods, MiP-seq enhances decoding capacity and reduces sequencing and imaging costs while maintaining the efficacy of detection of gene mutations, allele-specific expression and RNA modifications. MiP-seq can be integrated with in vivo calcium imaging and Raman imaging, which enabled us to generate a spatial multi-omics atlas of mouse brain tissues and to correlate gene expression with neuronal activity and cellular biochemical fingerprints. We also report a sequential dilution strategy for resolving optically crowded signals during in situ sequencing. High-throughput in situ pairwise sequencing may facilitate the multidimensional analysis of molecular and functional maps of tissues.

Development and Validation of a Prognostic Nomogram for Breast Cancer Patients With Multi-Organ Metastases: An Analysis of the Surveillance, Epidemiology, and End Results Program Database

BACKGROUND

Multi-organ metastases represent a substantial life-threatening risk for breast cancer (BC) patients. Nonetheless, the current dearth of assessment tools for patients with multi-organ metastatic BC adversely impacts their evaluation.

METHODS

We conducted a retrospective analysis of BC patients with multi-organ metastases using data from the SEER database from 2010 to 2019. The patients were randomly allocated into a training cohort and a validation cohort in a 7:3 ratio. Univariate COX regression analysis, the LASSO, and multivariate Cox regression analyses were performed to identify independent prognostic factors in the training set. Based on these factors, a nomogram was constructed to estimate overall survival (OS) probability for BC patients with multi-organ metastases. The performance of the nomogram was evaluated using C-indexes, ROC curves, calibration curves, decision curve analysis (DCA) curves, and the risk classification system for validation.

RESULTS

A total of 3626 BC patients with multi-organ metastases were included in the study, with 2538 patients in the training cohort and 1088 patients in the validation cohort. Age, grade, metastasis location, surgery, chemotherapy, and subtype were identified as significant independent prognostic factors for OS in BC patients with multi-organ metastases. A nomogram for predicting 1-year, 3-year, and 5-year OS was constructed. The evaluation metrics, including C-indexes, ROC curves, calibration curves, and DCA curves, demonstrated the excellent predictive performance of the nomogram. Additionally, the risk grouping system effectively stratified BC patients with multi-organ metastases into distinct prognostic categories.

CONCLUSION

The developed nomogram showed high accuracy in predicting the survival probability of BC patients with multi-organ metastases, providing valuable information for patient counseling and treatment decision making.

Chemoresistance to additive PARP/PI3K dual inhibition in triple-negative breast cancer cell lines is associated with adaptive stem cell-like prevalence

Cancer stem-like cells (CSCs) are posited to exhibit specialized oncogenic capacity to drive malignancies. CSCs are distinguished by enhanced hallmarks of cancer, including apoptosis avoidance, phenotypic plasticity and aberrant growth pathway signaling. Standard-of-care chemotherapies targeted to rapidly cycling cells routinely fail to eliminate this resistant subpopulation, leading to disease recurrence and metastasis. Triple-negative breast cancer (TNBC), a highly aggressive subtype of breast cancer, is enriched for tumor-propagating CD44+/CD24-/low CSCs, which are poorly ablated by chemotherapeutics and are associated with poor prognosis. CD44 governs sustained PI3K signaling in breast cancer, which is essential for CSC maintenance. PI3K inhibition can elicit DNA damage and down-regulate BRCA1 expression, which in turn enhance the synthetic lethality of PARP inhibitors. Here, we examined a dual chemotherapeutic approach targeting these pathways by combining a pan-PI3K inhibitor (Buparlisib) and a PARP1 inhibitor (Olaparib) on a panel of TNBC cell lines with distinct mutational profiles and proportions of CSCs. We observed differential sensitivity to this dual inhibition strategy and varying cellular stress and resistance responses across eight TNBC lines. The dual chemotherapeutic effect is associated with a reduction in S-phase cells, an increased in apoptotic cells and elevated expression of cleaved PARP, indicating a provoked replicative stress response. We observed that PARP/PI3K inhibition efficacy was potentiated by repeated administration in some TNBC lines and identified critical treatment schedules, which further potentiated the dual chemotherapeutic effect. Dual inhibition induced small but significant increases in CSC relative abundance as marked by CD44+/CD24-/low or ALDH1+ cells and increased stress and survival signaling in multiple TNBC cell lines, suggesting this sub-population contributes to TNBC chemoresistance. These results suggest the additive effects of PARP and PI3K inhibition against CSC phenotypes may be enhanced by temporally-staged administration in TNBC cells.

Context-dependent roles for ubiquitous mitochondrial creatine kinase CKMT1 in breast cancer progression

Metabolic reprogramming is a hallmark of cancer, enabling cancer cells to rapidly proliferate, invade, and metastasize. We show that creatine levels in metastatic breast cancer cell lines and secondary metastatic tumors are driven by the ubiquitous mitochondrial creatine kinase (CKMT1). We discover that, while CKMT1 is highly expressed in primary tumors and promotes cell viability, it is downregulated in metastasis. We further show that CKMT1 downregulation, as seen in breast cancer metastasis, drives up mitochondrial reactive oxygen species (ROS) levels. CKMT1 downregulation contributes to the migratory and invasive potential of cells by ROS-induced upregulation of adhesion and degradative factors, which can be reversed by antioxidant treatment. Our study thus reconciles conflicting evidence about the roles of metabolites in the creatine metabolic pathway in breast cancer progression and reveals that tight, context-dependent regulation of CKMT1 expression facilitates cell viability, cell migration, and cell invasion, which are hallmarks of metastatic spread.

Invasive papillary carcinoma of the breast

Invasive papillary carcinoma is a rare form of breast cancer that is more likely to occur in postmenopausal women. Previous studies have been limited to case reports and small retrospective studies, leading to low awareness of this type of tumor and difficult clinical management. According to the available literature, invasive papillary carcinoma exhibits unique pathological features and biological behaviors. Invasive papillary carcinoma is mostly luminal type, with a low rate of lymph node metastasis, which underlies its favorable prognosis. The effectiveness of adjuvant therapy in reducing tumor burden and improving prognosis in patients with invasive papillary carcinoma remains uncertain. Due to the rarity of the lesion, conducting prospective clinical trials is impractical. The use of biological models, such as organoids, can help alleviate the impact of the scarcity of this condition on research. In addition, invasive papillary carcinoma is affected by specific genomic events, and more extensive studies of gene expression profiling may provide molecular-level insights to make optimal therapeutic decisions.

Natural History of Germline BRCA1 Mutated and BRCA Wild-type Triple-negative Breast Cancer

We report a deep next-generation sequencing analysis of 13 sequentially obtained tumor samples, eight sequentially obtained circulating tumor DNA (ctDNA) samples and three germline DNA samples over the life history of 3 patients with triple-negative breast cancer (TNBC), 2 of whom had germline pathogenic BRCA1 mutation, to unravel tumor evolution. Tumor tissue from all timepoints and germline DNA was subjected to whole-exome sequencing (WES), custom amplicon deep sequencing (30,000X) of a WES-derived somatic mutation panel, and SNP arrays for copy-number variation (CNV), while whole transcriptome sequencing (RNA-seq) was performed only on somatic tumor.There was enrichment of homologous recombination deficiency signature in all tumors and widespread CNV, which remained largely stable over time. Somatic tumor mutation numbers varied between patients and within each patient (range: 70-216, one outlier). There was minimal mutational overlap between patients with TP53 being the sole commonly mutated gene, but there was substantial overlap in sequential samples in each patient. Each patient's tumor contained a founding ("stem") clone at diagnosis, which persisted over time, from which all other clones ("subclone") were derived ("branching evolution"), which contained mutations in well-characterized cancer-related genes like PDGFRB, ARID2, TP53 (Patient_02), TP53, BRAF, BRIP1, CSF3R (Patient_04), and TP53, APC, EZH2 (Patient_07). Including stem and subclones, tumors from all patients were polyclonal at diagnosis and during disease progression. ctDNA recapitulated most tissue-derived stem clonal and subclonal mutations while detecting some additional subclonal mutations. RNA-seq revealed a stable basal-like pattern, with most highly expressed variants belonging to stem clone.

SIGNIFICANCE

In germline BRCA1 mutated and BRCA wild-type patients, TNBC shows a branching evolutionary pattern of mutations with a single founding clone, are polyclonal throughout their disease course, and have widespread copy-number aberrations. This evolutionary pattern may be associated with treatment resistance or sensitivity and could be therapeutically exploited.

Extracting, filtering and simulating cellular barcodes using CellBarcode tools

Identifying true DNA cellular barcodes among polymerase chain reaction and sequencing errors is challenging. Current tools are restricted in the diversity of barcode types supported or the analysis strategies implemented. As such, there is a need for more versatile and efficient tools for barcode extraction, as well as for tools to investigate which factors impact barcode detection and which filtering strategies to best apply. Here we introduce the package CellBarcode and its barcode simulation kit, CellBarcodeSim, that allows efficient and versatile barcode extraction and filtering for a range of barcode types from bulk or single-cell sequencing data using a variety of filtering strategies. Using the barcode simulation kit and biological data, we explore the technical and biological factors influencing barcode identification and provide a decision tree on how to optimize barcode identification for different barcode settings. We believe that CellBarcode and CellBarcodeSim have the capability to enhance the reproducibility and interpretation of barcode results across studies.

Tumor-derived cell-free DNA and circulating tumor cells: partners or rivals in metastasis formation?

The origin of metastases is a topic that has sparked controversy. Despite recent advancements, metastatic disease continues to pose challenges. The first admitted model of how metastases develop revolves around cells breaking away from the primary tumor, known as circulating tumor cells (CTCs). These cells survive while circulating through the bloodstream and subsequently establish themselves in secondary organs, a process often referred to as the "metastatic cascade". This intricate and dynamic process involves various steps, but all the mechanisms behind metastatic dissemination are not yet comprehensively elucidated. The "seed and soil" theory has shed light on the phenomenon of metastatic organotropism and the existence of pre-metastatic niches. It is now established that these niches can be primed by factors secreted by the primary tumor before the arrival of CTCs. In particular, exosomes have been identified as important contributors to this priming. Another concept then emerged, i.e. the "genometastasis" theory, which challenged all other postulates. It emphasizes the intriguing but promising role of cell-free DNA (cfDNA) in metastasis formation through oncogenic formation of recipient cells. However, it cannot be ruled out that all these theories are intertwined. This review outlines the primary theories regarding the metastases formation that involve CTCs, and depicts cfDNA, a potential second player in the metastasis formation. We discuss the potential interrelationships between CTCs and cfDNA, and propose both in vitro and in vivo experimental strategies to explore all plausible theories.

Patient-Derived Xenografts: A Valuable Preclinical Model for Drug Development and Biomarker Discovery

Patient-derived xenografts (PDXs), established by implanting patient tumor cells into immunodeficient mice, offer a platform for faithfully replicating human tumors. They closely mimic the histopathology, genomics, and drug sensitivity of patient tumors. This chapter highlights the versatile applications of PDXs, including studying tumor biology, metastasis, and chemoresistance, as well as their use in biomarker identification, drug screening, and personalized medicine. It also addresses challenges in using PDXs in cancer research, including variations in metastatic potential, lengthy establishment timelines, stromal changes, and limitations in immunocompromised models. Despite these challenges, PDXs remain invaluable tools guiding patient treatment and advancing preclinical drug development.

New insights into the correlations between circulating tumor cells and target organ metastasis

Organ-specific metastasis is the primary cause of cancer patient death. The distant metastasis of tumor cells to specific organs depends on both the intrinsic characteristics of the tumor cells and extrinsic factors in their microenvironment. During an intermediate stage of metastasis, circulating tumor cells (CTCs) are released into the bloodstream from primary and metastatic tumors. CTCs harboring aggressive or metastatic features can extravasate to remote sites for continuous colonizing growth, leading to further lesions. In the past decade, numerous studies demonstrated that CTCs exhibited huge clinical value including predicting distant metastasis, assessing prognosis and monitoring treatment response et al. Furthermore, increasingly numerous experiments are dedicated to identifying the key molecules on or inside CTCs and exploring how they mediate CTC-related organ-specific metastasis. Based on the above molecules, more and more inhibitors are being developed to target CTCs and being utilized to completely clean CTCs, which should provide promising prospects to administer advanced tumor. Recently, the application of various nanomaterials and microfluidic technologies in CTCs enrichment technology has assisted to improve our deep insights into the phenotypic characteristics and biological functions of CTCs as a potential therapy target, which may pave the way for us to make practical clinical strategies. In the present review, we mainly focus on the role of CTCs being involved in targeted organ metastasis, especially the latest molecular mechanism research and clinical intervention strategies related to CTCs.

Patient Derived Xenografts (PDX) Models as an Avatar to Assess Personalized Therapy Options in Uveal Melanoma: A Feasibility Study

Uveal melanoma is the most common primary intraocular malignancy in adults. Up to 50% of UM patients develop metastatic disease, usually in the liver. When metastatic, the prognosis is poor, and few treatment options exist. Here, we investigated the feasibility of establishing patient-derived xenografts (PDXs) from a patient's tumor in order to screen for therapies that the patient could benefit from. Samples obtained from 29 primary tumors and liver metastases of uveal melanoma were grafted into SCID mice. PDX models were successfully established for 35% of primary patient tumors and 67% of liver metastases. The tumor take rate was proportional to the risk of metastases. PDXs showed the same morphology, the same GNAQ/11, BAP1, and SF3B1 mutations, and the same chromosome 3 and 8q status as the corresponding patient samples. Six PDX models were challenged with two compounds for 4 weeks. We show that, for 31% of patients with high or intermediate risk of metastasis, the timing to obtain efficacy results on PDX models derived from their primary tumors was compatible with the selection of the therapy to treat the patient after relapse. PDXs could thus be a valid tool ("avatar") to select the best personalized therapy for one third of patients that are most at risk of relapse.

The role of E3 ubiquitin ligases in bone homeostasis and related diseases

The ubiquitin-proteasome system (UPS) dedicates to degrade intracellular proteins to modulate demic homeostasis and functions of organisms. These enzymatic cascades mark and modifies target proteins diversly through covalently binding ubiquitin molecules. In the UPS, E3 ubiquitin ligases are the crucial constituents by the advantage of recognizing and presenting proteins to proteasomes for proteolysis. As the major regulators of protein homeostasis, E3 ligases are indispensable to proper cell manners in diverse systems, and they are well described in physiological bone growth and bone metabolism. Pathologically, classic bone-related diseases such as metabolic bone diseases, arthritis, bone neoplasms and bone metastasis of the tumor, etc., were also depicted in a UPS-dependent manner. Therefore, skeletal system is versatilely regulated by UPS and it is worthy to summarize the underlying mechanism. Furthermore, based on the current status of treatment, normal or pathological osteogenesis and tumorigenesis elaborated in this review highlight the clinical significance of UPS research. As a strategy possibly remedies the limitations of UPS treatment, emerging PROTAC was described comprehensively to illustrate its potential in clinical application. Altogether, the purpose of this review aims to provide more evidence for exploiting novel therapeutic strategies based on UPS for bone associated diseases.

Xenomake: a pipeline for processing and sorting xenograft reads from spatial transcriptomic experiments

Xenograft models are attractive models that mimic human tumor biology and permit one to perturb the tumor microenvironment and study its drug response. Spatially resolved transcriptomics (SRT) provide a powerful way to study the organization of xenograft models, but currently there is a lack of specialized pipeline for processing xenograft reads originated from SRT experiments. Xenomake is a standalone pipeline for the automated handling of spatial xenograft reads. Xenomake handles read processing, alignment, xenograft read sorting, quantification, and connects well with downstream spatial analysis packages. We additionally show that Xenomake can correctly assign organism specific reads, reduce sparsity of data by increasing gene counts, while maintaining biological relevance for studies.

Ethnic disparities in the immune microenvironment of triple negative breast cancer and its role in therapeutic outcomes

In 2020, newly diagnosed breast cancer (BC) cases surpassed that of lung cancer among women, making it the most common female cancer globally. In spite of recent increases in incidence rates, mortality due to BC has declined since 1989. These declines have been attributed to advancements in treatment modalities as well as increased mammography surveillance. Despite these advances, African American (AA) women are 40% more likely to die from BC than Caucasian women. Multifactorial etiology has been implicated in the disparity of BC mortality rates among AA women. As an example, AA women have a disproportionate incidence of triple negative breast cancer (TNBC), which has a poor prognosis and marginal treatment options. Increasingly, the tumor microenvironment (TME) has gained relevance as it relates to primary tumor progression, metastasis and treatment possibilities. The treatment outcomes or pathological complete response (pCR) in TNBC among AA women are affected by differences in TME. The TME of AA women exhibit several variances in acellular and cellular components associated with pro-tumorigenic effects. For example, increased levels of the adipocyte-related hormone, resistin, the pro-inflammatory cytokine, IL-6, and the CC chemokine, CCL2, within the TME of AA women gives rise to an increased density of M2 macrophages, also known as tumor-associated macrophages. Elevated levels of vascular endothelial growth factor in the TME of AA women increase the vascular density or vascularity, which facilitate aggressive tumor growth and metastasis. Furthermore, a pro-tumorigenic TME is supported by increased levels of the CXC chemokine, CXCL12 that results in the recruitment of regulatory T lymphocytes (Tregs ). Due to these and other differences in the TME of AA women, precision oncology can target specific aspects of the TME that may contribute to a poorer prognosis. In addition to the discrepancies in the TME, AA women face socio-economic barriers that limit their ability to access state-of-the-art, novel therapies against metastatic TNBC. In this review, we will provide a brief overview of the tumor immune microenvironment, immune-based treatment options for TNBC and their potential to decrease health disparities due to ethnicity.

Experimental and spontaneous metastasis assays can result in divergence in clonal architecture

Intratumoural heterogeneity is associated with poor outcomes in breast cancer. To understand how malignant clones survive and grow in metastatic niches, in vivo models using cell lines and patient-derived xenografts (PDX) have become the gold standard. Injections of cancer cells in orthotopic sites (spontaneous metastasis assays) or into the vasculature (experimental metastasis assays) have been used interchangeably to study the metastatic cascade from early events or post-intravasation, respectively. However, less is known about how these different routes of injection impact heterogeneity. Herein we directly compared the clonality of spontaneous and experimental metastatic assays using the human cell line MDA-MB-231 and a PDX model. Genetic barcoding was used to study the fitness of the subclones in primary and metastatic sites. Using spontaneous assays, we found that intraductal injections resulted in less diverse tumours compared to other routes of injections. Using experimental metastasis assays via tail vein injection of barcoded MDA-MB-231 cells, we also observed an asymmetry in metastatic heterogeneity between lung and liver that was not observed using spontaneous metastasis assays. These results demonstrate that these assays can result in divergent clonal outputs in terms of metastatic heterogeneity and provide a better understanding of the biases inherent to each technique.

Implications of tumour heterogeneity on cancer evolution and therapy resistance: lessons from breast cancer

Tumour heterogeneity is pervasive amongst many cancers and leads to disease progression, and therapy resistance. In this review, using breast cancer as an exemplar, we focus on the recent advances in understanding the interplay between tumour cells and their microenvironment using single cell sequencing and digital spatial profiling technologies. Further, we discuss the utility of lineage tracing methodologies in pre-clinical models of breast cancer, and how these are being used to unravel new therapeutic vulnerabilities and reveal biomarkers of breast cancer progression. © 2023 The Pathological Society of Great Britain and Ireland.

A transcriptome-informed QSP model of metastatic triple-negative breast cancer identifies predictive biomarkers for PD-1 inhibition

Triple-negative breast cancer (TNBC), a highly metastatic breast cancer subtype, has limited treatment options. While a small number of patients attain clinical benefit with single-agent checkpoint inhibitors, identifying these patients before the therapy remains challenging. Here, we developed a transcriptome-informed quantitative systems pharmacology model of metastatic TNBC by integrating heterogenous metastatic tumors. In silico clinical trial with an anti-PD-1 drug, pembrolizumab, predicted that several features, such as the density of antigen-presenting cells, the fraction of cytotoxic T cells in lymph nodes, and the richness of cancer clones in tumors, could serve individually as biomarkers but had a higher predictive power as combinations of two biomarkers. We showed that PD-1 inhibition neither consistently enhanced all antitumorigenic factors nor suppressed all protumorigenic factors but ultimately reduced the tumor carrying capacity. Collectively, our predictions suggest several candidate biomarkers that might effectively predict the response to pembrolizumab monotherapy and potential therapeutic targets to develop treatment strategies for metastatic TNBC.

Glycyrrhetinic acid suppresses breast cancer metastasis by inhibiting M2-like macrophage polarization via activating JNK1/2 signaling

BACKGROUND

Breast cancer metastasis is leading cause of cancer death among women worldwide. Tumor-associated macrophages (TAMs) have been considered as potential targets for treating breast cancer metastasis because they promote tumor growth and development. Glycyrrhetinic acid (GA) is one of the most important phytochemicals of licorice which has shown promising anti-cancer efficacies in pre-clinical trials. However, the regulatory effect of GA on the polarization of TAMs remains elusive.

PURPOSE

To investigate the role of GA in regulating the polarization of M2 macrophages and inhibiting breast cancer metastasis, and to further explore its underlying mechanisms of action.

STUDY DESIGN

IL-4 / IL-13-treated RAW 264.7 and THP-1 cells were used as the M2-polarized macrophages in vitro. A 4T1 mouse breast cancer model and the tail vein breast cancer metastasis model were applied to study the effect of GA on breast cancer growth and metastasis in vivo.

RESULTS

In vitro studies showed that GA significantly inhibited IL-4 / IL 13-induced M2-like polarization in RAW 264.7 and THP-1 macrophages without affecting M1-like polarization. GA strongly decreased the expression of M2 macrophage markers CD206 and Arg-1, and reduced the levels of the pro-angiogenic molecules VEGF, MMP9, MMP2 and IL-10 in M2 macrophages. GA also increased the phosphorylation of JNK1/2 in M2 macrophages. Moreover, GA significantly suppressed M2 macrophage-induced cell proliferation and migration in 4T1 cancer cells and HUVECs. Interestingly, the inhibitory effects of GA on M2 macrophages were abolished by a JNK inhibitor. Animal studies showed that GA significantly suppressed tumor growth, angiogenesis, and lung metastasis in BALB/c mice bearing breast tumor. In tumor tissues, GA reduced the number of M2 macrophages but elevated the proportion of M1 macrophages, accompanied by activation of JNK signaling. Similar results were found in the tail vein breast cancer metastasis model.

CONCLUSION

This study demonstrated for the first time that GA could effectively suppress breast cancer growth and metastasis by inhibiting macrophage M2 polarization via activating JNK1/2 signaling. These findings indicate that GA could be served as the lead compound for the future development of anti-breast cancer drug.

Mitochondrial structure and function adaptation in residual triple negative breast cancer cells surviving chemotherapy treatment

Neoadjuvant chemotherapy (NACT) used for triple negative breast cancer (TNBC) eradicates tumors in ~45% of patients. Unfortunately, TNBC patients with substantial residual cancer burden have poor metastasis free and overall survival rates. We previously demonstrated mitochondrial oxidative phosphorylation (OXPHOS) was elevated and was a unique therapeutic dependency of residual TNBC cells surviving NACT. We sought to investigate the mechanism underlying this enhanced reliance on mitochondrial metabolism. Mitochondria are morphologically plastic organelles that cycle between fission and fusion to maintain mitochondrial integrity and metabolic homeostasis. The functional impact of mitochondrial structure on metabolic output is highly context dependent. Several chemotherapy agents are conventionally used for neoadjuvant treatment of TNBC patients. Upon comparing mitochondrial effects of conventional chemotherapies, we found that DNA-damaging agents increased mitochondrial elongation, mitochondrial content, flux of glucose through the TCA cycle, and OXPHOS, whereas taxanes instead decreased mitochondrial elongation and OXPHOS. The mitochondrial effects of DNA-damaging chemotherapies were dependent on the mitochondrial inner membrane fusion protein optic atrophy 1 (OPA1). Further, we observed heightened OXPHOS, OPA1 protein levels, and mitochondrial elongation in an orthotopic patient-derived xenograft (PDX) model of residual TNBC. Pharmacologic or genetic disruption of mitochondrial fusion and fission resulted in decreased or increased OXPHOS, respectively, revealing longer mitochondria favor oxphos in TNBC cells. Using TNBC cell lines and an in vivo PDX model of residual TNBC, we found that sequential treatment with DNA-damaging chemotherapy, thus inducing mitochondrial fusion and OXPHOS, followed by MYLS22, a specific inhibitor of OPA1, was able to suppress mitochondrial fusion and OXPHOS and significantly inhibit regrowth of residual tumor cells. Our data suggest that TNBC mitochondria can optimize OXPHOS through OPA1-mediated mitochondrial fusion. These findings may provide an opportunity to overcome mitochondrial adaptations of chemoresistant TNBC.

Orchestration of Collective Migration and Metastasis by Tumor Cell Clusters

Metastatic dissemination has lethal consequences for cancer patients. Accruing evidence supports the hypothesis that tumor cells can migrate and metastasize as clusters of cells while maintaining contacts with one another. Collective metastasis enables tumor cells to colonize secondary sites more efficiently, resist cell death, and evade the immune system. On the other hand, tumor cell clusters face unique challenges for dissemination particularly during systemic dissemination. Here, we review recent progress toward understanding how tumor cell clusters overcome these disadvantages as well as mechanisms they utilize to gain advantages throughout the metastatic process. We consider useful models for studying collective metastasis and reflect on how the study of collective metastasis suggests new opportunities for eradicating and preventing metastatic disease.

Identifying biomarkers of differential chemotherapy response in TNBC patient-derived xenografts with a CTD/WGCNA approach

Although systemic chemotherapy remains the standard of care for TNBC, even combination chemotherapy is often ineffective. The identification of biomarkers for differential chemotherapy response would allow for the selection of responsive patients, thus maximizing efficacy and minimizing toxicities. Here, we leverage TNBC PDXs to identify biomarkers of response. To demonstrate their ability to function as a preclinical cohort, PDXs were characterized using DNA sequencing, transcriptomics, and proteomics to show consistency with clinical samples. We then developed a network-based approach (CTD/WGCNA) to identify biomarkers of response to carboplatin (MSI1, TMSB15A, ARHGDIB, GGT1, SV2A, SEC14L2, SERPINI1, ADAMTS20, DGKQ) and docetaxel (c, MAGED4, CERS1, ST8SIA2, KIF24, PARPBP). CTD/WGCNA multigene biomarkers are predictive in PDX datasets (RNAseq and Affymetrix) for both taxane- (docetaxel or paclitaxel) and platinum-based (carboplatin or cisplatin) response, thereby demonstrating cross-expression platform and cross-drug class robustness. These biomarkers were also predictive in clinical datasets, thus demonstrating translational potential.

Predictors of success in establishing orthotopic patient-derived xenograft models of triple negative breast cancer

Patient-derived xenograft (PDX) models of breast cancer are an effective discovery platform and tool for preclinical pharmacologic testing and biomarker identification. We established orthotopic PDX models of triple negative breast cancer (TNBC) from the primary breast tumors of patients prior to and following neoadjuvant chemotherapy (NACT) while they were enrolled in the ARTEMIS trial (NCT02276443). Serial biopsies were obtained from patients prior to treatment (pre-NACT), from poorly responsive disease after four cycles of Adriamycin and cyclophosphamide (AC, mid-NACT), and in cases of AC-resistance, after a 3-month course of different experimental therapies and/or additional chemotherapy (post-NACT). Our study cohort includes a total of 269 fine needle aspirates (FNAs) from 217 women, generating a total of 62 PDX models (overall success-rate = 23%). Success of PDX engraftment was generally higher from those cancers that proved to be treatment-resistant, whether poorly responsive to AC as determined by ultrasound measurements mid-NACT (p = 0.063), RCB II/III status after NACT (p = 0.046), or metastatic relapse within 2 years of surgery (p = 0.008). TNBC molecular subtype determined from gene expression microarrays of pre-NACT tumors revealed no significant association with PDX engraftment rate (p = 0.877). Finally, we developed a statistical model predictive of PDX engraftment using percent Ki67 positive cells in the patient's diagnostic biopsy, positive lymph node status at diagnosis, and low volumetric reduction of the patient's tumor following AC treatment. This novel bank of 62 PDX models of TNBC provides a valuable resource for biomarker discovery and preclinical therapeutic trials aimed at improving neoadjuvant response rates for patients with TNBC.

The complex nature of heterogeneity and its roles in breast cancer biology and therapeutic responsiveness

Heterogeneity is a complex feature of cells and tissues with many interacting components. Depending on the nature of the research context, interacting features of cellular, drug response, genetic, molecular, spatial, temporal, and vascular heterogeneity may be present. We describe the various forms of heterogeneity with examples of their interactions and how they play a role in affecting cellular phenotype and drug responses in breast cancer. While cellular heterogeneity may be the most widely described and invoked, many forms of heterogeneity are evident within the tumor microenvironment and affect responses to the endocrine and cytotoxic drugs widely used in standard clinical care. Drug response heterogeneity is a critical determinant of clinical response and curative potential and also is multifaceted when encountered. The interactive nature of some forms of heterogeneity is readily apparent. For example, the process of metastasis has the properties of both temporal and spatial heterogeneity within the host, whereas each individual metastatic deposit may exhibit cellular, genetic, molecular, and vascular heterogeneity. This review describes the many forms of heterogeneity, their integrated activities, and offers some insights into how heterogeneity may be understood and studied in the future.

Non-Invasive and Real-Time Monitoring of the Breast Cancer Metastasis Degree via Metabolomics

Breast cancer (BC) is a serious threat to women's health and metastasis is the major cause of BC-associated mortality. Various techniques are currently used to preoperatively describe the metastatic status of tumors, based on which a comprehensive treatment protocol was determined. However, accurately staging a tumor before surgery remains a challenge, which may lead to the miss of optimal treatment options. More severely, the failure to detect and remove occult micrometastases often causes tumor recurrences. There is an urgent need to develop a more precise and non-invasive strategy for the detection of the tumor metastasis in lymph nodes and distant organs. Based on the facts that tumor metastasis is closely related to the primary tumor microenvironment (TME) evolutions and that metabolomics profiling of the circulatory system can precisely reflect subtle changes within TME, we suppose whether metabolomic technology can be used to achieve non-invasive and real-time monitoring of BC metastatic status. In this study, the metastasis status of BC mouse models with different tumor-bearing times was firstly depicted to mimic clinical anatomic TNM staging system. Metabolomic profiling together with metastasis-related changes in TME among tumor-bearing mice with different metastatic status was conducted. A range of differential metabolites reflecting tumor metastatic states were screened and in vivo experiments proved that two main metastasis-driving factors in TME, TGF-β and hypoxia, were closely related to the regular changes of these metabolites. The differential metabolites level changes were also preliminarily confirmed in a limited number of clinical BC samples. Metabolite lysoPC (16:0) was found to be useful for clinical N stage diagnosis and the possible cause of its changes was analyzed by bioinformatics techniques.

DNA barcoding reveals ongoing immunoediting of clonal cancer populations during metastatic progression and immunotherapy response

Cancers evade the immune system through the process of cancer immunoediting. While immune checkpoint inhibitors are effective for reactivating tumour immunity in some cancer types, many other solid cancers, including breast cancer, remain largely non-responsive. Understanding how non-responsive cancers evade immunity and whether this occurs at the clonal level will improve immunotherapeutic design. Here we use DNA barcoding to track murine mammary cancer cell clones during immunoediting and determine clonal transcriptional profiles that allow immune evasion following anti-PD1 plus anti-CTLA4 immunotherapy. Clonal diversity is significantly restricted by immunotherapy treatment in both primary tumours and metastases, demonstrating selection for pre-existing breast cancer cell populations and ongoing immunoediting during metastasis and treatment. Immunotherapy resistant clones express a common gene signature associated with poor survival of basal-like breast cancer patient cohorts. At least one of these genes has an existing small molecule that can potentially be used to improve immunotherapy response.

Mastering the use of cellular barcoding to explore cancer heterogeneity

Tumours are often composed of a multitude of malignant clones that are genomically unique, and only a few of them may have the ability to escape cancer therapy and grow as symptomatic lesions. As a result, tumours with a large degree of genomic diversity have a higher chance of leading to patient death. However, clonal fate can be driven by non-genomic features. In this context, new technologies are emerging not only to track the spatiotemporal fate of individual cells and their progeny but also to study their molecular features using various omics analysis. In particular, the recent development of cellular barcoding facilitates the labelling of tens to millions of cancer clones and enables the identification of the complex mechanisms associated with clonal fate in different microenvironments and in response to therapy. In this Review, we highlight the recent discoveries made using lentiviral-based cellular barcoding techniques, namely genetic and optical barcoding. We also emphasize the strengths and limitations of each of these technologies and discuss some of the key concepts that must be taken into consideration when one is designing barcoding experiments. Finally, we suggest new directions to further improve the use of these technologies in cancer research.

Molecular Characterization and Prospective Evaluation of Pathologic Response and Outcomes with Neoadjuvant Therapy in Metaplastic Triple-Negative Breast Cancer

PURPOSE

Metaplastic breast cancer (MpBC) is a rare subtype of breast cancer that is commonly triple-negative and poorly responsive to neoadjuvant therapy in retrospective studies.

EXPERIMENTAL DESIGN

To better define clinical outcomes and correlates of response, we analyzed the rate of pathologic complete response (pCR) to neoadjuvant therapy, survival outcomes, and genomic and transcriptomic profiles of the pretreatment tumors in a prospective clinical trial (NCT02276443). A total of 211 patients with triple-negative breast cancer (TNBC), including 39 with MpBC, received doxorubicin-cyclophosphamide-based neoadjuvant therapy.

RESULTS

Although not meeting the threshold for statistical significance, patients with MpBCs were less likely to experience a pCR (23% vs. 40%; P = 0.07), had shorter event-free survival (29.4 vs. 32.2 months, P = 0.15), metastasis-free survival (30.3 vs. 32.4 months, P = 0.22); and overall survival (32.6 vs. 34.3 months, P = 0.21). This heterogeneity is mirrored in the molecular profiling. Mutations in PI3KCA (23% vs. 9%, P = 0.07) and its pathway (41% vs. 18%, P = 0.02) were frequently observed and enriched in MpBCs. The gene expression profiles of each histologically defined subtype were distinguishable and characterized by distinctive gene signatures. Among nonmetaplastic (non-Mp) TNBCs, 10% possessed a metaplastic-like gene expression signature and had pCR rates and survival outcomes similar to MpBC.

CONCLUSIONS

Further investigations will determine if metaplastic-like tumors should be treated more similarly to MpBC in the clinic. The 23% pCR rate in this study suggests that patients with MpBC should be considered for NAT. To improve this rate, a pathway analysis predicted enrichment of histone deacetylase (HDAC) and RTK/MAPK pathways in MpBC, which may serve as new targetable vulnerabilities.

DNA barcoded competitive clone-initiating cell analysis reveals novel features of metastatic growth in a cancer xenograft model

A problematic feature of many human cancers is a lack of understanding of mechanisms controlling organ-specific patterns of metastasis, despite recent progress in identifying many mutations and transcriptional programs shown to confer this potential. To address this gap, we developed a methodology that enables different aspects of the metastatic process to be comprehensively characterized at a clonal resolution. Our approach exploits the application of a computational pipeline to analyze and visualize clonal data obtained from transplant experiments in which a cellular DNA barcoding strategy is used to distinguish the separate clonal contributions of two or more competing cell populations. To illustrate the power of this methodology, we demonstrate its ability to discriminate the metastatic behavior in immunodeficient mice of a well-established human metastatic cancer cell line and its co-transplanted LRRC15 knockdown derivative. We also show how the use of machine learning to quantify clone-initiating cell (CIC) numbers and their subsequent metastatic progeny generated in different sites can reveal previously unknown relationships between different cellular genotypes and their initial sites of implantation with their subsequent respective dissemination patterns. These findings underscore the potential of such combined genomic and computational methodologies to identify new clonally-relevant drivers of site-specific patterns of metastasis.

Model-based inference of metastatic seeding rates in de novo metastatic breast cancer reveals the impact of secondary seeding and molecular subtype

We present a stochastic network model of metastasis spread for de novo metastatic breast cancer, composed of tumor to metastasis (primary seeding) and metastasis to metastasis spread (secondary seeding), parameterized using the SEER (Surveillance, Epidemiology, and End Results) database. The model provides a quantification of tumor cell dissemination rates between the tumor and metastasis sites. These rates were used to estimate the probability of developing a metastasis for untreated patients. The model was validated using tenfold cross-validation. We also investigated the effect of HER2 (Human Epidermal Growth Factor Receptor 2) status, estrogen receptor (ER) status and progesterone receptor (PR) status on the probability of metastatic spread. We found that dissemination rate through secondary seeding is up to 300 times higher than through primary seeding. Hormone receptor positivity promotes seeding to the bone and reduces seeding to the lungs and primary seeding to the liver, while HER2 expression increases dissemination to the bone, lungs and primary seeding to the liver. Secondary seeding from the lungs to the liver seems to be hormone receptor-independent, while that from the lungs to the brain appears HER2-independent.

Identification of two cancer stem cell-like populations in triple-negative breast cancer xenografts

Gene expression analysis at the single-cell level by next-generation sequencing has revealed the existence of clonal dissemination and microheterogeneity in cancer metastasis. The current spatial analysis technologies can elucidate the heterogeneity of cell–cell interactions in situ. To reveal the regional and expressional heterogeneity in primary tumors and metastases, we performed transcriptomic analysis of microtissues dissected from a triple-negative breast cancer (TNBC) cell line MDA-MB-231 xenograft model with our automated tissue microdissection punching technology. This multiple-microtissue transcriptome analysis revealed three cancer cell-type clusters in the primary tumor and axillary lymph node metastasis, two of which were cancer stem cell (CSC)-like clusters (CD44/MYC-high, HMGA1-high). Reanalysis of public single-cell RNA-sequencing datasets confirmed that the two CSC-like populations existed in TNBC xenograft models and in TNBC patients. The diversity of these multiple CSC-like populations could cause differential anticancer drug resistance, increasing the difficulty of curing this cancer.

Summary: We identified two types of cancer stem cell (CSC)-like populations in triple-negative breast cancer xenografts and patients. These CSC-like populations could potentially make tumors more drug resistant and thus more difficult to treat.

Identification of biomarkers of response to preoperative talazoparib monotherapy in treatment naïve gBRCA+ breast cancers

Germline mutations in BRCA1 or BRCA2 exist in ~2–7% of breast cancer patients, which has led to the approval of PARP inhibitors in the advanced setting. We have previously reported a phase II neoadjuvant trial of single agent talazoparib for patients with germline BRCA pathogenic variants with a pathologic complete response (pCR) rate of 53%. As nearly half of the patients treated did not have pCR, better strategies are needed to overcome treatment resistance. To this end, we conducted multi-omic analysis of 13 treatment naïve breast cancer tumors from patients that went on to receive single-agent neoadjuvant talazoparib. We looked for biomarkers that were predictive of response (assessed by residual cancer burden) after 6 months of therapy. We found that all resistant tumors exhibited either the loss of SHLD2, expression of a hypoxia signature, or expression of a stem cell signature. These results indicate that the deep analysis of pre-treatment tumors can identify biomarkers that are predictive of response to talazoparib and potentially other PARP inhibitors, and provides a framework that will allow for better selection of patients for treatment, as well as a roadmap for the development of novel combination therapies to prevent emergence of resistance.

A Cell Double-Barcoding System for Quantitative Evaluation of Primary Tumors and Metastasis in Animals That Uncovers Clonal-Specific Anti-Cancer Drug Effects

Simple Summary

The main problem in treating advanced cancers is a metastatic spread when individual cancer cells leave the primary tumor and colonize to distant organs. In drug development, it is important to quantitatively assess effects of novel drug candidates on both primary tumors and metastasis. Unfortunately, current methods of monitoring metastasis in mouse models have low sensitivity and are not quantitative. Here, we developed a methodology to monitor drug effects on metastasis that is quantitative and has a very high sensitivity and resolution. In fact, it allows monitoring effects of drugs on individual cancer cells in animals.

Abstract

Imaging in monitoring metastasis in mouse models has low sensitivity and is not quantitative. Cell DNA barcoding, demonstrating high sensitivity and resolution, allows monitoring effects of drugs on the number of tumor and metastatic clones. However, this technology is not suitable for comparison of sizes of metastatic clones in different animals, for example, drug treated and untreated, due to high biological and technical variability upon tumor and metastatic growth and isolation of barcodes from tissue DNA. However, both numbers of clones and their sizes are critical parameters for analysis of drug effects. Here we developed a modification of the barcoding approach for monitoring drug effects on tumors and metastasis that is quantitative, highly sensitive and highly reproducible. This novel cell double-barcoding system allows simultaneously following the fate of two or more cell variants or cell lines in xenograft models in vivo, and also following the fates of individual clones within each of these populations. This system allows comparing effects of drugs on different cell populations and thus normalizing drug effects by drug-resistant lines, which corrects for both biological and technical variabilities and significantly increases the reproducibility of results. Using this barcoding system, we uncovered that effects of a novel DYRK1B kinase inhibitor FX9847 on primary tumors and metastasis is clone-dependent, while a distinct drug osimertinib demonstrated clone-independent effects on cancer cell populations. Overall, a cell double-barcoding approach can significantly enrich our understanding of drug effects in basic research and preclinical studies.

Tissue architecture delineates field cancerization in BRAFV600E-induced tumor development

Cancer cells hijack developmental growth mechanisms but whether tissue morphogenesis and architecture modify tumorigenesis is unknown. Here, we characterized a new mouse model of sporadic thyroid carcinogenesis based on inducible expression of BRAF carrying a Val600 Glu (V600E) point mutation (BRAFV600E) from the thyroglobulin promoter (TgCreERT2). Spontaneous activation of this Braf-mutant allele due to leaky activity of the Cre recombinase revealed that intrinsic properties of thyroid follicles determined BRAF-mutant cell fate. Papillary thyroid carcinomas developed multicentrically within a normal microenvironment. Each tumor originated from a single follicle that provided a confined space for growth of a distinct tumor phenotype. Lineage tracing revealed oligoclonal tumor development in infancy and early selection of BRAFV600E kinase inhibitor-resistant clones. Somatic mutations were few, non-recurrent and limited to advanced tumors. Female mice developed larger tumors than males, reproducing the gender difference of human thyroid cancer. These data indicate that BRAFV600E-induced tumorigenesis is spatiotemporally regulated depending on the maturity and heterogeneity of follicles. Moreover, thyroid tissue organization seems to determine whether a BRAF-mutant lineage becomes a cancerized lineage. The TgCreERT2;BrafCA/+ sporadic thyroid cancer mouse model provides a new tool to evaluate drug therapy at different stages of tumor evolution.

Longitudinal Monitoring of Intra-Tumoural Heterogeneity Using Optical Barcoding of Patient-Derived Colorectal Tumour Models

Simple Summary

Colorectal cancer (CRC) is the second most common cancer worldwide. Despite improvements in the clinical management of CRC, outcomes for those with metastatic disease remain extremely poor. One reason for this is tumour heterogeneity, which refers to the observation that each cell within complex tumour cell populations displays different genetic features and biological behaviours. Such tumour heterogeneity is known to impact treatment efficacy and promote tumour recurrence. Here, we present a multi-colour barcoding methodology that allows for different lineages of colorectal cancer cells to be identified and monitored, thus allowing for tumour heterogeneity to be quantified in real-time. We show that discrete cell lineages can be quantified by both fluorescence microscopy and flow cytometry. Using this approach, we show that the cell culture models that are traditionally used in cancer research display limited heterogeneity, whereas patient-derived organoids—which are generated from fresh tumour resections—more faithfully represent the heterogeneity observed in cancer patients.

Abstract

Geno- and phenotypic heterogeneity amongst cancer cell subpopulations are established drivers of treatment resistance and tumour recurrence. However, due to the technical difficulty associated with studying such intra-tumoural heterogeneity, this phenomenon is seldom interrogated in conventional cell culture models. Here, we employ a fluorescent lineage technique termed “optical barcoding” (OBC) to perform simultaneous longitudinal tracking of spatio-temporal fate in 64 patient-derived colorectal cancer subclones. To do so, patient-derived cancer cell lines and organoids were labelled with discrete combinations of reporter constructs, stably integrated into the genome and thus passed on from the founder cell to all its clonal descendants. This strategy enables the longitudinal monitoring of individual cell lineages based upon their unique optical barcodes. By designing a novel panel of six fluorescent proteins, the maximum theoretical subpopulation resolution of 64 discriminable subpopulations was achieved, greatly improving throughput compared with previous studies. We demonstrate that all subpopulations can be purified from complex clonal mixtures via flow cytometry, permitting the downstream isolation and analysis of any lineages of interest. Moreover, we outline an optimized imaging protocol that can be used to image optical barcodes in real-time, allowing for clonal dynamics to be resolved in live cells. In contrast with the limited intra-tumour heterogeneity observed in conventional 2D cell lines, the OBC technique was successfully used to quantify dynamic clonal expansions and contractions in 3D patient-derived organoids, which were previously demonstrated to better recapitulate the heterogeneity of their parental tumour material. In summary, we present OBC as a user-friendly, inexpensive, and high-throughput technique for monitoring intra-tumoural heterogeneity in in vitro cell culture models.

Clonal populations of a human TNBC model display significant functional heterogeneity and divergent growth dynamics in distinct contexts

Intratumoral heterogeneity has been described for various tumor types and models of human cancer, and can have profound effects on tumor progression and drug resistance. This study describes an in-depth analysis of molecular and functional heterogeneity among subclonal populations (SCPs) derived from a single triple-negative breast cancer cell line, including copy number analysis, whole-exome and RNA sequencing, proteome analysis, and barcode analysis of clonal dynamics, as well as functional assays. The SCPs were found to have multiple unique genetic alterations and displayed significant variation in anchorage independent growth and tumor forming ability. Analyses of clonal dynamics in SCP mixtures using DNA barcode technology revealed selection for distinct clonal populations in different in vitro and in vivo environmental contexts, demonstrating that in vitro propagation of cancer cell lines using different culture conditions can contribute to the establishment of unique strains. These analyses also revealed strong enrichment of a single SCP during the development of xenograft tumors in immune-compromised mice. This SCP displayed attenuated interferon signaling in vivo and reduced sensitivity to the antiproliferative effects of type I interferons. Reduction in interferon signaling was found to provide a selective advantage within the xenograft microenvironment specifically. In concordance with the previously described role of interferon signaling as tumor suppressor, these findings suggest that similar selective pressures may be operative in human cancer and patient-derived xenograft models.

Subclonal heterogeneity and evolution in breast cancer

Subclonal heterogeneity and evolution are characteristics of breast cancer that play a fundamental role in tumour development, progression and resistance to current therapies. In this review, we focus on the recent advances in understanding the epigenetic and transcriptomic changes that occur within breast cancer and their importance in terms of cancer development, progression and therapy resistance with a particular focus on alterations at the single-cell level. Furthermore, we highlight the utility of using single-cell tracing and molecular barcoding methodologies in preclinical models to assess disease evolution and response to therapy. We discuss how the integration of single-cell profiling from patient samples can be used in conjunction with results from preclinical models to untangle the complexities of this disease and identify biomarkers of disease progression, including measures of intra-tumour heterogeneity themselves, and how enhancing this understanding has the potential to uncover new targetable vulnerabilities in breast cancer.

Evolution of metastasis: new tools and insights

Metastasis is an evolutionary process occurring across multiple organs and timescales. Due to its continuous and dynamic nature, this multifaceted process has been challenging to investigate and remains incompletely understood, in part due to the lack of tools capable of probing genomic evolution at high enough resolution. However, technological advances in genetic sequencing and editing have provided new and powerful methods to refine our understanding of the complex series of events that lead to metastatic dissemination. In this review, we summarize the latest genetic and lineage-tracing approaches developed to unravel the genetic evolution of metastasis. The findings that have emerged have enhanced our comprehension of the mechanistic trajectories and timescales of metastasis and could provide new strategies for therapy.

The Origins of Phenotypic Heterogeneity in Cancer

Heterogeneity is a pervasive feature of cancer, and understanding the sources and regulatory mechanisms underlying heterogeneity could provide key insights to help improve the diagnosis and treatment of cancer. In this review, we discuss the origin of heterogeneity in the phenotype of individual cancer cells. Genotype-phenotype (G-P) maps are widely used in evolutionary biology to represent the complex interactions of genes and the environment that lead to phenotypes that impact fitness. Here, we present the rationale of an extended G-P (eG-P) map with a cone structure in cancer. The eG-P cone is formed by cells that are similar at the genome layer but gradually increase variability in the epigenome, transcriptome, proteome, metabolome and signalome layers to produce large variability at the phenome layer. Experimental evidence from single-cell -omics analyses supporting the cancer eG-P cone concept is presented, and the impact of epimutations and the interaction of cancer and tumor microenvironmental eG-P cones are integrated with the current understanding of cancer biology. The eG-P cone concept uncovers potential therapeutic strategies to reduce cancer evolution and improve cancer treatment. More methods to study phenotypes in single cells will be key to better understand cancer cell fitness in tumor biology and therapeutics.

Single-cell RNA sequencing reveals cell heterogeneity and transcriptome profile of breast cancer lymph node metastasis

Molecular mechanisms underlying breast cancer lymph node metastasis remain unclear. Using single-cell sequencing, we investigated the transcriptome profile of 96,796 single cells from 15 paired samples of primary tumors and axillary lymph nodes. We identified nine cancer cell subclusters including CD44 + / ALDH2 + /ALDH6A1 + breast cancer stem cells (BCSCs), which had a copy-number variants profile similar to that of normal breast tissue. Importantly, BCSCs existed only in primary tumors and evolved into metastatic clusters infiltrating into lymph nodes. Furthermore, transcriptome data suggested that NECTIN2-TIGIT-mediated interactions between metastatic breast cancer cells and tumor microenvironment (TME) cells, which promoted immune escape and lymph node metastasis. This study is the first to delineate the transcriptome profile of breast cancer lymph node metastasis using single-cell RNA sequencing. Our findings offer novel insights into the mechanisms underlying breast cancer metastasis and have implications in developing novel therapies to inhibit the initiation of breast cancer metastasis.

Mammary-specific expression of Trim24 establishes a mouse model of human metaplastic breast cancer

Conditional overexpression of histone reader Tripartite motif containing protein 24 (TRIM24) in mouse mammary epithelia (Trim24COE) drives spontaneous development of mammary carcinosarcoma tumors, lacking ER, PR and HER2. Human carcinosarcomas or metaplastic breast cancers (MpBC) are a rare, chemorefractory subclass of triple-negative breast cancers (TNBC). Comparison of Trim24COE metaplastic carcinosarcoma morphology, TRIM24 protein levels and a derived Trim24COE gene signature reveals strong correlation with human MpBC tumors and MpBC patient-derived xenograft (PDX) models. Global and single-cell tumor profiling reveal Met as a direct oncogenic target of TRIM24, leading to aberrant PI3K/mTOR activation. Here, we find that pharmacological inhibition of these pathways in primary Trim24COE tumor cells and TRIM24-PROTAC treatment of MpBC TNBC PDX tumorspheres decreased cellular viability, suggesting potential in therapeutically targeting TRIM24 and its regulated pathways in TRIM24-expressing TNBC.

Aiduqing formula inhibits breast cancer metastasis by suppressing TAM/CXCL1-induced Treg differentiation and infiltration

Background

Metastasis represents the leading cause of death in patients with breast cancer. Traditional Chinese medicine is particularly appreciated for metastatic diseases in Asian countries due to its benefits for survival period prolongation and immune balance modulation. However, the underlying molecular mechanisms remain largely unknown. This study aimed to explore the antimetastatic effect and immunomodulatory function of a clinical formula Aiduqing (ADQ).

Methods

Naive CD4⁺ T cells, regulatory T cells (Tregs), and CD8⁺ T cells were sorted by flow cytometry. Then, breast cancer cells and these immune cells were co-cultured in vitro or co-injected into mice in vivo to simulate their coexistence. Flow cytometry, ELISA, qPCR, double luciferase reporter gene assay, and chromatin immunoprecipitation assay were conducted to investigate the immunomodulatory and antimetastatic mechanisms of ADQ.

Results

ADQ treatment by oral gavage significantly suppressed 4T1-Luc xenograft growth and lung metastasis in the orthotopic breast cancer mouse model, without noticeable hepatotoxicity, nephrotoxicity, or hematotoxicity. Meanwhile, ADQ remodeled the immunosuppressive tumor microenvironment (TME) by increasing the infiltration of tumor-infiltrating lymphocytes (TILs) and cytotoxic CD8⁺ T cells, and decreasing the infiltration of Tregs, naive CD4⁺ T cells, and tumor-associated macrophages (TAMs). Molecular mechanism studies revealed that ADQ remarkably inhibited CXCL1 expression and secretion from TAMs and thus suppressed the chemotaxis and differentiation of naive CD4⁺ T cells into Tregs, leading to the enhanced cytotoxic effects of CD8⁺ T cells. Mechanistically, TAM-derived CXCL1 promoted the differentiation of naive CD4⁺ T cells into Tregs by transcriptionally activating the NF-κB/FOXP3 signaling. Lastly, mouse 4T1-Luc xenograft experiments validated that ADQ formula inhibited breast cancer immune escape and lung metastasis by suppressing the TAM/CXCL1/Treg pathway.

Conclusions

This study not only provides preclinical evidence supporting the application of ADQ in inhibiting breast cancer metastasis but also sheds novel insights into TAM/CXCL1/NF-κB/FOXP3 signaling as a promising therapeutic target for Treg modulation and breast cancer immunotherapy.

Resveratrol attenuates TNBC lung metastasis by down-regulating PD-1 expression on pulmonary T cells and converting macrophages to M1 phenotype in a murine tumor model

Triple-negative breast cancer (TNBC) is an invasive breast cancer with the characteristics of easy to develop distant metastasis. Immune escape is one of the main reasons for TNBC growth and metastasis. Enhancement of T cell-mediated anti-tumor activity may benefit to inhibit tumor metastasis and improve the efficacy of cancer therapy. As a natural bioactive substance, resveratrol shows potential capability to prevent or suppress the development of a variety of cancers through direct or indirect effects, including immunoregulatory effect. However, whether resveratrol might affect lung metastasis of TNBC, and whether the effect of resveratrol might be associated with resveratrol-regulated immune responses in tumor microenvironment is still unknown. In this study, by using an experimental metastatic mouse 4 T1 tumor model, we identified that resveratrol may suppress TNBC lung metastasis by elevating local anti-tumor immunity. Indeed, an increase in the cytotoxic activity of CD8⁺T cells as well as the levels of type 1 cytokine IFN-γ and IL-2 in the lungs of resveratrol-treated tumor bearing mice were observed. The enhanced CD8⁺T cell activity and Th1 immune responses by resveratrol administration might be related to the down-regulated PD-1 expression on pulmonary CD8⁺T cells and CD4⁺T cells. Resveratrol may also convert macrophages to M1 phenotype in the lungs of tumor bearing mice. However, it seems likely resveratrol has no effect on pulmonary myeloid-derived suppressor cell activation. Our results provide an evidence that resveratrol might be a promising candidate agent for adjuvant therapy in the process of TNBC metastasis.

Spatial omics and multiplexed imaging to explore cancer biology

Understanding intratumoral heterogeneity-the molecular variation among cells within a tumor-promises to address outstanding questions in cancer biology and improve the diagnosis and treatment of specific cancer subtypes. Single-cell analyses, especially RNA sequencing and other genomics modalities, have been transformative in revealing novel biomarkers and molecular regulators associated with tumor growth, metastasis and drug resistance. However, these approaches fail to provide a complete picture of tumor biology, as information on cellular location within the tumor microenvironment is lost. New technologies leveraging multiplexed fluorescence, DNA, RNA and isotope labeling enable the detection of tens to thousands of cancer subclones or molecular biomarkers within their native spatial context. The expeditious growth in these techniques, along with methods for multiomics data integration, promises to yield a more comprehensive understanding of cell-to-cell variation within and between individual tumors. Here we provide the current state and future perspectives on the spatial technologies expected to drive the next generation of research and diagnostic and therapeutic strategies for cancer.

Collective metastasis: coordinating the multicellular voyage

The metastatic process is arduous. Cancer cells must escape the confines of the primary tumor, make their way into and travel through the circulation, then survive and proliferate in unfavorable microenvironments. A key question is how cancer cells overcome these multiple barriers to orchestrate distant organ colonization. Accumulating evidence in human patients and animal models supports the hypothesis that clusters of tumor cells can complete the entire metastatic journey in a process referred to as collective metastasis. Here we highlight recent studies unraveling how multicellular coordination, via both physical and biochemical coupling of cells, induces cooperative properties advantageous for the completion of metastasis. We discuss conceptual challenges and unique mechanisms arising from collective dissemination that are distinct from single cell-based metastasis. Finally, we consider how the dissection of molecular transitions regulating collective metastasis could offer potential insight into cancer therapy.

Breast cancer as an example of tumour heterogeneity and tumour cell plasticity during malignant progression

Heterogeneity within a tumour increases its ability to adapt to constantly changing constraints, but adversely affects a patient's prognosis, therapy response and clinical outcome. Intratumoural heterogeneity results from a combination of extrinsic factors from the tumour microenvironment and intrinsic parameters from the cancer cells themselves, including their genetic, epigenetic and transcriptomic traits, their ability to proliferate, migrate and invade, and their stemness and plasticity attributes. Cell plasticity constitutes the ability of cancer cells to rapidly reprogramme their gene expression repertoire, to change their behaviour and identities, and to adapt to microenvironmental cues. These features also directly contribute to tumour heterogeneity and are critical for malignant tumour progression. In this article, we use breast cancer as an example of the origins of tumour heterogeneity (in particular, the mutational spectrum and clonal evolution of progressing tumours) and of tumour cell plasticity (in particular, that shown by tumour cells undergoing epithelial-to-mesenchymal transition), as well as considering interclonal cooperativity and cell plasticity as sources of cancer cell heterogeneity. We review current knowledge on the functional contribution of cell plasticity and tumour heterogeneity to malignant tumour progression, metastasis formation and therapy resistance.

The site of breast cancer metastases dictates their clonal composition and reversible transcriptomic profile

Intratumoral heterogeneity is a driver of breast cancer progression, but the nature of the clonal interactive network involved in this process remains unclear. Here, we optimized the use of optical barcoding to visualize and characterize 31 cancer subclones in vivo. By mapping the clonal composition of thousands of metastases in two clinically relevant sites, the lungs and liver, we found that metastases were highly polyclonal in lungs but not in the liver. Furthermore, the transcriptome of the subclones varied according to their metastatic niche. We also identified a reversible niche-driven signature that was conserved in lung and liver metastases collected during patient autopsies. Among this signature, we found that the tumor necrosis factor-α pathway was up-regulated in lung compared to liver metastases, and inhibition of this pathway affected metastasis diversity. These results highlight that the cellular and molecular heterogeneity observed in metastases is largely dictated by the tumor microenvironment.

Treatment scheduling effects on the evolution of drug resistance in heterogeneous cancer cell populations

The effect of scheduling of targeted therapy combinations on drug resistance is underexplored in triple-negative breast cancer (TNBC). TNBC constitutes heterogeneous cancer cell populations the composition of which can change dynamically during treatment resulting in the selection of resistant clones with a fitness advantage. We evaluated crizotinib (ALK/MET inhibitor) and navitoclax (ABT-263; Bcl-2/Bcl-xL inhibitor) combinations in a large design consisting of 696 two-cycle sequential and concomitant treatment regimens with varying treatment dose, duration, and drug holiday length over a 26-day period in MDA-MB-231 TNBC cells and found that patterns of resistance depend on the schedule and sequence in which the drugs are given. Further, we tracked the clonal dynamics and mechanisms of resistance using DNA-integrated barcodes and single-cell RNA sequencing. Our study suggests that longer formats of treatment schedules in vitro screening assays are required to understand the effects of resistance and guide more realistically in vivo and clinical studies.

Long non-coding RNA CASC9/microRNA-590-3p axis participates in lutein-mediated suppression of breast cancer cell proliferation

Previous studies have shown that lutein can inhibit the proliferation of breast cancer cells. However, the mechanism of lutein inhibiting the proliferation of breast cancer cells remains unclear. The present study aimed to determine whether the long non-coding RNA (lncRNA) Cancer Susceptibility 9 (CASC9)/microRNA (miR)-590-3p axis participates in the antiproliferative effects of lutein via lncRNA microarray hybridization, reverse transcription-quantitative PCR, dual-luciferase reporter and MTT assays. The results demonstrated that CASC9 was the most significantly downregulated lncRNA in MCF7 cells treated with lutein. miR-590-3p was identified as the target of CASC9. In addition, lutein downregulated CASC9 expression and upregulated miR-590-3p expression in dose- and time-dependent manners, respectively. CASC9 knockdown or overexpression of miR-590-3p inhibited the proliferation of breast cancer cells. Notably, simultaneous transfection with miR-590-3p mimics and CASC9 small interfering RNA increased the potency of lutein in inhibiting the proliferation of breast cancer cells. Taken together, these results suggest that the CASC9/miR-590-3p axis participates in the antiproliferative effects of lutein on breast cancer.