hsa-miR-34a在肿瘤中的表达及生物信息学分析

doi:10.3969/j.issn.1673-8640.2022.07.012

摘要/Abstract

摘要：

目的分析人微小RNA-34a（hsa-miR-34a）在不同肿瘤中的表达情况,预测其靶基因,并分析其生物学功能、作用机制。方法通过miRBase数据库获得序列,分析hsa-miR-34a在人类基因组中的位置及其物种间的保守性。通过TargetScan、miRWalk、miRDB和mirDIP数据库预测hsa-miR-34a的靶基因并取交集,通过TiGER数据库分析所得靶基因集在人体不同组织中的表达情况。利用PANTHER和DAVID在线工具对hsa-miR-34a的靶基因进行基因本体论（GO）富集分析和REACTOME信号通路富集分析。采用STRING数据库构建hsa-miR-34a的靶基因对应蛋白质的互作网络。结果 hsa-miR-34a在乳腺癌、前列腺癌、膀胱癌等肿瘤中表达显著降低。编码hsa-miR-34a的基因序列在不同物种间具有高度保守性。共获得225个靶基因,采用Target Score功能排序后取前30个靶基因。hsa-miR-34a各靶基因在不同组织间的相对表达量差异均无统计学意义（P>0.05）。GO富集分析结果显示,hsa-miR-34a靶基因集中表达于网格蛋白囊泡和囊泡膜及其物质转运、多细胞生物过程和发育、神经系统发生和发育、解剖结构发育等。REACTOME生物通路富集分析结果显示,hsa-miR-34a靶基因集中富集于B型和G型肉毒杆菌毒素毒性、信号素4D（Sema4D）介导的细胞黏附和迁移抑制、神经递质的合成和释放等生物信号通路。根据蛋白质互作网络对靶基因重新进行排序,排在前5位的蛋白质对应的靶基因分别为SYT1、VAMP2、MET、RRAS、CACNA1E。结论 hsa-miR-34a与多种肿瘤的发生、发展有关,分析靶基因及其生物学功能,可为后续的临床和基础研究提供参考依据。

关键词: 微小RNA-34a, 生物信息学分析, 肿瘤, 靶基因, 功能富集, 信号通路富集

Abstract:

Objective To analyze the expression of homo sapiens microRNA（miR）-34a（hsa-miR-34a） in different tumors,to predict its target genes,and to analyze its biological function and mechanism. Methods Sequence was obtained from miRBase database to analyze the position of hsa-miR-34a in human genome and its conservation among species. TargetScan,miRWalk,miRDB and mirDIP were used to predict the target gene set of hsa-miR-34a and take the intersection,the expression of which in different human tissues was analyzed through TiGER database. Gene Ontology（GO） enrichment analysis and REACTOME signaling pathway enrichment analysis of hsa-miR-34a target genes were performed using PANTHER and DAVID online tools. The interaction network of hsa-miR-34a target genes corresponding to target proteins was constructed by STRING database. Results The expressions of hsa-miR-34a in breast cancer,prostate cancer,bladder cancer and other tumors were decreased. The sequence encoding hsa-miR-34a was highly conserved among different species. A total of 225 target genes were obtained,and the top 30 target genes were selected after Target Score functional sequencing. There was no statistical significance in the relative expression of target genes among different human tissues（P>0.05）. GO enrichment analysis showed that hsa-miR-34a target genes were mainly expressed in clathrin vesicles,vesicle membranes and their material transport,multicellular biological processes and development,neurogenesis and development,anatomical structure development and so on. REACTOME signaling pathway enrichment analysis showed that hsa-miR-34a target genes were concentrated in biological signaling pathways such as botulinum toxin B and G toxicities,cell adhesion and migration inhibition mediated by semaphorin 4D（Sema4D） and the synthesis and release of neurotransmitters. The target genes were reordered according to the protein interaction network. The target genes corresponding to the top 5 proteins were SYT1,VAMP2,MET,RRAS and CACNA1E. Conclusions The hsa-miR-34a is related to the occurrence and development of a variety of tumors. The analysis of its target genes and biological functions can provide a reference for subsequent clinical and basic research.

Key words: MicroRNA-34a, Bioinformatics analysis, Tumor, Target gene, Functional enrichment, Signaling pathway enrichment

中图分类号:

R446.7

文舒展, 付子乐, 陈淑英. hsa-miR-34a在肿瘤中的表达及生物信息学分析[J]. 检验医学, 2022, 37(7): 657-663.

WEN Shuzhan, FU Zile, CHEN Shuying. Expression and bioinformatics analysis of hsa-miR-34a in tumors[J]. Laboratory Medicine, 2022, 37(7): 657-663.

图/表 8

参考文献 27

[1]	RÍOS J M. The big world of the small RNAs[J]. Rev Latinoam Microbiol, 2006, 48(2):73-78.
[2]	MATRANGA C, ZAMORE P D. Small silencing RNAs[J]. Curr Biol, 2007, 17(18):R789-R793.
[3]	DU T, ZAMORE P D. MicroPrimer:the biogenesis and function of microRNA[J]. Development, 2005, 132(21):4645-4652. DOI URL
[4]	HERMEKING H. The miR-34 family in cancer and apoptosis[J]. Cell Death Differ, 2010, 17(2):193-199. DOI URL
[5]	TAZAWA H, TSUCHIYA N, IZUMIYA M, et al. Tumor-suppressive miR-34a induces senescence-like growth arrest through modulation of the E2F pathway in human colon cancer cells[J]. Proc Natl Acad Sci U S A, 2007, 104(39):15472-15477. DOI URL
[6]	WU J, WU G, LV L, et al. MicroRNA-34a inhibits migration and invasion of colon cancer cells via targeting to Fra-1[J]. Carcinogenesis, 2012, 33(3):519-528. DOI URL
[7]	PANG R T, LEUNG C O, YE T M, et al. MicroRNA-34a suppresses invasion through downregulation of Notch1 and Jagged1 in cervical carcinoma and choriocarcinoma cells[J]. Carcinogenesis, 2010, 31(6):1037-1044. DOI URL
[8]	PANG R T, LEUNG C O, LEE C L, et al. MicroRNA-34a is a tumor suppressor in choriocarcinoma via regulation of Delta-like1[J]. BMC Cancer, 2013, 13:25. DOI URL
[9]	王明丽, 李智, 徐晋豫, 等. miR-34a通过调控YY1抑制肾癌细胞的增殖及侵袭[J]. 检验医学, 2012, 27(8):635-640.
[10]	ZHANG C, YAO Z, ZHU M, et al. Inhibitory effects of microRNA-34a on cell migration and invasion of invasive urothelial bladder carcinoma by targeting Notch1[J]. J Huazhong Univ Sci Technolog Med Sci, 2012, 32(3):375-382. DOI URL
[11]	YU G, YAO W, XIAO W, et al. MicroRNA-34a functions as an anti-metastatic microRNA and suppresses angiogenesis in bladder cancer by directly targeting CD44[J]. J Exp Clin Cancer Res, 2014, 33(1):779. DOI URL
[12]	YAN K, GAO J, YANG T, et al. MicroRNA-34a inhibits the proliferation and metastasis of osteosarcoma cells both in vitro and in vivo[J]. PLoS One, 2012, 7(3):e33778. DOI URL
[13]	WU X, ZHONG D, GAO Q, et al. MicroRNA-34a inhibits human osteosarcoma proliferation by downregulating ether à go-go 1 expression[J]. Int J Med Sci, 2013, 10(6):676-682. DOI URL
[14]	LI L, YUAN L, LUO J, et al. MiR-34a inhibits proliferation and migration of breast cancer through down-regulation of Bcl-2 and SIRT1[J]. Clin Exp Med, 2013, 13(2):109-117. DOI URL
[15]	PENG Y, GUO J J, LIU Y M, et al. MicroRNA-34A inhibits the growth,invasion and metastasis of gastric cancer by targeting PDGFR and MET expression[J]. Biosci Rep, 2014, 34(3):e00112. DOI URL
[16]	OKADA N, LIN C P, RIBEIRO M C, et al. A positive feedback between p53 and miR-34 miRNAs mediates tumor suppression[J]. Genes Dev, 2014, 28(5):438-450. DOI URL
[17]	LIANG J, LI Y, DANIELS G, et al. LEF1 targeting EMT in prostate cancer invasion is regulated by miR-34a[J]. Mol Cancer Res, 2015, 13(4):681-688. DOI URL
[18]	ZHANG R, SU J, XUE S L, et al. HPV E6/p53 mediated down-regulation of miR-34a inhibits Warburg effect through targeting LDHA in cervical cancer[J]. Am J Cancer Res, 2016, 6(2):312-320.
[19]	CORTEZ M A, IVAN C, VALDECANAS D, et al. PDL1 regulation by p53 via miR-34[J]. J Natl Cancer Inst, 2015, 108(1):djv303.
[20]	SUN T Y, XIE H J, LI Z, et al. MiR-34a regulates HDAC1 expression to affect the proliferation and apoptosis of hepatocellular carcinoma[J]. Am J Transl Res, 2017, 9(1):103-114.
[21]	SLABÁKOVÁ E, CULIG Z, REMŠÍK J, et al. Alternative mechanisms of miR-34a regulation in cancer[J]. Cell Death Dis, 2017, 8(10):e3100. DOI URL
[22]	SOTILLO E, LAVER T, MELLERT H, et al. Myc overexpression brings out unexpected antiapoptotic effects of miR-34a[J]. Oncogene, 2011, 30(22):2587-2594. DOI URL
[23]	LODYGIN D, TARASOV V, EPANCHINTSEV A, et al. Inactivation of miR-34a by aberrant CpG methylation in multiple types of cancer[J]. Cell Cycle, 2008, 7(16):2591-2600. DOI URL
[24]	CHANG T C, WENTZEL E A, KENT O A, et al. Transactivation of miR-34a by p53 broadly influences gene expression and promotes apoptosis[J]. Mol Cell, 2007, 26(5):745-752. DOI URL
[25]	MANDKE P, WYATT N, FRASER J, et al. MicroRNA-34a modulates MDM4 expression via a target site in the open reading frame[J]. PLoS One, 2012, 7(8):e42034. DOI URL
[26]	KARNI-SCHMIDT O, LOKSHIN M, PRIVES C. The roles of MDM2 and MDMX in cancer[J]. Annu Rev Pathol, 2016, 11:617-644. DOI URL
[27]	ZHANG B, VOGELZANG A, MIYAJIMA M, et al. B cell-derived GABA elicits IL-10+ macrophages to limit anti-tumour immunity[J]. Nature, 2021, 599(7885):471-476. DOI URL

疾病种类	表达情况	检测样本	作用	文献	发表年份
结肠癌	降低	细胞系	调节E2F信号通路,发挥细胞增殖抑制作用	TAZAWA等^[5]	2007年
		组织	p53通过miR-34a依赖性方式调节Fra-1蛋白表达,并最终调节细胞迁移和侵袭	WU等^[6]	2012年
绒毛膜癌	降低	细胞系	抑制Notch1和Jagged1基因的表达	PANG等^[7]	2010年
		细胞系	抑制DLL1基因的表达	PANG等^[8]	2013年
肾癌	降低	组织	调节YY1基因的表达,抑制细胞增殖	王明丽等^[9]	2012年
膀胱癌	降低	组织	拮抗Notch1蛋白并抑制膀胱癌细胞的迁移和侵袭	ZHANG等^[10]	2012年
		组织、细胞系、血液	靶向作用于CD44基因,发挥抑制转移和血管生成作用	YU等^[11]	2014年
骨肉瘤	降低	组织	体内外均可通过下调c-Met^①来抑制骨肉瘤细胞的迁移和侵袭	YAN等^[12]	2012年
		组织和细胞系	可以通过下调Eag1^②表达来抑制骨肉瘤的生长	WU等^[13]	2013年
乳腺癌	降低	细胞系	抑制细胞增殖和侵袭,诱导细胞凋亡	LI等^[14]	2013年
胃癌	降低	组织	通过PI3K^③/Akt^④通路靶向调控PDGFR^⑤和EMT^⑥,抑制胃癌发生	PENG等^[15]	2014年
肺腺癌	降低	组织	抑制p53负调节因子基因MDM4的表达	OKADA等^[16]	2014年
前列腺癌	降低	组织	通过抑制基因LEF1表达来抑制EMT,进而抑制前列腺肿瘤细胞的迁移和侵袭	LIANG等^[17]	2015年
宫颈癌	降低	组织	抑制LDHA^⑦,从而抑制肿瘤的Warburg效应,起到抑制肿瘤生长和迁移的作用	ZHANG等^[18]	2016年
非小细胞肺癌	降低	组织	下调PD-L1^⑧表达,从而抑制肿瘤的免疫逃逸	CORTEZ等^[19]	2015年
肝癌	降低	组织	调控HDAC1基因,抑制肝癌细胞增殖	SUN等^[20]	2017年

疾病种类	表达情况	检测样本	作用	文献	发表年份
结肠癌	降低	细胞系	调节E2F信号通路,发挥细胞增殖抑制作用	TAZAWA等^[5]	2007年
		组织	p53通过miR-34a依赖性方式调节Fra-1蛋白表达,并最终调节细胞迁移和侵袭	WU等^[6]	2012年
绒毛膜癌	降低	细胞系	抑制Notch1和Jagged1基因的表达	PANG等^[7]	2010年
		细胞系	抑制DLL1基因的表达	PANG等^[8]	2013年
肾癌	降低	组织	调节YY1基因的表达,抑制细胞增殖	王明丽等^[9]	2012年
膀胱癌	降低	组织	拮抗Notch1蛋白并抑制膀胱癌细胞的迁移和侵袭	ZHANG等^[10]	2012年
		组织、细胞系、血液	靶向作用于CD44基因,发挥抑制转移和血管生成作用	YU等^[11]	2014年
骨肉瘤	降低	组织	体内外均可通过下调c-Met^①来抑制骨肉瘤细胞的迁移和侵袭	YAN等^[12]	2012年
		组织和细胞系	可以通过下调Eag1^②表达来抑制骨肉瘤的生长	WU等^[13]	2013年
乳腺癌	降低	细胞系	抑制细胞增殖和侵袭,诱导细胞凋亡	LI等^[14]	2013年
胃癌	降低	组织	通过PI3K^③/Akt^④通路靶向调控PDGFR^⑤和EMT^⑥,抑制胃癌发生	PENG等^[15]	2014年
肺腺癌	降低	组织	抑制p53负调节因子基因MDM4的表达	OKADA等^[16]	2014年
前列腺癌	降低	组织	通过抑制基因LEF1表达来抑制EMT,进而抑制前列腺肿瘤细胞的迁移和侵袭	LIANG等^[17]	2015年
宫颈癌	降低	组织	抑制LDHA^⑦,从而抑制肿瘤的Warburg效应,起到抑制肿瘤生长和迁移的作用	ZHANG等^[18]	2016年
非小细胞肺癌	降低	组织	下调PD-L1^⑧表达,从而抑制肿瘤的免疫逃逸	CORTEZ等^[19]	2015年
肝癌	降低	组织	调控HDAC1基因,抑制肝癌细胞增殖	SUN等^[20]	2017年

序列号	物种	名称	保守序列	染色体定位
MI0000268	人	hsa-miR-34a	22-UGGCAGUGU-30	chr1：9151668-9151777 [-]^①
MI0002798	倭黑猩猩	ppa-miR-34a	22-UGGCAGUGU-30	chr1：9141308-9141417 [-]^①
MI0002801	猕猴	mml-miR-34a	22-UGGCAGUGU-30	chr1：8320927-8321036 [-]^①
MI0002799	婆罗洲猩猩	ppy-miR-34a	22-UGGCAGUGU-30	chr1：221288574-221288683 [+]^②
MI0000584	小鼠	mmu-miR-34a	20-UGGCAGUGU-28	chr4：150068454-150068555 [+]^②
MI0001251	原鸡	gga-miR-34a	21-UGGCAGUGU-29	chr21：3262992-3263100 [-]^①
MI0008101	家犬	cfa-miR-34a	1-UGGCAGUGU-9	chr5：62485832-62485897 [-]^①
MI0010019	文昌鱼	bfl-miR-34a	21-UGGCAGUGU-29

序列号	物种	名称	保守序列	染色体定位
MI0000268	人	hsa-miR-34a	22-UGGCAGUGU-30	chr1：9151668-9151777 [-]^①
MI0002798	倭黑猩猩	ppa-miR-34a	22-UGGCAGUGU-30	chr1：9141308-9141417 [-]^①
MI0002801	猕猴	mml-miR-34a	22-UGGCAGUGU-30	chr1：8320927-8321036 [-]^①
MI0002799	婆罗洲猩猩	ppy-miR-34a	22-UGGCAGUGU-30	chr1：221288574-221288683 [+]^②
MI0000584	小鼠	mmu-miR-34a	20-UGGCAGUGU-28	chr4：150068454-150068555 [+]^②
MI0001251	原鸡	gga-miR-34a	21-UGGCAGUGU-29	chr21：3262992-3263100 [-]^①
MI0008101	家犬	cfa-miR-34a	1-UGGCAGUGU-9	chr5：62485832-62485897 [-]^①
MI0010019	文昌鱼	bfl-miR-34a	21-UGGCAGUGU-29

GO序列	细胞组分名称	P值	基因数量	基因名称
60203	网格蛋白修饰的谷氨酸转运囊泡膜	2.98 ×10^-5	2	SYT1、VAMP2
60199	网格蛋白修饰的谷氨酸转运囊泡	2.98×10^-5	2	SYT1、VAMP2
70083	网格蛋白修饰的单胺运输囊泡膜	4.17×10^-5	2	SYT1、VAMP2
70081	网格蛋白修饰的单胺运输囊泡	4.17×10^-5	2	SYT1、VAMP2
61202	网格蛋白修饰的γ-氨基丁酸转运囊泡膜	8.92×10^-5	2	SYT1、VAMP2