scholarly journals Single cell transcriptomics reveals molecular subtype and functional heterogeneity in models of breast cancer

2018 ◽  
Author(s):  
Daniel L. Roden ◽  
Laura A. Baker ◽  
Benjamin Elsworth ◽  
Chia-Ling Chan ◽  
Kate Harvey ◽  
...  

AbstractBreast cancer has long been classified into a number of molecular subtypes that predict prognosis and therefore influence clinical treatment decisions. Cellular heterogeneity is also evident in breast cancers and plays a key role in the development, evolution and metastatic progression of many cancers. How clinical heterogeneity relates to cellular heterogeneity is poorly understood, so we approached this question using single cell gene expression analysis of well established in vitro and in vivo models of disease.To explore the cellular heterogeneity in breast cancer we first examined a panel of genes that define the PAM50 classifier of molecular subtype. Five breast cancer cell line models (MCF7, BT474, SKBR3, MDA-MB-231, and MDA-MB-468) were selected as representatives of the intrinsic molecular subtypes (luminal A and B, basal-like, and Her2-enriched). Single cell multiplex RT-PCR was used to isolate and quantify the gene expression of single cells from each of these models, and the PAM50 classifier applied. Using this approach, we identified heterogeneity of intrinsic subtypes at single-cell level, indicating that cells with different subtypes exist within a cell line. Using the Chromium 10X system, this study was extended into thousands of cells from the MCF7 cell-line and an ER+ patient derived xenograft (PDX) model and again identified significant intra-tumour heterogeneity of molecular subtype.Estrogen Receptor (ER) is an important driver and therapeutic target in many breast cancers. It is heterogeneously expressed in a proportion of clinical cases but the significance of this to ER activity is unknown. Significant heterogeneity in the transcriptional activation of ER regulated genes was observed within tumours. This differential activation of the ER cistrome aligned with expression of two known transcriptional co-regulatory factors of ER (FOXA1 and PGR).To examine the degree of heterogeneity for other important phenotypic traits, we used an unsupervised clustering approach to identify cellular sub-populations with diverse cancer associated transcriptional properties, such as: proliferation; hypoxia; and treatment resistance. In particular, we show that we can identify two distinct sub-populations of cells that may have denovo resistance to endocrine therapies in a treatment naïve PDX model of ER+ breast cancer. One of these consists of cells with a non-proliferative transcriptional phenotype that is enriched for transcriptional properties of ERBB2 tumours. The other is heavily enriched for components of the primary cilia. Gene regulatory networks were used to identify transcription factor regulons that are active in each cell, leading us to identify potential transcriptional drivers (such as E2F7, MYB and RFX3) of the cilia associated endocrine resistant cells. This rare subpopulation of cells also has a highly heterogenous mix of intrinsic subtypes highlighting a potential role of intra-tumour subtype heterogeneity in endocrine resistance and metastatic potential.Overall, These results suggest a high degree of cellular heterogeneity within breast cancer models, even cell lines, that can be functionally dissected into sub-populations of cells with transcriptional phenotypes of potential clinical relevance.

2020 ◽  
pp. 848-859
Author(s):  
Max Klebe ◽  
Carlo Fremd ◽  
Mark Kriegsmann ◽  
Katharina Kriegsmann ◽  
Thomas Albrecht ◽  
...  

PURPOSE Conversion of tumor subtype frequently occurs in the course of metastatic breast cancer but is a poorly understood phenomenon. This study aims to compare molecular subtypes with subsequent lung or pleural metastasis. PATIENTS AND METHODS In a cohort of 57 patients with breast cancer and lung or pleural metastasis (BCLPM), we investigated paired primary and metastatic tissues for differential gene expression of 269 breast cancer genes. The PAM50 classifier was applied to identify intrinsic subtypes, and differential gene expression and cluster analysis were used to further characterize subtypes and tumors with subtype conversion. RESULTS In primary breast cancer, the most frequent molecular subtype was luminal A (lumA; 49.1%); it was luminal B (lumB) in BCLPM (38.6%). Subtype conversion occurred predominantly in lumA breast cancers compared with other molecular subtypes (57.1% v 27.6%). In lumA cancers, 62 genes were identified with differential expression in metastatic versus primary disease, compared with only 10 differentially expressed genes in lumB, human epidermal growth factor receptor 2 (HER2)–enriched, and basal subtypes combined. Gene expression changes in lumA cancers affected not only the repression of the estrogen receptor pathway and cell cycle–related genes but also the WNT pathway, proteinases ( MME, MMP11), and motility-associated cytoskeletal proteins (CK5, CK14, CK17). Subtype-switched lumA cancers were further characterized by cell proliferation and cell cycle checkpoint gene upregulation and dysregulation of the p53 pathway. This involved 83 notable gene expression changes. CONCLUSION Our results indicate that gene expression changes and subsequent subtype conversion occur on a large scale in metastatic luminal A–type breast cancer compared with other molecular subtypes. This underlines the significance of molecular changes in metastatic disease, especially in tumors of initially low aggressive potential.


2015 ◽  
Vol 33 (7_suppl) ◽  
pp. 305-305 ◽  
Author(s):  
Jordan Kardos ◽  
Jonathan J Melquist ◽  
David D. Chism ◽  
Woonyoung Choi ◽  
Katherine Cockerill ◽  
...  

305 Background: African American (AA) patients with urothelial carcinoma (UC) have been known to have a worse prognosis even when corrected for variables such as tumor stage and grade. Analysis of gene expression of several malignancies has resulted in the discovery of molecular subtypes with well-defined intrinsic biology. Recent studies in high grade (HG), muscle-invasive UC have led to the identification of two intrinsic, molecular subsets termed “luminal” and “basal” with characteristics of stages of urothelial differentiation, and that remarkably reflect the luminal and basal-like molecular subtypes of breast cancer. Patients with basal-like UC have a significantly worse overall survival. Methods: A total of 215 HG muscle-invasive UC tumors from the MDACC (n=75) and TCGA (n=140) were used to make intrinsic subtype calls using gene expression profiling (MDACC: DASL [cDNA-mediated Annealing, Selection, extension, and Ligation] and TCGA: RNA seq). Basal and luminal subtype calls were derived using previously published subtype classifiers (Damrauer et. al. PNAS, 2014 and Choi et. al. Cancer Cell, 2014). Patients were classified into AA and non-AA (white, Hispanic, or Asian) based upon self-reported race. Results: In total there were 16 and 199 tumors from AA and non-AA patients respectively. In non-AA patients, the proportion of tumors that were classified as basal and luminal were approximately equal (93 and 106 respectively), while in AA patients, there was enrichment of basal tumors (12 basal and 4 luminal) (p=0.03735, Fisher’s exact test). Conclusions: AA patients are enriched in the basal molecular subtype of UC. Similar findings have been previously documented in AA women with breast cancer. The enrichment of basal UC in AAs suggests that a biological explanation may in part underlie the poor outcomes seen in AA patients. Future studies will explore the prognostic and predictive implications of basal subtype in AA patients with UC.


2013 ◽  
Vol 20 (3) ◽  
pp. 339-348 ◽  
Author(s):  
Sewha Kim ◽  
Do Hee Kim ◽  
Woo-Hee Jung ◽  
Ja Seung Koo

The aim of this study was to investigate the expression of glutamine metabolism-related proteins to determine whether glutamine is metabolized differently according to breast cancer molecular subtype. We generated a tissue microarray of 702 breast cancer patients and performed immunohistochemical staining for glutamine metabolism-related proteins, including glutaminase 1 (GLS1 (GLS)), glutamate dehydrogenase (GDH (H6PD)), and amino acid transporter-2 (ASCT2 (SLC1A5)), which were separately evaluated in tumor and stroma compartments and then analyzed by breast cancer molecular subtypes. Breast cancers were classified as follows: 293 luminal A (41.7%), 166 luminal B (23.6%), 67 HER2 type (9.6%), and 176 TNBC (25.1%). HER2 type showed the highest stromal GLS1 (P=0.001), tumoral GDH (P=0.001), stromal GDH (P<0.001), and tumoral ASCT (P<0.001) expression. We identified differential expression of glutamine metabolism-related proteins according to molecular subtype of breast cancer. The highest glutamine metabolic activity was seen in HER2-type breast cancer.


Author(s):  
Player A ◽  
◽  
Abraham N ◽  
Abdulrahman N ◽  
Nsende E ◽  
...  

Purpose: The MYBL1 gene is a strong transcriptional activator, associated with cell cycle signaling and differentiation. Data show the gene is overexpressed in triple negative breast cancers. Considering the possibility that MYBL1 might be involved in events associated with the pathogenesis of these cancers, we sought to identify genes associated with MYBL1 expression in triple negative breast cancer. Methods: shRNA lentiviral knockdown was used to down-regulate the MYBL1 gene. Microarray analyses were used to identify genes either directly or indirectly affected by targeting MYBL1 knockdown. Data analyses was performed utilizing Affymetrix TAC 4.0, Chip X transcription factor analyses, Target Scan miRNA analyses, and STRING analyses was used to determine protein: protein interaction and pathway analyses. Web Gestalt and Gene Ontology were used to determine pathway and gene-set enrichments. Publicly available patient and cell line datasets were retrieved and processed using resources available in Gene Expression Omnibus and Oncomine. The polymerase chain reaction and western analyses were used to determine transcript and protein levels, respectively. Results: Knockdown of MYBL1 in a triple negative breast cell line led to down-regulation of MYBL2, TCF19, KIF18b along with an enrichment of cell cycle signaling genes. Gene expression analyses show that MYBL1, MYBL2, TCF19 and KIF18b display a similar pattern of expression in breast cell lines and many of the archival patient datasets examined. Conclusion: TNBC is a heterogeneous subtype, so these data suggest that cancers that over-express MYBL1, express MYBL2, TCF19 and KIF18b. Bioinformatic analyses suggest MYBL1 regulates MYBL2 which leads to regulation of TCF19 and KIF18b.


Animals ◽  
2021 ◽  
Vol 11 (2) ◽  
pp. 502
Author(s):  
Filipe Almeida ◽  
Andreia Gameiro ◽  
Jorge Correia ◽  
Fernando Ferreira

Feline mammary carcinoma (FMC) is the third most common type of neoplasia in cats, sharing similar epidemiological features with human breast cancer. In humans, histone deacetylases (HDACs) play an important role in the regulation of gene expression, with HDAC inhibitors (HDACis) disrupting gene expression and leading to cell death. In parallel, microtubules inhibitors (MTIs) interfere with the polymerization of microtubules, leading to cell cycle arrest and apoptosis. Although HDACis and MTIs are used in human cancer patients, in cats, data is scarce. In this study, we evaluated the antitumor properties of six HDACis (CI-994, panobinostat, SAHA, SBHA, scriptaid, and trichostatin A) and four MTIs (colchicine, nocodazole, paclitaxel, and vinblastine) using three FMC cell lines (CAT-MT, FMCp, and FMCm), and compared with the human breast cancer cell line (SK-BR-3). HDACis and MTIs exhibited dose-dependent antitumor effects in FMC cell lines, and for all inhibitors, the IC50 values were determined, with one feline cell line showing reduced susceptibility (FMCm). Immunoblot analysis confirmed an increase in the acetylation status of core histone protein HDAC3 and flow cytometry showed that HDACis and MTIs lead to cellular apoptosis. Overall, our study uncovers HDACis and MTIs as promising anti-cancer agents to treat FMCs.


2010 ◽  
Vol 28 (7) ◽  
pp. 1161-1167 ◽  
Author(s):  
Anita K. Dunbier ◽  
Helen Anderson ◽  
Zara Ghazoui ◽  
Elizabeth J. Folkerd ◽  
Roger A'Hern ◽  
...  

Purpose To determine whether plasma estradiol (E2) levels are related to gene expression in estrogen receptor (ER)–positive breast cancers in postmenopausal women. Materials and Methods Genome-wide RNA profiles were obtained from pretreatment core-cut tumor biopsies from 104 postmenopausal patients with primary ER-positive breast cancer treated with neoadjuvant anastrozole. Pretreatment plasma E2 levels were determined by highly sensitive radioimmunoassay. Genes were identified for which expression was correlated with pretreatment plasma E2 levels. Validation was performed in an independent set of 73 ER-positive breast cancers. Results The expression of many known estrogen-responsive genes and gene sets was highly significantly associated with plasma E2 levels (eg, TFF1/pS2, GREB1, PDZK1 and PGR; P < .005). Plasma E2 explained 27% of the average expression of these four average estrogen-responsive genes (ie, AvERG; r = 0.51; P < .0001), and a standardized mean of plasma E2 levels and ER transcript levels explained 37% (r, 0.61). These observations were validated in an independent set of 73 ER-positive tumors. Exploratory analysis suggested that addition of the nuclear coregulators in a multivariable analysis with ER and E2 levels might additionally improve the relationship with the AvERG. Plasma E2 and the standardized mean of E2 and ER were both significantly correlated with 2-week Ki67, a surrogate marker of clinical outcome (r = −0.179; P = .05; and r = −0.389; P = .0005, respectively). Conclusion Plasma E2 levels are significantly associated with gene expression of ER-positive breast cancers and should be considered in future genomic studies of ER-positive breast cancer. The AvERG is a new experimental tool for the study of putative estrogenic stimuli of breast cancer.


Sign in / Sign up

Export Citation Format

Share Document