- Review
- Open access
- Published:
Breast cancer heterogeneity and its implication in personalized precision therapy
Experimental Hematology & Oncology volume 12, Article number: 3 (2023)
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.
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.
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.
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.
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
05 January 2024
The original online version of this article was revised: Grant number is corrected
23 January 2024
A Correction to this paper has been published: https://doi.org/10.1186/s40164-024-00472-z
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
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.
Siegel RL, Miller KD, Fuchs HE, Jemal A. Cancer statistics, 2021. CA Cancer J Clin. 2021;71:7–33.
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.
Visvader JE. Cells of origin in cancer. Nature. 2011;469:314–22.
Dagogo-Jack I, Shaw AT. Tumour heterogeneity and resistance to cancer therapies. Nat Rev Clin Oncol. 2018;15:81–94.
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.
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.
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.
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.
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.
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.
Pott S. Simultaneous measurement of chromatin accessibility, DNA methylation, and nucleosome phasing in single cells. Elife. 2017;6: e23203.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Tang ZR, Zhang R, Lian ZX, Deng SL, Yu K. Estrogen-receptor expression and function in female reproductive disease. Cells. 2019;8:1123.
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.
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.
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.
Gao JJ, Swain SM. Luminal A breast cancer and molecular assays: a review. Oncologist. 2018;23:556–65.
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.
Loibl S, Gianni L. HER2-positive breast cancer. Lancet. 2017;389:2415–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.
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.
Fougner C, Bergholtz H, Norum JH, Sørlie T. Re-definition of claudin-low as a breast cancer phenotype. Nat Commun. 2020;11:1787.
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.
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.
Turner NH, Di Leo A. HER2 discordance between primary and metastatic breast cancer: assessing the clinical impact. Cancer Treat Rev. 2013;39:947–57.
Kroigard AB, Larsen MJ, Thomassen M, Kruse TA. Molecular concordance between primary breast cancer and matched metastases. Breast J. 2016;22:420–30.
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.
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.
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.
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.
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.
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.
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.
Dias AS, Almeida CR, Helguero LA, Duarte IF. Metabolic crosstalk in the breast cancer microenvironment. Eur J Cancer. 2019;121:154–71.
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.
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.
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.
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.
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.
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.
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.
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.
Gupta GP, Massagué J. Cancer metastasis: building a framework. Cell. 2006;127:679–95.
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.
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.
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.
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.
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.
Micalizzi DS, Maheswaran S, Haber DA. A conduit to metastasis: circulating tumor cell biology. Genes Dev. 2017;31:1827–40.
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.
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.
Folberg R, Maniotis AJ. Vasculogenic mimicry. APMIS. 2004;112:508–25.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Foulkes WD, Smith IE, Reis-Filho JS. Triple-negative breast cancer. N Engl J Med. 2010;363:1938–48.
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.
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.
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.
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.
Bomken S, Fiser K, Heidenreich O, Vormoor J. Understanding the cancer stem cell. Br J Cancer. 2010;103:439–45.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Marusyk A, Polyak K. Tumor heterogeneity: causes and consequences. Biochim Biophys Acta. 2010;1805:105–17.
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.
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.
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.
Pelham CJ, Nagane M, Madan E. Cell competition in tumor evolution and heterogeneity: merging past and present. Semin Cancer Biol. 2020;63:11–8.
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.
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.
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.
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.
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.
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.
Xing F, Saidou J, Watabe K. Cancer associated fibroblasts (CAFs) in tumor microenvironment. Front Biosci. 2010;15:166–79.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Togo S, Polanska UM, Horimoto Y, Orimo A. Carcinoma-associated fibroblasts are a promising therapeutic target. Cancers. 2013;5:149–69.
Kalluri R. The biology and function of fibroblasts in cancer. Nat Rev Cancer. 2016;16:582–98.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Wu Q, Li B, Li J, Sun S, Yuan J, Sun S. Cancer-associated adipocytes as immunomodulators in cancer. Biomark Res. 2021;9:2.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Lin Y, Xu J, Lan H. Tumor-associated macrophages in tumor metastasis: biological roles and clinical therapeutic applications. J Hematol Oncol. 2019;12:76.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Pardoll DM. The blockade of immune checkpoints in cancer immunotherapy. Nat Rev Cancer. 2012;12:252–64.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Johnston LA. Socializing with MYC: cell competition in development and as a model for premalignant cancer. Cold Spring Harb Perspect Med. 2014;4: a014274.
Patel MS, Shah HS, Shrivastava N. c-Myc-dependent cell competition in human cancer cells. J Cell Biochem. 2017;118:1782–91.
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.
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.
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.
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.
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.
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.
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.
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.
Hausser J, Alon U. Tumour heterogeneity and the evolutionary trade-offs of cancer. Nat Rev Cancer. 2020;20:247–57.
Lim B, Lin Y, Navin N. Advancing cancer research and medicine with single-cell genomics. Cancer Cell. 2020;37:456–70.
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.
Sanger F, Nicklen S, Coulson AR. DNA sequencing with chain-terminating inhibitors. Proc Natl Acad Sci USA. 1977;74:5463–7.
Heather JM, Chain B. The sequence of sequencers: the history of sequencing DNA. Genomics. 2016;107:1–8.
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.
Aslam B, Basit M, Nisar MA, Khurshid M, Rasool MH. Proteomics: technologies and their applications. J Chromatogr Sci. 2017;55:182–96.
Timp W, Timp G. Beyond mass spectrometry, the next step in proteomics. Sci Adv. 2020;6:eaax8978.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Aldhaeebi MA, Alzoubi K, Almoneef TS, Bamatraf SM, Attia H, Ramahi OM. Review of microwaves techniques for breast cancer detection. Sensors. 2020;20:2390.
Kjaer A. Molecular imaging of cancer using PET and SPECT. Adv Exp Med Biol. 2006;587:277–84.
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.
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.
Gillies RJ, Kinahan PE, Hricak H. Radiomics: images are more than pictures, they are data. Radiology. 2016;278:563–77.
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.
Crivelli P, Ledda RE, Parascandolo N, Fara A, Soro D, Conti M.