- Open Access
Down-regulation of microRNA-144-3p and its clinical value in non-small cell lung cancer: a comprehensive analysis based on microarray, miRNA-sequencing, and quantitative real-time PCR data
Respiratory Researchvolume 20, Article number: 48 (2019)
Previous studies have shown that miR-144-3p might be a potential biomarker in non-small cell lung cancer (NSCLC). Nevertheless, the comprehensive mechanism behind the effects of miR-144-3p on the origin, differentiation, and apoptosis of NSCLC, as well as the relationship between miR-144-3p and clinical parameters, has been rarely reported.
We investigated the correlations between miR-144-3p expression and clinical characteristics through data collected from Gene Expression Omnibus (GEO) microarrays, the relevant literature, The Cancer Genome Atlas (TCGA), and real-time quantitative real-time PCR (RT-qPCR) analyses to determine the clinical role of miR-144-3p in NSCLC. Furthermore, we investigated the biological function of miR-144-3p by Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. Protein-protein interaction (PPI) network was created to identify the hub genes.
From the comprehensive meta-analysis, the combined SMD of miR-144-3p was − 0.95 with 95% CI of (− 1.37, − 0.52), indicating that less miR-144-3p was expressed in the NSCLC tissue than in the normal tissue. MiR-144-3p expression was significantly correlated with stage, lymph node metastasis and vascular invasion (all P < 0.05). As for the bioinformatics analyses, a total of 37 genes were chosen as the potential targets of miR-144-3p in NSCLC. These promising target genes were highly enriched in various key pathways such as the protein digestion and absorption and the thyroid hormone signaling pathways. Additionally, PPI revealed five genes—C12orf5, CEP55, E2F8, STIL, and TOP2A—as hub genes with the threshold value of 6.
The current study validated that miR-144-3p was lowly expressed in NSCLC. More importantly, miR-144-3p might function as a latent tumor biomarker in the prognosis prediction for NSCLC. The results of bioinformatics analyses may present a new method for investigating the pathogenesis of NSCLC.
Lung cancer (LC) is recognized as a life-threatening malady as the incidence and mortality rates are ranked second among all neoplasms worldwide . Accounting for 85% of all diagnosed LC cases, non-small cell lung cancer (NSCLC) is divided mainly into lung squamous cell carcinoma (LUSC), lung adenocarcinoma (LUAD), and large cell carcinoma (LCC) . Currently, the main therapy for NSCLC is a combination of surgery and chemotherapy . Although great progress has been made in the early detection, diagnosis, and targeted treatments of NSCLC, the five-year survival rates are still low, varying from 4 to 17% depending on regional differences and the disease stage [3,4,5]. Patients with severer tumor and comorbidity burdens are at a higher risk of death from not receiving effective or specific treatments . Thus, the understanding of the molecular mechanisms in NSCLC and the identification of new therapeutic targets are crucial.
MiRNAs are single-stranded ncRNAs of approximately 20 nucleotides in length. They play an important role in the regulation of gene expression by post-transcriptionally binding with the mRNAs of target genes . Participating in cellular differentiation and homeostasis, miRNAs play pivotal roles in cancer . Tissue-specific miRNAs act as novel potential biomarkers in the diagnosis, treatment, and prognosis of cancer [9, 10]. Recently, the effects of miRNAs on oncogenesis and tumor progression have been receiving a great deal of attention.
MicroRNA-144-3p (miR-144-3p), having various functions in different types of cancers, is one of the miRNAs related to cancer. MiR-144-3p acts as a suppressive factor in laryngeal squamous cell carcinoma , gastric cancer , hepatocellular carcinoma , and pancreatic cancer . However, it is an oncogene in renal carcinoma , nasopharyngeal carcinoma , and colorectal cancer . Previous studies have also implicated that miR-144-3p is involved in cell proliferation, apoptosis, and autophagy by targeting the TP53-inducible glycolysis and apoptosis regulator (TIGAR) in LC . The down-regulation of miR-144-3p results in metabolic alterations of LC cells by regulating the glucose transporter 1(GLUT1) . Previous studies have shown that miR-144-3p might be a biomarker and target with great potential. Nevertheless, the comprehensive mechanism behind the effects of miR-144-3p on the origin, differentiation, and apoptosis of NSCLC, as well as the relationship between miR-144-3p and clinical parameters, has been rarely reported.
This study investigated the correlations between miR-144-3p expression and clinical characteristics through data collected from Gene Expression Omnibus (GEO) microarrays, the relevant literature, The Cancer Genome Atlas (TCGA), and real-time quantitative real-time PCR (RT-qPCR) analyses to determine the clinical role of miR-144-3p in NSCLC. The latent mechanism of action in NSCLC was subsequently examined by using 12 predictive programs to forecast the genes targeted by miR-144-3p. In addition, bioinformatics analyses, which included Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and protein–protein interaction (PPI) network analyses, were performed.
Materials and methods
A microarray search of miR-144-3p in NSCLC was conducted in the GEO database (http://www.ncbi.nlm.nih.gov/geo/) with the following keywords: (MicroRNA OR “Micro RNA” OR “non-coding RNA” OR ncRNA OR “small RNA” OR miRNA) AND (Lung OR pulmonary) AND (cancer OR tumor OR neoplasm OR malignancy OR carcinoma OR adenocarcinoma OR AC OR SCC OR NSCLC). The entry type was restricted to “series,” and the organism was filtered by “Homo sapiens.” The criteria for inclusion were as follows: (1) patients diagnosed with NSCLC and its subtypes were investigated; (2) cancerous and noncancerous samples were involved; (3) the healthy and malignant groups included at least three samples in the form of tissue, blood, or plasma; (4) and the expression profiling data for miR-144-3p were available. Related studies were retrieved from the PubMed, Google Scholar, China National Knowledge Infrastructure (CNKI), Chongqing VIP electronic (VIP) and Chinese Wanfang databases. Figure 1 shows the workflow for the study.
MicroRNA-144-3p expression data from the Cancer genome atlas
TCGA (https://cancergenome.nih.gov/) was used for obtaining detailed information about the expression value of miR-144-3p in NSCLC and noncancerous samples. The differences in the miR-144-3p expression in the NSCLC samples and the normal controls were calculated with IBM SPSS Statistics V22.0 software.
Quantitative real-time PCR
For the current study, 125 matched samples were provided by the Department of Pathology of the First Affiliated Hospital of Guangxi Medical University. Formalin fixation and paraffin embedding (FFPE) were performed to preserve the specimens. This aspect of the research was authorized by the hospital’s ethics committee. Next, the miR-144-3p expression in 125 paired clinical samples was detected by RT-qPCR with the Applied Biosystems 7900HT Fast Real-Time PCR System software. The miR-144-3p sequence was as follows: UACAGUAUAGAUGAUGUACU. The formula for 2-Δcq was used for the calculation of the miR-144-3p expression value.
Statistical analysis and comprehensive meta-analysis
After miR-144-3p was log2-transformed, the expression profiling information was used to calculate the number, mean (M) and standard deviation (SD) for each control and experimental group with IBM SPSS Statistics V22.0 software. In addition, Stata 12.0 software was used for performing a comprehensive meta-analysis of data aggregated from multiple sources (microarray, literature, miRNA sequencing, and RT-qPCR). The analysis of the miR-144-3p expression in the NSCLC and tumor-free specimens was displayed on forest plots that illustrate the standardized mean difference (SMD) and the 95% confidential interval (CI). The chi-squared test of Q and the I2 statistic were calculated to assess heterogeneity across the studies and to determine the appropriateness of applying either a random effects model or fixed effects model to the pooling process. To measure publication bias, Egger’s and Begg’s tests and a funnel plot, for which significance was p < 0.05, were performed.
Latent targets of microRNA-144-3p in non-small cell lung cancer
MiRWALK2.0, an online archive of data on miRNA-target interactions , was mined to forecast the genes targeted by miR-144-3p. In total, 12 servers with miRWalk, miRMap, MicroT4, miRNAMap, TargetScan PICTAR2, miRBridge, PITA, miRanda, RNAhybrid, miRDB, RNA22 were used. Only those genes projected by more than six of the servers were recognized as target genes. The high-expressed genes in LUAD and LUSC were acquired through Gene Expression Profiling Interactive Analysis (GEPIA). The overlapping genes among the up-regulated genes in LUAD and LUSC and the predicted target genes, were viewed as promising targets of miR-144-3p in NSCLC. A review of the literature on the specific target genes of miR-144-3p in NSCLC was conducted. The target genes determined by previous studies and the predicted target genes, namely promising target genes, were used in the functional analysis.
Functional analysis for promising target genes
The GO vocabularies, which include biological processes (BPs), cellular components (CCs), and molecular functions (MFs), were enriched by Metascape (http://metascape.org/gp/). The functional annotation of the underlying target genes was then elucidated by KEGG pathway analysis with Metascape tool. In addition, a PPI network was constructed to reveal the hub genes of the potential target genes on STRING, a web portal for undermining the integrated function of multiple genes .
Expression of hub genes from the Cancer genome atlas and the genotype-tissue expression database
To further confirm the function of hub genes in NSCLC and their relationship to miR-144-3p, a search of TCGA and the Genotype-Tissue Expression (GTEx) database was performed to determine the expression pattern of the hub genes in NSCLC. Box plots of the hub genes in the NSCLC and non-cancer samples were developed through GEPIA.
Confirmation of the expression and clinical value of microRNA-144-3p in non-small cell lung cancer, based on gene expression omnibus datasets
MicroRNA-144-3p expression in non-small cell lung cancer obtained through gene expression omnibus microarrays
A total of 19 microarrays from the GEO database met the entry criteria. The features of the included GEO datasets are depicted in Table 1. Of the microarrays, 14 were obtained from tissue, and 5 were derived from blood (GSE27486, GSE40738, GSE64951, GSE93300, and GSE114711). In addition, the expression data from the NSCLC and control groups were collected on the basis of the GEO database. With respect to the data from the tissue samples, the NSCLC groups had a significantly lower level of miR-144-3p expression than the control groups in GSE25508, GSE48414, GSE51853, GSE56036, GSE63805, GSE72526, GSE74190, and GSE102286 (p = 0.0202, p < 0.0001, p < 0.0001, p = 0.0011, p < 0.0001, p = 0.0102, p < 0.0001, and p < 0.0001, respectively (Fig. 2)). In contrast, no notable distinction in miR-144-3p expression was detected between the NSCLC and the control groups in the other microarrays (GSE14936, GSE29248, GSE36681, GSE47525, GSE53882, and GSE77380). Regarding the data from the blood samples, miR-144-3p expression in NSCLC was found to decrease significantly in GSE27486 and GSE40738 (p = 0.0196, p = 0.0036, respectively (Fig. 3)).
Results of meta-analysis of gene expression omnibus datasets
A meta-analysis was conducted on the basis of the 19 included microarrays from the GEO database. The results are demonstrated in Fig. 4a. Given the apparent heterogeneity (p < 0.05, I2 = 94.3%), a random effects model was applied, and remarkable down-regulation (SMD = − 0.89; 95% CI − 1.34, − 0.44; p = 0.000) of miR-144-3p was found in the NSCLC groups. A sensitivity analysis was later conducted to explore whether a particular microarray played a vital role in significant heterogeneity (Fig. 4b). After an individual study was removed each time of meta-analysis, the combined effect was compared to the previous one. No study was found to have played a crucial role in any of the enrolled studies.
A funnel plot was generated to estimate publication bias (Fig. 4c). To further clarify the heterogeneity source, a subgroup analysis was performed. It was based on multiple characteristics: sample source (tissue vs. blood), and cancer type (adenocarcinoma vs. squamous cell carcinoma). As is illustrated in Fig. 5, significant heterogeneity was observed in the tissue subgroup (I2 = 95.8%, p = 0.000). Significant heterogeneity was also found in adenocarcinoma (I2 = 92.9%, p = 0.000) and squamous cell carcinoma (I2 = 95.3%, p = 0.000). These results suggest that sample source and cancer type might be sources of heterogeneity.
A review of the literature relevant to miR-144-3p was conducted through searches in the PubMed, Google Scholar, CNKI, VIP, and Wanfang databases. Neither the M nor the SD of miR-144-3p in the NSCLC and normal groups was provided in the literature; therefore, no available data could be obtained from the existing studies.
Confirmation of the expression and clinical effects of microRNA-144-3p in non-small cell lung cancer, based on the Cancer genome atlas data
MicroRNA-144-3p expression and prognostic value in non-small cell lung cancer tissues
TCGA contained 376 samples for LUSC patients and 488 samples for LUAD patients. Regarding LUSC, miR-144-3p expression was remarkably downregulated in comparison with the normal controls (2.8193 ± 1.40600 vs. 5.5678 ± 1.27693, p < 0.001 (Fig. 6a and Table 2)). In terms of LUAD, the expression level of miR-144-3p was obviously less than that in the healthy tissue (2.8959 ± 1.35967 vs. 5.2775 ± 1.64708, p < 0.001 (Fig. 6b and Table 3)). The data from LUSC and LUAD were combined for further examination of the miR-144-3p expression in NSCLC. As is illustrated in Fig. 6c and Table 4, miR-144-3p was significantly reduced in the NSCLC tissue compared to the non-cancerous lung tissue (2.8632 ± 1.37928 vs. 5.4243 ± 1.4702, p < 0.0001). A Kaplan–Meier curve was later used to identify the effects of the expression of miR-144-3p on survival time. As is shown in Fig. 7, the p values for the three Kaplan–Meier curves were all greater than 0.05, thus indicating no significant difference in survival time between the group with low levels of miR-144-3p and the one with high levels.
Relationships between microRNA-144-3p and clinical pathology of non-small cell lung cancer, based on the Cancer genome atlas data
As can be seen in Tables 2 and 3, the clinical characteristics of 332 LUSC patients and 445 LUAD patients were downloaded from TCGA. Regarding LUSC, a significant difference in miR-144-3p was found for stage (p = 0.040) and primary tumor (T) (p = 0.035). LUSC patients in Stages III–IV (2.4469 ± 1.41079) had a lower expression of miR-144-3p than those in Stages I–II (2.8878 ± 1.39556). The miR-144-3p expression of LUSC in T3–T4 (2.5088 ± 1.36178) was more significantly decreased than T1–T2 (2.9023 ± 1.40854). In terms of LUAD, a significant difference of miR-144-3p expression was observed for stage (p = 0.031). Patients in Stages III–IV (3.1868 ± 1.45395) had higher expression values of miR-144-3p. The data on LUAD and LUSC, based on TCGA, were pooled for further validation. As is illustrated in Table 4, the significance in the statistics for the T stage was based on the lower miR-144-3p expression in patients in T3–T4 (p < 0.05).
Quantitative real-time PCR analysis
The microRNA-144-3p expression and its significance in non-small cell lung cancer prognoses
Using RT-qPCR, the clinical expression value of miR-144-3p in 125 matched tissues was evaluated. As is illustrated in Fig. 8a and Table 5, the NSCLC samples exhibited a significantly lower expression level of miR-144-3p than the non-cancerous samples (2.808 ± 1.303 vs. 4.813 ± 2.618, p < 0.001). Next, miR-144-3p expression was then analyzed for LUAD and for LUSC. As is shown in Fig. 8b and c, unlike what was found in the adjacent non-cancerous tissue, apparently lowly expressed miR-144-3p was observed in both LUSC and LUAD (p = 0.0004, p < 0.0001). A Kaplan–Meier curve was generated to assess whether miR-144-3p is appropriate for the prognosis prediction of LUAD. As is depicted in Fig. 9, the LUAD patients who exhibited lower expression values of miR-144-3p might have worse outcomes (p = 0.397).
Correlations between microRNA-144-3p expression and clinical characteristic for non-small cell lung cancer patients
The miR-144-3p expression in NSCLC cases was significantly different in lymph node metastasis and vascular invasion (Table 5). The lower value of miR-144-3p was found in patients with lymph node metastasis but not in those without it (Table 5). Patients with vascular invasion maintained a miR-144-3p value of 2.1400 ± 1.2263, and the expression level of those without vascular invasion was 3.0682 ± 1.24395 (Table 5). To further validate the correlation of miR-144-3p and clinicopathological characteristics, the NSCLC cases were divided into LUSC and LUAD groups. For the LUSC group (Table 6), no statistical differences were seen for smoking, vascular invasion, or lymph node metastasis. However, for LUAD, the statistical analyses indicated significant differences for smoking, vascular invasion, and lymph node metastasis. In comparison to patients without vascular invasion, those with vascular invasion had lower miR-144-3p expression values. The miR-144-3p expression in patients with a smoking habit was significantly down-regulated over that of patients without the habit (p = 0.027, Table 7). Besides, the miR-144-3p expression in NSCLC patients who were considered to have lymph node metastasis was markedly reduced.
Meta-analysis of combination of gene expression omnibus, the Cancer genome atlas, and quantitative real-time PCR data
Because no relevant studies were found, data from three sources (microarrays, miRNA-sequencing, and RT-qPCR) containing 2264 NSCLC samples and 968 non-cancerous samples were extracted for integrated meta-analyses with a random effects model because of high heterogeneity (I2 = 95.4%, p = 0.000). The differences in individuals, NSCLC subtype, and sample source were considered the sources of heterogeneity. As is shown in Fig. 10a, the combined SMD of miR-144-3p was − 0.95 with 95% CI of (− 1.37, − 0.52), indicating that less miR-144-3p was expressed in the NSCLC tissue than in the normal tissue. The sensitivity analysis (Fig. 10b) indicated significant differences among the studies; however, no specific study had a significant effect on high heterogeneity. The evaluation of publication bias was performed with Begg’s and Egger’s tests and a funnel plot (Fig. 10c). In general, the funnel plot was symmetrical, and the p values obtained from the Begg’s and Egger’s tests were 0.833 and 0.335, respectively. In sum, the results indicate that publication bias for the studies was controlled passably.
Promising target genes collection
From miRWALK2.0, 1635 genes targeted by miR-144-3p in NSCLC predicted by more than six algorithms were obtained. A total of 1109 overexpressed genes in LUAD were collected on the basis of GEPIA. In addition, 1922 overexpressed genes in LUSC were acquired from GEPIA. After intersection, 34 predicted target genes were selected. Four specific targets of miR-144-3p were verified in previous studies related to LC (Table 8). The TIGAR is also known as the C12orf5 gene. Accordingly, a total of 37 potential target genes were collected.
Gene ontology and Kyoto encyclopedia of genes and genomes analyses
For further interpretation of the function of the promising genes targeted by miR-144-3p in NSCLC, KEGG, and GO annotations were performed in Metascape. For the GO analysis, three categories were used: BP, CC, and MF. For the BP, renal system development (GO: 0072001) and the nucleobase-containing small molecule metabolic process (GO: 0055086) were the top two pathways (Fig. 11a). For the CC, the potential target differentially expressed genes (DEGs) were predominantly enriched in the centriole (GO: 0005814), Golgi membrane (GO: 0000139), and mitochondrial envelope (GO: 0005740) (Fig. 11c). For MF, the three significantly involved items were RNA polymerase II proximal promoter sequence-specific DNA binding (GO: 0000978); cofactor binding (GO: 0048037); and transferase activity, transferring glycosyl groups (GO: 0016757) (Fig. 11b). Regarding KEGG, the top two enriched pathways were the protein digestion and absorption (hsa04974) and the thyroid hormone signaling pathways (hsa04919) (Fig. 12). PPI revealed five genes—C12orf5, CEP55, E2F8, STIL, and TOP2A—as hub genes with the threshold value of 6 (Fig. 13).
Expression of hub genes from the Cancer genome atlas and the genotype-tissue expression datasets
Of the five hub genes, not including C12orf5, that occupied the central region of the PPI network, four genes (CEP55, E2F8, STIL, and TOP2A) were significantly up-regulated in the NSCLC group compared to the control group (Fig. 14).
Although previous studies have documented the expression of miR-144-3p in NSCLC, correlations of the clinical features with NSCLC and miR-144-3p have been seldom reported. This study used a larger number of samples for a systematic investigation of the relationships. A highlight of this study was the use of the computational biology method to explore the latent mechanism of miR-144-3p in NSCLC.
A decreased expression of miR-144-3p was found in the NSCLC cases, as presented by the GEO, TCGA, and RT-qPCR. A comprehensive meta-analysis was the focus of the current study. Data were obtained from multiple sources: namely, microarrays, previous studies, RT-qPCR, and TCGA. The meta-analysis results revealed that the decline of miR-144-3p in NSCLC was consistent across studies [19, 22]. It was therefore concluded that there was a considerable decrease of miR-144-3p expression in NSCLC. Based on the data from TCGA, the miR-144-3p expression level in LUSC and LUAD was related to stage, and the predominant difference for T was found in LUAD only.
The results revealed that miR-144-3p could be involved in the occurrence and development of NSCLC, and low miR-144-3p could indicate the promotion of NSCLC. Regarding RT-qPCR, miR-144-3p expression was related to lymph node metastasis and vascular invasion. Moreover, the patients with lymph node metastasis and vascular invasion had a tendency to down-regulate miR-144-3p expression. Regarding LUAD, the amount of miR-144-3p differed greatly depending on the smoking status of the patient, the existence of lymph node metastasis, and the presence of vascular invasion. However, no statistical differences were found for LUSC. In sum, miR-144-3p might serve as a marker to monitor the progression of NSCLC.
No differences were observed in the survival times for the low-miR-144-3p group and the high-miR-144-3p group, according to TCGA data. In contrast, the PCR results revealed a trend that suggests that LUAD patients with lower miR-144-3p levels might have worse outcomes although not to a significant extent. A study by Wu et al.  revealed that miR-144-3p was one of the independent prognostic risk factors for LUAD patients, with those with low miR-144-3p having poor prognoses. Further investigations are required for corroboration.
At the present time, the specific molecular mechanism of NSCLC is not widely understood. Therefore, bioinformatics analyses were performed to discover the inherent mode of NSCLC activity at the molecular level. Based on miRWALK2.0 and TCGA data, the candidate targets of miR-144-3p were projected. In an attempt to explore the roles of these genes further, KEGG and GO annotation analyses were performed. According to GO enrichment, the candidate targets of miR-144-3p might have an important effect on the progression of NSCLC by modulating several cellular biology processes, such as renal system development. Moreover, these genes could also have an important effect, such as cofactor binding, on MF. The results of the KEGG analysis also showed the roles of the candidate targets of miR-144-3p in NSCLC. The top two enriched pathways were the protein digestion and absorption and the thyroid hormone signaling pathways. This suggests that the promising targets of miR-144-3p could be involved in the aforementioned pathways to influence the occurrence and progression of NSCLC.
The protein digestion and absorption pathway is a key pathway in several human cancers. A study by Shi et al. showed that this pathway might contribute to the pulmonary metastasis of osteosarcoma patients . It might be involved in the up-regulation of differentially expressed genes in breast cancer , and it has already been associated with the down-regulated differentially expressed genes in pancreatic neuroendocrine tumors . The protein digestion and absorption pathway is connected mainly to differentially expressed genes, and this affects the occurrence and development of enchondromas . Additionally, B-cell malignancies are relevant to the protein digestion and absorption pathway . However, few studies have been conducted on the protein digestion and absorption pathway in NSCLC. Therefore, more studies are required to further validate the specific molecular mechanism in NSCLC.
According to the PPI network data, five genes (TIGAR, CEP55, E2F8, STIL, and TOP2A) have been recognized as hub genes in NSCLC, thus proving their roles as ideal candidates for miR-144-3p targets. The current study found that four of the 5 hub genes—CEP55, E2F8, STIL, and TOP2A—were more greatly upregulated in the NSCLC cases than in the normal cases. In addition, the decreased miR-144-3p expression in NSCLC was validated in this study. Therefore, the overexpression of the four hub genes in NSCLC was indirect proof that these genes might be the targets of miR-144-3p. By targeting the TIGAR, miR-144-3p suppresses proliferation, mediates programmed cell death, and increases autophagy in LC cells . Studies have found that CEP55 might have an important effect on the proliferation of LC [29, 30]. The role of the E2F8 gene has been studied in several cancers. Sun et al.  provided evidence that E2F8 was targeted by miR-144-3p in papillary thyroid cancer; however, a connection between E2F8 and miR-144-3P in NSCLC has not been reported so far.
The current study has limitations. Because a large-scale clinical sample was not used, more investigations with large clinical samples are required for further confirmation of the function of miR-144-3p in NSCLC prognoses. In addition, experiments were not conducted for the detection of the expression of hub genes. The expression of these genes needs to be verified by additional well-designed studies. Although bioinformatics analyses were performed, the specific molecular mechanisms were not identified.
In conclusion, the current study validated that miR-144-3p was lowly expressed in NSCLC. More importantly, miR-144-3p might function as a latent tumor biomarker in the prognosis prediction for NSCLC. The results of bioinformatics analyses may present a new method for investigating the pathogenesis of NSCLC.
Gene Expression Omnibus
Kyoto Encyclopedia of Genes and Genomes
Lung large cell carcinoma
Lung squamous cell carcinoma
Non-small cell lung cancer
Protein–protein interaction network
Quantitative real-time PCR
Standardized mean difference
The Cancer Genome Atlas
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2018. CA Cancer J Clin. 2018;68:7–30.
Inamura K. Lung cancer: understanding its molecular pathology and the 2015 who classification. Front Oncol. 2017;7:193.
Hirsch FR, Scagliotti GV, Mulshine JL, Kwon R, Curran WJ Jr, Wu YL, et al. Lung cancer: current therapies and new targeted treatments. Lancet. 2017;389:299–311.
Yun YH, Kim YA, Sim JA, Shin AS, Chang YJ, Lee J, et al. Prognostic value of quality of life score in disease-free survivors of surgically-treated lung cancer. BMC Cancer. 2016;16:505.
Tang Y, Qiao G, Xu E, Xuan Y, Liao M, Yin G. Biomarkers for early diagnosis, prognosis, prediction, and recurrence monitoring of non-small cell lung cancer. Onco Targets Ther. 2017;10:4527–34.
Rios J, Gosain R, Goulart BH, Huang B, Oechsli MN, McDowell JK, et al. Treatment and outcomes of non-small-cell lung cancer patients with high comorbidity. Cancer Manag Res. 2018;10:167–75.
Beermann J, Piccoli MT, Viereck J, Thum T. Non-coding RNAs in development and disease: background, mechanisms, and therapeutic approaches. Physiol Revs. 2016;96:1297.
Bracken CP, Scott HS, Goodall GJ. A network-biology perspective of microRNA function and dysfunction in cancer. Nat Rev Genet. 2016;17:719.
Zhou Q, Huang SX, Zhang F, Li SJ, Liu C, Xi YY, et al. MicroRNAs: a novel potential biomarker for diagnosis and therapy in patients with non-small cell lung cancer. Cell Prolif. 2017. https://doi.org/10.1111/cpr.12394.
Kang L, Tao X, He BS, Pan YQ, Sun HL, Peng HX, et al. MicroRNA expression profiles predict progression and clinical outcome in lung adenocarcinoma. Onco Targets Ther. 2016;9:5679–92.
Zhang SY, Lu ZM, Lin YF, Chen LS, Luo XN, Song XH, et al. miR-144-3p, a tumor suppressive microRNA targeting ETS-1 in laryngeal squamous cell carcinoma. Oncotarget. 2016;7:11637.
Li B, Zhang S, Hao S, Li C. MicroRNA-144-3p suppresses gastric cancer progression by inhibiting epithelial-to-mesenchymal transition through targeting PBX3. Biochem Biophys Res Commun. 2017;484:241–7.
Liang H, Ye Z, Yin S, Mo W, Wang H, Zhao J, et al. A comprehensive insight into the clinicopathologic significance of miR-144-3p in hepatocellular carcinoma. Onco Targets Ther. 2017;10:3405–19.
Li J, Sun P, Yue Z, Zhang D, You K, Wang J. miR-144-3p induces cell cycle arrest and apoptosis in pancreatic cancer cells by targeting proline-rich protein 11 expression via the mitogen-activated protein kinase signaling pathway. DNA Cell Biol. 2017;36:619.
Xiao W, Lou N, Ruan H, Bao L, Xiong Z, Yuan C, et al. Mir-144-3p promotes cell proliferation, metastasis, sunitinib resistance in clear cell renal cell carcinoma by downregulating ARID1A. Cell Physiol Biochem. 2017;43:2420–33.
Zhang LY, Hofun LV, Wong AM, Kwong DL, Zhu YH, Dong SS, et al. MicroRNA-144 promotes cell proliferation, migration and invasion in nasopharyngeal carcinoma through repression of PTEN. Carcinogenesis. 2013;34:454.
Ye JJ, Cao J. MicroRNAs in colorectal cancer as markers and targets: recent advances. World J Gastroenterol. 2014;20:4288.
Chen S, Li P, Li J, Wang Y, Du Y, Chen X, et al. MiR-144 inhibits proliferation and induces apoptosis and autophagy in lung cancer cells by targeting TIGAR. Cell Physiol Biochem. 2015;35:997–1007.
Liu M, Gao J, Huang Q, Jin Y, Wei Z. Downregulating microRNA-144 mediates a metabolic shift in lung cancer cells by regulating GLUT1 expression. Oncol Lett. 2016;11:3772.
Dweep H, Gretz N. miRWalk2.0: a comprehensive atlas of microRNA-target interactions. Nat Methods. 2015;12:697.
Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, et al. STRING v10: protein–protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 2015;43:D447–52.
Zhang G, An H, Fang X. MicroRNA-144 regulates proliferation, invasion, and apoptosis of cells in malignant solitary pulmonary nodule via zinc finger E-box-binding homeobox 1. Int J Clin Exp Pathol. 2015;8:5960–7.
Wu C, Xu B, Zhou Y, Ji M, Zhang D, Jiang J, Wu C. Correlation between serum IL-1β and miR-144-3p as well as their prognostic values in LUAD and LUSC patients. Oncotarget. 2016;7:85876–87.
Shi Z, Zhou H, Pan B, Lu L, Wei Z, Shi L, et al. Exploring the key genes and pathways of osteosarcoma with pulmonary metastasis using a gene expression microarray. Mol Med Rep. 2017;16:7423–31.
Dong LF, Xu SY, Long JP, Wan F, Chen YD. RNA-sequence analysis reveals differentially expressed genes (DEGs) in patients exhibiting different risks of tumor metastasis. Med Sci Monit. 2017;23:2842–9.
Wang DD, Liu ZW, Han MM, Zhu ZM, Tu YL, Dou CQ, et al. Microarray based analysis of gene expression patterns in pancreatic neuroendocrine tumors. Eur Rev Med Pharmacol Sci. 2015;19:3367–74.
Shi Z, Zhou H, Pan B, Lu L, Kang Y, Liu L, et al. Exploring the key genes and pathways in enchondromas using a gene expression microarray. Oncotarget. 2017;8:43967–77.
Hsu LI, Briggs F, Shao X, Metayer C, Wiemels JL, Chokkalingam AP, et al. Pathway analysis of genome-wide association study in childhood leukemia among Hispanics. Cancer Epidemiol Biomark Prev. 2016;25:815–22.
Ma XP, Zhang W, Wu BQ, Qin J. Correlations between mRNA levels of Centrosomal protein 55 (CEP55) and clinical features of patients with lung cancer. Med Sci Monit. 2018;24:3093–7.
Liu L, Mei Q, Zhao J, Dai Y, Fu Q. Suppression of CEP55 reduces cell viability and inducesapoptosis in human lung cancer. Oncol Rep. 2016;36:1939–45.
Sun J, Shi R, Zhao S, Li X, Lu S, Bu H, et al. E2F8, a direct target of miR-144, promotes papillary thyroid cancer progression via regulating cell cycle. J Exp Clin Cancer Res. 2017;36:40.
The current study was supported by the future academic star project of Guangxi Medical University (2017).
Availability of data and materials
One hundred twenty-five (125) paired non-small cell lung cancer specimens were obtained from the First Affiliated Hospital of Guangxi Medical University (Guangxi, China).
Ethics approval and consent to participate
This study was approved by the Ethics Committee of the First Affiliated Hospital of Guangxi Medical University.
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.