Identification of lymph node metastasis-related microRNAs in lung adenocarcinoma and analysis of the underlying mechanisms using a bioinformatics approach
Restricted accessResearch articleFirst published online April, 2017
Identification of lymph node metastasis-related microRNAs in lung adenocarcinoma and analysis of the underlying mechanisms using a bioinformatics approach
This study aimed to screen lymphatic metastasis-related microRNAs (miRNAs) in lung adenocarcinoma and explore their underlying mechanisms using bioinformatics. The miRNA expression in primary lung adenocarcinoma, matched adjacent non-tumorigenic and lymph node metastasis tissues of patients were profiled via microarray. The screened metastasis-related miRNAs were then validated using quantitative real-time PCR in a second cohort of lung adenocarcinoma patients with lymphatic metastasis. Significance was determined using a paired t-test. Target genes of the metastasis-related miRNAs were predicted using TargetScan, and transcription factors (TFs) were predicted based on the TRANSFAC and ENCODE databases. Furthermore, the related long non-coding RNAs (lncRNAs) were screened with starBase v2.0. The miRNA-TF-mRNA and lncRNA-miRNA-mRNA networks were constructed to determine the key interactions associated with lung adenocarcinoma metastasis. According to the miRNA microarray results, there were 10 miRNAs that were differentially expressed in metastatic tissues compared with primary tumor and adjacent non-tumorigenic tissues. Among them were increased levels of miR-146a-5p, miR-342-3p, and miR-150-5p, which were validated in the second cohort. Based on the miRNA-TF-mRNA network, vascular endothelial growth factor A and transcription factors (TFs) including TP53, SMAD4, and EP300 were recognized as critical targets of the three miRNAs. Interactions involving SNHG16–miR-146a-5p–SMAD4 and RP6-24A23.7–miR-342-3p/miR-150-5p–EP300 were highlighted according to the lncRNA-miRNA-mRNA network. miR-146a-5p, miR-342-3p, and miR-150-5p are lymphatic metastasis-related miRNAs in lung adenocarcinoma. Bioinformatics analyses demonstrated that SNHG16 might inhibit the interaction between miR-146a-5p and SMAD4, while RP6-24A23.7 might weaken miR-342-3p–EP300 and miR-150-5p–EP300 interactions in metastasis.
KamiyaKHayashiYDouguchiJHashiguchiAYamadaTIzumiYWatanabeMKawamuraMHorinouchiHShimadaNKobayashiKSakamotoM. Histopathological features and prognostic significance of the micropapillary pattern in lung adenocarcinoma. Mod Pathol2008; 21: 992–1001.
2.
YouldenDRCrambSMBaadePD. The international epidemiology of lung cancer: geographical distribution and secular trends. J Thorac Oncol2008; 3: 819–31.
3.
XieLYangZLiGShenLXiangXLiuXXuDXuLChenYTianZChenX. Genome-wide identification of bone metastasis-related microRNAs in lung adenocarcinoma by high-throughput sequencing. PloS One2013; 8: e61212–e61212.
4.
GaikwadASouzaCAInacioJRGuptaASekhonHSSeelyJMDennieCGomesMM. Aerogenous metastases: a potential game changer in the diagnosis and management of primary lung adenocarcinoma. AJR Am J Roentgenol2014; 203: W570–82.
5.
WeiLZhaoSChenMChenZChenX. [Clinical research on surgical treatment for non-small cell lung cancer of diameter less than 2 cm]. Zhonghua yi xue za zhi2013; 93: 2502–4.
6.
LiuSLiYQiWZhaoYHuangAShengWLeiBLinPZhuHLiWShenH. Expression of Tiam1 predicts lymph node metastasis and poor survival of lung adenocarcinoma patients. Diagn Pathol2014; 9: 69–69.
7.
ChaoLYi-ShengHYuCLi-XuYXin-LanLDong-LanLJieCYi-LonWHuiLY. Relevance of EGFR mutation with micropapillary pattern according to the novel IASLC/ATS/ERS lung adenocarcinoma classification and correlation with prognosis in Chinese patients. Lung Cancer2014; 86: 164–9.
8.
BernardsRWeinbergRA. A progression puzzle. Nature2002; 418: 823–823.
9.
NguyenDXBosPDMassagueJ. Metastasis: from dissemination to organ-specific colonization. Nat Rev Cancer2009; 9: 274–84.
10.
WongKYHuangXChimCS. DNA methylation of microRNA genes in multiple myeloma. Carcinogenesis2012; 33: 1629–38.
11.
SongQXuYYangCChenZJiaCChenJZhangYLaiPFanXZhouXLinJLiMMaWLuoSBaiX. miR-483-5p promotes invasion and metastasis of lung adenocarcinoma by targeting RhoGDI1 and ALCAM. Cancer Res2014; 74: 3031–42.
12.
RoybalJDZangYAhnYHYangYGibbonsDLBairdBNAlvarezCThilaganathanNLiuDDSaintignyPHeymachJVCreightonCJKurieJM. miR-200 Inhibits lung adenocarcinoma cell invasion and metastasis by targeting Flt1/VEGFR1. Mol Cancer Res2011; 9: 25–35.
13.
MengWYeZCuiRPerryJDedousi-HuebnerVHuebnerAWangYLiBVoliniaSNakanishiHKimTSuhSSAyersLWRossPCroceCMChakravartiAJinVXLautenschlaegerT. MicroRNA-31 predicts the presence of lymph node metastases and survival in patients with lung adenocarcinoma. Clin Cancer Res2013; 19: 5423–33.
14.
WapinskiOChangHY. Long noncoding RNAs and human disease. Trends Cell Biol2011; 21: 354–61.
15.
ThaiPStattSChenCHLiangECampbellCWuR. Characterization of a novel long noncoding RNA, SCAL1, induced by cigarette smoke and elevated in lung cancer cell lines. Am J Respir Cell Mol Biol2013; 49: 204–11.
16.
SalmenaLPolisenoLTayYKatsLPandolfiPP. A ceRNA hypothesis: the Rosetta Stone of a hidden RNA language?Cell2011; 146: 353–8.
17.
JiPDiederichsSWangWBöingSMetzgerRSchneiderPMTidowNBrandtBBuergerHBulkE. MALAT-1, a novel noncoding RNA, and thymosin β4 predict metastasis and survival in early-stage non-small cell lung cancer. Oncogene2003; 22: 8031–41.
18.
BaskervilleSBartelDP. Microarray profiling of microRNAs reveals frequent coexpression with neighboring miRNAs and host genes. RNA2005; 11: 241–7.
19.
ZhangYSzustakowskiJSchinkeM. Bioinformatics analysis of microarray data. Cardiovasc Genom2009; 573: 259–84.
20.
TravisWDBrambillaENicholsonAGYatabeYAustinJHBeasleyMBChirieacLRDacicSDuhigEFliederDB. The 2015 World Health Organization classification of lung tumors: impact of genetic, clinical and radiologic advances since the 2004 classification. J Thorac Oncol2015; 10: 1243–60.
21.
LinFLiRPanZXZhouBYu deBWangXGMaXSHanJShenMLiuHL. miR-26b promotes granulosa cell apoptosis by targeting ATM during follicular atresia in porcine ovary. PloS One2012; 7: e38640–e38640.
22.
Ma Y, Zhang P, Wang F, Zhang H, Yang J, Peng J, Liu W, Qin H. miR-150 as a potential biomarker associated with prognosis and therapeutic outcome in colorectal cancer. Gut 2012;61:1447–53.
23.
GaoXGulariEZhouX. In situ synthesis of oligonucleotide microarrays. Biopolymers2004; 73: 579–96.
24.
BolstadBMIrizarryRAAstrandMSpeedTP. A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics2003; 19: 185–93.
25.
MarkouASourvinouIVorkasPYousefGLianidouE. Clinical evaluation of microRNA expression profiling in non small cell lung cancer. Lung Cancer2013; 81: 388–96.
26.
KanehisaMGotoS. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res2000; 28: 27–30.
27.
LewisBPBurgeCBBartelDP. Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell2005; 120: 15–20.
28.
Huang daWShermanBTTanQCollinsJRAlvordWGRoayaeiJStephensRBaselerMWLaneHCLempickiRA. The DAVID Gene Functional Classification Tool: a novel biological module-centric algorithm to functionally analyze large gene lists. Genome Biol2007; 8: R183–R183.
29.
JiaAYCastillo-MartinMBonalDMSanchez-CarbayoMSilvaJMCordon-CardoC. MicroRNA-126 inhibits invasion in bladder cancer via regulation of ADAM9. Br J Cancer2014; 110: 2945–54.
30.
MatysVFrickeEGeffersRGosslingEHaubrockMHehlRHornischerKKarasDKelAEKel-MargoulisOVKloosDULandSLewicki-PotapovBMichaelHMunchRReuterIRotertSSaxelHScheerMThieleSWingenderE. TRANSFAC: transcriptional regulation, from patterns to profiles. Nucleic Acids Res2003; 31: 374–8.
31.
FranceschiniASzklarczykDFrankildSKuhnMSimonovicMRothALinJMinguezPBorkPvon MeringCJensenLJ. STRING v9.1: protein–protein interaction networks, with increased coverage and integration. Nucleic Acids Res2012; 41: D808–D15.
32.
HurstDREdmondsMDScottGKBenzCCVaidyaKSWelchDR. Breast cancer metastasis suppressor 1 up-regulates miR-146, which suppresses breast cancer metastasis. Cancer Res2009; 69: 1279–83.
33.
ZhouLZhaoXHanYLuYShangYLiuCLiTJinZFanDWuK. Regulation of UHRF1 by miR-146a/b modulates gastric cancer invasion and metastasis. FASEB J2013; 27: 4929–39.
34.
LauvrakSMuntheEKresseSStratfordENamløsHMeza-ZepedaLMyklebostO. Functional characterisation of osteosarcoma cell lines and identification of mRNAs and miRNAs associated with aggressive cancer phenotypes. Br J Cancer2013; 109: 2228–36.
35.
WangHWuJMengXYingXZuoYLiuRPanZKangTHuangW. MicroRNA-342 inhibits colorectal cancer cell proliferation and invasion by directly targeting DNA methyltransferase 1. Carcinogenesis2011; 32: 1033–42.
36.
LiX-rChuH-jLvTWangLKongS-fDaiS-z. miR-342-3p suppresses proliferation, migration and invasion by targeting FOXM1 in human cervical cancer. FEBS Lett2014; 588: 3298–307.
37.
Xie X, Liu H, Wang M, Ding F, Xiao H, Hu F, Hu R, Mei J. miR-342-3p targets RAP2B to suppress proliferation and invasion of non-small cell lung cancer cells. Tumour Biol 2015;36:5031–8.
38.
ItoMTeshimaKIkedaSKitadateAWatanabeANaraMYamashitaJOhshimaKSawadaKTagawaH. MicroRNA-150 inhibits tumor invasion and metastasis by targeting the chemokine receptor CCR6, in advanced cutaneous T-cell lymphoma. Blood2014; 123: 1499–511.
39.
WangWHChenJZhaoFZhangBRYuHSJinHYDaiJH. MiR-150-5p suppresses colorectal cancer cell migration and invasion through targeting MUC4. Asian Pac J Cancer Prev2014; 15: 6269–73.
40.
SrivastavaSKBhardwajASinghSAroraSWangBGrizzleWESinghAP. MicroRNA-150 directly targets MUC4 and suppresses growth and malignant behavior of pancreatic cancer cells. Carcinogenesis2011; 32: 1832–9.
41.
BaffaRFassanMVoliniaSO'HaraBLiuCGPalazzoJPGardimanMRuggeMGomellaLGCroceCM. MicroRNA expression profiling of human metastatic cancers identifies cancer gene targets. J Pathol2009; 219: 214–21.
42.
JiangHHeCGengSShengHShenXZhangXLiHZhuSChenXYangC. RhoT1 and Smad4 are correlated with lymph node metastasis and overall survival in pancreatic cancer. PloS One2012; 7: e42234–e42234.
43.
NatsugoeSXiangmingCMatsumotoMOkumuraHNakashimaSSakitaHIshigamiSBabaMTakaoSAikouT. Smad4 and transforming growth factor beta1 expression in patients with squamous cell carcinoma of the esophagus. Clin Cancer Res2002; 8: 1838–42.
44.
TanakaTWatanabeTKazamaYTanakaJKanazawaTKazamaSNagawaH. Loss of Smad4 protein expression and 18qLOH as molecular markers indicating lymph node metastasis in colorectal cancer—a study matched for tumor depth and pathology. J Surg Oncol2008; 97: 69–73.
45.
LiuJChoS-NAkkantiBJinNMaoJLongWChenTZhangYTangXWistubII. ErbB2 pathway activation upon Smad4 loss promotes lung tumor growth and metastasis. Cell Rep2015; 10: 1599–613.
46.
BianCLiZXuYWangJXuLShenH. Clinical outcome and expression of mutant P53, P16, and Smad4 in lung adenocarcinoma: a prospective study. World J Surg Oncol2015; 13: 1–1.
47.
XuPLiuJDerynckR. Post-translational regulation of TGF-β receptor and Smad signaling. FEBS Lett2012; 586: 1871–84.
48.
KatsunoYLamouilleSDerynckR. TGF-β signaling and epithelial–mesenchymal transition in cancer progression. Curr Opin Oncol2013; 25: 76–84.
ZhuYYuMLiZKongCBiJLiJGaoZLiZ. ncRAN, a newly identified long noncoding RNA, enhances human bladder tumor growth, invasion, and survival. Urology2011; 77: 510:e1–e5, .
51.
GaytherSABatleySJLingerLBannisterAThorpeKChinSFDaigoYRussellPWilsonASowterHMDelhantyJDPonderBAKouzaridesTCaldasC. Mutations truncating the EP300 acetylase in human cancers. Nat Genet2000; 24: 300–3.
52.
MeesSTMardinWAWendelCBaeumerNWillscherESenningerNSchleicherCColombo-BenkmannMHaierJ. EP300ndel C, regulated metastasis suppressor gene in ductal adenocarcinomas of the pancreas. Int J Cancer2010; 126: 114–24.
53.
NikiTIbaSTokunouMYamadaTMatsunoYHirohashiS. Expression of vascular endothelial growth factors A, B, C, and D and their relationships to lymph node status in lung adenocarcinoma. Clin Cancer Res2000; 6: 2431–9.
54.
ChenC-HLaiJ-MChouT-YChenC-YSuL-JLeeY-CChengT-SHongY-RChouC-KWhang-PengJ. VEGFA upregulates FLJ10540 and modulates migration and invasion of lung cancer via PI3K/AKT pathway. PLoS One2009; 4: e5052–e5052.
55.
QiuTZhouXWangJDuYXuJHuangZZhuWShuYLiuP. MiR-145, miR-133a and miR-133b inhibit proliferation, migration, invasion and cell cycle progression via targeting transcription factor Sp1 in gastric cancer. FEBS Lett2014; 588: 1168–77.
56.
KongL-MLiaoC-GZhangYXuJLiYHuangWZhangYBianHChenZ-N. A regulatory loop involving miR-22, Sp1, and c-Myc modulates CD147 expression in breast cancer invasion and metastasis. Cancer Res2014; 74: 3764–78.
