Next Article in Journal
Structural Determinants of Substrate Specificity of SplF Protease from Staphylococcus aureus
Previous Article in Journal
Chemistry and Toxicology of Major Bioactive Substances in Inocybe Mushrooms
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 22
Issue 4
10.3390/ijms22042219
ijms-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

Runx3 Induces a Cell Shape Change and Suppresses Migration and Metastasis of Melanoma Cells by Altering a Transcriptional Profile

by
Ning Wang
1,
Haiying Zhang
2,
Xiulin Cui
1,
Chao Ma
1,
Linghui Wang
1 and
Wenguang Liu
1,*
1
Institute of Genetics and Cell Biology, School of Life Sciences, Northeast Normal University, No. 5268, Renmin St., Changchun 130024, China
2
Key Laboratory of Pathobiology of Ministry of Education, Norman Bethune College of Medicine, Jilin University, No. 126, Xinmin St., Changchun 130021, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2021, 22(4), 2219; https://doi.org/10.3390/ijms22042219
Submission received: 13 January 2021 / Revised: 8 February 2021 / Accepted: 18 February 2021 / Published: 23 February 2021
(This article belongs to the Section Molecular Oncology)
Browse Figures
Versions Notes

Abstract

:
Runt-related transcription factor-3 (Runx3) is a tumor suppressor, and its contribution to melanoma progression remains unclear. We previously demonstrated that Runx3 re-expression in B16-F10 melanoma cells changed their shape and attenuated their migration. In this study, we found that Runx3 re-expression in B16-F10 cells also suppressed their pulmonary metastasis. We performed microarray analysis and uncovered an altered transcriptional profile underlying the cell shape change and the suppression of migration and metastasis. This altered transcriptional profile was rich in Gene Ontology/Kyoto Encyclopedia of Genes and Genomes (GO/KEGG) annotations relevant to adhesion and the actin cytoskeleton and included differentially expressed genes for some major extracellular matrix (ECM) proteins as well as genes that were inversely associated with the increase in the metastatic potential of B16-F10 cells compared to B16-F0 melanoma cells. Further, we found that this altered transcriptional profile could have prognostic value, as evidenced by myelin and lymphocyte protein (MAL) and vilin-like (VILL). Finally, Mal gene expression was correlated with metastatic potential among the cells and was targeted by histone deacetylase (HDAC) inhibitors in B16-F10 cells, and the knockdown of Mal gene expression in B16-F0 cells changed their shape and enhanced the migratory and invasive traits of their metastasis. Our study suggests that self-entrapping of metastatic Runx3-negative melanoma cells via adhesion and the actin cytoskeleton could be a powerful therapeutic strategy.

1. Introduction

Cancer results in more than 8 million deaths every year worldwide, and cancer metastasis is the most frequent cause of these deaths, accounting for over 90% of mortalities [1,2]. Metastasis is hard to cure through current therapeutic regimens due to its systemic and complex nature and drug resistance [1,3,4,5,6]. As the molecular and cellular mechanisms involved in metastasis are still poorly understood, the discovery and development of effective treatments of cancer have been obstructed [3,5,6,7,8].
Cancer metastasis is a multistep cellular process; it involves cancer cell invasion, intravasation, survival in the circulation, arrest within the capillary bed, extravasation, and colonization at a distant site [1,3,7,8]. Along with this process, the local and foreign microenvironment, angiogenesis, platelet interaction, immunosurveillance, etc., together with many intrinsic characteristics of cancer cells, such as the epithelial–mesenchymal transition (EMT), energy metabolism, protease secretion, migratory capacity, exosome production, stemness, dormancy, and epi-/genetic drivers, affect successful metastasis [1,5,7,8]. Therefore, this process could provide unlimited candidates for various therapeutic targets [1,7,8,9].
Metastatic cancer cells can adapt their shape to be suitable for migration and invasion according to both environmental and genetic cues [5,10,11,12,13,14,15]. This cell shape change is a reflection of the dynamics in adhesion and the cytoskeleton, as well as other molecular and cellular processes [16]. Both of these, as well as the interplay between them, determine cancer cells’ migratory capacity [1,17,18,19], so these two evolutionarily conserved cellular processes could be therapeutically targeted for the treatment of metastasis [7,20]. The specific mechanisms that determine favorable cell shape changes for cancer cell migration and invasion during metastasis are not completely understood. Thus, the identification and characterization of key factors that can regulate cancer cells’ morphology would provide new insights into the pathology of metastasis [13].
The extracellular matrix (ECM) is a complex architecture that is capable of remodeling, and it supports cancer cell proliferation, survival, adhesion, migration, and invasion [21,22]. Many early and recent studies have thoroughly elucidated the cancer biology of individual ECM proteins or their genes, and a number of candidates for metastasis markers and therapeutic targets have been identified [23,24,25,26,27,28,29,30,31,32,33]. Although ECM proteins are very important in cancer metastasis, they remain underexplored [29,34].
The transcription factor runt-related transcription factor-3 (Runx3) plays a pivotal role in neurogenesis and thymopoiesis [35]. The depletion of Runx3 is a cause of gastric cancer development, and evidence that inactivation of Runx3 expression or function is a key event in gastric cancer and many other types of cancer has accumulated [36,37,38,39]. Loss of Runx3 expression has also frequently been found in melanoma cells as compared to melanocytes, and there is a good negative correlation between Runx3 expression and melanoma malignancy [40]. In addition, many functional studies have demonstrated that Runx3 is a tumor suppressor in many types of cancers but is a tumor promoter in other types [37,38,39,41,42,43,44].
The mode of action of Runx3 in the suppression of tumor progression is complex, but evidence shows that Runx3 governs the EMT during tumor progression [45,46,47,48]. The morphological changes in breast cancer cells are reminiscent of back–epithelial transitions and are closely accompanied by predominantly decreased Matrigel invasion after an increase in exogenous Runx3 expression [36]. In metastatic melanoma cells, a number of genes, including the tumor suppressor PTEN, were recently identified with a screen for morphological complexity as functionally conserved cell shape regulators that are relevant to disease progression [13]. Like PTEN, Runx3 may also function as a cell shape regulator when it suppresses melanoma cell migration and invasion in vitro [49]. Recently, additional roles of Runx3 in metastasis have been implicated in processes such as regulation of the cancer microenvironment, angiogenesis, stemness, immunosurveillance, and inflammation [42,48,50,51,52].
Myelin and lymphocyte protein (Mal) is a four-time transmembrane protein, and it plays a role in adhesion and the actin cytoskeleton [53,54]. Mal is also a tumor suppressor [53,55], and it may be involved in suppressing metastasis [56]. As a transcription factor, Runx3 has to fulfill its anti-metastasis function via its downstream genes. Whether Mal can be a therapeutic target for the treatment of metastasis is potentially interesting.
In this study, we used a metastatic cell model of mouse melanoma to further investigate the role of Runx3 in metastasis and to uncover the underlying transcriptional profile. We hope that this study can help in the discovery of therapeutic targets for the treatment of metastasis.

2. Materials and Methods

2.1. Generation of Stable Cell Lines by Lentiviral Infection

For the mock control cell line, transfection was performed with the empty lentiviral construct PLKO.1 and the packaging constructs psPAX2 and pMD.2G using the Lipofectamine 2000 transfection reagent (11668019, Invitrogen, Carlsbad, CA, USA). The viral supernatant was harvested at 48 h and filtered through a 0.45 µm filter. For the stable knockdown cell line for Mal, lentiviral particles were purchased from Santa Cruz (sc-44786-v, Mal shRNA (m)). B16-F0 mouse melanoma cells were infected with the lentiviral supernatants with the aid of 5 µg/mL of polybrene. All stable colonies were selected through treatment with 5 µg/mL of puromycin. The success of the knockdown of Mal gene expression in puromycin-resistant colonies was confirmed by RT-qPCR and Western blotting.

2.2. Cell Culture and Treatment, Cell Proliferation Assay, and Cell Cycle Assay

The generation of the Runx3-re-expressing B16-F10 melanoma cell line and the cultures of all the different melanoma cell lines in this study (B16-F10, mock control B16-F10, B16-F0, mock control B16-F0, B16-F10/Runx3, and B16-F0/shMal) have been previously described [49]. B16-F10 has high pulmonary metastatic ability, whereas B16-F0 is the parent cell line of B16-F10 and its pulmonary metastatic ability is low [57,58]. Histone deacetylase (HDAC) inhibitors (S1095, S2132, S2818, S2244, and S8464) were purchased from Selleck, and treatment of cells with an individual inhibitor was maintained for 18 h at 10 μM. The cell proliferation assay and the cell cycle assay were performed, as previously described, using the MTT method and flow cytometry for PI (Propidium Iodide) (Cat. No. BD5012, Bioworld)-stained cells, respectively [49,59].

2.3. Wound-Healing Assay and Transwell Migration/Invasion Assay

To address their migratory potential, cells were grown to 100% confluence in a culture medium containing a normal concentration of serum (10%). At hour 0, a 0.5-mm-wide wound was made with a rubber scrapper, and the detached cells were washed off with phosphate-buffered saline (PBS) buffer. The remaining cells were maintained in a medium supplemented with 1% serum for 48 h. During this period, cells were allowed to migrate into the wounded area. Transwell migration and invasion assays were performed, as previously described [49], using Costar Transwell chambers (Cat. No. 3422, Corning) with or without a Matrigel (Cat. No. 356234, Corning) coating. In brief, 10e5 cells were seeded on the upper part of the chamber with a medium containing 1% FBS, and the lower part of the chamber was incubated with a medium containing 10% FBS. Migrated or invaded cells on the membrane were stained with crystal violet.

2.4. Subcutaneous Tumor Formation and Pulmonary Metastasis Assay

C57BL/6J mice were bought from the Animal Experimental Center of Jilin University (Changchun, China). All mice were housed in a pathogen-free and air-conditioned environment with light cycles of 12 h on and 12 h off and free access to tap water and food. The starting age of experimental mice was 6–10 weeks old. Where applicable, mice were anesthetized. When assessing the tumor volume, a diameter of less than or around 1.0 cm per tumor was used. Tumor formation was generated by subcutaneous injection of 10e6 cells per mouse, and the tumors were allowed to grow for 3 weeks. The mice were sacrificed, and the tumors were carefully harvested. To detect pulmonary metastasis, 10e5 cells were intravenously injected per mouse, and the mice were sacrificed to record the metastasis foci in the lungs at 3 or 4 weeks after injection.

2.5. Microarray

SmartArray microarray chips (CapitalBio, Beijing, China) were prepared from 32K Mouse Genome Array version 4.0 (http://www.Operon.com, accessed on 25 April 2011). Total RNA was prepared using Trizol reagent (15596026, Invitrogen, ThermoFisher, Waltham, MA, USA) and was purified using the NucleoSpin RNA Clean-up kit (Macherey-Nagel, Düren, Germany). cRNA was amplified and labeled with Cy5- or Cy3-dCTP using the Jingxin cRNA amplification and labeling kit (CapitalBio, Beijing, China) after the synthesis of double-stranded cDNA. Hybridization on chips was performed in a hybridization buffer (3× SSC, 0.2% SDS, 5× Denhart’s, 25% formamide) at 42 °C overnight. After washing, the chips were scanned using a LuxScan 10 KA scanner, and the signals were analyzed using LuxScan 3.0 software (CapitalBio, Beijing, China). Differentially expressed genes from the comparative analysis were picked with a 0.5- or 2-fold threshold. Pathway and Gene Ontology (GO) annotations were performed using Molecule Annotation System version 3.0 (CapitalBio, Beijing, China).

2.6. RNA Preparation and RT-qPCR

Total RNA was extracted using Trizol reagent. Reverse transcription from the RNA was performed using the Takara PrimeScript RT reagent kit with gDNA Eraser (RR047A, Takara, Dalian, China). qPCR was performed using the 2× SYBR Green qPCR Master Mix (B21202, Bimake, Shanghai, China), and the PCR protocol was 95 °C for 15 s and 60 °C for 1 min for 40 cycles. The primer sequences used in this study were as follows: actin-forward, 5′-AACAGTCCGCCTAGAAGCAC-3′; actin-reverse, 5′-CGTTGACATCCGTAAAGACC-3′; Mal-forward, 5′-ATGCAGCCTACCACTGTGTG-3′; Mal-reverse, 5′-CAGCATGGACCACGTAGATC-3′; vilin-like (Vill)-forward, 5′-GAGGACCTAGGAGACCAGAC-3′; and Vill-reverse, 5′-ATTCTCCGCATCTACCACAG-3′.

2.7. Western Blotting

Western blotting was performed, as previously described [49]. In brief, cells were lysed with extraction buffer. Total proteins were separated after electrophoresis and were transferred onto an NC (nitrocellulose) membrane (pore size: 0.22 um; GE Healthcare). The intensities of the bands were recorded by a Tanon chemiluminescence detector (S5200, Tanon, Shanghai, China) after the signals were developed using the ECL (enhanced chemiluminescence) system (GE Healthcare). The protein levels of Mal and β-actin were indirectly detected using anti-Mal antibodies (ab15418, Abcam, Shanghai, China) and anti-β-actin antibodies (No. 3700, CST), respectively.

2.8. Histology

For histological analysis, subcutaneous tumors were fixed in PBS buffer containing 10% formaldehyde and were embedded in paraffin. Sections were stained with hematoxylin and eosin. The cell shape was observed under a microscope.

2.9. Actin Cytoskeleton Staining

Cells were seeded onto cover glasses and were allowed to grow to near 100% confluence. Cells were also seeded, allowed to grow to 100% confluence, and then subjected to a wound-healing assay. Cells were washed with PBS, fixed with 10% formaldehyde for 10 min, and permeabilized in 0.1% Triton X-100 for 3 min. Cells were then stained for 15 min with phalloidin-FITC (ab235137, Abcam). Photos were captured with an Olympus fluorescence microscope.

2.10. The Cancer Genome Atlas (TCGA) Data Mining and Gene Expression Correlation Analysis

For MAL and VILL, we obtained the fragments per kilobase of exon per million reads (FPKM) expression data and the related clinical data of patients with various types of cancer from The Cancer Genome Atlas (TCGA) database (https://www.proteinatlas.org, accessed on 1 September 2020). Patients were stratified into MAL- or VILL-high and MAL- or VILL-low groups. Survival curves were made with GraphPad Prism 7 software. Correlation analysis of the gene expression of MAL with that of other genes was performed using an online tool (http://www.linkedomics.org/login.php) [60].

2.11. Statistical Analysis

The t-test and log-rank test were performed using GraphPad Prism 7 software. The Spearman correlation test was performed using an online tool (http://www.linkedomics.org/login.php, accessed on 1 October 2020) [60]. Data were shown as the mean ± SEM. p-Values equal to or less than 0.05 were considered statistically significant. n denotes the number of experiments or samples.

3. Results

3.1. Runx3 Re-Expression Induces a Cell Shape Change and Suppresses the Migration and Metastasis of B16-F10 Melanoma Cells

Melanocytes are derived from neural crest cells, and they are generally regarded as branched/dendritic pigment cells [61,62]. Melanoma has an origin in melanocytes [63]. Therefore, the cell shape of B16-F10 melanoma cells was described as elongated and branched when they were cultured in vitro [49]. Upon restoration of Runx3 expression, B16-F10 cells underwent a significant cell shape change from elongated and branched to spread and unbranched [49] (Figure 1A). In the meantime, stress fiber formation was considerably enhanced [49] (Figure 1A). In association with the cell shape change, migration was significantly attenuated [49] (Figure 1B and Figure S1). Interestingly, stress fibers in the migrating B16-F10/Runx3 cells degenerated, with the exception of their ends (Figure 1B as compared to Figure 1A). In contrast, along the migration front, some B16-F10 mock control cells showed microspike-like structures (Figure 1B as compared to Figure 1A), which were never formed in the Runx3-expressing B16-F10 cells (Table S1). Thus, it seemed that B16-F10/Runx3 cells entrapped themselves to at least slow down their migration. Further, although the B16-F10 mock control cells demonstrated a high potential of pulmonary metastasis, as expected [57], the metastatic potential was abolished by Runx3, whereas subcutaneous tumor formation was reduced by Runx3 to 61.5% compared to the mock control cells (Figure 1C,D). Meanwhile, the morphological differences between B16-F10/Runx3 and B16-F10 cells remained in vivo (Figure 1E).

3.2. Runx3 Re-Expression Alters a Specific Transcriptional Profile in B16-F10 Melanoma Cells

We performed microarray analysis for the identification of differentially expressed genes (DEGs) regulated by Runx3 re-expression in B16-F10 (high-metastatic) melanoma cells [57]. We also performed microarray analysis for the identification of DEGs in the comparison of B16-F0 (low-metastatic) melanoma cells with B16-F10 cells, since B16-F10 cells gained a high potential for pulmonary metastasis from B16-F0 cells via unknown genetic and/or epigenetic mechanisms [58]. In general, 783 DEGs and 499 DEGs were regulated in the case of Runx3 re-expression in B16-F10 cells and B16-F0 vs. B16-F10 cells, respectively. The differential expression of a few DEGs was confirmed previously and in this study for the validation of microarray data [49,64] (Figure S2). With an emphasis on the aforementioned cell shape change and its associated attenuation of the migratory and metastatic potential, we performed molecular annotation analysis of the microarray data using the GO database and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. According to the GO_cellular component/biological process/molecular function annotations, we found in both cases that the various repertoires of DEGs were assigned to the terms of GO:0005886-plasma membrane, GO:0005856-cytoskeleton, GO:0007155-cell adhesion, GO:0003779/0051015-actin/actin filament binding, and GO:0030054-cell junction (Figure 2A). According to the KEGG_pathway annotation, we found in both cases that the different repertoires of DEGs were assigned to the terms cell adhesion molecules (CAMs), focal adhesion, and regulation of the actin cytoskeleton (Figure 2A).
Notably, we found that Runx3 specifically upregulated a number of DEGs that encode major ECM proteins in B16-F10 cells, some of which have been reported to play key roles in focal adhesion formation [23,24,65] (Figure 2B).
Interestingly, among the upregulated DEGs in the case of B16-F0 vs. B16-F10 cells, 61 were also found to be upregulated by Runx3 in B16-F10 cells; among the downregulated DEGs in the case of B16-F0 vs. B16-F10 cells, 27 were also downregulated by Runx3 in B16-F10 cells, indicating that there was a commonly altered transcriptional profile in both cases (Figure 2C). In addition, 10 CAM-encoding DEGs (H2-K1, H2-Bl, H2-Q6, H2-Q7, H2-Q1, H2-D1, H2-T3, H2-T9, H2-T23, and Mcam) were upregulated in both cases (Figure 2C); indeed, downregulation of murine CAMs has frequently been found in metastatic tumors [66].

3.3. The Runx3-Altered Transcriptional Profile May Have Prognostic Value for Various Cancers, as Evidenced by MAL and/or VILL

To address the clinical relevance of the Runx3-altered transcriptional profile, we selected Mal and Vill from the list of commonly altered DEGs to perform data mining using the gene expression data and the corresponding clinical parameters that were available from The Cancer Genome Atlas (TCGA). In Kaplan–Meier survival analysis, patients were stratified into MAL-/VILL-high and MAL-/VILL-low groups for each type of cancer. MAL has previously been reported to be a metastasis suppressor in esophageal cancer, gastric cancer, and head and neck squamous cell carcinoma and seems to confer survival advantages for relevant patients [53,56,67,68]. We demonstrated here that high-MAL gene expression seems to confer survival advantages for patients with melanoma, lung adenocarcinoma, pancreatic adenocarcinoma, renal cancer, or colorectal adenocarcinoma (Figure 3A). VILL’s role in metastasis is unclear. We demonstrated that low-VILL gene expression seems to confer survival advantages for patients with melanoma, lung adenocarcinoma, pancreatic adenocarcinoma, renal cancer, or glioma (Figure 3A). Notably, we demonstrated that both high-MAL gene expression and low-VILL gene expression seem to synergistically confer survival advantages for patients with melanoma, lung adenocarcinoma, pancreatic adenocarcinoma, or renal cancer (Figure 3B). We also performed Kaplan–Meier survival analysis for RUNX3 within the same types of cancer. We demonstrated that RUNX3 does not have prognostic value or a positive prognostic value, even though it is usually regarded as a tumor suppressor gene (Figure S3).

3.4. Mal Gene Expression Is Correlated with Metastatic Potential among Cellular Models and Is Subject to Epigenetic Regulation

Based on our microarray data and the results of our survival analysis, we selected Mal for further investigation. We validated the gene expression of Mal in B16-F10, B16-F0, and B16-F10/Runx3 cells. As shown in Figure 4A, both B16-F10/Runx3 cells and B16-F0 cells had an increased level of Mal mRNA compared to B16-F10 cells. As Runx3 suppressed metastasis of B16-F10 cells and B16-F0 cells are known to have a lower metastatic potential than B16-F10 cells, we believe that the level of Mal gene expression is correlated with the metastatic potential among the melanoma cells used in this study.
To address the pharmaceutical regulation of Mal gene expression, we randomly selected five HDAC inhibitors; two of these molecules are pan-HDAC inhibitors, and the others are HDAC2/8-selective inhibitors. We found that these HDAC inhibitors dramatically induced Mal gene expression in B16-F10 cells (Figure 4B), suggesting that Mal gene expression is subject to epigenetic regulation.
We also performed correlation analysis between the gene expression of MAL and that of other genes in melanoma (Figure 4C). We found that MAL was positively correlated with a known metastasis suppressor, BRMS1 [7] (Figure 4C). We also found that MAL was negatively correlated with a few genes involved in de novo methylation of genomic DNA or deacetylation of histones (Figure 4C).

3.5. Mal Is Necessary for Suppressing the Migratory and Invasive Traits of Metastasis in Melanoma Cells

To investigate the functional roles of Mal in melanoma, we lentivirally reduced the gene expression of Mal in B16-F0 cells (Figure 5A). We found that the knockdown of Mal gene expression apparently resulted in a cell shape change from spread to elongated cells (Figure 5B, Figure S4). In line with the cell shape change, the knockdown of Mal gene expression altered the actin cytoskeleton and, notably, induced the formation of microspike-like structures in some cells along the migration front (Figure 5B). Meanwhile, the knockdown of Mal gene expression slightly reduced the proliferation rate, as shown by the altered progress of the cell cycle, but had no effect on the subcutaneous tumor formation rate (Figure 5C–E). However, it significantly enhanced the migratory and invasive potential (Figure 5F–H). Unfortunately, we found that the knockdown of Mal gene expression did not dramatically increase the number of metastasis foci in the lungs (Figure 5I). Instead, it increased the frequency of metastasis-focus-positive lungs (Figure 5I).

4. Discussion

The search for therapeutic targets for metastasis has been a continuous journey [1,7,8,9,14,20,22,69,70]. Targeting cancer cell motility can be an effective strategy for metastasis therapy because cell motility is a key driving trait of metastasis and is required for multiple stages of metastasis at least via adhesion and the actin cytoskeleton [1,20]. In addition, increases in the metastatic potential of melanoma cells could be associated with cell shape change [13], a phenomenon that has already been noticed [26,57]. In this study, by using cellular models of mouse melanoma, we uncovered a specific Runx3-regulated transcriptional profile that contained DEGs that encoded some ECM proteins, as well as DEGs that were inversely associated with increases in the metastatic potential. This profile was correlated with the phenotype of the cell shape change and the attenuation of migration and metastasis, with comprehensive annotations referring to the actin cytoskeleton and adhesion. The DEGs were identified using microarray data of mouse melanoma, but they can be interpreted in future clinical studies, as evidenced by others [71]. Thus, we suggest that the self-entrapping of cells via adhesion and the actin cytoskeleton could be a powerful strategy for metastasis therapy for melanoma.
Whether Runx3 is a metastasis suppressor in melanoma is an important issue. First, it was reported that Runx3 expression is frequently lost in melanoma and that there is a good negative correlation between Runx3 expression and melanoma malignancy [40]. Second, we found that Runx3 re-expression abolishes the pulmonary metastasis of B16-F10 cells, but it slightly affects the subcutaneous tumor formation rate (Figure 1C,D). This evidence is in accordance with the functional definition of a metastasis suppressor by Khan and Steeg [7]. Third, and most importantly, we found that Runx3 can regulate a specific profile of DEGs that is inversely associated with an increase in the metastatic potential of B16-F10 cells compared to B16-F0 cells (Figure 2C).
It is a distinct feature that Runx3 re-expression in B16-F10 cells upregulates fibronectin (Fn) and the other major ECM genes without downregulating any other ECM genes (Figure 2B). This feature may implicate a metastasis-suppressing role played by Fn, even though Fn is usually regarded as an EMT marker in many cancer studies. In fact, early studies suggested that the loss of Fn is correlated with the acquisition of metastasis [72]. Recently, the role of Fn in protecting against the EMT process was confirmed in vivo [73]. In addition, the metastatic potential seems to be strongly related to the binding capacity of the specific constituent combination of Fn with other ECM molecules [34]. Therefore, the roles of ECM proteins in metastasis should be considered carefully.
Several previous reports have demonstrated that Mal is a tumor suppressor [53,55]. Mal is a four-time transmembrane protein; in particular, it is a lipid raft constituent [53], and it could be a therapeutic target because of its features as a cell surface protein [14]. Mal is thus a participant in RhoA signaling and has a role in stress fiber formation and focal adhesion formation [54]. Mal expression in esophageal cancer cells suppresses migration and invasion [67]. Suppression of Mal expression in gastric cancer is correlated with metastasis [68]. Downregulation of Mal expression may contribute to the metastatic potential of head and neck squamous cell carcinoma [53,56]. In this study, it was demonstrated that Mal is necessary for suppressing the migratory and invasive traits of metastasis in mouse melanoma cells; its expression can be of prognostic value and can be subject to epigenetic regulation. Since transcription factors are hard to target in cancers, we suggest that Mal can be a therapeutic target for the treatment of metastasis of melanoma cells.
Like MAL, VILL gene expression also has a prognostic value, and an increase in MAL gene expression and/or a decrease in VILL gene expression may confer survival advantages in several types of cancer. In view of making use of MAL and VILL for the treatment of metastasis of these types of cancer, one can use small molecules to enhance MAL gene expression or function and/or to attenuate VILL gene expression or function. Our result shows that Mal gene expression would be pharmaceutically enhanced by HDAC inhibitors in metastatic mouse melanoma cells. Moreover, this result implicates that there could be an HDAC-mediated epigenetic mechanism when B16-F10 cells gain high potential for pulmonary metastasis compared to B16-F0 cells. Further studies are needed to address whether there are HDAC-mediated epigenetic mechanisms in human melanoma. If there exists such a mechanism that could be implicated in melanoma metastasis, a specific HDAC inhibitor could be useful in the future for the treatment of metastasis of melanoma cells.

Supplementary Materials

The following are available online at https://www.mdpi.com/1422-0067/22/4/2219/s1: Figure S1: Runx3 expression in B16-F10 melanoma cells delayed the migration rate in wound healing (quantification); Figure S2: The Vill gene expression was validated in the cell lines used in this study with qRT-PCR; Figure S3: The gene expression of RUNX3 had no prognostic/positive prognostic value for various cancers; Figure S4: The knockdown of Mal gene expression changed the cell shape (at subconfluence); and Table S1: Frequency of microspike-like structures per 10 random views.

Author Contributions

N.W., H.Z., X.C.: investigation, visualization, and writing—review and editing; C.M.: data curation; L.W.: investigation; W.L.: conceptualization, project administration, funding acquisition, data curation, writing—original draft preparation, and writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Natural Science Foundation of Jilin Province (20180101257JC).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki. All mouse experiments were approved by the Animal Care Committee of Northeast Normal University (no. 2013-05; date of approval: 2013-5-1).

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw microarray data and the processed microarray data can be viewed at https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE161991 and https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE161992.

Acknowledgments

We would like to thank Qi Chen for her early assistance in the comparison of microarray data.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

DEGdifferentially expressed gene
ECMextracellular matrix
EMTepithelial–mesenchymal transition
FPKMfragments per kilobase of exon per million reads
GOGene Ontology
HDAChistone deacetylase
KEGGKyoto Encyclopedia of Genes and Genomes
Malmyelin and lymphocyte protein
PBSphosphate-buffered saline
qRT-PCRquantitative reverse transcription–polymerase chain reaction
Runx3runt-related transcription factor-3
TCGAThe Cancer Genome Atlas
Villvilin-like

References

  1. Stoletov, K.; Beatty, P.H.; Lewis, J.D. Novel therapeutic targets for cancer metastasis. Expert Rev. Anticancer 2020, 20, 97–109. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. The ICGC/TCGA Pan-Cancer Analysis of Whole Genomes Consortium. Pan-cancer analysis of whole genomes. Nature 2020, 578, 82–93. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Valastyan, S.; Weinberg, R.A. Tumor metastasis: Molecular insights and evolving paradigms. Cell 2011, 147, 275–292. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Fidler, I.J.; Kripke, M.L. The challenge of targeting metastasis. Cancer Metastasis Rev. 2015, 34, 635–641. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Lambert, A.W.; Pattabiraman, D.R.; Weinberg, R.A. Emerging Biological Principles of Metastasis. Cell 2017, 168, 670–691. [Google Scholar] [CrossRef] [Green Version]
  6. Suhail, Y.; Cain, M.P.; Vanaja, K.; Kurywchak, P.A.; Levchenko, A.; Kalluri, R. Systems Biology of Cancer Metastasis. Cell Syst. 2019, 9, 109–127. [Google Scholar] [CrossRef] [Green Version]
  7. Khan, I.; Steeg, P.S. Metastasis suppressors: Functional pathways. Lab. Investig. 2018, 98, 198–210. [Google Scholar] [CrossRef] [Green Version]
  8. Anderson, R.L.; Balasas, T.; Callaghan, J.; Coombes, R.C.; Evans, J.; Hall, J.A.; Kinrade, S.; Jones, D.; Jones, P.S.; Jones, R.; et al. A framework for the development of effective anti-metastatic agents. Nat. Rev. Clin. Oncol. 2019, 16, 185–204. [Google Scholar] [CrossRef] [Green Version]
  9. Steeg, P.S. Targeting metastasis. Nat. Rev. Cancer. 2016, 16, 201–218. [Google Scholar] [CrossRef]
  10. Sanz-Moreno, V.; Gadea, G.; Ahn, J.; Paterson, H.; Marra, P.; Pinner, S.; Sahai, E.; Marshall, C.J. Rac activation and inactivation control plasticity of tumor cell movement. Cell 2008, 135, 510–523. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Keren, K.; Pincus, Z.; Allen, G.M.; Barnhart, E.L.; Marriott, G.; Mogilner, A.; Theriot, J.A. Mechanism of shape determination in motile cells. Nature 2008, 453, 475–480. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Barnhart, E.L.; Lee, K.C.; Keren, K.; Mogilner, A.; Theriot, J.A. An adhesion-dependent switch between mechanisms that determine motile cell shape. PLoS Biol. 2011, 9, e1001059. [Google Scholar] [CrossRef]
  13. Yin, Z.; Sadok, A.; Sailem, H.; McCarthy, A.; Xia, X.; Li, F.; Garcia, M.A.; Evans, L.; Barr, A.R.; Perrimon, N.; et al. A screen for morphological complexity identifies regulators of switch-like transitions between discrete cell shapes. Nat. Cell Biol. 2013, 15, 860–871. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Chakraborty, S.; Lakshmanan, M.; Swa, H.L.F.; Chen, J.; Zhang, X.; Ong, Y.S.; Loo, L.S.; Akincilar, S.C.; Gunaratne, J.; Tergaonkar, V.; et al. An oncogenic role of Agrin in regulating focal adhesion integrity in hepatocellular carcinoma. Nat. Commun. 2015, 6, 6184. [Google Scholar] [CrossRef]
  15. Roszkowska, K.A.; Gizinski, S.; Sady, M.; Gajewski, Z.; Olszewski, M.B. Gain-of-Function Mutations in p53 in Cancer Invasiveness and Metastasis. Int. J. Mol. Sci. 2020, 21, 1334. [Google Scholar]
  16. Mogilner, A.; Keren, K. The shape of motile cells. Curr. Biol. 2009, 19, R762–R771. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Gardel, M.L.; Schneider, I.C.; Aratyn-Schaus, Y.; Waterman, C.M. Mechanical integration of actin and adhesion dynamics in cell migration. Annu. Rev. Cell Dev. Biol. 2010, 26, 315–333. [Google Scholar]
  18. Chang, C.W.; Kumar, S. Vinculin tension distributions of individual stress fibers within cell-matrix adhesions. J. Cell Sci. 2013, 126, 3021–3030. [Google Scholar] [CrossRef] [Green Version]
  19. Paul, C.D.; Mistriotis, P.; Konstantopoulos, K. Cancer cell motility: Lessons from migration in confined spaces. Nat. Rev. Cancer 2017, 17, 131–140. [Google Scholar] [PubMed] [Green Version]
  20. Wangpu, X.; Lu, J.; Xi, R.; Yue, F.; Sahni, S.; Park, K.C.; Menezes, S.; Huang, M.L.H.; Zheng, M.; Kovacevic, Z.; et al. Targeting the Metastasis Suppressor, N-Myc Downstream Regulated Gene-1, with Novel Di-2-Pyridylketone Thiosemicarbazones: Suppression of Tumor Cell Migration and Cell-Collagen Adhesion by Inhibiting Focal Adhesion Kinase/Paxillin Signaling. Mol. Pharm. 2016, 89, 521–540. [Google Scholar] [CrossRef] [Green Version]
  21. Hynes, R.O. The extracellular matrix: Not just pretty fibrils. Science 2009, 326, 1216–1219. [Google Scholar] [CrossRef] [Green Version]
  22. Huang, H. Matrix Metalloproteinase-9 (MMP-9) as a Cancer Biomarker and MMP-9 Biosensors: Recent Advances. Sensors 2018, 18, 3249. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Hynes, R.O.; Destree, A.T. Relationships between fibronectin (LETS protein) and actin. Cell 1978, 15, 875–886. [Google Scholar] [CrossRef]
  24. Lawler, J.; Weinstein, R.; Hynes, R.O. Cell attachment to thrombospondin: The role of ARG-GLY-ASP, calcium, and integrin receptors. J. Cell Biol. 1988, 107, 2351–2361. [Google Scholar] [CrossRef]
  25. Roberts, D.D. Regulation of tumor growth and metastasis by thrombospondin-1. FASEB J. 1996, 10, 1183–1191. [Google Scholar] [CrossRef]
  26. Clark, E.A.; Golub, T.R.; Lander, E.S.; Hynes, R.O. Genomic analysis of metastasis reveals an essential role for RhoC. Nature 2000, 406, 532–535. [Google Scholar] [CrossRef] [PubMed]
  27. Armstrong, L.C.; Bornstein, P. Thrombospondins 1 and 2 function as inhibitors of angiogenesis. Matrix Biol. 2003, 22, 63–71. [Google Scholar] [CrossRef]
  28. Murphy-Ullrich, J.E.; Iozzo, R.V. Thrombospondins in physiology and disease: New tricks for old dogs. Matrix Biol. 2012, 31, 152–154. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Naba, A.; Clauser, K.R.; Lamar, J.M.; Carr, S.A.; Hynes, R.O. Extracellular matrix signatures of human mammary carcinoma identify novel metastasis promoters. eLife 2014, 3, e01308. [Google Scholar] [CrossRef] [Green Version]
  30. Venning, F.A.; Wullkopf, L.; Erler, J.T. Targeting ECM Disrupts Cancer Progression. Front Oncol. 2015, 5, 224. [Google Scholar] [CrossRef] [Green Version]
  31. Yamauchi, M.; Barker, T.H.; Gibbons, D.L.; Kurie, J.M. The fibrotic tumor stroma. J. Clin. Investig. 2018, 128, 16–25. [Google Scholar] [CrossRef]
  32. Paolillo, M.; Schinelli, S. Extracellular Matrix Alterations in Metastatic Processes. Int. J. Mol. Sci. 2019, 20, 4947. [Google Scholar] [CrossRef] [Green Version]
  33. Ma, H.P.; Chang, H.L.; Bamodu, O.A.; Yadav, V.K.; Huang, T.Y.; Wu, A.T.H.; Yeh, C.T.; Tsai, S.H.; Lee, W.H. Collagen 1A1 (COL1A1) Is a Reliable Biomarker and Putative Therapeutic Target for Hepatocellular Carcinogenesis and Metastasis. Cancers 2019, 11, 786. [Google Scholar] [CrossRef] [Green Version]
  34. Reticker-Flynn, N.E.; Malta, D.F.B.; Winslow, M.M.; Lamar, J.M.; Xu, M.J.; Underhill, G.H.; Hynes, R.O.; Jacks, T.E.; Bhatia, S.N. A combinatorial extracellular matrix platform identifies cell-extracellular matrix interactions that correlate with metastasis. Nat. Commun. 2012, 3, 1122. [Google Scholar] [CrossRef] [Green Version]
  35. Mevel, R.; Draper, J.E.; Lie-A-Ling, M.; Kouskoff, V.; Lacaud, G. RUNX transcription factors: Orchestrators of development. Development 2019, 146, dev148296. [Google Scholar] [CrossRef] [Green Version]
  36. Lau, Q.C.; Raja, E.; Salto-Tellez, M.; Liu, Q.; Ito, K.; Inoue, M.; Putti, T.C.; Loh, M.; Ko, T.K.; Huang, C.; et al. RUNX3 is frequently inactivated by dual mechanisms of protein mislocalization and promoter hypermethylation in breast cancer. Cancer Res. 2006, 66, 6512–6520. [Google Scholar] [CrossRef] [Green Version]
  37. Kudo, Y.; Tsunematsu, T.; Takata, T. Oncogenic role of RUNX3 in head and neck cancer. J. Cell Biochem. 2011, 112, 387–393. [Google Scholar] [CrossRef]
  38. Lee, Y.S.; Lee, J.W.; Jang, J.W.; Chi, X.Z.; Kim, J.H.; Li, Y.H.; Kim, M.K.; Kim, D.M.; Choi, B.S.; Kim, E.G.; et al. Runx3 inactivation is a crucial early event in the development of lung adenocarcinoma. Cancer Cell 2013, 24, 603–616. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Ito, Y.; Bae, S.C.; Chuang, L.S.H. The RUNX family: Developmental regulators in cancer. Nat. Rev. Cancer 2015, 15, 81–95. [Google Scholar] [CrossRef] [PubMed]
  40. Kitago, M.; Martinez, S.R.; Nakamura, T.; Sim, M.S.; Hoon, D.S. Regulation of RUNX3 tumor suppressor gene expression in cutaneous melanoma. Clin. Cancer Res. 2009, 15, 2988–2994. [Google Scholar] [CrossRef] [Green Version]
  41. Chen, L.F. Tumor suppressor function of RUNX3 in breast cancer. J. Cell Biochem. 2012, 113, 1470–1477. [Google Scholar] [CrossRef] [Green Version]
  42. Manandhar, S.; Lee, Y.M. Emerging role of RUNX3 in the regulation of tumor microenvironment. BMB Rep. 2018, 51, 174–181. [Google Scholar] [CrossRef] [Green Version]
  43. Chen, H.; Crosley, P.; Azad, A.K.; Gupta, N.; Gokul, N.; Xu, Z.; Weinfeld, M.; Postovit, L.M.; Pangas, S.A.; Hitt, M.M.; et al. RUNX3 Promotes the Tumorigenic Phenotype in KGN, a Human Granulosa Cell Tumor-Derived Cell Line. Int. J. Mol. Sci. 2019, 20, 3471. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Date, Y.; Ito, K. Oncogenic RUNX3: A Link between p53 Deficiency and MYC Dysregulation. Mol. Cells 2020, 43, 176–181. [Google Scholar] [PubMed]
  45. Chang, T.L.; Ito, K.; Ko, T.K.; Liu, Q.; Salto-Tellez, M.; Yeoh, K.G.; Fukamachi, H.; Ito, Y. Claudin-1 has tumor suppressive activity and is a direct target of RUNX3 in gastric epithelial cells. Gastroenterology 2010, 138, 255–265.e3. [Google Scholar] [CrossRef] [Green Version]
  46. Voon, D.C.; Wang, H.; Koo, J.K.W.; Nguyen, T.A.P.; Hor, Y.T.; Chu, Y.S.; Ito, K.; Fukamachi, H.; Chan, S.L.; Thiery, J.P.; et al. Runx3 protects gastric epithelial cells against epithelial-mesenchymal transition-induced cellular plasticity and tumorigenicity. Stem Cells 2012, 30, 2088–2099. [Google Scholar] [CrossRef] [PubMed]
  47. Tanaka, S.; Shiraha, H.; Nakanishi, Y.; Nishina, S.; Matsubara, M.; Horiguchi, S.; Takaoka, N.; Iwamuro, M.; Kataoka, J.; Kuwaki, K.; et al. Runt-related transcription factor 3 reverses epithelial-mesenchymal transition in hepatocellular carcinoma. Int. J. Cancer 2012, 131, 2537–2546. [Google Scholar] [CrossRef] [Green Version]
  48. Kulkarni, M.; Tan, T.Z.; Syed Sulaiman, N.B.; Lamar, J.M.; Bansal, P.; Cui, J.; Qiao, Y.; Ito, Y. RUNX1 and RUNX3 protect against YAP-mediated EMT, stem-ness and shorter survival outcomes in breast cancer. Oncotarget 2018, 9, 14175–14192. [Google Scholar] [CrossRef] [Green Version]
  49. Zhang, X.; Wang, L.; Zeng, X.; Fujita, T.; Liu, W. Runx3 inhibits melanoma cell migration through regulation of cell shape change. Cell Biol. Int. 2017, 41, 1048–1055. [Google Scholar] [CrossRef]
  50. Chen, F.; Liu, X.; Bai, J.; Pei, D.; Zheng, J. The emerging role of RUNX3 in cancer metastasis (Review). Oncol Rep. 2016, 35, 1227–1236. [Google Scholar] [CrossRef] [Green Version]
  51. Lotem, J.; Levanon, D.; Negreanu, V.; Bauer, O.; Hantisteanu, S.; Dicken, J.; Groner, Y. Runx3 in Immunity, Inflammation and Cancer. Adv. Exp. Med. Biol. 2017, 962, 369–393. [Google Scholar]
  52. Kim, B.R.; Na, Y.J.; Kim, J.L.; Jeong, Y.A.; Park, S.H.; Jo, M.J.; Jeong, S.; Kang, S.; Oh, S.C.; Lee, D.H. RUNX3 suppresses metastasis and stemness by inhibiting Hedgehog signaling in colorectal cancer. Cell Death Differ. 2020, 27, 676–694. [Google Scholar] [CrossRef]
  53. Cao, W.; Zhang, Z.Y.; Xu, Q.; Sun, Q.; Yan, M.; Zhang, J.; Zhang, P.; Han, Z.G.; Chen, W.T. Epigenetic silencing of MAL, a putative tumor suppressor gene, can contribute to human epithelium cell carcinoma. Mol. Cancer 2010, 9, 296. [Google Scholar] [CrossRef] [Green Version]
  54. Dubash, A.D.; Koetsier, J.L.; Amargo, E.V.; Najor, N.A.; Harmon, R.M.; Green, K.J. The GEF Bcr activates RhoA/MAL signaling to promote keratinocyte differentiation via desmoglein-1. J. Cell Biol. 2013, 202, 653–666. [Google Scholar] [CrossRef] [Green Version]
  55. Kalmar, A.; Peterfia, B.; Hollosi, P.; Galamb, O.; Spisak, S.; Wichmann, B.; Bodor, A.; Toth, K.; Patai, A.V.; Valcz, G.; et al. DNA hypermethylation and decreased mRNA expression of MAL, PRIMA1, PTGDR and SFRP1 in colorectal adenoma and cancer. BMC Cancer 2015, 15, 736. [Google Scholar] [CrossRef] [Green Version]
  56. Beder, L.B.; Gunduz, M.; Hotomi, M.; Fujihara, K.; Shimada, J.; Tamura, S.; Gunduz, E.; Fukushima, K.; Yaykasli, K.; Grenman, R.; et al. T-lymphocyte maturation-associated protein gene as a candidate metastasis suppressor for head and neck squamous cell carcinomas. Cancer Sci. 2009, 100, 873–880. [Google Scholar] [CrossRef]
  57. Fidler, I.J. Selection of successive tumour lines for metastasis. Nat. New Biol. 1973, 242, 148–149. [Google Scholar] [CrossRef]
  58. Nakamura, K.; Yoshikawa, N.; Yamaguchi, Y.; Kagota, S.; Shinozuka, K.; Kunitomo, M. Characterization of mouse melanoma cell lines by their mortal malignancy using an experimental metastatic model. Life Sci. 2002, 70, 791–798. [Google Scholar] [CrossRef]
  59. Chen, H.; Shan, J.; Chen, D.; Wang, R.; Qi, W.; Wang, H.; Ke, Y.; Liu, W.; Zeng, X. CtIP promotes G2/M arrest in etoposide-treated HCT116 cells in a p53-independent manner. J. Cell Physiol. 2019, 234, 11871–11881. [Google Scholar] [CrossRef]
  60. Vasaikar, S.V.; Straub, P.; Wang, J.; Zhang, B. LinkedOmics: Analyzing multi-omics data within and across 32 cancer types. Nucleic Acids Res. 2018, 46, D956–D963. [Google Scholar] [CrossRef] [Green Version]
  61. Nichols, S.E., Jr.; Reams, W.M., Jr. The occurrence and morphogenesis of melanocytes in the connective tissues of the PET/MCV mouse strain. J. Embryol. Exp. Morphol. 1960, 8, 24–32. [Google Scholar]
  62. Cramer, S.F. The origin of epidermal melanocytes. Implications for the histogenesis of nevi and melanomas. Arch. Pathol. Lab. Med. 1991, 115, 115–119. [Google Scholar]
  63. Gong, F.; Guo, Y.; Niu, Y.; Jin, J.; Zhang, X.; Shi, X.; Zhang, L.; Li, R.; Chen, L.; Ma, R.Z. Epigenetic silencing of TET2 and TET3 induces an EMT-like process in melanoma. Oncotarget 2017, 8, 315–328. [Google Scholar] [CrossRef] [Green Version]
  64. Shi, X.; Deepak, V.; Wang, L.; Ba, X.; Komori, T.; Zeng, X.; Liu, W. Thrombospondin-1 is a putative target gene of Runx2 and Runx3. Int. J. Mol. Sci. 2013, 14, 14321–14332. [Google Scholar] [CrossRef] [Green Version]
  65. Murphy-Ullrich, J.E.; Hook, M. Thrombospondin modulates focal adhesions in endothelial cells. J. Cell Biol. 1989, 109, 1309–1319. [Google Scholar] [CrossRef]
  66. Goodenow, R.S.; Vogel, J.M.; Linsk, R.L. Histocompatibility antigens on murine tumors. Science 1985, 230, 777–783. [Google Scholar] [CrossRef]
  67. Mimori, K.; Shiraishi, T.; Mashino, K.; Sonoda, H.; Yamashita, K.; Yoshinaga, K.; Masuda, T.; Utsunomiya, T.; Alonso, M.A.; Inoue, H.; et al. MAL gene expression in esophageal cancer suppresses motility, invasion and tumorigenicity and enhances apoptosis through the Fas pathway. Oncogene 2003, 22, 3463–3471. [Google Scholar] [CrossRef] [Green Version]
  68. Kurashige, J.; Sawada, G.; Takahashi, Y.; Eguchi, H.; Sudo, T.; Ikegami, T.; Yoshizumi, T.; Soejima, Y.; Ikeda, T.; Kawanaka, H.; et al. Suppression of MAL gene expression in gastric cancer correlates with metastasis and mortality. Fukuoka Igaku Zasshi 2013, 104, 344–349. [Google Scholar]
  69. Fidler, I.J.; Radinsky, R. Search for genes that suppress cancer metastasis. J. Natl. Cancer Inst. 1996, 88, 1700–1703. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Pascual, G.; Avgustinova, A.; Mejetta, S.; Martin, M.; Castellanos, A.; Attolini, C.S.O.; Berenguer, A.; Prats, N.; Toll, A.; Hueto, J.A.; et al. Targeting metastasis-initiating cells through the fatty acid receptor CD36. Nature 2017, 541, 41–45. [Google Scholar] [CrossRef] [PubMed]
  71. Xu, L.; Shen, S.S.; Hoshida, Y.; Subramanian, A.; Ross, K.; Brunet, J.P.; Wagner, S.N.; Ramaswamy, S.; Mesirov, J.P.; Hynes, R.O. Gene expression changes in an animal melanoma model correlate with aggressiveness of human melanoma metastases. Mol. Cancer Res. 2008, 6, 760–769. [Google Scholar] [CrossRef] [Green Version]
  72. Hynes, R.O.; Yamada, K.M. Fibronectins: Multifunctional modular glycoproteins. J. Cell Biol. 1982, 95, 369–377. [Google Scholar] [CrossRef] [Green Version]
  73. Yoshino, T.; Saito, D.; Atsuta, Y.; Uchiyama, C.; Ueda, S.; Sekiguchi, K.; Takahashi, Y. Interepithelial signaling with nephric duct is required for the formation of overlying coelomic epithelial cell sheet. Proc. Natl. Acad. Sci. USA 2014, 111, 6660–6665. [Google Scholar] [CrossRef] [Green Version]
Ijms 22 02219 g001 550
Figure 1. Runt-related transcription factor-3 (Runx3) re-expression in B16-F10 melanoma cells resulted in a cell shape change and suppressed the pulmonary metastasis of the cells. (A) Runx3 changed the cell shape and altered the actin cytoskeleton. The green fluorescent signals indicate actin filaments. (B) Runx3 delayed cell migration, and altered the dynamics of stress fiber formation in the migrating cells. Start denotes when the wound is just made, and end denotes when wound healing is terminated. The green fluorescent signals indicate actin filaments. The arrows indicate the directions of migration fronts. The arrowheads indicate the ends of stress fibers. The thin arrows indicate microspike-like structures. (C) Runx3 completely suppressed pulmonary metastasis. n = 4. The black dots on the surface of the lungs indicate metastasis foci. The windows-arrows representatively indicate metastasis foci. (D) Runx3 slightly reduced the tumor formation rate. The averaged tumor volumes are compared to each other (mean ± SEM, n = 3). * p < 0.05 is generated by the t-test. A representative image is shown. Bar: 1 cm. (E) Runx3 changed the cell shape in vivo. A representative section is shown. Ctrl: mock control B16-F10 cells; Runx3: B16-F10/Runx3 cells. All bars: 22 μm.
Figure 1. Runt-related transcription factor-3 (Runx3) re-expression in B16-F10 melanoma cells resulted in a cell shape change and suppressed the pulmonary metastasis of the cells. (A) Runx3 changed the cell shape and altered the actin cytoskeleton. The green fluorescent signals indicate actin filaments. (B) Runx3 delayed cell migration, and altered the dynamics of stress fiber formation in the migrating cells. Start denotes when the wound is just made, and end denotes when wound healing is terminated. The green fluorescent signals indicate actin filaments. The arrows indicate the directions of migration fronts. The arrowheads indicate the ends of stress fibers. The thin arrows indicate microspike-like structures. (C) Runx3 completely suppressed pulmonary metastasis. n = 4. The black dots on the surface of the lungs indicate metastasis foci. The windows-arrows representatively indicate metastasis foci. (D) Runx3 slightly reduced the tumor formation rate. The averaged tumor volumes are compared to each other (mean ± SEM, n = 3). * p < 0.05 is generated by the t-test. A representative image is shown. Bar: 1 cm. (E) Runx3 changed the cell shape in vivo. A representative section is shown. Ctrl: mock control B16-F10 cells; Runx3: B16-F10/Runx3 cells. All bars: 22 μm.
Ijms 22 02219 g001
Ijms 22 02219 g002 550
Figure 2. Microarray analysis uncovered an altered and specific transcriptional profile underlying the cell shape change and the suppression of the metastatic potential by Runx3. (A) Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotations assigned different repertoires of differentially expressed genes (DEGs) that were regulated by Runx3 in B16-F10 cells and, in the case of B16-F0 cells vs. B16-F10 cells, to the terms relevant to the actin cytoskeleton and adhesion. Each GO/KEGG annotation term is inscribed as a caption. The number of DEGs within a repertoire is described at x axis. The columns for the downregulated and the upregulated DEGs are shown in different colors as explained at the top of the figure., Down/Up by Runx3 indicates the appearance of DEGs in the case of B16-F10/Runx3 cells vs. mock control B16-F10 cells. Down/Up in B16-F0 indicates the appearance of DEGs in the case of B16-F0 cells vs. B16-F10 cells. Down/Up in common indicates the appearance of same DEGs in both cases. (B) Runx3 upregulated the expression of extracellular matrix (ECM) genes. The summary table indicates the gene names and their corresponding fold changes. (C) Runx3 regulated the expression of a list of DEGs that were inversely associated with an increase in the metastatic potential of B16-F10 cells compared to B16-F0 cells. The heatmap shows the gene names and their fold changes in the cases of B16-F10/Runx3 cells vs. mock control B16-F10 cells and B16-F0 cells vs. B16-F10 cells, respectively.
Figure 2. Microarray analysis uncovered an altered and specific transcriptional profile underlying the cell shape change and the suppression of the metastatic potential by Runx3. (A) Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotations assigned different repertoires of differentially expressed genes (DEGs) that were regulated by Runx3 in B16-F10 cells and, in the case of B16-F0 cells vs. B16-F10 cells, to the terms relevant to the actin cytoskeleton and adhesion. Each GO/KEGG annotation term is inscribed as a caption. The number of DEGs within a repertoire is described at x axis. The columns for the downregulated and the upregulated DEGs are shown in different colors as explained at the top of the figure., Down/Up by Runx3 indicates the appearance of DEGs in the case of B16-F10/Runx3 cells vs. mock control B16-F10 cells. Down/Up in B16-F0 indicates the appearance of DEGs in the case of B16-F0 cells vs. B16-F10 cells. Down/Up in common indicates the appearance of same DEGs in both cases. (B) Runx3 upregulated the expression of extracellular matrix (ECM) genes. The summary table indicates the gene names and their corresponding fold changes. (C) Runx3 regulated the expression of a list of DEGs that were inversely associated with an increase in the metastatic potential of B16-F10 cells compared to B16-F0 cells. The heatmap shows the gene names and their fold changes in the cases of B16-F10/Runx3 cells vs. mock control B16-F10 cells and B16-F0 cells vs. B16-F10 cells, respectively.
Ijms 22 02219 g002
Ijms 22 02219 g003 550
Figure 3. The gene expression of myelin and lymphocyte protein (MAL) and/or vilin-like (VILL) had prognostic value for various cancers. (A) Kaplan–Meier survival analysis implicated MAL or VILL in various cancers. (B) Kaplan–Meier survival analysis implicated both MAL and VILL in various cancers. The information about cancer type, patient number, and p-value is inscribed at each panel. The gene level was defined as fragments per kilobase of exon per million reads (FPKM). The median cutoff was used to group patients into low (L) and high (H) expression of MAL and/or VILL. The p-values are generated by the log-rank test.
Figure 3. The gene expression of myelin and lymphocyte protein (MAL) and/or vilin-like (VILL) had prognostic value for various cancers. (A) Kaplan–Meier survival analysis implicated MAL or VILL in various cancers. (B) Kaplan–Meier survival analysis implicated both MAL and VILL in various cancers. The information about cancer type, patient number, and p-value is inscribed at each panel. The gene level was defined as fragments per kilobase of exon per million reads (FPKM). The median cutoff was used to group patients into low (L) and high (H) expression of MAL and/or VILL. The p-values are generated by the log-rank test.
Ijms 22 02219 g003
Ijms 22 02219 g004 550
Figure 4. Mal gene expression was correlated with metastatic potential among the cell lines and was subject to epigenetic regulation. (A) Mal gene expression was validated in the cells used in this study with qRT-PCR. The relative Mal gene expression levels were compared among the indicated cell lines (means ± SEM, n = 3). (B) Treatment with histone deacetylase (HDAC) inhibitors induced the gene expression of Mal in B16-F10 cells. The relative Mal gene expression levels of inhibitor-treated cells were compared to those of untreated (DMSO) cells (means ± SEM, n = 3). (C) The gene expression of MAL was correlated with that of other genes in human melanoma. The top left panel indicates the overall MAL correlation. The other panels indicate the correlation between MAL and the other genes. The information about patient number and p-value is inscribed at each panel. The gene level at x/y axis is defined from FPKM. The p-values are generated by the Spearman correlation test. * p < 0.05 and ** p < 0.01 are generated by the t-test.
Figure 4. Mal gene expression was correlated with metastatic potential among the cell lines and was subject to epigenetic regulation. (A) Mal gene expression was validated in the cells used in this study with qRT-PCR. The relative Mal gene expression levels were compared among the indicated cell lines (means ± SEM, n = 3). (B) Treatment with histone deacetylase (HDAC) inhibitors induced the gene expression of Mal in B16-F10 cells. The relative Mal gene expression levels of inhibitor-treated cells were compared to those of untreated (DMSO) cells (means ± SEM, n = 3). (C) The gene expression of MAL was correlated with that of other genes in human melanoma. The top left panel indicates the overall MAL correlation. The other panels indicate the correlation between MAL and the other genes. The information about patient number and p-value is inscribed at each panel. The gene level at x/y axis is defined from FPKM. The p-values are generated by the Spearman correlation test. * p < 0.05 and ** p < 0.01 are generated by the t-test.
Ijms 22 02219 g004
Ijms 22 02219 g005 550
Figure 5. Mal was necessary for suppressing the migratory and invasive traits of metastatic melanoma cells. (A) The knockdown of Mal gene expression was validated with qRT-PCR and Western blotting. The relative Mal gene expression levels are compared to each other (means ± SEM, n = 3). Mal’s molecular weight: 18 kDa; actin’s molecular weight: 43 kDa. (B) The knockdown of Mal gene expression changed the cell shape and altered the actin cytoskeleton of steady or migrating cells. The green fluorescent signals indicate actin filaments. The arrows indicate the directions of migration fronts. The thin arrows indicate microspike-like structures. Bars: 22 μm. (C) The knockdown of Mal gene expression slightly reduced the proliferation rate. The relative proliferation rates were determined by 3-(4,5-dimethyl-2-thiazolyl)-2,5-diphenyl-2-H-tetrazolium bromide (MTT) and were compared to each other (means ± SEM, n = 3). (D) The knockdown of Mal gene expression slightly delayed the progress of the cell cycle. The percentages of cells at various phases (G0/G1, S, and G2/M) were determined by flow cytometry (means ± SEM, n = 3, ns: not significant). (E) The knockdown of Mal gene expression had no effect on the tumor formation rate. The averaged tumor volumes were compared to each other (means ± SEM, n = 4, ns: not significant). (F) The knockdown of Mal gene expression increased the migration rate in wound healing. The relative scratch open areas were compared to each other (means ± SEM, n = 3). (G) The knockdown of Mal gene expression increased the migration rate through the Transwell chambers. A representative image is shown on the left of the panel, and the averaged cell numbers were compared to each other, as shown on the right (means ± SEM, n = 6). (H) The knockdown of Mal gene expression increased the invasion rate through Matrigel-coated Transwell chambers. A representative image is shown on the left of the panel, and the averaged cell numbers were compared to each other, as shown on the right (means ± SEM, n = 6). (I) The knockdown of Mal gene expression increased the frequency of metastasis-focus-positive lungs. A representative image is shown. The arrowheads indicate metastasis foci. The frequencies are inscribed at the top of the image. Ctrl: mock control; shMal: short hairpin RNA interference of Mal. * p < 0.05 and ** p < 0.01 are generated by the t-test.
Figure 5. Mal was necessary for suppressing the migratory and invasive traits of metastatic melanoma cells. (A) The knockdown of Mal gene expression was validated with qRT-PCR and Western blotting. The relative Mal gene expression levels are compared to each other (means ± SEM, n = 3). Mal’s molecular weight: 18 kDa; actin’s molecular weight: 43 kDa. (B) The knockdown of Mal gene expression changed the cell shape and altered the actin cytoskeleton of steady or migrating cells. The green fluorescent signals indicate actin filaments. The arrows indicate the directions of migration fronts. The thin arrows indicate microspike-like structures. Bars: 22 μm. (C) The knockdown of Mal gene expression slightly reduced the proliferation rate. The relative proliferation rates were determined by 3-(4,5-dimethyl-2-thiazolyl)-2,5-diphenyl-2-H-tetrazolium bromide (MTT) and were compared to each other (means ± SEM, n = 3). (D) The knockdown of Mal gene expression slightly delayed the progress of the cell cycle. The percentages of cells at various phases (G0/G1, S, and G2/M) were determined by flow cytometry (means ± SEM, n = 3, ns: not significant). (E) The knockdown of Mal gene expression had no effect on the tumor formation rate. The averaged tumor volumes were compared to each other (means ± SEM, n = 4, ns: not significant). (F) The knockdown of Mal gene expression increased the migration rate in wound healing. The relative scratch open areas were compared to each other (means ± SEM, n = 3). (G) The knockdown of Mal gene expression increased the migration rate through the Transwell chambers. A representative image is shown on the left of the panel, and the averaged cell numbers were compared to each other, as shown on the right (means ± SEM, n = 6). (H) The knockdown of Mal gene expression increased the invasion rate through Matrigel-coated Transwell chambers. A representative image is shown on the left of the panel, and the averaged cell numbers were compared to each other, as shown on the right (means ± SEM, n = 6). (I) The knockdown of Mal gene expression increased the frequency of metastasis-focus-positive lungs. A representative image is shown. The arrowheads indicate metastasis foci. The frequencies are inscribed at the top of the image. Ctrl: mock control; shMal: short hairpin RNA interference of Mal. * p < 0.05 and ** p < 0.01 are generated by the t-test.
Ijms 22 02219 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, N.; Zhang, H.; Cui, X.; Ma, C.; Wang, L.; Liu, W. Runx3 Induces a Cell Shape Change and Suppresses Migration and Metastasis of Melanoma Cells by Altering a Transcriptional Profile. Int. J. Mol. Sci. 2021, 22, 2219. https://doi.org/10.3390/ijms22042219

AMA Style

Wang N, Zhang H, Cui X, Ma C, Wang L, Liu W. Runx3 Induces a Cell Shape Change and Suppresses Migration and Metastasis of Melanoma Cells by Altering a Transcriptional Profile. International Journal of Molecular Sciences. 2021; 22(4):2219. https://doi.org/10.3390/ijms22042219

Chicago/Turabian Style

Wang, Ning, Haiying Zhang, Xiulin Cui, Chao Ma, Linghui Wang, and Wenguang Liu. 2021. "Runx3 Induces a Cell Shape Change and Suppresses Migration and Metastasis of Melanoma Cells by Altering a Transcriptional Profile" International Journal of Molecular Sciences 22, no. 4: 2219. https://doi.org/10.3390/ijms22042219

APA Style

Wang, N., Zhang, H., Cui, X., Ma, C., Wang, L., & Liu, W. (2021). Runx3 Induces a Cell Shape Change and Suppresses Migration and Metastasis of Melanoma Cells by Altering a Transcriptional Profile. International Journal of Molecular Sciences, 22(4), 2219. https://doi.org/10.3390/ijms22042219

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