Research | Open | Published:
Analysis of microRNA (miRNA) expression profiles reveals 11 key biomarkers associated with non-small cell lung cancer
World Journal of Surgical Oncologyvolume 15, Article number: 175 (2017)
Non-small cell lung cancer (NSCLC) accounts for more than 85% of lung cancer cases which cause most of cancer-related deaths globally. However, the results vary largely in different studies due to different platforms and sample sizes. Here, we aim to identify the key miRNAs in the carcinogenesis of NSCLC that might be potential biomarkers for this cancer.
Meta-analysis was performed on miRNA profile using seven datasets of NSCLC studies. Furthermore, we predicted and investigated the functions of genes regulated by key miRNAs.
Eleven key miRNAs were identified, including 2 significantly upregulated ones (hsa-miR-21-5p and hsa-miR-233-3p) and 9 downregulated ones (hsa-miR-126-3p, hsa-miR-133a-3p, hsa-miR-140-5p, hsa-miR-143-5p, hsa-miR-145-5p, hsa-miR-30a-5p, hsa-miR-30d-3p, hsa-miR-328-3pn, and hsa-miR-451). The functional enrichment analysis revealed that both up- and downregulated miRNAs were proportionally associated with regulation of transcription from RNA polymerase II promoter. According to transcription factor analysis, there were 65 (43.9%) transcription factors influenced by both up- and downregulated miRNAs.
In this study, 11 meta-signature miRNAs, as well as their target genes and transcription factors, were found to play significant role in carcinogenesis of NSCLC. These target genes identified in our study may be profitable to diagnosis and prognostic prediction of NSCLC as biomarkers.
Although great effort has been put in research and prevention strategies, lung cancer is still the most common cancer globally . To date, it even shows increasing tendency in both annual incidence and mortality . Non-small cell lung cancer (NSCLC) accounts for almost 85–90% of all lung cancers . Additionally, lung cancer remains to be the most common cause giving rise to cancer-related deaths worldwide. In the USA, one 5-year survival rate is shown to be less than 20% in patients . Usually, the occurrence and development of NSCLC results from multistep carcinogenesis, and these involve in various changes of gene expression and signal transduction pathways [5, 6]. Nevertheless, there are still unclear mechanisms, which help to promote carcinogenesis of NSCLC, that need to be exploited. Hence, a search of novel biomarkers for the diagnosis and prognostic prediction of NSCLC is critical for the patients in order to receive proper therapies as early as possible.
MicroRNAs (miRNAs) are a type of small non-coding RNA molecules, which acts in RNA silencing and post-transcriptionally regulating gene expressions [7, 8]. Kentaro Inamura and Yuichi Ishikawa summarized the characteristic nature of miRNAs which can be used for diagnostics, prognostics, and targeted therapeutics in tumor : (1) the easy availability of miRNAs makes them promising biomarkers for non-invasive liquid biopsy in cancer diagnosis; (2) the stability of miRNAs in formalin-fixed paraffin-embedded (FFPE) samples guarantees the accuracy of measurements using miRNAs which can be performed from FFPE samples collected and preserved in hospitals. Indeed, miRNAs are novel and suitable diagnostic markers to be applied in histological classification or genetic alterations .
Recently, increasing studies reporting various expression levels of miRNAs in different cancers (including NSCLC) suggests important roles of certain miRNAs in carcinogenesis [5, 10, 11]. For instance, several miRNAs, such as miRNA-224 and miRNA-30d-5p, have been reported to function either as cancer promoters or as suppressors in NSCLCs [12, 13]. Additionally, previous studies have revealed that the key biological functions of some miRNAs, such as miRNA-107, miRNA-93, and miRNA-148b, are closely associated with clinicopathological variables and to display a diagnostic and prognostic property in NSCLC [14,15,16]. Even though more and more miRNAs are identified, there are still abundant undegraded miRNAs in tissues as well as in body fluids (e.g., blood, plasma, serum, and sputum) [17, 18]. However, several factors usually give rise to common drawbacks of recent miRNA expression profiling studies resulting from a lack of agreement: inconsistent annotation, application of different technological variations, different methods on data processing and analysis, different sample sizes, and ongoing discovery of novel miRNAs [19, 20].
In order to minimize the drawbacks, meta-analysis approach was applied in this study, which combined several independent studies to increase the statistical power, to identify the key miRNAs in carcinogenesis of NSCLC. Furthermore, we investigated the functions of the target genes and the transcription factors regulated by these miRNAs. Our study provided potential target genes that may be profitable to early detection and precise diagnosis of NSCLC by using miRNAs identified here.
Selection of studies and datasets
NSCLC miRNA expression profiling studies were searched and collected using PubMed databases with terms “non-small cell lung cancer” and “miRNA or microRNA.” Articles for meta-analysis since 2000 were selected as the following criteria. Full text of every study was carefully evaluated. Original experimental literatures that provide detailed miRNA expressions between NSCLC and non-cancerous control in human were selected for further analysis. Expression profile studies of individual preselected candidate genes and studies performed with only cell lines were removed. Studies focusing on different histologic subtypes without using non-cancerous tissue were dislodged as well.
The miRNAs listed by statistically significant expression variations between NSCLC and non-cancerous control in human were extracted from the published papers. Meanwhile, in case the gene list was not available, we also contacted the authors directly. All miRNA names were standardized by miRBase version 21. Many traditional “major” miRNA annotations throughout the main text were re-designated according to miRBase database version 21. Viral non-miRNA and miRNA probes were excluded from the analysis. After the standardization of precursor names, the pre-miRNAs, which are reported in certain studies, were used in the analyses.
Cluster analysis of datasets
To evaluate the correlations between the outcomes of individual studies, we performed hierarchical cluster analysis. Briefly, overall rank matrix was constructed according to rank matrixes generated from separate analyses for upregulated and downregulated gene with RobustRankAggreg of R package. In the matrix, value 0.5 indicates that this miRNA was not reported in that study, while value above 0.5 indicates that it is upregulated (one minus the normalized rank of miRNA on the upregulated gene lists) and value below 0.5 indicates that this miRNA is downregulated in that study (normalized rank on the downregulated gene lists). Spearman rank correlation with average linkage method was used in cluster analysis.
The extracted miRNAs were prioritized basing on the p values in t test. p < 0.05 was considered as threshold for significance identification. For a better consistence of the identified miRNAs, we applied a novel robust rank aggregation method—RobustRankAggreg of R package . In this method, one p value to each element was assigned to indicate how prior it is when compared to a null model. In order to evaluate the stability of acquired p values, we left one out cross-validation on the robust rank aggregation algorithm and repeated analyses for 10,000 times. Then, one random gene list was taken out from the analysis for each time. At the end, average of acquired p values from each repeat for each miRNA was calculated.
miRNA target prediction
The potential targets of meta-signature miRNAs were predicted from databases by utilizing target prediction algorithms, TargetScan v6.1 . Considering the experimental threshold of target genes below which protein production is highly repressed together with the prediction algorithms in TargetScan v6.1, the maximum number of target genes for each miRNA with significant expression variation was set to 300. But the other parameters were as default.
A functional enrichment analysis was performed to examine the enrichment of annotated terms. Gene ontology (GO) enrichment analysis was conducted via using the Database for Annotation Visualization and Integrated Discovery (DAVID). Predicted target genes for each miRNA were used as input and followed by FDR-corrected p values (p < 0.05) and then visualized as a heatmap. Clustering of the heatmap was according to Pearson correlation and average linkage.
Transcription factor analysis
The key up- and downregulated miRNAs were uploaded to Tfacts database (http://www.tfacts.org/) [22, 23]. The transcription factors were predicted by p value, q value, E value, and FDR (each value is less than 0.05). The transcription factors regulated by up- and downregulated miRNAs were compared to find out the overlapped ones.
Here, we applied meta-analysis to identify the key significantly up- and downregulated miRNAs from the seven available NSCLC miRNA datasets (Table 1) [24,25,26,27,28,29]. Furthermore, we also predicted potential target genes of the selected miRNAs, of which function annotation and transcription factor (TF) analysis were performed as well (Fig. 1).
Selection of miRNAs from the datasets
According to search criteria, 7 NSCLC miRNA expression profiling datasets from 7 published literatures were gained from public databases. A short description of these studies, such as the assay type, sample numbers, and number of miRNA probes, are provided in Table 1. Although the number of patients investigated was ranging from 6 to 434 across the studies, qRT-PCR was applied on each dataset to verify the reliability of expression difference of miRNAs. Expansion of studied miRnome from 7 datasets was reflected by Scalable Vector Graphics (SVG) of the differentially expressed miRNAs in all datasets. We found that generally, the miRNA expression profiles from different studies had huge variations (Fig. 2).
The size of significantly deregulated miRNA varies largely across the 7 studies. Three datasets showed more than 20 deregulated miRNAs, namely, dataset 1 (38 deregulated miRNAs), dataset 2 (29 deregulated miRNAs), and dataset 6 (40 deregulated miRNAs). Dataset 6 has the most upregulated miRNAs (27 upregulated miRNAs), while dataset 1 has the most downregulated miRNAs (18 downregulated miRNAs) (Fig. 3a).
NSCLC miRNA meta-signature
In total, based on the 7 independent datasets we selected, 65 (63 in 1 study, 1 and 1 in 2 and 3 study) were reported as significantly upregulated while 52 (42 in 1 study, 9 and 1 in 2 and 3 study) as significantly downregulated (Fig. 3b). However, the dataset-specific upregulated miRNAs were 63 (96.9% of non-redundant upregulated miRNAs) and dataset-specific downregulated miRNAs were 43 (82.7% of non-redundant downregulated miRNAs). The high percentages emphasized the importance of meta-analysis on expression difference of miRNAs in NSCLC. Using robust rank aggregation, we found a statistically significant meta-signature consisting of 2 upregulated and 9 downregulated miRNAs in NSCLC samples compared to non-cancerous tissue (Table 2). These 11 meta-signature miRNAs that displayed highly statistical significance were reported by at least 2 datasets and were considered as key miRNAs in carcinogenesis of NSCLC.
According to miRBase, meta-signature miRNA genes were located on different chromosomal locations. Chromosome 5, chromosome 16, and chromosome 17 all harbored 2 different miRNAs; meanwhile, chromosomes 6, 8, 9, and 18 severally harbored one miRNAs (Table 2).
Target prediction and functional enrichment of meta-signature miRNAs
The high-stringency target gene prediction for meta-signature miRNAs was conducted by TargetScan. We obtained target genes from experimentally supported databases and prediction algorithms. The distribution of target genes is summarized in Additional file 1: Table S1. All the meta-signature miRNAs had the same number of target genes (300). Briefly, 527 non-redundant target genes were affected by the 2 upregulated miRNAs, and 1882 non-redundant target genes were affected by the 9 downregulated genes. The number of target genes affected by each miRNA varied largely. The numbers of target genes affected by 3 miRNAs, namely, hsa-miR-126-3p, hsa-miR-328-3p, and hsa-miR-451, were less than 300, which means these miRNAs may be less conserved in different species .
To elucidate the biological function of meta-signature miRNAs, enrichment analyses were conducted via predicted target genes by using the gene ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) signal pathway analysis on The Database for Annotation, Visualization and Integrated Discovery (DAVID; v6.8). KEGG pathway enrichment analyses were concordant for all the 11 meta-signature miRNAs (Fig. 4). Beside hsa-miR-126-3p, other 10 meta-signature miRNAs generated rich GO terms, especially hsa-miR-21-5p which had the richest GO terms.
Enrichment analysis of the 11 miRNA target sets showed that target genes of upregulated miRNAs were most frequently associated with regulation of transcription from RNA polymerase II promoter, while the targets of downregulated were enriched in positive regulation of transcription from RNA polymerase II promoter as well. Pathways in small GTPase-mediated signal transduction and regulation of branching involved in ureteric bud morphogenesis were also influenced by downregulated miRNAs (Table 3).
We further analyzed the pathway enrichment of target genes influenced by the miRNAs and found that the target genes regulated by upregulated miRNAs were mainly enriched in pathways in cancer, proteoglycans in cancer, mitogen-activated protein kinases (MAPK) signaling pathway, Ras signaling pathway, and signaling pathways regulating pluripotency of stem cells (Fig. 5a). Additionally, the target genes regulated by downregulated miRNAs were mainly enriched in four groups of pathway, which are endocytosis, actin cytoskeleton, Hippo signaling (Salvador/Warts/Hippo (SWH) pathway), and bacterial invasion of epithelial cells (Fig. 5b).
Transcription factor analysis
We further investigated the transcription factors regulated by meta-signature miRNAs and compared the difference of transcription factors between up- and downregulated miRNAs (Additional file 2: Table S2). Transcription factor analysis revealed that 195 interactions were formed between 83 transcription factors and 2 upregulated miRNAs, while 633 interactions were formed between 130 transcription factors and 9 downregulated miRNAs. As shown in the Venn diagram, 65 (43.9%) transcription factors were influenced by both up- and downregulated miRNAs (Fig. 6a). Then, the E value scores and intersection percentage for significantly regulated TFs were displayed in detail (Fig. 6b).
Studies analyzing miRNA expressions depending on the evolution of lung cancer have been carried out and reviewed . It is obvious that miRNAs may be potential biomarkers due to their increased specificity and selectivity. However, there is still a need to assign and validate the key miRNAs for different types of lung cancer. In this study, the preferred method we used for gene expression meta-analysis is involved in analysis of raw miRNA expression datasets; however, the unavailability of the raw data often limits the possibility of such rigorous approach. Moreover, changes in the number of miRNAs known at different moment and the technological platform applied in any individual study complicate the integration of raw dataset.
In addition, noisiness of microarray data and the relatively small sample size have caused inconsistency of biological conclusions. Therefore, to get rid of these limitations above, we directly analyzed 7 prioritized miRNA datasets from the 7 published studies and then followed by robust rank aggregation method.
In total, 11 meta-signature miRNAs were identified to play essential roles in carcinogenesis of NSCLC, namely, hsa-miR-21-5p, hsa-miR-223-3p, hsa-miR-126-3p, hsa-miR-140-5p, hsa-miR-133a-3p, hsa-miR-143-5p, hsa-miR-145-5p hsa-miR-30a-5p hsa-miR-30d-3p hsa-miR-328-3p hsa-miR-451, and hsa-miR-140-5p, of which the first 2 showed upregulation and the others displayed downregulated expression. Interestingly, besides the publications we applied in this study, several studies have also reported the importance of these miRNAs not only in studies regarding NSCLC but also in other kinds of cancers. For instance, the overexpression of hsa-miR-21-5p, one of the best studied miRNAs, not only was suggested to be an important prognostic factor for the overall survival of NSCLC patients [31, 32] but also was deregulated in other cancers, such as breast cancer and cervical cancer; it seems to be a common oncogenic miRNA of no tissue specificity [33, 34]. Hsa-miR-126-3p is also reported to display a tumor-suppressive effect on breast cancer cells and be associated with aggressive disease phenotype . Hsa-miR-143-5p exhibits a strong tumor suppressive effect and has been shown to act as anti-oncomirs in gastric cancer and directly target prostaglandin-endoperoxide synthase 2 (COX-2) . Hsa-miR-145-5p also functions as tumor suppressor in gall bladder cancer . Therefore, these meta-signature miRNAs showed great potential to benefit not only NSCLC but also other cancer diagnoses as biomarkers.
The biological function of each meta-signature miRNA was also investigated in our study. One single miRNA may affect expression of multiple target genes; on the other hand, a specific mRNA may be regulated by certain different miRNAs. Subsequently, the miRNAs may induce changes in various processes and pathways and present a further level of mechanism by which NSCLC may be induced. Enriched KEGG showed that target genes of upregulated miRNAs were most frequently associated with regulation of transcription from RNA polymerase II promoter and MAPK pathway. For example, hsa-miR-21, a well-known onco-miRNA, acts on various target genes such as programmed cell death protein 4 (PDCD4)  and PTEN . Has-miR-21 has also been reported to influence the tumor formation by regulating Ras-MEK/ERK pathway. The MAPK pathway plays an essential role in lung cancer; the inhibition of MAPK by flavanone was reported to inhibit the metastasis of lung cancer cells . Meanwhile, the targets of downregulated were enriched in positive regulation of transcription from RNA polymerase II promoter, as well as pathways in small GTPase-mediated signal transduction and regulation of branching involved in ureteric bud morphogenesis and Hippo pathway. The Hippo pathway has an important role in tumorigenesis; in mammals, the regulated Hippo pathway can even promote the lung cancer metastasis . In our result, both up- and downregulated miRNAs were found associated with regulation of transcription from the RNA polymerase II promoter. As polymerase II promoter catalyzes the transcription of DNA to synthesize precursors of mRNA and most snRNA and microRNA in the nucleus of eukaryotic cells , it is still not easy to draw a conclusion clearly that one of the miRNAs will regulate another one by influencing this promoter. Still, there is a large part of studies regarding the roles of the pathways we identified in NSCLC that needs to be further explored.
Currently, our analysis is restricted to comparison of cancerous and non-cancerous tissue only; however, the 11 most frequently and significantly reported differentially expressed miRNAs could be considered as potential diagnostic or/and prognostic biomarkers. Furthermore, functional enrichment analysis and target predictions may also provide a thread for clarifying the functions of miRNAs in the precise underlying mechanisms and the carcinogenesis of NSCLC. Additionally, from clinical perspectives, sufficient sensitivity and specificity of miRNAs should be further detected in the well-designed clinical studies. Therefore, the findings of this study may provide substantially clinical implications or significance.
In conclusion, we utilized a meta-analysis method to critically converge and dissect the heterogeneous miRNA expression profiling datasets regarding NSCLC. In total, 11 meta-signature miRNAs were identified, which are consistently up- and downregulated, and highly significant across 7 different studies. Due to the development of miRNA-based tests, our analysis achieves certain highlights on identifying new biomarkers for NSCLC, but still, it is in need that rigorous evaluation of the results should be carried out prior to proceeding to clinical treatments.
Jemal A, Thun M, Yu XQ, Hartman AM, Cokkinides V, Center MM, Ross H, Ward EM. Changes in smoking prevalence among U.S. adults by state and region: estimates from the tobacco use supplement to the current population survey, 1992-2007. BMC Public Health. 2011;11:512.
Field JK, Hansell DM, Duffy SW, Baldwin DR. CT screening for lung cancer: countdown to implementation. Lancet Oncol. 2013;14:e591–600.
Jemal A, Bray F, Center MM, Ferlay J, Ward E, Forman D. Global cancer statistics. CA Cancer J Clin. 2011;61:69–90.
Molina JR, Yang P, Cassivi SD, Schild SE, Adjei AA. Non-small cell lung cancer: epidemiology, risk factors, treatment, and survivorship. Mayo Clin Proc. 2008;83:584.
Markou A, Sourvinou I, Vorkas PA, Yousef GM, Lianidou E. Clinical evaluation of microRNA expression profiling in non small cell lung cancer. Lung Cancer. 2013;81:388–96.
Kalari KR, Rossell D, Necela BM, Asmann YW, Nair A, Baheti S, Kachergus JM, Younkin CS, Baker T, Carr JM, et al. Deep sequence analysis of non-small cell lung cancer: integrated analysis of gene expression, alternative splicing, and single nucleotide variations in lung adenocarcinomas with and without oncogenic KRAS mutations. Front Oncol. 2012;2:12.
Yanaihara N, Caplen N, Bowman E, Seike M, Kumamoto K, Yi M, Stephens RM, Okamoto A, Yokota J, Tanaka T, et al. Unique microRNA molecular profiles in lung cancer diagnosis and prognosis. Cancer Cell. 2006;9:189–98.
Landi MT, Zhao Y, Rotunno M, Koshiol J, Liu H, Bergen AW, Rubagotti M, Goldstein AM, Linnoila I, Marincola FM, et al. MicroRNA expression differentiates histology and predicts survival of lung cancer. Clin Cancer Res. 2010;16:430–41.
Inamura K, Ishikawa Y. MicroRNA in lung cancer: novel biomarkers and potential tools for treatment. J Clin Med. 2016;5:36.
Liu ZL, Wang H, Liu J, Wang ZX. MicroRNA-21 (miR-21) expression promotes growth, metastasis, and chemo- or radioresistance in non-small cell lung cancer cells by targeting PTEN. Mol Cell Biochem. 2013;372:35–45.
Vosa U, Vooder T, Kolde R, Vilo J, Metspalu A, Annilo T. Meta-analysis of microRNA expression in lung cancer. Int J Cancer. 2013;132:2884–93.
Cui R, Meng W, Sun HL, Kim T, Ye Z, Fassan M, Jeon YJ, Li B, Vicentini C, Peng Y, et al. MicroRNA-224 promotes tumor progression in nonsmall cell lung cancer. Proc Natl Acad Sci U S A. 2015;112:E4288–97.
Chen D, Guo W, Qiu Z, Wang Q, Li Y, Liang L, Liu L, Huang S, Zhao Y, He X. MicroRNA-30d-5p inhibits tumour cell proliferation and motility by directly targeting CCNE2 in non-small cell lung cancer. Cancer Lett. 2015;362:208–17.
Zhong KZ, Chen WW, Hu XY, Jiang AL, Zhao J. Clinicopathological and prognostic significance of microRNA-107 in human non small cell lung cancer. Int J Clin Exp Pathol. 2014;7:4545–51.
Chen X, Yang H, Zhou X, Zhang L, Lu X. MiR-93 targeting EphA4 promotes neurite outgrowth from spinal cord neurons. J Mol Neurosci. 2016;58:517–24.
Ghasemkhani N, Shadvar S, Masoudi Y, Talaei AJ, Yahaghi E, Goudarzi PK, Shakiba E. Down-regulated MicroRNA 148b expression as predictive biomarker and its prognostic significance associated with clinicopathological features in non-small-cell lung cancer patients. Diagn Pathol. 2015;10:164.
Gilad S, Meiri E, Yogev Y, Benjamin S, Lebanony D, Yerushalmi N, Benjamin H, Kushnir M, Cholakh H, Melamed N, et al. Serum microRNAs are promising novel biomarkers. PLoS One. 2008;3:e3148.
Nitu R, Rogobete AF, Gundogdu F, Tanasescu S, Boruga O, Sas A, Popovici SE, Hutanu D, Pilut C, Sarau CA. MicroRNAs expression as novel genetic biomarker for early prediction and continuous monitoring in pulmonary cancer. Biochem Genet. 2017:1–10.
Grimson A, Farh KK, Johnston WK, Garrett-Engele P, Lim LP, Bartel DP. MicroRNA targeting specificity in mammals: determinants beyond seed pairing. Mol Cell. 2007;27:91–105.
Maragkakis M, Alexiou P, Papadopoulos GL, Reczko M, Dalamagas T, Giannopoulos G, Goumas G, Koukis E, Kourtis K, Simossis VA, et al. Accurate microRNA target prediction correlates with protein repression levels. BMC Bioinformatics. 2009;10:295.
Kolde R, Laur S, Adler P, Vilo J. Robust rank aggregation for gene list integration and meta-analysis. Bioinformatics. 2012;28:573–80.
Essaghir A, Toffalini F, Knoops L, Kallin A, van Helden J, Demoulin JB. Transcription factor regulation can be accurately predicted from the presence of target gene signatures in microarray gene expression data. Nucleic Acids Res. 2010;38:e120.
Essaghir A, Demoulin JB. A minimal connected network of transcription factors regulated in human tumors and its application to the quest for universal cancer biomarkers. PLoS One. 2012;7:e39666.
Wu C, Cao Y, He Z, He J, Hu C, Duan H, Jiang J. Serum levels of miR-19b and miR-146a as prognostic biomarkers for non-small cell lung cancer. Tohoku J Exp Med. 2014;232:85–95.
Donnem T, Fenton CG, Lonvik K, Berg T, Eklo K, Andersen S, Stenvold H, Al-Shibli K, Al-Saad S, Bremnes RM. MicroRNA signatures in tumor tissue related to angiogenesis in non-small cell lung cancer. PLoS One. 2012;7:e29671.
Liu X-G, Zhu W-Y, Huang Y-Y, Ma L-N, Zhou S-Q, Wang Y-K, Zeng F, Zhou J-H, Zhang Y-K. High expression of serum miR-21 and tumor miR-200c associated with poor prognosis in patients with lung cancer. Med Oncol. 2012;29:618–26.
Wang L-K, Hsiao T-H, Hong T-M, Chen H-Y, Kao S-H, Wang W-L, Yu S-L, Lin C-W, Yang P-C. MicroRNA-133a suppresses multiple oncogenic membrane receptors and cell invasion in non-small cell lung carcinoma. PLoS One. 2014;9:e96765.
Wang R, Wang Z, Yang J, Pan X, De W, Chen L. MicroRNA-451 functions as a tumor suppressor in human non-small cell lung cancer by targeting ras-related protein 14 (RAB14). Oncogene. 2011;30:2644–58.
Xu C, Zheng Y, Lian D, Ye S, Yang J, Zeng Z. Analysis of microRNA expression profile identifies novel biomarkers for non-small cell lung cancer. Tumori. 2015;101:104–10.
Lewis BP, Shih I-h, Jones-Rhoades MW, Bartel DP, Burge CB. Prediction of mammalian microRNA targets. Cell. 2003;115:787–98.
Markou A, Tsaroucha EG, Kaklamanis L, Fotinou M, Georgoulias V, Lianidou ES. Prognostic value of mature microRNA-21 and microRNA-205 overexpression in non-small cell lung cancer by quantitative real-time RT-PCR. Clin Chem. 2008;54:1696–704.
Gao W, Shen H, Liu L, Xu J, Xu J, Shu Y. MiR-21 overexpression in human primary squamous cell lung carcinoma is associated with poor patient prognosis. J Cancer Res Clin Oncol. 2011;137:557–66.
Meng F, Henson R, Wehbe-Janek H, Ghoshal K, Jacob ST, Patel T. MicroRNA-21 regulates expression of the PTEN tumor suppressor gene in human hepatocellular cancer. Gastroenterology. 2007;133:647–58.
Zhou L, Yang ZX, Song WJ, Li QJ, Yang F, Wang DS, Zhang N, Dou KF. MicroRNA-21 regulates the migration and invasion of a stem-like population in hepatocellular carcinoma. Int J Oncol. 2013;43:661–9.
Han M, Wang F, Gu Y, Pei X, Guo G, Yu C, Li L, Zhu M, Xiong Y, Wang Y. MicroRNA-21 induces breast cancer cell invasion and migration by suppressing smad7 via EGF and TGF-beta pathways. Oncol Rep. 2016;35:73–80.
Wu XL, Cheng B, Li PY, Huang HJ, Zhao Q, Dan ZL, Tian DA, Zhang P. MicroRNA-143 suppresses gastric cancer cell growth and induces apoptosis by targeting COX-2. World J Gastroenterol. 2013;19:7758–65.
Letelier P, Garcia P, Leal P, Alvarez H, Ili C, Lopez J, Castillo J, Brebi P, Roa JC. miR-1 and miR-145 act as tumor suppressor microRNAs in gallbladder cancer. Int J Clin Exp Pathol. 2014;7:1849–67.
Bao L, Yan Y, Xu C, Ji W, Shen S, Xu G, Zeng Y, Sun B, Qian H, Chen L, et al. MicroRNA-21 suppresses PTEN and hSulf-1 expression and promotes hepatocellular carcinoma progression through AKT/ERK pathways. Cancer Lett. 2013;337:226–36.
Hsiao YC, Kuo WH, Chen PN, Chang HR, Lin TH, Yang WE, Hsieh YS, Chu SC. Flavanone and 2′-OH flavanone inhibit metastasis of lung cancer cells via down-regulation of proteinases activities and MAPK pathway. Chem Biol Interact. 2007;167:193–206.
Lin CW, Chang YL, Chang YC, Lin JC, Chen CC, Pan SH, Wu CT, Chen HY, Yang SC, Hong TM, Yang PC. MicroRNA-135b promotes lung cancer metastasis by regulating multiple targets in the Hippo pathway and LZTS1. Nat Commun. 2013;4:1877.
Young RA. RNA polymerase II. Annu Rev Biochem. 1991;60:689–715.
This study was supported by Technologies R&D Program of Shaanxi Province, the Diagnosis and Micrometastases Research in Early Stage Non-small Cell Lung Cancer (2012K13-02-16).
Availability of data and materials
All data and material were available in GEO database.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.