吉林大学学报(医学版)  2016, Vol. 42 Issue (01): 85-88

扩展功能

文章信息

于飞飞, 田雨溪, 刘晓冬, 马淑梅
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
宫颈鳞状细胞癌组织中miRNA差异性表达及其靶基因作为诊断标志物的意义
于飞飞, 田雨溪, 刘晓冬, 马淑梅     
吉林大学公共卫生学院卫生部放射生物学重点实验室, 吉林长春 130021
摘要: 目的: 采用生物信息学方法研究宫颈鳞状细胞癌(CESC)组织中差异性表达的miRNA并预测其靶基因,以期找到影响CESC发生发展及可用于肿瘤标志物的miRNA。 方法: 癌症基因组图谱 (TCGA)数据库中下载3例CESC组织及3例配对正常组织,用统计学方法对差异性表达的基因及miRNA进行筛选,通过靶基因预测网站对差异性表达的miRNA进行靶基因预测,运用KEGG数据库分析靶基因参与的肿瘤相关信号通路。 结果: 与同源正常组织比较,CESC组织有18个miRNA和1180个基因有差异性表达(P<0.05),其中有15个miRNA和411个基因上调,3个miRNA和770个基因下调,采用靶基因预测网站有7个miRNA与预测到的靶基因呈负向调控关系,6个上调miRNA的靶基因下调,1个下调miRNA的靶基因上调,靶基因集中在Wnt、MAPK、P53和cAMP等肿瘤相关信号通路。 结论: 差异性表达的miRNA包括hsa-miR-27a、hsa-miR-148b、hsa-miR-185、hsa-miR-200a、hsa-miR-200b、hsa-miR-221和hsa-mir-133b,差异性表达的miRNA及其靶基因在CESC的发生发展过程中起重要作用,有可能成为诊断CESC的肿瘤标志物。
关键词: 宫颈鳞状细胞癌    微小RNA    肿瘤标志物    
Differential expression of miRNA in cervical squamous cell carcinoma tissue and significance of its target genes as diagnostic markers
YU Feifei, TIAN Yuxi, LIU Xiaodong, MA Shumei    
Key Laboratory of Radiobiology, Ministry of Health, School of Public Health, Jilin University, Changchun 130021, China
Abstract: Objective: To study the differential expression of miRNA in cervical squamous cell carcinoma(CESC)tissue and predict its target genes using bioinformatics methods,and to find the miRNA as tumor markers in the pathological process of CESC.Methods: The CESC samples and their matching normal CESC tissue in TCGA database were download,and the differential expression of genes and miRNAs were screened with statistical method,and the target gene prediction was performed for the differently expressed miRNA with target gene prediction websites.KEGG database was used to analyze the target genes involved in tumor related signaling pathway.Results: Compared with their matching normal CESC tissue,21 miRNA and 1 180 genes were found to be differentially expressed(P<0.05),18 miRNA and 411 genes were up-regulated,and 3 miRNA and 770 genes were down-regulated.7 miRNAs had negative correlationship with the predicted target genes;6 up-regulated miRNA were accompanied by their target genes down-regulation,and down-regulated 1 miRNA were accompanied by their target genes up-regulation;their target genes were involved in Wnt,MAPK,P53 and cAMP tumor related signaling pathways.Conclusion: The differentially expressed miRNAs including hsa-miR-27a,hsa-miR-148b,hsa-miR-185,hsa-miR-200a,hsa-miR-200b,hsa-miR-221 and hsa-mir-133b and their target genes play important roles in the development of CESC,and they are possible to become the tumor markers in the diagnosis of CESC.
Key words: cervical squamous cell carcinoma    microRNA    tumor marker    

生物信息学是将分子生物学与信息处理技术结合,以计算机为工具对生物信息进行储存、检索和分析的交叉学科,其目的是利用各种数据库,分析整理其数据的意义而揭示大量复杂的生物数据所赋有的生物学奥秘。癌症基因组图谱(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

表 1 CESC组织中差异性表达的miRNA Tab. 1 Differential expressions of miRNAs in CESC tissue
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
2.2 差异性表达miRNA作为肿瘤标志物能力判断

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)。

表 2 肿瘤相关信号通路中异常表达的miRNAs及其靶基因 Tab. 2 Abnormal expressions of miRNAs and target genes in tumor-related signal pathway
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
3 讨论

宫颈癌是妇女生殖系统中最常见的恶性肿瘤,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.