扩展功能
文章信息
- 于飞飞, 田雨溪, 刘晓冬, 马淑梅
- YU Feifei, TIAN Yuxi, LIU Xiaodong, MA Shumei
- 宫颈鳞状细胞癌组织中miRNA差异性表达及其靶基因作为诊断标志物的意义
- Differential expression of miRNA in cervical squamous cell carcinoma tissue and significance of its target genes as diagnostic markers
- 吉林大学学报(医学版), 2016, 42(01): 85-88
- Journal of Jilin University (Medicine Edition), 2016, 42(01): 85-88
- 10.13481/j.1671-587x.20160117
-
文章历史
- 收稿日期: 2015-06-21
生物信息学是将分子生物学与信息处理技术结合,以计算机为工具对生物信息进行储存、检索和分析的交叉学科,其目的是利用各种数据库,分析整理其数据的意义而揭示大量复杂的生物数据所赋有的生物学奥秘。癌症基因组图谱(the cancer genome atlas,TCGA)数据库是由美国国家癌症研究所(NCI)及国家人类基因组研究所(NHGRI)联合建立,其中包含丰富的数据类型和肿瘤类型,不需任何费用即可获得大量数据,其次在TCGA下载的数据已经经过了前期繁琐的标准化处理,节省了人工处理时间和资源。 Grieb等[1]使用TCGA数据库成功筛选出治疗三阴性乳腺癌的新靶分子双微体2癌基因结合蛋白(MDM2 binding protein,MTBP)。Yang等[2] 通过分析TCGA数据库为卵巢癌的治疗提供新的思路。本研究利用TCGA数据库下载和分析宫颈鳞状细胞癌(cervical squamous cell carcinoma,CESC)组织相对于配对正常组织中miRNA和基因的表达是否有显著差异,寻找潜在的可用于诊断CESC的肿瘤标志物。miRNAs是一类长度为18~25个核苷酸的单链非编码RNA,在细胞的生存、增殖和肿瘤形成过程中发挥重要的作用[3, 4]。在动物细胞中成熟的miRNA主要通过与目的基因调控序列的3′端 非翻译区(3′-UTR)结合而使mRNA降解或抑制其转录后翻译过程使蛋白水平下降,从而起到调节基因表达的作用[5]。研究[6]显示:miRNA还可以受到其上游基因、转录因子及表观遗传因子的调控参与肿瘤的发生发展过程。在宫颈癌的研究[7, 8, 9]中已证实:miR-182通过调控DNMT3a促进癌细胞凋亡,miR-21的表达水平可作为肿瘤的预后标志物,miR-133a作为治疗肿瘤的潜在靶基因。本研究通过生物信息学方法预测差异性miRNA及其靶基因,正确认识miRNA与靶基因的相互作用,为进一步研究CESC诊断标志物及其发生机制奠定基础。
1 资料与方法 1.1 资料来源使用TCGA数据库获取CESC肿瘤和正常配对组织的基因数据及miRNA数据。
1.2 下载肿瘤正常组织配对样本数据库进入TCGA数据门户主页,选择对应的肿瘤进入下载页面。在下载页面中选miRNAseq和RNAseqV2,选中肿瘤组织和正常组织的配对样本,点击“Build Archive”输入邮箱地址,根据所提供的链接开始下载数据。
1.3 整理肿瘤正常组织配对样本数据库在TCGA下载的数据中,每个样本的基因和miRNA均为有独立编号的文件,肿瘤组织样本和正常组织样本一一对应,将每种肿瘤标准化的基因或miRNA数据汇总到Excel表格里,建立原始的数据库。
1.4 在CESC中筛选差异性表达的基因和miRNA①为提高筛选的准确性和简化筛选过程,首先剔除肿瘤或正常组织中荧光信号值为0的基因或miRNA数据,这些数据是芯片在检测过程中由本底的荧光值决定的;②对肿瘤组织数据和正常组织数据进行统计学分析,P < 0.05为差异有统计学意义,纳入筛选范畴;③肿瘤组织荧光均值/正常组织荧光均值≥2被认为是上调的基因或miRNA,肿瘤组织荧光均值/正常组织荧光均值≤0.05为下调的基因或miRNA,纳入后续的分析。
1.5 受试者工作特征(ROC)曲线为评价本研究方法特异性和敏感性的关系,用Medcalc软件绘制ROC曲线将灵敏度与特异性以图示方法结合。曲线下面积(AUC)越接近于1,说明诊断效果越好。AUC在 0.5~0.7时诊断价值较弱,AUC在0.8~0.9时有一定诊断价值,AUC在0.9以上时有较高诊断价值。AUC=0.5时,说明诊断方法完全不起作用,无诊断价值。
1.6 约登指数(Youden’ s index,YI)YI是灵敏度与特异度之和减去1,表示筛检方法发现真正的患者与非患者的总能力,YI越大说明筛查的效果越好,真实性越大。
1.7 miRNA靶基因预测使用TargetScan、PICTAR、miRWalk、miRanda和DIANAmT 5个常用的miRNA靶基因预测网站对差异表达的miRNA分别进行靶基因预测,5个网站全部预测出的靶基因为有效靶基因。根据miRNA与靶基因呈负调控关系,将上调miRNA的预测靶基因与差异性表达下调基因汇总,将下调miRNA的预测靶基因与差异性表达上调基因汇总,重合基因即为差异miRNA在CESC中的靶基因。
1.8 miRNA靶基因功能分析利用KEGG数据库筛选出差异表达的基因参与的肿瘤相关信号通路,确定影响CESC发生发展的基因。
1.9 统计学分析采用SPSS 13.0软件包进行统计学处理。对CESC组织与配对正常组织基因和miRNA的差异进行正态性检验,符合正态性分布的样本组间比较采用配对样本t检验,否则采用配对秩和检验。以P < 0.05为差异有统计学意义。
2 结果 2.1 CESC组织中差异性表达基因和miRNA与同源正常组织比较,在CESC组织中共有18个miRNAs(P < 0.05)和1 180个基因(P < 0.05)差异性表达,其中miRNA有15个上调,3个下调;基因有411个上调,770个下调。见表1。
miRNA | Tumor | Normal | Fold | P |
hsa-mir-200a | 956 | 72.34 | 13.22 | <0.05 |
hsa-mir-942 | 9.68 | 1.20 | 11.40 | <0.05 |
hsa-mir-200b | 823 | 109 | 8.08 | <0.05 |
hsa-mir-942 | 3.45 | 0.30 | 11.40 | <0.05 |
hsa-mir-21 | 381 542 | 59 372 | 6.43 | <0.05 |
hsa-mir-877 | 0.94 | 0.18 | 5.35 | <0.05 |
hsa-mir-766 | 5.57 | 1.58 | 3.52 | <0.05 |
hsa-mir-27a | 2 164 | 692 | 3.13 | <0.05 |
hsa-mir-3615 | 5.31 | 1.45 | 3.65 | <0.05 |
hsa-mir-3610 | 1.99 | 0.60 | 3.33 | <0.05 |
hsa-mir-24-2 | 1 789 | 854 | 2.09 | <0.05 |
hsa-mir-185 | 51.65 | 18.16 | 2.85 | <0.05 |
hsa-mir-221 | 572 | 201 | 2.84 | <0.05 |
hsa-mir-148b | 107 | 38.92 | 2.75 | <0.05 |
hsa-mir-653 | 4.54 | 1.70 | 2.67 | <0.05 |
hsa-mir-143 | 96 795 | 351 649 | 0.28 | <0.05 |
hsa-mir-133b | 2.24 | 11.18 | 0.20 | <0.05 |
hsa-mir-139 | 18.45 | 106 | 0.17 | <0.05 |
ROC曲线结果显示:筛选出的18种差异性表达miRNA的灵敏度、特异度均为100%,AUC均为1,约登指数均为1,提示筛选出来的miRNA作为肿瘤标志物的效力强。
2.3 差异性表达miRNA靶基因预测通过TargetScan、PICTAR、miRWalk、miRanda和DIANAmT 5个靶基因网站的预测,预测出靶基因的miRNA共9个,其中7个上调,2个下调;根据靶基因预测网站共筛选出靶基因97个,其中91个下调,6个上调。
2.4 CESC中靶基因的相关信号通路分析采用KEGG信号通路数据库对miRNA调控基因参与信号通路研究显示:由miRNA调控的基因可能广泛参与肿瘤相关信号通路,包括MAPK、Wnt、P53、cAMP及AMPK信号通路。在肿瘤相关通路中有7个由miRNA调控的基因参与,6个基因下调,1个基因上调(表2)。
miRNA | Up/Down | Target gene |
hsa-miR-27a | Up | FOXO1,PRICKLE2 |
hsa-miR-148b | Up | PRICKLE2 |
hsa-miR-185 | Up | SOX13 |
hsa-miR-200a | Up | THBS1 |
hsa-miR-200b | Up | EDNRA |
hsa-miR-221 | Up | MAPK10 |
hsa-mir-133b | Down | TCF7 |
宫颈癌是妇女生殖系统中最常见的恶性肿瘤,CESC占宫颈癌的80%~85%[10],高发年龄为50~55岁。近年国外研究[11]报道:在普通人群中宫颈癌的发病率有所下降,但在发展中国家,尤其在我国,由于人乳头瘤病毒(HPV)感染增多、性传播疾病和筛查程序的混乱等原因,宫颈癌的发病率和死亡率反而增加并呈年轻化倾向,其中99.7%的宫颈癌与HPV感染相关。虽然子宫颈细胞学的早期检查联合手术和放化疗的综合治疗方案已明显改善宫颈癌患者的生活质量,但肿瘤的侵袭和转移却难以根治,这也是绝大多数宫颈癌患者死亡的主要因素。因此,病因预防、早期发现和控制肿瘤的转移对提高宫颈癌患者的生存率,降低死亡率有重要意义。
目前miRNA已是基因表达调控网络中重要的组成部分,可能参与调控多条信号传导通路。不同肿瘤组织中含有不同的miRNA表达谱,对不同肿瘤组织特定miRNA表达水平的研究不仅有利于阐明肿瘤的发生发展机制,更为肿瘤的早期诊断和治疗提供了理论依据。近年来研究报道宫颈癌与miRNA的异常表达有密切关系。Shi等[12]研究发现:miR-218多态性(rs11134527)与中国人子宫颈癌的发病有关。Wilting等[13]运用直接定量甲基化特异性PCR发现宫颈癌组miRNA甲基化程度与对照组miRNA甲基化程度存在明显差异,证实宫颈癌的发生发展与miRNA的甲基化有关。HPV感染者中miRNA-218下调参与宫颈癌的发病机制也已经被证实[14, 15]。本研究通过比较CESC与宫颈正常组织筛选出差异表达的miRNA,并通过靶基因预测网站找出CESC中差异表达的miRNA可能的靶基因,并且将靶基因进行肿瘤相关信号通路分析,找出参与肿瘤发生发展的基因,进而推测miRNA在CESC中的作用。
本研究通过比较CESC组织与配对正常组织,共发现有18个miRNA和1180个基因差异表达,通过对差异表达的miRNA行ROC曲线分析发现其特异度和灵敏度较高,提示筛选出来的差异表达的miRNA可作为潜在的肿瘤标志物。运用生物信息学靶基因预测网站发现:hsa-miR-27a、hsa-miR-148b、hsa-miR-185、hsa-miR-200a、hsa-miR-200b、hsa-miR-221和hsa-mir-133b可以预测到靶基因;通过KEGG数据库对信号通路分析显示:靶基因参与肿瘤相关信号通路包括Wnt信号通路、MAPK信号通路、P53信号通路和cAMP信号通路,提示差异表达的miRNA及其靶基因可能在CESC的发生发展中起重要作用。由于本研究的样本数量有限,hsa-miR-27a、hsa-miR-148b、hsa-miR-185、hsa-miR-200a、hsa-miR-200b、hsa-miR-221和hsa-mir-133b做为CESC诊断标志物的准确性和特异性仍需进一步验证。
[1] | Grieb BC,Chen X,Eischen CM. MTBP is overexpressed in triple-negative breast cancer and contributes to its growth and survival[J].Mol Cancer Res,2014,12(9):1216-1224. |
[2] | Yang D,Sun Y,Hu L,et al.Integrated analyses identify a master microRNA regulatory network for the mesenchymal subtype in serous ovarian cancer[J].Cancer Cell,2013,23(2):186-199. |
[3] | Calin GA,Croce CM.MicroRNA-cancer connection:the beginning of a new tale[J].Cancer Res,2006,66(15):7390-7394. |
[4] | Michael MZ,O'Connor SM,van Holst Pellekaan NG,et al. Reduced accumulation of specific microRNAs in colorectal neoplasia[J].Mol Cancer Res,2003,1(12):882-891. |
[5] | Wilting SM,Snijders PJ,Verlaat W,et al.Altered microRNA expression associated with chromosomal changes contributes to cervical carcinogenesis[J].Oncogene,2013,32(1):106-116. |
[6] | Rao Q,Shen Q,Zhou H,et al.Aberrant microRNA expression in human cervical carcinomas[J].Med Oncol,2012,29(2):1242-1248. |
[7] | Sun J,Ji J,Huo G,et al.miR-182 induces cervical cancer cell apoptosis through inhibiting the expression of DNMT3a[J].Int J Clin Exp Pathol,2015,8(5):4755-4763. |
[8] | Han Y,Xu GX,Lu H,et al.Dysregulation of miRNA-21 and their potential as biomarkers for the diagnosis of cervical cancer[J].Int J Clin Exp Pathol,2015,8(6):7131-7139. |
[9] | Song X,Shi B,Huang K,et al.miR-133a inhibits cervical cancer growth by targeting EGFR[J].Oncol Rep,2015,34(3):1573-1580. |
[10] | 丰有吉,沈 铿.妇产科学[M].2版.北京: 人民卫生出版社,2010. |
[11] | Saavedra KP,Brebi PM,Roa JC.Epigenetic alterations in preneoplastic and neoplastic lesions of the cervix[J].Clin Epigenetics,2012,4(1):13. |
[12] | Shi TY,Chen XJ,Zhu ML,et al.A pri-miR-218 variant and risk of cervical carcinoma in Chinese women[J].BMC Cancer,2013,13:19. |
[13] | Wilting SM,Verlaat W,Jaspers A,et al.Methylation-mediated transcriptional repression of microRNAs during cervical carcinogenesis[J].Epigenetics,2013,8(2):220-228. |
[14] | Li Y,Liu J,Yuan C,et al.High-risk human papillomavirus reduces the expression of microRNA-218 in women with cervical intraepithelial neoplasia[J].J Int Med Res,2010,38(5):1730-1736. |
[15] | 莫祥兰,韦荣干,卢军营,等.Sirt1和ProEx C在宫颈高度鳞状上皮内病变细胞学诊断中的价值[J].中国实验诊断学,2015,19(11):1886-1889. |