BackgroundCancer stem-like cells (CSCs) have been extensively researched as the primary drivers of therapy resistance and tumor relapse in patients with breast cancer. However, due to lack of specific molecular markers, increased phenotypic plasticity and no clear clinicopathological features, the assessment of CSCs presence and functionality in solid tumors is challenging. While several potential markers, such as CD24/CD44, have been proposed, the extent to which they truly represent the stem cell potential of tumors or merely provide static snapshots is still a subject of controversy. Recent studies have highlighted the crucial role of the tumor microenvironment (TME) in influencing the CSC phenotype in breast cancer. The interplay between the tumor and TME induces significant changes in the cancer cell phenotype, leading to the acquisition of CSC characteristics, therapeutic resistance, and metastatic spread. Simultaneously, CSCs actively shape their microenvironment by evading immune surveillance and attracting stromal cells that support tumor progression.MethodsIn this study, we associated in vitro mammosphere formation assays with bulk tumor microarray profiling and deconvolution algorithms to map CSC functionality and the microenvironmental landscape in a large cohort of 125 breast tumors.ResultsWe found that the TME score was a significant factor associated with CSC functionality. CSC-rich tumors were characterized by an immune-suppressed TME, while tumors devoid of CSC potential exhibited high immune infiltration and activation of pathways involved in the immune response. Gene expression analysis revealed IFNG, CXCR5, CD40LG, TBX21 and IL2RG to be associated with the CSC phenotype and also displayed prognostic value for patients with breast cancer.ConclusionThese results suggest that the characterization of CSCs content and functionality in tumors can be used as an attractive strategy to fine-tune treatments and guide clinical decisions to improve patients therapy response.