Figures
Abstract
Tumor-promoting cytokines are a cause of tumor progression; therefore, identifying key regulatory microRNAs (miRNAs) for controlling their production is important. The aim of this study is to identify promising miRNAs associated with tumor-promoting cytokines in non-small cell lung cancer (NSCLC). We identified circulating miRNAs from 16 published miRNA profiles. The selected miRNAs were validated in the serum of 32 NSCLC patients and compared with 33 patients with other lung diseases and 23 healthy persons using quantitative real-time PCR. The cytokine concentration was investigated using the enzyme-linked immunoassay in the same sample set, with clinical validation of the miRNAs. The correlation between miRNA expression and cytokine concentration was evaluated by Spearman’s rank correlation. For consistent direction, one up-regulated miRNA (miR-145) was found in four studies, and seven miRNAs were reported in three studies. One miRNA (miR-20a) and four miRNAs (miR-25-3p, miR-223, let-7f, and miR-20b) were reported in six and five studies. However, their expression was inconsistent. In the clinical validation, serum miR-145 was significantly down-regulated, whereas serum miR-20a was significantly up-regulated in NSCLC, compared with controls. Regarding serum cytokine, all cytokines [vascular endothelial growth factor (VEGF), interleukin-6 (IL-6), and transforming growth factor β (TGF-β)], except tumor necrosis factor-α (TNF-α), had a higher level in NSCLC patients than controls. In addition, we found a moderate correlation between the TGF-β concentration and miR-20a (r = −0.537, p = 0.002) and miR-223 (r = 0.428, p = 0.015) and a weak correlation between the VEGF concentration with miR-20a (r = 0.376, p = 0.037) and miR-223 (r = −0.355, p = 0.046). MiR-145 and miR-20a are potential biomarkers for NSCLC. In addition, the regulation of tumor-promoting cytokine, through miR-20a and miR-223, might be a new therapeutic approach for lung cancer.
Citation: Chaniad P, Trakunran K, Geater SL, Keeratichananont W, Thongsuksai P, Raungrut P (2020) Serum miRNAs associated with tumor-promoting cytokines in non-small cell lung cancer. PLoS ONE 15(10): e0241593. https://doi.org/10.1371/journal.pone.0241593
Editor: Esra Bozgeyik, Tekirdag Namik Kemal University, TURKEY
Received: March 12, 2020; Accepted: October 16, 2020; Published: October 30, 2020
Copyright: © 2020 Chaniad et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the manuscript and its Supporting information files.
Funding: This study is funded by a grant from the Prince of Songkla University [grant number MED590692S], the Research Center for Cancer Control in Thailand [grant number MEDRC59036], and the Faculty of Medicine, Prince of Songkla University [grant number REC60350042]. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Lung cancer is the most common and leading cause of cancer death worldwide, accounting for 11.6% of the total cases and 18.4% of the total cancer deaths in 2018, respectively [1]. Approximately 84% of lung cancer is non-small cell lung cancer (NSCLC), which is usually diagnosed in the advanced stage (30%–79% of all NSCLC cases) and accompanied by lymph node and distant metastasis [2, 3]. Tumor progression is supported by tumor-promoting cytokines, in which high serum levels of several cytokines, such as vascular endothelial growth factor (VEGF), interleukin-6 (IL-6), transforming growth factor β (TGF-β), and tumor necrosis factor-α (TNF-α), are correlated with advanced stages of lung cancer [4, 5]. Currently, there is accumulating evidence showing that these cytokines are regulated by microRNAs (miRNAs). Therefore, the identification of miRNAs associated with cytokine production is desired.
MiRNAs are short noncoding RNAs, 20–22 nucleotides in length. They suppress gene expression at the post-transcriptional level by binding to the 3’-untranslated regions of mRNAs [6]. Numerous studies have shown that miRNAs play a part in cell differentiation, proliferation, apoptosis, and cytokine production [6, 7]. In addition, the remarkable stability of serum miRNAs under various conditions was found in our previous report [8]. Several studies have revealed that the expression levels display oscillatory changes in response to TNF-α in macrophage cells [9], and VEGF and IL-6 in bone marrow mesenchymal stromal cells [10]. Moreover, changes in miRNA expression have been demonstrated to be involved with cytokine signaling pathways in several types of cancers, including colorectal [11], breast [12], and lung [13] cancers. Although there has been accumulating evidence to support the regulation of cytokines by miRNAs, studies conducted to determine the association of miRNAs and tumor-promoting cytokines in clinical samples in lung cancer are limited.
In this study, we identified differentially expressed (DE) miRNAs through the systematic review of previously published studies on miRNA profiling in NSCLC. A vote-counting procedure and bioinformatics analysis were used to obtain miRNAs associated with cytokine signaling pathways. We then validated the selected miRNAs using quantitative real-time PCR in serum samples of NSCLC compared with controls, which included patients with other lung diseases (OL) and healthy persons (HP). We also determined the serum concentration of the common tumor-promoting cytokines (IL-6, VEGF, TGF-β, and TNF-α) using the enzyme-linked immunoassay (ELISA) in samples from the same patients with miRNA validation. Finally, we investigated the association of miRNA expression and cytokine concentration in the serum of NSCLC patients.
Materials and methods
Literature search
The literature on miRNA profiling of lung cancer, published between 2006 and 2016, was retrieved from “PubMed” and “Scopus” databases. The search terms “miRNA or miR or microRNA,” “lung cancer,” and “profiling” were used.
Study selection criteria
Eligible studies were required to meet the following inclusion criteria: 1) miRNA expression profiling of patients with NSCLC; 2) conducted on blood, serum, or plasma; 3) used samples from patients with OL or HP for comparison; 4) were full-text articles in English. The exclusion criteria were: 1) patients received any previous treatment; 2) studies that did not report p-values, or reported a false discovery rate (FDR); 3) studies were review articles; 4) studies that did not report sample size.
Data extraction
Lists of DE miRNAs were extracted from each selected article. Related information was also retrieved, including authors, year of publication, country, specimen type, histological subtypes, number of samples (both cases and controls), stage of cancer, array platforms, cut-off criteria (p-value or FDR), and numbers of DE miRNAs. All articles were independently assessed by two researchers (Chaniad P and Raungrut P).
Ranking
The vote-counting strategy reported in Griffith’s and Chan’s studies [14, 15] was used. The DE miRNAs were ranked for importance in the following order: 1) the number of articles that consistently reported as differentially expressed; 2) direction of change of DE miRNA, and 3) the total number of samples. Moreover, we considered the number of target genes supported by reliable and experimental evidence from the miRTarBase 7.0 [16].
Bioinformatics analysis
The miRTarbase 7.0 containing experimentally validated miRNA-target interactions was used to identify potential gene targets [16]. Enrichment analyses of Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were analyzed by DIANA-miRPath v3.0. [17]. The integrated genes regulatory network was generated by Cytoscape 3.5.1 [18]. Cytokine-related genes were retrieved from the Mouse Genome Informatics Web Site [19].
Study subjects and sample collection
This study was approved by the Ethics Committee on Human Research, Faculty of Medicine, Prince of Songkla University (REC 59-210-04-2 and REC 59-211-04-2). All subjects, including 32 NSCLC patients, 33 OL patients, and 23 HP, were recruited from Songklanagarind Hospital, Songkhla, Thailand, from May 2016 to April 2017. NSCLC was diagnosed according to 2015 WHO classification of lung and pleural tumors [20]. The staging of cancer was based on the Tumor Node Metastasis cancer staging system of the AJCC Cancer Staging Manual (7th Edition) [21]. No patients received any treatment before sample collection.
For the control groups, subjects who were age- and gender-matched with the NSCLC group were included. The OL group were patients who had symptoms similar to lung cancer, such as tuberculosis, bronchiectasis, interstitial lung disease, pneumonia, and chronic obstructive pulmonary disease. The HP group was recruited from people who had an annual checkup for two consecutive years. Written informed consent was obtained from all subjects. Peripheral blood (5 ml) was collected in a clotting tube (Greiner Bio-One, Kremsmünster, Austria). Serum was isolated by centrifugation and prepared for miRNA isolation, as previously described [8].
Quantitative real-time PCR (qRT-PCR)
The total RNA from serum was extracted using Trizol® LS (Invitrogen, California, USA) according to the manufacturer’s protocol. The quantity and quality of the total RNA were measured using a NanoDrop® ND-1000 UV-Vis spectrophotometer (Thermo Scientific, Massachusetts, USA) at an optical density ratio of A260/280 nm and A260/230 nm. The total RNA (50 ng) was reverse transcribed into complementary DNA (cDNA) by the miScript II RT kit (Qiagen, Hilden, Germany). Reverse transcription was performed by Thermal cycler (Bio-Rad, California, USA) at 37 °C for 60 mins, and 95 °C for 5 mins. The cDNA was diluted by adding 200 μl of RNase-free water. Subsequently, miRNA amplification was performed using a Bio-Rad CFX96 qPCR system (Bio-Rad, California, USA) with the miScript SYBR® green PCR kit (Qiagen, Hilden, Germany) according to the manufacturer’s protocol. Primers for miR-145-5p, miR-206-5p, miR-20a-5p, miR-223-5p, miR-25-3p, let-7f-5p, miR-20b-5p, and U6 small nuclear RNA (RNU6), were used as internal controls, and obtained from Qiagen (Qiagen, Hilden, Germany). Data were presented as the cycle threshold (Ct) value, which was determined using the default threshold settings. The difference between the Ct value (ΔCt) of the miRNA of interest and RNU6 was calculated, and the relative expression was presented using the 2−ΔCT method [22].
Enzyme-linked immunoassay
Concentrations of VEGF, IL-6, TGF-β, and TNF-α were measured quantitatively using 100 μl of serum with a sandwich human ELISA kit (Prepotech, Rehovot, Israel) according to the manufacturer’s instructions. Absorbance was read at 405 nm with a reference wavelength of 650 nm by a microplate reader (Molecular Devices, California, USA). Samples were assayed in duplicate, and concentrations of each cytokine were calculated by the constructed standard curve.
Statistical analysis
Data distribution of relative miRNA expression and cytokine concentration was assessed by the Shapiro–Wilk normality test. Outliers of the data were detected using the Grubbs’ test. Differences in the data between the groups of subjects were analyzed by the Student’s t-test if the data had a normal distribution, and by the Wilcoxon–Mann–Whitney test in cases of non-normal distribution. The correlation between miRNA expression and cytokine concentration was analyzed by the Spearman’s rank correlation. A p-value of less than 0.05 was considered statistically significant. All plots and statistical analyses were performed using the R-statistical software version 3.4.4.
Results
Included articles
Fig 1 depicts the number of selected articles from each selection process. A total of 725 studies were identified using our search terms. There were 536 remaining articles, after removing 189 duplicates. After screening the titles and abstracts, 47 eligible articles were included. According to the exclusion criteria, only 16 articles were finally included in our analysis (Fig 1). All characteristics of these articles are shown in Table 1.
Differentially expressed miRNAs
In total, 229 DE miRNAs were obtained from the 16 miRNA profiles studied, of which 181 and 48 miRNAs were reported with consistent and inconsistent direction (both up- and down-regulation), respectively. Among the 181 consistent miRNAs, 97 miRNAs were up-regulated, and 84 miRNAs were down-regulated.
Only miR-145 was reported to be consistently up-regulated in four studies. Six miRNAs (miR-320, miR-151a, miR-16, miR-200b-3p, miR-205, and miR-574) and one miRNA (miR-1285-3p) were most consistently reported to be up-regulated and down-regulated in three studies, respectively (Table 2). For inconsistent direction, miR-20a was reported in six studies (4 up-/2 down-regulated expression), and miRNAs (miR-25-3p, miR-223, let-7f, and miR-20b) were reported in five studies.
List of potential total genes and targeted genes that are supported by strong and experimental evidences was presented in S1 Table for the consistently up-regulated miRNA, S2 Table for the consistently down-regulated miRNA, and S3 Table for the inconsistency reported miRNA.
Functional enrichment analysis and gene interaction network
We performed a functional enrichment analysis to find the molecular networks and their target genes of seven DE miRNAs, including miR-145, miR-206, miR-20a, miR-223, miR-25-3p, let-7f, and miR-20b. The most significantly enriched GO terms on molecular functions, cellular components, and biological functions are iron binding (1662 genes), organelles (2973 genes), and cellular nitrogen metabolic processes (1499 genes) (Fig 2A). Also, the KEGG pathway showed that target genes were enriched in pathways for cancer (155 genes), PI3K-Akt signaling pathways (110 genes), proteoglycans in cancer (97 genes), MAPK signaling pathways (88 genes), and protein processing in the endoplasmic reticulum (Fig 2B). The complexity of the miRNA-mRNA interaction network is shown in Fig 2C. The results showed that a single miRNA regulates numerous genes. The highest number of target genes were identified for miR-145 (143 genes), followed by miR-20a (71 genes). In addition, these miRNAs work synergistically. We also revealed that our selected miRNAs regulated cytokine-associated genes, such as VEGFA, IGF1R, TGFBR1, and HIF1A.
Significantly enriched GO term (A) and KEGG pathway (B) of seven DE miRNAs. (C) Visualization of miRNAs and their associated target genes with Cytoscape. Red: miRNAs; Green: mRNAs; and yellow star: cytokine-associated genes.
Validation of selected differentially expressed miRNAs
The expression level of seven miRNAs, including miR-145, miR-206, miR-20a, miR-223, miR-25-3p, let-7f, and miR-20b, was validated the serum of 32 patients with NSCLC, 33 patients with OL, and 23 HP using qRT-PCR analysis. NSCLC patients had a mean age of 59.9 years. Eighteen cases were male, and 14 were female. Histological subtype was 84.4% (27/32 of cases) for adenocarcinoma (ADC) and 15.6% (5/32 of cases) for squamous cell carcinoma (SCC). The majority of patients exhibited stage IV lung cancer (90.6%, 29/32 of cases), whereas one and two cases were stage II and I lung cancer, respectively. The OL group included 22 patients with chronic obstructive pulmonary disease, 6 with bronchiectasis, 3 with tuberculosis, 1 with pneumonia, and 1 with idiopathic pulmonary fibrosis.
miR-145 was significantly down-regulated in NSCLC patients compared with OL patients (p = 0.003), HP (p = 0.038), and all controls (p = 0.004) (Fig 3). miR-20a was significantly up-regulated in NSCLC patients compared with OL patients (p = 0.045) and all controls (p = 0.039), but not with those in the HP groups. However, no significant difference was observed for other miRNAs.
(*) indicates significant difference of p < 0.05.
Serum cytokine concentration in NSCLC patients and controls
Concentrations of tumor-promoting cytokines, including IL-6, VEGF, TGF-β, and TNF-α were determined in the serum of NSCLC patients and controls using ELISA. The results revealed that NSCLC patients had significantly higher VEGF concentrations than those of OL and HP (p < 0.001) (Fig 4A). A significantly elevated concentration was also found in NSCLC patients compared with the HP group for IL-6 (p = 0.023) and TGF-β (p < 0.001), and all controls for IL-6 (p = 0.025) (Fig 4B and 4C). There was no significant difference in serum TNF-α concentration between NSCLC patients and controls (Fig 4D).
(*) indicates a significant difference of p < 0.05.
Correlation between miRNA expression and cytokine level
We evaluated the correlation between miRNA expression and cytokine concentration in the serum of NSCLC patients by Spearman’s correlation. A moderate correlation was shown between the concentration of TGF-β and the expression level of miR-20a (r = −0.537; p = 0.002) and miR-223 (r = 0.428; p = 0.015), in both negative and positive directions, respectively (Fig 5A and 5B). We also detected a weak, positive correlation of VEGF concentration with miR-20a expression (r = 0.376; p = 0.037) (Fig 5C), whereas miR-223 expression was weakly, negatively correlated with serum VEGF concentration (r = −0.355; p = 0.046) (Fig 5D). However, the other miRNAs did not exhibit any correlation with cytokines (Table 3).
(A) miR-20a versus TGF-β, (B) miR-223 versus TGF-β, (C) miR-20a versus VEGF, and (D) miR-223 versus VEGF. (*) indicates significant difference of p < 0.05 from Spearman’s rank correlation.
Discussion
In the present study, the lists of DE miRNAs were identified and then validated in clinical samples. Among 16 lung cancer miRNA profiling studies, eight miRNAs were consistently reported to be differentially expressed in at least three studies. In contrast, seven miRNAs were reported to be differentially expressed with inconsistent direction (either up- or down-regulation) in at least four studies. In the validation study, miR-145 and miR-20a showed a significant difference between NSCLC patients and controls. We also found that the cytokine concentrations of VEGF, IL-6, and TGF-β, but not TNF-α, were higher in NSCLC patients than in the control groups. In addition, the expression of miR-20a and miR-223 were correlated with the concentrations of TGF-β and VEGF in NSCLC patients.
Although several miRNA profiling studies to identify potential biomarkers for lung cancer have been conducted, the lists of DE miRNAs are not consistent. To deal with this issue, identification based on the vote-counting strategy, which has not yet been reported for lung cancer, was used. We demonstrated that the number of consistent DE miRNAs (79%, 181/229 miRNAs) was higher than those of inconsistent DE miRNAs (21%, 48/229 miRNAs). This result indicated a high degree of concordance among the studies. Inconsistent DE miRNAs could potentially be explained by several factors, such as the diversity of specimens, histological subtypes, clinical stages, and profiling platforms.
From the identification, miR-145 was consistently found to be up-regulated in four studies [26, 27, 29, 34]. However, our validation results in the clinical samples demonstrated the down-regulation of miR-145 in NSCLC patients compared with controls. Previous studies have shown that overexpression of miR-145 increased cell death of prostate cancer cells [39], and inhibited the invasion and metastasis of osteosarcoma cells [40]. Moreover, low expression levels of serum miR-145 predicted poor survival in patients with ovarian cancer [41]. In lung cancer, a decreased expression of miR-145 in tissues was found to have a significantly worse prognosis [42] and associated with a shorter time to relapse in NSCLC [43]. Our finding in clinical validation was consistent with these previous studies, indicating that the miR-145 might play a functional role as tumor-suppressive miRNA.
MiR-20a was found to be up-regulated in four studies [23, 27, 29, 34] and down-regulated in two studies [26, 32]. Our validation study supported the evidence for the up-regulation of miR-20a in NSCLC compared with controls. Although several studies have revealed the aberrant expression of miR-20a in various types of cancer, data on its function is inconsistent. For example, miR-20a overexpression inhibited the invasion of endometrial cancer cells [44] and the proliferation of neuroblastoma cells [45]. In contrast, miR-20a has been found to promote proliferation, invasion, and metastasis in hepatocellular carcinoma [46] and colorectal cancer cells [47]. Regarding clinical studies, increased expression of plasma/serum miR-20a was found in glioblastoma [48], gastric cancer [49], and NSCLC [50]. It was also associated with overall patient survival and disease-free survival. In addition, our bioinformatics analysis revealed that miR-20a regulated several cancer-related genes, particularly tumor suppressor genes, such as PTEN, BCL2L11, FBXO31, and DUSP2. This evidence was in accordance with our result in clinical samples that indicated an oncogenic property of miR-20a in lung cancer.
One important discrepancy between the results in our identification and clinical validation is the stage distribution of NSCLC patients. The studies of Gao et al. (2016) [24], Geng et al. (2014) [29], and Foss et al. (2011) [35] include early-stage NSCLC, while most of the patients in our study were in an advanced stage. This explanation was supported by the study of Aiso et al. [51], demonstrating that the serum expressions of miR-145 and miR-20a were higher in the early stages than the advanced stages of NSCLC.
Cytokines have been proposed as blood-based biomarkers and therapeutic targets, because they play essential roles in cancer development, including lung cancer [5]. The increased serum concentration of several cytokines, including VEGF, IL-6, TGF-β, and TNF-α, have been proposed as lung cancer biomarkers [5]. In the present study, we confirmed the increased serum concentration of all cytokines, except TNF-α in NSCLC compared with HP. Since the overestimation of the test could be found if only healthy persons are used as controls [52], both HP and patients with OL patients who may have symptoms similar to lung cancers were also included in this study for comparison.
There are accumulating studies showing that miRNAs associate with tumorigenesis through the regulation of cytokine production [11–13]. We found a moderate correlation between TGF-β concentration and expression levels of miR-20a and miR-223, and a weak correlation between VEGF concentration and expression levels of miR-20a and miR-223. Several studies have demonstrated that miR-20a regulates related genes in TGF-β-signaling, such as TGIF2, E2F5, MYC, ALK5, and TGFBR2 [11, 53]. In addition, some evidence has revealed that miR-20a is associated with VEGF production in breast cancer [12]. For miR-223, Berenstein et al. (2016) [10] have shown that decreased expression of miR-223 increases VEGF expression in bone marrow mesenchymal stromal cells. Liu et al. (2018) [54] have observed that miR-223 is increased in TGF-β1-stimulated cardiac fibroblasts. Our current results support the findings of theses previous studies, showing that miR-20a and miR-223 may regulate genes associated with TGF-β and VEGF.
There are some limitations to our present study. First, we included advanced-stage NSCLC for clinical validation. Since the DE miRNAs were identified from miRNA profiling conducted from various stages of NSCLC, the samples from patients with different stages should be more suitable than those from advanced stages of disease in our study. Second, the number of clinical samples, both NSCLC, and controls, is small. This low number may be subject to a type II error (false negative). Therefore, a large clinical sample should be used to assess the diagnostic performance of miR-145 and miR-20a for clinical utility.
Conclusions
We identified miR-145 and miR-20a as potential biomarkers for NSCLC. Our study provides supporting evidence of increased tumor-promoting cytokines in NSCLC, which are correlated with miR-20a and miR-223. Our results suggest that the regulation of tumor-promoting cytokines, through miR-20a and miR-223, might be a new therapeutic approach for lung cancer.
Supporting information
S1 Table. List of the genes associated with the consistently up-regulated miRNA.
https://doi.org/10.1371/journal.pone.0241593.s001
(PDF)
S2 Table. List of the genes associated with the consistently down-regulated miRNA.
https://doi.org/10.1371/journal.pone.0241593.s002
(PDF)
S3 Table. List of the genes associated with the inconsistency reported miRNA.
https://doi.org/10.1371/journal.pone.0241593.s003
(PDF)
Acknowledgments
Research center for cancer control in Thailand are acknowledged for their research facilities. We thank the office of International Affairs, Faculty of Medicine, Prince of Songkla University for English language editing.
References
- 1. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68: 394–424. pmid:30207593
- 2.
Howlader N, Noone AM, Krapcho M, Miller D, Brest A, Yu M, et al. SEER Cancer Statistics Review, 1975–2016 Bethesda, MD: National Cancer Institute; April 2019 [https://seer.cancer.gov/csr/1975_2016/].
- 3. Nakamura K, Ukawa S, Okada E, Hirata M, Nagai A, Yamagata Z, et al. Characteristics and prognosis of Japanese male and female lung cancer patients: The BioBank Japan Project. J Epidemiol. 2017;27: S49–s57. pmid:28202209
- 4. Song XY, Zhou SJ, Xiao N, Li YS, Zhen DZ, Su CY, et al. Research on the relationship between serum levels of inflammatory cytokines and non-small cell lung cancer. Asian Pac J Cancer Prev. 2013;14: 4765–4768. pmid:24083740
- 5. Marrugal A, Ojeda L, Paz-Ares L, Molina-Pinelo S, Ferrer I. Proteomic-Based Approaches for the Study of Cytokines in Lung Cancer. Dis Markers. 2016;2016: 2138627. pmid:27445423
- 6. Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116: 281–297. pmid:14744438
- 7. Fan J, Kou X, Yang Y, Chen N. MicroRNA-Regulated Proinflammatory Cytokines in Sarcopenia. Mediators Inflamm. 2016;2016: 1438686. pmid:27382188
- 8. Trakunram K, Champoochana N, Chaniad P, Thongsuksai P, Raungrut P. MicroRNA Isolation by Trizol-Based Method and Its Stability in Stored Serum and cDNA Derivatives. Asian Pac J Cancer Prev. 2019;20: 1641–1647. pmid:31244282
- 9. Tili E, Michaille JJ, Cimino A, Costinean S, Dumitru CD, Adair B, et al. Modulation of miR-155 and miR-125b levels following lipopolysaccharide/TNF-alpha stimulation and their possible roles in regulating the response to endotoxin shock. J Immunol. 2007;179: 5082–5089. pmid:17911593
- 10. Berenstein R, Nogai A, Waechter M, Blau O, Kuehnel A, Schmidt-Hieber M, et al. Multiple myeloma cells modify VEGF/IL-6 levels and osteogenic potential of bone marrow stromal cells via Notch/miR-223. Mol Carcinog. 2016;55: 1927–1939. pmid:27023728
- 11. Pellatt AJ, Mullany LE, Herrick JS, Sakoda LC, Wolff RK, Samowitz WS, et al. The TGFβ-signaling pathway and colorectal cancer: associations between dysregulated genes and miRNAs. J Transl Med. 2018;16: 191. pmid:29986714
- 12. Luengo-Gil G, Gonzalez-Billalabeitia E, Perez-Henarejos SA, Navarro Manzano E, Chaves-Benito A, Garcia-Martinez E, et al. Angiogenic role of miR-20a in breast cancer. PLoS One. 2018;13: e0194638. pmid:29617404
- 13. Yang Y, Ding L, Hu Q, Xia J, Sun J, Wang X, et al. MicroRNA-218 functions as a tumor suppressor in lung cancer by targeting IL-6/STAT3 and negatively correlates with poor prognosis. Mol Cancer. 2017;16: 141. pmid:28830450
- 14. Chan SK, Griffith OL, Tai IT, Jones SJ. Meta-analysis of colorectal cancer gene expression profiling studies identifies consistently reported candidate biomarkers. Cancer Epidemiol Biomarkers Prev. 2008;17: 543–552. pmid:18349271
- 15. Griffith OL, Melck A, Jones SJ, Wiseman SM. Meta-analysis and meta-review of thyroid cancer gene expression profiling studies identifies important diagnostic biomarkers. J Clin Oncol. 2006;24: 5043–5051. pmid:17075124
- 16. Chou CH, Shrestha S, Yang CD, Chang NW, Lin YL, Liao KW, et al. miRTarBase update 2018: a resource for experimentally validated microRNA-target interactions. Nucleic Acids Res. 2018;46: D296–d302. pmid:29126174
- 17. Vlachos IS, Paraskevopoulou MD, Karagkouni D, Georgakilas G, Vergoulis T, Kanellos I, et al. DIANA-TarBase v7.0: indexing more than half a million experimentally supported miRNA:mRNA interactions. Nucleic Acids Res. 2015;43: D153–159. pmid:25416803
- 18. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13: 2498–2504. pmid:14597658
- 19. Krupke DM, Begley DA, Sundberg JP, Richardson JE, Neuhauser SB, Bult CJ. The Mouse Tumor Biology Database: A Comprehensive Resource for Mouse Models of Human Cancer. Cancer Res. 2017;77: e67–e70. pmid:29092943
- 20. Travis WD, Brambilla E, Nicholson AG, Yatabe Y, Austin JHM, Beasley MB, et al. The 2015 World Health Organization Classification of Lung Tumors: Impact of Genetic, Clinical and Radiologic Advances Since the 2004 Classification. J Thorac Oncol. 2015;10: 1243–1260. pmid:26291008
- 21. Edge SB, Compton CC. The American Joint Committee on Cancer: the 7th edition of the AJCC cancer staging manual and the future of TNM. Ann Surg Oncol. 2010;17: 1471–1474. pmid:20180029
- 22. Rao X, Huang X, Zhou Z, Lin X. An improvement of the 2^(-delta delta CT) method for quantitative real-time polymerase chain reaction data analysis. Biostat Bioinforma Biomath. 2013;3: 71–85. pmid:25558171
- 23. Fan L, Qi H, Teng J, Su B, Chen H, Wang C, et al. Identification of serum miRNAs by nano-quantum dots microarray as diagnostic biomarkers for early detection of non-small cell lung cancer. Tumour Biol. 2016;37: 7777–7784. pmid:26695145
- 24. Gao X, Wang Y, Zhao H, Wei F, Zhang X, Su Y, et al. Plasma miR-324-3p and miR-1285 as diagnostic and prognostic biomarkers for early stage lung squamous cell carcinoma. Oncotarget. 2016;7: 59664–59675. pmid:27517633
- 25. Halvorsen AR, Bjaanaes M, LeBlanc M, Holm AM, Bolstad N, Rubio L, et al. A unique set of 6 circulating microRNAs for early detection of non-small cell lung cancer. Oncotarget. 2016;7: 37250–37259. pmid:27191990
- 26. Leidinger P, Brefort T, Backes C, Krapp M, Galata V, Beier M, et al. High-throughput qRT-PCR validation of blood microRNAs in non-small cell lung cancer. Oncotarget. 2016;7: 4611–4623. pmid:26672767
- 27. Nadal E, Truini A, Nakata A, Lin J, Reddy RM, Chang AC, et al. A Novel Serum 4-microRNA Signature for Lung Cancer Detection. Sci Rep. 2015;5: 12464. pmid:26202143
- 28. Wozniak MB, Scelo G, Muller DC, Mukeria A, Zaridze D, Brennan P. Circulating MicroRNAs as Non-Invasive Biomarkers for Early Detection of Non-Small-Cell Lung Cancer. PLoS One. 2015;10: e0125026. pmid:25965386
- 29. Geng Q, Fan T, Zhang B, Wang W, Xu Y, Hu H. Five microRNAs in plasma as novel biomarkers for screening of early-stage non-small cell lung cancer. Respir Res. 2014;15: 149. pmid:25421010
- 30. Rani S, Gately K, Crown J, O’Byrne K, O’Driscoll L. Global analysis of serum microRNAs as potential biomarkers for lung adenocarcinoma. Cancer Biol Ther. 2013;14: 1104–1112. pmid:24025412
- 31. Heegaard NH, Schetter AJ, Welsh JA, Yoneda M, Bowman ED, Harris CC. Circulating micro-RNA expression profiles in early stage nonsmall cell lung cancer. Int J Cancer. 2012;130: 1378–1386. pmid:21544802
- 32. Patnaik SK, Yendamuri S, Kannisto E, Kucharczuk JC, Singhal S, Vachani A. MicroRNA expression profiles of whole blood in lung adenocarcinoma. PLoS One. 2012;7: e46045. pmid:23029380
- 33. Roth C, Stuckrath I, Pantel K, Izbicki JR, Tachezy M, Schwarzenbach H. Low levels of cell-free circulating miR-361-3p and miR-625* as blood-based markers for discriminating malignant from benign lung tumors. PLoS One. 2012;7: e38248. pmid:22675530
- 34. Chen X, Hu Z, Wang W, Ba Y, Ma L, Zhang C, et al. Identification of ten serum microRNAs from a genome-wide serum microRNA expression profile as novel noninvasive biomarkers for nonsmall cell lung cancer diagnosis. Int J Cancer. 2012;130: 1620–1628. pmid:21557218
- 35. Foss KM, Sima C, Ugolini D, Neri M, Allen KE, Weiss GJ. miR-1254 and miR-574-5p: serum-based microRNA biomarkers for early-stage non-small cell lung cancer. J Thorac Oncol. 2011;6: 482–488. pmid:21258252
- 36. Silva J, Garcia V, Zaballos A, Provencio M, Lombardia L, Almonacid L, et al. Vesicle-related microRNAs in plasma of nonsmall cell lung cancer patients and correlation with survival. Eur Respir J. 2011;37: 617–623. pmid:20595154
- 37. Wang ZX, Bian HB, Wang JR, Cheng ZX, Wang KM, De W. Prognostic significance of serum miRNA-21 expression in human non-small cell lung cancer. J Surg Oncol. 2011;104: 847–851. pmid:21721011
- 38. Keller A, Leidinger P, Borries A, Wendschlag A, Wucherpfennig F, Scheffler M, et al. miRNAs in lung cancer—Studying complex fingerprints in patient’s blood cells by microarray experiments. BMC Cancer. 2009;9: 353. pmid:19807914
- 39. Chen X, Gong J, Zeng H, Chen N, Huang R, Huang Y, et al. MicroRNA145 targets BNIP3 and suppresses prostate cancer progression. Cancer Res. 2010;70: 2728–2738. pmid:20332243
- 40. Fan L, Wu Q, Xing X, Wei Y, Shao Z. MicroRNA-145 targets vascular endothelial growth factor and inhibits invasion and metastasis of osteosarcoma cells. Acta Biochim Biophys Sin (Shanghai). 2012;44: 407–414. pmid:22472569
- 41. Liang H, Jiang Z, Xie G, Lu Y. Serum microRNA-145 as a novel biomarker in human ovarian cancer. Tumour Biol. 2015;36: 5305–5313. pmid:25722112
- 42. Yanaihara N, Caplen N, Bowman E, Seike M, Kumamoto K, Yi M, et al. Unique microRNA molecular profiles in lung cancer diagnosis and prognosis. Cancer Cell. 2006;9: 189–198. pmid:16530703
- 43. Campayo M, Navarro A, Vinolas N, Diaz T, Tejero R, Gimferrer JM, et al. Low miR-145 and high miR-367 are associated with unfavourable prognosis in resected nonsmall cell lung cancer. Eur Respir J. 2013;41: 1172–1178. pmid:22835608
- 44. Huang Y, Yang N. MicroRNA-20a-5p inhibits epithelial to mesenchymal transition and invasion of endometrial cancer cells by targeting STAT3. Int J Clin Exp Pathol 2018;11: 5715–5724. pmid:31949657
- 45. Yu Y, Zhang J, Jin Y, Yang Y, Shi J, Chen F, et al. MiR-20a-5p suppresses tumor proliferation by targeting autophagy-related gene 7 in neuroblastoma. Cancer Cell Int. 2018;18: 5. pmid:29311760
- 46. Chen Y, Wang X, Cheng J, Wang Z, Jiang T, Hou N, et al. MicroRNA-20a-5p targets RUNX3 to regulate proliferation and migration of human hepatocellular cancer cells. Oncol Rep. 2016;36: 3379–3386. pmid:27748919
- 47. Cheng D, Zhao S, Tang H, Zhang D, Sun H, Yu F, et al. MicroRNA-20a-5p promotes colorectal cancer invasion and metastasis by downregulating Smad4. Oncotarget. 2016;7: 45199–45213. pmid:27286257
- 48. Zhao H, Shen J, Hodges TR, Song R, Fuller GN, Heimberger AB. Serum microRNA profiling in patients with glioblastoma: a survival analysis. Mol Cancer. 2017;16: 59. pmid:28284220
- 49. Yang R, Fu Y, Zeng Y, Xiang M, Yin Y, Li L, et al. Serum miR-20a is a promising biomarker for gastric cancer. Biomed Rep. 2017;6: 429–434. pmid:28413641
- 50. Xu X, Zhu S, Tao Z, Ye S. High circulating miR-18a, miR-20a, and miR-92a expression correlates with poor prognosis in patients with non-small cell lung cancer. Cancer Med. 2018;7: 21–31. pmid:29266846
- 51. Aiso T, Ohtsuka K, Ueda M, Karita S, Yokoyama T, Takata S, et al. Serum levels of candidate microRNA diagnostic markers differ among the stages of non-small-cell lung cancer. Oncol Lett. 2018;16: 6643–6651. pmid:30405804
- 52. Rutjes AW, Reitsma JB, Vandenbroucke JP, Glas AS, Bossuyt PM. Case-control and two-gate designs in diagnostic accuracy studies. Clin Chem. 2005;51: 1335–1341. pmid:15961549
- 53. Correia AC, Moonen JR, Brinker MG, Krenning G. FGF2 inhibits endothelial-mesenchymal transition through microRNA-20a-mediated repression of canonical TGF-beta signaling. J Cell Sci. 2016;129: 569–579. pmid:26729221
- 54. Liu X, Xu Y, Deng Y, Li H. MicroRNA-223 Regulates Cardiac Fibrosis After Myocardial Infarction by Targeting RASA1. Cell Physiol Biochem. 2018;46: 1439–1454. pmid:29689569