Identification of a core EMT signature that separates basal-like breast cancers into partial- and post-EMT subtypes.

Epithelial-mesenchymal transition (EMT) is a cellular plasticity program critical for embryonic development and tissue regeneration, and aberrant EMT is associated with disease including cancer. The high degree of plasticity in the mammary epithelium is reflected in extensive heterogeneity among breast cancers. Here, we have analyzed RNA-sequencing data from three different mammary epithelial cell line-derived EMT models and identified a robust mammary EMT gene expression signature that separates breast cancers into distinct subgroups. Most strikingly, the basal-like breast cancers form two subgroups displaying partial-EMT and post-EMT gene expression patterns. We present evidence that key EMT-associated transcription factors play distinct roles at different stages of EMT in mammary epithelial cells.

ALYREF, a novel factor involved in breast carcinogenesis, acts through transcriptional and post-transcriptional mechanisms selectively regulating the short NEAT1 isoform.

The RNA-binding protein ALYREF (THOC4) is involved in transcriptional regulation and nuclear mRNA export, though its role and molecular mode of action in breast carcinogenesis are completely unknown. Here, we identified high ALYREF expression as a factor for poor survival in breast cancer patients. ALYREF significantly influenced cellular growth, apoptosis and mitochondrial energy metabolism in breast cancer cells as well as breast tumorigenesis in orthotopic mouse models. Transcriptional profiling, phenocopy and rescue experiments identified the short isoform of the lncRNA NEAT1 as a molecular trigger for ALYREF effects in breast cancer. Mechanistically, we found that ALYREF binds to the NEAT1 promoter region to enhance the global NEAT1 transcriptional activity. Importantly, by stabilizing CPSF6, a protein that selectively activates the post-transcriptional generation of the short isoform of NEAT1, as well as by direct binding and stabilization of the short isoform of NEAT1, ALYREF selectively fine-tunes the expression of the short NEAT1 isoform. Overall, our study describes ALYREF as a novel factor contributing to breast carcinogenesis and identifies novel molecular mechanisms of regulation the two isoforms of NEAT1.

Serglycin Is Involved in TGF-ß Induced Epithelial-Mesenchymal Transition and Is Highly Expressed by Immune Cells in Breast Cancer Tissue.

Serglycin is a proteoglycan highly expressed by immune cells, in which its functions are linked to storage, secretion, transport, and protection of chemokines, proteases, histamine, growth factors, and other bioactive molecules. In recent years, it has been demonstrated that serglycin is also expressed by several other cell types, such as endothelial cells, muscle cells, and multiple types of cancer cells. Here, we show that serglycin expression is upregulated in transforming growth factor beta (TGF-ß) induced epithelial-mesenchymal transition (EMT). Functional studies provide evidence that serglycin plays an important role in the regulation of the transition between the epithelial and mesenchymal phenotypes, and it is a significant EMT marker gene. We further find that serglycin is more expressed by breast cancer cell lines with a mesenchymal phenotype as well as the basal-like subtype of breast cancers. By examining immune staining and single cell sequencing data of breast cancer tissue, we show that serglycin is highly expressed by infiltrating immune cells in breast tumor tissue.

Expression and functions of long non-coding RNA NEAT1 and isoforms in breast cancer.

NEAT1 is a highly abundant nuclear architectural long non-coding RNA. There are two overlapping NEAT1 isoforms, NEAT1_1 and NEAT1_2, of which the latter is an essential scaffold for the assembly of a class of nuclear ribonucleoprotein bodies called paraspeckles. Paraspeckle formation is elevated by a wide variety of cellular stressors and in certain developmental processes, either through transcriptional upregulation of the NEAT1 gene or through a switch from NEAT1_1 to NEAT1_2 isoform production. In such conditions, paraspeckles modulate cellular processes by sequestering proteins or RNA molecules. NEAT1 is abnormally expressed in many cancers and a growing body of evidence suggests that, in many cases, high NEAT1 levels are associated with therapy resistance and poor clinical outcome. Here we review the current knowledge of NEAT1 expression and functions in breast cancer, highlighting its established role in postnatal mammary gland development. We will discuss possible isoform-specific roles of NEAT1_1 and NEAT1_2 in different breast cancer subtypes, which critically needs to be considered when studying NEAT1 and breast cancer.

Current Status of Circulating Tumor Cells, Circulating Tumor DNA, and Exosomes in Breast Cancer Liquid Biopsies.

Breast cancer is the most common cancer among women worldwide. Although the five-, ten- and fifteen-year survival rates are good for breast cancer patients diagnosed with early-stage disease, some cancers recur many years after completion of primary therapy. Tumor heterogeneity and clonal evolution may lead to distant metastasis and therapy resistance, which are the main causes of breast cancer-associated deaths. In the clinic today, imaging techniques like mammography and tissue biopsies are used to diagnose breast cancer. Even though these methods are important in primary diagnosis, they have limitations when it comes to longitudinal monitoring of residual disease after treatment, disease progression, therapy responses, and disease recurrence. Over the last few years, there has been an increasing interest in the diagnostic, prognostic, and predictive potential of circulating cancer-derived material acquired through liquid biopsies in breast cancer. Thanks to the development of sensitive devices and platforms, a variety of tumor-derived material, including circulating cancer cells (CTCs), circulating DNA (ctDNA), and biomolecules encapsulated in extracellular vesicles, can now be extracted and analyzed from body fluids. Here we will review the most recent studies on breast cancer, demonstrating the clinical potential and utility of CTCs and ctDNA. We will also review literature illustrating the potential of circulating exosomal RNA and proteins as future biomarkers in breast cancer. Finally, we will discuss some of the advantages and limitations of liquid biopsies and the future perspectives of this field in breast cancer management.

The expression of the long NEAT1_2 isoform is associated with human epidermal growth factor receptor 2-positive breast cancers.

The long non-coding RNA NEAT1 locus is transcribed into two overlapping isoforms, NEAT1_1 and NEAT1_2, of which the latter is essential for the assembly of nuclear paraspeckles. NEAT1 is abnormally expressed in a wide variety of human cancers. Emerging evidence suggests that the two isoforms have distinct functions in gene expression regulation, and recently it was shown that NEAT1_2, but not NEAT1_1, expression predicts poor clinical outcome in cancer. Here, we report that NEAT1_2 expression correlates with HER2-positive breast cancers and high-grade disease. We provide evidence that NEAT1_1 and NEAT1_2 have distinct expression pattern among different intrinsic breast cancer subtypes. Finally, we show that NEAT1_2 expression and paraspeckle formation increase upon lactation in humans, confirming what has previously been demonstrated in mice.

Cancer-specific SNPs originate from low-level heteroplasmic variants in human mitochondrial genomes of a matched cell line pair.

Low-level mitochondrial heteroplasmy is a common phenomenon in both normal and cancer cells. Here, we investigate the link between low-level heteroplasmy and mitogenome mutations in a human breast cancer matched cell line by high-throughput sequencing. We identified 23 heteroplasmic sites, of which 15 were common between normal cells (Hs578Bst) and cancer cells (Hs578T). Most sites were clustered within the highly conserved Complex IV and ribosomal RNA genes. Two heteroplasmic variants in normal cells were found as fixed mutations in cancer cells. This indicates a positive selection of these variants in cancer cells. RNA-Seq analysis identified upregulated L-strand specific transcripts in cancer cells, which include three mitochondrial long non-coding RNA molecules. We hypothesize that this is due to two cancer cell-specific mutations in the control region.

The long noncoding RNA NEAT1 and nuclear paraspeckles are up-regulated by the transcription factor HSF1 in the heat shock response.

The long noncoding RNA (lncRNA) NEAT1 (nuclear enriched abundant transcript 1) is the architectural component of nuclear paraspeckles, and it has recently gained considerable attention as it is abnormally expressed in pathological conditions such as cancer and neurodegenerative diseases. NEAT1 and paraspeckle formation are increased in cells upon exposure to a variety of environmental stressors and believed to play an important role in cell survival. The present study was undertaken to further investigate the role of NEAT1 in cellular stress response pathways. We show that NEAT1 is a novel target gene of heat shock transcription factor 1 (HSF1) and is up-regulated when the heat shock response pathway is activated by sulforaphane (SFN) or elevated temperature. HSF1 binds specifically to a newly identified conserved heat shock element in the NEAT1 promoter. In line with this, SFN induced the formation of NEAT1-containing paraspeckles via an HSF1-dependent mechanism. HSF1 plays a key role in the cellular response to proteotoxic stress by promoting the expression of a series of genes, including those encoding molecular chaperones. We have found that the expression of HSP70, HSP90, and HSP27 is amplified and sustained during heat shock in NEAT1-depleted cells compared with control cells, indicating that NEAT1 feeds back via an unknown mechanism to regulate HSF1 activity. This interrelationship is potentially significant in human diseases such as cancer and neurodegenerative disorders.

ICAM1 expression is induced by proinflammatory cytokines and associated with TLS formation in aggressive breast cancer subtypes.

Intratumoral formation of tertiary lymphoid structures (TLS) within the tumor microenvironment is considered to be a consequence of antigen challenge during anti-tumor responses. Intracellular adhesion molecule 1 (ICAM1) has been implicated in a variety of immune and inflammatory responses, in addition to associate with triple negative breast cancer (TNBC). In this study, we detected TLS in the aggressive tumor phenotypes TNBC, HER2+ and luminal B, whereas the TLS negative group contained solely tumors of the luminal A subtype. We show that ICAM1 is exclusively expressed in TNBC and HER2 enriched subtypes known to be associated with inflammation and the formation of TLS. Furthermore, cell from normal mammary epithelium and breast cancer cell lines expressed ICAM1 upon stimulation with the proinflammatory cytokines TNFα, IL1ß and IFNγ. ICAM1 overexpression was induced in MCF7, MDA-MB-468 and SK-BR-3 cells regardless of hormone receptor status. Taken together, our findings show that ICAM1 is expressed in aggressive subtypes of breast cancer and its expression is inducible by well-known proinflammatory cytokines. ICAM1 may be an attractive molecular target for TNBC, but further investigations elucidating the role of ICAM1 in targeted therapies have to take into consideration selective subtypes of breast cancer.

Regulation of atypical MAP kinases ERK3 and ERK4 by the phosphatase DUSP2.

The atypical MAP kinases ERK3 and ERK4 are activated by phosphorylation of a serine residue lying within the activation loop signature sequence S-E-G. However, the regulation of ERK3 and ERK4 phosphorylation and activity is poorly understood. Here we report that the inducible nuclear dual-specificity MAP kinase phosphatase (MKP) DUSP2, a known regulator of the ERK and p38 MAPKs, is unique amongst the MKP family in being able to bind to both ERK3 and ERK4. This interaction is mediated by a conserved common docking (CD) domain within the carboxyl-terminal domains of ERK3 and ERK4 and the conserved kinase interaction motif (KIM) located within the non-catalytic amino terminus of DUSP2. This interaction is direct and results in the dephosphorylation of ERK3 and ERK4 and the stabilization of DUSP2. In the case of ERK4 its ability to stabilize DUSP2 requires its kinase activity. Finally, we demonstrate that expression of DUSP2 inhibits ERK3 and ERK4-mediated activation of its downstream substrate MK5. We conclude that the activity of DUSP2 is not restricted to the classical MAPK pathways and that DUSP2 can also regulate the atypical ERK3/4-MK5 signalling pathway in mammalian cells.

Metabolic re-wiring of isogenic breast epithelial cell lines following epithelial to mesenchymal transition.

Epithelial to mesenchymal transition (EMT) has implications in tumor progression and metastasis. Metabolic alterations have been described in cancer development but studies focused on the metabolic re-wiring that takes place during EMT are still limited. We performed metabolomics profiling of a breast epithelial cell line and its EMT derived mesenchymal phenotype to create genome-scale metabolic models descriptive of both cell lines. Glycolysis and OXPHOS were higher in the epithelial phenotype while amino acid anaplerosis and fatty acid oxidation fueled the mesenchymal phenotype. Through comparative bioinformatics analysis, PPAR-γ1, PPAR- γ2 and AP-1 were found to be the most influential transcription factors associated with metabolic re-wiring. In silico gene essentiality analysis predicts that the LAT1 neutral amino acid transporter is essential for mesenchymal cell survival. Our results define metabolic traits that distinguish an EMT derived mesenchymal cell line from its epithelial progenitor and may have implications in cancer progression and metastasis. Furthermore, the tools presented here can aid in identifying critical metabolic nodes that may serve as therapeutic targets aiming to prevent EMT and inhibit metastatic dissemination.

Distinct Small RNA Signatures in Extracellular Vesicles Derived from Breast Cancer Cell Lines.

Breast cancer is a heterogeneous disease, and different subtypes of breast cancer show distinct cellular morphology, gene expression, metabolism, motility, proliferation, and metastatic potential. Understanding the molecular features responsible for this heterogeneity is important for correct diagnosis and better treatment strategies. Extracellular vesicles (EVs) and their associated molecules have gained much attention as players in intercellular communication, ability to precondition specific organs for metastatic invasion, and for their potential role as circulating cancer biomarkers. EVs are released from the cells and contain proteins, DNA, and long and small RNA species. Here we show by high-throughput small RNA-sequencing that EVs from nine different breast cancer cell lines share common characteristics in terms of small RNA content that are distinct from their originating cells. Most strikingly, a highly abundant small RNA molecule derived from the nuclear 28S rRNA is vastly enriched in EVs. The miRNA profiles in EVs correlate with the cellular miRNA expression pattern, but with a few exceptions that includes miR-21. This cancer-associated miRNA is retained in breast cancer cell lines. Finally, we report that EVs from breast cancer cell lines cluster together based on their small RNA signature when compared to EVs derived from other cancer cell lines. Altogether, our data demonstrate that breast cancer cell lines manifest a specific small RNA signature in their released EVs. This opens up for further evaluation of EVs as breast cancer biomarkers.

New insights into the activation, interaction partners and possible functions of MK5/PRAK.

MAP kinase-activated protein kinase 5 (MK5) was first described as a downstream target of the p38 MAP kinase pathway leading to its alternative acronym of p38-regulated/activated protein kinase (PRAK). However, since the discovery that MK5 is a bona fide interaction partner of the atypical MAP kinases ERK3 and ERK4 and that this interaction leads to both the activation and subcellular relocalisation of MK5, there has been considerable debate as to the relative roles of these MAPK pathways in mediating the activation and biological functions of MK5. Here we discuss recent progress in defining novel upstream components of the ERK3/ERK4 signalling pathway, our increased understanding of the mechanism by which MK5 interacts with and is activated by ERK3 and ERK4, and the discovery of novel interaction partners for MK5. Finally, we review recent literature that suggests novel biological functions for MK5 in a range of physiological and pathophysiological conditions including neuronal function and cancer.

Performance comparison of digital microRNA profiling technologies applied on human breast cancer cell lines.

MicroRNA profiling represents an important first-step in deducting individual RNA-based regulatory function in a cell, tissue, or at a specific developmental stage. Currently there are several different platforms to choose from in order to make the initial miRNA profiles. In this study we investigate recently developed digital microRNA high-throughput technologies. Four different platforms were compared including next generation SOLiD ligation sequencing and Illumina HiSeq sequencing, hybridization-based NanoString nCounter, and miRCURY locked nucleic acid RT-qPCR. For all four technologies, full microRNA profiles were generated from human cell lines that represent noninvasive and invasive tumorigenic breast cancer. This study reports the correlation between platforms, as well as a more extensive analysis of the accuracy and sensitivity of data generated when using different platforms and important consideration when verifying results by the use of additional technologies. We found all the platforms to be highly capable for microRNA analysis. Furthermore, the two NGS platforms and RT-qPCR all have equally high sensitivity, and the fold change accuracy is independent of individual miRNA concentration for NGS and RT-qPCR. Based on these findings we propose new guidelines and considerations when performing microRNA profiling.

Docking of PRAK/MK5 to the atypical MAPKs ERK3 and ERK4 defines a novel MAPK interaction motif.

ERK3 and ERK4 are atypical MAPKs in which the canonical TXY motif within the activation loop of the classical MAPKs is replaced by SEG. Both ERK3 and ERK4 bind, translocate, and activate the MAPK-activated protein kinase (MK) 5. The classical MAPKs ERK1/2 and p38 interact with downstream MKs (RSK1-3 and MK2-3, respectively) through conserved clusters of acidic amino acids, which constitute the common docking (CD) domain. In contrast to the classical MAPKs, the interaction between ERK3/4 and MK5 is strictly dependent on phosphorylation of the SEG motif of these kinases. Here we report that the conserved CD domain is dispensable for the interaction of ERK3 and ERK4 with MK5. Using peptide overlay assays, we have defined a novel MK5 interaction motif (FRIEDE) within both ERK4 and ERK3 that is essential for binding to the C-terminal region of MK5. This motif is located within the L16 extension lying C-terminal to the CD domain in ERK3 and ERK4 and a single isoleucine to lysine substitution in FRIEDE totally abrogates binding, activation, and translocation of MK5 by both ERK3 and ERK4. These findings are the first to demonstrate binding of a physiological substrate via this region of the L16 loop in a MAPK. Furthermore, the link between activation loop phosphorylation and accessibility of the FRIEDE interaction motif suggests a switch mechanism for these atypical MAPKs in which the phosphorylation status of the activation loop regulates the ability of both ERK3 and ERK4 to bind to a downstream effector.

Does MK5 reconcile classical and atypical MAP kinases?

MAP kinase-activated protein kinase 5 (MK5) was originally described as a protein kinase activated downstream of the p38 MAP kinase and is also named p38-regulated/activated protein kinase (PRAK). However, while MK5 is most similar in sequence to the two p38 regulated MAPKAP kinases MK2 and MK3, recent data has shown that in contrast to these enzymes MK5 is not activated in response to either cellular stress or pro-inflammatory cytokines. This lack of response to stimuli which cause robust activation of p38 MAP kinase in vivo is supported by data obtained using transgenic mice lacking MK5. Unlike animals lacking MK2 and MK3, MK5 null mice respond normally to endotoxic shock and display an unchanged pattern of cytokine expression in response to LPS. Clues as to the physiological function of MK5 have come from the recent observation that MK5 is uniquely regulated and activated following complex formation with the atypical MAP kinases ERK3 and ERK4. Thus, it is possible that MK5 is unique amongst the MAPKAP kinases in being regulated downstream of signaling pathways other than the classical MAP kinases p38 and ERK1/2.

The Ser(186) phospho-acceptor site within ERK4 is essential for its ability to interact with and activate PRAK/MK5.

ERK (extracellular-signal-regulated kinase) 4 [MAPK (mitogen-activated protein kinase) 4] and ERK3 (MAPK6) are atypical MAPKs. One major difference between these proteins and the classical MAPKs is substitution of the conserved T-X-Y motif within the activation loop by a single phospho-acceptor site within an S-E-G motif. In the present study we report that Ser(186) of the S-E-G motif in ERK4 is phosphorylated in vivo. Kinase-dead ERK4 is also phosphorylated on Ser(186), indicating that an ERK4 kinase, rather than autophosphorylation, is responsible. Co-expression of MK5 [MAPK-activated protein kinase 5; also known as PRAK (p38-regulated/activated kinase)], a physiological target of ERK4, increases phosphorylation of Ser(186). This is not dependent on MK5 activity, but does require interaction between ERK4 and MK5 suggesting that MK5 binding either prevents ERK4 dephosphorylation or facilitates ERK4 kinase activity. ERK4 mutants in which Ser(186) is replaced with either an alanine residue or a phospho-mimetic residue (glutamate) are unable to activate MK5 and Ser(186) is also required for cytoplasmic anchoring of MK5. Both defects seem to reflect an impaired ability of the ERK4 mutants to interact with MK5. We find that there are at least two endogenous pools of wild-type ERK4. One form exhibits reduced mobility when analysed using SDS/PAGE. This is due to MK5-dependent phosphorylation and only this retarded ERK4 species is both phosphorylated on Ser(186) and co-immunoprecipitates with wild-type MK5. We conclude that binding between ERK4 and MK5 facilitates phosphorylation of Ser(186) and stabilization of the ERK4-MK5 complex. This results in phosphorylation and activation of MK5, which in turn phosphorylates ERK4 on sites other than Ser(186) resulting in the observed mobility shift.

Regulation of MAPK-activated protein kinase 5 activity and subcellular localization by the atypical MAPK ERK4/MAPK4.

MAPK-activated protein kinase 5 (MK5) was recently identified as a physiological substrate of the atypical MAPK ERK3. Complex formation between ERK3 and MK5 results in phosphorylation and activation of MK5, concomitant stabilization of ERK3, and the nuclear exclusion of both proteins. However, ablation of ERK3 in HeLa cells using small interfering RNA or in fibroblasts derived from ERK3 null mice reduces the activity of endogenous MK5 by only 50%, suggesting additional mechanisms of MK5 regulation. Here we identify the ERK3-related kinase ERK4 as a bona fide interaction partner of MK5. Binding of ERK4 to MK5 is accompanied by phosphorylation and activation of MK5. Furthermore, complex formation also results in the relocalization of MK5 from nucleus to cytoplasm. However unlike ERK3, ERK4 is a stable protein, and its half-life is not modified by the presence or absence of MK5. Finally, although knock-down of ERK4 protein in HeLa cells reduces endogenous MK5 activity by approximately 50%, a combination of small interfering RNAs targeting both ERK4 and ERK3 causes a further reduction in the MK5 activity by more than 80%. We conclude that MK5 activation is dependent on both ERK3 and ERK4 in these cells and that these atypical MAPKs are both physiological regulators of MK5 activity.

p62/SQSTM1 forms protein aggregates degraded by autophagy and has a protective effect on huntingtin-induced cell death.

Autophagic degradation of ubiquitinated protein aggregates is important for cell survival, but it is not known how the autophagic machinery recognizes such aggregates. In this study, we report that polymerization of the polyubiquitin-binding protein p62/SQSTM1 yields protein bodies that either reside free in the cytosol and nucleus or occur within autophagosomes and lysosomal structures. Inhibition of autophagy led to an increase in the size and number of p62 bodies and p62 protein levels. The autophagic marker light chain 3 (LC3) colocalized with p62 bodies and co-immunoprecipitated with p62, suggesting that these two proteins participate in the same complexes. The depletion of p62 inhibited recruitment of LC3 to autophagosomes under starvation conditions. Strikingly, p62 and LC3 formed a shell surrounding aggregates of mutant huntingtin. Reduction of p62 protein levels or interference with p62 function significantly increased cell death that was induced by the expression of mutant huntingtin. We suggest that p62 may, via LC3, be involved in linking polyubiquitinated protein aggregates to the autophagy machinery.

Interaction codes within the family of mammalian Phox and Bem1p domain-containing proteins.

The Phox and Bem1p (PB1) domain constitutes a recently recognized protein-protein interaction domain found in the atypical protein kinase C (aPKC) isoenzymes, lambda/iota- and zeta PKC; members of mitogen-activated protein kinase (MAPK) modules like MEK5, MEKK2, and MEKK3; and in several scaffold proteins involved in cellular signaling. Among the last group, p62 and Par6 (partitioning-defective 6) are involved in coupling the aPKCs to signaling pathways involved in cell survival, growth control, and cell polarity. By mutation analyses and molecular modeling, we have identified critical residues at the interaction surfaces of the PB1 domains of aPKCs and p62. A basic charge cluster interacts with an acidic loop and helix both in p62 oligomerization and in the aPKC-p62 interaction. Subsequently, we determined the abilities of mammalian PB1 domain proteins to form heteromeric and homomeric complexes mediated by this domain. We report several novel interactions within this family. An interaction between the cell polarity scaffold protein Par6 and MEK5 was found. Furthermore, p62 interacts both with MEK5 and NBR1 in addition to the aPKCs. Evidence for involvement of p62 in MEK5-ERK5 signaling is presented.