唾液腺腺样囊性癌中相关微小RNAs及其靶基因的生物信息学分析

刘发煇; 侯婉云; 梁家东; 徐晶晶; 罗春英

doi:10.16098/j.issn.0529-1356.2021.04.015

PDF(11507 KB)

解剖学报 ›› 2021, Vol. 52 ›› Issue (4) : 601-608. DOI: 10.16098/j.issn.0529-1356.2021.04.015

肿瘤生物学

唾液腺腺样囊性癌中相关微小RNAs及其靶基因的生物信息学分析

刘发煇^1,2 侯婉云^1,2 梁家东² 徐晶晶1 罗春英^1,2*

作者信息 +

Bioinformatics analysis of the microRNAs and target genes of microRNAs in salivary adenoid cystic carcinoma

LIU Fa-hui^1,2 HOU Wan-yun^1,2 LIANG Jia-dong² XU Jin-jing¹ LUO Chun-ying^{1, 2*}

Author information +

文章历史 +

摘要

目的探讨唾液腺腺样囊性癌(SACC)潜在的微小RNAs(miRNAs）分子标志物，构建miRNA-mRNA调控网络，并阐明其潜在的分子机制。方法从Gene Expression Omnibus(GEO)数据库下载2个SACC的微阵列芯片数据，通过R语言进行分析差异的miRNAs与mRNA。应用FunRich 3.1.3软件对差异miRNA进行转录因子富集分析以及预测差异miRNAs的靶基因。对SACC中差异miRNAs的靶基因进行基因本体论（GO）富集分析与京都基因与基因组百科全书（KEGG）通路富集分析以及蛋白互作分析。使用Cytoscape 3.7.0构建miRNA-mRNA调控网络。结果 1. 共筛选出144个差异表达的miRNAs (DEMs)和1216个差异表达的mRNAs(DEGs)；2. KEGG信号通路富集分析发现，靶基因主要参与Rap1信号通路、丝裂原活化蛋白激酶（MAPK）信号通路、磷脂酰肌醇-3激酶/蛋白激酶B（PI3K/Akt）信号通路以及肌动蛋白细胞骨架的调节；3. STRING蛋白相互作用分析发现，ACSL1、SCD、MGLL、FABP4在蛋白相互作用网络中处于核心地位。结论我们筛选出SACC与相对正常组织间差异明显的miRNA和mRNA，并分析发现了这些差异分子主要参与的信号通路和功能。

Abstract

Objective To identify potential microRNAs(miRNAs) in salivary adenoid cystic carcinoma and to construct a miRNA-mRNA regulatory network to better understand its potential molecular mechanisms. Methods Two microarray datasets of SACC were downloaded from the database Gene Expression Omnibus(GEO), and the differentially expressed miRNAs and mRNA were analyzed by the R language. FunRich 3.1.3 software was used to enrich and analyze the transcription factors of differential miRNAs and to predict the target genes of differentially expressed miRNAs. The target genes of differential miRNAs in SACC were utilized to perform Gene Onotology(GO) and Kyoto Encyclopedia of Gene Genomes(KEGG) pathway enrichment analyses, and protein-protein interaction. The miRNA-mRNA regulatory network was constructed in Cytoscape 3.7.0. Results A total of 144 differentially expressed miRNA (DEMs) and 1216 differentially expressed mRNA (DEGs) were screened. The enrichment analysis of KEGG signaling pathway revealed that target genes were mainly involved in the regulation of Rap1 signaling pathway, mitogen active protein kinase(MAPK) signaling pathway, phosphatidylinositol 3-kinase/protein kinase B(PI3K/Akt) signaling pathway, and regulation of actin cytoskeleton. STRING protein interaction analysis shows that ACSL1, SCD, MGLL, FABP4 may be the key proteins in the protein interaction network. Conclusion Differentially expressed miRNA and mRNA between SACC tissues and normal tissues were screened out and the signaling pathways and functions of these differential molecules were found in our research.

导出引用

刘发煇侯婉云梁家东徐晶晶罗春英. 唾液腺腺样囊性癌中相关微小RNAs及其靶基因的生物信息学分析[J]. 解剖学报. 2021, 52(4): 601-608 https://doi.org/10.16098/j.issn.0529-1356.2021.04.015

LIU Fa-hui HOU Wan-yun LIANG Jia-dong XU Jin-jing LUO Chun-ying. Bioinformatics analysis of the microRNAs and target genes of microRNAs in salivary adenoid cystic carcinoma[J]. Acta Anatomica Sinica. 2021, 52(4): 601-608 https://doi.org/10.16098/j.issn.0529-1356.2021.04.015

中图分类号： R739

参考文献

［1］ Fan J, Yang J, Qiao W, et al. LAPTM4B-35 expression is associated with pathological grades and clinical stages in salivary adenoid cystic carcinoma［J］. Oncol Lett, 2020, 19(1):317-322.

［2］ Jiang Y, Tang Y, Wang S, et al. PRRX1-induced epithelial-to-mesenchymal transition in salivary adenoid cystic carcinoma activates the metabolic reprogramming of free fatty acids to promote invasion and metastasis［J］. Cell Prolif, 2020, 53(1):e12705.

［3］ Qiao Z, Zou Y, Zhao H. MicroRNA-140-5p inhibits salivary adenoid cystic carcinoma progression and metastasis via targeting survivin［J］. Cancer Cell Int, 2019, 19(1): 301-312.

［4］ Xia Y, Li B, Gao N, et al. Expression of tumor-associated calcium signal transducer 2 in patients with salivary adenoid cystic carcinoma: Correlation with clinicopathological features and prognosis［J］. Oncol Lett, 2014, 8(4): 1670-1674.

［5］ Hu XM, Lin T, Huang XY, et al. ID1 contributes to cell growth invasion and migration in salivary adenoid cystic carcinoma［J］. Mol Med Rep, 2017, 16(6): 8907-8915.

［6］ Liang Y, Ye J, Jiao J, et al. Down-regulation of miR-125a-5p is associated with salivary adenoid cystic carcinoma progression via targeting p38/JNK/ERK signal pathway［J］. Am J Transl Res, 2017, 9(3): 1101-1113.

［7］ Qin SL, Zhou Y, Hu W, et al. Bioinformatics analysis and functional analysis of bladder cancer gene microarray［J］. Acta Anatomica Sinica, 2019, 50(3): 310-316. (in Chinese)

秦松林,周洋,胡威,等.膀胱癌基因芯片生物信息分析及功能探讨［J］.解剖学报,2019,50(3):310-316.

［8］ Liu ZhJ, Dou LT, Zhuang Zh, et al. Screening and survival analysis of differentially expressed genes in papillary thyroid ［J］. Acta Anatomica Sinica, 2020, 51(1): 51-57. (in Chinese)

刘中娇, 窦龙涛, 庄振, 等. 乳头状甲状腺癌差异表达基因的筛选及生存分析［J］. 解剖学报, 2020, 51(1): 51-57.

［9］ Kiss O, Tkés AM, Vranic S, et al. Expression of miRNAs in adenoid cystic carcinomas of the breast and salivary glands［J］. Virchows Arch, 2015, 467(5): 551-562.

［10］ Qiao Z, Zou Y, Zhao H. MicroRNA-140-5p inhibits salivary adenoid cystic carcinoma progression and metastasis via targeting survivin［J］. Cancer Cell Int, 2019, 19(1): 301-312.

［11］ Brown AL, Al-Samadi A, Sperandio M, et al. MiR-455-3p, miR-150 and miR-375 are aberrantly expressed in salivary gland adenoid cystic carcinoma and polymorphous adenocarcinoma［J］. J Oral Pathol Med, 2019, 48(9): 840-845.

［12］ Li J, Gao H, Meng L, et al. Mithramycin inhibits epithelial-to-mesenchymal transition and invasion by downregulating SP1 and SNAI1 in salivary adenoid cystic carcinoma［J］. Tumor Biol, 2017, 39(6): 1010428317708697.

［13］ Tian H, Ge C, Li H, et al. Ribonucleotide reductase M2B inhibits cell migration and spreading by early growth response protein 1-mediated phosphatase and tensin homolog/Akt1 pathway in hepatocellular carcinoma［J］. Hepatology, 2014, 59(4): 1459-1470.

［14］ Wang W, Xiong Y, Ding X, et al. Cathepsin L activated by mutant p53 and Egr-1 promotes ionizing radiation-induced EMT in human NSCLC［J］. J Exp Clin Cancer Res, 2019, 38(1): 61.

［15］ Wu Y, Shi W, Tang T, et al. miR-29a contributes to breast cancer cells epithelial-mesenchymal transition, migration, and invasion via down-regulating histone H4K20 trimethylation through directly targeting SUV420H2［J］. Cell Death Dis, 2019, 10(3): 1-16.

［16］ Zhong G, Lou W, Yao M, et al. Identification of novel mRNA-miRNA-lncRNA competing endogenous RNA network associated with prognosis of breast cancer［J］. Epigenomics, 2019, 11(13): 1501-1518.

［17］ Ding LC, Huang XY, Zheng FF, et al. FZD2 inhibits the cell growth and migration of salivary adenoid cystic carcinomas［J］. Oncol Rep, 2016, 35(2): 1006-1012.

［18］ Xie S, Yu X, Li Y, et al. Upregulation of lncRNA ADAMTS9-AS2 promotes salivary adenoid cystic carcinoma metastasis via PI3K/Akt and MEK/Erk signaling［J］. Mol Ther, 2018, 26(12): 2766-2778.

［19］ He Q, Zhou X, Li S, et al. MicroRNA-181a suppresses salivary adenoid cystic carcinoma metastasis by targeting MAPK-Snai2 pathway［J］. Biochim Biophys Acta, 2013, 1830(11): 5258-5266.