hsa-miR-216a-5p的靶基因及調控網絡的生物信息學分析

馬 亮，張運峰，董建悅，胡 芬

生命科學與技術

hsa-miR-216a-5p的靶基因及調控網絡的生物信息學分析

馬 亮1，張運峰2，董建悅2，胡 芬1

（1. 華北理工大學 生命科學學院，河北 唐山 063210；2. 唐山師范學院 生命科學系，河北 唐山 063000）

利用miRBase、生物信息工具bioinfo和TargetScan7.2、miRDB、miRDIP和miTarBase四個數據庫獲取了八個物種的miR-216a-5p成熟序列、表達譜數據和129個靶基因。對靶基因集合進行GO、KEGG分析，結果顯示hsa-miR-216a-5p主要富集在轉錄DNA模板、轉錄生長因子beta受體信號通路、蛋白穩定化等10個生物過程（P<0.01），泛素蛋白連接酶活性、SMAD蛋白結合、蛋白泛素化等4個分子功能（P<0.01）以及糖尿病、內吞作用、TRF通道的炎性介質調節等5個KEGG通路中（P<0.05）。借助構建PPI蛋白互作網絡，從129個靶基因中篩選出了PTEN、SIRT1、TGFBR2等10個關鍵基因。

hsa-miR-216a-5p；靶基因；生物信息

miRNA（microRNA）是一類廣泛存在于動植物和某些病毒體內長度為22nt左右的單鏈非編碼RNA分子，通過結合靶mRNA的3’UTR區阻止其翻譯或促進其降解而發揮功能[1]。miRNA具有高度的序列保守性，從而能夠廣泛調節多種生命活動。miR-216a-5p在多種癌癥的發生發展中具有調控作用。在肝細胞癌（HCC）的發生過程中，mR-216a-5p可通過JAK2靶向抑制癌細胞的增殖和侵襲[2]，通過下調MMP16抑制肺癌細胞的侵襲[3]，也能以PAK2為靶基因實現對膀胱癌細胞增殖和凋亡的調控[4]，另外還可抑制HMGB1的表達從而抑制胃癌細胞的增殖和遷移[5]。

到目前為止，hsa-miR-216a-5p眾多靶基因尚不十分清楚。本研究擬采用生物信息學的方法利用公共數據庫對hsa-miR-216a-5p進行靶基因預測，并采用DAVID網站對獲得靶基因集合進行GO功能注釋和KEGG通路富集分析，利用STRING數據庫和Cystocape 3.7.1軟件構建蛋白互作網絡并篩選出關鍵基因，為進一步深入分析hsa-miR-216a-5p的功能提供參考依據。

1 材料與方法

1.1 miR-216a-5p序列信息獲取

通過在線公共數據庫miRBase（http://www. mirbase.org）下載人、獼猴等八個物種的miR-216a-5p的成熟序列，利用MEGA-X軟件分析其序列保守性。

1.2 hsa-miR-216a-5p表達譜分析

采用癌癥基因組圖譜計劃（TCGA）數據庫信息，利用生物信息學在線分析工具（http:// bioinfo.life.hust.edu.cn/miR_path/index）獲得hsa-miR-216a-5p在人類各個癌組織和癌旁組織的表達譜。

1.3 靶基因預測

分別采用數據庫TargetScan7.2、miRDB、miRDIP進行hsa-miR-216a-5p的靶基因預測，使用venny2.1進行韋恩圖繪制。利用miTarBase數據庫獲取經過實驗驗證的靶基因，綜合四個數據庫整合靶基因集合，用于后續分析。

1.4 GO和KEGG分析

利用在線網站DAVID（https://david.ncifcrf. gov）對得到的hsa-miR-216a-5p靶基因集合，分別進行生物過程、細胞組分、分子功能GO功能富集分析和KEGG通路富集分析。

1.5 靶基因PPI網絡構建及關鍵基因篩選

將上述獲得的靶基因集合輸入STRING（https://string-db.org/cgi/input.pl）數據庫，構建PPI蛋白互作網絡，并結合Cystocape 3.7.1軟件進行數據可視化，通過使用CytoHubba插件進行關鍵基因的篩選。

2 結果

2.1 hsa-miR-216a-5p基因序列保守性分析

如圖1，在人、獼猴、小鼠、大鼠、安大略鮭中miR-216a-5p的成熟序列為22個堿基，在兔、犰狳、豚鼠中為23個堿基，八個物種miR- 216a-5p成熟序列基本一致。

圖1 miR-216a-5p 在8個物種中成熟microRNA序列比對結果

2.2 hsa-miR-216a-5p表達譜分析

如圖2，相較于正常組織，在膀胱尿路上皮癌（BLCA）、宮頸癌（CESC）、食管癌（ESCA）等中呈現高表達，miR-216a-5p在這些癌癥中可能充當癌基因。另外在腎嫌色細胞癌（KICH）、腎乳頭狀細胞癌（KIRP）中呈現低表達，推測其可能作為抑癌基因存在。

圖2 TCGA分析hsa-miR-216a-5p的表達譜。

2.3 hsa-miR-216a-5p靶基因預測

采用數據庫TargetScan、miRDB、miRDIP預測hsa-miR-216a-5p的靶基因，分別獲得372個、577個、845個靶基因，對數據取交集共獲得122個靶基因，詳見圖3a。結合miTarBase數據庫經過雙熒光素酶報告實驗、western blot和qPCR驗證的靶基因10個（PTEN, SIRT1, CDC42, CD44, SMAD7, BECN1, HNF4A, JAK2, CEMIP, CBL），共獲得129個靶基因，詳見圖3b。

2.4 hsa-miR-216a-5p靶基因的GO功能注釋和KEGG通路富集

采用DAVID對hsa-miR-216a-5p靶基因進行GO分析，結果表明靶基因主要富集在轉錄DNA模板、轉錄和轉化生長因子beta受體信號通路、表皮生長因子激活受體的隱性調節等10個生物過程（P<0.01），泛素蛋白連接酶活性、SMAD蛋白結合、蛋白泛素化和蛋白結合4個分子功能（P<0.01），詳見圖4。

圖4 hsa-miR-216a-5p靶基因GO功能顯著性富集分析結果

使用在線網站DAVID進行了靶基因KEGG通路富集分析，結果發現共有47個靶基因富集到KEGG通路上，其中主要富集在年輕的成年發病型糖尿病、內吞作用、TRF通道的炎性介質調節、細胞黏著和癌癥microRNA等5個通路（P<0.05）。結果見表1。

表1 hsa-miR-216a-5p靶基因KEGG通路富集分析結果

2.5 構建靶基因PPI網絡并篩選關鍵基因

通過上述實驗獲得的hsa-miR-216a-5p靶基因集結合STRING數據庫構建其蛋白互作網絡圖（PPI），得到一個有76個node，123個edge的蛋白互作網絡圖。下載節點數據文件導入Cytoscape 3.7.1，使其數據可視化，利用CytoHubba插件計算PPI中各個基因的節點度（degree），篩選核心蛋白質的編碼基因。進而得到10個關鍵基因：PTEN、SIRT1、TGFBR2、CDC42、XIAP、CD44、SMAD7、CBL、ZEB1和HNF4A，詳見圖5。

圖5 has-miR-216a-5p關鍵靶基因篩選

3 討論

miRNA與眾多生理過程有密切關系，如分化、增殖、凋亡與發育等[6]。目前，miR-216a-5p的研究主要集中于癌癥的發生發展中。CHEN等[7]研究發現，在腎細胞癌中miR-216a-5p顯著上調表達，并促進癌細胞增殖抑制其凋亡，可作為腎細胞癌的標志物，但尚未給出其靶向調控基因。然而在肺細胞癌中miR-216a-5p卻通過靶向調控Bcl-2、Bax、Bad蛋白（family proteins），從而抑制細胞增殖和轉移并影響細胞周期活動[8]；Zhang Y等人[9]研究證實，在乳腺癌的發生發展中，miR-216a-5p通過靶向調控PAK2蛋白的表達進而抑制癌細胞增殖和轉移。通過這些研究可知，miR-216a-5p是調控癌癥細胞增殖、轉移和凋亡的關鍵因子，由此可推測，其可作為癌癥早期診斷的生物標志物。

miRNA作為一種小型單鏈RNA分子越來越多地受到關注。研究表明，其主要通過和靶mRNA的3’UTR區結合抑制其翻譯或促進其降解而發揮功能，所以往往伴隨著眾多靶基因的異常表達。如何快速準確地預測其靶基因顯得尤為重要。生物信息學方法預測mRNA靶基因的工具很多，目前為止使用較多的有TargetScan7.2、miRDB、miRDIP這三種[10,11]。其中miRDB為第一代預測軟件，其不受物種限制，支持在線和本地，支持各種平臺，產出多但精度低。TargetScan7.2為第二代預測軟件，預測精度大為提高，假陽性較低；miRDIP為最新預測軟件，搜索功能強大，整合了多達30個數據庫的資源但精度低。這三者所使用的算法各有優勢，能互為補充，還可降低假陽性率和擴大搜索范圍。本研究整合經過實驗驗證數據庫miTarBase結果，共獲得了129個靶基，為進一步實驗工作指明了方向。

[1] Iorio M, Croce C. MicroRNAs in cancer: small molecules with a huge impact[J]. Journal of Clinical Oncology, 2009, 27(34): 5848-5856.

[2] 周維,譚曉宇,李慧芬,等. miR-216a-5p通過JAK2靶向抑制肝癌細胞的增殖和侵襲[J].中國醫藥生物技術,2019,14(2):155-163.

[3] 崔艷麗,李小亮,何利珍.miR-216a-5p通過抑制HMGB1降低胃癌細胞系MGC-803的增殖和遷移[J].基礎醫學與臨床,2019,39(4):546-551.

[4] 邵煥軍,趙振伶,郝麗娜,等. miR-216a-5p靶向作用于PAK2對膀胱癌細胞增殖和凋亡影響的體外研究[J].醫學研究雜志,2018,47(6):151-155.

[5] 安寧,李宏敏,于瑞蓮,等.miR-216a-5p通過下調MMP16表達抑制肺癌細胞的侵襲[J].中國癌癥雜志,2015, 25(8):588-594.

[6] He L, Hannon G J. MicroRNAs: Small RNAs with abig role in gene regulation[J]. Nature Reviews Genetics, 2004, 5(7): 522- 531.

[7] Chen P, Quan J, Jin L, et al. miR-216a-5p acts as an oncogene in renal cell carcinoma[J]. Exp Ther Med, 2018, 15(4): 4039-4046.

[8] Sun Y, Hu B, Wang Y, et al. miR-216a-5p inhibits malignant progression in small cell lung cancer: involvement of the Bcl-2 family proteins[J]. Cancer Manag Res, 2018, 18(10): 4735-4745.

[9] Zhang Y, Lin P, Zou J Y, et al. MiR-216a-5p act as a tumor suppressor, regulating the cell proliferation and metastasis bytargeting PAK2 in breast cancer[J]. Eur Rev Med Pharmacol Sci, 2019, 23(6): 2469-2475.

[10] Nathan W, Wang X W.miRDB: an online resource for microRNA target prediction and functional annotations [J]. Nucleic Acids Research, 2015, 43(D1): D146-152.

[11] Tokar T, Pastrello C, Rossos AEM, et al. mirDIP 4. 1-integrative database of human microRNA target predictions[J]. Nucleic Acids Res, 2018, 46(D1): D360-D370.

Bioinformatics Analysis of Target Gene and Regulatory Network of hsa-miR-216a-5p

MA Liang1, ZHANG Yun-feng2, DONG Jian-yue2, HU Fen1

(1. College of Life Sciences, North China University of Science and Technology, Tangshan 063210, China; 2. Department of Life Science, Tangshan Normal University, Tangshan063000, China)

The mature sequences of miR-216a-5p were obtained by miRBase, and its conservation was analyzed by MEGA-X. The expression profile of has-miR-216a-5p was get from bio info web. A total of 129 target genes were obtained integrating four database target genes of TargetScan7.2、miRDB、miRDIP and miTarBase. The GO function annotation, KEGG pathway analysis and PPI were performed on the target genes. The results show that they are mainly enriched in 10 biological processes such as transcription and DNA-templated (p<0.01), 4 molecular functions such as ubiquitin protein ligase activity (p<0.01), and 5 KEGG pathways such as maturity onset diabetes of the young (p<0.05). Ten key genes such as PTEN, SIRT1, TGFBR2 were screened out from 129 target genes by constructing PPI protein interaction network.

hsa-miR-216a-5p; target genes; bioinformatics

Q279

A

1009-9115(2020)03-0067-04

10.3969/j.issn.1009-9115.2020.03.016

河北省自然科學基金項目（H2018209140，C2019105055）

2019-12-16

2020-04-06

馬亮（1994-），男，河北邯鄲人，碩士研究生，研究方向為腫瘤遺傳學。

胡芬（1983-），女，河北唐山人，博士，副教授，研究方向為生物化學與分子生物學。

（責任編輯、校對：李春香）