Skip to main content
Download PDF
  • Primary research
  • Open access
  • Published:

Re-expression of miR-200s in claudin‐low mammary tumor cells alters cell shape and reduces proliferation and invasion potentially through modulating other miRNAs and SUZ12 regulated genes

Cancer Cell International volume 21, Article number: 89 (2021) Cite this article

Abstract

Background

MicroRNAs are a class of non-coding RNAs that regulate gene expression through binding to mRNAs and preventing their translation. One family of microRNAs known as the miR-200 family is an important regulator of epithelial identity. The miR-200 family consists of five members expressed in two distinct clusters; the miR-200c/141 cluster and the miR-200b/200a/429 cluster. We have found that murine and human mammary tumor cells with claudin-low characteristics are associated with very low levels of all five miR-200s.

Methods

To determine the impact of miR-200s on claudin-low mammary tumor cells, the miR-200c/141 cluster and the miR-200b/200a/429 cluster were stably re-expressed in murine (RJ423) and human (MDA-MB-231) claudin-low mammary tumor cells. Cell proliferation and migration were assessed using BrdU incorporation and transwell migration across Matrigel coated inserts, respectively. miRNA sequencing and RNA sequencing were performed to explore miRNAs and mRNAs regulated by miR-200 re-expression while Enrichr-based pathway analysis was utilized to identify cellular functions modified by miR-200s.

Results

Re-expression of the miR-200s in murine and human claudin-low mammary tumor cells partially restored an epithelial cell morphology and significantly inhibited proliferation and cell invasion in vitro. miRNA sequencing and mRNA sequencing revealed that re-expression of miR-200s altered the expression of other microRNAs and genes regulated by SUZ12 providing insight into the complexity of miR-200 function. SUZ12 is a member of the polycomb repressor complex 2 that suppresses gene expression through methylating histone H3 at lysine 27. Flow cytometry confirmed that re-expression of miR-200s increased histone H3 methylation at lysine 27.

Conclusions

Re-expression of miR-200s in claudin-low mammary tumor cells alters cell morphology and reduces proliferation and invasion, an effect potentially mediated by SUZ12-regulated genes and other microRNAs.

Background

MicroRNAs (miRNAs) are small, non-coding RNA molecules that can be found in both introns and exons [1]. MiRNAs regulate mRNA translation by binding to mRNAs and in turn, target the mRNA for degradation or impair ribosome mediated translation of the mRNA [2, 3]. The interaction between miRNA and mRNA is primarily regulated by nucleotides 1–8 of the miRNA. This region, known as the seed sequence, binds to complementary sequences in the 3’ untranslated regions of mRNAs [4,5,6,7,8,9,10,11,12]. However, more recent studies have shown that non-canonical sites can also influence miRNA-mRNA interaction [13, 14] and miRNAs have been found to directly regulate gene transcription, activating Toll-like receptors and upregulating protein expression [15, 16], further increasing their complexity.

A number of miRNAs have been implicated in regulating breast cancer including the miR-200 family. The miR-200 family consists of 5 members organized into two clusters, the miR-200c/141 cluster and the miR-200b/200a/429 cluster [17,18,19]. miR-200b, miR-200c, and miR-429 share a common seed sequence (AAUACUG) while miR-200a and miR-141 share the same seed sequence (AACACUG) that differs from the seed sequence of miR-200b, miR-200c and miR-429 by one nucleotide [20]. Therefore, both miR-200 clusters express miRNAs with the AAUACUG and AACACUG seed sequences and thus are predicted to share mRNA targets. The most extensively characterized function of the miR-200 family is the regulation of epithelial and mesenchymal phenotypes. miR-200 family members negatively regulate mesenchymal transcription factors such as Zeb1/2, Twist1/2, Snai1/2, [19, 21,22,23] while inducing the expression of epithelial genes [24, 25]. Therefore, loss of miR-200s facilitates epithelial to mesenchymal transition (EMT).

Consistent with their roles in epithelial identity, most studies found that miR-200s are expressed in human luminal breast cancer but lost in triple negative breast cancer (TNBC) [26,27,28,29,30]. TNBCs lack expression of the estrogen receptor (ER), progesterone receptor (PR), and HER2, and have a poorer prognosis than luminal breast cancer [31,32,33,34]. Claudin-low breast cancers share features with TNBC including frequent loss of ER, PR and HER2 expression, high expression of mesenchymal genes and features of an epithelial–mesenchymal transition (EMT) phenotype [35]. In addition, claudin-low tumors express low levels of cell-cell adhesion and tight junction genes including claudins 3, 4, 7.

Our previous work with miR-200s in a murine claudin-low mammary tumor cell line, RJ423, revealed that re-expression of the miR-200b/200a/429 cluster reduced mesenchymal gene expression, restored a more epithelial-like morphology, reduced migration and inhibited tumor growth and metastasis in vivo [36]. In this current study, both the miR-200b/200a/429 and miR-200c/141 clusters were re-expressed in RJ423 cells as well as the human claudin-low breast cancer cell line, MDA-MB-231. Our data showed that when miR-200s are expressed at sufficiently high levels, they can restore a more epithelial like morphology and impair cell migration. miRNA sequencing revealed that re-expression of miR-200s significantly influenced the expression of > 50 miRNAs and RNA sequencing revealed significant alterations in > 1000 mRNAs. One of the molecular features associated with the differentially expressed mRNAs was an alteration in SUZ12 regulated genes. SUZ12 is a core protein of the polycomb repressor complex 2 PRC2 that regulates gene transcription through methylating histone H3 at lysine 27 (H3K27). Therefore, miR-200s appear to influence cell morphology and migration through (i) directly targeting mRNAs, (ii) altering the expression of other miRNAs, and (iii) regulating methylation of histone H3.

Methods

Cell lines and transfections

MCF-7 (cat #HTB-22) and MDA-MB-231 (cat #HTB-26) cells were purchased from ATCC (Manassas, VA). MDA-200c141 and MDA-200ba429 cells were created by transfecting MDA-MB-231 cells with pLenti 4.1 Ex miR200c-141 (cat #35,534, Addgene, Watertown, MA) or pLenti 4.1 Ex miR200b-200a-429 (cat #35,533, Addgene, Watertown, MA), respectively. MDA-231EV cells were created by infecting MDA-MB-231 cells with copGFP control lentiviral particles (cat #sc-108,084, Santa Cruz Biotechnology Inc, Dallas, TX). MCF-7, MDA-MB-231 and variants of MDA-MB-231 cells were maintained in DMEM media supplemented with 10% FBS, 1 mM sodium pyruvate, 4 mM glutamine, 1% antibiotic/antimycotic and 60 µg/ml puromycin (InvivoGen, San Diego, CA). RJ345 and RJ423 cells were created from mammary tumors that developed in distinct MTB-IGFIR transgenic mice and have previously been described [36,37,38]. RJ423-200c/141 cells were created by transfecting RJ423 cells with pLenti 4.1 Ex miR200c-141 (cat #35,534, Addgene, Watertown, MA). RJ423-200c/141 cells were cultured in DMEM media containing 10% FBS, 1 mM sodium pyruvate, 10 mM 4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid (HEPES), 4 mM glutamine, 2 mM hydrocortisone, 5 µg/ml estrogen, 5 µg/ml prolactin, 10 µg/ml EGF, 10 µg/ml insulin, and 1% antibiotic-antimycotic (all supplements purchased from ThermoFisher Scientific, Waltham, MA) as well as 300 µg/ml of puromycin (InvivoGen, San Diego, CA). RJ345, RJ423EV and RJ423-200ba429 cells have been described previously as have the culture conditions for these cells [36].

Crystal violet staining and transwell invasion assays

Crystal violet staining and transwell invasion assays were performed as described in [36] using 25,000 human tumor cells and 30,000 murine tumor cells for the transwell assays.

RNA extraction, taqman qRT-PCR for microRNA expression and qRT-PCR for gene expression

RNA extraction, Taqman qRT-PCR and qRT-PCR for gene expression were performed as described in [37]. All miR-200 Taqman probes were obtained from ThermoFisher Scientific (Waltham, MA). miR-200 expression for the murine mammary tumor cell lines was normalized to sno202 and sno234 while the miR-200 expression of the human breast cancer cell lines was normalized to RNU44 and RNU48. All gene primers were obtained from Bio-Rad Laboratories (Mississauga, ON); Cdh1 (qMmuCID0005843), Hprt (qMmuCED0045738), Snai1 (qMmuCID0024342), Snai2 (qMmuCED0046072), Twist1 (qMmuCED0004065), Twist2 (qMmuCID0009652) Vim (qMmuCID0005527), Zeb1 (qMmuCID0009095), and Zeb2 (qMmuCID0014662). Hprt was used as the housekeeping gene.

miRNA sequencing

miRNA sequencing libraries were generated using NEB Multiplex small RNA library Prep Set for Illumina and sequencing quality was determined using an Agilent 2100 Bioanalyzer. Libraries were sequenced using an Illumina NextSeq 500 instrument. The Q30 scores for all samples were above 93%. Reads were then 3′-adaptor trimmed and filtered ≤ 15 bp reads with cutadapt software (v1.14). Trimmed reads were aligned to the reference genome with bowtie software. miRNA expression levels were calculated using mirdeep2 (v0.0.8) and differentially expressed miRNAs were performed with edgeR (v3.18.1). Library preparation, sequencing and data analysis were performed by Arraystar Inc. (Rockville, MD). Four independent samples were sequenced.

RNA sequencing

RNA sequencing for one set of RJ423EV samples and the RJ423ba429 samples was performed at the Genome Quebec Innovation Centre at McGill University using the Illumina Hiseq 2500 v4 PE125 as previously described [37]. RNA sequencing for a second set of RJ423EV samples as well as RJ423-200c/141, MDAEV, MDA-200c/141 and MDA-200ba429 were performed by Arraystar Inc (Arraystar Inc., Rockville MD). All Fastq files were processed using Genialis software (Genialis Inc, Houston, TX) following the standard RNA-seq pipeline which uses BBDuk to remove adapters and trim reads, STAR to align the reads, and feature counts to generate gene level counts. RNA sequencing has been uploaded to GEO under accession number GSE150107. Note that our original data for RJ423EV and RJ423-200ba429 samples found at GSE113162 [36] were analyzed by the Genome Quebec Innovation Centre at McGill University and thus might differ from the data in this manuscript that was analyzed using Genialis software. Three independent samples were sequenced for the RJ423 variants and four independent samples were sequenced for the MDA-231 variants.

BrdU and H3K27me3 flow cytometry

For the murine cell lines, a FITC BrdU flow kit (BD Biosciences, San Jose, CA, cat #559,619) and for the human cells an APC BrdU flow kit (BD Biosciences, San Jose, CA, cat #552,598) were used following the manufacturer’s protocol. The APC kit was required for the human cell lines as MDA-231EV cell lines express GFP. Briefly, cells were incubated with 1 mM BrdU in fully supplemented media for 45 min. Cells were then fixed, washed and analyzed on an Accuri C6 cytometer (BD Biosciences, San Jose, CA) using a flow rate of 35 μl/min and 20,000 events were collected. For the H3K27me3 analysis, cells were fixed in 4% paraformaldehyde for 15min and then permeabilized in 90% methanol for 10 min. Cells were then incubated with a 1:200 dilution of anti-H3K27me3 (Cell Signaling cat #9733S) or IgG KP isotype control (Cell Signaling cat #3900s) for 1 h followed by a 30 min incubation with the appropriate secondary antibody (Alexa Fluor 488 for the mouse cells and APC for the human cells). Cells were then resuspended in 5% 7-AAD solution for 5 min. Cells were analyzed on an Accuri C6 cytometer (BD Biosciences, San Jose, CA) using a flow rate of 35 μl/min and 20,000 events were collected.

Venn diagrams

Venn diagrams were generated using Venny 2.1.0 [38].

KEGG pathway, gene ontology and ChEA

For RNA sequencing data Enhrichr [39, 40] was used to identify KEGG and Gene Ontology Pathways affected as well as transcription factors implicated in regulating the differentially expressed genes (ENCODE and ChEA Consensus TFs from ChIP-X).

miRNA binding site analysis

Mouse genes with a target score ≥ 50 for either the miR-200c/141 cluster (miR-200c-3p, miR-200c-5p, miR-141-3p, miR-141-5p) or the miR-200b/200a/429 cluster (miR-200a-3p, miR-200a-5p, miR-200b-3p, miR-200b-5p, miR-429-3p, miR-429-5p) were identified in miRDB [41, 42] and compared to the mRNAs differentially expressed in the RNA-seq data of RJ423c141 or RJ423ba429 cells. The same approach was used for the human genes.

Statistics

Statistical significance was determined using an ANOVA followed by a Dunnett’s test using GraphPad Prism 9 (San Diego, CA). For the MCF-7 and MDA-231 variants, all comparisons were made using MDA-231EV cells as the control and for the RJ345 and RJ423 variants, all comparisons were made using RJ423EV cells as the control.

Results

miR-200 expression and the impact on mesenchymal gene expression and cell shape

As shown in Fig. 1, the murine (RJ345, Fig. 1a) and human (MCF-7, Fig. 1c) mammary tumor cell lines with epithelial characteristics expressed high levels of all 5 miR-200 family members compared to the murine (RJ423/RJ423EV) and human (MDA-MB-231/MDA-231EV) claudin-low mammary tumor cell lines, respectively. Transfection of either the miR-200c/141 cluster or miR-200b/200a/429 cluster in RJ423 or MDA-MB-231 cells significantly increased the expression of one or more of the appropriate miR-200s. For example, RJ423ba429 cells expressed significantly higher levels of miR-200b, miR-200a and miR-429 compared to RJ423EV and MDA-231c141 cells expressed significantly higher levels of miR-200c and miR-141 than MDA-231EV cells.

Fig. 1
figure 1

MicroRNA and gene expression. Expression of the miR-200 family (a, c) and epithelial and mesenchymal genes (b, d) in RJ345, RJ423, RJ423EV, RJ423c141 and RJ423ba429 cells (a, b) and MCF-7, MDA-MB-231, MDA-231EV, MDA-231c141 and MDA-231ba429 cells (c, d). miR-200 expression was determined using Taqman qRT-PCR while gene expression was determined using qRT-PCR. The graphs represent the mean ± SEM (n = 3) and the * indicates as significant difference (p < 0.05) from RJ423EV cells (a, b) or MDA-231EV cells (c, d)

Full size image

Since the miR-200 family has been implicated in regulating epithelial to mesenchymal transition (EMT) [19, 21,22,23,24,25], the expression of the epithelial gene Cdh1, and several mesenchymal genes were evaluated. The altered miR-200 expression in RJ423 and MDA-MB-231 cells did not significantly increase the expression of the epithelial gene, Cdh1 and had inconsistent effects on mesenchymal genes (Fig. 1b, d). Only Twist2 and Vim were significantly reduced by miR-200 expression in RJ423 cells (Fig. 1c) and only Snai1 was significantly reduced by miR-200 expression in MDA-MB-231 cells (Fig. 1d) however, the expression of Snai1 was extremely low in all the cell lines. Somewhat surprisingly, re-expression of miR-200c/141 in both RJ423 and MDA-MB-231 cells induced a significant increase in Snai2 expression (Fig. 1b, d). Zeb1, Zeb2 are predicted targets of all five miR-200s while Snai2 is a predicted target of miR-200b, miR-200c and miR-429 according to miRDB (mirdb.org). Zeb1 was reduced approximately twofold in both RJ423ba429 cells compared to RJ423EV cells and MDA-231c141 cells compared to MDA-231EV cells but these differences did not reach statistical significance. The mesenchymal genes evaluated in this study are regulated by other miRNAs in addition to the miR-200 family and several transcription factors and thus it is the balance of transcriptional induction and repression as well as post-transcriptional regulation that would determine gene expression levels.

During routine culture of these cells it was observed that RJ345 (Fig. 2a) and MCF-7 (Fig. 2f) cells displayed a more rounded cell shape typically associated with epithelial cells while RJ423 (Fig. 2b) and MDA-MB-231 (Fig. 2g) cells have a more spindle shaped morphology typically associated with mesenchymal cells. Re-expression of miR-200b/200a/429 in RJ423 cells (Fig. 2e and [36]) and re-expression of miR-200c/141 in MDA-MB-231 cells (Fig. 2i) reverted these cells to a more epithelial morphology. This morphological change was observed in all RJ423ba429 (Fig. 2e) and MDA-231c141 (Fig. 2i) cells and was maintained when the cells were passaged.

Fig. 2
figure 2

Cell morphology. Crystal violet staining of murine (a–e) and human (f–j) mammary tumor cells. Mammary tumor cells with a luminal phenotype (a, f) display a more rounded/cuboidal morphology while mammary tumor cells with a claudin-low phenotype (b, g) display a more spindle-like morphology. Re-expression of the miR-200b/200a/429 cluster in RJ423 cells (e) and the miR-200c/141 cluster in MDA-MB-231 cells (i) reverted these cells to a more luminal morphology. Scale bars are 25 µm

Full size image

To evaluate whether miR-200 re-expression altered cell proliferation, BrdU incorporation was assessed using flow cytometry. Re-expression of miR-200c/141 or miR-200b/200a/429 in murine tumor cells significantly reduced proliferation (Fig. 3a). In the MDA-MB-231 cells re-expression of the miR-200c/141 cluster but not the miR-200b/200a/429 cluster significantly reduced proliferation in MDA-MB-231 cells (Fig. 3b).

Fig. 3
figure 3

Cell proliferation and invasion. Proliferation (a, b) and transwell migration across Matrigel coated inserts (c–l) in murine (a, c–g) and human (b, h–l) mammary tumor cells. Proliferation was determined using BrdU incorporation while transwell migration was determined by calculating the percentage of the lower transwell membrane covered with crystal violet stained cells. Representative images of the crystal violet stained cells on the bottom of the transwell insert have been provided for murine (c–f) and human (h–k) mammary tumor cells. Graphs (a, b, g, l) represent mean ± SEM with n = 3 for the BrdU assay and n = 4 for the transwell assay. * Indicate values significantly different (p < 0.05) from either RJ423EV or MDA-231EV cells

Full size image

To determine the impact of miR-200s on migration, matrigel-coated, transwell assays were performed. In RJ423 cells, re-expression of either the miR-200c/141 cluster or the miR-200b/200a/429 cluster significantly reduced cell migration and invasion however, these cells still migrated and invaded more efficiently than the RJ345 cells (Fig. 3c–g). Re-expression of the miR-200c/141 cluster but not the miR-200b/200a/429 cluster in MDA-MB-231 cells significantly reduced migration and invasion and re-expression of the miR-200c/141 cluster returned migration and invasion capacity to a level similar to MCF-7 cells (Fig. 3h–l).

Re-expression of miR-200s influence the expression of other miRNAs

miRNA sequencing was performed to confirm that the appropriate miR-200 family members had been upregulated in RJ423c141, RJ423ba429, MDA-231c141 and MDA-231ba429 cells and to determine whether re-expression of miR-200s influenced the expression of other miRNAs. The complete miRNA sequencing files for each cell line compared to their respective control line are presented as Additional file 1 and the top 10 upregulated and downregulated miRNAs based on p-values are presented in Table 1. In each cell line, the miR-200s expressed by the transfected plasmid were found within the top 10 upregulated miRNAs confirming appropriate miR-200 expression. In contrast to the qPCR data (Fig. 1a, c), miRNA-seq found a significant increase in miR-200a and miR-429 in the MDA-231ba429 cells compared to the MDA-231EV cells. The miRNA sequencing data also found a significant increase in miR-200c in RJ423ba429 cells and significant increase in miR-200b and miR-429 in MDA-231c141 cells (Additional file 1).

Table 1 Top 10 upregulated and downregulated miRNAs as determined by miRNA-seq
Full size table

To further explore miRNAs that were unique and shared in RJ423ba429 cells compared to RJ423c141 cells and MDA-231c141 cells compared to MDA-231ba429 cells, Venn diagrams were created (Additional file 2). These diagrams show miRNAs that were significantly, differentially expressed (log2FC ≥ 1, p < 0.01) in RJ423c141 and RJ423ba429 cells compared to RJ423EV cells (Additional file 2A) and MDA-231c131 and MDA-231ba429 cells compared to MDA-231EV cells (Additional file 2B). Re-expression of the miR-200b/200a/429 cluster altered the expression of more miRNAs than re-expression of the miR-200c/141 cluster in RJ423 cells and only a small number of miRNAs were shared (Additional file 2A). In the human cancer cells, re-expression of the miR-200c/141 cluster altered more miRNAs than re-expression of the miR-200b/200a/429 cluster and like the murine cells, a relatively small proportion of miRNAs were shared (Additional file 2B). Additional file 2C illustrates the number shared and unique miRNAs in the mouse and human miR-200 modified cells.

Re-expression of miR-200s altered the expression of genes regulated by SUZ12

RNA sequencing was also performed to identify genes altered by miR-200s. Hierarchical clustering revealed that the cell types segregated into discrete clusters and RJ423c141 cells clustered more closely with RJ423EV cells than RJ423ba429 cells (Fig. 4a). The analysis of this data was complicated by the fact that the RNA sequencing for the RJ423c141 and RJ423ba429 cells was performed on different days. In each RNA sequencing run, the same RNA from RJ423EV cells was used. However, the hierarchical clustering reveals some batch effects from the different runs as the two sets of RJ423EV samples (RJ423EVB samples that were utilized with the RJ423ba49 cells and RJ423EVC samples that were utilized with the RJ423c141 cells) showed distinct RNA expression profiles (Fig. 4b). Using a of log2FC ≥ 1 and FDR ≤ 0.01, 494 genes were upregulated and 1002 genes were downregulated in RJ423c141 cells compared to RJ423EV cells and 186 (~ 12%) of these genes had predicted miR-200c or miR-141 binding sites (using miRDB, Fig. 4b). When RJ423ba429 cells were compared to RJ423EV cells, 1691 genes were upregulated and 1899 were downregulated and 437 (~ 12%) of these genes had predicted miR-200b, 200a or miR-429 binding sites (Fig. 4b). Of the genes differentially regulated in RJ423c141 vs. RJ423EV and RJ423ba429 vs. RJ423EV, 913 were shared and 583 and 2677 were unique to RJ423c141 and RJ423ba429 cells, respectively (Fig. 4b).

Fig. 4
figure 4

Hierarchical clustering of RNA sequencing data (n = 4). Hierarchical clustering (a, c) and Venn diagram (b, d) analysis of mRNAs differentially expressed in RJ423c141 and RJ423ba429 cells relative to RJ423EV cells (a, b) as well as in, MDA-231c141 and MDA-231ba429 cells relative to MDA-231EV cells (c, d). The Venn diagrams also indicate the percentage of differentially expressed genes with predicted miR-200 binding sites as determined by the miRDB database

Full size image

RNA sequencing was also performed on MDA-231EV, MDA-231c141 and MDA-231ba429 cells. Again, hierarchical clustering revealed that the 3 cell types segregated into discrete clusters and MDA-231b429 cells were more similar to MDA-231EV cells than MDA-231c141 cells (Fig. 4c). Using a of log2FC ≥ 1 and FDR ≤ 0.01, 1733 genes were upregulated, and 1663 genes were downregulated in MDA-231c141 cells compared to MDA-231EV cells and 409 (~ 12%) have predicted miR-200c or miR-141 binding sites (Fig. 3d). When MDA-231ba429 cells were compared to MDA-231EV cells, 506 genes were upregulated and 948 were downregulated and 159 (~ 11%) of these genes had predicted miR-200b, miR-200a or miR-429 binding sites (Fig. 4d). Of the genes differentially regulated in MDA-231c141 vs. MDA-231EV and MDA-231ba429 vs. MDA-231EV, 766 were shared and 2630 and 688 were unique to MDA-231c141 and MDA-231ba429 cells, respectively (Fig. 4d).

Since RJ423ba429 cells and MDA-231c141 cells display a change in cell morphology (Fig. 2), the genes unique to these two cell lines were further examined using Enrichr [40, 41] to identify pathways associated with the unique genes in RJ423ba429 and MDA-231c141 cells (Fig. 5a, b). Although there were only a small number of pathways shared between RJ423ba429 and MDA-231c141 cells, the top ChEA category was the same in both cell lines, SUZ12 (Fig. 5c). The Encode and ChEA Consensus TFs from ChIP-X database from Enrichr lists 4439 predicted Suz12 targets and it was found that 863 of the 2677 (~ 32%) unique genes in RJ423ba429 cells and 706 of the 2630 (~ 27%) unique genes in MDA-231c141 cells were predicted targets of SUZ12.

Fig. 5
figure 5

Pathway analysis of genes regulated by miR-200s. KEGG Pathways (a), Gene Ontology (b), and ENCODE/ChEA analysis (c) of genes differentially expressed in RJ423ba429 cells (white bars) or MDA-231c141 cells (grey bars)

Full size image

When comparing the RJ423ba429 unique gene profile to the MDA-231c141 unique gene profile, it was determined that 370 genes were shared between the unique gene profiles of the two cell lines (Fig. 6a). The KEGG and Gene Ontology Pathways as well as the ChEA pathways of these shared genes are presented in Fig. 5b, c and suggest that altering miR-200 expression in these cells influences pathways such as hippo signaling, endothelial migration and metallopeptidase activity (Fig. 6d). The top ChEA pathways remained SUZ12 and 194 of the 370 shared genes (~ 52%) were predicted SUZ12 targets.

Fig. 6
figure 6

Venn diagrams and pathway analysis of shared genes. Venn diagram (a), KEGG Pathways (b), Gene Ontology (c) and ENCODE/ChEA analysis (d) of genes differentially expressed in both RJ423ba429 and MDA-231c141 cells

Full size image

H3K27me3 levels are elevated in RJ423ba429 and MDA-231c141 cells

SUZ12 is a member of the polycomb repressor complex 2 (PRC2) which mono-, di- and tri-methylates lysine 27 on histone H3 [43,44,45,46,47,48,49,50]. To further investigate differences in SUZ12 activity, flow cytometry was used to quantify trimethylated histone H3 (H3K27me3). As shown in Fig. 7, RJ423ba429 cells had significantly higher levels of H3K27me3 than RJ423EV cells (Fig. 7a) and MDA-231c141 cells had significantly higher levels of H3K27me3 than MDA-231EV cells (Fig. 7b). The luminal breast cancer cell line, MCF-7, also had significantly elevated levels of H3K27me3 compared to the claudin-low breast cancer cell line MDA-231EV while RJ345 cells that have both luminal and basal-like characteristics, did not have higher levels of H3K27me3 compared to the claudin-low RJ423EV cells.

Fig. 7
figure 7

Quantification of H3K27me3 levels using flow cytometry. The levels of H3K27me3 in RJ345, RJ423c141, and RJ423ba429 cells were expressed relative to the levels of H3K27me3 in RJ423EV cells (a) while the levels of H3K27me3 in MCF-7, MDA-231c141 and MDA-231ba429 cells were expressed relative to the levels of H3K27me3 in MDA-231EV cells (b). The graphs represent the mean ± SEM (n = 5) and the * indicates p < 0.05 relative to RJ423EV cells (a) and MDA-231EV cells (b)

Full size image

Discussion

Re-expression of miR-200s in human and murine claudin-low mammary tumor cells is sufficient to induce a partial reversion to an epithelial morphology and to reduce cell migration. The miR-200b/200a/429 cluster but not the miR-200c/141 cluster was able to induce these alterations in cell shape and migration in the murine mammary tumor cell line, RJ423, while it was the miR-200c/141 cluster, and not the miR-200b/200a/429 cluster, that was capable of inducing these alterations in MDA-MB-231 cells. While cluster-specific effects in murine and human cell lines are possible, it is also possible that miR-200 expression levels are an important factor. In RJ423 cells, the miR-200b/200a/429 cluster was expressed at higher levels than the miR-200c/141 cluster while in the MDA-MB-231 cells the miR-200c/141 cluster was expressed at higher levels than the miR-200b/200a/429 cluster. Therefore, miR-200s may require expression levels beyond a certain threshold in order to have significant cellular effects. It also appears that if one miR-200 cluster is expressed at sufficiently high levels, it can positively influence the expression from the other miR-200 cluster. This was demonstrated in both RJ423ba429 where re-expression of the miR-200b/200a/429 cluster also induced a significant increase in miR-200c while in MDA-231c141 cells re-expression of the miR-200c/141 cluster induced a significant increase in miR-200b and miR-429 expression as determined by miRNA sequencing.

While the impact of miRNAs regulating the expression of mRNAs has been extensively investigated, the regulation of miRNAs by other miRNAs has largely been ignored. A small number of papers have shown that miRNAs can regulate other miRNAs [51,52,53], however, no paper describing a comprehensive examination of miRNAs regulated by the miR-200 family could be found. The miRNA sequencing results from our study demonstrate that simply re-expressing miR-200s induced alterations in the expression of a number of other miRNAs. Not only was the number of miRNAs regulated by miR-200 re-expression unexpected (> 100 miRNAs regulated in RJ423ba429 and MDA-231c141 cells), so was the small overlap (16–18%) of miRNAs regulated by the miR-200c/141 and miR-200b/200a/429 clusters in both the murine and human cancer cell lines. The small overlap in expressed miRNAs suggest that the miR-200c/141 and miR-200b/200a/429 clusters regulate the expression of different miRNAs.

Our observations that re-rexpression of miR-200s in claudin-low mammary tumor cells can inhibit proliferation and migration/invasion as well as altering cell shape is consistent with several published studies. These studies showed that re-expression of miR-200c in MDA-MB-231 cells significantly suppressed cell migration and promoted a more rounded cell morphology [54, 55], while treatment of MDA-231 cells with a miR-200c mimic was shown to significantly inhibit transwell invasion [56,57,58,59] or proliferation [57, 60]. Although our current study did not find that re-expression of the miR-200b/200a/429 cluster in MDA-MB-231 cells significantly impacted proliferation or invasion, other groups have shown that treating MDA-MB-231 cells with a miR-200a mimic [61], a miR-200b mimic [62,63,64,65], a miR-429 mimic [63], or overexpressing miR-200b [66], significantly reduced migration/invasion and treating MDA-MB-231 cells with a miR-200b or miR-429 mimic significantly reduced proliferation [63, 64]. It should be noted that two studies found miR-200 effects that contracted our findings. The first study stably expressed either the miR-200b/200a/429 or miR-200c/141 clusters in MDA-MB-231 cells and found that both clusters reduced cell growth (the miR-200c/141 cluster was more effective than the miR-200b/200a/429 cluster) and both clusters promoted migration and invasion [67]. The second study showed that re-expression of miR-200c significantly promoted invasion in a transwell assay [68]. Studies on murine claudin-low mammary tumor cells are far rarer and the only study outside our lab on murine claudin-low mammary tumor cells was performed in vivo and found that re-expression of the miR-200c/141 cluster in p53null mammary tumors reduced tumor growth and metastasis to the lungs [69].

RNA sequencing revealed a similar pattern with RJ423ba429 and MDA-231c141 cells expressing more differentially expressed mRNAs than RJ423c141 or MDA-231ba429 cells, respectively. It also revealed that the miR-200c/141 and miR-200b/200a/429 clusters targeted distinct mRNAs despite both clusters sharing identical seed sequences. Cluster specific mRNA (and miRNA) expression profiles were observed in both murine and human mammary tumor cell lines reducing the possibility of this being an artifact. Another interesting observation from the RNA sequencing study was that no more that 12% of the differentially expressed mRNAs in any of the cell lines were predicted targets of the miR-200 family. This is not completely unexpected given the large number of miRNAs also regulated by the miR-200 family as described above. However, this data suggests that studies focusing only on predicted mRNA targets may miss a majority of mRNAs regulated by their miRNA of interest .

To reduce the complexity, differentially expressed mRNAs unique to RJ423ba429 and MDA-231c141 cells (the two lines showing changes in cell shape) were further examined. In both cell lines, the top ChEA pathway identified was SUZ12. SUZ12 is a component of the Polycomb repressive complex 2 (PRC2) that mono-, di- and trimethylates lysine 27 of histone H3 (H3K27) [43,44,45,46,47,48,49,50]. Trimethylation of H3K27 promotes chromosome compaction preventing transcription factors from accessing transcription start sites in promoters [70,71,72]. miR-200b and miR-200c have been reported to target SUZ12 [73,74,75], downregulating its expression and leading to hypomethylation of H3K27 and loss of H3K27 trimethylation. However, re-expression of the miR-200b/200a/429 cluster in RJ423 cells and re-expression of the miR-200c/141 cluster in MDA-MB-231 cells resulted in a significant increase in H3K27me3 levels. In normal mammary epithelial cells, higher levels of H3K27me3 are observed in mature luminal cells compared to mammary stem cells and H3K27me3 levels have been observed to increase during mammary epithelial differentiation during pregnancy [76]. Therefore, re-expression of the miR-200b/200a/429 cluster in RJ423 cells and re-expression of miR-200c/141 cluster in MDA-MB-231 cells may reduce the stem cell features of these cells and drive them to a more differentiated phenotype. Evaluation of gene specific H3K27 methylation patterns are required to confirm this and identify the genes affected.

Conclusions

Understanding how miRNAs influence cellular function is already complex due to the large number of predicted targets and false positives in existing databases, however, our study demonstrates that miR-200s can regulate other microRNAs as well as epigenetic modulators thus adding additional complexity to this system. We have also shown that re-expression of miR-200s in claudin-low mammary tumor cells can revert claudin-low cells to a more epithelial phenotype provided the level of miR-200 re-expression is sufficiently high and this reversion is associated with alterations in SUZ12 regulated genes and reduced cell migration.

Availability of data and materials

RNA sequencing has been uploaded to GEO under accession number GSE150107 and the miRNA sequencing files are presented as Additional file 1.

Abbreviations

APC:

Allophycocyanin

BrdU:

5-Bromo-2′-deoxyuridine

EMT:

Epithelial to mesenchymal transition

ER:

Estrogen receptor

FITC:

Fluorescein isothiocyanate

H3K27:

Lysine 27 of histone H3

H3K27me3:

Trimethylated lysine 27 of histone H3

HER2:

Human epidermal growth factor receptor 2

MDA-231ba429:

MDA-MB-231 cells stably expressing the miR-200b/200a/429 cluster

MDA-231c141:

MDA-MB-231 cells stably expressing the miR-200c/141 cluster

MDA-231EV:

MDA-MB-231 cells stably expressing an empty vector

miRNA:

microRNA

PCR2:

Polycomb repressor complex 2

PR:

Progesterone receptor

RJ423ba429:

RJ423 cells stably expressing the miR-200b/200a/429 cluster

RJ423c141:

RJ423 cells stably expressing the miR-200c/141 cluster

RJ423EV:

RJ423 cells stably expressing an empty vector

TNBC:

Triple negative breast cancer

References

  1. Ha M, Kim VN. Regulation of microRNA biogenesis. Nat Rev Mol Cell Biol. 2014;15(8):509–24.

    Article  CAS  PubMed  Google Scholar 

  2. Ambros V. The functions of animal microRNAs. Nature. 2004;431(7006):350–5.

    Article  CAS  PubMed  Google Scholar 

  3. Mei J, Hao L, Wang H, Xu R, Liu Y, Zhu Y, Liu C. Systematic characterization of non-coding RNAs in triple-negative breast cancer. Cell Prolif. 2020;53:e12801.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116(2):281–97.

    Article  CAS  PubMed  Google Scholar 

  5. Di LG, Calin GA, Croce CM. MicroRNAs: fundamental facts and involvement in human diseases. Birth Defects Res C Embryo Today. 2006;78(2):180–9.

    Article  CAS  Google Scholar 

  6. Czech B, Hannon GJ. Small RNA sorting: matchmaking for Argonautes. Nat Rev Genet. 2011;12(1):19–31.

    Article  CAS  PubMed  Google Scholar 

  7. Hock J, Meister G. The Argonaute protein family. Genome Biol. 2008;9(2):210.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  8. Filipowicz W, Bhattacharyya SN, Sonenberg N. Mechanisms of post-transcriptional regulation by microRNAs: are the answers in sight? Nat Rev Genet. 2008;9(2):102–14.

    Article  CAS  PubMed  Google Scholar 

  9. Bartel DP. MicroRNAs: target recognition and regulatory functions. Cell. 2009;136(2):215–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Lewis BP, Shih IH, Jones-Rhoades MW, Bartel DP, Burge CB. Prediction of mammalian microRNA targets. Cell. 2003;115(7):787–98.

    Article  CAS  PubMed  Google Scholar 

  11. Brennecke J, Stark A, Russell RB, Cohen SM. Principles of microRNA-target recognition. PLoS Biol. 2005;3(3):e85.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  12. Jonas S, Izaurralde E. Towards a molecular understanding of microRNA-mediated gene silencing. Nat Rev Genet. 2015;16(7):421–33.

    Article  CAS  PubMed  Google Scholar 

  13. McGeary SE, Lin KS, Shi CY, Pham TM, Bisaria N, Kelley GM, Bartel DP. The biochemical basis of microRNA targeting efficacy. Science. 2019;366(6472):eaav1741.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Grimson A, Farh KK, Johnston WK, Garrett-Engele P, Lim LP, Bartel DP. MicroRNA targeting specificity in mammals: determinants beyond seed pairing. Mol Cell. 2007;27(1):91–105.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Matsui M, Chu Y, Zhang H, Gagnon KT, Shaikh S, Kuchimanchi S, Manoharan M, Corey DR, Janowski BA. Promoter RNA links transcriptional regulation of inflammatory pathway genes. Nucleic Acids Res. 2013;41(22):10086–109.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Dragomir MP, Knutsen E, Calin GA. SnapShot: unconventional miRNA functions. Cell. 2018;174(4):1038.

    Article  CAS  PubMed  Google Scholar 

  17. Saini HK, Enright AJ, Griffiths-Jones S. Annotation of mammalian primary microRNAs. BMC Genom. 2008;9:564.

    Article  CAS  Google Scholar 

  18. Park SM, Gaur AB, Lengyel E, Peter ME. The miR-200 family determines the epithelial phenotype of cancer cells by targeting the E-cadherin repressors ZEB1 and ZEB2. Genes Dev. 2008;22(7):894–907.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Hill L, Browne G, Tulchinsky E. ZEB/miR-200 feedback loop: at the crossroads of signal transduction in cancer. Int J Cancer. 2013;134(4):745–54.

    Article  CAS  Google Scholar 

  20. Trumbach D, Prakash N. The conserved miR-8/miR-200 microRNA family and their role in invertebrate and vertebrate neurogenesis. Cell Tissue Res. 2015;359(1):161–77.

    Article  PubMed  CAS  Google Scholar 

  21. Howe EN, Cochrane DR, Richer JK. The miR-200 and miR-221/222 microRNA families: opposing effects on epithelial identity. J Mammary Gland Biol Neoplasia. 2012;17(1):65–77.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Radisky DC. miR-200c at the nexus of epithelial–mesenchymal transition, resistance to apoptosis, and the breast cancer stem cell phenotype. Breast Cancer Res. 2011;13(3):110.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Bracken CP, Gregory PA, Kolesnikoff N, Bert AG, Wang J, Shannon MF, Goodall GJ. A double-negative feedback loop between ZEB1-SIP1 and the microRNA-200 family regulates epithelial-–mesenchymal transition. Cancer Res. 2008;68(19):7846–54.

    Article  CAS  PubMed  Google Scholar 

  24. Peinado H, Olmeda D, Cano A. Snail, Zeb and bHLH factors in tumour progression: an alliance against the epithelial phenotype? Nat Rev Cancer. 2007;7(6):415–28.

    Article  CAS  PubMed  Google Scholar 

  25. Yang J, Weinberg RA. Epithelial–mesenchymal transition: at the crossroads of development and tumor metastasis. Dev Cell. 2008;14(6):818–29.

    Article  CAS  PubMed  Google Scholar 

  26. Bockmeyer CL, Christgen M, Muller M, Fischer S, Ahrens P, Langer F, Kreipe H, Lehmann U. MicroRNA profiles of healthy basal and luminal mammary epithelial cells are distinct and reflected in different breast cancer subtypes. Breast Cancer Res Treat. 2011;130(3):735–45.

    Article  CAS  PubMed  Google Scholar 

  27. Cochrane DR, Cittelly DM, Howe EN, Spoelstra NS, McKinsey EL, LaPara K, Elias A, Yee D, Richer JK. MicroRNAs link estrogen receptor alpha status and Dicer levels in breast cancer. Horm Cancer. 2010;1(6):306–19.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Korpal M, Lee ES, Hu G, Kang Y. The miR-200 family inhibits epithelial–mesenchymal transition and cancer cell migration by direct targeting of E-cadherin transcriptional repressors ZEB1 and ZEB2. J Biol Chem. 2008;283(22):14910–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Mekala JR, Naushad SM, Ponnusamy L, Arivazhagan G, Sakthiprasad V, Pal-Bhadra M. Epigenetic regulation of miR-200 as the potential strategy for the therapy against triple-negative breast cancer. Gene. 2018;641:248–58.

    Article  CAS  PubMed  Google Scholar 

  30. Humphries B, Wang Z, Oom AL, Fisher T, Tan D, Cui Y, Jiang Y, Yang C. MicroRNA-200b targets protein kinase Calpha and suppresses triple-negative breast cancer metastasis. Carcinogenesis. 2014;35(10):2254–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Sorlie T, Perou CM, Tibshirani R, Aas T, Geisler S, Johnsen H, Hastie T, Eisen MB, van de Jeffrey RM, et al. Gene expression patterns of breast carcinomas distinguish tumor subclasses with clinical implications. Proc Natl Acad Sci USA. 2001;98(19):10869–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Perou CM, Sorlie T, Eisen MB, van de Jeffrey RM, Rees SS, Pollack CA, Ross JR, Johnsen DT, Akslen H, et al. Molecular portraits of human breast tumours. Nature. 2000;406(6797):747–52.

    Article  CAS  PubMed  Google Scholar 

  33. Banerjee S, Reis-Filho JS, Ashley S, Steele D, Ashworth A, Lakhani SR, Smith IE. Basal-like breast carcinomas: clinical outcome and response to chemotherapy. J Clin Pathol. 2006;59(7):729–35.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Dent R, Trudeau M, Pritchard KI, Hanna WM, Kahn HK, Sawka CA, Lickley LA, Rawlinson E, Sun P, Narod SA. Triple-negative breast cancer: clinical features and patterns of recurrence. Clin Cancer Res. 2007;13(15 Pt 1):4429–34.

    Article  PubMed  Google Scholar 

  35. Prat A, Parker JS, Karginova O, Fan C, Livasy C, Herschkowitz JI, He X, Perou CM. Phenotypic and molecular characterization of the claudin-low intrinsic subtype of breast cancer. Breast Cancer Res. 2010;12(5):R68.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  36. Watson KL, Jones RA, Bruce A, Moorehead RA. The miR-200b/200a/429 cluster prevents metastasis and induces dormancy in a murine claudin-low mammary tumor cell line. Exp Cell Res. 2018;369(1):17–26.

    Article  CAS  PubMed  Google Scholar 

  37. Jones R, Watson K, Bruce A, Nersesian S, Kitz J, Moorehead R. Re-expression of miR-200c suppresses proliferation, colony formation and in vivo tumor growth of murine claudin-low mammary tumor cells. Oncotarget. 2017;8(14):23727–49.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Oliveros JC. VENNY. An interactive tool for comparing lists with Venn Diagrams. 2007.

  39. Chen EY, Tan CM, Kou Y, Duan Q, Wang Z, Meirelles GV, Clark NR, Ma’ayan A. Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinform. 2013;14:128.

    Article  Google Scholar 

  40. Kuleshov MV, Jones MR, Rouillard AD, Fernandez NF, Duan Q, Wang Z, Koplev S, Jenkins SL, Jagodnik KM, Lachmann A, et al. Enrichr: a comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res. 2016;44(W1):W90–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Chen Y, Wang X. miRDB: an online database for prediction of functional microRNA targets. Nucleic Acids Res. 2020;48(D1):D127–31.

    Article  CAS  PubMed  Google Scholar 

  42. Liu W, Wang X. Prediction of functional microRNA targets by integrative modeling of microRNA binding and target expression data. Genome Biol. 2019;20(1):18.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Martin C, Moorehead RA. Polycomb repressor complex 2 function in breast cancer. Int J Oncol. 2020;57(5):1085–94.

    CAS  PubMed Central  PubMed  Google Scholar 

  44. Vizan P, Beringer M, Ballare C, Di Croce L. Role of PRC2-associated factors in stem cells and disease. FEBS J. 2015;282(9):1723–35.

    Article  CAS  PubMed  Google Scholar 

  45. Thornton SR, Butty VL, Levine SS, Boyer LA. Polycomb repressive complex 2 regulates lineage fidelity during embryonic stem cell differentiation. PLoS ONE. 2014;9(10):e110498.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  46. Cholewa-Waclaw J, Bird A, von Schimmelmann M, Schaefer A, Yu H, Song H, Madabhushi R, Tsai LH. The role of epigenetic mechanisms in the regulation of gene expression in the nervous system. J Neurosci. 2016;36(45):11427–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Moritz LE, Trievel RC. Structure, mechanism, and regulation of polycomb-repressive complex 2. J Biol Chem. 2018;293(36):13805–14.

    Article  CAS  PubMed  Google Scholar 

  48. Hyun K, Jeon J, Park K, Kim J. Writing, erasing and reading histone lysine methylations. Exp Mol Med. 2017;49(4):e324.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Holoch D, Margueron R. Mechanisms regulating PRC2 recruitment and enzymatic activity. Trends Biochem Sci. 2017;42(7):531–42.

    Article  CAS  PubMed  Google Scholar 

  50. Laugesen A, Hojfeldt JW, Helin K. Molecular mechanisms directing PRC2 recruitment and H3K27 methylation. Mol Cell. 2019;74(1):8–18.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Kim JS, Kim EJ, Lee S, Tan X, Liu X, Park S, Kang K, Yoon JS, Ko YH, Kurie JM, et al. MiR-34a and miR-34b/c have distinct effects on the suppression of lung adenocarcinomas. Exp Mol Med. 2019;51(1):1–10.

    PubMed  PubMed Central  Google Scholar 

  52. Tang R, Li L, Zhu D, Hou D, Cao T, Gu H, Zhang J, Chen J, Zhang CY, Zen K. Mouse miRNA-709 directly regulates miRNA-15a/16-1 biogenesis at the posttranscriptional level in the nucleus: evidence for a microRNA hierarchy system. Cell Res. 2012;22(3):504–15.

    Article  CAS  PubMed  Google Scholar 

  53. Wang D, Sun X, Wei Y, Liang H, Yuan M, Jin F, Chen X, Liu Y, Zhang CY, Li L, et al. Nuclear miR-122 directly regulates the biogenesis of cell survival oncomiR miR-21 at the posttranscriptional level. Nucleic Acids Res. 2018;46(4):2012–29.

    Article  CAS  PubMed  Google Scholar 

  54. Ljepoja B, Schreiber C, Gegenfurtner FA, Garcia-Roman J, Kohler B, Zahler S, Radler JO, Wagner E, Roidl A. Inducible microRNA-200c decreases motility of breast cancer cells and reduces filamin A. PLoS ONE. 2019;14(11):e0224314.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Hurteau GJ, Carlson JA, Spivack SD, Brock GJ. Overexpression of the microRNA hsa-miR-200c leads to reduced expression of transcription factor 8 and increased expression of E-cadherin. Cancer Res. 2007;67(17):7972–6.

    Article  CAS  PubMed  Google Scholar 

  56. Zhang T, Wan JG, Liu JB, Deng M. MiR-200c inhibits metastasis of breast tumor via the downregulation of Foxf2. Genet Mol Res. 2017. https://doi.org/10.4238/gmr16038971.

    Article  PubMed  Google Scholar 

  57. Zhang J, Li G, Chen Y, Fang L, Guan C, Bai F, Ma M, Lyu J, Meng QH. Metformin inhibits tumorigenesis and tumor growth of breast cancer cells by upregulating miR-200c but downregulating AKT2 expression. J Cancer. 2017;8(10):1849–64.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  58. Chang BP, Wang DS, Xing JW, Yang SH, Chu Q, Yu SY. miR-200c inhibits metastasis of breast cancer cells by targeting HMGB1. J Huazhong Univ Sci Technol Med Sci. 2014;34(2):201–6.

    Article  CAS  Google Scholar 

  59. Luo D, Wilson JM, Harvel N, Liu J, Pei L, Huang S, Hawthorn L, Shi H. A systematic evaluation of miRNA: mRNA interactions involved in the migration and invasion of breast cancer cells. J Transl Med. 2013;11:57.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Ren Y, Han X, Yu K, Sun S, Zhen L, Li Z, Wang S. microRNA-200c downregulates XIAP expression to suppress proliferation and promote apoptosis of triple-negative breast cancer cells. Mol Med Rep. 2014;10(1):315–21.

    Article  CAS  PubMed  Google Scholar 

  61. Kim HK, Park JD, Choi SH, Shin DJ, Hwang S, Jung HY, Park KS. Functional link between miR-200a and ELK3 regulates the metastatic nature of breast cancer. Cancers. 2020;12(5):1225.

    Article  CAS  PubMed Central  Google Scholar 

  62. Yuan J, Xiao C, Lu H, Yu H, Hong H, Guo C, Wu Z. miR-200b regulates breast cancer cell proliferation and invasion by targeting radixin. Exp Ther Med. 2020;19(4):2741–50.

    CAS  PubMed  PubMed Central  Google Scholar 

  63. Li D, Wang H, Song H, Xu H, Zhao B, Wu C, Hu J, Wu T, Xie D, Zhao J, et al. The microRNAs miR-200b-3p and miR-429-5p target the LIMK1/CFL1 pathway to inhibit growth and motility of breast cancer cells. Oncotarget. 2017;8(49):85276–89.

    Article  PubMed  PubMed Central  Google Scholar 

  64. Zheng Q, Cui X, Zhang D, Yang Y, Yan X, Liu M, Niang B, Aziz F, Liu S, Yan Q, et al. miR-200b inhibits proliferation and metastasis of breast cancer by targeting fucosyltransferase IV and alpha1,3-fucosylated glycans. Oncogenesis. 2017;6(7):e358.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Hong H, Yu H, Yuan J, Guo C, Cao H, Li W, Xiao C. MicroRNA-200b impacts breast cancer cell migration and invasion by regulating Ezrin-Radixin-Moesin. Med Sci Monit. 2016;22:1946–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Pang Y, Liu J, Li X, Xiao G, Wang H, Yang G, Li Y, Tang SC, Qin S, Du N, et al. MYC and DNMT3A-mediated DNA methylation represses microRNA-200b in triple negative breast cancer. J Cell Mol Med. 2018;22(12):6262–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  67. Choi SK, Kim HS, Jin T, Hwang EH, Jung M, Moon WK. Overexpression of the miR-141/200c cluster promotes the migratory and invasive ability of triple-negative breast cancer cells through the activation of the FAK and PI3K/AKT signaling pathways by secreting VEGF-A. BMC Cancer. 2016;16(1):570.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  68. Meng Z, Zhang R, Wang Y, Zhu G, Jin T, Li C, Zhang S. miR-200c/PAI-2 promotes the progression of triple negative breast cancer via M1/M2 polarization induction of macrophage. Int Immunopharmacol. 2020;81:106028.

    Article  CAS  PubMed  Google Scholar 

  69. Knezevic J, Pfefferle AD, Petrovic I, Greene SB, Perou CM, Rosen JM. Expression of miR-200c in claudin-low breast cancer alters stem cell functionality, enhances chemosensitivity and reduces metastatic potential. Oncogene. 2015;34(49):5997–6006.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Aranda S, Mas G, Di Croce L. Regulation of gene transcription by Polycomb proteins. Sci Adv. 2015;1(11):e1500737.

    Article  PubMed  PubMed Central  Google Scholar 

  71. Kouznetsova VL, Tchekanov A, Li X, Yan X, Tsigelny IF. Polycomb repressive 2 complex-molecular mechanisms of function. Protein Sci. 2019;28(8):1387–99.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  72. Bracken AP, Dietrich N, Pasini D, Hansen KH, Helin K. Genome-wide mapping of Polycomb target genes unravels their roles in cell fate transitions. Genes Dev. 2006;20(9):1123–36.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Peng F, Jiang J, Yu Y, Tian R, Guo X, Li X, Shen M, Xu M, Zhu F, Shi C, et al. Direct targeting of SUZ12/ROCK2 by miR-200b/c inhibits cholangiocarcinoma tumourigenesis and metastasis. Br J Cancer. 2013;109(12):3092–104.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  74. Iliopoulos D, Lindahl-Allen M, Polytarchou C, Hirsch HA, Tsichlis PN, Struhl K. Loss of miR-200&nbsp;inhibition of Suz12 leads to polycomb-mediated repression required for the formation and maintenance of cancer stem cells. Mol Cell. 2010;39(5):761–72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  75. Polytarchou C, Iliopoulos D, Struhl K. An integrated transcriptional regulatory circuit that reinforces the breast cancer stem cell state. Proc Natl Acad Sci USA. 2012;109(36):14470–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  76. Lim YY, Wright JA, Attema JL, Gregory PA, Bert AG, Smith E, Thomas D, Lopez AF, Drew PA, Khew-Goodall Y, et al. Epigenetic modulation of the miR-200 family is associated with transition to a breast cancer stem-cell-like state. J Cell Sci. 2013;126(Pt 10):2256–66.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was funded by a CIHR operating grant (MOP-136970) and a CHIR project grant (PJT-162218) awarded to RAM.

Author information

Authors and Affiliations

  1. Department of Biomedical Sciences, Ontario Veterinary College, University of Guelph, Guelph, ON, Canada

    K. Simpson, G. Conquer-van Heumen, K. L. Watson, M. Roth, C. J. Martin & R. A. Moorehead

Authors
  1. K. Simpson
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. G. Conquer-van Heumen
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. K. L. Watson
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. M. Roth
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. C. J. Martin
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. R. A. Moorehead
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

KS performed the crystal violet staining, transwell migration assay, flow cytometry for BrdU and H3K27me3 and qRT-PCR  for the human breast cancer cells as well as editing the manuscript. GC-VH performed the crystal violet staining, transwell migration assay and flow cytometry for BrdU and H3K27me3 for the murine mammary cells as well as editing the manuscript. CJM performed the qRT-PCR for the mRNAs in the control and miR-200 transfected murine cells. MR performed the qRT-PCR for the microRNAs in the control and miR-200 transfected murine cells. KLW assisted with the qRT-PCR for both the miRNAs and the mRNAs. RAM performed the analysis of the miRNA sequencing and mRNA sequencing data, wrote the manuscript and ran the project. All authors read and approved the final manuscript.

Corresponding author

Correspondence to R. A. Moorehead.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

miRNA sequencing files for RJ423c141 vs. RJ423EV, RJ423ba429 vs. RJ423EV, MDA-231c141 vs. MDA-231ba429 and MDA-231ba429 vs. MDA-231EV cell lines.

Additional file 2.

Venn diagrams showing unique and shared miRNAs in the mouse cell lines (A), the human cell lines (B), and a combination of the mouse and human cell lines (C).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Simpson, K., Conquer-van Heumen, G., Watson, K.L. et al. Re-expression of miR-200s in claudin‐low mammary tumor cells alters cell shape and reduces proliferation and invasion potentially through modulating other miRNAs and SUZ12 regulated genes. Cancer Cell Int 21, 89 (2021). https://doi.org/10.1186/s12935-021-01784-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12935-021-01784-4

Keywords

Cancer Cell International

ISSN: 1475-2867

Contact us