Next Article in Journal
MED12 Regulates Smooth Muscle Cell Functions and Participates in the Development of Aortic Dissection
Previous Article in Journal
Comparative Genomic Hybridization to Microarrays in Fetuses with High-Risk Prenatal Indications: Polish Experience with 7400 Pregnancies
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 13
Issue 4
10.3390/genes13040691
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Elevated Expression of miR-200c/141 in MDA-MB-231 Cells Suppresses MXRA8 Levels and Impairs Breast Cancer Growth and Metastasis In Vivo

by
Kaitlyn E. Simpson
,
Katrina L. Watson
and
Roger A. Moorehead
*
Department of Biomedical Sciences, Ontario Veterinary College, University of Guelph, Guelph, ON N1G 2W1, Canada
*
Author to whom correspondence should be addressed.
Genes 2022, 13(4), 691; https://0-doi-org.brum.beds.ac.uk/10.3390/genes13040691
Submission received: 9 March 2022 / Revised: 28 March 2022 / Accepted: 9 April 2022 / Published: 14 April 2022
(This article belongs to the Section Molecular Genetics and Genomics)
Browse Figures
Review Reports Versions Notes

Abstract

:
Breast cancer cells with mesenchymal characteristics, particularly the claudin-low subtype, express extremely low levels of miR-200s. Therefore, this study examined the functional impact of restoring miR-200 expression in a human claudin-low breast cancer cell line MDA-MB-231. MDA-MB-231 cells were stably transfected with a control vector (MDA-231EV) or the miR-200c/141 cluster (MDA-231c141). Injection of MDA-231c141 cells into the 4th mammary gland of NCG mice produced tumors that developed significantly slower than tumors produced by MDA-231EV cells. Spontaneous metastasis to the lungs was also significantly reduced in MDA-231c141 cells compared to MDA-231EV cells. RNA sequencing of MDA-231EV and MDA-231c141 tumors identified genes including MXRA8 as being downregulated in the MDA-231c141 tumors. MXRA8 was further investigated as elevated levels of MXRA8 were associated with reduced distant metastasis free survival in breast cancer patients. Quantitative RT-PCR and Western blotting confirmed that MXRA8 expression was significantly higher in mammary tumors induced by MDA-231EV cells compared to those induced by MDA-231c141 cells. In addition, MXRA8 protein was present at high levels in metastatic tumor cells found in the lungs. This is the first study to implicate MXRA8 in human breast cancer, and our data suggests that miR-200s inhibit growth and metastasis of claudin-low mammary tumor cells in vivo through downregulating MXRA8 expression.

1. Introduction

MicroRNAs (miRNAs) are small, non-coding RNAs that regulate mRNA expression by binding to mRNAs and in turn, prevent mRNA translation [1,2]. The interaction between miRNAs and mRNAs is primarily regulated by a region known as the seed sequence that binds to complementary sequences in the 3′ untranslated regions of mRNAs [3,4,5,6,7,8,9,10,11].
MicroRNAs have been implicated in regulating a number of cancers, including breast cancer. One family of miRNAs implicating in breast cancer is the miR-200 family. This family consists of five members organized into two clusters, the miR-200b/200a/429 cluster and the miR-200c/141 cluster [12,13,14]. These five miRNAs have highly similar seed sequences and the 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), which differs from the seed sequence of miR-200b, miR-200c and miR-429 by one nucleotide [15]. The most completely characterized function of the miR-200 family is their role in maintaining epithelial identity; miR-200 family members reduce the expression of mesenchymal transcription factors, such as ZEB1/2, TWIST1/2, SNAI1/2 [14,16,17,18], and increase the expression of epithelial genes such as E-cadherin [19,20].
Given their roles in epithelial identity, it is not surprising that miR-200s are expressed in human luminal breast cancer but are only expressed at very low levels in cancers with mesenchymal characteristics, such as claudin-low breast cancer [21,22,23,24,25]. Claudin-low breast cancers share features with triple negative breast cancer (TNBC), including low levels of the estrogen receptor (ER), progesterone receptor (PR) and human epidermal growth factor receptor 2 (HER2) [26]. Claudin-low breast cancers are differentiated from other TNBCs by the low levels of cell-cell adhesion and tight junction genes such as claudins 3, 4, 7 [26,27,28].
Our previous work investigating the function of miR-200s in human claudin-low breast cancer cell lines found that re-expression of the miR-200c/141 cluster in human MDA-MB-231 cells reverted these cells to a more epithelial morphology, significantly reduced proliferation, and impaired migration and invasion in vitro [29]. Thus, the goal of this study was to evaluate the impact of miR-200c/141 overexpression on in vivo mammary tumor growth and metastasis. In the current study, it was shown that re-expression of the miR-200c/141 cluster in MDA-MB-231 cells significantly impaired mammary tumor growth in vivo. Re-expression of the miR-200c/141 cluster in MDA-MB-231 cells was also capable of significantly inhibiting spontaneous metastasis to the lungs. The reduced levels of MXRA8 in the mammary tumors induced by MDA-231c141 cells implicate this gene in regulating mammary tumor growth and metastasis.

2. Materials and Methods

2.1. Cell Lines and Culture Conditions

MDA-231EV cells were created by infecting MDA-MB-231 cells with copGFP lentiviral particles (cat sc-108,084, Santa Cruz Biotechnology Inc., Dallas, TX, USA), and stable cells were selected by continual culture in 5 µg/mL puromycin (InvivoGen, San Diego, CA, USA). MDA-231c141 cells were created by transfecting MDA-MB-231 cells with pLenti 4.1 Ex miR200c-141 (cat #35,534, Addgene, Watertown, MA, USA) and stable cells were selected by continual culture in 5 µg/mL puromycin (InvivoGen, San Diego, CA, USA). The in vitro characterization of MDA-231EV and MDA-231c141 cells have previously been described [29].

2.2. Animals and Ethics

Mice were housed and cared for following guidelines established by the Central Animal Facility at the University of Guelph, and the guidelines established by the Canadian Council of Animal Care. This study was approved by the Animal Care Committee at the University of Guelph (AUP #3995).
MDA-231EV and MDA-231c141 cells were collected from logarithmically growing cultures, and 1 × 106 cells were injected into the mammary glands of female NCG (NOD-Prkdcem26Cd52Il2rgem26Cd22/NjuCrl) mice (Charles River, Wilmington MA, USA). Ten mice were injected with each cell line and injections occurred at three different times, separated by several weeks. For each independent round of injections, both MDA-231EV and MDA-231c141 cells were injected. All mice were monitored twice per week by palpating the mammary glands. Once a palpable mammary tumor was identified, tumor size was measured using digital calipers. MDA-231EV and MDA-231c141 tumor bearing mice were collected once the MDA-231EV mammary tumors reached ~10% of the mouse’s body weight (maximum size allowed by the Canadian Council on Animal Care).

2.3. Tumor Specific Growth Rates (SGR)

Tumor specific growth rates were calculated using the formula SGR = ln (V2/V1)/(t2 − t1), where V1 and V2 are the tumor volumes at the time of palpation (t1) and euthanasia (t2), respectively.

2.4. Histology and Immunohistochemistry

Mammary tumors and major organs were fixed in 10% formalin overnight and embedded in paraffin. Sections were cut and stained with hematoxylin and eosin for histologic analysis. Immunohistochemistry was performed as previously described [30], using a 1:200 dilution of an anti-vimentin antibody (ab16700, Abcam, Toronto, ON, Canada) or a 1:100 dilution of an anti-MXRA8 antibody (ab185444, Abcam). Tissue sections were scanned using a Motic Easyscan digital slide scanner (Motic, Richmond, BC, Canada).
Lung tumor burden was calculated using lung tissue sections from MDA-231EV and MDA-231c141 tumor bearing mice stained with the human anti-vimentin antibody. Comparison of the H&E stained sections with vimentin stained sections confirmed that metastatic cells stained positive for vimentin. Vimentin stained sections were analyzed using ImageJ software [31], and the area of the section occupied by vimentin positive cells was determined using the color threshold function of ImageJ software [31]. Lung area was determined using the freehand drawing tool to outline each piece of lung tissue and total lung area on each section was calculated using Image J software [31].

2.5. RNA Extraction and RNA Sequencing

RNA was isolated from four independent MDA-231EV and four independent MDA-231c141 tumors. All samples had a 260/280 ratio > 1.95 (Nanodrop) and ~2 μg of RNA was sent to Genome Quebec. Genome Quebec analyzed the RNA on an Agilent Bioanalyzer 2100 (Agilent, Santa Clara, CA, USA), and the samples had an RIN of 9.5 or greater. Libraries were generated from 250 ng of total RNA which included mRNA enrichment using the NEBNext Poly(A) Magnetic Isolation Module (New England BioLabs, Whitby, ON, Canada), cDNA synthesis using NEBNext RNA First Strand Synthesis and NEBNext Ultra Directional RNA Second Strand Synthesis Modules (New England BioLabs, Whitby, ON, Canada). The libraries were normalized, pooled and denatured in 0.05N NaOH and neutralized using HT1 buffer. The pool was loaded at 225 pM on an Illumina NovaSeq S4 lane using Xp protocol as per the manufacturer’s recommendations. The run was performed for 2 × 100 cycles (paired-end mode). A phiX library was used as a control and mixed with libraries at 1% level. Base calling was performed with RTA v3. Program bcl2fastq2 v2.20 was then used to demultiplex samples and generate fastq reads.
Fastq reads were analyzed using Genialis software (Genialis Inc., Houston, TX, USA) following the software’s standard RNA-seq pipeline, as previously described [29]. Hierarchical clustering was performed using Genialis software (Genialis Inc., Houston, TX, USA) with all genes and Pearson distance measure and average linkage clustering. Pathway analysis was performed using Enrichr software [32,33]. The RNA sequencing data has been uploaded to the GEO database as accession number GSE193479.

2.6. Real-Time PCR

For TaqMan real-time PCR, the same RNA that was extracted for sequencing was used. 100 ng of RNA was reversed transcribed using the TaqMan microRNA reverse transcription kit (Thermo Fisher Scientific, Burlington, ON, Canada) following the manufacturer’s protocol. Primers for reverse transcribing and amplifying cDNA for miR-141 (ID 000463), miR-200a (ID 000502), miR-200b (ID 001800), miR-200c (ID 002300), miR-429, RNU44 (ID 001094) and RNU48 (ID 001006) were obtained from Thermo Fisher Scientific (Burlington, ON, Canada). qPCR was performed using a CFX96 real-time PCR machine (Bio-Rad Laboratories, Mississauga, ON, Canada) using 2 μL of each cDNA reaction, 10 μL of TaqMan universal master mix II no UNG (Thermo Fisher Scientific, Burlington, ON, Canada), 1 μL of 20x primer, and 7 μL of water, under the following incubation conditions; 50 °C for 2 min, 95 °C for 10min, and then 40 cycles of 95 °C for 15 s and 60 °C for 1 min. miRNA levels were calculated using Bio-Rad CFX Manager 3.1 (Bio-Rad Laboratories, Mississauga, ON, Canada) and were normalized to the levels of RNU44 and RNU48.
For gene expression, 1 ug of total RNA (same RNA as used for sequencing) was reverse transcribed using 4 μL of qScript cDNA supermix (Quanta Biosciences, Beverly, MA, USA) per 20 μL reaction and the following incubation conditions on a thermal cycler; 25 °C for 5 min, 42 °C for 30 min, 85 °C for 5 min, and then held at 4 °C until used for qPCR. For qPCR, 2 μL of cDNA was mixed with 10 μL of SensiFast SYBR No-Rox (FroggaBio Inc, Concord, ON) and 7 μL of water and 1 µL of MXRA8 primer (qHsaCED0045992) or HPRT primer (qHsaCID0016375). Both primers were obtained from Bio-Rad Laboratories (Mississauga, ON, Canada). Gene expression was quantified using a CFX96 real-time PCR machine (Bio-Rad Laboratories, Mississauga, ON, Canada) and the following program; 95 °C for 2 min and then 40 cycles of 95 °C for 5 s and 60 °C for 30 s. HPRT was used as the housekeeping gene and gene expression was quantified using CFX Manager software (Bio-Rad Laboratories, Mississauga, ON, Canada).

2.7. Western Blotting

Western blotting was performed as previously described [34]. Briefly, protein was isolated from four independent MDA-231EV and MDA-231c141 tumors. Forty micrograms of protein for each sample were separated on an 8% polyacrylamide and transferred to nitrocellulose membrane. Membranes were incubated with a 1:500 dilution of an anti-MXRA8 antibody (ab185444; Abcam, Toronto, ON, Canada) overnight at 4 °C. Proteins were detected with HRP-linked anti-rabbit IgG secondary antibody (1:2000) (New England Biolabs, Ltd., Whitby, ON, Canada) and Clarity Western ECL substrate (Bio-Rad Canada, Mississauga, ON, Canada). The protein was then imaged using the ChemiDocTMXRS+ System (Bio-Rad Canada, Mississauga, ON, Canada), and quantification was performed using Image Lab software (Bio-Rad Canada, Mississauga, ON, Canada). Following detection of MXRA8, the membranes were blocked and incubated with a 1:5000 dilution of an anti-β-actin antibody (cat. no. 8457; New England Biolabs, Whitby, ON, Canada).

2.8. Kaplan–Meier Plotter

The online survival tool, Kaplan–Meier Plotter, was used to determine whether gene expression was associated with distant metastases-free survival (DMFS) [35]. This analysis was restricted to basal-like breast cancers as MDA-MB-231 cells represent a class of basal-like breast cancer. Patients were split by the median value [35].

2.9. Statistics

For analysis comparing the means of two different groups, a Student’s t-test was performed. Means were considered statistically different when p < 0.05.

3. Results

3.1. miR-200s Impair In Vivo Tumor Growth and Metastasis

Previously published work from our lab characterized the MDA-231EV and MDA-231c141 cells in vitro [29]. MDA-231EV cells were transfected with a control vector, while MDA-231c141 cells were transfected with a vector containing the miR-200c/141 cluster. To gain a better understanding of the function of miR-200s in tumorigenesis, MDA-231EV (1 × 106 cells) or MDA-231c141 (1 × 106 cells) were injected into the mammary fat pads of NCG (NOD-Prkdcem26Cd52Il2rgem26Cd22/NjuCrl) mice. As shown in Figure 1A, MDA-231c141 cells displayed reduced tumor growth compared to MDA-231EV cells (Figure 1A). Specific growth rates were calculated for each tumor, and the specific growth rate for MDA-231c141 tumors (1.48 ± 0.13) was significantly lower than MDA-231EV (3.24 ± 0.05) tumors.
Histological examination of the mammary tumors produced from intramammary injections of MDA-231EV or MDA-231c141 cells revealed that both tumors presented as solid sheets of tumor cells with variable amounts of stroma and necrosis (Figure 1B,C).
To assess whether miR-200 expression observed in the cell lines were maintained in vivo, Taqman RT-PCR was performed. As shown in Figure 1D–H, tumors produced by the injection of MDA-231c141 cells had significantly higher levels of miR-200c (Figure 1D) and miR-141 (Figure 1E) than MDA-231EV tumors.
MDA-231EV tumors spontaneously metastasized to the lungs in 100% of the mice (Figure 2A–C, Table 1), while MDA-231c141 tumors only metastasized to the lungs in 70% of the mice (Figure 2D–F, Table 1). Average tumor burden (based on positive staining for human vimentin) in the lungs of mice harboring MDA-231EV tumors were significantly higher than the average tumor burden of the lungs of mice harboring MDA-231c141 tumors (Figure 2G).
MDA-231EV tumors also frequently metastasized to the pancreas and spleen (Table 1, Supplemental File S1A–F) and less frequently to the liver (Table 1, Supplemental File S1G,H). The metastatic lesions that developed near the pancreas and spleen were unusual in that they appeared to grow just outside these tissues with sporadic regions of tumor cell invasion (Supplemental File S1A–F).

3.2. Elevated MXRA8 Is Associated with Tumor Growth, Metastasis, and Patient Prognosis

RNA sequencing was performed on tumors derived from MDA-231EV and MDA-231c141 cells. Hierarchical clustering of the tumors revealed that the MDA-231EV and MDA-231c141 tumors segregated into discrete clusters (Supplemental File S2). Using cut-offs of logFC ≥ 1 and FDR < 0.00001 it was determined that MDA-231c141 tumors had 1273 transcripts that were upregulated, and 1570 transcripts were downregulated compared to MDA-231EV tumors. Pathway analysis of the 2843 differentially expressed transcripts in the MDA-231c141 vs. MDA-231EV tumors using Enrichr software revealed the top ENCODE and ChEA Consensus TFS from ChIP-X was SUZ12, while the top MSigDB Hallmark pathway was epithelial mesenchymal transition, and the top GO Biological Process was axon guidance (Figure 3).
The top five differentially expressed transcripts in the MDA-231c141 vs. MDA-231EV comparison (based on FDR) are presented in Table 2, and all five were significantly downregulated in the MDA-231c141 tumors. Since all five transcripts were downregulated in the MDA-231c141 tumors, three databases (miRDB, microT-CDS, and miRWalk) were investigated to determine whether these transcripts contained potential miR-200c or miR-141 binding sites. The database miRDB uses the target prediction tool MirTarget [36,37] and microT-CDS identifies miRNA recognition elements in the 3′UTR and coding sequences [38], while miRWalk searches for potential binding sites in 3′-UTR, 5′-UTR and coding regions using TarPmiR [39,40]. As shown in Supplemental File Table S1, only miRwalk but not miRDB or MicroT-CDS indicated that these five transcripts had potential miR-200c or miR-141 binding sites.
To evaluate whether CPE, GNG2, MXRA8, SELENBP1 or AGR2 were associated with survival in breast cancer patients, these genes were examined using kmplot.com. Tumor types were restricted to the basal-like breast cancer subtype as the basal-like subtype is the closest subtype to claudin-low (MDA-MB-231 cells have characteristics of claudin-low breast cancer [26], and patients were split by the median expression values. Only elevated MXRA8 (Figure 4) was significantly associated with a decrease in distant metastasis free survival in basal-like breast cancer patients.
The levels of MXRA8 were then evaluated in the primary tumors using qRT-PCR (Figure 5A) and Western blotting (Figure 5B,C). At both the mRNA and protein level, MXRA8 was significantly higher in the MDA-231EV tumors compared to the MDA-231c141 tumors. To evaluate the levels of MXRA8 protein in lung metastasis, immunohistochemistry was performed. The metastatic tumor cells expressed high levels of MXRA8 compared to the adjacent normal lung tissue (Figure 6A,B).

4. Discussion

Although it is generally accepted that miR-200s negatively regulate breast cancer growth and metastasis, studies directly assessing miR-200c/miR-141 function in vivo are limited, and the findings are not always consistent. Only four studies have investigated the impact of miR-200c overexpression on in vivo growth of human breast cancer cell lines. All four studies found that elevated levels of miR-200s inhibited primary mammary tumor growth following intramammary [41,42,43] or subcutaneous injection [44]. Similarly, four of the five studies investigating miR-200c or miR-200c/141 overexpression in murine mammary tumor cells found that elevated miR-200s significantly reduce primary mammary tumor growth [44,45,46,47]. One study found that elevated expression of miR-200c and miR-141 in the murine mammary tumor cell line, 4TO7, did not significantly affect tumor growth [48]. The majority of these studies are consistent with our findings and suggest that miR-200s negatively regulate growth of primary mammary tumors.
Only one other study has evaluated the impact of miR-200c overexpression on spontaneous metastasis of human breast cancer cells and, like our study, showed that elevated miR-200c could significantly reduce spontaneous metastasis to the lungs [49]. Interestingly, two studies evaluating lung metastasis of MDA-MB-231 cells overexpressing miR-200c or miR-200c/141 following tail vein injection found that elevated miR-200 levels were associated with an increase in lung metastases [50,51]. There have also been two studies on metastasis of murine mammary tumor cells, with one study using an experimental metastasis assay [47], and the other evaluating spontaneous metastasis from a primary mammary tumor [46]. Both studies found that miR-200c/141 expression was negatively associated with lung metastases [46,47]. Since spontaneous metastasis assays more accurately reflect metastasis in the patient, our study and the studies by Liu et al. [46] and Li et al. [49] suggest that elevated miR-200c/141 expression will reduce metastatic dissemination.
Transcriptome profiling of MDA-231EV and MDA-231c141 tumors revealed a number of genes regulated by the miR-200c/141 cluster that could influence mammary tumor growth and metastasis including MXRA8. MXRA8 or matrix remodeling associated 8 is a transmembrane protein that can influence integrin signaling and can regulate cell–cell interactions. MXRA8 also serves as the receptor for arthritogenic alphaviruses [52]. While MXRA8 has been associated with esophageal, kidney, gingivobuccal cancer [53,54,55], there are no publications on MXRA8 in breast cancer. Therefore, our work demonstrates for the first time that MXRA8 is elevated in fast growing mammary tumors with a high metastatic potential. Our work also suggests that MXRA8 is associated with breast cancer metastasis, as lung metastases expressed high levels of MXRA8 and high MXRA8 was associated with poor distant metastasis free survival in patients with basal-like breast cancer. Further support demonstrating that MXRA8 plays a role in breast cancer stems from our work with murine mammary tumor cell lines. RJ423 cells, a murine mammary tumor cell line with claudin-low characteristics, express low levels of miR-200s and are highly tumorigenic and metastatic compared to RJ345 cells, which have some luminal characteristics [29,56,57,58]. A search of our RNA sequencing data revealed that RJ423 expressed significantly higher levels of MXRA8 (logFC 6.9, adjusted pval 5.4 × 10−66) than RJ345 cells. Thus, in two different cell lines, MXRA8 is elevated in aggressive, metastatic tumors expressing low miR-200s.
How MXRA8 is regulated by miR-200c/141 will require additional work. While MXRA8 has potential miR-200c or miR-141 binding sites, these sites were only predicted by miRWalk, and thus additional studies are required to determine whether MXRA8 is a direct miR-200 target. In addition, miR-200c and miR-141 may also regulate the expression of genes like MXRA8 indirectly through influencing the expression of transcription factors or altering the methylation status of DNA or histones. It is interesting to note that the top ENCODE and ChEA Consensus TFs from the ChIP-X pathway identified by Enrichr software was SUZ12. Our previously published transcriptome profiling of MDA-231c141 cells in 2D culture also found the top ENCODE and ChEA Consensus TFs from ChIP-X pathway to be SUZ12. SUZ12 is a member of the polycomb repressor complex 2 (PRC2) that mono-, di-, and tri-methylates histone H3 at lysine 27 (H3K27) [59,60]. Methylation of H3K27 is associated with decreased transcription. MXRA8 is a predicted target of SUZ12 (Gene Set-SUZ12 (maayanlab.cloud)), and thus suppression of MXRA8 expression in MDA-231c141 cells may be mediated by increased H3K27 methylation.
In summary, this study provides further evidence that elevated expression of the miR-200c/141 cluster is associated with reduced primary mammary tumor growth and metastasis. Moreover, this is the first study implicating MXRA8 as a potential regulator of mammary tumor growth and metastasis.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/genes13040691/s1, Supplemental File S1, Vimentin stained metastatic MDA-231EV cells in (A–D) pancreas, (E,F) spleen, and (G,H) liver. Metastatic tumor cells were frequently found near the (A,B) pancreas and (E,F) spleen and (C,D) only rarely invaded the pancreas. Scale bars for A,C,E,G are 600 µm and for B,D,F,H are 100 µm: Supplemental File S2, Hierarchical clustering of MDA-231EV, MDA-231c141 and MDA-231ba429 tumors; Supplemental File Table S1, Transcripts with miR-200c or miR-141 binding sites.

Author Contributions

K.E.S. helped with the human tumor cell injections, performed the vimentin and MXRA8 immunohistochemistry, performed RT-PCR on the tumor tissue, performed some of the miR-200 qRT-PCR, and edited the manuscript. K.L.W. maintained the animal colony, measured the tumors, collected tissue, extracted RNA, performed some of the miR-200 qRT-PCR, and edited the manuscript. R.A.M. performed the RNA-seq data analysis, wrote the manuscript, and directed project. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by a CHIR project grant (PJT-162218) awarded to R.A.M.

Institutional Review Board Statement

Animals were housed and cared for following guidelines established by the Central Animal Facility at the University of Guelph and the guidelines established by the Canadian Council of Animal Care. This study was approved by the Animal Care Committee at the University of Guelph (AUP #3995).

Informed Consent Statement

Not Applicable.

Data Availability Statement

RNA sequencing has been uploaded to GEO under accession number GSE193479.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ambros, V. The functions of animal microRNAs. Nature 2004, 431, 350–355. [Google Scholar] [CrossRef] [PubMed]
  2. 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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Bartel, D.P. MicroRNAs: Genomics, biogenesis, mechanism, and function. Cell 2004, 116, 281–297. [Google Scholar] [CrossRef] [Green Version]
  4. Di, L.G.; Calin, G.A.; Croce, C.M. MicroRNAs: Fundamental facts and involvement in human diseases. Birth Defects Res. C Embryo Today 2006, 78, 180–189. [Google Scholar] [CrossRef]
  5. Czech, B.; Hannon, G.J. Small RNA sorting: Matchmaking for Argonautes. Nat. Rev. Genet. 2011, 12, 19–31. [Google Scholar] [CrossRef] [PubMed]
  6. Hock, J.; Meister, G. The Argonaute protein family. Genome Biol. 2008, 9, 210. [Google Scholar] [CrossRef] [PubMed]
  7. Filipowicz, W.; Bhattacharyya, S.N.; Sonenberg, N. Mechanisms of post-transcriptional regulation by microRNAs: Are the answers in sight? Nat. Rev. Genet. 2008, 9, 102–114. [Google Scholar] [CrossRef]
  8. Bartel, D.P. MicroRNAs: Target recognition and regulatory functions. Cell 2009, 136, 215–233. [Google Scholar] [CrossRef] [Green Version]
  9. Lewis, B.P.; Shih, I.H.; Jones-Rhoades, M.W.; Bartel, D.P.; Burge, C.B. Prediction of mammalian microRNA targets. Cell 2003, 115, 787–798. [Google Scholar] [CrossRef] [Green Version]
  10. Brennecke, J.; Stark, A.; Russell, R.B.; Cohen, S.M. Principles of microRNA-target recognition. PLoS Biol. 2005, 3, e85. [Google Scholar] [CrossRef]
  11. Jonas, S.; Izaurralde, E. Towards a molecular understanding of microRNA-mediated gene silencing. Nat. Rev. Genet. 2015, 16, 421–433. [Google Scholar] [CrossRef] [PubMed]
  12. Saini, H.K.; Enright, A.J.; Griffiths-Jones, S. Annotation of mammalian primary microRNAs. BMC Genom. 2008, 9, 564. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Park, S.M.; Gaur, A.B.; Lengyel, E.; Peter, M.E. The miR-200 family determines the epithelial phenotype of cancer cells by targeting the E-cadherin repressors ZEB1 and ZEB2. Genes Dev. 2008, 22, 894–907. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Hill, L.; Browne, G.; Tulchinsky, E. ZEB/miR-200 feedback loop: At the crossroads of signal transduction in cancer. Int. J. Cancer 2012, 132, 745–754. [Google Scholar] [CrossRef] [PubMed]
  15. 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, 161–177. [Google Scholar] [CrossRef]
  16. Howe, E.N.; Cochrane, D.R.; Richer, J.K. The miR-200 and miR-221/222 microRNA families: Opposing effects on epithelial identity. J. Mammary Gland Biol. Neoplasia 2012, 17, 65–77. [Google Scholar] [CrossRef] [Green Version]
  17. Radisky, D.C. miR-200c at the nexus of epithelial-mesenchymal transition, resistance to apoptosis, and the breast cancer stem cell phenotype. Breast Cancer Res. 2011, 13, 110. [Google Scholar] [CrossRef]
  18. Bracken, C.P.; Gregory, P.A.; Kolesnikoff, N.; Bert, A.G.; Wang, J.; Shannon, M.F.; Goodall, G.J. A double-negative feedback loop between ZEB1-SIP1 and the microRNA-200 family regulates epithelial-mesenchymal transition. Cancer Res. 2008, 68, 7846–7854. [Google Scholar] [CrossRef] [Green Version]
  19. 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, 415–428. [Google Scholar] [CrossRef]
  20. Yang, J.; Weinberg, R.A. Epithelial-mesenchymal transition: At the crossroads of development and tumor metastasis. Dev. Cell 2008, 14, 818–829. [Google Scholar] [CrossRef] [Green Version]
  21. Bockmeyer, C.L.; 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, 735–745. [Google Scholar] [CrossRef] [PubMed]
  22. Cochrane, D.R.; Cittelly, D.M.; Howe, E.N.; Spoelstra, N.S.; McKinsey, E.L.; LaPara, K.; Elias, A.; Yee, D.; Richer, J.K. MicroRNAs link estrogen receptor α status and Dicer levels in breast cancer. Horm. Cancer 2010, 1, 306–319. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Korpal, M.; Lee, E.S.; 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, 14910–14914. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Mekala, J.R.; Naushad, S.M.; 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–258. [Google Scholar] [CrossRef] [PubMed]
  25. Humphries, B.; Wang, Z.; Oom, A.L.; 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, 2254–2263. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Prat, A.; Parker, J.S.; Karginova, O.; Fan, C.; Livasy, C.; Herschkowitz, J.I.; He, X.; Perou, C.M. Phenotypic and molecular characterization of the claudin-low intrinsic subtype of breast cancer. Breast Cancer Res. 2010, 12, R68. [Google Scholar] [CrossRef] [Green Version]
  27. Hennessy, B.T.; Gonzalez-Angulo, A.M.; Stemke-Hale, K.; Gilcrease, M.Z.; Krishnamurthy, S.; Lee, J.S.; Fridlyand, J.; Sahin, A.; Agarwal, R.; Joy, C.; et al. Characterization of a naturally occurring breast cancer subset enriched in epithelial-to-mesenchymal transition and stem cell characteristics. Cancer Res. 2009, 69, 4116–4124. [Google Scholar] [CrossRef] [Green Version]
  28. Herschkowitz, J.I.; Simin, K.; Weigman, V.J.; Mikaelian, I.; Usary, J.; Hu, Z.; Rasmussen, K.E.; Jones, L.P.; Assefnia, S.; Chandrasekharan, S.; et al. Identification of conserved gene expression features between murine mammary carcinoma models and human breast tumors. Genome Biol. 2007, 8, R76. [Google Scholar] [CrossRef] [Green Version]
  29. Simpson, K.; Conquer-van Heumen, G.; Watson, K.L.; Roth, M.; Martin, C.J.; Moorehead, R.A. 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. 2021, 21, 89. [Google Scholar] [CrossRef]
  30. Linnerth, N.M.; Sirbovan, K.; Moorehead, R.A. Use of a transgenic mouse model to identify markers of human lung tumors. Int. J. Cancer 2005, 114, 977–982. [Google Scholar] [CrossRef]
  31. Schneider, C.A.; Rasband, W.S.; Eliceiri, K.W. NIH Image to ImageJ: 25 years of image analysis. Nat. Methods 2012, 9, 671–675. [Google Scholar] [CrossRef] [PubMed]
  32. Chen, E.Y.; Tan, C.M.; Kou, Y.; Duan, Q.; Wang, Z.; Meirelles, G.V.; Clark, N.R.; Ma’ayan, A. Enrichr: Interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinform. 2013, 14, 128. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Kuleshov, M.V.; Jones, M.R.; Rouillard, A.D.; Fernandez, N.F.; Duan, Q.; Wang, Z.; Koplev, S.; Jenkins, S.L.; Jagodnik, K.M.; Lachmann, A.; et al. Enrichr: A comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res. 2016, 44, W90–W97. [Google Scholar] [CrossRef] [Green Version]
  34. Chorner, P.M.; Moorehead, R.A. A-674563, a putative AKT1 inhibitor that also suppresses CDK2 activity, inhibits human NSCLC cell growth more effectively than the pan-AKT inhibitor, MK-2206. PLoS ONE 2018, 13, e0193344. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Gyorffy, B.; Lanczky, A.; Eklund, A.C.; Denkert, C.; Budczies, J.; Li, Q.; Szallasi, Z. An online survival analysis tool to rapidly assess the effect of 22,277 genes on breast cancer prognosis using microarray data of 1809 patients. Breast Cancer Res. Treat. 2010, 123, 725–731. [Google Scholar] [CrossRef] [Green Version]
  36. Chen, Y.; Wang, X. miRDB: An online database for prediction of functional microRNA targets. Nucleic Acids Res. 2020, 48, D127–D131. [Google Scholar] [CrossRef] [Green Version]
  37. Liu, W.; Wang, X. Prediction of functional microRNA targets by integrative modeling of microRNA binding and target expression data. Genome Biol. 2019, 20, 18. [Google Scholar] [CrossRef]
  38. Paraskevopoulou, M.D.; Georgakilas, G.; Kostoulas, N.; Vlachos, I.S.; Vergoulis, T.; Reczko, M.; Filippidis, C.; Dalamagas, T.; Hatzigeorgiou, A.G. DIANA-microT web server v5.0: Service integration into miRNA functional analysis workflows. Nucleic Acids Res. 2013, 41, W169–W173. [Google Scholar] [CrossRef] [Green Version]
  39. Sticht, C.; De La Torre, C.; Parveen, A.; Gretz, N. miRWalk: An online resource for prediction of microRNA binding sites. PLoS ONE 2018, 13, e0206239. [Google Scholar] [CrossRef]
  40. Ding, J.; Li, X.; Hu, H. TarPmiR: A new approach for microRNA target site prediction. Bioinformatics 2016, 32, 2768–2775. [Google Scholar] [CrossRef] [Green Version]
  41. Zhang, D.D.; Li, Y.; Xu, Y.; Kim, J.; Huang, S. Phosphodiesterase 7B/microRNA-200c relationship regulates triple-negative breast cancer cell growth. Oncogene 2019, 38, 1106–1120. [Google Scholar] [CrossRef] [PubMed]
  42. D’Ippolito, E.; Plantamura, I.; Bongiovanni, L.; Casalini, P.; Baroni, S.; Piovan, C.; Orlandi, R.; Gualeni, A.V.; Gloghini, A.; Rossini, A.; et al. miR-9 and miR-200 Regulate PDGFRbeta-Mediated Endothelial Differentiation of Tumor Cells in Triple-Negative Breast Cancer. Cancer Res. 2016, 76, 5562–5572. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Shimono, Y.; Zabala, M.; Cho, R.W.; Lobo, N.; Dalerba, P.; Qian, D.; Diehn, M.; Liu, H.; Panula, S.P.; Chiao, E.; et al. Downregulation of miRNA-200c links breast cancer stem cells with normal stem cells. Cell 2009, 138, 592–603. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Song, C.; Liu, L.Z.; Pei, X.Q.; Liu, X.; Yang, L.; Ye, F.; Xie, X.; Chen, J.; Tang, H.; Xie, X. miR-200c inhibits breast cancer proliferation by targeting KRAS. Oncotarget 2015, 6, 34968–34978. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Williams, M.M.; Christenson, J.L.; O’Neill, K.I.; Hafeez, S.A.; Ihle, C.L.; Spoelstra, N.S.; Slansky, J.E.; Richer, J.K. MicroRNA-200c restoration reveals a cytokine profile to enhance M1 macrophage polarization in breast cancer. NPJ Breast Cancer 2021, 7, 64. [Google Scholar] [CrossRef] [PubMed]
  46. Liu, B.; Du, R.; Zhou, L.; Xu, J.; Chen, S.; Chen, J.; Yang, X.; Liu, D.X.; Shao, Z.M.; Zhang, L.; et al. miR-200c/141 Regulates Breast Cancer Stem Cell Heterogeneity via Targeting HIPK1/β-Catenin Axis. Theranostics 2018, 8, 5801–5813. [Google Scholar] [CrossRef]
  47. Knezevic, J.; Pfefferle, A.D.; Petrovic, I.; Greene, S.B.; Perou, C.M.; Rosen, J.M. Expression of miR-200c in claudin-low breast cancer alters stem cell functionality, enhances chemosensitivity and reduces metastatic potential. Oncogene 2015, 34, 5997–6006. [Google Scholar] [CrossRef] [Green Version]
  48. Korpal, M.; Ell, B.J.; Buffa, F.M.; Ibrahim, T.; Blanco, M.A.; Celia-Terrassa, T.; Mercatali, L.; Khan, Z.; Goodarzi, H.; Hua, Y.; et al. Direct targeting of Sec23a by miR-200s influences cancer cell secretome and promotes metastatic colonization. Nat. Med. 2011, 17, 1101–1108. [Google Scholar] [CrossRef] [Green Version]
  49. Li, X.; Roslan, S.; Johnstone, C.N.; Wright, J.A.; Bracken, C.P.; Anderson, M.; Bert, A.G.; Selth, L.A.; Anderson, R.L.; Goodall, G.J.; et al. MiR-200 can repress breast cancer metastasis through ZEB1-independent but moesin-dependent pathways. Oncogene 2014, 33, 4077–4088. [Google Scholar] [CrossRef]
  50. Jin, T.; Suk Kim, H.; Ki Choi, S.; Hye Hwang, E.; Woo, J.; Suk Ryu, H.; Kim, K.; Moon, A.; Kyung Moon, W. microRNA-200c/141 upregulates SerpinB2 to promote breast cancer cell metastasis and reduce patient survival. Oncotarget 2017, 8, 32769–32782. [Google Scholar] [CrossRef] [Green Version]
  51. 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. [Google Scholar] [CrossRef] [PubMed]
  52. Zhang, R.; Kim, A.S.; Fox, J.M.; Nair, S.; Basore, K.; Klimstra, W.B.; Rimkunas, R.; Fong, R.H.; Lin, H.; Poddar, S.; et al. Mxra8 is a receptor for multiple arthritogenic alphaviruses. Nature 2018, 557, 570–574. [Google Scholar] [CrossRef] [PubMed]
  53. Zhang, D.; Qian, C.; Wei, H.; Qian, X. Identification of the Prognostic Value of Tumor Microenvironment-Related Genes in Esophageal Squamous Cell Carcinoma. Front. Mol. Biosci. 2020, 7, 599475. [Google Scholar] [CrossRef] [PubMed]
  54. Li, S.; Xu, W. Mining TCGA database for screening and identification of hub genes in kidney renal clear cell carcinoma microenvironment. J. Cell Biochem. 2019, 121, 3952–3960. [Google Scholar] [CrossRef]
  55. Bhosale, P.G.; Cristea, S.; Ambatipudi, S.; Desai, R.S.; Kumar, R.; Patil, A.; Kane, S.; Borges, A.M.; Schaffer, A.A.; Beerenwinkel, N.; et al. Chromosomal Alterations and Gene Expression Changes Associated with the Progression of Leukoplakia to Advanced Gingivobuccal Cancer. Transl. Oncol. 2017, 10, 396–409. [Google Scholar] [CrossRef] [Green Version]
  56. Campbell, C.I.; Thompson, D.E.; Siwicky, M.D.; Moorehead, R.A. Murine mammary tumor cells with a claudin-low genotype. Cancer Cell Int. 2011, 11, 28. [Google Scholar] [CrossRef] [Green Version]
  57. Saleh, S.; Thompson, D.E.; McConkey, J.; Murray, P.; Moorehead, R.A. Osteopontin regulates proliferation, apoptosis, and migration of murine claudin-low mammary tumor cells. BMC Cancer 2016, 16, 359. [Google Scholar] [CrossRef] [Green Version]
  58. Watson, K.L.; Jones, R.A.; Bruce, A.; Moorehead, R.A. 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, 17–26. [Google Scholar] [CrossRef]
  59. Martin, C.J.; Moorehead, R.A. Polycomb repressor complex 2 function in breast cancer (Review). Int. J. Oncol. 2020, 57, 1085–1094. [Google Scholar] [CrossRef]
  60. Ferrari, K.J.; Scelfo, A.; Jammula, S.; Cuomo, A.; Barozzi, I.; Stutzer, A.; Fischle, W.; Bonaldi, T.; Pasini, D. Polycomb-dependent H3K27me1 and H3K27me2 regulate active transcription and enhancer fidelity. Mol. Cell 2014, 53, 49–62. [Google Scholar] [CrossRef] [Green Version]
Genes 13 00691 g001 550
Figure 1. Mammary tumor growth curves for (A) MDA-MB-231 cells transfected with an empty vector (MDA-231EV) or the miR-200c/141 cluster (MDA-231c141). Hematoxylin and eosin stained sections of tumors from (B) MDA-231EV or (C) MDA-231c141 injections. Scale bars are 60 µm. Expression of (D) miR-200c, (E) miR-141, (F) miR-200b, (G) miR-200a, and (H) miR-429 in the tumors induced by MDA-231EV or MDA-231c141 injections. miR-200 values were normalized to RNU44 and RNU48. The individual values have been indicated along with the mean ± SEM. A Student’s t-test was performed and significance is indicated as * p < 0.05, ** p < 0.01, ns, not significant.
Figure 1. Mammary tumor growth curves for (A) MDA-MB-231 cells transfected with an empty vector (MDA-231EV) or the miR-200c/141 cluster (MDA-231c141). Hematoxylin and eosin stained sections of tumors from (B) MDA-231EV or (C) MDA-231c141 injections. Scale bars are 60 µm. Expression of (D) miR-200c, (E) miR-141, (F) miR-200b, (G) miR-200a, and (H) miR-429 in the tumors induced by MDA-231EV or MDA-231c141 injections. miR-200 values were normalized to RNU44 and RNU48. The individual values have been indicated along with the mean ± SEM. A Student’s t-test was performed and significance is indicated as * p < 0.05, ** p < 0.01, ns, not significant.
Genes 13 00691 g001
Genes 13 00691 g002 550
Figure 2. Lung sections from mice harbouring (A–C) MDA-231EV or (D–F) MDA-231c141 primary mammary tumors. Panels (A,D) show hematoxylin and eosin stained lung metastases at lower magnification while, panels (C,F) show representative lung sections at higher magnification. Panels (B,E) show staining of the lung metastases with a human, anti-vimentin antibody and the percentage of the lungs containing cells staining positive for human vimentin are quantified in (G). Scale bars of (A,B,D,E) are 300 µm while scale bars for (C,F) are 60 µm. (G) Quantification of lung tumor burden. To calculate lung tumor burden, the area of the lung with positive vimentin staining was expressed relative to the total lung area. The individual values have been indicated along with the mean ± SEM. A Student’s t-test was performed and significance is indicated as ** p < 0.01.
Figure 2. Lung sections from mice harbouring (A–C) MDA-231EV or (D–F) MDA-231c141 primary mammary tumors. Panels (A,D) show hematoxylin and eosin stained lung metastases at lower magnification while, panels (C,F) show representative lung sections at higher magnification. Panels (B,E) show staining of the lung metastases with a human, anti-vimentin antibody and the percentage of the lungs containing cells staining positive for human vimentin are quantified in (G). Scale bars of (A,B,D,E) are 300 µm while scale bars for (C,F) are 60 µm. (G) Quantification of lung tumor burden. To calculate lung tumor burden, the area of the lung with positive vimentin staining was expressed relative to the total lung area. The individual values have been indicated along with the mean ± SEM. A Student’s t-test was performed and significance is indicated as ** p < 0.01.
Genes 13 00691 g002
Genes 13 00691 g003 550
Figure 3. The top (A) ENCODE and ChEA Consensus TFs from ChIP-X, (B) MSigDB Hallmark 2020, and (C) GO Biological Process 2021 identified in Enrichr using the 2843 differentially expressed genes in MDA-231EV tumors compared to MDA-231c141 tumors.
Figure 3. The top (A) ENCODE and ChEA Consensus TFs from ChIP-X, (B) MSigDB Hallmark 2020, and (C) GO Biological Process 2021 identified in Enrichr using the 2843 differentially expressed genes in MDA-231EV tumors compared to MDA-231c141 tumors.
Genes 13 00691 g003
Genes 13 00691 g004 550
Figure 4. Distant metastasis free survival (DMFS) curves from Kaplan-Meier plotter using basal-like breast cancers and median cut-offs for (A) CPE, (B) GNG2, (C) MXRA8, (D) SELENBP1, and (E) AGR2. The red line indicates elevated gene expression while the black line indicates reduced gene expression. The hazard ratio (HR) and significance (logrank P) are indicated in the upper right hand corner of each graph.
Figure 4. Distant metastasis free survival (DMFS) curves from Kaplan-Meier plotter using basal-like breast cancers and median cut-offs for (A) CPE, (B) GNG2, (C) MXRA8, (D) SELENBP1, and (E) AGR2. The red line indicates elevated gene expression while the black line indicates reduced gene expression. The hazard ratio (HR) and significance (logrank P) are indicated in the upper right hand corner of each graph.
Genes 13 00691 g004
Genes 13 00691 g005 550
Figure 5. The expression of (A) MXRA8 mRNA levels as determined by qRT-PCR and (B,C) MXRA8 protein levels as determined by Western blotting. MXRA8 mRNA levels were normalized to HPRT, while MXRA8 protein levels were normalized using GAPDH. (C) Quantitation of the MXRA8 band intensities from the Western blot relative to the GAPDH band intensities. In both A and C, the individual values have been indicated along with the mean ± SEM. A Student’s t-test was performed and significance is indicated as * p < 0.05; *** p < 0.001.
Figure 5. The expression of (A) MXRA8 mRNA levels as determined by qRT-PCR and (B,C) MXRA8 protein levels as determined by Western blotting. MXRA8 mRNA levels were normalized to HPRT, while MXRA8 protein levels were normalized using GAPDH. (C) Quantitation of the MXRA8 band intensities from the Western blot relative to the GAPDH band intensities. In both A and C, the individual values have been indicated along with the mean ± SEM. A Student’s t-test was performed and significance is indicated as * p < 0.05; *** p < 0.001.
Genes 13 00691 g005
Genes 13 00691 g006 550
Figure 6. Immunohistochemistry at (A) low and (B) high magnification for MXRA8 in lung metastases from MDA-231EV primary tumors. Scales bars are 600 µm for (A) and 100 µm for (B).
Figure 6. Immunohistochemistry at (A) low and (B) high magnification for MXRA8 in lung metastases from MDA-231EV primary tumors. Scales bars are 600 µm for (A) and 100 µm for (B).
Genes 13 00691 g006
Table
Table 1. Metastatic frequency of MDA-231EV and MDA-231c141 tumors.
Table 1. Metastatic frequency of MDA-231EV and MDA-231c141 tumors.
LungPancreas/SpleenLiver
MDA-231EV10/10 (100%)8/10 (80%)6/10 (60%)
MDA-231c1417/10 (70%)2/10 (20%)0/10 (0%)
Table
Table 2. Top five genes differentially expressed in MDA-231c141 vs. MDA-231EV tumors.
Table 2. Top five genes differentially expressed in MDA-231c141 vs. MDA-231EV tumors.
TranscriptlogFCFDR
CPE−8.81.42 × 10−223
GNG2−5.21.12 × 10−210
MXRA8−7.92.50 × 10−209
SELENBP1−7.44.88 × 10−160
AGR2−10.03.44 × 10−148
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Simpson, K.E.; Watson, K.L.; Moorehead, R.A. Elevated Expression of miR-200c/141 in MDA-MB-231 Cells Suppresses MXRA8 Levels and Impairs Breast Cancer Growth and Metastasis In Vivo. Genes 2022, 13, 691. https://0-doi-org.brum.beds.ac.uk/10.3390/genes13040691

AMA Style

Simpson KE, Watson KL, Moorehead RA. Elevated Expression of miR-200c/141 in MDA-MB-231 Cells Suppresses MXRA8 Levels and Impairs Breast Cancer Growth and Metastasis In Vivo. Genes. 2022; 13(4):691. https://0-doi-org.brum.beds.ac.uk/10.3390/genes13040691

Chicago/Turabian Style

Simpson, Kaitlyn E., Katrina L. Watson, and Roger A. Moorehead. 2022. "Elevated Expression of miR-200c/141 in MDA-MB-231 Cells Suppresses MXRA8 Levels and Impairs Breast Cancer Growth and Metastasis In Vivo" Genes 13, no. 4: 691. https://0-doi-org.brum.beds.ac.uk/10.3390/genes13040691

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop