Metabolomics at the tumor microenvironment interface: Decoding cellular conversations.

Cancer heterogeneity remains a significant challenge for effective cancer treatments. Altered energetics is one of the hallmarks of cancer and influences tumor growth and drug resistance. Studies have shown that heterogeneity exists within the metabolic profile of tumors, and personalized-combination therapy with relevant metabolic interventions could improve patient response. Metabolomic studies are identifying novel biomarkers and therapeutic targets that have improved treatment response. The spatial location of elements in the tumor microenvironment are becoming increasingly important for understanding disease progression. The evolution of spatial metabolomics analysis now allows scientists to deeply understand how metabolite distribution contributes to cancer biology. Recently, these techniques have spatially resolved metabolite distribution to a subcellular level. It has been proposed that metabolite mapping could improve patient outcomes by improving precision medicine, enabling earlier diagnosis and intraoperatively identifying tumor margins. This review will discuss how altered metabolic pathways contribute to cancer progression and drug resistance and will explore the current capabilities of spatial metabolomics technologies and how these could be integrated into clinical practice to improve patient outcomes.

A robust platform for integrative spatial multi-omics analysis to map immune responses to SARS-CoV-2 infection in lung tissues.

The SARS-CoV-2 (COVID-19) virus has caused a devastating global pandemic of respiratory illness. To understand viral pathogenesis, methods are available for studying dissociated cells in blood, nasal samples, bronchoalveolar lavage fluid and similar, but a robust platform for deep tissue characterization of molecular and cellular responses to virus infection in the lungs is still lacking. We developed an innovative spatial multi-omics platform to investigate COVID-19-infected lung tissues. Five tissue-profiling technologies were combined by a novel computational mapping methodology to comprehensively characterize and compare the transcriptome and targeted proteome of virus infected and uninfected tissues. By integrating spatial transcriptomics data (Visium, GeoMx and RNAScope) and proteomics data (CODEX and PhenoImager HT) at different cellular resolutions across lung tissues, we found strong evidence for macrophage infiltration and defined the broader microenvironment surrounding these cells. By comparing infected and uninfected samples, we found an increase in cytokine signalling and interferon responses at different sites in the lung and showed spatial heterogeneity in the expression level of these pathways. These data demonstrate that integrative spatial multi-omics platforms can be broadly applied to gain a deeper understanding of viral effects on cellular environments at the site of infection and to increase our understanding of the impact of SARS-CoV-2 on the lungs.

Overexpression of miRNA-9 enhances galectin-3 levels in oral cavity cancers.

Oral cavity cancer (OCC) is the predominant subtype of head and neck cancer (HNC) and has up to 50% mortality. Genome-wide microRNA (miR) sequencing data indicates overexpression of miR-9-5p in HNC tumours, however, the biological role of miR-9-5p in OCC is complex; it can either act as a tumour suppressor or an oncomir, regulating many target genes at the post-transcriptional level. We have investigated the overexpression of miR-9-5p in three OCC cell lines. We have evaluated its expression levels and Galectin-3 as potential biomarkers in saliva samples collected from controls and OCC patients. We found that over expression of miR-9-5p in OCC cell lines resulted in a significant reduction in cell proliferation and migration, and an increase in apoptosis, which was paralleled by an increase in Galectin-3 secretion and export of Galectin-3 protein. Our data are consistent with miR-9-5p being a modulator of Galectin-3 via the AKT/γ-catenin pathway. In addition, the positive correlation between the levels of miR-9-5p expression and secreted Galectin-3 in saliva reflects a similar relationship in vivo, and supports the utility of their integrative evaluation in OCC. Our findings indicate that both miR-9-5p and Galectin-3 are critical biomolecules in the progression of OCC.

Integrin alpha-2 and beta-1 expression increases through multiple generations of the EDW01 patient-derived xenograft model of breast cancer-insight into their role in epithelial mesenchymal transition in vivo gained from an in vitro model system.

BACKGROUND: Breast cancers acquire aggressive capabilities via epithelial to mesenchymal transition (EMT), in which various integrins/integrin-linked kinase signalling are upregulated. METHODS: We investigated this in two patient-derived xenografts (PDXs) developed from breast-to-bone metastases, and its functional significance in a breast cancer cell line system. ED03 and EDW01 PDXs were grown subcutaneously in immunocompromised SCID mice through 11 passages and 7 passages, respectively. Tumour tissue was assessed using immunohistochemistry (IHC) for oestrogen receptor (ER)-alpha, E-cadherin, vimentin, Twist1, beta-catenin, P120-RasGAP, CD44, CD24 and Ki67, and RT-qPCR of EMT-related factors (CDH1, VIM, CD44, CD24), integrins beta 1 (ITGB1), alpha 2 (ITGA2) and ILK. Integrin and ILK expression in epidermal growth factor (EGF)-induced EMT of the PMC42-ET breast cancer cell line was assessed by RT-qPCR and Western blotting, as were the effects of their transient knockdown via small interfering RNA +/- EGF. Cell migration, changes in cell morphology and adhesion of siRNA-transfected PMC42-ET cells to various extracellular matrix (ECM) substrates was assessed. RESULTS: The ED03 (ER+/PR-/HER2-/lobular) and EDW01 (ER+/PR-/HER2-/ductal) PDXs were both classified as molecular subtype luminal A. ED03 xenografts exhibited mutated E-cadherin with minimal expression, but remained vimentin-negative across all passages. In EDW01, the hypoxic indicator gene CAIX and Twist1 were co-ordinately upregulated at passages 4-5, corresponding with a decrease in E-cadherin. At passages 6-7, VIM was upregulated along with ITGB1 and ITGA2, consistent with an increasing EMT. The ED03 PDX displayed minimal change over passages in mice, for all genes examined. ILK, ITGB1 and ITGA2 mRNAs were also increased in the EGF-induced EMT of PMC42-ET cells (in which CDH1 was downregulated) although siRNA against these targets revealed that this induction was not necessary for the observed EMT. However, their knockdown significantly reduced EMT-associated adhesion and Transwell migration. CONCLUSION: Our data suggest that despite an increase in ITGA2 and ITGB1 gene expression in the EMT exhibited by EDW01 PDX over multiple generations, this pathway may not necessarily drive the EMT process.

High threshold of ß1 integrin inhibition required to block collagen I-induced membrane type-1 matrix metalloproteinase (MT1-MMP) activation of matrix metalloproteinase 2 (MMP-2).

BACKGROUND: Matrix metalloproteinase-2 (MMP-2) is an endopeptidase that facilitates extracellular matrix remodeling and molecular regulation, and is implicated in tumor metastasis. Type I collagen (Col I) regulates the activation of MMP-2 through both transcriptional and post-transcriptional means; however gaps remain in our understanding of the involvement of collagen-binding ß1 integrins in collagen-stimulated MMP-2 activation. METHODS: Three ß1 integrin siRNAs were used to elucidate the involvement of ß1 integrins in the Col I-induced MMP-2 activation mechanism. ß1 integrin knockdown was analyzed by quantitative RT-PCR, Western Blot and FACS analysis. Adhesion assay and collagen gel contraction were used to test the biological effects of ß1 integrin abrogation. MMP-2 activation levels were monitored by gelatin zymography. RESULTS: All three ß1 integrin siRNAs were efficient at ß1 integrin knockdown and FACS analysis revealed commensurate reductions of integrins α2 and α3, which are heterodimeric partners of ß1, but not αV, which is not. All three ß1 integrin siRNAs inhibited adhesion and collagen gel contraction, however only the siRNA showing the greatest magnitude of ß1 knockdown inhibited Col I-induced MMP-2 activation and reduced the accompanying upregulation of MT1-MMP, suggesting a dose response threshold effect. Re-transfection with codon-swapped ß1 integrin overcame the reduction in MMP-2 activation induced by Col-1, confirming the ß1 integrin target specificity. MMP-2 activation induced by TPA or Concanavalin A (Con A) was not inhibited by ß1 integrin siRNA knockdown. CONCLUSION: Together, the data reveals that strong abrogation of ß1 integrin is required to block MMP-2 activation induced by Col I, which may have implications for the therapeutic targeting of ß1 integrin.

The in situ transcriptomic landscape of breast tumour-associated and normal adjacent endothelial cells.

BACKGROUND AND AIMS: Triple Negative Breast Cancer (TNBC) is associated with increased angiogenesis, which is known to aid tumour growth and metastasis. Anti-angiogenic therapies that have been developed to target this feature have mostly generated disappointing clinical results. Further research into targeted approaches is limited by a lack of understanding of the in situ molecular profile of tumour-associated vasculature. In this study, we aimed to understand the differences in the molecular profiles of tumour endothelial cells vs normal-adjacent endothelial cells in TNBC tissues. METHOD: We have applied unbiased whole transcriptome spatial profiling of in situ gene expressions of endothelial cells localized in full-face patient TNBC tissues (n = 4) and normal-adjacent regions of the same patient breast tissues. RESULTS: Our comparative analysis revealed that 2412 genes were differentially expressed (padj < 0.05) between the tumour endothelial cells and normal-adjacent endothelial cells. Pathway enrichment showed the enrichment of gene sets related to cell-cell, cell-ECM adhesion, chromatin organization and remodeling, and protein-DNA complex subunit organization. CONCLUSION: Overall, the results revealed unique molecular profiles and signalling pathways of tumour-associated vasculature, which is a critical step towards larger cohort studies investigating potential targets for TNBC prognosis and anti-angiogenic treatments.

In situ single-cell profiling sheds light on IFI27 localisation during SARS-CoV-2 infection.

The utilization of single-cell resolved spatial transcriptomics to delineate immune responses during SARS-CoV-2 infection was able to identify M1 macrophages to have elevated expression of IFI27 in areas of infection.

Circulating Tumour Cells Indicate the Presence of Residual Disease Post-Castration in Prostate Cancer Patient-Derived Xenograft Models.

Castrate-resistant prostate cancer (CRPC) is the lethal form of prostate cancer. Epithelial mesenchymal plasticity (EMP) has been associated with disease progression to CRPC, and prostate cancer therapies targeting the androgen signalling axis, including androgen deprivation therapy (ADT), promote EMP. We explored effects of castration on EMP in the tumours and circulating tumour cells (CTCs) of patient-derived xenograft (PDX)-bearing castrated mice using human-specific RT-qPCR assays and immunocytochemistry. Expression of prostate epithelial cell marker KLK3 was below detection in most tumours from castrated mice (62%, 23/37 mice), consistent with its known up-regulation by androgens. Endpoint tumour size after castration varied significantly in a PDX model-specific pattern; while most tumours were castration-sensitive (BM18, LuCaP70), the majority of LuCaP105 tumours continued to grow following castration. By contrast, LuCaP96 PDX showed a mixed response to castration. CTCs were detected in 33% of LuCaP105, 43% of BM18, 47% of LuCaP70, and 54% of LuCaP96 castrated mice using RPL32 mRNA measurement in plasma. When present, CTC numbers estimated using human RPL32 expression ranged from 1 to 458 CTCs per ml blood, similar to our previous observations in non-castrated mice. In contrast to their non-castrated counterparts, there was no relationship between tumour size and CTC burden in castrated mice. Unsupervised hierarchical clustering of the gene expression profiles of CTCs collected from castrated and non-castrated mice revealed distinct CTC sub-groups within the pooled population that were classified as having mesenchymal, epithelial, or EMP hybrid gene expression profiles. The epithelial signature was only found in CTCs from non-castrated mice. Hybrid and mesenchymal signatures were detected in CTCs from both castrated and non-castrated mice, with an emphasis towards mesenchymal phenotypes in castrated mice. Post-castration serum PSA levels were either below detection or very low for all the CTC positive samples highlighting the potential usefulness of CTCs for disease monitoring after androgen ablation therapy. In summary, our study of castration effects on prostate cancer PDX CTCs showed that CTCs were often detected in the castrate setting, even in mice with no palpable tumours, and demonstrated the superior ability of CTCs to reveal residual disease over the conventional clinical biomarker serum PSA.

Portable NMR for quantification of breast density in vivo: Proof-of-concept measurements and comparison with quantitative MRI.

Mammographic Density (MD) is the degree of radio-opacity of the breast in an X-ray mammogram. It is determined by the Fibroglandular: Adipose tissue ratio. MD has major implications in breast cancer risk and breast cancer chemoprevention. This study aimed to investigate the feasibility of accurate, low-cost quantification of MD in vivo without ionising radiation. We used single-sided portable nuclear magnetic resonance ("Portable NMR") due to its low cost and the absence of radiation-related safety concerns. Fifteen (N = 15) healthy female volunteers were selected for the study and underwent an imaging routine consisting of 2D X-ray mammography, quantitative breast 3T MRI (Dixon and T1-based 3D compositional breast imaging), and 1D compositional depth profiling of the right breast using Portable NMR. For each participant, all the measurements were made within 3-4 h of each other. MRI-determined tissue water content was used as the MD-equivalent quantity. Portable NMR depth profiles of tissue water were compared with the equivalent depth profiles reconstructed from Dixon and T1-based MR images, which were used as the MD-equivalent reference standard. The agreement between the depth profiles acquired using Portable NMR and the reconstructed reference-standard profiles was variable but overall encouraging. The agreement was somewhat inferior to that seen in breast tissue explant measurements conducted in vitro, where quantitative micro-CT was used as the reference standard. The lower agreement in vivo can be attributed to an uncertainty in the positioning of the Portable NMR sensor on the breast surface and breast compression in Portable NMR measurements. The degree of agreement between Portable NMR and quantitative MRI is encouraging. While the results call for further development of quantitative Portable NMR, they demonstrate the in-principle feasibility of Portable NMR-based quantitative compositional imaging in vivo and show promise for the development of safe and low-cost protocols for quantification of MD suitable for clinical applications.

Defining the E-cadherin repressor interactome in epithelial-mesenchymal transition: the PMC42 model as a case study.

Epithelial-mesenchymal transition (EMT) is a feature of migratory cellular processes in all stages of life, including embryonic development and wound healing. Importantly, EMT features cluster with disease states such as chronic fibrosis and cancer. The dissolution of the E-cadherin-mediated adherens junction (AJ) is a key preliminary step in EMT and may occur early or late in the growing epithelial tumour. This is a first step for tumour cells towards stromal invasion, intravasation, extravasation and distant metastasis. The AJ may be inactivated in EMT by directed E-cadherin cleavage; however, it is increasingly evident that the majority of AJ changes are transcriptional and mediated by an expanding group of transcription factors acting directly or indirectly to repress E-cadherin expression. A review of the current literature has revealed that these factors may regulate each other in a hierarchical pattern where Snail1 (formerly Snail) and Snail2 (formerly Slug) are initially induced, leading to the activation of Zeb family members, TCF3, TCF4, Twist, Goosecoid and FOXC2. Within this general pathway, many inter-regulatory relationships have been defined which may be important in maintaining the EMT phenotype. This may be important given the short half-life of Snail1 protein. We have investigated these inter-regulatory relationships in the mesenchymal breast carcinoma cell line PMC42 (also known as PMC42ET) and its epithelial derivative, PMC42LA. This review also discusses several newly described regulators of E-cadherin repressors including oestrogen receptor-α and new discoveries in hypoxia- and growth factor-induced EMT. Finally, we evaluated how these findings may influence approaches to current cancer treatment.

Epithelial mesenchymal transition traits in human breast cancer cell lines parallel the CD44(hi/)CD24 (lo/-) stem cell phenotype in human breast cancer.

We review here the recently emerging relationship between epithelial-mesenchymal transition (EMT) and breast cancer stem cells (BCSC), and provide analyses of published data on human breast cancer cell lines, supporting their utility as a model for the EMT/BCSC state. Genome-wide transcriptional profiling of these cell lines has confirmed the existence of a subgroup with mesenchymal tendencies and enhanced invasive properties ('Basal B'/Mesenchymal), distinct from subgroups with either predominantly luminal ('Luminal') or mixed basal/luminal ('Basal A') features (Neve et al. Cancer Cell, 2006). A literature-derived EMT gene signature has shown specific enrichment within the Basal B subgroup of cell lines, consistent with their over-expression of various EMT transcriptional drivers. Basal B cell lines are found to resemble BCSC, being CD44(high)CD24(low). Moreover, gene products that distinguish Basal B from Basal A and Luminal cell lines (Basal B Discriminators) showed close concordance with those that define BCSC isolated from clinical material, as reported by Shipitsin et al. (Cancer Cell, 2007). CD24 mRNA levels varied across Basal B cell lines, correlating with other Basal B Discriminators. Many gene products correlating with CD24 status in Basal B cell lines were also differentially expressed in isolated BCSC. These findings confirm and extend the importance of the cellular product of the EMT with Basal B cell lines, and illustrate the value of analysing these cell lines for new leads that may improve breast cancer outcomes. Gene products specific to Basal B cell lines may serve as tools for the detection, quantification, and analysis of BCSC/EMT attributes.

Diversity of Epithelial-Mesenchymal Phenotypes in Circulating Tumour Cells from Prostate Cancer Patient-Derived Xenograft Models.

Metastasis is the leading cause of cancer-related deaths worldwide. The epithelial-mesenchymal plasticity (EMP) status of primary tumours has relevance to metastatic potential and therapy resistance. Circulating tumour cells (CTCs) provide a window into the metastatic process, and molecular characterisation of CTCs in comparison to their primary tumours could lead to a better understanding of the mechanisms involved in the metastatic cascade. In this study, paired blood and tumour samples were collected from four prostate cancer patient-derived xenograft (PDX) models (BM18, LuCaP70, LuCaP96, LuCaP105) and assessed using an EMP-focused, 42 gene human-specific, nested quantitative RT-PCR assay. CTC burden varied amongst the various xenograft models with LuCaP96 having the highest number of CTCs per mouse (mean: 704; median: 31) followed by BM18 (mean: 101; median: 21), LuCaP70 (mean: 73; median: 16) and LuCaP105 (mean: 57; median: 6). A significant relationship was observed between tumour size and CTC number (p = 0.0058). Decreased levels of kallikrein-related peptidase 3 (KLK3) mRNA (which encodes prostate-specific antigen; PSA) were observed in CTC samples from all four models compared to their primary tumours. Both epithelial- and mesenchymal-associated genes were commonly expressed at higher levels in CTCs compared to the bulk primary tumour, although some common EMT-associated genes (CDH1, VIM, EGFR, EPCAM) remained unchanged. Immunofluorescence co-staining for pan-cytokeratin (KRT) and vimentin (VIM) indicated variable proportions of CTCs across the full EMP axis, even in the same model. EMP hybrids predominated in the BM18 and LuCaP96 models, but were not detected in the LuCaP105 model, and variable numbers of KRT+ and human VIM+ cells were observed in each model. SERPINE1, which encodes plasminogen activator inhibitor-1 (PAI-1), was enriched at the RNA level in CTCs compared to primary tumours and was the most commonly expressed mesenchymal gene in the CTCs. Co-staining for SERPINE1 and KRT revealed SERPINE1+ cells in 7/11 samples, six of which had SERPINE+KRT+ CTCs. Cell size variation was observed in CTCs. The majority of samples (8/11) contained larger CTCs ranging from 15.3 to 37.8 µm, whilst smaller cells (10.7 ± 4.1 µm, similar in size to peripheral blood mononuclear cells (PBMCs)) were identified in 6 of 11 samples. CTC clusters were also identified in 9/11 samples, containing 2-100 CTCs per cluster. Where CTC heterogeneity was observed in the clusters, epithelial-like cells (KRT+VIM-) were located on the periphery of the cluster, forming a layer around hybrid (KRT+VIM+) or mesenchymal-like (KRT-VIM+) cells. The CTC heterogeneity observed in these models emphasises the complexity in CTC isolation and classification and supports the increasingly recognised importance of the epithelial-mesenchymal hybrid state in cancer progression and metastasis.

RASSF1A Suppression as a Potential Regulator of Mechano-Pathobiology Associated with Mammographic Density in BRCA Mutation Carriers.

High mammographic density (MD) increases breast cancer (BC) risk and creates a stiff tissue environment. BC risk is also increased in BRCA1/2 gene mutation carriers, which may be in part due to genetic disruption of the tumour suppressor gene Ras association domain family member 1 (RASSF1A), a gene that is also directly regulated by tissue stiffness. High MD combined with BRCA1/2 mutations further increase breast cancer risk, yet BRCA1/2 mutations alone or in combination do not increase MD. The molecular basis for this additive effect therefore remains unclear. We studied the interplay between MD, stiffness, and BRCA1/2 mutation status in human mammary tissue obtained after prophylactic mastectomy from women at risk of developing BC. Our results demonstrate that RASSF1A expression increased in MCF10DCIS.com cell cultures with matrix stiffness up until ranges corresponding with BiRADs 4 stiffnesses (~16 kPa), but decreased in higher stiffnesses approaching malignancy levels (>50 kPa). Similarly, higher RASSF1A protein was seen in these cells when co-cultivated with high MD tissue in murine biochambers. Conversely, local stiffness, as measured by collagen I versus III abundance, repressed RASSF1A protein expression in BRCA1, but not BRCA2 gene mutated tissues; regional density as measured radiographically repressed RASSF1A in both BRCA1/2 mutated tissues. The combinatory effect of high MD and BRCA mutations on breast cancer risk may be due to RASSF1A gene repression in regions of increased tissue stiffness.

Corrigendum: Heparanase Promotes Syndecan-1 Expression to Mediate Fibrillar Collagen and Mammographic Density in Human Breast Tissue Cultured ex vivo.

[This corrects the article DOI: 10.3389/fcell.2020.00599.].

Applications of RNA from circulating tumor cells.

Circulating tumour cells (CTCs) are shed into the bloodstream from both primary and secondary tumours and provide a non-invasive means to study tumor progression and response to treatment. Assessment of ribonucleic acid (RNA) and monitoring dynamic changes in gene expression profiles of CTCs extends their clinical and prognostic power and establish their role in guiding treatment. Among these methods, droplet digital (RT-ddPCR) technique provides a high sensitivity and detectibility of CTCs. RNA-sequencing (RNAseq) is the most comprehensive method, that would allow the simultaneous measurement of a large number of genes and theoretically the whole transcriptome. Since CTCs are heterogeneous in nature, single cell RNAseq methods are very valuable in assessing population dynamics and functional states of CTCs. While RNA in situ hybridization (RNA-ISH) is used relatively less frequently, it also allows for the assessment of expression of multiple genes within individual CTCs. Epithelial to Mesenchymal Transition (EMT) is a major contributor to metastasis, providing a mechanism to allow cells to become migratory and invasive, and to survive in the bloodstream. Monitoring CTCs undergoing EMT may lead to improvement in their prognostic and predictive power. Here, we review various RNA analysis of CTCs and those that undergo EMT and their application in diagnosis, prognosis and treatment of cancers.

Identifying Therapies to Combat Epithelial Mesenchymal Plasticity-Associated Chemoresistance to Conventional Breast Cancer Therapies Using An shRNA Library Screen.

BACKGROUND: Breast cancer (BC) is a heterogeneous disease for which the commonly used chemotherapeutic agents primarily include the anthracyclines (doxorubicin, epirubicin), microtubule inhibitors (paclitaxel, docetaxel, eribulin), and alkylating agents (cyclophosphamide). While these drugs can be highly effective, metastatic tumours are frequently refractory to treatment or become resistant upon tumour relapse. METHODS: We undertook a cell polarity/epithelial mesenchymal plasticity (EMP)-enriched short hairpin RNA (shRNA) screen in MDA-MB-468 breast cancer cells to identify factors underpinning heterogeneous responses to three chemotherapeutic agents used clinically in breast cancer: Doxorubicin, docetaxel, and eribulin. shRNA-transduced cells were treated for 6 weeks with the EC10 of each drug, and shRNA representation assessed by deep sequencing. We first identified candidate genes with depleted shRNA, implying that their silencing could promote a response. Using the Broad Institute's Connectivity Map (CMap), we identified partner inhibitors targeting the identified gene families that may induce cell death in combination with doxorubicin, and tested them with all three drug treatments. RESULTS: In total, 259 shRNAs were depleted with doxorubicin treatment (at p < 0.01), 66 with docetaxel, and 25 with eribulin. Twenty-four depleted hairpins overlapped between doxorubicin and docetaxel, and shRNAs for TGFB2, RUNX1, CCDC80, and HYOU1 were depleted across all the three drug treatments. Inhibitors of MDM/TP53, TGFBR, and FGFR were identified by CMap as the top pharmaceutical perturbagens and we validated the combinatorial benefits of the TGFBR inhibitor (SB525334) and MDM inhibitor (RITA) with doxorubicin treatment, and also observed synergy between the inhibitor SB525334 and eribulin in MDA-MB-468 cells. CONCLUSIONS: Taken together, a cell polarity/EMP-enriched shRNA library screen identified relevant gene products that could be targeted alongside current chemotherapeutic agents for the treatment of invasive BC.

Heparanase Promotes Syndecan-1 Expression to Mediate Fibrillar Collagen and Mammographic Density in Human Breast Tissue Cultured ex vivo.

Mammographic density (MD) is a strong and independent factor for breast cancer (BC) risk and is increasingly associated with BC progression. We have previously shown in mice that high MD, which is characterized by the preponderance of a fibrous stroma, facilitates BC xenograft growth and metastasis. This stroma is rich in extracellular matrix (ECM) factors, including heparan sulfate proteoglycans (HSPGs), such as the BC-associated syndecan-1 (SDC1). These proteoglycans tether growth factors, which are released by heparanase (HPSE). MD is positively associated with estrogen exposure and, in cell models, estrogen has been implicated in the upregulation of HPSE, the activity of which promotes SDC expression. Herein we describe a novel measurement approach (single-sided NMR) using a patient-derived explant (PDE) model of normal human (female) mammary tissue cultured ex vivo to investigate the role(s) of HPSE and SDC1 on MD. Relative HSPG gene and protein analyses determined in patient-paired high vs. low MD tissues identified SDC1 and SDC4 as potential mediators of MD. Using the PDE model we demonstrate that HPSE promotes SDC1 rather than SDC4 expression and cleavage, leading to increased MD. In this model system, synstatin (SSTN), an SDC1 inhibitory peptide designed to decouple SDC1-ITGαvß3 parallel collagen alignment, reduced the abundance of fibrillar collagen as assessed by picrosirius red viewed under polarized light, and reduced MD. Our results reveal a potential role for HPSE in maintaining MD via its direct regulation of SDC1, which in turn physically tethers collagen into aligned fibers characteristic of MD. We propose that inhibitors of HPSE and/or SDC1 may afford an opportunity to reduce MD in high BC risk individuals and reduce MD-associated BC progression in conjunction with established BC therapies.

Staurosporine augments EGF-mediated EMT in PMC42-LA cells through actin depolymerisation, focal contact size reduction and Snail1 induction - a model for cross-modulation.

BACKGROUND: A feature of epithelial to mesenchymal transition (EMT) relevant to tumour dissemination is the reorganization of actin cytoskeleton/focal contacts, influencing cellular ECM adherence and motility. This is coupled with the transcriptional repression of E-cadherin, often mediated by Snail1, Snail2 and Zeb1/deltaEF1. These genes, overexpressed in breast carcinomas, are known targets of growth factor-initiated pathways, however it is less clear how alterations in ECM attachment cross-modulate to regulate these pathways. EGF induces EMT in the breast cancer cell line PMC42-LA and the kinase inhibitor staurosporine (ST) induces EMT in embryonic neural epithelial cells, with F-actin de-bundling and disruption of cell-cell adhesion, via inhibition of aPKC. METHODS: PMC42-LA cells were treated for 72 h with 10 ng/ml EGF, 40 nM ST, or both, and assessed for expression of E-cadherin repressor genes (Snail1, Snail2, Zeb1/deltaEF1) and EMT-related genes by QRT-PCR, multiplex tandem PCR (MT-PCR) and immunofluorescence +/- cycloheximide. Actin and focal contacts (paxillin) were visualized by confocal microscopy. A public database of human breast cancers was assessed for expression of Snail1 and Snail2 in relation to outcome. RESULTS: When PMC42-LA were treated with EGF, Snail2 was the principal E-cadherin repressor induced. With ST or ST+EGF this shifted to Snail1, with more extreme EMT and Zeb1/deltaEF1 induction seen with ST+EGF. ST reduced stress fibres and focal contact size rapidly and independently of gene transcription. Gene expression analysis by MT-PCR indicated that ST repressed many genes which were induced by EGF (EGFR, CAV1, CTGF, CYR61, CD44, S100A4) and induced genes which alter the actin cytoskeleton (NLF1, NLF2, EPHB4). Examination of the public database of breast cancers revealed tumours exhibiting higher Snail1 expression have an increased risk of disease-recurrence. This was not seen for Snail2, and Zeb1/deltaEF1 showed a reverse correlation with lower expression values being predictive of increased risk. CONCLUSION: ST in combination with EGF directed a greater EMT via actin depolymerisation and focal contact size reduction, resulting in a loosening of cell-ECM attachment along with Snail1-Zeb1/deltaEF1 induction. This appeared fundamentally different to the EGF-induced EMT, highlighting the multiple pathways which can regulate EMT. Our findings add support for a functional role for Snail1 in invasive breast cancer.

Multi-Omics Characterization of the Spontaneous Mesenchymal-Epithelial Transition in the PMC42 Breast Cancer Cell Lines.

Epithelial-mesenchymal plasticity (EMP), encompassing epithelial-mesenchymal transition (EMT) and mesenchymal-epithelial transition (MET), are considered critical events for cancer metastasis. We investigated chromosomal heterogeneity and chromosomal instability (CIN) profiles of two sister PMC42 breast cancer (BC) cell lines to assess the relationship between their karyotypes and EMP phenotypic plasticity. Karyotyping by GTG banding and exome sequencing were aligned with SWATH quantitative proteomics and existing RNA-sequencing data from the two PMC42 cell lines; the mesenchymal, parental PMC42-ET cell line and the spontaneously epithelially shifted PMC42-LA daughter cell line. These morphologically distinct PMC42 cell lines were also compared with five other BC cell lines (MDA-MB-231, SUM-159, T47D, MCF-7 and MDA-MB-468) for their expression of EMP and cell surface markers, and stemness and metabolic profiles. The findings suggest that the epithelially shifted cell line has a significantly altered ploidy of chromosomes 3 and 13, which is reflected in their transcriptomic and proteomic expression profiles. Loss of the TGFßR2 gene from chromosome 3 in the epithelial daughter cell line inhibits its EMT induction by TGF-ß stimulus. Thus, integrative 'omics' characterization established that the PMC42 system is a relevant MET model and provides insights into the regulation of phenotypic plasticity in breast cancer.

Interrogation of Phenotypic Plasticity between Epithelial and Mesenchymal States in Breast Cancer.

Dynamic interconversions between transitional epithelial and mesenchymal states underpin the epithelial mesenchymal plasticity (EMP) seen in some carcinoma cell systems. We have delineated epithelial and mesenchymal subpopulations existing within the PMC42-LA breast cancer cell line by their EpCAM expression. These purified but phenotypically plastic states, EpCAMHigh (epithelial) and EpCAMLow (mesenchymal), have the ability to regain the phenotypic equilibrium of the parental population (i.e., 80% epithelial and 20% mesenchymal) over time, although the rate of reversion in the mesenchymal direction (epithelial-mesenchymal transition; EMT) is higher than that in the epithelial direction (mesenchymal-epithelial transition; MET). Single-cell clonal propagation was implemented to delineate the molecular and cellular features of this intrinsic heterogeneity with respect to EMP flux. The dynamics of the phenotypic proportions of epithelial and mesenchymal states in single-cell generated clones revealed clonal diversity and intrinsic plasticity. Single cell-derived clonal progenies displayed differences in their functional attributes of proliferation, stemness marker (CD44/CD24), migration, invasion and chemo-sensitivity. Interrogation of genomic copy number variations (CNV) with whole exome sequencing (WES) in the context of chromosome count from metaphase spread indicated that chromosomal instability was not influential in driving intrinsic phenotypic plasticity. Overall, these findings reveal the stochastic nature of both the epithelial and mesenchymal subpopulations, and the single cell-derived clones for differential functional attributes.