白雪峰，劉躍娟，齊 靖，張 建

（哈爾濱醫(yī)科大學大慶校區(qū)醫(yī)學信息學院，黑龍江大慶，163319）

胰腺癌中miRNAs、mRNAs基因表達譜分析

白雪峰，劉躍娟，齊 靖，張 建

（哈爾濱醫(yī)科大學大慶校區(qū)醫(yī)學信息學院，黑龍江大慶，163319）

本文研究miRNAs、mRNAs在胰腺癌中的調(diào)控機制。通過對胰腺癌表達譜進行分析，基于生物信息學方法，建立miRNAs-mRNAs網(wǎng)絡(luò)，通過模塊挖掘獲得5個功能顯著性模塊，網(wǎng)絡(luò)中結(jié)點度較大的HUB節(jié)點miRNAs，這些結(jié)點與胰腺癌的發(fā)生、發(fā)展密切相關(guān)。對每個模塊進行聚類分析，得到度數(shù)小卻聚集的集合，這些基因也和胰腺癌密切相關(guān)。對這些模塊進行功能注釋，多個模塊注釋到代謝相關(guān)的功能上，表明胰腺癌與代謝密切相關(guān)；另外，多個模塊中具有多個相同功能注釋，多個調(diào)控相反的功能在同一模塊當中出現(xiàn)。這表明在組織自身原來就存在著平衡調(diào)控機制。通過本文的研究得到miRNAs與mRNAs協(xié)同作用的一些結(jié)果，為今后更深入的研究胰腺癌提供指導和借鑒。

胰腺癌；miRNAs；mRNAs；協(xié)同

0 引言

胰腺癌（pancreatic cancer，PC）因其起病隱匿、早期診斷難、手術(shù)切除困難、放化療不敏感以及預后差、死亡率高等原因，成為危害最大的惡性腫瘤之一，其中90%以上為胰腺導管腺癌。因此，深入了解胰腺癌發(fā)生發(fā)展的分子機制，發(fā)現(xiàn)新的潛在治療靶點尤為重要［1-2］。

近幾年的研究表明，miRNA的表達水平與人類重大疾病如癌癥等密切相關(guān)。這使得miRNA可以作為新的生物標志物用于癌癥等重大疾病的早期診斷并可作為新的基因藥物作用靶點。

最近，研究發(fā)現(xiàn)miRNA在胰腺癌中存在多個異常表達，有的是在胰腺癌癌前病變中研究中發(fā)現(xiàn)的［7-8］；有的是在胰腺癌組織標本中研究發(fā)現(xiàn)的；還有的是在胰腺癌病人血清中研究中發(fā)現(xiàn)的［9-10］。因此研究者認為這些差異表達的miRNA或許會為研究胰腺癌提供新的線索。

目前，雖然識別出一些與胰腺癌發(fā)生發(fā)展關(guān)系密切的miRNAs和mRNAs，但大部分miRNAs和mRNAs在胰腺癌中的作用還是未知的。本研究從miRNA、mRNA表達譜數(shù)據(jù)，構(gòu)建胰腺癌疾病網(wǎng)絡(luò)，進行失調(diào)網(wǎng)絡(luò)的模塊挖掘，初步探索其作用機制，為揭示胰腺癌發(fā)生發(fā)展過程中miRNA與mRNAs的關(guān)系和調(diào)控作用機制提供新的線索，為胰腺癌的生物治療提供理論依據(jù)。

1 miRNA與動脈粥樣硬化

1.1 數(shù)據(jù)準備

從 NCBI上下載mRNA芯片數(shù)據(jù)（GSE43795）和miRNA芯片數(shù)據(jù)（GSE43796），兩種芯片包括相同的11病人樣本（5 normal/6 PCA），數(shù)據(jù)處理得到30500個mRNA和1205miRNA。差異基因的選取通過SAM算法獲得，通過R代碼進行數(shù)據(jù)處理，選取差異表達的mRNA和miRNA（q-value>0.001和q-value>0.05），得到1585mRNAs和106個miRNAs，這些差異表達的mRNAs和miRNAs作為胰腺癌的潛在致病因素。miRNAs與mRNA的互作對從HMDD數(shù)據(jù)庫下載，與疾病相關(guān)的miRNAs為117，與胰腺癌miRNA芯片取交集共有miRNAs92個，這些數(shù)據(jù)將用做以后分析使用。對獲得的疾病表達譜數(shù)據(jù)為獲得更多關(guān)于胰腺癌生物學信息，通過以下步驟行處理（圖1）。利用差異表達的mRNAs和在HMDD數(shù)據(jù)庫miRNA調(diào)控的mRNAs映射到KEGG通路當中，通過Subpathway-GM進行子通路挖掘，最終形成mRNAs-miRNAs互作網(wǎng)絡(luò)。用Cytoscape插件對網(wǎng)絡(luò)進行功能模塊挖掘，獲得5個功能顯著性模塊，最后通過David進行基因的功能注釋。

1.2 表達譜網(wǎng)絡(luò)分析

利用Subpathwaymine 進行子通路重構(gòu)，合并子通路形成miRNA-mRNA互作網(wǎng)絡(luò)，結(jié)點數(shù)731個，其中包括mRNAs和mRNAs（91/640），3030條邊。點表示基因，邊表示兩個基因間的聯(lián)系，在網(wǎng)絡(luò)中，度越大的基因連接的邊越多，表明在網(wǎng)絡(luò)中可能起到HUB結(jié)點的作用。另外，還有一些結(jié)點雖然度小，但在結(jié)點間傳輸起到重要作用的橋梁作用，這些結(jié)點的關(guān)系，可能作為PCA中的潛在的致病基因。

圖1 數(shù)據(jù)分析流程圖

1.3 網(wǎng)絡(luò)模塊挖掘

利用Cytoscape中的JActivemodule插件進行模塊挖掘。獲得5個功能性顯著模塊（module1～module5）。用ClusterOne進行對模塊進行聚類，尋找顯著性miRNA-mRNA互作對（p-value<0.05）。

1.4 Gene ontology（GO）分析

通過David在線平臺完成對模塊中重要結(jié)點進行功能分析，篩選出具有診斷意義的顯示miRNA-mRNA關(guān)系對。

miRNAs在許多疾病的發(fā)生與發(fā)展中扮演了重要的角色，有研究證明，在疾病中miRNAs對mRNAs的調(diào)控行使著重要的生物學功能。通過以上步驟進一步對胰腺癌中miRNAs對mRNAs的調(diào)控關(guān)系做出探索。

2 結(jié)果

2.1 差異表達miRNAs、mRNAs

對PCA表達譜數(shù)據(jù)進行處理，得到差異表達的mRNAs、miRNAs數(shù)量分別為1585、106個。取表達譜差異的mRNAs、miRNAs前5%做熱圖，結(jié)果如圖2所示。在miRNAs表達譜的熱圖中，module1中顯著表達的有mir-320c；module2包括:mir-222、mir-92b、mir-320c；moudule3則有mir-320c、mir-320b、mir-654、mir-92a、mir-92b；module4含有mir-222、mir-320c、mir-92b；module5未出現(xiàn)上述熱圖中的miRNAs。通過比較可知，表達譜中mir-222在兩個模塊中出現(xiàn)，mir-92b在三個模塊中出現(xiàn)，而mir-320c在module1～module4模塊中全部出現(xiàn)。研究表明［11］mir-222通過靶向p57控制癌細胞的增殖。而mir-92b在胰腺癌的研究中未有文獻報道，mir-92b在膀胱癌中僅對癌細胞轉(zhuǎn)移和浸襲起作用，卻不會促使細胞的增殖［12］。miR-222可能通過MAPK信號轉(zhuǎn)導通路參與胰腺癌細胞侵襲轉(zhuǎn)移的調(diào)控［13］。文獻［14］報道m(xù)ir-320c可能是通過另外一種模式lncRNA-mi RNA-mRNA進行調(diào)控。

圖2 mRNAs和miRNAs表達譜熱圖

在mRNAs表達譜的熱圖中，module2中顯著表達的有YWHAH；module3中含有GLS2、YWHAH；module4含有IFNA16；module5出現(xiàn)熱圖中ACADL、YWHAH。通過比較可知，表達譜中YWHAH在moudule2、moudule3、moudule5中出現(xiàn)。 促血管形成Gremlin 1 通過綁定YWHAH可能發(fā)揮著致癌作用，這一結(jié)果在宮頸癌肺癌、卵巢、腎、乳腺癌、結(jié)腸、胰腺癌中得到證實，Gremlin 1結(jié)合目標蛋白YWHAH可能為人類癌癥的提供新的診斷和治療方案［15］。ACADL顯著表達與胰腺固狀瘤有關(guān)［16］。而GLS2蛋白質(zhì)對腫瘤的生長有抑制作用［17］。

2.2 miRNAs-mRNAs網(wǎng)絡(luò)模塊

癌癥屬于復雜疾病，對建立的疾病網(wǎng)絡(luò)已經(jīng)證明顯現(xiàn)冪律分布［18］，從網(wǎng)絡(luò)分析來看，網(wǎng)絡(luò)中結(jié)點度越大的點，表明可能會影響更多的鄰居結(jié)點，因此，這些點被考慮扮演hubs結(jié)點的角色。圖3顯示miRNAs-mRNAs網(wǎng)絡(luò)結(jié)點的度分布是冪律分布。通過模塊挖掘來研究miRNAs、mRNAs在癌癥當中的共調(diào)控方式、以及基于局部策略來分析miRNAs的在胰腺癌中的功能。

圖3 網(wǎng)絡(luò)度分布

根據(jù)得到的5個功能模塊當中（圖4），大部分miRNAs具有較大的度（表1），如:mir-615、mir-222、mir-93-5p、mir-415-3p、mir-877等。在各個模塊中有些結(jié)點的度較大，mir-93-5p在5個模塊中度分別為:19、19、23、26、20；mir-222在模塊二、四、五的度數(shù)為:15、12、15；mir-24在模塊二、三、四、五分別為:18、19、18、13；mir-98在模塊四、五都為45。mir-615-3p在模塊一、二結(jié)點度為61、和65。有意思的是，研究證明hsa-miR-615-5p其異常表達可能與胰腺癌的發(fā)生、發(fā)展、浸潤、轉(zhuǎn)移有關(guān)，mir-615-5p對胰腺癌有抑制作用［19］，而在前面表達譜數(shù)據(jù)中mir-615-5p是下調(diào)的，mir-615-3p是上調(diào)的；mir-615-5p、mir-615-3p分別通過不同的方式來調(diào)節(jié)mRNAs表達。

表1 模塊中的邊和點

在五個模塊中，對模塊進一步分析，通過clusterone插件對模塊進行挖掘，通過挖掘每個模塊中重要節(jié)點用紅色標出，這表明除了度大的點起到了HUB結(jié)點的作用，在網(wǎng)絡(luò)中還有一些結(jié)點度雖然小，在局部聯(lián)系較為緊密，成為一個團隊發(fā)揮著作用，這些合作的結(jié)點可能在疾病發(fā)生過程中起到作用。在模塊一中RPL32、mir-186-5p、RPL27、RPL29、mir-652-3p聚在一個集合中，在這個小集合中，在狗身上檢測到RPL32在各組織上表達都是穩(wěn)定的［20］，在胰腺癌中mir-186 -5p影響細胞增殖和轉(zhuǎn)移［21］，mir-652與ZEB1結(jié)合抑制胰腺癌腫瘤細胞上皮間質(zhì)（EMT）的酸性微環(huán)境，酸性/ mir-655/ZEB1/EMT可能是腫瘤發(fā)生的新機制［22］。模塊二中YWHAH、FOXO3、CDC25C、YWHAB、YWHAZ構(gòu)成的集合中，YWHAH、FOXO3、CDC25C與胰腺癌密切相關(guān)［15］［23］［24］。從模塊二中還可以看到Y(jié)WHAH和mir-186直接相連，F(xiàn)OXO3和mir-615-3p直接相連，mir-615-3p在模塊中起到了HUB結(jié)點的作用。另外mir-744、RAN、EIF4G1作為一組聯(lián)系緊密的集合被發(fā)現(xiàn)，miR-744 發(fā)現(xiàn)成為新的預后標記物［25］，而RAN通過AR和CXCR4調(diào)控胰腺癌的侵襲和轉(zhuǎn)移［26］。RAN也與mir-615-3p直接相連。模塊三中mir-10a、mir-320c、LDHA、MT-CO3被聚在一起，DNMT3A綁定HK1和LDHA，而這兩個基因作為代謝開關(guān)，動物實驗表明敲降這兩個基因會引起dnmt3a缺失，進一步證明了DNA甲基化與胰腺β細胞的功能直接作用關(guān)系［27］。MT-CO3與糖尿病有關(guān)［28］。過表達mir-10a會抑制HOXA1并引起胰腺癌局部侵襲［29］。通過研究表明，mir-320c通過SMARCC1在治療胰腺癌中對藥物吉西他賓（gemcitabine） 產(chǎn)生抗藥性［30］。模塊四當中發(fā)現(xiàn)新的集合mir-122-5p、EGLN3、EPAS1，其中mir-122與多個miRNAs篩選作胰腺癌血液中的標記物［31］；EGLN3與胰腺癌腫瘤細胞的分化、增殖有關(guān)［32］；干擾EPAS1用作治療胰腺癌疾病［33］。模塊五中YWHAH、FOXO3、CDC25C、YWHAB、YWHAZ與模塊二中的集合相同，mir-125a-5p、TCEB2、EPAS1單獨聚在一起，降低mir-125a-5p的表達，會抑制胰腺腫瘤的增殖和細胞周期，并促進早期癌細胞的凋亡［34］。TCEB2m則可能用于卵巢癌的治療［35］。

2.3 模塊功能分析

五個模塊分別進行GO功能注釋（圖5），以上五個模塊注釋來看，在多個模塊或同一模塊當中多次出現(xiàn)注釋到一對相反的功能上，如:positive regulation of biological process和process、negative regulation of biological process；positive regulation of metabolic process、 negative regulation of metabolic process等，這表明在組織自身原來就存在著平衡調(diào)控機制，這些功能的失衡可能會引起疾病發(fā)生。多個模塊注釋到代謝相關(guān)的功能上，從功能上說明疾病也與代謝密切相關(guān)。另外除了細胞周期的功能之外，還有甲基化功能也被注釋。通過功能注釋表明多個基因通過協(xié)同在多個模塊中來行使某一功能，而多個相同的功能可能通過多個基因在不同的途徑下實現(xiàn)。

圖4 miRNAs-mRNAs模塊圖

圖5 模塊功能注釋

3 結(jié)語

為了明確胰腺癌中miRNAs與mRNAs如何行使調(diào)控作用，通過對胰腺癌芯片數(shù)據(jù)進行分析，得到miRNA-mRNAs網(wǎng)絡(luò)，對網(wǎng)絡(luò)進行模塊挖掘獲得5個顯著性模塊，對模塊進行分析，mir-615-3p、mir-222、mir-93-5p、mir-415-3p、mir-877在網(wǎng)絡(luò)中起到了HUB節(jié)點的作用，其中mir-222在表達譜中差異表達。在模塊當中還有一些節(jié)點度雖然比較小，通過聚類方式發(fā)現(xiàn)這些結(jié)點以“團隊”合作的方式行使功能，如:YWHAH、FOXO3、CDC25C、YWHAB、YWHAZ，通過文獻也發(fā)現(xiàn)這些成員與胰腺癌的發(fā)生都有密切關(guān)系。從功能上多個模塊中出現(xiàn)相反的功能注釋，如:positive regulation of cellular process、negative regulation of cellular process。一方面表明這些模塊中存在著的行使相同功能的基因，另一方面也表明，在同一模塊當中一些基因在平衡著這些功能起到重要作用。從功能注釋到甲基化，說明癌癥的發(fā)生機制十分復雜，是多階段、多因素復雜作用的結(jié)果。

最后，通過本文的研究，發(fā)現(xiàn)miRNAs與mRNAs協(xié)同作用的一些結(jié)果，為今后更深入的研究胰腺癌提供指導和借鑒。

［1］ Mayo SC，Austin DF，Sheppard BC，et al.Evolving pre—operative evaluation of patients with pancreatic cancer:does laparoscopy have a mle in the cuⅡ℃nt era？［J］.J Am Coll Surg，2009，208（1）:87-95.

［2］ Galasso D，Camuccio A，Larghi A，et al.Pancreatic cancer:diagnosis and endoscopic staging［J］.Eur Rev Med Pharmacol Sci，20lO，14（4）:375-385.

［3］ Calin G A，Corce C M. Nat. Rev. Cancer，2006，6:857-866.

［4］ Gilad S，Meiri E，Yogev Y，Benjanmin S，Lebanony D，Yerushalmi N，Bentwich Z ，Hod M，Goren Y，Goren Y，Chajut A. PloS ONE，2008，3:e314.

［5］ Mitchell P S，Parkin R K，Kroh E M，F(xiàn)ritz B R，Wyman S K，Pogosova-Agadjanyan E L，Peterson A，Noteboom J，O’Briant K C，Allen A，Lin D W，Urban N，Drescher C W，Knudser B S，Stirewlt D L，Gentleman R ，Vessella R L，Nelson P S，Martin D B，Tewari M. Proc. Natl. Acad. Sci.，2008，105:10513-10518.

［6］ Ryu JK，Hong SM，Karikari CA，et al.Abe玎ant microRNA-155 expression is an early event in the multistep progression of pancreatic adenocarcinoma［J［.Pancreato-logy，2010，10（1）:66-73.

［7］ Habbe N，Koorstra JB，Mendell JT，et al.MicroRNA miR一155 is a biomarker of early pancreatic neoplasia［J］.Cancer Biol Ther，2009，8（4）:340-346.

［8］ Wang J，Chen J，Chang P，et al.MicroRNAs in plasma of pancreatic ductal adenocarcinoma patients as novel blood—based biomarkers of disease［J］.Cancer Prev Res（Phila），2009，2（9）:807—813.

［9］ Li A，0mura N，Hong SM，et al.Pancreatic cancers epigenetically silence SIPl and hypomethylate and overexDress miR.200a/ 200b in association with elevated circulating miR一200a and miR一200b levels［J］.Cancer Res，2010，70（13）:5226—5237.

［10］ Bloomston M，F(xiàn)rankel WL，Petmcca F，et al.MicroRNA expression pattems to diffrentiate pancreatic adenocarci nomal pancreas and chronic pancreatitis［J］.JAMA，2007，297（17）:1901—1908.

［11］ Zhao Y，Wang Y，Yang Y，Liu J，Song Y，Cao Y，Chen X，Yang W，Wang F，Gao J，Li Z，Yang C. MicroRNA-222 Controls Human Pancreatic Cancer Cell Line Capan-2 Proliferation by P57 Targeting. J Cancer. 2015 Oct 16；6（12）:1230-5.

［12］ Huang J，Wang B，Hui K，Zeng J，F(xiàn)an J，Wang X，Hsieh JT，He D，Wu K. miR-92b targets DAB2IP to promote EMT in bladder cancer migration and invasion. Oncol Rep. 2016 Sep；36（3）:1693-701.

［13］ 李捍司.miRNA-222參與調(diào)控胰腺癌細胞侵襲轉(zhuǎn)移分子機制的研究［D］.中國醫(yī)科大學，2014.

［14］ Ye S，Yang L，Zhao X，Song W，Wang W，Zheng S. Bioinformatics method to predict two regulation mechanism:TF-miRNA-mRNA and lncRNA-miRNA-mRNA in pancreatic cancer. Cell Biochem Biophys. 2014 Dec；70（3）:1849-58.

［15］ Namkoong H，Shin SM，Kim HK，Ha SA，Cho GW，Hur SY，Kim TE，Kim JW. The bone morphogenetic protein antagonist gremlin 1 is overexpressed in human cancers and interacts with YWHAH protein. BMC Cancer. 2006 Mar 18；6:74.

［16］ Zhu Y，Xu H，Chen H，Xie J，Shi M，Shen B，Deng X，Liu C，Zhan X，Peng C.Proteomic analysis of solid pseudopapillary tumor of the pancreas reveals dysfunction of the endoplasmic reticulum protein processing pathway. Mol Cell Proteomics. 2014 Oct；13（10）: 2593-603.

［17］ Szeliga M，?wik?a J，Obara-Michlewska M，Cichocki A，Albrecht J. Glutaminases in slowly proliferating gastroenteropancreatic neuroendocrine neoplasms/tumors（GEP-NETs）:Selective overexpression of mRNA coding for the KGA isoform. Exp MolPathol. 2016 Feb；100（1）:74-8.

［18］ Hu G，Agarwal P. Human disease-drug network based on genomic expression profiles. PLoS One. 2009；4（8）:e6536.

［19］ Jiang Y，Zhang Y，Li F，Du X，Zhang J. CDX2 inhibits pancreatic adenocarcinoma cell proliferation via promoting tumor suppressor miR-615-5p. Tumour Biol. 2016 Jan；37（1）:1041-9.

［20］ Peters IR，Peeters D，Helps CR，Day MJ. Development and application of multiple internal reference （housekeeper） gene assays for accurate normalisation of canine gene expression studies. Vet Immunol Immunopathol. 2007 May 15；117（1-2）:55-66.

［21］ Zhang ZL，Bai ZH，Wang XB，Bai L，Miao F，Pei HH. miR-186 and 326 predict the prognosis of pancreatic ductal adenocarcinoma and affect the proliferation and migration of cancer cells. PLoS One. 2015 Mar 5；10（3）:e0118814.

［22］ Deng S，Li X，Niu Y，Zhu S，Jin Y，Deng S，Chen J，Liu Y，He C，Yin T，Yang Z，Tao J，Xiong J，Wu H，Wang C，Zhao G. MiR-652 inhibits acidic microenvironment-induced epithelial- mesenchymal transition of pancreatic cancer cells by targeting ZEB1. Oncotarget. 2015 Nov 24；6（37）:39661-75.

［23］ Wu JX，Hong YH，Yang XG. Bis（acetylacetonato）-oxidovanadium（IV） and sodium metavanadate inhibit cell proliferation via ROS-induced sustained MAPK/ERK activation but with elevated AKT activity in human pancreatic cancer AsPC-1 cells. J Biol Inorg Chem. 2016 Sep 10.

［24］ Wang J，He P，Gaida M，Yang S，Schetter AJ，Gaedcke J，Ghadimi BM，Ried T，Yfantis H，Lee D，Weiss JM，Stauffer J，Hanna N，Alexander HR，Hussain SP.Inducible nitric oxide synthase enhances disease aggressiveness in pancreatic cancer. Oncotarget. 2016 Jun 29.

［25］ Zhou W，Li Y，Gou S，Xiong J，Wu H，Wang C，Yan H，Liu T. MiR-744 increases tumorigenicity of pancreatic cancer by activating Wnt/β-catenin pathway.Oncotarget. 2015 Nov 10；6（35）:37557-69.

［26］ Deng L，Shang Y，Guo S，Liu C，Zhou L，Sun Y，Nie Y，F(xiàn)an D，Lu Y，Guo X. Ran GTPase protein promotes metastasis and invasion in pancreatic cancer by deregulating the expression of AR and CXCR4. Cancer Biol Ther. 2014 Aug；15（8）:1087-93.

［27］ Dhawan S，Tschen SI，Zeng C，Guo T，Hebrok M，Matveyenko A，Bhushan A. DN methylation directs functional maturation of pancreatic β cells. J Clin Invest.2015 Jul 1；125（7）:2851-60.

［28］ Tabebi M，Mkaouar-Rebai E，Mnif M，Kallabi F，Ben Mahmoud A，Ben Saad W，Charfi N，Keskes-Ammar L，Kamoun H，Abid M，F(xiàn)akhfakh F. A novel mutation MT-COIII m.9267G>C and MT-COI m.5913G>A mutation in mitochondrial genes in a Tunisian family with maternally inherited diabetes and deafness （MIDD） associated with severe nephropathy. Biochem Biophys Res Commun. 2015 Apr 10；459（3）:353-60.

［29］ Ohuchida K，Mizumoto K，Lin C，Yamaguchi H，Ohtsuka T，Sato N，Toma H，Nakamura M，Nagai E，Hashizume M，Tanaka M. MicroRNA-10a is overexpressed in human pancreatic cancer and involved in its invasiveness partially via suppression of the HOXA1 gene. Ann Surg Oncol. 2012 Jul；19（7）:2394-402.

［30］ Iwagami Y，Eguchi H，Nagano H，Akita H，Hama N，Wada H，Kawamoto K，Kobayashi S，Tomokuni A，Tomimaru Y，Mori M，Doki Y. miR-320c regulates gemcitabine-resistance in pancreatic cancer via SMARCC1. Br J Cancer. 2013 Jul 23；109（2）:502-11.

［31］ Schultz NA，Dehlendorff C，Jensen BV，Bjerregaard JK，Nielsen KR，Bojesen SE，Calatayud D，Nielsen SE，Yilmaz M，Holl?nder NH，Andersen KK，Johansen JS. MicroRNA biomarkers in whole blood for detection of pancreatic cancer. JAMA. 2014 Jan 22-29；311（4）:392-404.

［32］ Su Y，Loos M，Giese N，Hines OJ，Diebold I，G?rlach A，Metzen E，Pastorekova S，F(xiàn)riess H，Büchler P. PHD3 regulates differentiation，tumour growth and Angio genesis in pancreatic cancer. Br J Cancer. 2010 Nov 9；103（10）:1571-9.

［33］ Pan X，Zhu Q，Sun Y，Li L，Zhu Y，Zhao Z，Zuo J，F(xiàn)ang W，Li K. PLGA/poloxamer nanoparticles loaded with EPAS1 siRNA for the treatment of pancreatic cancer in vitro and in vivo. Int J Mol Med. 2015 Apr；35（4）:995-1002.

［34］ Jia CW，Sun Y，Zhang TT，Lu ZH，Chen J. Effects of miR-125a-5p on Cell Proliferation，Apoptosis and Cell Cycle of Pancreatic Cancer Cells. Zhongguo Yi Xue Ke Xue Yuan Xue Bao. 2016 Aug；38（4）:415-21.

miRNAs and mRNAs gene expression profile analysis in pancreatic cancer

BAI Xue-feng， LIU Yue-juan， QI jing， ZHANG jian
（Department of Medical Informatics， Harbin Medical University， Daqing， 163319， China）

In this paper，we study miRNAs-mRNAs regulatory mechanism about pancreatic cancer。Application of bioinformatics method，we set up the miRNAs-mRNA disease network by the pancreatic cancer expression profile，and obtain five significance functional module by mining module of plugins .The larger degree nodes of miRNAs play a role HUB node in this network，while these are closely related to the initial and development of pancreatic cancer. whereafter use of clustering method to each module. the small degree nodes are gathered a collection，of these are more closely related to pancreatic cancer. The GO annotation of genes showed that pancreatic cancer was closely correlated to the metabolism function； moreover，many of same functional annotation appear in multiple modules，while many contrary of functional annotation appear in the single module. This suggests that there is a balance in its original regulatory mechanism.In this paper，achieved some synergy results about of miRNAs-mRNAs，these results may provide reference and help for future disease research on pancreatic cancer。

Pancreatic cancer；mRNAs；mRNAs；Synergy

黑龍江省衛(wèi)生廳科研課題（編號:2013116）

白雪峰，男，講師，主要研究方向生物信息學。