Dataset Information

Transcriptome analysis of heterogeneity in mouse model of metastatic breast cancer.

ABSTRACT:

Background

Cancer metastasis is a complex process involving the spread of malignant cells from a primary tumor to distal organs. Understanding this cascade at a mechanistic level could provide critical new insights into the disease and potentially reveal new avenues for treatment. Transcriptome profiling of spontaneous cancer models is an attractive method to examine the dynamic changes accompanying tumor cell spread. However, such studies are complicated by the underlying heterogeneity of the cell types involved. The purpose of this study was to examine the transcriptomes of metastatic breast cancer cells using the well-established MMTV-PyMT mouse model.

Methods

Organ-derived metastatic cell lines were harvested from 10 female MMTV-PyMT mice. Cancer cells were isolated and sorted based on the expression of CD44^low/EpCAM^high or CD44^high/EpCAM^high surface markers. RNA from each cell line was extracted and sequenced using the NextSeq 500 Illumina platform. Tissue-specific genes were compared across the different metastatic and primary tumor samples. Reads were mapped to the mouse genome using STAR, and gene expression was quantified using RSEM. Single-cell RNA-seq (scRNA-seq) was performed on select samples using the ddSeq platform by BioRad and analyzed using Seurat v3.2.3. Monocle2 was used to infer pseudo-time progression.

Results

Comparison of RNA sequencing data across all cell populations produced distinct gene clusters. Differential gene expression patterns related to CD44 expression, organ tropism, and immunomodulatory signatures were observed. scRNA-seq identified expression profiles based on tissue-dependent niches and clonal heterogeneity. These cohorts of data were narrowed down to identify subsets of genes with high expression and known metastatic propensity. Dot plot analyses further revealed clusters expressing cancer stem cell and cancer dormancy markers. Changes in relevant genes were investigated across pseudo-time and tissue origin using Monocle2. These data revealed transcriptomes that may contribute to sub-clonal evolution and treatment evasion during cancer progression.

Conclusions

We performed a comprehensive transcriptome analysis of tumor heterogeneity and organ tropism during breast cancer metastasis. These data add to our understanding of metastatic progression and highlight targets for breast cancer treatment. These markers could also be used to image the impact of tumor heterogeneity on metastases.

SUBMITTER: Ionkina AA

PROVIDER: S-EPMC8477508 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Breast cancer (BC) is a heterogeneous disease, which is primarily classified according to hormone receptors and HER2 expression. Despite the many advances in BC diagnosis and management, the identification of novel actionable therapeutic targets expressed by cancerous cells has always been a daunting task due to the large heterogeneity of the disease and the presence of non-cancerous cells (i.e., immune cells and stromal cells) within the tumor microenvironment. In the current study, we employed computational algorithms to decipher the cellular composition of estrogen receptor-positive (ER+), HER2+, ER+HER2+, and triple-negative BC (TNBC) subtypes from a total of 49,899 single cells' publicly available transcriptomic data derived from 26 BC patients. Restricting the analysis to EPCAM+Lin- tumor epithelial cells, we identified the enriched gene sets in each BC molecular subtype. Integration of single-cell transcriptomic with CRISPR-Cas9 functional screen data identified 13 potential therapeutic targets for ER+, 44 potential therapeutic targets for HER2+, and 29 potential therapeutic targets for TNBC. Interestingly, several of the identified therapeutic targets outperformed the current standard of care for each BC subtype. Given the aggressive nature and lack of targeted therapies for TNBC, elevated expression of ENO1, FDPS, CCT6A, TUBB2A, and PGK1 predicted worse relapse-free survival (RFS) in basal BC (n = 442), while elevated expression of ENO1, FDPS, CCT6A, and PGK1 was observed in the most aggressive BLIS TNBC subtype. Mechanistically, targeted depletion of ENO1 and FDPS halted TNBC cell proliferation, colony formation, and organoid tumor growth under 3-dimensional conditions and increased cell death, suggesting their potential use as novel therapeutic targets for TNBC. Differential expression and gene set enrichment analysis in TNBC revealed enrichment in the cycle and mitosis functional categories in FDPShigh, while ENO1high was associated with numerous functional categories, including cell cycle, glycolysis, and ATP metabolic processes. Taken together, our data are the first to unravel the unique gene signatures and to identify novel dependencies and therapeutic vulnerabilities for each BC molecular subtype, thus setting the foundation for the future development of more effective targeted therapies for BC.