Skip to main content

Breast cancer heterogeneity and its implication in personalized precision therapy

A Correction to this article was published on 23 January 2024

This article has been updated

Abstract

Breast cancer heterogeneity determines cancer progression, treatment effects, and prognosis. However, the precise mechanism for this heterogeneity remains unknown owing to its complexity. Here, we summarize the origins of breast cancer heterogeneity and its influence on disease progression, recurrence, and therapeutic resistance. We review the possible mechanisms of heterogeneity and the research methods used to analyze it. We also highlight the importance of cell interactions for the origins of breast cancer heterogeneity, which can be further categorized into cooperative and competitive interactions. Finally, we provide new insights into precise individual treatments based on heterogeneity.

Introduction

Breast cancer is the most common cancer among women worldwide, with increasing incidence. Owing to its heterogeneity [1], it is also the second leading cause of cancer-related deaths in women [2, 3]. The heterogeneity of breast cancer is attributed to differences in the genomic, epigenomic, transcriptomic, and proteomic characteristics of the cancer cells. These factors affect tumor properties such as proliferation, apoptosis, metastasis, and therapeutic response. This heterogeneity is also observed in tumor tissues among different patients or different metastases (intertumor heterogeneity) and within the individual tumor from the same patient (intratumor heterogeneity) (Fig. 1) [4]. The heterogeneity is also spatiotemporal. Temporal heterogeneity refers to the dynamic variations in the molecular makeup of cancer cells during tumor progression [5]. Spatial heterogeneity describes the distribution and interactions of different cell populations in complex tumor structures. Several scholars believe breast cancer is the most heterogeneous cancer type in women worldwide [6]. Breast cancer heterogeneity increases the difficulty of early diagnosis, treatment selection, and prognosis prediction. For instance, in situ recurrence, distant metastasis, and many clinical problems in breast cancer treatment arise from the heterogeneity [1, 6,7,8]. In this review, we summarize the different levels of heterogeneity, their influence on breast cancer progression, the underlying mechanisms, new research methods to measure breast cancer heterogeneity, and treatment strategies in terms of heterogeneity.

Fig. 1
figure 1

Inter- and intratumor heterogeneity of breast cancer. A Intertumoral heterogeneity: a Variation in different breast cancer subtypes among different patients; b Variation in different breast cancer subtypes in patients with primary breast cancer and its metastasis in the same patient. B Intratumor heterogeneity: variation in different cell types within a tumor

Different levels of heterogeneity

The most fundamental manifestations and causes of heterogeneity of breast cancer are genomic, epigenomic, transcriptomic, and proteomic alterations expressed in spatial and temporal dimensions. At the genomic level, the heterogeneity comprises of mutations [9] and copy number aberrations (CNAs) [10]. At the epigenomic level, the heterogeneity comprises of transposase-accessible chromatin [10], DNase [11], nucleosome occupancy [12], methylome [12], and so on. At the transcriptomic level, the heterogeneity refers to the stochastic gene expression [9]. At the proteomic level, the heterogeneity comprises of protein modifications and signal transduction [9]. The Cancer Genome Atlas [13], Molecular Taxonomy of Breast Cancer International Consortium [14], and more recently, single-cell profiling [15,16,17,18,19] have revealed intertumor and intratumor heterogeneity also includes temporal and spatial variations. Changes in the genome, epigenome, transcriptome, and proteome heterogeneity can be reflected in various aspects, such as biomarkers, metabolism, cell cycle, tumor microenvironment (TME), epithelial–mesenchymal transition (EMT), circulating tumor cells (CTCs), and clinical pathology.

Biomarkers

Regarding biomarker heterogeneity, the American Society of Clinical Oncology/College of American Pathologists recommends detection of the expression levels of estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth factor receptor 2 (HER2) in all invasive breast carcinomas by immunohistochemistry [20, 21]. ER and PR are often referred to collectively as hormone receptors (HRs). ER is produced by three genes, named ERα, ERβ, and membrane receptor G protein-coupled receptor 30 (GPR30) [22]. PR is produced by two genes, PRA and PRB [23]. HER2 is produced by ERBB2 [24]. Based on these biomarkers, breast cancers are typically classified into four subtypes: luminal A (ER+, PR+, and HER2), luminal B (ER+, PR+, and HER2+), HER2-enriched (ER, PR, and HER2+), and triple-negative breast cancers (TNBC) (ER, PR, and HER2) [25,26,27,28]. According to Fudan University Shanghai Cancer Center classification criteria, TNBC can further be categorized into immunomodulatory (IM), luminal androgen receptor (LAR), mesenchymal-like (MES), basal-like and immune-suppressed (BLIS) subtypes [29]. Recently, the claudin-low subtype has been defined as mesenchymal cell-like, stem cell-like, and has high expression of immune response genes [30]. It was demonstrated that the claudin-low subtype is an additional complex phenotype that may permeate breast tumors of different subtypes [31]. Furthermore, HR+/HER2+ tumors behave differently, with some resembling more luminal A subtype and others resembling HER2-enriched subtype [32]. Some scholars believe that it should be considered an individual subtype. Breast cancer treatment is often based on the ER, PR, and HER2 status in primary cancer. However, increasing data have shown substantial discordance between the primary sites or the receptor profiles of primary breast cancers and their distant metastases [33,34,35,36,37,38,39,40,41]. This discordance can also be observed between different metastatic sites [42]. This conversion may reflect sampling or treatment effects; however, they may also indicate an alteration in cancer biology and heterogeneity.

Metabolism

Recent studies have demonstrated that metabolic levels differ among cells in a tumor mass, resulting from differences in various metabolic pathways [43]. The degree of biomarker positivity is negatively correlated with the intensity of metabolism and aggressive clinical behavior [44, 45]. For example, TNBC, with a more aggressive clinical course, showed increased demand for nicotinamide, 1-ribosyl-nicotinamide, and NAD+ than other types of breast cancer. Compared with primary cancer cells, metastatic cells shift their metabolic levels with lower glucose metabolism [46]. In contrast, metastatic cells have higher levels of fumarate, malate, and succinyl-carnitine [46]. Furthermore, metastatic cells have higher fatty acid and amino acid metabolism than primary cancer cells [46].

Cell cycle

An uncontrolled cell cycle is a hallmark of breast cancer [47]. Nevertheless, the cell cycle progression within breast cancer cells varies considerably [48]. Using cell cycle phase analysis, researchers have found that most cancer cells in tumors do not undergo the conventional cell cycle phases [48] and are resistant to treatment. Compared with other subtypes of breast cancer, TNBC show a higher proportion of mitoses in the overall cycling cells [49]. Furthermore, tumor size was demonstrated to positively affect the proportion of proliferated cells.

TME

The TME plays a significant role in tumor growth and metastasis and has always been an important obstacle during tumor therapy. The TME has many heterogeneous cell populations, including fibroblasts, adipocytes, immune cells, epithelial cells, pericytes, and extracellular matrix components [50]. Interactions between tumor cells and their microenvironment contribute to tumor heterogeneity, thus promoting tumor progression [51].

EMT, CSCs, and CTCs

Generally, EMT is regarded as the first and foremost step in the cancer metastasis [52, 53]. It is considered a heterogeneous phenomenon [54] that can also be involved in embryonic development and wound healing. EMT correlates with cancer stem cell (CSC) plasticity and therapeutic resistance [55, 56]. Increasing evidence indicates that EMT promotes breast cancer metastasis [57]. CTCs represent the intermediate stage of metastasis [58]. CTCs are found in the peripheral blood of breast cancer patients with metastatic lesions [58]. They can be derived from primary and metastatic sites [58]. CTCs exhibit high inter- and intra-patient heterogeneity compared with normal breast cancer cells [59]. Therefore, the analysis of CTCs can be an essential tool for studying the heterogeneity between cells from the primary and metastatic sites. Because CTCs can be obtained by minimally invasive peripheral blood extraction, they can be collected continuously to monitor molecular changes as the disease progresses and, significantly, under therapeutic pressure [60]. Thus, analysis of a single CTC enhances the assessment of heterogeneity.

Microcirculation

Heterogeneity in tumor microcirculation is caused by the combination or addition of vessels, angiogenic vessels, mosaic vessels, and vasogenic mimicry [61]. As such, antiangiogenic agents alone can lead to poor efficacy and resistance. Vasogenic mimicry (VM) is a vascular structure formed by cancer cells that can transport blood and cancer cells to efficiently obtain oxygen and nutrients independent of endothelial cells [62, 63]. The presence of VM in breast cancer is a heterogeneous phenomenon. It was demonstrated that, compared with other subtypes, VM tends to form in TNBC [64, 65]. Furthermore, VM is associated with high tumor grade, tumor progression, invasion, metastasis, and poor prognosis [63]. Tubular type and patterned matrix type are two types of VM [63]. The tubular type may be morphologically confused with vessels lined with endothelial cells. Patterned matrix type is not morphologically or topologically similar to vessels at all [61], and is characteristic of highly invasive tumors in which therapeutic regimens that target angiogenesis may be ineffective, resulting in recurrence and metastasis [61].

Clinical pathology

In a clinical and histopathologic evaluation of heterogeneity, the American Joint Committee on Cancer incorporated a tumor, lymph node, metastasis (TNM) staging system [66]. TNM is scored based on tumor volume, the status of nearby lymph nodes, and the presence of distant metastasis to judge the heterogeneity and guide the prognosis and treatment of patients [67]. The grading of breast cancer also highlights its heterogeneity. Grading is assessed according to three morphological parameters: the proportion of carcinoma in glands and tubular structures, degree of nuclear pleomorphism, and mitosis rate in three hierarchical (low, medium, and high) systems [68].

Breast cancer heterogeneity affects disease progression

Therapeutic resistance, recurrence, and metastasis are the most severe problems in the breast cancer treatment [69]. These problems cannot be completely separated from each other. Chemotherapy can kill sensitive cancer cells; however, the resistant cells left behind can cause tumor recurrence (or relapse) and even cancer metastasis [70].

Therapeutic resistance

In cancer, genomic instability results in a significant degree of intercellular genetic heterogeneity. According to the tumor stem cell theory, the original genetic variation can lead to changes in epigenetic factors [70] and the cell cycle. Gene products that are directly involved in the development of drug resistance may be increased or diminished following changes in these factors [70]. Therefore, the mechanisms underlying the heterogeneity must be studied to address therapeutic resistance. Neoadjuvant chemotherapy for many patients with TNBC, is effective in approximately half of all cases [71, 72]. The genomic and molecular basis of chemotherapy resistance in patients with TNBC is still poorly understood. Kim et al. applied single-cell DNA and RNA sequencing in addition to bulk exome sequencing of samples from 20 patients with TNBC during NAC [16] and discovered a new chemoresistance model containing two modes of evolution. They found the preexistence of therapy-resistant genotypes in these patients prior to NAC. Furthermore, during treatment with NAC, a small portion of genotypes selected by NAC were subjected to transcriptional reprogramming. In this model, two modes of evolution (adaptive and acquired) were in operation, resulting in therapy resistance and recurrence. Similar findings have been observed in ER+ breast cancer. Marta et al. performed extensive single-cell gene expression profiling of breast cancer cell line MCF7 and its derivatives, including docetaxel-resistant cells [73]. They found untreated parental cells already contained a small portion of cells in which EMT and stemness-related genes were upregulated and cell cycle-related genes were downregulated, similar to drug-resistant cells. During chemotherapy, this population may be positively selected, leading to treatment failure. In addition, epigenetic reprogramming, aberrant co-factor activity, hypoxia, stromal factors, and inflammatory and immune components are thought to contribute to the resistance to endocrine therapies.

Radiation therapy is an adjunctive therapy for many primary cancers and one of the most commonly used methods for breast cancer treatment [74]. Despite the use of advanced radiotherapy methods, radiotherapy accuracy, and fractionation schemes, radiation therapy resistance is frequently observed [75, 76] in the form of tumor progression and recurrence. Solid tumors, such as breast and lung tumors, contain markedly heterogeneous CSC populations [77]. Radiation therapy resistance in breast tumors is highly correlated with the proportion of CSCs [78]. A higher proportion of breast CSCs (BCSCs) leads to higher radioresistance and lower biological efficiency of the tumor.

As mentioned above, immune cells as a part of the TME, its subtype, percentage, and distribution vary significantly in breast cancer masses [79]. “Immune subtypes” of breast cancer can be defined as neutrophil-enriched subtypes (NES) and macrophage-enriched subtypes (MES) [79]. MES mainly consists of C–C chemokine receptor type 2 (CCR2)-dependent macrophages that respond to immune checkpoint blockade (ICB) therapy. NES contains widespread or localized immunosuppressive neutrophils (or granulocytic myeloid-derived suppressor cells) resistant to ICB therapy and includes a small number of macrophages unaffected by CCR2 knockout. Converting MES to NES enables an initially sensitive MES model to gain ICB resistance [79].

Recurrence and metastasis

Owing to heterogeneity, small portions of breast cancer cells resistant to chemotherapy, radiation therapy, and immune therapy are selected and survive during treatment, as discussed above. The survival of resistant cancer cells leads to treatment failure and recurrence after treatment. BCSCs are highly correlated to the recurrence of breast cancer patients even after surgery [80]. BCSCs are highly heterogeneous, and patients with higher marker expression of 74 genes, which have since been revealed by single-cell RNA sequencing, in BCSCs were more likely to show recurrence than patients with lower marker expression [81, 82]. However, different subtypes of breast cancer exhibit various recurrence risks. TNBC has a higher recurrence rate than other types of breast cancer, partially due to a lack of effective targeting strategy, especially in the first five years [83]. Tumor-infiltrating lymphocytes (TILs) are associated with survival benefits in HER2-positive and TNBC; however, in luminal/HER2-negative tumors, high TILs are negatively associated with survival benefits. Thus, it is challenging to determine the prognosis of treatment and the recurrence risk owing to heterogeneity.

More than 90% of breast cancer-related deaths result from metastasis [13, 84, 85]. The 5-year overall survival rate of patients with breast cancer without metastasis is greater than 80% [86, 87]. Circulating breast cancer cells can be divided into various cell types based on differences in biomarker status, epithelial/mesenchymal phenotype, aggregation status, and other factors [88, 89]. RNA sequencing has revealed that weakly migratory subpopulations are largely epithelial and highly migratory subpopulations are largely mesenchymal [88]. Interestingly, the metastases of weakly migratory subpopulations are significantly more frequent than those of highly migratory subpopulations when injected in situ [88]. Hapach et al. found that weakly migratory cells gather as CTC clusters in the circulation rather than remain as a single CTC. The higher expression of E-cadherin also contributed to increased metastasis [88]. CTCs may exploit their genetic potential and communicate with the surrounding environment, such as chemokine systems, blood cells, and extracellular matrix components, to regulate organ-specific metastasis of breast cancer [89]. VM may also play an essential role in metastasis. It was found that patients with VM tended to have a higher rate of hematogenous recurrence and a lower 5-year survival rate. By inhibiting VM, primary tumor growth and lung metastasis are reduced [90].

Cell interactions contribute to breast cancer heterogeneity

Previous research has mainly focused on the intracellular autonomic changes that lead to abnormal activity of cancer cells, which is a hallmark of cancer. However, there is increasing recognition of how non-cell-autonomous factors influence cancer cell survival, growth, and proliferation and their function and contribution to cancer. An increasing number of studies have reported competitive or cooperative interactions between cells in recent years.

In some cases, breast cancer subclones may interact to gain a selective growth advantage, and interclonal cooperation is essential for cancer maintenance [91]. Cooperation between different cancer cell subtypes may also result in a malignant phenotype [92]. In contrast, more aggressive clones and subclones can eliminate less aggressive and normal host cells through competitive interactions. Cellular competition refers to the survival of the fittest at the cellular level. The winner of the cell competition can identify and eliminate defective or damaged cells with low fitness characteristics (called "losers"), thus creating growth space and nutrients for the growth of cancer cells. Cellular competition may inhibit early carcinogenesis by eliminating mutant cells from standard wild-type cells [93,94,95]. Competition between cancer cells and the TME is one of the mechanisms of tumor invasion, diffusion, evolution, and formation of the intratumor [96,97,98] and intertumor heterogeneity [97]. This phenomenon has been linked to several cancer-related genes and thus may play an essential role in cancer. Cancer hypoxia, clonal stem cell selection, and immune cell response are thought to have a potential connection with the mechanisms that occur during cell competition [99].

Notably, cancer initiation and progression are not static but highly dynamic processes; therefore, subclones of cancer cells and cell interactions change over time [100]. Changing cancer cell interactions eventually produce heterogeneous cancers with various genotypes and phenotypes [101, 102] (Fig. 2). Therefore, how cell-to-cell interactions shape cancer populations in both space and time and to what extent cell interactions actively shape cancer dynamics require further exploration.

Fig. 2
figure 2

Schematic of breast cancer microenvironment. Breast cancer cells infiltrate tissues and interact with various cell types

Breast cancer cells—fibroblasts

Cancer-associated fibroblasts (CAFs) are specialized fibroblasts in the cancer microenvironment [103] derived from various cells. CAFs can stimulate cancer growth and progression by secreting stromal cell-derived factor 1 (SDF-1)/C-X-C motif ligand 12 (CXCL12) [104, 105], transforming growth factor-β (TGF-β) [105], interleukin (IL)-8 [106], and IL-32 [107]. CAFs harboring activated hypoxia-inducible factor-α (HIF-α) can promote cancer growth and lymph node metastasis via the paracrine production of recycled nutrients to feed cancer cells [108]. After coculture with CAFs, the expression levels of S100 calcium-binding protein A4 (S100A4), TGF-β, fibroblast growth factor (FGF)2, FGF7, matrix metalloproteinase (MMP)-2, MMP9, MMP11, platelet-derived growth factor (PDGF) A, PDGFB, vascular endothelial growth factor A (VEGFA), IL-6, IL-8, urokinase-type plasminogen activator (uPA), and tissue inhibitor matrix metalloproteinase 1 (TIMP1) in breast cancer cells are increased [109]. TGFβ secreted by breast cancer cells can stimulate bone morphogenetic protein (BMP) antagonist GREMLIN1 (GREM1) in CAFs [110]. GREM1 further promotes mesenchymal phenotype, stemness, and invasion by disrupting BMP/SMAD signaling in breast cancer cells [110]. CAFs can secrete exosomes containing microRNA(miR)-21, -378e, -143 [111], -16, -148a [112], and -500a-5p [113], which are transferred to breast cancer cells and subsequently promote the proliferation and metastasis of cancer cells. Interestingly, four CAF subpopulations in metastatic lymph nodes play different roles.

CAF-S1 activates the EMT of breast cancer cells through CXCL12 and TGFβ pathways, mainly playing a role in stimulating cancer cell migration. At the same time, CAF-S4 induces the invasion of breast cancer cells through NOTCH signaling in a three-dimensional space [105]. Furthermore, collagen type I secreted by CAF is related to a decrease in intratumor chemotherapeutic drug uptake, which results in resistance to chemotherapy [114, 115]. CAFs can function as a niche, allowing diffuse cancer cells to metastasize and colonize distant organs [116]. Furthermore, CAFs can remodel the extracellular matrix and immune microenvironment, inducing cancer cell resistance to immune checkpoint inhibitors (ICIs) and leading to refractory cancer [117,118,119,120]. Interestingly, breast cancer cells can induce the transformation of fibroblasts into CAFs, which in turn recruit breast cancer cells to promote cancer progression [121]. Thus, cancer cells interacting with CAFs could lead to resistance to anticancer chemotherapy and disease progression, contributing to breast cancer heterogeneity.

High-throughput live-cell imaging has revealed that normal human fibroblasts can suppress cancer cells to varying degrees [122]. The expression of phosphatase and tensin homolog (PTEN) in normal fibroblasts can suppress epithelial mammary tumors [123] through the Pten–miR-320–Ets2 axis [123, 124]. Furthermore, transplanting normal fibroblasts from a normal breast into a xenograft model of progressive proliferative breast disease can inhibit cellular invasiveness and tumor progression [125] through enhanced apoptosis. TNBC cells can provide integrin β4 (ITGB4) proteins via exosomes to induce mitophagy in the CAFs [126]. ITGB4-induced mitophagy in CAFs can promote the conversion of pyruvate to lactate to feed breast cancer cells, contributing to TNBC progression [126].

Breast cancer cells—adipocytes

Cancer-associated adipocytes (CAAs), a main component of the TME, modify cancer cell phenotype by interacting with tumor cells [127, 128]. Moreover, leptin and adiponectin produced by adipocytes can increase the proliferation and invasion of cancer cells [129]. Adipocyte-derived IL-6 and leptin promote breast cancer metastasis by activating the Janus kinase (JAK)/signal transducer and activator of transcription 3 (STAT3) and phosphoinositide 3-kinase (PI3K)/protein kinase B (AKT) signaling pathways, thus upregulating lysyl hydroxylase (PLOD2) expression [130]. After the activation of PI3K/AKT in the breast cancer cells, the expression of IL-6, IL-1β, and tumor necrosis factor-α (TNF-α) in breast cancer cells is increased [131]. Leptin can also promote plasminogen activator inhibitor-1 expression, which is related to EMT, by activating STAT3 to promote metastasis [132]. Free fatty acids (FFAs) [133], monocyte chemoattractant protein-1 (MCP-1) [134], CC-chemokine ligand 5(CCL5) [135], and insulin-like growth factor 1 (IGF-1) [136] secreted by CAAs can induce breast cancer cells to be more invasive or more proliferative. Interestingly, higher expression of HIF-1α in adipocytes causes the loss of the ERα protein [137], which may play a role in receptor conversion and therapy resistance. In addition, a coculture study showed that adipocytes could secrete adipocyte soluble factors, which can promote higher major vault protein (MVP) expression of breast cancer cells to reduce drug accumulation in the nuclei, and efflux into the extracellular medium through extracellular vesicle secretion, thus inducing a multidrug-resistant phenotype [115]. Moreover, higher MVP expression is found at the invasive front in cancer cells closer to adipocytes, suggesting that MVP-mediated drug resistance may be clinically relevant [115]. The interactions between cancer cells and adipocytes affect cancer progression and heterogeneity.

Breast cancer cells—immune/inflammatory cells

Inflammatory processes are associated with the initiation and progression of breast cancer, which often presents with many immune cells [138]. The infiltration of immune cells in tumors may be heterogeneous, with significant differences in immune cell types and protein composition among patients and tumor stages [139]. Tumor-associated macrophages (TAMs) are essential in TME and are involved in carcinogenesis, progression, and therapy resistance in the cancer [140,141,142,143]. IL-1β, IL-4, IL-10, IL-13, and TGF-β1, which are secreted by TAMs, are of significance in promoting tumor growth, EMT, or CSC formation [144,145,146]. TGF-β derived from breast cancer cells can stimulate the uPA gene, increase uPA mRNA stability, and activate uPA expression in TAMs [147]. TAMs can concentrate uPA at the cell surface by expressing urokinase receptor (uPAR), which can reshape the extracellular matrix through proteolysis, leading to the movement of breast cancer cells [148]. Breast cancer cells can downregulate the expression of epigenetic factor lysine demethylase 6B (KDM6B) in TAMs by delivering miR-138-5p via exosomes to TAMs [149]. Then, the decrease of KDM6B stimulates M2 polarization of the TAMs [149]. M2 macrophages can secrete chitinase 3-like protein 1 (CHI3L1) to interact with interleukin-13 receptor α2 (IL-13Rα2) chain on the plasma membranes of cancer cells to activate the mitogen-activated protein kinase (MAPK) signaling pathway, by which the expression of MMP genes are upregulated, leading to the promotion of the metastasis of breast cancer cells [150]. Breast cancer cells can downregulate the expression of protein phosphatase 2 catalytic subunit alpha (PPP2CA) in TAMs by delivering miR-183-5p via exosomes to TAMs [151]. PPP2CA downregulation enhances NF-κB signaling and promotes macrophage expression of IL-1β, IL-6, and TNF-α, which contributed to tumor progression [151].

In contrast, immune/inflammatory cells are recruited into tumor tissues during tumor growth to facilitate the overgrowth of cancer cells and suppress them. Activated T cells, B cells, natural killer (NK) cells, activated monocytes, dendritic cells, myeloid cells, and thymocyte subsets can express programmed cell death protein 1 (PD-1), an inhibitory immune checkpoint in breast cancer [152]. PD-L1, a PD-1 ligand, is expressed on cancer and immune cells and is associated with histological grade, hormone status, and prognosis of breast cancer [153, 154]. The expression of PD-L1 is regulated by diverse signaling in TME [155, 156]. Combinations of ICIs with chemotherapy have demonstrated efficacy in treating advanced/metastatic TNBC. They thus have been approved by the United States Food and Drug Administration (FDA) [157,158,159,160]. Progranulin (PGRN), a multifunctional growth factor expressed by breast cancer cells, can induce immune escape via upregulating PD-L1 expression on TAMs and promoting CD8+ T-cell exclusion [161]. The secretion of miR‐27a‐3p‐loaded exosomes from breast cancer cells, induced by endoplasmic reticulum stress, stimulates the PD‐L1 in macrophages and promotes immune evasion of breast cancer cells by activating the PTEN‐AKT/PI3K pathway [162]. The decreased expression of N-Myc downstream-regulated gene 2 (NDRG2) in breast cancer can upregulate PD-L1 expression of breast cancer cells and inhibit the proliferation of T cells [163]. Mouse UL16 binding protein-like transcript 1 (Mult1) is highly expressed in breast cancer cells carrying p53, which can bind to NK group 2 member D receptors (NKG2DRs) on NK cells, leading to NK activation and increased IFN-γ production, thus enabling the immune-eradication of cancer cells [164]. Previous studies mainly focused on the composition of immune cells in tumors. Still, more scholars believe that the spatial architecture of immune cells can affect tumor immunity and response to treatment, such as the spatial distribution of immune cells in tumors, the distance between immune cells, and other cells in the TME. [165].

Breast cancer cells—normal epithelial cells

Breast cancer originates from epithelial progenitor cells in the breast [166]. Therefore, breast epithelial cells are inextricably associated with breast cancer cells. Normal mammary epithelial cells are closest to breast cancer cells in the early stages of carcinogenesis. The interaction between cancer cells and epithelial cells is dynamic; normal epithelial cells initially inhibit tumor cells. Coculture verified that normal epithelial cells could induce transcriptional activation p53 of breast cancer cells, resulting in apoptosis of these cells in which tyrosine phosphatases may be involved [167]. Furthermore, normal epithelial cells could induce p53-independent cell cycle arrest of breast cancer cells [167].

But the tumor-suppressive effect may not last throughout solid tumor development. Breast cancer cells can transform normal breast epithelial cells to exhibit cancerous properties, including enhanced proliferation and migration, EMT-related features, and loss of apical-basal polarity [168]. S100A8/A9 can mediate dynamic interaction between normal breast epithelial cells and adjacent breast cancer cells to induce reprogramming [169]. Coculture with breast cancer cells or induction of overexpression of S100A8/A9 revealed that normal epithelial cells exhibit features of EMT, and the capacity of cell proliferation, migration, colony formation, and three-dimensional sphere formation is increased [169]. Angiostatin, amphiregulin, CD14, ETL, insulin-like growth factor (IGFBP)-2, IGFBP-6, IGFBP-7, IL-6, latent TGF-β BP1, osteoprotegerin, PDGF-AA, secreted protein acidic and rich in cysteine (SPARC), thrombospondin (TSP), TIMP1, and VEGF are also involved in this process [168]. Exosomes from cancer cells containing miRNAs play a significant role in intercellular communication [170]. Breast cancer cell-derived exosomes, miR-21 and -10b, can silence the expression of PTEN and homeobox D10 (HOXD10) in normal epithelial cells to induce tumor formation in a Dicer-dependent manner [171].

The interaction between epithelial and breast cancer cells changes over time, transforming from inhibition to promotion of tumor cells, affecting tumor growth, proliferation, invasion, and metastasis, thus leading to heterogeneity in breast cancer.

Breast cancer cells—breast cancer cells

Experiments have shown that the cooperation between breast cancer cells plays a vital role in the occurrence and development of breast cancer. This kind of cooperation can confer an advantage to breast cancer cell proliferation, invasion, and metastasis. Interestingly, neither luminal nor basal cells can generate cancers by themselves when transplanted into secondary recipient mice, whereas the combined transplantation of both cell types is highly carcinogenic [91]. The Wnt signaling pathway may play an essential role in this process. When this bi-clonal carcinoma is disturbed by Wnt blocking, basal subclones recruit allogeneic Wnt-producing cells to restore tumor growth [91]. EMT occurs in a subset of cells within primary breast cancers and has been shown to increase the growth of neighboring cancer cells [172]. For example, the DACH-Eyes absent (EYA)–Sineoculis homeobox homolog (SIX) network is involved in the initiation and progression of breast cancer, and their expression determines the prognosis of breast cancer patients [173,174,175,176]. When SIX1 is expressed in EMT cancer cells, VEGF-C is activated, and through neuropilin 2 (NRP2)/Fms-related tyrosine kinase 4 (FLT4), GLI signals in adjacent epithelial cancer cells are activated. The activation of GLI signals further promotes the growth, invasion, and metastasis of adjacent epithelial cancer cells, ultimately contributing to the progression of breast cancer [177]. Glutamine is a metabolite involved in the metabolic symbiosis between breast cancer cells. Glutamine-independent (luminal) breast cancer cells could synthesize glutamine to rescue glutamine-dependent (basal) BC cells in an in vitro coculture experiment, suggesting that the cell cooperation between basal and luminal cells in the mammary ducts may be achieved through glutamine symbiosis [178]. Hypoxic breast cancer cells can take up glucose and produce lactate, which can be taken up by aerobic breast cancer cells and catabolized to alanine and glutamate, which can then be exported to cells. Aerobic breast cancer cells have heterogeneity in the consumption capacity of lactate. Among them, those with high consumption of lactate express high monocarboxylate transporter (MCT)1/low MCT4 and reduce the use of glucose, which can leave more glucose to be used by other breast cancer cells [179]. This kind of cooperation can enable tumors to adapt to environmental changes and promote the generation of heterogeneity through the symbiosis of nutrients under different oxygen concentrations or nutrient states.

However, aerobic breast cancer cells with low consumption of lactate, with low MCT1/high MCT4 expression, consume less lactate and use relatively more glucose, which leads to the reduction of glucose available to hypoxic tumor cells [179]. Nutritional competition is a way in which cells compete. There are other aspects of cell competition, and many genes are thought to play essential roles in this process. MYC is one of the most commonly mutated genes in breast cancer [180], showing a strong correlation with cell competition [181]. MYC contributes to the generation of more competitive phenotypes; higher c-MYC breast cancer cells have higher numbers when cocultured with lower c-MYC breast cancer cells [182]. Higher expression of JNK, related to apoptosis, was observed in loser cells [182]. The alteration of human giant larvae homolog 1 (HUGL1) on breast cancer cell membrane can relieve the inhibition of Yes-Associated Protein (YAP), thus activating the transcription and expression of c-MYC and promoting the occurrence of Caspase-3 dependent apoptosis of adjacent cells with lower c-MYC expression, which leads to promotion of breast cancer cells with high expression of c-MYC to outcompete other cells in cell competition [183]. When breast cancer cells express high octamer-binding transcription factor 4 (Oct4), they also show inhibition to breast cancer cells with low Oct4 expression, similar to c-MYC [184]. This inhibition is even more potent when cells overexpress Oct4 and c-MYC simultaneously [184]. During this process, overexpression of Oct4 and c-MYC can induce higher expression of enolase 1 (Eno1), heat shock protein 90 α family class b member 1 (HSP90AB1), eukaryotic elongation factor 2 (Eef2), vinculin (VCL), Trail, and p53 to suppress other breast cancer cells [184]. During this process, the expression of PD-L1 and lysine demethylase 3A (Kdm3a) is downregulated and that of CD44 is upregulated in suppressed cells [184]. Breast cancer cells with higher MYC expression are more stem-like, which is regulated by MAPK/ERK. A higher proportion of MYC-positive cells in tumors significantly correlated with a short patient survival [185]. Thus, this c-MYC-mediated cellular competition helps eliminate the weaker cells to provide space and nutrition for the stronger ones, which may promote the emergence of more malignant phenotypes, heterogeneity, and progression of tumors, either naturally or during treatment. It was verified that the plasmacytoma variant translocation 1 (PVT1) promoter could inhibit MYC expression via promoter competition [186]. Furthermore, disruption of the PVT1 promoter by CRISPR can enhance breast cancer cell competition promoted by MYC [186]. The Flower is a downstream factor of MYC [96]. Breast cancer cells expressing Win human Flower isoforms (hFWE2 and hFWE4) have a competitive advantage over cells expressing Lose human isoforms (hFWE1 and hFWE3), contributing to the proliferation of winners and the cell death of losers [187]. Furthermore, breast cancer cells predominantly express the Win isoforms, whereas the stroma predominantly expresses the Lose isoforms, which may result in tumor growth and metastasis [187]. The activation of KRAS leads to the activation of Rac1 signaling, which can transmit the winner state of cells by downregulating contractile myosin, resulting in an increased mechanical deformability [188]. After an increase in mechanical deformability, cancer cells (winners) can engulf other cells (losers), and the engulfed cells experience cell death [188]. The competitive status of breast cancer cells is not invariable and will change or even reverse with the change of tumor living environment, such as pH and glucose concentration [100].

The interaction of breast cancer cells with TME was illustrated (Fig. 3). Indeed, many other mechanisms can explain the heterogeneity of breast cancer, such as the differentiation state of the cell of origin [98], cell plasticity, and CSC theory [189], genetic evolution of cancer, and multitask evolution [190, 191]. These mechanisms affect heterogeneity in different ways, each with advantages and shortcomings. Cell interaction provides a more comprehensive explanation of mechanisms, especially cell competition. Other mechanisms are involved in this process; for example, competition based on glycometabolism and lipid metabolism need to be explored intensively.

Fig. 3
figure 3

Interactions among diverse cells in the breast cancer microenvironment. The breast cancer microenvironment includes fibroblasts, adipocytes, immune cells, epithelial cells, and breast cancer cells. The interactions can be promoted (cooperative) or inhibited (competitive), leading to the generation of two cell types. The two types of cells undergo different fates, one more competitive, with more diverse and intense biological behavior, and the other gradually dying off

Research methods to measure heterogeneity

Given that breast cancer cells exhibit heterogeneity during the evolution of mechanisms of disease progression and therapeutic evasion routes, it is critical to monitor and quantify tumor heterogeneity during diagnosis, treatment, and possibly disease resistance. Heterogeneity is present at all tumor tissue levels, so its quantification must be measured at multiple scales.

Measuring heterogeneity in mixed cellular populations

Sequencing is one of the most used tools to measure heterogeneity [192, 193]. Since sequencing is an innovative technology in a constantly evolving field [9, 194], understanding of tumor heterogeneity is growing. Sanger sequencing is a typical example of first-generation sequencing at the genomic level. Sanger sequencing involves amplification of the DNA sequence, random termination by di-deoxy nucleotides, and fragments by their molecular weight [194]. At the transcriptomic level, RNA sequencing is often applied by the same sequencing platforms used to sequence DNA [195]. At the epigenomic level, nucleosome occupancy and methylome sequencing (NOMe–seq) based on nanopore sequencing, using an exogenous M. CviPI GpC methyltransferase, enables determination of phased patterns of native CpG methylation and chromatin accessibility [196]. As sequencing of the genome, epigenome, and transcriptome does not indicate fundamental biological processes, there was a pressing need for the sequencing of a whole protein.

At the proteomic level, in addition to traditional techniques, such as chromatography-based techniques, ELISA, and western blotting, many advanced methods have appeared, such as protein microarray, gel-based approaches, mass spectrometry (MS), Edman sequencing, and other high-throughput techniques [197, 198]. In the current development trend, the method of multi-omics combination will be a hot spot. The combination of third-generation DNA sequencing and proteomics is a prime example; the approximate shape, volume, charge, rotational diffusion coefficient, and dipole moment of individual proteins can be determined by zeptolitre sensing volume of bilayer-coated solid-state nanopores [199].

Measuring heterogeneity in single cells

Though it is insufficient, most of the above analyses are based on the analysis of tissue lysates, which hinders the cellular origin of the analyzed molecules, and some of the information of the smaller cell subsets is lost, affecting the diagnosis, treatment, and prognosis of tumors. A better understanding of genetic differences within cells is urgent and will lead to a better understanding of the heterogeneity of breast cancer. Based on second-generation sequencing, single-cell sequencing has progressed rapidly in recent years, revolutionizing our ability to drill down into thousands of individual cells' genomic, transcriptomic, epigenomic, and metabolic properties to enable an unbiased analysis of cells in tumoral lesions [200]. At the genomic level, the highly multiplexed single nucleus sequencing (HM-SNS) has been demonstrated as a good tool to explore the heterogeneity of copy number profiling in the breast cancer [201], which has been applied to explain the chemoresistance revolution in TNBC [16]. Deger et al. reported a modular pipeline constructed of five individually adjustable steps, including enrichment, isolation, whole-genome amplification, low pass sequencing, and bioinformatical analysis, to perform single CTCs genomics at a high resolution by copy number profiling [60]. Each step can be adjusted to use different methods according to specific research needs. Furthermore, it is a minimally invasive method to monitor molecular changes occurring during disease progression. At the epigenomic level, detecting multiallelic copy number events and copy-neutral loss-of-heterozygosity by single-cell transposase-accessible chromatin sequencing (scATAC-seq) enables the analysis of the contribution of chromosome instability and chromatin remodeling to tumor evolution [10]. Single-cell deoxyribonuclease I digestion sequencing (scDNase-seq) [11] and single-cell NOMe–seq (scNOMe-seq) [12] have emerged but are still rarely used at present. Due to its advantages in the field of epigenetics, it is believed to bring great convenience and breakthrough in the study of the heterogeneity of breast cancer shortly. Single-cell chromatin immunoprecipitation followed by sequencing approach (scChIP-seq), a high-throughput droplet microfluidics platform, paves the way to study the role of chromatin heterogeneity [202]. At the transcriptomic level, single-cell RNA sequencing (scRNA-seq) involves the capture of mRNA from compartmentalized cells, reverse transcription into cDNA, amplification, and identification by oligomeric indexes and molecular markers enabling high-throughput analysis of transcriptome [203]; scRNA-seq has been used widely in exploring heterogeneity of breast cancer, including metabolism [204], characterizing heterogeneous subgroups of different subtypes of breast cancer [205, 206], cell cycle [207], EMT [208], CTCs [209], and cell interaction [209]. At the proteomic level, the development of microfluidic imaging flow cytometry, which utilizes high-speed imaging technology, has facilitated the high spatial and temporal resolution of protein expression in individual cells. For example, an optofluidic flow cytometer (OFCM) combines a multistage microfluidic chip and a four-color fluorescence detection system. Single-cell phenotypic analysis can be done on 1.2 mL of whole blood per hour with OFCM and CTC separation, 3D focusing in the microchannel, and counting in real time [210]. Liquid chromatography - quadrupole time-of-flight tandem mass spectrometer (LC-Q-TOF-MS/MS) can achieve an extremely high degree of separation by LC; in addition, it provides information on the structure of flavonoids by Q-TOF-MS/MS [211], which has been applied to reveal the mechanism of treatment of TNBC [212].

Measuring heterogeneity in tissue slices preserving spatial information

Single-cell analysis can result in a loss of spatial information. To overcome this obstacle, the appearance of spatial omics has improved the imaging capabilities to produce subcellular spatial visualization and quantitative analysis of breast cancer tissue. At the genomic level, by combining laser-capture micro-dissection, laser catapulting, whole-genome amplification, and single-cell DNA sequencing, topographic single-cell sequencing (TSCS) can be used to measure genomic copy number profiles while preserving spatial information in the tissue slices. TSCS has revealed that one or more clones escape the ducts and migrate into the adjacent tissues to establish invasive carcinomas [213]. At the epigenomic level, the transposome-assisted single nucleus barcoding approach for ATAC-seq (SNuBar-ATAC) can easily label and multiplex a large number of samples together for parallel sequencing in a single microdroplet experiment with spatial information in breast cancer [214]. It has also been used in studying chromatin accessibility changes induced by drug treatment combinations. At the transcriptomic level, spatial transcriptomics provides quantitative gene expression data and visualization of mRNA distribution in tissue slices by locating tissue slices on aligned reverse transcription primers using unique location barcodes [215]. For miRNA analysis, fluorophore-encoded error-corrected labels (FluoELs) have low cytotoxicity and are capable of quantifying and spatially resolving breast cancer-related miRNAs and evaluating their coordination. For error-corrected quantification, FluoELs combine proportional dual fluorophores with a constant quantitative single fluorophore [216]. Cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq) can spatially characterize high-resolution immunophenotypes, including new PD-L1/PD-L2+ macrophage populations associated with clinical outcomes [19]. At the proteomic level, multiplexed ion beam imaging by the time-of-flight (MIBI-TOF) couples bright ion sources and orthogonal TOF-MS to image metal-tagged antibodies at the subcellular level [217] and has been used to explore the TME structure. MIBI-TOF verified that myoepithelial disruption could be a protective factor against recurrence [218]. Imaging mass cytometry, which couples immunohistochemical and immunocytochemical methods with high-resolution laser ablation to single-cell cytometry by TOF, reserving spatial information and delineating cell subpopulations and cell interactions in breast cancer [219]. High-resolution atlas of breast cancer has been well established, such as the atlas of CAFs [220], immune ecosystem [221, 222], and the interaction of these components with spatial information [223].

Measuring heterogeneity in patients through medical imaging

In addition to the level of tumor cells or tissues, the study of tumor heterogeneity on the patient level is important. Imaging techniques, including mammography, ultrasound (US), and magnetic resonance imaging (MRI) are conventional and effective tools for measuring the heterogeneity of breast cancer patients. Nevertheless, several limitations are associated with these techniques, such as high cost, harmful radiation, lack of sensibility, and inconvenience to the patients [224]. Positron-emission tomography (PET) and single-photon emission computed tomography (SPECT) are useful imaging techniques for the detection of metabolism and metastasis in breast cancer patients, but there is a lack of anatomical details in the imaging [225]. Thus, the combination of several techniques is recommended [226], such as PET/MRI and SPECT/MRI. However, these classical imaging methods, including mammography, US, MRI, PET, and SPECT, are primarily qualitative and subjective in tumor evaluation. Analyses based on tumor biopsies provide limited data on tumor characteristics because the sample extracted may not always represent the heterogeneity of the entire patient's tumor [227]. Due to the need for a quantitative and whole heterogeneity evaluation of breast cancer, radiomics has been gradually used in this field. Radiomics begins with obtaining high-quality images from mammography, US, MRI, PET, etc., and from these images, identifies an area of interest. They are segmented with operator edits and eventually rendered in three dimensions. In addition to clinical and genomic data, the quantitative features are extracted to generate reports for storage in the database. Then, these data are mined, which can be done using artificial intelligence, machine learning, or statistical approaches. Then, diagnostic, prediction, or prognostic models are established to obtain results [228]. The application of radiomics in breast cancer has been well reviewed [229, 230], and it will be possible to improve the diagnostic accuracy of breast imaging and increase the accuracy of heterogeneous studies by adding radiomics to the standard radiology workflow.

In conclusion, the quantification of heterogeneity of breast cancer needs to be measured at multiple scales. Each method has its advantages and disadvantages (Table 1). Depending on the actual situation, choosing the appropriate single or combined approach is essential.

Table 1 Advantages and disadvantages of various research methods to measure heterogeneity

Breast cancer heterogeneity and precise treatment

As mentioned above, considering heterogeneity, precise breast cancer treatments should be based on genetic, transcriptomic, epigenetic, proteomic, biomarkers, metabolism, cell cycle, TME, EMT, CTCs, and clinical and histopathological aspects. There are many standard treatment schemes that have been widely adopted in clinical practice based on breast cancer heterogeneity. Surgery, radiotherapy, and chemotherapy are standard therapies for all subtypes of breast cancer [232]. Standard breast surgery is either mastectomy, which means a complete excision of the breast, usually followed by reconstruction, or lumpectomy, which preserves the breast [233]. Axillary is the removal of the axillary lymph nodes to help determine the spread of cancer cells and remove the cancer cells [233]. Radiotherapy refers to high-energy radiation applied to the whole or part of the breast (after lumpectomy), the chest wall (after mastectomy), and regional lymph nodes [234]. A short course of radiotherapy of 3–4 weeks is usually as effective as a long course. Chemotherapy generally includes anthracyclines, taxanes, antimetabolites, alkylating agents, platinum drugs, and vinca alkaloids [235]. Chemotherapy can be further defined into NAC (before surgery) and adjuvant chemotherapy (after surgery). NAC, which often includes targeted agents, is a standard of care for non-metastatic but inoperable breast cancer and can help to reduce and adjust risk before or after breast and axillary surgery [232]. Different subtypes of tumors have varied sensitivities to neoadjuvant chemotherapy, possibly due to different immunological infiltrate biology. Increased TIL concentrations are positively associated with better response to neoadjuvant chemotherapy in HER2+ breast cancer and TNBC. However, increased TILs is negatively associated with better response to NAC in luminal-HER2 breast cancer [236]. Adjuvant chemotherapy is used in breast cancer patients with lymph node metastases or at high risk for recurrence [237]. As with neoadjuvant chemotherapy, ER breast cancer patients have lower recurrence and mortality rates than ER+ breast cancer patients after adjuvant chemotherapy [238]. Surgery, radiotherapy, and chemotherapy are common strategies for treating patients with BC. However, because patients respond differently to radiotherapy or chemotherapy, they are not effective enough to treat all BC molecular subtypes, so personalized treatment is essential.

Cyclin-dependent kinase (CDK) 4 and CDK6 inhibitors, PI3K inhibitors, polyadenosine-diphosphate-ribose polymerase (PARP) inhibitors, and anti-PD-L1 immunotherapy are examples of standard treatment options for metastatic breast cancer [232].

Personalized precision therapy for HER2+ breast cancer

As mentioned above, patients are usually managed based on the ER/PR/HER2 status of the cancer [233, 239]. As for early HER2+ breast cancer treatment, NAC includes trastuzumab plus pertuzumab plus chemotherapy (including taxanes, with or without anthracyclines or platinum drugs such as carboplatin) and is recommended before surgery [232]. Trastuzumab (Herceptin), a humanized monoclonal antibody targeting HER2, was the first biological drug-approved therapy and gold standard for HER2+ breast cancer [240]. Trastuzumab binds to subdomain IV of the HER2 extracellular domain [241]. Pertuzumab (Perjeta) binds to subdomain II of the HER2 extracellular domain to function [242]. After surgery, if patients have achieved pathological complete response (pCR), trastuzumab with or without pertuzumab is recommended for 12 months. If not, trastuzumab emtansine (T-DM1) is recommended for 14 cycles [232]. T-DM1 is administered because HER2 heterogeneity makes anti-HER2 therapy alone insufficient to eradicate HER2+ cancer completely. It is an antibody–drug conjugate composed of an anti-HER2 antibody and a microtubule inhibitor [243], showing high targeting ability and an excellent therapeutic effect on metastasis in several clinical trials [244, 245].

For metastatic/advanced HER2+ breast cancer, trastuzumab plus pertuzumab plus docetaxel or paclitaxel are recommended as the first-line therapy [232]. Second-line therapy consists of T-DMI. However, its role in eliminating HER2 heterogeneity is limited, with no pCR observed among patients with HER2 heterogeneity [246]. Trastuzumab deruxtecan (DS-8201) is another antibody–drug conjugate composed of an anti-HER2 antibody and a potent topoisomerase I inhibitor [247]. Its drug-antibody ratio is higher than T-DM1 while maintaining good pharmacokinetic characteristics [247]. DS-8201 has also been shown to play a therapeutic role in HER2+ metastatic breast cancer in several clinical trials [248, 249]. Furthermore, DS-8201 has been demonstrated to have a bystander-killing effect, which means deruxtecan is released into the spaces between cells and eliminates HER2 tumor cells after DS-8201 kills HER2-positive cancer cells to overcome intratumor heterogeneity, both in vitro and in vivo [250]. Trastuzumab duocarmazine (SYD985), which contains an anti-HER2 antibody, and duocarmycin [251], which can reduce the number of cancer cells through DNA damage and bystander-killing effects [251, 252], can be used to overcome intratumor heterogeneity. Its efficacy and safety have been demonstrated in a phase I study [253].

Personalized precision therapy for HR+ breast cancer

For the treatment of HR+/HER2 breast cancer, endocrine therapy, such as selective ER modulators (SERMs), selective ER deregulators (SERDs), and aromatase inhibitors (AIs), are mainstream. SERMs, such as tamoxifen, toremifene, bazedoxifene, and raloxifene, are antiestrogens that can competitively bind to the ER and render estrogen inactive [233]. SERDs, such as fulvestrant, once bound to the ER, inhibit receptor dimerization and prevent the ER's translocation to the nucleus, leading to its degradation and, in contrast to SERMs, completely block the ER signaling pathway [254]. AIs can be divided into two categories: steroidal AIs and non-steroidal AIs, which can block aromatase activity, inhibiting estrogen synthesis [255]. After surgery or radiotherapy, endocrine therapy is recommended for early-stage HR+ breast cancer. For postmenopausal women, tamoxifen is recommended, or AIs if contraindicated. For premenopausal women, ovarian function inhibitors are added [232, 233].

For metastatic/advanced HR+ breast cancer, a CDK4/6 inhibitor combined with endocrine therapy should be considered a standard of care [256]. For example, the combination of a CDK4/6 inhibitor named ribociclib and AIs, or tamoxifen, has better efficacy than endocrine therapy alone [257]. Another first-line treatment for advanced ER+ breast cancer is fulvestrant monotherapy [258]. Fulvestrant alone has no significant disadvantage compared with the combination of other endocrine therapies [233]. When treating patients with PIK3CA-mutant advanced HR+ tumors, a PIK3 inhibitor named alpelisib combined with fulvestrant is recommended, which can improve median overall survival by about eight months than fulvestrant alone [259]. PARP is a damage recognition repair protein of single-strand break, which plays a vital role in initiating the repair of the single-strand break. Inhibition of PARP leads to the accumulation of single-strand breaks and the formation of double-strand breaks.

BRCA-mutated cells cannot repair double-strand breaks, eventually leading to cell death [260]. When treating patients with BRCA1 or BRCA2-mutant advanced HR+ tumors, PARP inhibitors, such as olaparib or talazoparib, are recommended [261, 262]. Antibody-drug conjugates targeting the ER receptor are less well studied, among which Proteolysis-Targeting Chimeras (PROTACs) targeting ER are currently attracting attention in medicinal chemistry [263]. PROTACs are degrader-antibody conjugates that bind both an E3 ligase and a protein of interest and work by forming a ternary complex that initiates protein of interest ubiquitination and degradation by hijacking the ubiquitin–proteasome system [264]. ARV-471, a PROTAC targeting ERα, has entered the clinical phase II trial [265], has prominent antitumor activity, and reduces ER protein expression by more than 90% in cell lines and 62% in tumor tissues. Furthermore, several new PROTACs are undergoing experimental assessment [263]. The studies of PROTACs will make a breakthrough contribution to overcoming drug resistance in ER+ breast cancer in the future. As for HR+/HER2+ patients, a combination of hormone therapy and targeted anti-HER2 approaches, which have been shown to control the disease over the long term, is preferable [32].

Personalized precision therapy for triple-negative breast cancer

The conventional treatment mode of TNBC, known for its challenges in therapy with intratumor heterogeneity, involves a combination of surgery, radiotherapy, and chemotherapy. Despite the emergence of new biological and targeted agents, the primary treatment for TNBC remains cytotoxic chemotherapy. As for the treatment of early TNBC breast cancer, chemotherapy (including adriamycin, cyclophosphamide, and paclitaxel, with or without carbo platinum) is recommended before surgery [232, 266]. After surgery, if patients have achieved pCR, they need only be followed up regularly. If not, these patients should consider capecitabine, an oral prodrug of fluorouracil, for 6–8 courses [232]. There are clinical trials that show that adjuvant chemotherapy with capecitabine after standard neoadjuvant chemotherapy containing anthracyclines, taxanes, or both is safe and effective in prolonging disease-free survival and overall survival in patients with TNBC who have residual invasive disease on pathological testing [267].

For metastatic/advanced TNBC breast cancer, patients should complete testing for a germline BRCA mutation and PD-L1 before treatment [232]. Because of remarkable genomic instability and increased immune infiltration, certain TNBCs display high PD-L1 expression compared with other subtypes [268]. These properties make patients with TNBC good candidates for ICIs such as atezolizumab and pembrolizumab. Therapeutic blockade of PD-L1 with atezolizumab [269] or therapeutic blockade of PD-1 with pembrolizumab [270] can activate and enhance tumor-specific T-cell responses, resulting in increased antitumor activity. In a phase I study of 111 patients with TNBC, the PD-1 inhibitor pembrolizumab [157] demonstrated therapeutic effects and acceptable safety. If patients are PD-L1+, nab-paclitaxel plus atezolizumab or paclitaxel plus pembrolizumab are recommended [232]. If patients are BRCA mutation-positive, PARPi are recommended [232] as examples of PARPi, olaparib and veliparib are being evaluated in clinical experiments [271, 272]. If patients are PD-L1 or BRCA wild-type, combination chemotherapy or sequential single-agent chemotherapy should be considered, including anthracycline or taxanes. If these drugs are not available, carboplatin, eribulin, vinorelbine, capecitabine, and sacituzumab govitecan are recommended [232]. Sacituzumab govitecan is a type of antibody-drug conjugate that combines a topoisomerase I inhibitor SN-38 and an antibody targeting trophoblast antigen 2 (Trop-2) [273]. In a randomized, phase III trial, among patients with metastatic TNBC, progression-free survival and overall survival were about twice as long in those treated with sacituzumab govitecan as in those treated with single-agent chemotherapy [274]. Vaccination is an emerging approach to stopping the recurrence of TNBC in high-risk patients. The individual selection of vaccine antigens from a set of different candidate peptides based on the pre-existing host immunity is known as personalized peptide vaccination [275]. In a phase II study of personalized peptide vaccination for metastatic recurrent breast cancer patients, personalized peptide vaccination was reported to have limited adverse events and boost cytotoxic T lymphocyte or IgG response in patients, especially TNBC patients [275].

The precise treatment of each subtype of TNBC is still in its infancy. LAR is characterized by high expression of androgen receptors, which can be a target for treatment. Bicalutamide, an androgen receptor (AR) agonist, was evaluated in a multicenter phase II trial. The results showed that the 6-month clinical benefit rate was 19% for bicalutamide, and the median progression-free survival was 12 weeks with no grade 4/5 treatment-related adverse events observed [276]. Enzalutamide and abiraterone are androgen receptor inhibitors, showing significant clinical activity and tolerability in patients with locally advanced or metastatic AR+ TNBC [277, 278]. Furthermore, 70% of LAR tumors exhibit somatic mutations in the PI3K signaling pathway; thus, PI3K and AKT inhibitors may be beneficial [279]. LAR retained retinoblastoma tumor suppressor and showed frequent CDK inhibitor 2A gene alterations, both of which are associated with CDK4/6 inhibitors. Therefore, patients with LAR tumors may benefit from CDK4/6 inhibitors [279]. IM subtype is characterized by elevated immune cell signaling and TILs. Additionally, ICI genes such as PD1 and PD-L1 showed high expression in IM [279]. Thus, ICIs may be a potential strategy for IM. BLIS subtype is characterized by high genomic instability, suggesting that PARP inhibitors and other DNA-damaging chemotherapy can be a potential strategy for BLIS [279]. The MES subtype is characterized by overexpression of stem cell-related genes such as the STAT3 signaling pathway, suggesting that targeting the STAT3 pathway can be a potential strategy for MES [279].

Personalized precision therapy for claudin-low breast cancer

Research on claudin-low breast cancer has gradually increased in popularity recently. The high rate of metastasis and mortality in patients with claudin-low tumors remains troubling [280] because of its aggression, chemoresistance, and lack of targeted therapies. Claudin-low tumors are characterized by high enrichment for EMT markers, immune response genes, and cancer stem cell-like features; these are potential therapeutic targets. Topsentinol L trisulfate can inhibit the activation of AMP-activated protein kinase and checkpoint kinase 1 and promote activation of p38 in claudin-low breast cancer. Topsentinol L trisulfate has higher efficacy against claudin-low breast cancer than other types [281]. Most claudin-low breast tumors are reported to be ER, PR, and HER2; therefore, the treatment of TNBC can play a guiding role in the treatment of claudin-low breast cancer to a certain extent [30, 282].

Personalized precision therapy targeting cell interactions

Targeting cell interactions as a therapeutic breakthrough is a promising approach. In theory, all the genes or proteins mentioned above that interact with cells have potential therapeutic implications. Currently, related drugs are gradually being used in clinical practice. AVI-4126 is a phosphorodiamidate morpholino oligomer that can inhibit MYC expression by preventing ribosomal assembly [283]. In a phase I study, the findings suggest that AVI-4126 is concentrated in breast cancer as expected and can inhibit the target [283]. TAS-119 has the potential to target the MYC and Wnt/β-catenin pathways [284], with no complete response or partial response after the use of TAS-119 in patients with MYC-amplified/β-catenin mutations [285]. Furthermore, short-term doxorubicin and cisplatin administration may induce a more favorable TME and increase the likelihood of response to PD-1 blockade in the TNBC [158]. Bevacizumab, a type of VEGF inhibitor [286, 287], and cediranib, a VEGFR1-3 inhibitor [288], have also achieved good results in breast cancer treatment. Bispecific antibodies (BsAbs) are another candidate for the effective treatment of cancer and are currently in various stages of clinical trials [289, 290]. A BsAb named M802 consists of a monovalent unit against HER2 and a single chain unit against CD3. M802 could recruit CD3+ immune cells to eliminate breast cancer cells. It was more cytotoxic than trastuzumab in cells with high expression of HER2, low expression of HER2, and trastuzumab resistance [291]. Ertumaxomab is a bispecific, trifunctional antibody that binds to tumor cells expressing HER2/neu, T cells expressing CD3, and accessory cells expressing Fcγ, forming a tri-cell complex [292]. In a phase I trial, Ertumaxomab showed promising efficacy and transient and reversible side effects of HER2 overexpressing cancers, including breast cancer [292]. Anti-CD3 × anti-HER2 bispecific antibody-armed activated T cells (HER2 BATs) showed promising efficacy in killing breast cancer cells expressing high levels of HER2 and low HER2 receptor expressing breast cancer cells in a non-MHC restricted manner by perforin and granzyme B [293, 294]. A phase I trial suggested that HER2 BATs can target HER2+ and HER2 tumors and increase Th1 cytokines and IL12 against breast cancer [295]. In a phase II trial, immune consolidation (IL-2 and granulocyte–macrophage-colony stimulating factor) with HER2 BATs after chemotherapy was demonstrated to increase the proportion of patients who were stable at four months and the median overall survival in HER2 metastatic breast cancer patients [293]. Unfortunately, immunotherapy for ER+ breast cancer has historically fared poorly, thus making ER+ breast cancer the most difficult subtype for BsAbs to target [296]. This will be one of the development directions of BsAbs in the future.

Personalized precision therapy guided by single-cell sequencing

Despite significant advances in heterogeneous individualized therapy for breast cancer, with a growing list of targeted drugs in adjuvant, neoadjuvant, and metastatic therapies, the need to finetune precision drug therapy to improve survival and eliminate late relapse remains unmet. Single-cell sequencing is a powerful tool for exploring tumor heterogeneity and guiding the improvement of traditional therapies. As mentioned above, single-cell sequencing can detect cell resistance to conventional therapies to apply further novel drugs targeting these cells [73]. Single-cell RNA sequencing can be a powerful tool for discovering novel therapies associated with the immune checkpoint crosstalk [297]. Radioresistant tumor cells have a higher rate of PD-L1 positivity and tumor mutation burden, as indicated by the single-cell RNA sequencing [297]. Thus, combining immune checkpoint therapy and radiation therapy guided by single-cell RNA sequencing may be beneficial. Single-cell sequencing can also monitor changes in the genome, epigenome, transcriptome, and proteome after treatment and adjust the treatment regimen. Single-cell RNA sequencing revealed that intratumoral plasmid IL12 can expand CD8+ T cells, induce a CXCR3 gene signature, and sensitize patients to anti-PD-L1 therapy [298]. A patient who had not previously responded to anti-PD-L1 therapy had increased CXCR3 after intratumoral injection of plasmid IL12. Immediately, she received additional anti-PD-1 therapy and showed a significant clinical response. Minimally invasive thermal therapy has been attempted for breast cancer and can induce an immune response. By performing single-cell RNA sequencing on peripheral blood mononuclear cells from patients before and after microwave ablation, researchers found that B cells are critical antigen-presenting cells that initiate CD4+ T cells in the microwave ablation-induced immune response [299]. This study provides a comprehensive picture of microwave ablation-induced systemic immune responses and opens a new way to identify and improve potential targets of immune responses. Combining single-cell sequencing with other technologies further promotes the development of personalized therapies, such as hydrosequencing, which can analyze ER, PR, HER2 expression, and other clinical markers of CTCs from 10 mL of blood [300].

Conclusion

We summarized the current state of knowledge on the origins of breast cancer heterogeneity and highlighted cell interactions in this process. We have reviewed different levels of breast cancer heterogeneity and described its enormous complexity. In particular, we have defined new research methods, a possible mechanism of cell interaction, and clinical treatment progress. Cell interactions have a unique advantage in explaining the genesis of heterogeneity and further treatment of breast cancer. This theory treats cells as individuals and endows them with characteristics and behaviors such as cooperation and competition between cells, just like individual organisms in nature, which is appropriate and reasonable. Targeting the heterogeneity of breast cancer at different levels, intercellular interactions, supported by single-cell sequencing and spatial transcriptome, will be the focus of future research and treatment.

Availability of data and materials

Not applicable.

Change history

Abbreviations

TME:

Tumor microenvironment

EMT:

Epithelial–mesenchymal transition

CTCs:

Circulating tumor cells

ER:

Estrogen receptor

PR:

Progesterone receptor

HRs:

Hormone receptors 

HER2:

Human epidermal growth factor receptor 2

TNBC:

Triple-negative breast cancer

LAR:

Luminal androgen receptor

MES:

Mesenchymal-like

BLIS:

Basal-like and immune suppressed

CSCs:

Cancer stem cells

BCSCs:

Breast cancer stem cells

TNM:

Tumor, lymph node, metastasis

VM:

Vasogenic mimicry

NAC:

Neoadjuvant chemotherapy

SERMs:

Selective ER modulators

SERDs:

Selective ER deregulators

AIs:

Aromatase inhibitors

TILs:

Tumor-infiltrating lymphocytes

NES:

Neutrophil-enriched subtype

MES:

Macrophage-enriched subtype

CCR2:

C–C chemokine receptor type 2

ICB:

Checkpoint blockade

CAFs:

Cancer-associated fibroblasts

SDF-1:

Stromal-derived factors

CXCL12:

C-X-C motif ligand 12

TGFβ:

Transforming growth factor-beta

IL-6:

Interleukin-6

IL-8:

Interleukin-8

HIF-α:

Hypoxia-inducible factor-α

S100A4:

S100 calcium-binding protein A4

FGF2:

Fibroblast growth factor 2

FGF7:

Fibroblast growth factor 7

MMP2:

Matrix metalloproteinase 2

MMP9:

Matrix metalloproteinase 9

MMP11:

Matrix metalloproteinase 11

PDGFA:

Platelet-derived growth factor A

PDGFB:

Platelet-derived growth factor B

VEGFA:

Vascular endothelial growth factor A

uPA:

Urokinase-type plasminogen activator

TIMP1:

Tissue inhibitor matrix metalloproteinase 1

BMP:

Bone morphogenetic protein

GREM1:

GREMLIN1

miR:

microRNA

PTEN:

Phosphatase and tensin homolog

ITGB4:

Integrin β 4

CAAs:

Cancer-associated adipocytes

JAK:

Janus kinase

STAT3:

Signal transducer and activator of transcription 3

PI3K:

Phosphoinositide 3-kinases

AKT:

Protein kinase B

PLOD2:

Lysyl hydroxylase

IL-1β:

Interleukin-1 beta

TNF-α:

Tumor necrosis factor  α

PAI-1:

Plasminogen activator inhibitor-1

FFAs:

Free fatty acids

MCP-1:

Monocyte chemoattractant protein-1

CCL5:

CC-chemokine ligand 5

IGF-1:

Insulin-like growth factor 1

MVP:

Major vault protein

TAMs:

Tumor-associated macrophages

IL-4:

Interleukin-4

IL-10:

Interleukin-10

IL-13:

Interleukin-13

uPAR:

Urokinase receptor

KDM6B:

Lysine demethylase 6B

CHI3L1:

Chitinase 3-like protein 1

IL-13Rα2:

Interleukin-13 receptor α2 chain

MAPK:

Mitogen-activated protein kinase

PPP2CA:

Protein phosphatase 2 catalytic subunit alpha

NK:

Natural killer

PD-1:

Programmed death  protein 1

ICIs:

Immune checkpoint inhibitors

FDA:

Food and Drug Administration

PGRN:

Progranulin

NDRG2:

N-myc downstream-regulated gene 2

Mult1:

Mouse UL16 binding protein-like transcript 1

NKG2DRs:

Natural killer group 2 member D receptors

S100A8:

S100 calcium-binding protein A8

S100A9:

S100 calcium-binding protein A9

IGFBP:

Insulin-like growth factor

SPARC:

Secreted protein acidic and rich in cysteine

TSP:

Thrombospondin

HOXD10:

Homeobox D10

DACH:

Diaminocyclohexane

EYA:

Eyes absent

SIX1:

Sineoculis homeobox homolog 1

NRP2:

Neuropilin 2

FLT4:

Fms-related tyrosine kinase 4

MCT:

Monocarboxylate transporter

HUGL1:

Human giant larvae homolog 1

YAP:

Yes-associated protein 1

Oct-4:

Octamer-binding transcription factor 4

Eno1:

Enolase 1

HSP90AB1:

Heat Shock Protein 90 Alpha Family Class B Member 1

Eef2:

Eukaryotic elongation factor 2

VCL:

Vinculin

Kdm3a:

Lysine demethylase 3A

PVT1:

Plasmacytoma variant translocation 1

hFWE:

Human Flower isoform

NOMe–seq:

Nucleosome occupancy and methylome sequencing

MS:

Mass spectrometry

HM-SNS:

highly multiplexed single nucleus sequencing

scChIP-seq:

Single-cell chromatin immunoprecipitation followed by sequencing

scATAC-seq:

Single-cell transposase-accessible chromatin-sequencing

scDNase-seq:

Single-cell DNase sequencing

scNOMe–seq:

Single-cell NOMe–seq

scRNA-seq:

Single-cell RNA sequencing

OFCM:

Optofluidic flow cytometer

LC-Q-TOF-MS/MS:

Liquid chromatography - quadrupole time-of-flight tandem mass spectrometer

TSCS:

Topographic single cell sequencing

SNuBar-ATAC:

Single nucleus barcoding approach for transposase-accessible chromatin-sequencing

FluoELs:

Fluorophores encoded error-corrected labels

CITE-seq:

Cellular indexing of transcriptomes and epitopes by sequencing

MIBI-TOF:

Multiplexed ion beam imaging by time of flight

IMC:

Imaging mass cytometry

US:

Ultrasound

MRI:

Magnetic resonance imaging

PET:

Positron-emission tomography

SPECT:

Single-photon emission computed tomography

T-DM1:

Trastuzumab emtansine

pCR:

Pathologic complete response

CDK:

Cyclin-dependent kinase

PARP:

Polyadenosine-diphosphate-ribose polymerase

DS-8201:

Trastuzumab deruxtecan

PROTACs:

Proteolysis-targeting chimera

SYD985:

Trastuzumab duocarmazine

Trop-2:

Trophoblast antigen 2

AR:

androgen receptor

BsAb:

Bispecific antibody

References

  1. Kimbung S, Johansson I, Danielsson A, Veerla S, Egyhazi Brage S, Frostvik Stolt M, et al. Transcriptional profiling of breast cancer metastases identifies liver metastasis-selective genes associated with adverse outcome in luminal a primary breast cancer. Clin Cancer Res. 2016;22:146–57.

    Article  CAS  PubMed  Google Scholar 

  2. Siegel RL, Miller KD, Fuchs HE, Jemal A. Cancer statistics, 2021. CA Cancer J Clin. 2021;71:7–33.

    Article  PubMed  Google Scholar 

  3. Yi M, Li T, Niu M, Luo S, Chu Q, Wu K. Epidemiological trends of women’s cancers from 1990 to 2019 at the global, regional, and national levels: a population-based study. Biomark Res. 2021;9:55.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Visvader JE. Cells of origin in cancer. Nature. 2011;469:314–22.

    Article  CAS  PubMed  Google Scholar 

  5. Dagogo-Jack I, Shaw AT. Tumour heterogeneity and resistance to cancer therapies. Nat Rev Clin Oncol. 2018;15:81–94.

    Article  CAS  PubMed  Google Scholar 

  6. Saleem MZ, Alshwmi M, Zhang H, Din SRU, Nisar MA, Khan M, et al. Inhibition of JNK-mediated autophagy promotes proscillaridin A- induced apoptosis via ROS generation, intracellular Ca(+2) oscillation and inhibiting STAT3 signaling in breast cancer cells. Front Pharmacol. 2020;11:01055.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Torre LA, Sauer AM, Chen MS Jr, Kagawa-Singer M, Jemal A, Siegel RL. Cancer statistics for Asian Americans, Native Hawaiians, and Pacific Islanders, 2016: Converging incidence in males and females. CA Cancer J Clin. 2016;66:182–202.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Chen W, Zheng R, Baade PD, Zhang S, Zeng H, Bray F, et al. Cancer statistics in China, 2015. CA Cancer J Clin. 2016;66:115–32.

    Article  PubMed  Google Scholar 

  9. Kashyap A, Rapsomaniki MA, Barros V, Fomitcheva-Khartchenko A, Martinelli AL, Rodriguez AF, et al. Quantification of tumor heterogeneity: from data acquisition to metric generation. Trends Biotechnol. 2022;40:647–76.

    Article  CAS  PubMed  Google Scholar 

  10. Wu CY, Lau BT, Kim HS, Sathe A, Grimes SM, Ji HP, et al. Integrative single-cell analysis of allele-specific copy number alterations and chromatin accessibility in cancer. Nat Biotechnol. 2021;39:1259–69.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Gao W, Ku WL, Pan L, Perrie J, Zhao T, Hu G, et al. Multiplex indexing approach for the detection of DNase I hypersensitive sites in single cells. Nucleic Acids Res. 2021;49: e56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Pott S. Simultaneous measurement of chromatin accessibility, DNA methylation, and nucleosome phasing in single cells. Elife. 2017;6: e23203.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Koboldt DC, Fulton RS, McLellan MD, Schmidt H, Kalicki-Veizer J, McMichael JF, et al. Comprehensive molecular portraits of human breast tumours. Nature. 2012;490:61–70.

    Article  CAS  Google Scholar 

  14. Curtis C, Shah SP, Chin SF, Turashvili G, Rueda OM, Dunning MJ, et al. The genomic and transcriptomic architecture of 2,000 breast tumours reveals novel subgroups. Nature. 2012;486:346–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Eirew P, Steif A, Khattra J, Ha G, Yap D, Farahani H, et al. Dynamics of genomic clones in breast cancer patient xenografts at single-cell resolution. Nature. 2015;518:422–6.

    Article  CAS  PubMed  Google Scholar 

  16. Kim C, Gao R, Sei E, Brandt R, Hartman J, Hatschek T, et al. Chemoresistance evolution in triple-negative breast cancer delineated by single-cell sequencing. Cell. 2018;173:879-93.e13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Wang Y, Waters J, Leung ML, Unruh A, Roh W, Shi XQ, et al. Clonal evolution in breast cancer revealed by single nucleus genome sequencing. Nature. 2014;512:155–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Navin N, Kendall J, Troge J, Andrews P, Rodgers L, McIndoo J, et al. Tumour evolution inferred by single-cell sequencing. Nature. 2011;472:90-U119.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Wu SZ, Al-Eryani G, Roden DL, Junankar S, Harvey K, Andersson A, et al. A single-cell and spatially resolved atlas of human breast cancers. Nat Genet. 2021;53:1334–47.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Hammond MEH, Hayes DF, Dowsett M, Allred DC, Hagerty KL, Badve S, et al. American Society of Clinical Oncology/College of American Pathologists Guideline recommendations for immunohistochemical testing of estrogen and progesterone receptors in breast cancer. J Clin Oncol. 2010;28:2784–95.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Wolff AC, Hammond MEH, Hicks DG, Dowsett M, McShane LM, Allison KH, et al. Recommendations for human epidermal growth factor receptor 2 testing in breast cancer: American Society of Clinical Oncology/College of American Pathologists clinical practice guideline update. J Clin Oncol. 2013;31:3997–4013.

    Article  PubMed  Google Scholar 

  22. Tang ZR, Zhang R, Lian ZX, Deng SL, Yu K. Estrogen-receptor expression and function in female reproductive disease. Cells. 2019;8:1123.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Graham JD, Yeates C, Balleine RL, Harvey SS, Milliken JS, Bilous AM, et al. Characterization of progesterone receptor A and B expression in human breast cancer. Cancer Res. 1995;55:5063–8.

    CAS  PubMed  Google Scholar 

  24. Petrelli F, Tomasello G, Barni S, Lonati V, Passalacqua R, Ghidini M. Clinical and pathological characterization of HER2 mutations in human breast cancer: a systematic review of the literature. Breast Cancer Res Treat. 2017;166:339–49.

    Article  CAS  PubMed  Google Scholar 

  25. Dietze EC, Sistrunk C, Miranda-Carboni G, O’Regan R, Seewaldt VL. Triple-negative breast cancer in African–American women: disparities versus biology. Nat Rev Cancer. 2015;15:248–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Gao JJ, Swain SM. Luminal A breast cancer and molecular assays: a review. Oncologist. 2018;23:556–65.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Ades F, Zardavas D, Bozovic-Spasojevic I, Pugliano L, Fumagalli D, de Azambuja E, et al. Luminal B breast cancer: molecular characterization, clinical management, and future perspectives. J Clin Oncol. 2014;32:2794–803.

    Article  PubMed  Google Scholar 

  28. Loibl S, Gianni L. HER2-positive breast cancer. Lancet. 2017;389:2415–29.

    Article  CAS  PubMed  Google Scholar 

  29. Liu YR, Jiang YZ, Xu XE, Yu KD, Jin X, Hu X, et al. Comprehensive transcriptome analysis identifies novel molecular subtypes and subtype-specific RNAs of triple-negative breast cancer. Breast Cancer Res. 2016;18:33.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Prat A, Parker JS, Karginova O, Fan C, Livasy C, Herschkowitz JI, et al. Phenotypic and molecular characterization of the claudin-low intrinsic subtype of breast cancer. Breast Cancer Res. 2010;12:18.

    Article  Google Scholar 

  31. Fougner C, Bergholtz H, Norum JH, Sørlie T. Re-definition of claudin-low as a breast cancer phenotype. Nat Commun. 2020;11:1787.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Kay C, Martínez-Pérez C, Meehan J, Gray M, Webber V, Dixon JM, et al. Current trends in the treatment of HR+/HER2+ breast cancer. Future Oncol. 2021;17:1665–81.

    Article  CAS  PubMed  Google Scholar 

  33. Schrijver W, Suijkerbuijk KPM, van Gils CH, van der Wall E, Moelans CB, van Diest PJ. Receptor conversion in distant breast cancer metastases: a systematic review and meta-analysis. J Natl Cancer Inst. 2018;110:568–80.

    Article  PubMed  Google Scholar 

  34. Turner NH, Di Leo A. HER2 discordance between primary and metastatic breast cancer: assessing the clinical impact. Cancer Treat Rev. 2013;39:947–57.

    Article  CAS  PubMed  Google Scholar 

  35. Kroigard AB, Larsen MJ, Thomassen M, Kruse TA. Molecular concordance between primary breast cancer and matched metastases. Breast J. 2016;22:420–30.

    Article  PubMed  Google Scholar 

  36. Cui X, Schiff R, Arpino G, Osborne CK, Lee AV. Biology of progesterone receptor loss in breast cancer and its implications for endocrine therapy. J Clin Oncol. 2005;23:7721–35.

    Article  CAS  PubMed  Google Scholar 

  37. Aurilio G, Disalvatore D, Pruneri G, Bagnardi V, Viale G, Curigliano G, et al. A meta-analysis of oestrogen receptor, progesterone receptor and human epidermal growth factor receptor 2 discordance between primary breast cancer and metastases. Eur J Cancer. 2014;50:277–89.

    Article  CAS  PubMed  Google Scholar 

  38. Thomson AH, McGrane J, Mathew J, Palmer J, Hilton DA, Purvis G, et al. Changing molecular profile of brain metastases compared with matched breast primary cancers and impact on clinical outcomes. Br J Cancer. 2016;114:793–800.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Szekely B, Nagy ZI, Farago Z, Kiss O, Lotz G, Kovacs KA, et al. Comparison of immunophenotypes of primary breast carcinomas and multiple corresponding distant metastases: an autopsy study of 25 patients. Clin Exp Metastasis. 2017;34:103–13.

    Article  CAS  PubMed  Google Scholar 

  40. Yao ZX, Lu LJ, Wang RJ, Jin LB, Liu SC, Li HY, et al. Discordance and clinical significance of ER, PR, and HER2 status between primary breast cancer and synchronous axillary lymph node metastasis. Med Oncol. 2014;31:798.

    Article  PubMed  Google Scholar 

  41. Yeung C, Hilton J, Clemons M, Mazzarello S, Hutton B, Haggar F, et al. Estrogen, progesterone, and HER2/neu receptor discordance between primary and metastatic breast tumours-a review. Cancer Metastasis Rev. 2016;35:427–37.

    Article  CAS  PubMed  Google Scholar 

  42. Chen R, Qarmali M, Siegal GP, Wei S. Receptor conversion in metastatic breast cancer: analysis of 390 cases from a single institution. Mod Pathol. 2020;33:2499–506.

    Article  CAS  PubMed  Google Scholar 

  43. Dias AS, Almeida CR, Helguero LA, Duarte IF. Metabolic crosstalk in the breast cancer microenvironment. Eur J Cancer. 2019;121:154–71.

    Article  CAS  PubMed  Google Scholar 

  44. Willmann L, Schlimpert M, Halbach S, Erbes T, Stickeler E, Kammerer B. Metabolic profiling of breast cancer: differences in central metabolism between subtypes of breast cancer cell lines. J Chromatogr B Analyt Technol Biomed Life Sci. 2015;1000:95–104.

    Article  CAS  PubMed  Google Scholar 

  45. Willmann L, Schlimpert M, Hirschfeld M, Erbes T, Neubauer H, Stickeler E, et al. Alterations of the exo- and endometabolite profiles in breast cancer cell lines: a mass spectrometry-based metabolomics approach. Anal Chim Acta. 2016;925:34–42.

    Article  CAS  PubMed  Google Scholar 

  46. Davis RT, Blake K, Ma D, Gabra MBI, Hernandez GA, Phung AT, et al. Transcriptional diversity and bioenergetic shift in human breast cancer metastasis revealed by single-cell RNA sequencing. Nat Cell Biol. 2020;22:310–20.

    Article  CAS  PubMed  Google Scholar 

  47. Mason JM, Wei X, Fletcher GC, Kiarash R, Brokx R, Hodgson R, et al. Functional characterization of CFI-402257, a potent and selective Mps1/TTK kinase inhibitor, for the treatment of cancer. Proc Natl Acad Sci U S A. 2017;114:3127–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Yano S, Hoffman RM. Real-time determination of the cell-cycle position of individual cells within live tumors using FUCCI cell-cycle imaging. Cells. 2018;7:168.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Lashen AG, Toss MS, Katayama A, Gogna R, Mongan NP, Rakha EA. Assessment of proliferation in breast cancer: cell cycle or mitosis? An observational study. Histopathology. 2021;79:1087–98.

    Article  PubMed  Google Scholar 

  50. Wang M, Zhao J, Zhang L, Wei F, Lian Y, Wu Y, et al. Role of tumor microenvironment in tumorigenesis. J Cancer. 2017;8:761–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Ma B, Cheng H, Mu C, Geng G, Zhao T, Luo Q, et al. The SIAH2-NRF1 axis spatially regulates tumor microenvironment remodeling for tumor progression. Nat Commun. 2019;10:1034.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Gupta GP, Massagué J. Cancer metastasis: building a framework. Cell. 2006;127:679–95.

    Article  CAS  PubMed  Google Scholar 

  53. Davis FM, Azimi I, Faville RA, Peters AA, Jalink K, Putney JW Jr, et al. Induction of epithelial–mesenchymal transition (EMT) in breast cancer cells is calcium signal dependent. Oncogene. 2014;33:2307–16.

    Article  CAS  PubMed  Google Scholar 

  54. Imaoka H, Ikeda M, Maehara K, Umemoto K, Ozaka M, Kobayashi S, et al. Clinical outcomes of chemotherapy in patients with undifferentiated carcinoma of the pancreas: a retrospective multicenter cohort study. BMC Cancer. 2020;20:946.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Tiran V, Lindenmann J, Brcic L, Heitzer E, Stanzer S, Tabrizi-Wizsy NG, et al. Primary patient-derived lung adenocarcinoma cell culture challenges the association of cancer stem cells with epithelial-to-mesenchymal transition. Sci Rep. 2017;7:10040.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Dong B, Li S, Zhu S, Yi M, Luo S, Wu K. MiRNA-mediated EMT and CSCs in cancer chemoresistance. Exp Hematol Oncol. 2021;10:12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Liu X, Li J, Cadilha BL, Markota A, Voigt C, Huang Z, et al. Epithelial-type systemic breast carcinoma cells with a restricted mesenchymal transition are a major source of metastasis. Sci Adv. 2019;5:eaav4275.

    Article  PubMed  PubMed Central  Google Scholar 

  58. Micalizzi DS, Maheswaran S, Haber DA. A conduit to metastasis: circulating tumor cell biology. Genes Dev. 2017;31:1827–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Marrinucci D, Bethel K, Kolatkar A, Luttgen MS, Malchiodi M, Baehring F, et al. Fluid biopsy in patients with metastatic prostate, pancreatic and breast cancers. Phys Biol. 2012;9: 016003.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Deger T, Mendelaar PAJ, Kraan J, van der Smissen WJCP, van der Vlugt-Daane M, Bindels EMJ, et al. A pipeline for copy number profiling of single circulating tumor cells to assess intra-patient tumor heterogeneity. Mol Oncol. 2021. https://doi.org/10.1002/1878-0261.13174.

    Article  Google Scholar 

  61. Folberg R, Maniotis AJ. Vasculogenic mimicry. APMIS. 2004;112:508–25.

    Article  PubMed  Google Scholar 

  62. Morales-Guadarrama G, García-Becerra R, Méndez-Pérez EA, García-Quiroz J, Avila E, Díaz L. Vasculogenic mimicry in breast cancer: clinical relevance and drivers. Cells. 2021;10:1758.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Luo Q, Wang J, Zhao W, Peng Z, Liu X, Li B, et al. Vasculogenic mimicry in carcinogenesis and clinical applications. J Hematol Oncol. 2020;13:19.

    Article  PubMed  PubMed Central  Google Scholar 

  64. Quiros-Gonzalez I, Tomaszewski MR, Aitken SJ, Ansel-Bollepalli L, McDuffus LA, Gill M, et al. Optoacoustics delineates murine breast cancer models displaying angiogenesis and vascular mimicry. Br J Cancer. 2018;118:1098–106.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Luan YY, Liu ZM, Zhong JY, Yao RY, Yu HS. Effect of grape seed proanthocyanidins on tumor vasculogenic mimicry in human triple-negative breast cancer cells. Asian Pac J Cancer Prev. 2015;16:531–5.

    Article  PubMed  Google Scholar 

  66. Giuliano AE, Connolly JL, Edge SB, Mittendorf EA, Rugo HS, Solin LJ, et al. Breast cancer-major changes in the American Joint Committee on Cancer eighth edition cancer staging manual. CA Cancer J Clin. 2017;67:290–303.

    Article  PubMed  Google Scholar 

  67. Amin MB, Greene FL, Edge SB, Compton CC, Gershenwald JE, Brookland RK, et al. The eighth edition AJCC cancer staging manual: continuing to build a bridge from a population-based to a more “personalized” approach to cancer staging. CA Cancer J Clin. 2017;67:93–9.

    Article  PubMed  Google Scholar 

  68. Elston CW, Ellis IO. Pathological prognostic factors in breast cancer. I. The value of histological grade in breast cancer: experience from a large study with long-term follow-up. Histopathology. 1991;19:403–10.

    Article  CAS  PubMed  Google Scholar 

  69. Basu A, Ramamoorthi G, Jia Y, Faughn J, Wiener D, Awshah S, et al. Immunotherapy in breast cancer: current status and future directions. Adv Cancer Res. 2019;143:295–349.

    Article  PubMed  Google Scholar 

  70. Mansoori B, Mohammadi A, Davudian S, Shirjang S, Baradaran B. The different mechanisms of cancer drug resistance: a brief review. Adv Pharm Bull. 2017;7:339–48.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  71. Liedtke C, Mazouni C, Hess KR, André F, Tordai A, Mejia JA, et al. Response to neoadjuvant therapy and long-term survival in patients with triple-negative breast cancer. J Clin Oncol. 2008;26:1275–81.

    Article  PubMed  Google Scholar 

  72. Foulkes WD, Smith IE, Reis-Filho JS. Triple-negative breast cancer. N Engl J Med. 2010;363:1938–48.

    Article  CAS  PubMed  Google Scholar 

  73. Prieto-Vila M, Usuba W, Takahashi R-u, Shimomura I, Sasaki H, Ochiya T, et al. Single-cell analysis reveals a preexisting drug-resistant subpopulation in the luminal breast cancer subtype. Cancer Res. 2019;79:4412–25.

    Article  CAS  PubMed  Google Scholar 

  74. Gilreath C, Boerma M, Qin Z, Hudson MK, Wang S. The hypoxic microenvironment of breast cancer cells promotes resistance in radiation therapy. Front Oncol. 2021;10: 629422.

    Article  PubMed  PubMed Central  Google Scholar 

  75. Agrawal RK, Alhasso A, Barrett-Lee PJ, Bliss JM, Bliss P, Bloomfield D, et al. First results of the randomised UK FAST Trial of radiotherapy hypofractionation for treatment of early breast cancer (CRUKE/04/015). Radiother Oncol. 2011;100:93–100.

    Article  PubMed  Google Scholar 

  76. Mast ME, Vredeveld EJ, Credoe HM, van Egmond J, Heijenbrok MW, Hug EB, et al. Whole breast proton irradiation for maximal reduction of heart dose in breast cancer patients. Breast Cancer Res Treat. 2014;148:33–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  77. Bomken S, Fiser K, Heidenreich O, Vormoor J. Understanding the cancer stem cell. Br J Cancer. 2010;103:439–45.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  78. Qi XS, Pajonk F, McCloskey S, Low DA, Kupelian P, Steinberg M, et al. Radioresistance of the breast tumor is highly correlated to its level of cancer stem cell and its clinical implication for breast irradiation. Radiother Oncol. 2017;124:455–61.

    Article  PubMed  PubMed Central  Google Scholar 

  79. Kim IS, Gao Y, Welte T, Wang H, Liu J, Janghorban M, et al. Immuno-subtyping of breast cancer reveals distinct myeloid cell profiles and immunotherapy resistance mechanisms. Nat Cell Biol. 2019;21:1113–26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  80. Dandawate PR, Subramaniam D, Jensen RA, Anant S. Targeting cancer stem cells and signaling pathways by phytochemicals: novel approach for breast cancer therapy. Semin Cancer Biol. 2016;40–41:192–208.

    Article  PubMed  PubMed Central  Google Scholar 

  81. Tong M, Deng Z, Yang M, Xu C, Zhang X, Zhang Q, et al. Transcriptomic but not genomic variability confers phenotype of breast cancer stem cells. Cancer Commun. 2018;38:56.

    Article  Google Scholar 

  82. Zeng X, Liu C, Yao J, Wan H, Wan G, Li Y, et al. Breast cancer stem cells, heterogeneity, targeting therapies and therapeutic implications. Pharmacol Res. 2021;163: 105320.

    Article  CAS  PubMed  Google Scholar 

  83. Dent R, Trudeau M, Pritchard KI, Hanna WM, Kahn HK, Sawka CA, et al. Triple-negative breast cancer: clinical features and patterns of recurrence. Clin Cancer Res. 2007;13:4429–34.

    Article  PubMed  Google Scholar 

  84. Arnedos M, Vicier C, Loi S, Lefebvre C, Michiels S, Bonnefoi H, et al. Precision medicine for metastatic breast cancer—limitations and solutions. Nat Rev Clin Oncol. 2015;12:693–704.

    Article  CAS  PubMed  Google Scholar 

  85. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71:209–49.

    Article  PubMed  Google Scholar 

  86. Allemani C, Matsuda T, Di Carlo V, Harewood R, Matz M, Nikšić M, et al. Global surveillance of trends in cancer survival 2000–14 (CONCORD-3): analysis of individual records for 37 513 025 patients diagnosed with one of 18 cancers from 322 population-based registries in 71 countries. Lancet. 2018;391:1023–75.

    Article  PubMed  PubMed Central  Google Scholar 

  87. Piccart M, van’t Veer LJ, Poncet C, Lopes Cardozo JMN, Delaloge S, Pierga JY, et al. 70-gene signature as an aid for treatment decisions in early breast cancer: updated results of the phase 3 randomised MINDACT trial with an exploratory analysis by age. Lancet Oncol. 2021;22:476–88.

    Article  CAS  PubMed  Google Scholar 

  88. Hapach LA, Carey SP, Schwager SC, Taufalele PV, Wang W, Mosier JA, et al. Phenotypic heterogeneity and metastasis of breast cancer cells. Cancer Res. 2021;81:3649.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  89. Yu T, Wang C, Xie M, Zhu C, Shu Y, Tang J, et al. Heterogeneity of CTC contributes to the organotropism of breast cancer. Biomed Pharmacother. 2021;137: 111314.

    Article  CAS  PubMed  Google Scholar 

  90. Shin SU, Cho HM, Das R, Gil-Henn H, Ramakrishnan S, Al Bayati A, et al. Inhibition of vasculogenic mimicry and angiogenesis by an anti-EGFR IgG1-human endostatin-P125A fusion protein reduces triple negative breast cancer metastases. Cells. 2021;10:2904.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  91. Cleary AS, Leonard TL, Gestl SA, Gunther EJ. Tumour cell heterogeneity maintained by cooperating subclones in Wnt-driven mammary cancers. Nature. 2014;508:113–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  92. Marusyk A, Polyak K. Tumor heterogeneity: causes and consequences. Biochim Biophys Acta. 2010;1805:105–17.

    CAS  PubMed  Google Scholar 

  93. Kajita M, Fujita Y. EDAC: Epithelial defence against cancer-cell competition between normal and transformed epithelial cells in mammals. J Biochem. 2015;158:15–23.

    Article  CAS  PubMed  Google Scholar 

  94. Chen CL, Schroeder MC, Kango-Singh M, Tao C, Halder G. Tumor suppression by cell competition through regulation of the Hippo pathway. Proc Natl Acad Sci USA. 2012;109:484–9.

    Article  CAS  PubMed  Google Scholar 

  95. Sasaki A, Nagatake T, Egami R, Gu G, Takigawa I, Ikeda W, et al. Obesity suppresses cell-competition-mediated apical elimination of RasV12-transformed cells from epithelial tissues. Cell Rep. 2018;23:974–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  96. Pelham CJ, Nagane M, Madan E. Cell competition in tumor evolution and heterogeneity: merging past and present. Semin Cancer Biol. 2020;63:11–8.

    Article  CAS  PubMed  Google Scholar 

  97. Marusyk A, Tabassum DP, Altrock PM, Almendro V, Michor F, Polyak K. Non-cell-autonomous driving of tumour growth supports sub-clonal heterogeneity. Nature. 2014;514:54–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  98. Esparza-López J, Escobar-Arriaga E, Soto-Germes S, Ibarra-Sánchez MDJ. Breast cancer intra-tumor heterogeneity: one tumor, different entities. Rev Invest Clin. 2017;69:66–76.

    PubMed  Google Scholar 

  99. Madan E, Peixoto ML, Dimitrion P, Eubank TD, Yekelchyk M, Talukdar S, et al. Cell competition boosts clonal evolution and hypoxic selection in cancer. Trends Cell Biol. 2020;30:967–78.

    Article  CAS  PubMed  Google Scholar 

  100. Freischel AR, Damaghi M, Cunningham JJ, Ibrahim-Hashim A, Gillies RJ, Gatenby RA, et al. Frequency-dependent interactions determine outcome of competition between two breast cancer cell lines. Sci Rep. 2021;11:4908.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  101. Kaminska K, Szczylik C, Bielecka ZF, Bartnik E, Porta C, Lian F, et al. The role of the cell–cell interactions in cancer progression. J Cell Mol Med. 2015;19:283–96.

    Article  PubMed  PubMed Central  Google Scholar 

  102. Brady-Kalnay SM. Molecular mechanisms of cancer cell-cell interactions: cell–cell adhesion-dependent signaling in the tumor microenvironment. Cell Adh Migr. 2012;6:344–5.

    Article  PubMed  PubMed Central  Google Scholar 

  103. Xing F, Saidou J, Watabe K. Cancer associated fibroblasts (CAFs) in tumor microenvironment. Front Biosci. 2010;15:166–79.

    Article  CAS  PubMed Central  Google Scholar 

  104. Orimo A, Gupta PB, Sgroi DC, Arenzana-Seisdedos F, Delaunay T, Naeem R, et al. Stromal fibroblasts present in invasive human breast carcinomas promote tumor growth and angiogenesis through elevated SDF-1/CXCL12 secretion. Cell. 2005;121:335–48.

    Article  CAS  PubMed  Google Scholar 

  105. Pelon F, Bourachot B, Kieffer Y, Magagna I, Mermet-Meillon F, Bonnet I, et al. Cancer-associated fibroblast heterogeneity in axillary lymph nodes drives metastases in breast cancer through complementary mechanisms. Nat Commun. 2020;11:404.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  106. Lim H, Koh M, Jin H, Bae M, Lee SY, Kim KM, et al. Cancer-associated fibroblasts induce an aggressive phenotypic shift in non-malignant breast epithelial cells via interleukin-8 and S100A8. J Cell Physiol. 2021;236:7014–32.

    Article  CAS  PubMed  Google Scholar 

  107. Wen S, Hou Y, Fu L, Xi L, Yang D, Zhao M, et al. Cancer-associated fibroblast (CAF)-derived IL32 promotes breast cancer cell invasion and metastasis via integrin β3-p38 MAPK signalling. Cancer Lett. 2019;442:320–32.

    Article  CAS  PubMed  Google Scholar 

  108. Chiavarina B, Whitaker-Menezes D, Migneco G, Martinez-Outschoorn UE, Pavlides S, Howell A, et al. HIF1-alpha functions as a tumor promoter in cancer associated fibroblasts, and as a tumor suppressor in breast cancer cells Autophagy drives compartment-specific oncogenesis. Cell Cycle. 2010;9:3534–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  109. Eiro N, González L, Martínez-Ordoñez A, Fernandez-Garcia B, González LO, Cid S, et al. Cancer-associated fibroblasts affect breast cancer cell gene expression, invasion and angiogenesis. Cell Oncol (Dordr). 2018;41:369–78.

    Article  CAS  PubMed  Google Scholar 

  110. Ren J, Smid M, Iaria J, Salvatori DCF, van Dam H, Zhu HJ, et al. Cancer-associated fibroblast-derived Gremlin 1 promotes breast cancer progression. Breast Cancer Res. 2019;21:109.

    Article  PubMed  PubMed Central  Google Scholar 

  111. Donnarumma E, Fiore D, Nappa M, Roscigno G, Adamo A, Iaboni M, et al. Cancer-associated fibroblasts release exosomal microRNAs that dictate an aggressive phenotype in breast cancer. Oncotarget. 2017;8:19592–608.

    Article  PubMed  PubMed Central  Google Scholar 

  112. Wu HJ, Hao M, Yeo SK, Guan JL. FAK signaling in cancer-associated fibroblasts promotes breast cancer cell migration and metastasis by exosomal miRNAs-mediated intercellular communication. Oncogene. 2020;39:2539–49.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  113. Chen B, Sang Y, Song X, Zhang D, Wang L, Zhao W, et al. Exosomal miR-500a-5p derived from cancer-associated fibroblasts promotes breast cancer cell proliferation and metastasis through targeting USP28. Theranostics. 2021;11:3932–47.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  114. Loeffler M, Kruger JA, Niethammer AG, Reisfeld RA. Targeting tumor-associated fibroblasts improves cancer chemotherapy by increasing intratumoral drug uptake. J Clin Invest. 2006;116:1955–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  115. Lehuédé C, Li X, Dauvillier S, Vaysse C, Franchet C, Clement E, et al. Adipocytes promote breast cancer resistance to chemotherapy, a process amplified by obesity: role of the major vault protein (MVP). Breast Cancer Res. 2019;21:7.

    Article  PubMed  PubMed Central  Google Scholar 

  116. Togo S, Polanska UM, Horimoto Y, Orimo A. Carcinoma-associated fibroblasts are a promising therapeutic target. Cancers. 2013;5:149–69.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  117. Kalluri R. The biology and function of fibroblasts in cancer. Nat Rev Cancer. 2016;16:582–98.

    Article  CAS  PubMed  Google Scholar 

  118. Bai X, Yi M, Jiao Y, Chu Q, Wu K. Blocking TGF-β signaling to enhance the efficacy of immune checkpoint inhibitor. Onco Targets Ther. 2019;12:9527–38.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  119. Yi M, Zhang J, Li A, Niu M, Yan Y, Jiao Y, et al. The construction, expression, and enhanced anti-tumor activity of YM101: a bispecific antibody simultaneously targeting TGF-β and PD-L1. J Hematol Oncol. 2021;14:27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  120. Yi M, Niu M, Zhang J, Li S, Zhu S, Yan Y, et al. Combine and conquer: manganese synergizing anti-TGF-β/PD-L1 bispecific antibody YM101 to overcome immunotherapy resistance in non-inflamed cancers. J Hematol Oncol. 2021;14:146.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  121. Lin D, Chen X, Lin Z, Lin J, Liu Y, Liu D. Paper-supported co-culture system for dynamic investigations of the lung-tropic migration of breast cancer cells. Biomed Mater. 2021;16: 025028.

    Article  CAS  PubMed  Google Scholar 

  122. Flaberg E, Markasz L, Petranyi G, Stuber G, Dicső F, Alchihabi N, et al. High-throughput live-cell imaging reveals differential inhibition of tumor cell proliferation by human fibroblasts. Int J Cancer. 2011;128:2793–802.

    Article  CAS  PubMed  Google Scholar 

  123. Trimboli AJ, Cantemir-Stone CZ, Li F, Wallace JA, Merchant A, Creasap N, et al. Pten in stromal fibroblasts suppresses mammary epithelial tumours. Nature. 2009;461:1084–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  124. Bronisz A, Godlewski J, Wallace JA, Merchant AS, Nowicki MO, Mathsyaraja H, et al. Reprogramming of the tumour microenvironment by stromal PTEN-regulated miR-320. Nat Cell Biol. 2012;14:159–67.

    Article  CAS  Google Scholar 

  125. Sadlonova A, Mukherjee S, Bowe DB, Gault SR, Dumas NA, Van Tine BA, et al. Human breast fibroblasts inhibit growth of the MCF10AT xenograft model of proliferative breast disease. Am J Pathol. 2007;170:1064–76.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  126. Sung JS, Kang CW, Kang S, Jang Y, Chae YC, Kim BG, et al. ITGB4-mediated metabolic reprogramming of cancer-associated fibroblasts. Oncogene. 2020;39:664–76.

    Article  CAS  PubMed  Google Scholar 

  127. Iyengar P, Combs TP, Shah SJ, Gouon-Evans V, Pollard JW, Albanese C, et al. Adipocyte-secreted factors synergistically promote mammary tumorigenesis through induction of anti-apoptotic transcriptional programs and proto-oncogene stabilization. Oncogene. 2003;22:6408–23.

    Article  CAS  PubMed  Google Scholar 

  128. Wu Q, Li B, Li J, Sun S, Yuan J, Sun S. Cancer-associated adipocytes as immunomodulators in cancer. Biomark Res. 2021;9:2.

    Article  PubMed  PubMed Central  Google Scholar 

  129. Nalabolu MR, Palasamudram K, Jamil K. Adiponectin and leptin molecular actions and clinical significance in breast cancer. Int J Hematol Oncol Stem Cell Res. 2014;8:31–40.

    PubMed  PubMed Central  Google Scholar 

  130. He JY, Wei XH, Li SJ, Liu Y, Hu HL, Li ZZ, et al. Adipocyte-derived IL-6 and leptin promote breast cancer metastasis via upregulation of lysyl hydroxylase-2 expression. Cell Commun Signal. 2018;16:100.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  131. Park JY, Kang SE, Ahn KS, Um JY, Yang WM, Yun M, et al. Inhibition of the PI3K-AKT-mTOR pathway suppresses the adipocyte-mediated proliferation and migration of breast cancer cells. J Cancer. 2020;11:2552–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  132. Li SJ, Wei XH, Zhan XM, He JY, Zeng YQ, Tian XM, et al. Adipocyte-derived leptin promotes PAI-1 -mediated breast cancer metastasis in a STAT3/miR-34a dependent manner. Cancers. 2020;12:3864.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  133. Wang YY, Attané C, Milhas D, Dirat B, Dauvillier S, Guerard A, et al. Mammary adipocytes stimulate breast cancer invasion through metabolic remodeling of tumor cells. JCI Insight. 2017;2: e87489.

    Article  PubMed  PubMed Central  Google Scholar 

  134. Fujisaki K, Fujimoto H, Sangai T, Nagashima T, Sakakibara M, Shiina N, et al. Cancer-mediated adipose reversion promotes cancer cell migration via IL-6 and MCP-1. Breast Cancer Res Treat. 2015;150:255–63.

    Article  CAS  PubMed  Google Scholar 

  135. D’Esposito V, Liguoro D, Ambrosio MR, Collina F, Cantile M, Spinelli R, et al. Adipose microenvironment promotes triple negative breast cancer cell invasiveness and dissemination by producing CCL5. Oncotarget. 2016;7:24495–509.

    Article  PubMed  PubMed Central  Google Scholar 

  136. D’Esposito V, Passaretti F, Hammarstedt A, Liguoro D, Terracciano D, Molea G, et al. Adipocyte-released insulin-like growth factor-1 is regulated by glucose and fatty acids and controls breast cancer cell growth in vitro. Diabetologia. 2012;55:2811–22.

    Article  PubMed  PubMed Central  Google Scholar 

  137. Yao-Borengasser A, Monzavi-Karbassi B, Hedges RA, Rogers LJ, Kadlubar SA, Kieber-Emmons T. Adipocyte hypoxia promotes epithelial–mesenchymal transition-related gene expression and estrogen receptor-negative phenotype in breast cancer cells. Oncol Rep. 2015;33:2689–94.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  138. DeNardo D, Andreu P, Coussens LM. Interactions between lymphocytes and myeloid cells regulate pro- versus anti-tumor immunity. Cancer Metastasis Rev. 2010;29:309–16.

    Article  PubMed  PubMed Central  Google Scholar 

  139. Perrone M, Talarico G, Chiodoni C, Sangaletti S. Impact of immune cell heterogeneity on HER2+ breast cancer prognosis and response to therapy. Cancers. 2021;13:6352.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  140. Wang L, He T, Liu J, Tai J, Wang B, Chen Z, et al. Pan-cancer analysis reveals tumor-associated macrophage communication in the tumor microenvironment. Exp Hematol Oncol. 2021;10:31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  141. Lin Y, Xu J, Lan H. Tumor-associated macrophages in tumor metastasis: biological roles and clinical therapeutic applications. J Hematol Oncol. 2019;12:76.

    Article  PubMed  PubMed Central  Google Scholar 

  142. Xu T, Yu S, Zhang J, Wu S. Dysregulated tumor-associated macrophages in carcinogenesis, progression and targeted therapy of gynecological and breast cancers. J Hematol Oncol. 2021;14:181.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  143. Zhu S, Yi M, Wu Y, Dong B, Wu K. Roles of tumor-associated macrophages in tumor progression: implications on therapeutic strategies. Exp Hematol Oncol. 2021;10:60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  144. Mantovani A, Sozzani S, Locati M, Allavena P, Sica A. Macrophage polarization: tumor-associated macrophages as a paradigm for polarized M2 mononuclear phagocytes. Trends Immunol. 2002;23:549–55.

    Article  CAS  PubMed  Google Scholar 

  145. Jang J-H, Kim D-H, Lim JM, Lee JW, Jeong SJ, Kim KP, et al. Breast cancer cell-derived soluble CD44 promotes tumor progression by triggering macrophage IL1β production. Cancer Res. 2020;80:1342–56.

    Article  CAS  PubMed  Google Scholar 

  146. Liu Q, Hodge J, Wang J, Wang Y, Wang L, Singh U, et al. Emodin reduces breast cancer lung metastasis by suppressing macrophage-induced breast cancer cell epithelial–mesenchymal transition and cancer stem cell formation. Theranostics. 2020;10:8365–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  147. Hildenbrand R, Jansen C, Wolf G, Böhme B, Berger S, von Minckwitz G, et al. Transforming growth factor-beta stimulates urokinase expression in tumor-associated macrophages of the breast. Lab Invest. 1998;78:59–71.

    CAS  PubMed  Google Scholar 

  148. Hildenbrand R, Wolf G, Böhme B, Bleyl U, Steinborn A. Urokinase plasminogen activator receptor (CD87) expression of tumor-associated macrophages in ductal carcinoma in situ, breast cancer, and resident macrophages of normal breast tissue. J Leukoc Biol. 1999;66:40–9.

    Article  CAS  PubMed  Google Scholar 

  149. Xun J, Du L, Gao R, Shen L, Wang D, Kang L, et al. Cancer-derived exosomal miR-138-5p modulates polarization of tumor-associated macrophages through inhibition of KDM6B. Theranostics. 2021;11:6847–59.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  150. Chen Y, Zhang S, Wang Q, Zhang X. Tumor-recruited M2 macrophages promote gastric and breast cancer metastasis via M2 macrophage-secreted CHI3L1 protein. J Hematol Oncol. 2017;10:36.

    Article  PubMed  PubMed Central  Google Scholar 

  151. Guo J, Duan Z, Zhang C, Wang W, He H, Liu Y, et al. Mouse 4T1 breast cancer cell-derived exosomes induce proinflammatory cytokine production in macrophages via miR-183. J Immunol. 2020;205:2916–25.

    Article  CAS  PubMed  Google Scholar 

  152. Pardoll DM. The blockade of immune checkpoints in cancer immunotherapy. Nat Rev Cancer. 2012;12:252–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  153. Stovgaard ES, Dyhl-Polk A, Roslind A, Balslev E, Nielsen D. PD-L1 expression in breast cancer: expression in subtypes and prognostic significance: a systematic review. Breast Cancer Res Treat. 2019;174:571–84.

    Article  CAS  PubMed  Google Scholar 

  154. Rozenblit M, Huang R, Danziger N, Hegde P, Alexander B, Ramkissoon S, et al. Comparison of PD-L1 protein expression between primary tumors and metastatic lesions in triple negative breast cancers. J Immunother Cancer. 2020;8: e001558.

    Article  PubMed  PubMed Central  Google Scholar 

  155. Yi M, Niu M, Xu L, Luo S, Wu K. Regulation of PD-L1 expression in the tumor microenvironment. J Hematol Oncol. 2021;14:10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  156. Yi M, Xu L, Jiao Y, Luo S, Li A, Wu K. The role of cancer-derived microRNAs in cancer immune escape. J Hematol Oncol. 2020;13:25.

    Article  PubMed  PubMed Central  Google Scholar 

  157. Nanda R, Chow LQM, Dees EC, Berger R, Gupta S, Geva R, et al. Pembrolizumab in patients with advanced triple-negative breast cancer: phase Ib KEYNOTE-012 study. J Clin Oncol. 2016;34:2460–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  158. Voorwerk L, Slagter M, Horlings HM, Sikorska K, van de Vijver KK, de Maaker M, et al. Immune induction strategies in metastatic triple-negative breast cancer to enhance the sensitivity to PD-1 blockade: the TONIC trial. Nat Med. 2019;25:920–8.

    Article  CAS  PubMed  Google Scholar 

  159. Emens LA, Adams S, Cimino-Mathews A, Disis ML, Gatti-Mays ME, Ho AY, et al. Society for Immunotherapy of Cancer (SITC) clinical practice guideline on immunotherapy for the treatment of breast cancer. J Immunother Cancer. 2021;9: e002597.

    Article  PubMed  PubMed Central  Google Scholar 

  160. Yi M, Zheng X, Niu M, Zhu S, Ge H, Wu K. Combination strategies with PD-1/PD-L1 blockade: current advances and future directions. Mol Cancer. 2022;21:28.

    Article  PubMed  PubMed Central  Google Scholar 

  161. Fang W, Zhou T, Shi H, Yao M, Zhang D, Qian H, et al. Progranulin induces immune escape in breast cancer via up-regulating PD-L1 expression on tumor-associated macrophages (TAMs) and promoting CD8(+) T cell exclusion. J Exp Clin Cancer Res. 2021;40:4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  162. Yao X, Tu Y, Xu Y, Guo Y, Yao F, Zhang X. Endoplasmic reticulum stress-induced exosomal miR-27a-3p promotes immune escape in breast cancer via regulating PD-L1 expression in macrophages. J Cell Mol Med. 2020;24:9560–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  163. Lee A, Lim S, Oh J, Lim J, Yang Y, Lee M-S, et al. NDRG2 expression in breast cancer cells downregulates PD-L1 expression and restores T cell proliferation in tumor-coculture. Cancers. 2021;13:6112.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  164. Uddin MB, Roy KR, Hill RA, Roy SC, Gu X, Li L, et al. p53 missense mutant G242A subverts natural killer cells in sheltering mouse breast cancer cells against immune rejection. Exp Cell Res. 2022;417: 113210.

    Article  CAS  PubMed  Google Scholar 

  165. Fu T, Dai L-J, Wu S-Y, Xiao Y, Ma D, Jiang Y-Z, et al. Spatial architecture of the immune microenvironment orchestrates tumor immunity and therapeutic response. J Hematol Oncol. 2021;14:98.

    Article  PubMed  PubMed Central  Google Scholar 

  166. Ding L, Shunkwiler LB, Harper NW, Zhao Y, Hinohara K, Huh SJ, et al. Deletion of Cdkn1b in ACI rats leads to increased proliferation and pregnancy-associated changes in the mammary gland due to perturbed systemic endocrine environment. PLoS Genet. 2019;15: e1008002.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  167. Toillon RA, Chopin V, Jouy N, Fauquette W, Boilly B, Le Bourhis X. Normal breast epithelial cells induce p53-dependent apoptosis and p53-independent cell cycle arrest of breast cancer cells. Breast Cancer Res Treat. 2002;71:269–80.

    Article  CAS  PubMed  Google Scholar 

  168. Li S, Zhou J, Wu H, Lu Q, Tai Y, Liu Q, et al. Oncogenic transformation of normal breast epithelial cells co-cultured with cancer cells. Cell Cycle. 2018;17:2027–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  169. Jo SH, Heo WH, Son H-Y, Quan M, Hong BS, Kim JH, et al. S100A8/A9 mediate the reprograming of normal mammary epithelial cells induced by dynamic cell–cell interactions with adjacent breast cancer cells. Sci Rep. 2021;11:1337.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  170. Li S, Yi M, Dong B, Jiao Y, Luo S, Wu K. The roles of exosomes in cancer drug resistance and its therapeutic application. Clin Transl Med. 2020;10: e257.

    Article  PubMed  PubMed Central  Google Scholar 

  171. Melo SA, Sugimoto H, O’Connell JT, Kato N, Villanueva A, Vidal A, et al. Cancer exosomes perform cell-independent microRNA biogenesis and promote tumorigenesis. Cancer Cell. 2014;26:707–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  172. Neelakantan D, Zhou H, Oliphant MUJ, Zhang X, Simon LM, Henke DM, et al. EMT cells increase breast cancer metastasis via paracrine GLI activation in neighbouring tumour cells. Nat Commun. 2017;8:15773.

    Article  PubMed  PubMed Central  Google Scholar 

  173. Liu Q, Li A, Yu S, Qin S, Han N, Pestell RG, et al. DACH1 antagonizes CXCL8 to repress tumorigenesis of lung adenocarcinoma and improve prognosis. J Hematol Oncol. 2018;11:53.

    Article  PubMed  PubMed Central  Google Scholar 

  174. Wu K, Katiyar S, Li A, Liu M, Ju X, Popov VM, et al. Dachshund inhibits oncogene-induced breast cancer cellular migration and invasion through suppression of interleukin-8. Proc Natl Acad Sci USA. 2008;105:6924–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  175. Dong B, Yi M, Luo S, Li A, Wu K. RDGN-based predictive model for the prognosis of breast cancer. Exp Hematol Oncol. 2020;9:13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  176. Liu Y, Han N, Zhou S, Zhou R, Yuan X, Xu H, et al. The DACH/EYA/SIX gene network and its role in tumor initiation and progression. Int J Cancer. 2016;138:1067–75.

    Article  CAS  PubMed  Google Scholar 

  177. Kong D, Zhou H, Neelakantan D, Hughes CJ, Hsu JY, Srinivasan RR, et al. VEGF-C mediates tumor growth and metastasis through promoting EMT-epithelial breast cancer cell crosstalk. Oncogene. 2021;40:964–79.

    Article  CAS  PubMed  Google Scholar 

  178. Kung HN, Marks JR, Chi JT. Glutamine synthetase is a genetic determinant of cell type-specific glutamine independence in breast epithelia. PLoS Genet. 2011;7: e1002229.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  179. Kennedy KM, Scarbrough PM, Ribeiro A, Richardson R, Yuan H, Sonveaux P, et al. Catabolism of exogenous lactate reveals it as a legitimate metabolic substrate in breast cancer. PLoS ONE. 2013;8: e75154.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  180. Deming SL, Nass SJ, Dickson RB, Trock BJ. C-myc amplification in breast cancer: a meta-analysis of its occurrence and prognostic relevance. Br J Cancer. 2000;83:1688–95.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  181. Johnston LA. Socializing with MYC: cell competition in development and as a model for premalignant cancer. Cold Spring Harb Perspect Med. 2014;4: a014274.

    Article  PubMed  PubMed Central  Google Scholar 

  182. Patel MS, Shah HS, Shrivastava N. c-Myc-dependent cell competition in human cancer cells. J Cell Biochem. 2017;118:1782–91.

    Article  CAS  PubMed  Google Scholar 

  183. Di Giacomo S, Sollazzo M, de Biase D, Ragazzi M, Bellosta P, Pession A, et al. Human cancer cells signal their competitive fitness through MYC activity. Sci Rep. 2017;7:12.

    Article  Google Scholar 

  184. Li K, Sun X, Zha R, Liu S, Feng Y, Sano T, et al. Counterintuitive production of tumor-suppressive secretomes from Oct4- and c-Myc-overexpressing tumor cells and MSCs. Theranostics. 2022;12:3084–103.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  185. Gupta N, Jung K, Wu C, Alshareef A, Alqahtani H, Damaraju S, et al. High Myc expression and transcription activity underlies intra-tumoral heterogeneity in triple-negative breast cancer. Oncotarget. 2017;8:28101–15.

    Article  PubMed  PubMed Central  Google Scholar 

  186. Cho SW, Xu J, Sun R, Mumbach MR, Carter AC, Chen YG, et al. Promoter of lncRNA gene PVT1 is a tumor-suppressor DNA boundary element. Cell. 2018;173:1398-412.e22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  187. Madan E, Pelham CJ, Nagane M, Parker TM, Canas-Marques R, Fazio K, et al. Flower isoforms promote competitive growth in cancer. Nature. 2019;572:260–4.

    Article  CAS  PubMed  Google Scholar 

  188. Sun Q, Luo TZ, Ren YX, Florey O, Shirasawa S, Sasazuki T, et al. Competition between human cells by entosis. Cell Res. 2014;24:1299–310.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  189. Gupta PB, Fillmore CM, Jiang GZ, Shapira SD, Tao K, Kuperwasser C, et al. Stochastic state transitions give rise to phenotypic equilibrium in populations of cancer cells. Cell. 2011;146:633–44.

    Article  CAS  PubMed  Google Scholar 

  190. Hausser J, Szekely P, Bar N, Zimmer A, Sheftel H, Caldas C, et al. Tumor diversity and the trade-off between universal cancer tasks. Nat Commun. 2019;10:5423.

    Article  PubMed  PubMed Central  Google Scholar 

  191. Hausser J, Alon U. Tumour heterogeneity and the evolutionary trade-offs of cancer. Nat Rev Cancer. 2020;20:247–57.

    Article  CAS  PubMed  Google Scholar 

  192. Lim B, Lin Y, Navin N. Advancing cancer research and medicine with single-cell genomics. Cancer Cell. 2020;37:456–70.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  193. Gerlinger M, Rowan AJ, Horswell S, Larkin J, Endesfelder D, Gronroos E, et al. Intratumor heterogeneity and branched evolution revealed by multiregion sequencing. N Engl J Med. 2012;366:883–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  194. Sanger F, Nicklen S, Coulson AR. DNA sequencing with chain-terminating inhibitors. Proc Natl Acad Sci USA. 1977;74:5463–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  195. Heather JM, Chain B. The sequence of sequencers: the history of sequencing DNA. Genomics. 2016;107:1–8.

    Article  CAS  PubMed  Google Scholar 

  196. Lee I, Razaghi R, Gilpatrick T, Molnar M, Gershman A, Sadowski N, et al. Simultaneous profiling of chromatin accessibility and methylation on human cell lines with nanopore sequencing. Nat Methods. 2020;17:1191–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  197. Aslam B, Basit M, Nisar MA, Khurshid M, Rasool MH. Proteomics: technologies and their applications. J Chromatogr Sci. 2017;55:182–96.

    Article  CAS  PubMed  Google Scholar 

  198. Timp W, Timp G. Beyond mass spectrometry, the next step in proteomics. Sci Adv. 2020;6:eaax8978.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  199. Yusko EC, Bruhn BR, Eggenberger OM, Houghtaling J, Rollings RC, Walsh NC, et al. Real-time shape approximation and fingerprinting of single proteins using a nanopore. Nat Nanotechnol. 2017;12:360–7.

    Article  CAS  PubMed  Google Scholar 

  200. Lei Y, Tang R, Xu J, Wang W, Zhang B, Liu J, et al. Applications of single-cell sequencing in cancer research: progress and perspectives. J Hematol Oncol. 2021;14:91.

    Article  PubMed  PubMed Central  Google Scholar 

  201. Gao R, Davis A, McDonald TO, Sei E, Shi X, Wang Y, et al. Punctuated copy number evolution and clonal stasis in triple-negative breast cancer. Nat Genet. 2016;48:1119–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  202. Grosselin K, Durand A, Marsolier J, Poitou A, Marangoni E, Nemati F, et al. High-throughput single-cell ChIP-seq identifies heterogeneity of chromatin states in breast cancer. Nat Genet. 2019;51:1060–6.

    Article  CAS  PubMed  Google Scholar 

  203. Zheng GXY, Terry JM, Belgrader P, Ryvkin P, Bent ZW, Wilson R, et al. Massively parallel digital transcriptional profiling of single cells. Nat Commun. 2017;8:14049.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  204. Yu TJ, Ma D, Liu YY, Xiao Y, Gong Y, Jiang YZ, et al. Bulk and single-cell transcriptome profiling reveal the metabolic heterogeneity in human breast cancers. Mol Ther. 2021;29:2350–65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  205. Karaayvaz M, Cristea S, Gillespie SM, Patel AP, Mylvaganam R, Luo CC, et al. Unravelling subclonal heterogeneity and aggressive disease states in TNBC through single-cell RNA-seq. Nat Commun. 2018;9:3588.

    Article  PubMed  PubMed Central  Google Scholar 

  206. Chung W, Eum HH, Lee HO, Lee KM, Lee HB, Kim KT, et al. Single-cell RNA-seq enables comprehensive tumour and immune cell profiling in primary breast cancer. Nat Commun. 2017;8:12.

    Article  Google Scholar 

  207. Jiang K, Dong M, Li C, Sheng J. Unraveling heterogeneity of tumor cells and microenvironment and its clinical implications for triple negative breast cancer. Front Oncol. 2021;11: 557477.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  208. Li X, Liu L, Goodall GJ, Schreiber A, Xu T, Li J, et al. A novel single-cell based method for breast cancer prognosis. PLoS Comput Biol. 2020;16: e1008133.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  209. Brechbuhl HM, Vinod-Paul K, Gillen AE, Kopin EG, Gibney K, Elias AD, et al. Analysis of circulating breast cancer cell heterogeneity and interactions with peripheral blood mononuclear cells. Mol Carcinog. 2020;59:1129–39.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  210. Li Q, Cui S, Xu Y, Wang Y, Jin F, Si H, et al. Consecutive sorting and phenotypic counting of CTCs by an optofluidic flow cytometer. Anal Chem. 2019;91:14133–40.

    Article  CAS  PubMed  Google Scholar 

  211. Yang M, Li J, Zhao C, Xiao H, Fang X, Zheng J. LC-Q-TOF-MS/MS detection of food flavonoids: principle, methodology, and applications. Crit Rev Food Sci Nutr. 2021. https://doi.org/10.1080/10408398.2021.1993128.

    Article  PubMed  Google Scholar 

  212. Qu J, Ke F, Liu Z, Yang X, Li X, Xu H, et al. Uncovering the mechanisms of dandelion against triple-negative breast cancer using a combined network pharmacology, molecular pharmacology and metabolomics approach. Phytomedicine. 2022;99: 153986.

    Article  CAS  PubMed  Google Scholar 

  213. Casasent AK, Schalck A, Gao R, Sei E, Long A, Pangburn W, et al. Multiclonal invasion in breast tumors identified by topographic single cell sequencing. Cell. 2018;172:205-17.e12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  214. Wang K, Xiao Z, Yan Y, Ye R, Hu M, Bai S, et al. Simple oligonucleotide-based multiplexing of single-cell chromatin accessibility. Mol Cell. 2021;81:4319-32.e10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  215. Stahl PL, Salmen F, Vickovic S, Lundmark A, Navarro JF, Magnusson J, et al. Visualization and analysis of gene expression in tissue sections by spatial transcriptomics. Science. 2016;353:78–82.

    Article  CAS  PubMed  Google Scholar 

  216. Wei W, Dai W, Yang F, Lu H, Zhang K, Xing Y, et al. Spatially resolved, error-robust multiplexed microRNA profiling in single living cells. Angew Chem Int Ed Engl. 2022;61: e202116909.

    Article  CAS  PubMed  Google Scholar 

  217. Keren L, Bosse M, Thompson S, Risom T, Vijayaragavan K, McCaffrey E, et al. MIBI-TOF: a multiplexed imaging platform relates cellular phenotypes and tissue structure. Sci Adv. 2019;5:eaax5851.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  218. Risom T, Glass DR, Averbukh I, Liu CC, Baranski A, Kagel A, et al. Transition to invasive breast cancer is associated with progressive changes in the structure and composition of tumor stroma. Cell. 2022;185:299-310.e18.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  219. Giesen C, Wang HA, Schapiro D, Zivanovic N, Jacobs A, Hattendorf B, et al. Highly multiplexed imaging of tumor tissues with subcellular resolution by mass cytometry. Nat Methods. 2014;11:417–22.

    Article  CAS  PubMed  Google Scholar 

  220. Bartoschek M, Oskolkov N, Bocci M, Lövrot J, Larsson C, Sommarin M, et al. Spatially and functionally distinct subclasses of breast cancer-associated fibroblasts revealed by single cell RNA sequencing. Nat Commun. 2018;9:5150.

    Article  PubMed  PubMed Central  Google Scholar 

  221. Wagner J, Rapsomaniki MA, Chevrier S, Anzeneder T, Langwieder C, Dykgers A, et al. A single-cell atlas of the tumor and immune ecosystem of human breast cancer. Cell. 2019;177:1330-45.e18.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  222. Hu Q, Hong Y, Qi P, Lu G, Mai X, Xu S, et al. Atlas of breast cancer infiltrated B-lymphocytes revealed by paired single-cell RNA-sequencing and antigen receptor profiling. Nat Commun. 2021;12:2186.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  223. Andersson A, Larsson L, Stenbeck L, Salmén F, Ehinger A, Wu SZ, et al. Spatial deconvolution of HER2-positive breast cancer delineates tumor-associated cell type interactions. Nat Commun. 2021;12:6012.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  224. Aldhaeebi MA, Alzoubi K, Almoneef TS, Bamatraf SM, Attia H, Ramahi OM. Review of microwaves techniques for breast cancer detection. Sensors. 2020;20:2390.

    Article  PubMed  PubMed Central  Google Scholar 

  225. Kjaer A. Molecular imaging of cancer using PET and SPECT. Adv Exp Med Biol. 2006;587:277–84.

    Article  PubMed  Google Scholar 

  226. Yang CT, Ghosh KK, Padmanabhan P, Langer O, Liu J, Eng DNC, et al. PET-MR and SPECT-MR multimodality probes: development and challenges. Theranostics. 2018;8:6210–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  227. Lambin P, Rios-Velazquez E, Leijenaar R, Carvalho S, van Stiphout RG, Granton P, et al. Radiomics: extracting more information from medical images using advanced feature analysis. Eur J Cancer. 2012;48:441–6.

    Article  PubMed  PubMed Central  Google Scholar 

  228. Gillies RJ, Kinahan PE, Hricak H. Radiomics: images are more than pictures, they are data. Radiology. 2016;278:563–77.

    Article  PubMed  Google Scholar 

  229. Tagliafico AS, Piana M, Schenone D, Lai R, Massone AM, Houssami N. Overview of radiomics in breast cancer diagnosis and prognostication. Breast. 2020;49:74–80.

    Article  PubMed  Google Scholar 

  230. Crivelli P, Ledda RE, Parascandolo N, Fara A, Soro D, Conti M.