A transcriptional miRNA-gene network associated with lung adenocarcinoma metastasis based on the TCGA database
- Authors:
- Published online on: January 14, 2016 https://doi.org/10.3892/or.2016.4560
- Pages: 2257-2269
Abstract
Introduction
Lung adenocarcinoma, the most common subtype of non-small cell lung cancer (NSCLC), leads to one million deaths each year, affecting an increasing percentage of the population over the past few years (1). Despite advances in suitable therapies, the 5-year survival of lung adenocarcinoma patients remains low. The high mortality rate may be attributed to the delay of detection, since patients are asymptomatic in early stages. Local and distant metastases occur in most cases by the time symptoms are obvious, resulting in treatment failure in advanced NSCLC (2). Therefore, it is of great importance to explore the mechanism of NSCLC metastasis and strategies to block metastasis of the disease.
As post-transcriptional modulators, microRNAs (miRNAs) are a class of endogenous, non-coding, single-stranded RNAs with 21–24 nucleotides (3). Dysregulation of miRNAs contributes to many pathological conditions, such as the initiation and progression of lung cancer. A number of studies have assessed the potential role of miRNA signatures to classify histological subtypes (4) or to predict diagnosis (5), metastasis, recurrence or survival of NSCLC patients (6–9). Wang et al revealed a set of 24 differentially expressed miRNAs between NSCLC metastatic primary loci and non-metastatic primary loci in a mouse model (9). Larzabal et al showed that miR-205 overexpression inhibited metastasis of NSCLC by targeting integrin α5 (a pro-invasive protein) upon TMPRSS4 blockade, which was confirmed by in vivo and in vitro experiments (10).
High throughput data, such as gene expression data from RNA sequencing or microarrays, could be widely used in the exploration of molecular mechanisms that drive tumor behavior. The Cancer Genome Atlas (TCGA) database, a publically available database, offers a multilayered view of genomic and epigenomic data of approximately 10,000 patient samples together with clinicopathological information across more than 30 human cancer types, which is a rich resource for data mining and biological discovery.
Based on the fact that the collection of lung adenocarcinoma specimens that have metastasized to other organs is difficult, investigation of the early molecular events underlying lung adenocarcinoma metastasis is difficult. Toward this end, we analyzed mRNA and miRNA expression data and clinical data derived from TCGA to identify differences in mRNA and miRNA expression in primary tumor tissues from lung adenocarcinoma that did and did not metastasize. In addition, combined with bioinformatic prediction, we constructed an miRNA-target gene regulatory network, suggesting the regulation of the beginning of lung adenocarcinoma metastasis.
Materials and methods
TCGA gene expression profiles
Lung adenocarcinoma level 3 mRNA and miRNA expression data, and the corresponding clinical information, were downloaded from the TCGA data portal (https://tcga-data.nci.nih.gov/tcga/). The tumor staging information for all patient samples was derived from TCGA clinical information. Gene expression data were available for 291 lung adenocarcinoma samples without metastases, and 248 lung adenocarcinoma samples with lymph node metastasis or distant metastases. The miRNA expression data were available for 186 lung adenocarcinoma samples without metastases, and 144 lung adenocarcinoma with metastases.
Ranking of differentially expressed genes and miRNAs
The raw expression data of all lung adenocarcinoma patients was downloaded, and transformed into log2 scale. Z-score normalization was also employed. The Limma (Linear Models for Microarray Data) package in R was used to identify the differentially expressed probe sets between the metastatic primary loci and non-metastatic primary loci by two-tailed Student's t-test, and p-values from the same platform were obtained. MetaMA package in R was used to combine p-values from different platforms, and the false discovery rate (FDR) was calculated for multiple comparisons using the Benjamini and Hochberg method. We selected differentially expressed mRNAs and miRNAs with criterion of FDR <0.05. Hierarchical clustering of differentially expressed genes was performed using the 'heatmap.2' function of the R/Bioconductor package 'gplots' (11).
Target gene prediction of differentially expressed miRNAs
To understand the potential association between differentially expressed mRNAs and miRNAs obtained in the study, we predicted the transcriptional targets of the identified miRNAs using the online tools of miRWalk (http://www.umm.uni-heidelberg.de/apps/zmf/mirwalk/) (12) based on six bioinformatic algorithms (DIANAmT, miRanda, miRDB, miRWalk, PicTar and TargetScan). Target genes that were commonly predicted by ≥4 algorithms or experimentally validated based on miRWalk database, were considered as putative targets.
Construction of the regulatory network of miRNA-target genes
Given that miRNAs tend to decrease the expression of their target genes, we matched putative target genes with the list of differentially expressed genes between the metastatic primary loci and non-metastatic primary loci to increase the accuracy of target prediction. We selected miRNA-target pairs whose expression was inversely correlated, to subject to further investigation (13–15). We conducted miRNA-target gene interaction networks with miRNA-target gene interacting pairs, whose expression levels are inversely correlated, and the miRNA regulation networks were visualized by Cytoscape (16).
Functional annotation
Functional enrichment analysis is essential to uncover biological functions of miRNA target genes. To gain insight into the biological function of the miRNA target genes, we performed Gene Ontology (GO) classification and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis based on the online software GENECODIS (17). GO which includes three categories such as biological process, molecular function and cellular component, provides a common descriptive framework of gene annotation and classification for analyzing gene set data. The KEGG pathway enrichment analysis was also performed to detect the potential pathway of miRNA target genes based on the KEGG pathway database, which is a recognized and comprehensive database including all types of biochemistry pathways (18). FDR <0.05 was set as the cut-off for selecting significantly enriched functional GO terms and KEGG pathways.
Validating the expression of lung adenocarcinoma metastasis-associated miRNA target genes
Based on the published studies, hsa-mir-133a-1 and hsa-let-7d were closely linked to the metastatic ability of lung adenocarcinoma (19,20). In our study, SurvMicro (http://bioinformatica.mty.itesm.mx:8080/Biomatec/Survmicro.jsp), an online validation tool, was used to evaluate the association of miRNA expression with survival (21) in human cancer datasets.
Furthermore, the expression of putative targets of hsa-mir-133a-1 (PKM2, PTBP1, FSCN1 and STRAP) and hsa-let-7d (KIT, BCL2, CNTN3, CISH, FLT3 and ESR2) was detected in primary tumor tissues from 10 lung adenocarcinoma patients by RT-qPCR, to validate the differential expression observed in the bioinformatic analysis. The 2−ΔΔCt method was used to analyze the data of RT-qPCR. SPSS version 13.0 was used to perform all statistical analyses. The sequences of the primers used are provided in Table I.
Among the 10 patients included, 5 presented with metastasis and 5 were without metastasis. The clinical specimens were provided by Daping Hospital. Written informed consent forms were received from the patients or legal guardians of the patients. All protocols and documents were approved by the Medical Ethics Committee of Daping Hospital. Frozen sections were prepared from the tumor tissues for the cytological or histological diagnosis. The resected tumor tissues were stored in liquid nitrogen until RNA extraction.
Results
Differences in miRNA and mRNA expression in lung adenocarcinoma with or without metastasis
Based on clinical information 'AJCC pathologic tumor stage' in the TCGA lung adenocarcinoma, we classified these data into lung adenocarcinoma with or without metastasis. After the gene expression data were downloaded, we performed differential expression analysis of the genes and miRNAs between the metastatic primary loci and non-metastatic primary loci. Genes (1,257) were identified to be differentially expressed under a threshold of FDR <0.05, with 585 upregulated and 672 downregulated genes in the lung adenocarcinoma with metastasis. The hierarchical clustering of differentially expressed genes is shown in Fig. 1. Fifty-eight miRNAs were identified as differentially expressed under the threshold of FDR <0.05, with 44 upregulated and 14 downregulated miRNAs in the lung adenocarcinoma samples with metastasis (Table II).
The regulatory network of differentially expressed mRNAs and miRNAs associated with lung adenocarcinoma metastasis
To find the potential link between the differentially expressed mRNAs and miRNAs observed in the bioinformatic analysis, we predicted the putative targets of the identified miRNAs in the miRWalk database. Matching predicted putative targets with those found to be disregulated genes in metastatic primary tumors, miRNA-target gene pairs whose expression levels were inversely correlated were selected. As a result, we identified 1,154 miRNA-target gene pairs for the upregulated miRNAs with 61 pairs validated by previous experiments, and 474 miRNA-target gene pairs for the downregulated miRNAs with 57 pairs validated by previous experiments (Table III).
Using the 770 miRNA-target gene pairs, an miRNA-target gene regulatory network was constructed (Fig. 2). In the miRNA-target gene regulatory network, we identified the top 10 miRNAs which regulated the most target genes, such as hsa-miR-7, hsa-miR-182, hsa-miR-324-3p, hsa-miR-139-5p, hsa-miR-130b, hsa-let-7f, hsa-miR-18a, hsa-miR-188-5p, hsa-let-7d, and hsa-miR-590-5p, and the target genes such as RPS6KA3, TSC1, AIM1, GAS7, GFOD1, GGA2, IGF1, IL28RA, and INSR were regulated by the most miRNAs.
GO classification and KEGG pathways of the miRNA target genes
We performed the GO classification and KEGG pathway enrichment analysis for miRNA target genes whose expression was differentially expressed. We found that lymphoid progenitor cell differentiation (GO:0002320, FDR=4.75E-04) and hematopoietic progenitor cell differentiation (GO:0002244, FDR=2.08E-03) were significantly enriched for biological processes. While for molecular functions, signaling receptor activity (GO:0038023, FDR=9.24E-04) and signal transducer activity (GO:0004871, FDR=8.77E-04) were significantly enriched (Table IV). The most significant pathway in our KEGG analysis was DNA replication (FDR=1.78E-04). Furthermore, T cell receptor signaling pathway (FDR=1.12E-03) and endocytosis (FDR=1.17E-03) were also highly enriched (Table V).
Validating the expression of lung adenocarcinoma metastasis-associated miRNA target genes
By the online tool of SurvMicro, survival analysis was performed to evaluate the correlation between miRNA expression level (hsa-miR-133a-1 and hsa-let-7d) and the overall survival time of the lung adenocarcinoma patients. With the data from TCGA database, Kaplan-Meier curves indicated that the expression of hsa-mir-133a-1 was significantly correlated with the overall survival time of lung adenocarcinoma patients (P=0.03654), while the expression of hsa-let-7d was not significantly correlated with the overall survival time of the lung adenocarcinoma patients (P=0.114) (Fig. 3).
We performed RT-qPCR for the putative target genes of hsa-mir-133a-1 and hsa-let-7d, which exhibited significantly differential expression between the primary tumor tissues from lung adenocarcinoma with metastasis and that from the lung adenocarcinoma patiemts without metastasis in the bioinformatic analysis. In the RT-qPCR assay, most of the target genes showed similar expression patterns to what was observed in the bioinformatic analysis. The expression levels of the selected upregulated and downregulated genes are shown in Fig. 4. Among the upregulated genes in the lung adenocarcinoma with metastasis, PKM2 and STRAP mRNA expression was significantly increased (P=0.04, P=0.05), and PTBP1 mRNA expression was also increased, but not obviously, while FSCN1 mRNA expression was decreased. Among the downregulated genes in the lung adenocarcinoma with metastasis, FLT3 mRNA expression was significantly decreased (P=0.023), and BCL2, ESR2, and KIT mRNA expression was also mildly decreased. In contrast, CISH mRNA expression was significantly increased (P=0.027). CNTN3 mRNA expression appeared unchanged between the metastatic primary loci and the non-metastatic primary loci.
Discussion
Due to a substantially lower cost, high-throughput gene expression data are becoming widely available. A key question in the field is how to use the data to identify both biological drivers and strong metastatic markers. In this study, we compared miRNA and mRNA expression data differences between lung adenocarcinoma with metastasis and lung adenocarcinoma without metastasis from the TCGA data portal with the aim to screen the miRNAs and genes associated with lung adenocarcinoma metastasis.
In line with previous findings, some of the identified differentially expressed miRNAs were found to be implicated in lung adenocarcinoma metastasis, such as hsa-mir-98, hsa-let-7d, hsa-mir-134, hsa-mir-196b, hsa-mir-381, hsa-mir-133a, hsa-mir-552, hsa-mir-944, hsa-mir-550, and hsa-mir-655. miR-98 was also found to be overexpressed in lung cancer cell lines, and it binds the 3′UTR of Fus1, a tumor suppressor, to inhibit protein expression (22). Despite no reports related to lung adenocarcinoma metastasis, a recent study uncovered the regulatory roles of miR-98 in melanoma metastasis (23). miR-134 and miR-655, belonging to the same cluster on chromosome 14q32, were demonstrated to be involved in TGF-β1-induced EMT which is recognized as a key element of cell invasion, migration, and metastasis, by directly targeting MAGI2, a scaffold protein required for PTEN (24). Directly targeting oncogenic receptors, such as IGF-1R, TGFBR1, and EGFR, miR-133a inhibits cell invasiveness and cell growth in lung cancer cells. Furthermore, an in vivo animal model showed that miR-133a markedly suppressed the metastatic ability of lung cancer cells (19). An miRNA microarray revealed that the expression level of miR-552 in colorectal cancer metastases in the lung was 39 times higher than that of primary lung adenocarcinomas, suggesting its possible roles in cancer metastases (25). miR-944 affected NSCLC cell growth, proliferation, and invasion by targeting a tumor suppressor, SOCS4. miR-944 was found to be associated with lymph node metastasis by determining its expression in 52 formalin-fixed paraffin-embedded SCC tissues (26).
More importantly, two miRNAs, hsa-mir-133a and hsa-let-7d, were found to be involved in lung adenocarcinoma metastasis by previous studies. Consequently, we performed survival analysis to evaluate the correlation between miRNA expression level and the overall survival time of lung adenocarcinoma patients, and we found that the expression of hsa-mir-133a was significantly correlated with the overall survival of lung adenocarcinoma, while the expression of hsa-let-7d was not significantly correlated with the survival of lung adenocarcinoma.
miR-133a was identified to function as a tumor suppressor in NSCLCs directly targeting several membrane receptors including IGF-1R, TGFBR1 and EGFR. In addition, by using the in vivo animal model, ectopically expressing miR-133a markedly suppressed the metastatic ability of the lung cancer cells, which may be applied in the clinical therapy of lung cancer in the future. Accordingly, NSCLC patients with higher expression levels of miR-133a had longer survival rates compared with those with lower miR-133a expression levels (19).
Let-7d, one of the members of the let-7 family, is deregulated in many types of cancers (27–30). Let-7d is involved in cellular differentiation, epithelial-to-mesenchymal transition (EMT) as a switch between EMT and MET (mesenchymal-epithelial transition), and regulates tumor-initiating cell (TIC) formation (31). Mairinger et al showed that let-7d expression was significantly associated with overall survival in pulmonary neuroendocrine tumors, suggesting its important role in metastasis and tumor progression (20).
The putative targets of hsa-mir-133a and hsa-let-7d were subjected to RT-qPCR validation in primary tumor tissues from 5 lung adenocarcinoma patients with metastasis and 5 lung adenocarcinoma patients without metastasis. The result of the RT-qPCR assay revealed that most of the genes showed similar expression patterns to what were observed in the bioinformatic analysis, suggesting that the findings in the bioinformatic analysis of the differential gene expression data were credible. PKM2, as a common putative target of both hsa-mir-133a and hsa-mir-552, was significantly upregulated at the mRNA level in the lung adenocarcinoma with metastasis. Interesting, the protein level of PKM2 was upregulated in the lung adenocarcinoma tissues compared with the paired surrounding normal tissue, which was also correlated with chemotherapy resistance, the severity of epithelial dysplasia, as well as a relatively poor prognosis (32,33), indicating that PKM2 could be used as a tumor marker for diagnosis and, in particular, as a potential target for cancer therapy.
STRAP, a putative target of hsa-mir-133a, was significantly increased at the mRNA level in lung adenocarcinoma with metastasis. As a WD domain-containing protein, STRAP inhibits TGF-β signaling through interaction with receptors and Smad7, and promotes growth and enhances tumorigenicity. Similarly, strong upregulation of STRAP was also observed in lung tumors by immunoblot analyses (34).
FLT3, a common putative target of hsa-let-7d, hsa-mir-196b and hsa-mir-376c, was significantly upregulated at the mRNA level among the downregulated genes in the lung adenocarcinoma with metastasis. It was found that in a mouse model of lung metastases treatment with the Flt3 ligand significantly reduced the number of lung metastases after laparotomy or radiation therapy, and thus prolonged survival (35,36). Unlike the findings in the bioinformatic analysis, CISH mRNA expression was significantly increased. As negative feedback regulators of JAK-STAT and several other signalling pathways, CISH, a putative common target of hsa-let-7d, hsa-mir-382, hsa-mir-9-3, hsa-mir-9-2 and hsa-mir-9-1, was recently found to be involved in the development of many solid organ and haematological malignancies (37). Prostate cancer LNCaP-S17 cells were resistant to exogenous IL-6-induced neuroendocrine differentiation due to increased levels of CISH and SOCS7 that blocked activation of the JAK2-STAT3 pathways (38).
Functional enrichment analysis of miRNA target genes showed that those genes were highly correlated with carcinogenesis. The most significantly enriched pathway based on the KEGG database was DNA replication. Based on the fact that cancer cells are characterized by uncontrolled cell proliferation properties, deregulation of DNA replication may promote the process of carcinogenesis. In addition, in cancer cells several DNA replication-initiation proteins, such as CDC6 and minichromosome maintenance (MCM) proteins, were overexpressed (39,40).
There are two limitations in this study. Firstly, the expression levels of the miRNAs were not determined in the metastatic primary loci compared with non-metastatic primary loci of lung adenocarcinomas. Secondly, the small size of the clinical samples for RT-qPCR validation was also a limitation of our study, which may be the reason that the RT-qPCR results were not significant. Given that surgery is not recommended for the traditional stage IIIB-IV patients with lung adenocarcinoma classified according to the taxonomy of cancer staging (TNM) system, it is difficult to obtain specimens of metastatic lung adenocarcinoma. Thus, detailed further studies are needed to investigate the selected miRNAs and the corresponding target genes in lung adenocarcinomas metastasis.
In summary, by comparing the difference in the miRNA and mRNA expression profiling in lung adenocarcinomas with and without metastases, an miRNA-regulated network involved in lung adenocarcinoma metastasis was identified combined with the bioinformatic analysis. RT-qPCR results indicated that the bioinformatic approach in our study was acceptable. Our results may contribute to the identification of markers which could be used to detect lung adenocarcinoma metastasis and assess prognosis.
References
Alberg AJ, Ford JG and Samet JM; American College of Chest Physicians: Epidemiology of lung cancer: ACCP evidence-based clinical practice guidelines (2nd edition). Chest. 132(Suppl 3): 29S–55S. 2007. View Article : Google Scholar : PubMed/NCBI | |
Xie L, Yang Z, Li G, Shen L, Xiang X, Liu X, Xu D, Xu L, Chen Y, Tian Z, et al: Genome-wide identification of bone metastasis-related microRNAs in lung adenocarcinoma by high-throughput sequencing. PLoS One. 8:e612122013. View Article : Google Scholar : PubMed/NCBI | |
Bartel DP: MicroRNAs: Genomics, biogenesis, mechanism, and function. Cell. 116:281–297. 2004. View Article : Google Scholar : PubMed/NCBI | |
Zhang YK, Zhu WY, He JY, Chen DD, Huang YY, Le HB and Liu XG: miRNAs expression profiling to distinguish lung squamous-cell carcinoma from adenocarcinoma subtypes. J Cancer Res Clin Oncol. 138:1641–1650. 2012. View Article : Google Scholar : PubMed/NCBI | |
Zheng D, Haddadin S, Wang Y, Gu LQ, Perry MC, Freter CE and Wang MX: Plasma microRNAs as novel biomarkers for early detection of lung cancer. Int J Clin Exp Pathol. 4:575–586. 2011.PubMed/NCBI | |
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. 16:430–441. 2010. View Article : Google Scholar : PubMed/NCBI | |
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. 9:189–198. 2006. View Article : Google Scholar : PubMed/NCBI | |
Yu SL, Chen HY, Chang GC, Chen CY, Chen HW, Singh S, Cheng CL, Yu CJ, Lee YC, Chen HS, et al: MicroRNA signature predicts survival and relapse in lung cancer. Cancer Cell. 13:48–57. 2008. View Article : Google Scholar : PubMed/NCBI | |
Wang R, Chen XF and Shu YQ: Prediction of non-small cell lung cancer metastasis-associated microRNAs using bioinformatics. Am J Cancer Res. 5:32–51. 2015.PubMed/NCBI | |
Larzabal L, de Aberasturi AL, Redrado M, Rueda P, Rodriguez MJ, Bodegas ME, Montuenga LM and Calvo A: TMPRSS4 regulates levels of integrin α5 in NSCLC through miR-205 activity to promote metastasis. Br J Cancer. 110:764–774. 2014. View Article : Google Scholar : PubMed/NCBI | |
Reimers M and Carey VJ: Bioconductor: An open source framework for bioinformatics and computational biology. Methods Enzymol. 411:119–134. 2006. View Article : Google Scholar : PubMed/NCBI | |
Dweep H, Sticht C, Pandey P and Gretz N: miRWalk - database: Prediction of possible miRNA binding sites by 'walking' the genes of three genomes. J Biomed Inform. 44:839–847. 2011. View Article : Google Scholar : PubMed/NCBI | |
Luo Z, Zhang L, Li Z, Li X, Li G, Yu H, Jiang C, Dai Y, Guo X, Xiang J, et al: An in silico analysis of dynamic changes in microRNA expression profiles in stepwise development of naso-pharyngeal carcinoma. BMC Med Genomics. 5:32012. View Article : Google Scholar | |
Lionetti M, Biasiolo M, Agnelli L, Todoerti K, Mosca L, Fabris S, Sales G, Deliliers GL, Bicciato S, Lombardi L, et al: Identification of microRNA expression patterns and definition of a microRNA/mRNA regulatory network in distinct molecular groups of multiple myeloma. Blood. 114:e20–e26. 2009. View Article : Google Scholar : PubMed/NCBI | |
Enerly E, Steinfeld I, Kleivi K, Leivonen SK, Aure MR, Russnes HG, Rønneberg JA, Johnsen H, Navon R, Rødland E, et al: miRNA-mRNA integrated analysis reveals roles for miRNAs in primary breast tumors. PLoS One. 6:e169152011. View Article : Google Scholar : PubMed/NCBI | |
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B and Ideker T: Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 13:2498–2504. 2003. View Article : Google Scholar : PubMed/NCBI | |
Tabas-Madrid D, Nogales-Cadenas R and Pascual-Montano A: GeneCodis3: A non-redundant and modular enrichment analysis tool for functional genomics. Nucleic Acids Res. 40:W478–W483. 2012. View Article : Google Scholar : PubMed/NCBI | |
Altermann E and Klaenhammer TR: PathwayVoyager: Pathway mapping using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. BMC Genomics. 6:602005. View Article : Google Scholar : PubMed/NCBI | |
Wang LK, Hsiao TH, Hong TM, Chen HY, Kao SH, Wang WL, Yu SL, Lin CW and Yang PC: MicroRNA-133a suppresses multiple oncogenic membrane receptors and cell invasion in non-small cell lung carcinoma. PLoS One. 9:e967652014. View Article : Google Scholar : PubMed/NCBI | |
Mairinger FD, Ting S, Werner R, Walter RF, Hager T, Vollbrecht C, Christoph D, Worm K, Mairinger T, Sheu-Grabellus SY, et al: Different micro-RNA expression profiles distinguish subtypes of neuroendocrine tumors of the lung: results of a profiling study. Mod Pathol. 27:1632–1640. 2014. View Article : Google Scholar : PubMed/NCBI | |
Aguirre-Gamboa R, Gomez-Rueda H, Martínez-Ledesma E, Martínez-Torteya A, Chacolla-Huaringa R, Rodriguez-Barrientos A, Tamez-Peña JG and Treviño V: SurvExpress: An online biomarker validation tool and database for cancer gene expression data using survival analysis. PLoS One. 8:e742502013. View Article : Google Scholar : PubMed/NCBI | |
Du L, Schageman JJ, Subauste MC, Saber B, Hammond SM, Prudkin L, Wistuba II, Ji L, Roth JA, Minna JD, et al: miR-93, miR-98, and miR-197 regulate expression of tumor suppressor gene FUS1. Mol Cancer Res. 7:1234–1243. 2009. View Article : Google Scholar : PubMed/NCBI | |
Li F, Li XJ, Qiao L, Shi F and Liu W, Li Y, Dang YP, Gu WJ, Wang XG and Liu W: miR-98 suppresses melanoma metastasis through a negative feedback loop with its target gene IL-6. Exp Mol Med. 46:e1162014. View Article : Google Scholar : PubMed/NCBI | |
Kitamura K, Seike M, Okano T, Matsuda K, Miyanaga A, Mizutani H, Noro R, Minegishi Y, Kubota K and Gemma A: miR-134/487b/655 cluster regulates TGF-β-induced epithelial-mesenchymal transition and drug resistance to gefitinib by targeting MAGI2 in lung adenocarcinoma cells. Mol Cancer Ther. 13:444–453. 2014. View Article : Google Scholar | |
Kim J, Lim NJ, Jang SG, Kim HK and Lee GK: miR-592 and miR-552 can distinguish between primary lung adenocarcinoma and colorectal cancer metastases in the lung. Anticancer Res. 34:2297–2302. 2014.PubMed/NCBI | |
Ma J, Mannoor K, Gao L, Tan A, Guarnera MA, Zhan M, Shetty A, Stass SA, Xing L and Jiang F: Characterization of microRNA transcriptome in lung cancer by next-generation deep sequencing. Mol Oncol. 8:1208–1219. 2014. View Article : Google Scholar : PubMed/NCBI | |
Boyerinas B, Park SM, Hau A, Murmann AE and Peter ME: The role of let-7 in cell differentiation and cancer. Endocr Relat Cancer. 17:F19–F36. 2010. View Article : Google Scholar | |
Ali S, Saleh H, Sethi S, Sarkar FH and Philip PA: MicroRNA profiling of diagnostic needle aspirates from patients with pancreatic cancer. Br J Cancer. 107:1354–1360. 2012. View Article : Google Scholar : PubMed/NCBI | |
Nagadia R, Pandit P, Coman WB, Cooper-White J and Punyadeera C: miRNAs in head and neck cancer revisited. Cell Oncol (Dordr). 36:1–7. 2013. View Article : Google Scholar | |
Liu X, Chen Z, Yu J, Xia J and Zhou X: MicroRNA profiling and head and neck cancer. Comp Funct Genomics. 2009:8375142009. View Article : Google Scholar : | |
Kolenda T, Przybyła W, Teresiak A, Mackiewicz A and Lamperska KM: The mystery of let-7d - a small RNA with great power. Contemp Oncol (Pozn). 18:293–301. 2014. | |
Peng XC, Gong FM, Zhao YW, Zhou LX, Xie YW, Liao HL, Lin HJ, Li ZY, Tang MH and Tong AP: Comparative proteomic approach identifies PKM2 and cofilin-1 as potential diagnostic, prognostic and therapeutic targets for pulmonary adenocarcinoma. PLoS One. 6:e273092011. View Article : Google Scholar : PubMed/NCBI | |
Shi HS, Li D, Zhang J, Wang YS, Yang L, Zhang HL, Wang XH, Mu B, Wang W, Ma Y, et al: Silencing of pkm2 increases the efficacy of docetaxel in human lung cancer xenografts in mice. Cancer Sci. 101:1447–1453. 2010. View Article : Google Scholar : PubMed/NCBI | |
Anumanthan G, Halder SK, Friedman DB and Datta PK: Oncogenic serine-threonine kinase receptor-associated protein modulates the function of Ewing sarcoma protein through a novel mechanism. Cancer Res. 66:10824–10832. 2006. View Article : Google Scholar : PubMed/NCBI | |
Carter JJ, Feingold DL, Oh A, Kirman I, Wildbrett P, Stapleton G, Asi Z, Fowler R, Bhagat G, Huang EH, et al: Perioperative immunomodulation with Flt3 kinase ligand or a whole tumor cell vaccine is associated with a reduction in lung metastasis formation after laparotomy in mice. Surg Innov. 13:41–47. 2006. View Article : Google Scholar : PubMed/NCBI | |
Chakravarty PK, Alfieri A, Thomas EK, Beri V, Tanaka KE, Vikram B and Guha C: Flt3-ligand administration after radiation therapy prolongs survival in a murine model of metastatic lung cancer. Cancer Res. 59:6028–6032. 1999. | |
Sasi W, Sharma AK and Mokbel K: The role of suppressors of cytokine signalling in human neoplasms. Mol Biol Int. 2014:6307972014. View Article : Google Scholar : PubMed/NCBI | |
Ge D, Gao AC, Zhang Q, Liu S, Xue Y and You Z: LNCaP prostate cancer cells with autocrine interleukin-6 expression are resistant to IL-6-induced neuroendocrine differentiation due to increased expression of suppressors of cytokine signaling. Prostate. 72:1306–1316. 2012. View Article : Google Scholar : PubMed/NCBI | |
Karakaidos P, Taraviras S, Vassiliou LV, Zacharatos P, Kastrinakis NG, Kougiou D, Kouloukoussa M, Nishitani H, Papavassiliou AG, Lygerou Z, et al: Overexpression of the replication licensing regulators hCdt1 and hCdc6 characterizes a subset of non-small-cell lung carcinomas: Synergistic effect with mutant p53 on tumor growth and chromosomal instability - evidence of E2F-1 transcriptional control over hCdt1. Am J Pathol. 165:1351–1365. 2004. View Article : Google Scholar : PubMed/NCBI | |
Li JN, Feng CJ, Lu YJ, Li HJ, Tu Z, Liao GQ and Liang C: mRNA expression of the DNA replication-initiation proteins in epithelial dysplasia and squamous cell carcinoma of the tongue. BMC Cancer. 8:3952008. View Article : Google Scholar |