Abstract Breast cancer is a heterogenous disease that can be classified into multiple subtypes including the most aggressive basal-like and triple-negative subtypes. Understanding the heterogeneity within the normal mammary basal epithelial cells holds the key to inform us about basal-like cancer cell differentiation dynamics as well as potential cells of origin. Although it is known that the mammary basal compartment contains small pools of stem cells that fuel normal tissue morphogenesis and regeneration, a comprehensive yet focused analysis of the transcriptional makeup of the basal cells is lacking. We used single-cell RNA-sequencing and multiplexed RNA in-situ hybridization to characterize mammary basal cell heterogeneity. We used bioinformatic and computational pipelines to characterize the molecular features as well as predict differentiation dynamics and cell–cell communications of the newly identified basal cell states. We used genetic cell labeling to map the in vivo fates of cells in one of these states. We identified four major distinct transcriptional states within the mammary basal cells that exhibit gene expression signatures suggestive of different functional activity and metabolic preference. Our in vivo labeling and ex vivo organoid culture data suggest that one of these states, marked by Egr2 expression, represents a dynamic transcriptional state that all basal cells transit through during pubertal mammary morphogenesis. Our study provides a systematic approach to understanding the molecular heterogeneity of mammary basal cells and identifies previously unknown dynamics of basal cell transcriptional states.
more »
« less
Positional influence on cellular transcriptional identity revealed through spatially segmented single-cell transcriptomics
Single-cell RNA sequencing (scRNA-seq) is a powerful technique for describing cell states. Identifying the spatial arrangement of these states in tissues remains challenging, with the existing methods requiring niche methodologies and expertise. Here, we describe segmentation by exogenous perfusion (SEEP), a rapid and integrated method to link surface proximity and environment accessibility to transcriptional identity within three-dimensional (3D) disease models. The method utilizes the steady-state diffusion kinetics of a fluores- cent dye to establish a gradient along the radial axis of disease models. Classification of sample layers based on dye accessibility enables dissociated and sorted cells to be characterized by transcriptomic and regional identities. Using SEEP, we analyze spheroid, organoid, and in vivo tumor models of high-grade serous ovarian cancer (HGSOC). The results validate long-standing beliefs about the relationship between cell state and position while revealing new concepts regarding how spatially unique microenvironments influence the identity of individual cells within tumors.
more »
« less
- Award ID(s):
- 2011754
- PAR ID:
- 10500615
- Publisher / Repository:
- Cell Systems
- Date Published:
- Journal Name:
- Cell Systems
- Volume:
- 14
- Issue:
- 6
- ISSN:
- 2405-4712
- Page Range / eLocation ID:
- 464 to 481.e7
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Breast and mammary epithelial cells experience different local environments during tissue development and tumorigenesis. Microenvironmental heterogeneity gives rise to distinct cell regulatory states whose identity and importance are just beginning to be appreciated. Cellular states diversify when clonal three-dimensional (3D) spheroids are cultured in basement membrane, and one such state is associated with stress tolerance and poor response to anticancer therapeutics. Here, we found that this state was jointly coordinated by the NRF2 and p53 pathways, which were costabilized by spontaneous oxidative stress within 3D cultures. Inhibition of NRF2 or p53 individually disrupted some of the transcripts defining the regulatory state but did not yield a notable phenotype in nontransformed breast epithelial cells. In contrast, combined perturbation prevented 3D growth in an oxidative stress–dependent manner. By integrating systems models of NRF2 and p53 signaling in a single oxidative stress network, we recapitulated these observations and made predictions about oxidative stress profiles during 3D growth. NRF2 and p53 signaling were similarly coordinated in normal breast epithelial tissue and hormone-negative ductal carcinoma in situ lesions but were uncoupled in triple-negative breast cancer (TNBC), a subtype in which p53 is usually mutated. Using the integrated model, we correlated the extent of this uncoupling in TNBC cell lines with the importance of NRF2 in the 3D growth of these cell lines and their predicted handling of oxidative stress. Our results point to an oxidative stress tolerance network that is important for single cells during glandular development and the early stages of breast cancer.more » « less
-
One of the fundamental questions in developmental biology is how a cell is specified to differentiate as a specialized cell type. Traditionally, plant cell types were defined based on their function, location, morphology, and lineage. Currently, in the age of single-cell biology, researchers typically attempt to assign plant cells to cell types by clustering them based on their transcriptomes. However, because cells are dynamic entities that progress through the cell cycle and respond to signals, the transcriptome also reflects the state of the cell at a particular moment in time, raising questions about how to define a cell type. We suggest that these complexities and dynamics of cell states are of interest and further consider the roles signaling, stochasticity, cell cycle, and mechanical forces play in plant cell fate specification. Once established, cell identity must also be maintained. With the wealth of single-cell data coming out, the field is poised to elucidate both the complexity and dynamics of cell states.more » « less
-
Abstract Advances in single-cell technologies allow scrutinizing of heterogeneous cell states, however, detecting cell-state transitions from snap-shot single-cell transcriptome data remains challenging. To investigate cells with transient properties or mixed identities, we present MuTrans, a method based on multiscale reduction technique to identify the underlying stochastic dynamics that prescribes cell-fate transitions. By iteratively unifying transition dynamics across multiple scales, MuTrans constructs the cell-fate dynamical manifold that depicts progression of cell-state transitions, and distinguishes stable and transition cells. In addition, MuTrans quantifies the likelihood of all possible transition trajectories between cell states using coarse-grained transition path theory. Downstream analysis identifies distinct genes that mark the transient states or drive the transitions. The method is consistent with the well-established Langevin equation and transition rate theory. Applying MuTrans to datasets collected from five different single-cell experimental platforms, we show its capability and scalability to robustly unravel complex cell fate dynamics induced by transition cells in systems such as tumor EMT, iPSC differentiation and blood cell differentiation. Overall, our method bridges data-driven and model-based approaches on cell-fate transitions at single-cell resolution.more » « less
-
Hematopoiesis has been studied via stem cell labeling using barcodes, viral integration sites (VISs), or in situ methods. Subsequent proliferation and differentiation preserve the tag identity, thus defining a clone of mature cells across multiple cell type or lineages. By tracking the population of clones, measured within samples taken at discrete time points, we infer physiological parameters associated with a hybrid stochastic-deterministic mathematical model of hematopoiesis. We analyze clone population data from Koelle et al. ( Koelle et al., 2017 ) and compare the states of clones (mean and variance of their abundances) and the state-space density of clones with the corresponding quantities predicted from our model. Comparing our model to the tagged granulocyte populations, we find parameters (stem cell carrying capacity, stem cell differentiation rates, and the proliferative potential of progenitor cells, and sample sizes) that provide reasonable fits in three out of four animals. Even though some observed features cannot be quantitatively reproduced by our model, our analyses provides insight into how model parameters influence the underlying mechanisms in hematopoiesis. We discuss additional mechanisms not incorporated in our model.more » « less
An official website of the United States government

