Phenotypic plasticity is a major factor in tumor heterogeneity and treatment resistance. In particular, cancer stem cells (CSCs) represent a small subpopulation within tumors that possesses self-renewal and tumor-forming capabilities. Understanding reprogramming, maintenance, and lineage properties of CSCs requires dedicated tools to disentangle the respective influences of phenotypic inheritance and cell-cell interactions. Here, we set up ultra-wide field microscopy to image breast cancer cell lines expressing a stemness fluorescent reporter over several days. The fluorescent reporter distinguishes three phenotypes: CSCs, cancer differentiated cells (CDCs), and intermediate/transiting cancer cells (iCCs). Spatial statistics indicate significant zonation in which CSCs cluster together and are spatially separated from CDCs, forming patterns resembling niches. Surprisingly, single-cell time series reveal spontaneous reprogramming events from CDC to CSC even in unperturbed populations. We identify that such transitions are prone to arise during the cell cycle. Moreover, lineage analysis shows that the phenotype is partially inherited from ancestor cells. However, such heredity is not sufficient to explain the spatial properties of the cell population, which also depend on cell-cell interactions. Indeed, we find that phenotypic transitions of cancer cells are influenced by the phenotypic state of neighboring cells. Reprogramming into CSCs is respectively promoted and inhibited by the presence of CSCs and CDCs in the neighborhood. Altogether, our results disentangle how phenotypic inheritance and intercellular interactions orchestrate the spatio-temporal self-organization of cancer cell heterogeneity, maintaining a subpopulation of CSCs within niches.

Cancer stem cell (CSC) has paved the way to many fundamental and translational studies. Recent studies have highlighted differentiated breast cancer cells (non-CSCs) switching phenotype to CSCs in response to various stimuli, depicting the existence of cancer stem cell plasticity. Although strategies to reduce the phenotypic plasticity of non-CSCs into CSCs are likely to prevent treatment-resilient cancer cells driving recurrence, most phenotypic plasticity mechanisms involve Notch, Wnt or MAPK signaling pathways. In this study, breast cancer cells were irradiated to identify soluble reprogramming factors. Using conditioned medias, protein arrays analyses, flow cytometry and in cellulo/in vivo functional assays, we demonstrated, for the first time, that radiation-induced chemokine expression, especially CXCL1 and CCL5 and their receptors CXCR2, CCR1 and CCR5, stimulates reprogramming of breast non-CSCs into CSCs. Treatment of non-CSCs with recombinant CXCL1 and CCL5 is sufficient to induce cell reprogramming, while their inhibition can be used to prevent reprogramming and sensitize tumor to radiation. Moreover, analysis of gene expression profiles from 38 public merged databases demonstrated that combined over-expression of CXCL1/CXCR2, CCL5/CCR1 or CCL5/CCR5 has a poorer prognosis in patients treated with radiotherapy, suggesting a promising way for patient stratification, where individuals with elevated cytokine levels could benefit from radiotherapy in conjunction with cytokine inhibitors. Taken together, our findings provide a rationale to consider these axes as potential targets and predictive biomarkers in breast cancer patients.

MOTIVATION:Nowadays, epigenetic gene regulations are studied in each part of the biology, from embryonic development to diseases such as cancers and neurodegenerative disorders. Currently, to quantify and compare CpG methylation levels of a specific region of interest, the most accessible technique is the bisulfite sequencing PCR (BSP). However, no existing user-friendly tool is able to analyze data from all approaches of BSP. Therefore, the most convenient way to process results from the direct sequencing of PCR products (direct-BSP) is to manually analyze the chromatogram traces, which is a repetitive and prone to error task.RESULTS:Here, we implement a new R-based tool, called ABSP for analysis of bisulfite sequencing PCR, providing a complete analytic process of both direct-BSP and cloning-BSP data. It uses the raw sequencing trace files (.ab1) as input to compute and compare CpG methylation percentages. It is fully automated and includes a user-friendly interface as a built-in R shiny app, quality control steps and generates publication-ready graphics.AVAILABILITY AND IMPLEMENTATION:The ABSP tool and associated data are available on GitHub at https://github.com/ABSP-methylation-tool/ABSP.SUPPLEMENTARY INFORMATION:Supplementary data are available at Bioinformatics online.

ABSTRACTWe developed an experimental methodology to monitor response dynamics of single human cells to hydrogen peroxide. Our approach is based on fluidic control of both magnitude and time-evolution of the external perturbation, and on high-throughput imaging of intracellular fluorescent redox reporters. We applied step stimulus to MCF7 cells with hydrogen peroxide concentration in the range of 10 to 500μM. First, our data highlights dynamic adaptation of Reactive Oxygen Species (ROS) scavenging system at several time scales. Nicotinamide Adenine Dinucleotide Phosphate (NADPH) level is rapidly restored within 10 min after a transient decrease, while glutathione (GSH) redox potential is slowly driven back toward pre-stimulus level (within one hour). Extra-cellular glucose is necessary for adaptation of both NADPH level and GSH system. Second, our results also reveal large cell-to-cell variability in the dynamic response to external ROS. Our experimental approach is readily usable to monitor other cellular redox systems such as thioredoxins. As response-heterogeneity can lead to fractional killing, we finally discuss how our methodology can be an opportunity to link dynamics of ROS scavenging and cellular fate such as cell death.