- Primary research
- Open Access
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)
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.
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.
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.
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.
MicroRNAs (miRNAs) are small, non-coding RNA molecules that can be found in both introns and exons . 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 . 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 . 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 . 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.
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 .
Crystal violet staining and transwell invasion assays
Crystal violet staining and transwell invasion assays were performed as described in  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 . 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 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 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 . 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  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 were generated using Venny 2.1.0 .
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.
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.
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.
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 ) 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.
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).
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).
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).
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.
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.
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.
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 , a miR-200b mimic [62,63,64,65], a miR-429 mimic , or overexpressing miR-200b , 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 . The second study showed that re-expression of miR-200c significantly promoted invasion in a transwell assay . 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 .
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 . 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.
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.
Epithelial to mesenchymal transition
Lysine 27 of histone H3
Trimethylated lysine 27 of histone H3
Human epidermal growth factor receptor 2
MDA-MB-231 cells stably expressing the miR-200b/200a/429 cluster
MDA-MB-231 cells stably expressing the miR-200c/141 cluster
MDA-MB-231 cells stably expressing an empty vector
Polycomb repressor complex 2
RJ423 cells stably expressing the miR-200b/200a/429 cluster
RJ423 cells stably expressing the miR-200c/141 cluster
RJ423 cells stably expressing an empty vector
Triple negative breast cancer
Ha M, Kim VN. Regulation of microRNA biogenesis. Nat Rev Mol Cell Biol. 2014;15(8):509–24.
Ambros V. The functions of animal microRNAs. Nature. 2004;431(7006):350–5.
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.
Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116(2):281–97.
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.
Czech B, Hannon GJ. Small RNA sorting: matchmaking for Argonautes. Nat Rev Genet. 2011;12(1):19–31.
Hock J, Meister G. The Argonaute protein family. Genome Biol. 2008;9(2):210.
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.
Bartel DP. MicroRNAs: target recognition and regulatory functions. Cell. 2009;136(2):215–33.
Lewis BP, Shih IH, Jones-Rhoades MW, Bartel DP, Burge CB. Prediction of mammalian microRNA targets. Cell. 2003;115(7):787–98.
Brennecke J, Stark A, Russell RB, Cohen SM. Principles of microRNA-target recognition. PLoS Biol. 2005;3(3):e85.
Jonas S, Izaurralde E. Towards a molecular understanding of microRNA-mediated gene silencing. Nat Rev Genet. 2015;16(7):421–33.
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.
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.
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.
Dragomir MP, Knutsen E, Calin GA. SnapShot: unconventional miRNA functions. Cell. 2018;174(4):1038.
Saini HK, Enright AJ, Griffiths-Jones S. Annotation of mammalian primary microRNAs. BMC Genom. 2008;9:564.
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.
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.
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.
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.
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.
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.
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.
Yang J, Weinberg RA. Epithelial–mesenchymal transition: at the crossroads of development and tumor metastasis. Dev Cell. 2008;14(6):818–29.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Oliveros JC. VENNY. An interactive tool for comparing lists with Venn Diagrams. 2007.
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.
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.
Chen Y, Wang X. miRDB: an online database for prediction of functional microRNA targets. Nucleic Acids Res. 2020;48(D1):D127–31.
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.
Martin C, Moorehead RA. Polycomb repressor complex 2 function in breast cancer. Int J Oncol. 2020;57(5):1085–94.
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.
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.
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.
Moritz LE, Trievel RC. Structure, mechanism, and regulation of polycomb-repressive complex 2. J Biol Chem. 2018;293(36):13805–14.
Hyun K, Jeon J, Park K, Kim J. Writing, erasing and reading histone lysine methylations. Exp Mol Med. 2017;49(4):e324.
Holoch D, Margueron R. Mechanisms regulating PRC2 recruitment and enzymatic activity. Trends Biochem Sci. 2017;42(7):531–42.
Laugesen A, Hojfeldt JW, Helin K. Molecular mechanisms directing PRC2 recruitment and H3K27 methylation. Mol Cell. 2019;74(1):8–18.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Aranda S, Mas G, Di Croce L. Regulation of gene transcription by Polycomb proteins. Sci Adv. 2015;1(11):e1500737.
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.
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.
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.
Iliopoulos D, Lindahl-Allen M, Polytarchou C, Hirsch HA, Tsichlis PN, Struhl K. Loss of miR-200 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.
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.
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.
This work was funded by a CIHR operating grant (MOP-136970) and a CHIR project grant (PJT-162218) awarded to RAM.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
miRNA sequencing files for RJ423c141 vs. RJ423EV, RJ423ba429 vs. RJ423EV, MDA-231c141 vs. MDA-231ba429 and MDA-231ba429 vs. MDA-231EV cell lines.
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).
About this article
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
- Claudin‐low breast cancer
- RNA sequencing
- miRNA sequencing