Hybrid Stem Cell States: Insights Into the Relationship Between Mammary Development and Breast Cancer Using Single-Cell Transcriptomics

High-Throughput Transcriptomics of Nontumorigenic Breast Cells Exposed to Environmentally Relevant Chemicals

BACKGROUND

There is a suite of chemicals, including metals, pesticides, and personal care product compounds, which are commonly detected at high levels in US Center for Disease Control's National Health and Nutrition Examination Survey (NHANES) chemical biomarker screens. Whether these chemicals influence development of breast cancer is not well understood.

OBJECTIVES

The objectives were to perform an unbiased concentration-dependent assessment of these chemicals, to quantify differences in cancer-specific genes and pathways, to describe if these differences occur at human population-relevant concentrations, and to specifically test for differences in markers of stemness and cellular plasticity.

METHODS

We treated nontumorigenic mammary epithelial cells, MCF10A, with 21 chemicals at four concentrations (25  nM , 250  nM , 2.5 μ M , and 25 μ M ) for 48 h. We conducted RNA-sequencing for these 408 samples, adapting the plexWell plate-based RNA-sequencing method to analyze differences in gene expression. We calculated gene and biological pathway-specific benchmark concentrations (BMCs) using BMDExpress3, identifying differentially expressed genes and generating the best fit benchmark concentration models for each chemical across all genes. We identified enriched biological processes and pathways for each chemical and tested whether chemical exposures change predicted cell type distributions. We contextualized benchmark concentrations relative to human population biomarker concentrations in NHANES.

RESULTS

We detected chemical concentration-dependent differences in gene expression for thousands of genes. Enrichment and cell type distribution analyses showed benchmark concentration responses correlated with differences in breast cancer-related pathways, including induction of basal-like characteristics for some chemicals, including arsenic, lead, copper, and methyl paraben. Comparison of benchmark data to NHANES chemical biomarker (urine or blood) concentrations indicated an overlap between exposure levels and levels sufficient to cause a gene expression response.

DISCUSSION

These analyses revealed that many of these 21 chemicals resulted in differences in genes and pathways involved in breast cancer in vitro at human exposure-relevant concentrations. https://doi.org/10.1289/EHP12886.

Redox signalling regulates breast cancer metastasis via phenotypic and metabolic reprogramming due to p63 activation by HIF1α

BACKGROUND

Redox signaling caused by knockdown (KD) of Glutathione Peroxidase 2 (GPx2) in the PyMT mammary tumour model promotes metastasis via phenotypic and metabolic reprogramming. However, the tumour cell subpopulations and transcriptional regulators governing these processes remained unknown.

METHODS

We used single-cell transcriptomics to decipher the tumour cell subpopulations stimulated by GPx2 KD in the PyMT mammary tumour and paired pulmonary metastases. We analyzed the EMT spectrum across the various tumour cell clusters using pseudotime trajectory analysis and elucidated the transcriptional and metabolic regulation of the hybrid EMT state.

RESULTS

Integration of single-cell transcriptomics between the PyMT/GPx2 KD primary tumour and paired lung metastases unraveled a basal/mesenchymal-like cluster and several luminal-like clusters spanning an EMT spectrum. Interestingly, the luminal clusters at the primary tumour gained mesenchymal gene expression, resulting in epithelial/mesenchymal subpopulations fueled by oxidative phosphorylation (OXPHOS) and glycolysis. By contrast, at distant metastasis, the basal/mesenchymal-like cluster gained luminal and mesenchymal gene expression, resulting in a hybrid subpopulation using OXPHOS, supporting adaptive plasticity. Furthermore, p63 was dramatically upregulated in all hybrid clusters, implying a role in regulating partial EMT and MET at primary and distant sites, respectively. Importantly, these effects were reversed by HIF1α loss or GPx2 gain of function, resulting in metastasis suppression.

CONCLUSIONS

Collectively, these results underscored a dramatic effect of redox signaling on p63 activation by HIF1α, underlying phenotypic and metabolic plasticity leading to mammary tumour metastasis.

TLE3 Sustains Luminal Breast Cancer Lineage Fidelity to Suppress Metastasis

Breast cancer subtypes and their phenotypes parallel different stages of the mammary epithelial cell developmental hierarchy. Discovering mechanisms that control lineage identity could provide novel avenues for mitigating disease progression. Here we report that the transcriptional corepressor TLE3 is a guardian of luminal cell fate in breast cancer and operates independently of the estrogen receptor. In luminal breast cancer, TLE3 actively repressed the gene-expression signature associated with highly aggressive basal-like breast cancers (BLBC). Moreover, maintenance of the luminal lineage depended on the appropriate localization of TLE3 to its transcriptional targets, a process mediated by interactions with FOXA1. By repressing genes that drive BLBC phenotypes, including SOX9 and TGFβ2, TLE3 prevented the acquisition of a hybrid epithelial-mesenchymal state and reduced metastatic capacity and aggressive cellular behaviors. These results establish TLE3 as an essential transcriptional repressor that sustains the more differentiated and less metastatic nature of luminal breast cancers. Approaches to induce TLE3 expression could promote the acquisition of less aggressive, more treatable disease states to extend patient survival.

SIGNIFICANCE

Transcriptional corepressor TLE3 actively suppresses SOX9 and TGFβ transcriptional programs to sustain the luminal lineage identity of breast cancer cells and to inhibit metastatic progression.

Single-cell transcriptome sequencing reveals potential novel combination of biomarkers for antibody-based cancer therapeutics in hepatocellular carcinoma

Background: Antibody-based cancer therapeutics is developing rapidly in recent years for its advantages in precisely targeting the tumor cells. However, tumor-specific cell surface antigens are still lacking, and the heterogeneity of tumor mass greatly impeded the development of effective drugs.

Methods: In the present study, single-cell RNA sequencing was used to dissect tumor heterogeneity in human hepatocellular carcinoma (HCC). Tissues from different spatial regions including the tumor, para-tumor, and distant normal liver tissues were dissociated into single cells, and the gene expressions were compared in a different subpopulation of cells from these regions and validated in independent cohorts.

Results: A total of 28 cell clusters with different distribution patterns and gene expression profiles were identified within a heterogenous tumor and its paired liver tissues. Differentially expressed genes encoding the plasma membrane in cells with hepatic lineage were further extracted from single-cell transcriptome sequencing and validated in TCGA database. A 3-gene signature was identified to be significantly upregulated in dominant HCC tumor cell subpopulations with prognostic significance and validated in multiple independent patient cohorts.

Conclusion: The composition of the three plasma membrane proteins on the surface of HCC tumor cells within a heterogenous tumor might indicate poor prognostic tumor subpopulations during cancer evolution and potential therapeutic targets.

Embryonic Programs in Cancer and Metastasis—Insights From the Mammary Gland

Cancer is characterized as a reversion of a differentiated cell to a primitive cell state that recapitulates, in many aspects, features of embryonic cells. This review explores the current knowledge of developmental mechanisms that are essential for embryonic mouse mammary gland development, with a particular focus on genes and signaling pathway components that are essential for the induction, morphogenesis, and lineage specification of the mammary gland. The roles of these same genes and signaling pathways in mammary gland or breast tumorigenesis and metastasis are then summarized. Strikingly, key embryonic developmental pathways are often reactivated or dysregulated during tumorigenesis and metastasis in processes such as aberrant proliferation, epithelial-to-mesenchymal transition (EMT), and stem cell potency which affects cellular lineage hierarchy. These observations are in line with findings from recent studies using lineage tracing as well as bulk- and single-cell transcriptomics that have uncovered features of embryonic cells in cancer and metastasis through the identification of cell types, cell states and characterisation of their dynamic changes. Given the many overlapping features and similarities of the molecular signatures of normal development and cancer, embryonic molecular signatures could be useful prognostic markers for cancer. In this way, the study of embryonic development will continue to complement the understanding of the mechanisms of cancer and aid in the discovery of novel therapeutic targets and strategies.

Transcriptional changes in the mammary gland during lactation revealed by single cell sequencing of cells from human milk

Under normal conditions, the most significant expansion and differentiation of the adult mammary gland occurs in response to systemic reproductive hormones during pregnancy and lactation to enable milk synthesis and secretion to sustain the offspring. However, human mammary tissue remodelling that takes place during pregnancy and lactation remains poorly understood due to the challenge of acquiring samples. We report here single-cell transcriptomic analysis of 110,744 viable breast cells isolated from human milk or non-lactating breast tissue, isolated from nine and seven donors, respectively. We found that human milk largely contains epithelial cells belonging to the luminal lineage and a repertoire of immune cells. Further transcriptomic analysis of the milk cells identified two distinct secretory cell types that shared similarities with luminal progenitors, but no populations comparable to hormone-responsive cells. Taken together, our data offers a reference map and a window into the cellular dynamics that occur during human lactation and may provide further insights on the interplay between pregnancy, lactation and breast cancer.

Human mammary tissue remodelling that takes place during pregnancy and lactation remains poorly understood. Here the authors characterize cells in human milk, identifying epithelial cells resembling luminal progenitors and immune cells, contributing insights into this process.

Application of Single-Cell Multi-Omics in Dissecting Cancer Cell Plasticity and Tumor Heterogeneity

Tumor heterogeneity, a hallmark of cancer, impairs the efficacy of cancer therapy and drives tumor progression. Exploring inter- and intra-tumoral heterogeneity not only provides insights into tumor development and progression, but also guides the design of personalized therapies. Previously, high-throughput sequencing techniques have been used to investigate the heterogeneity of tumor ecosystems. However, they could not provide a high-resolution landscape of cellular components in tumor ecosystem. Recently, advance in single-cell technologies has provided an unprecedented resolution to uncover the intra-tumoral heterogeneity by profiling the transcriptomes, genomes, proteomes and epigenomes of the cellular components and also their spatial distribution, which greatly accelerated the process of basic and translational cancer research. Importantly, it has been demonstrated that some cancer cells are able to transit between different states in order to adapt to the changing tumor microenvironment, which led to increased cellular plasticity and tumor heterogeneity. Understanding the molecular mechanisms driving cancer cell plasticity is critical for developing precision therapies. In this review, we summarize the recent progress in dissecting the cancer cell plasticity and tumor heterogeneity by use of single-cell multi-omics techniques.

Single-cell evaluation reveals shifts in the tumor-immune niches that shape and maintain aggressive lesions in the breast

There is an unmet clinical need for stratification of breast lesions as indolent or aggressive to tailor treatment. Here, single-cell transcriptomics and multiparametric imaging applied to a mouse model of breast cancer reveals that the aggressive tumor niche is characterized by an expanded basal-like population, specialization of tumor subpopulations, and mixed-lineage tumor cells potentially serving as a transition state between luminal and basal phenotypes. Despite vast tumor cell-intrinsic differences, aggressive and indolent tumor cells are functionally indistinguishable once isolated from their local niche, suggesting a role for non-tumor collaborators in determining aggressiveness. Aggressive lesions harbor fewer total but more suppressed-like T cells, and elevated tumor-promoting neutrophils and IL-17 signaling, disruption of which increase tumor latency and reduce the number of aggressive lesions. Our study provides insight into tumor-immune features distinguishing indolent from aggressive lesions, identifies heterogeneous populations comprising these lesions, and supports a role for IL-17 signaling in aggressive progression.

Partial EMT in Squamous Cell Carcinoma: A Snapshot

In the process of cancer EMT, some subgroups of cancer cells simultaneously exhibit both mesenchymal and epithelial characteristics, a phenomenon termed partial EMT (pEMT). pEMT is a plastic state in which cells coexpress epithelial and mesenchymal markers. In squamous cell carcinoma (SCC), pEMT is regulated, and the phenotype is maintained via the HIPPO pathway, NOTCH pathway and TGF-β pathways and by microRNAs, lncRNAs and the cancer microenvironment (CME); thus, SCC exhibits aggressive tumorigenic properties and high stemness, which leads collective migration and therapy resistance. Few studies have reported therapeutic interventions to address cells that have undergone pEMT, and this approach may be an effective way to inhibit the plasticity, drug resistance and metastatic potential of SCC.

Single-cell RNA sequencing in cancer: Applications, advances, and emerging challenges

Cancer has become one of the greatest threats to human health, and new technologies are urgently needed to further clarify the mechanisms of cancer so that better detection and treatment strategies can be developed. At present, extensive genomic analysis and testing of clinical specimens shape the insights into carcinoma. Nevertheless, carcinoma of humans is a complex ecosystem of cells, including carcinoma cells and immunity-related and stroma-related subsets, with accurate characteristics obscured by extensive genome-related approaches. A growing body of research shows that sequencing of single-cell RNA (scRNA-seq) is emerging to be an effective way for dissecting human tumor tissue at single-cell resolution, presenting one prominent way for explaining carcinoma biology. This review summarizes the research progress of scRNA-seq in the field of tumors, focusing on the application of scRNA-seq in tumor circulating cells, tumor stem cells, tumor drug resistance, the tumor microenvironment, and so on, which provides a new perspective for tumor research.

Mammary gland development from a single cell 'omics view

Understanding the complexity and heterogeneity of mammary cell subpopulations is vital to delineate the mechanisms behind breast cancer development, progression and prevention. Increasingly sophisticated tools for investigating these cell subtypes has led to the development of a greater understanding of these cell subtypes, complex interplay of certain subtypes and their developmental potential. Of note, increasing accessibility and affordability of single cell technologies has led to a plethora of studies being published containing data from mammary cell subtypes and their differentiation potential in both mice and human data sets. Here, we review the different types of single cell technologies and how they have been used to improve our understanding of mammary gland development.

Characterization of Gene Expression Signatures for the Identification of Cellular Heterogeneity in the Developing Mammary Gland

The developing mammary gland depends on several transcription-dependent networks to define cellular identities and differentiation trajectories. Recent technological advancements that allow for single-cell profiling of gene expression have provided an initial picture into the epithelial cellular heterogeneity across the diverse stages of gland maturation. Still, a deeper dive into expanded molecular signatures would improve our understanding of the diversity of mammary epithelial and non-epithelial cellular populations across different tissue developmental stages, mouse strains and mammalian species. Here, we combined differential mammary gland fractionation approaches and transcriptional profiles obtained from FACS-isolated mammary cells to improve our definitions of mammary-resident, cellular identities at the single-cell level. Our approach yielded a series of expression signatures that illustrate the heterogeneity of mammary epithelial cells, specifically those of the luminal fate, and uncovered transcriptional changes to their lineage-defined, cellular states that are induced during gland development. Our analysis also provided molecular signatures that identified non-epithelial mammary cells, including adipocytes, fibroblasts and rare immune cells. Lastly, we extended our study to elucidate expression signatures of human, breast-resident cells, a strategy that allowed for the cross-species comparison of mammary epithelial identities. Collectively, our approach improved the existing signatures of normal mammary epithelial cells, as well as elucidated the diversity of non-epithelial cells in murine and human breast tissue. Our study provides a useful resource for future studies that use single-cell molecular profiling strategies to understand normal and malignant breast development.

Hybrid E/M Phenotype(s) and Stemness: A Mechanistic Connection Embedded in Network Topology

Metastasis remains an unsolved clinical challenge. Two crucial features of metastasizing cancer cells are (a) their ability to dynamically move along the epithelial-hybrid-mesenchymal spectrum and (b) their tumor initiation potential or stemness. With increasing functional characterization of hybrid epithelial/mesenchymal (E/M) phenotypes along the spectrum, recent in vitro and in vivo studies have suggested an increasing association of hybrid E/M phenotypes with stemness. However, the mechanistic underpinnings enabling this association remain unclear. Here, we develop a mechanism-based mathematical modeling framework that interrogates the emergent nonlinear dynamics of the coupled network modules regulating E/M plasticity (miR-200/ZEB) and stemness (LIN28/let-7). Simulating the dynamics of this coupled network across a large ensemble of parameter sets, we observe that hybrid E/M phenotype(s) are more likely to acquire stemness relative to "pure" epithelial or mesenchymal states. We also integrate multiple "phenotypic stability factors" (PSFs) that have been shown to stabilize hybrid E/M phenotypes both in silico and in vitro-such as OVOL1/2, GRHL2, and NRF2-with this network, and demonstrate that the enrichment of hybrid E/M phenotype(s) with stemness is largely conserved in the presence of these PSFs. Thus, our results offer mechanistic insights into recent experimental observations of hybrid E/M phenotype(s) that are essential for tumor initiation and highlight how this feature is embedded in the underlying topology of interconnected EMT (Epithelial-Mesenchymal Transition) and stemness networks.

Low-Level Mouse DNA in Conditioned Medium Generates False Positive Cross-Species Contamination Results in Human Organoid Cultures

Cell line authentication is critical for preventing the use of mixed or misidentified cell lines in research. Current efforts include short tandem repeat (STR) analysis and PCR-based assays to detect mixed species cultures. Using PCR analysis with mouse-specific primers, we identified contaminating mouse DNA in growth factor conditioned medium (CM) derived from the L-WRN cell line (L-WRN CM), as well as in human organoid cultures maintained in the L-WRN CM. DNA isolated from L-WRN CM matched the L-WRN cell signature by STR analysis. Organoid lines that were positive for murine DNA by PCR were further analyzed via bulk RNA-sequencing and transcripts were aligned to the human and mouse genomes. RNA analysis failed to detect mouse-specific gene expression above background levels, suggesting no viable murine cells were present in the organoid cultures. We interpret our data to show conclusive evidence that mouse cell-derived CM can be a source of contaminating murine DNA detected in human organoid cultures, even though live, transcriptionally-active murine cells are not present. Together, our findings suggest that multiple methods may be required to authenticate human organoid or cell lines and urges cautious interpretation of DNA-based PCR cell line authentication results.

Hypoxia, partial EMT and collective migration: Emerging culprits in metastasis

Epithelial-mesenchymal transition (EMT) is a cellular biological process involved in migration of primary cancer cells to secondary sites facilitating metastasis. Besides, EMT also confers properties such as stemness, drug resistance and immune evasion which can aid a successful colonization at the distant site. EMT is not a binary process; recent evidence suggests that cells in partial EMT or hybrid E/M phenotype(s) can have enhanced stemness and drug resistance as compared to those undergoing a complete EMT. Moreover, partial EMT enables collective migration of cells as clusters of circulating tumor cells or emboli, further endorsing that cells in hybrid E/M phenotypes may be the 'fittest' for metastasis. Here, we review mechanisms and implications of hybrid E/M phenotypes, including their reported association with hypoxia. Hypoxia-driven activation of HIF-1α can drive EMT. In addition, cyclic hypoxia, as compared to acute or chronic hypoxia, shows the highest levels of active HIF-1α and can augment cancer aggressiveness to a greater extent, including enriching for a partial EMT phenotype. We also discuss how metastasis is influenced by hypoxia, partial EMT and collective cell migration, and call for a better understanding of interconnections among these mechanisms. We discuss the known regulators of hypoxia, hybrid EMT and collective cell migration and highlight the gaps which needs to be filled for connecting these three axes which will increase our understanding of dynamics of metastasis and help control it more effectively.