王 涛，潘 鑫

    (锦州医科大学医疗学院医学系，辽宁 锦州 121013)

    miRNA 作为一类小的非编码RNA，其长度大约在22 个核苷酸左右，能够与靶mRNA 3'端非翻译区互补，进而降解或阻止该mRNA 翻译出蛋白质，实现转录后水平的基因表达调控[1]。糖尿病作为一种慢性疾病，部分是由胰岛β 细胞受损而引起，胰腺或胰岛移植对于该类疾病治疗效果较好，但是由于供体缺乏以及免疫排斥导致临床应用受限[2，3]，而干细胞分化为胰岛素分泌细胞(insulin-producing cells，IPCs)再行移植，有望彻底根治[4，5]。生物信息学通过整合多个学科领域信息和知识，能够从整体层面揭示疾病或生物过程的复杂分子机制[6-8]。本研究运用生物信息学方法分析人胚胎干细胞(human embryonic stem cells，hESCs)向IPCs 诱导分化过程中mRNA 的表达谱，预测靶向作用于该mRNA 的miRNA，构建miRNA-mRNA 调控网络。通过对调控网络的分析，拟找到IPCs 诱导分化的关键节点，为进一步开展分子靶向干预提供新思路。

    1 资料与方法

    1.1 资料来源 选择GEO(https://www.ncbi.nlm.nih.gov/geo/)数据库内的GSE42094 数据集作为研究对象，该数据集包含23 个样本数据(GSM1032319-41)，涉及未分化hESCs、IPCs 分化的5 个时期、胰腺内胰岛和胎胰。选择未分化hESCs(hESCs 组)和IPCs 分化的5 个时期(Diff1、Diff2、Diff3、Diff4 和IPCs 组)研究mRNA 的表达谱。应用数据库自带软件GEO2R对6 组数据进行分析，获取差异表达的mRNA。

    1.2 GO 功能和KEGG 路径分析 选用DAVID(https://david.ncifcrf.gov/tools.jsp)数据库进行功能富集分析，主要分析GO 功能以及KEGG 路径，以P<0.05为差异有统计学意义。

    1.3 蛋白互作网络分析 选用STRING(https://stringdb.org/)数据库进行蛋白互作网络(protein-protein interaction，PPI)分析，选择Homo sapiens 作为“Organism”的条件，其余参数均为系统默认设置。

    1.4 miRNA 预测 针对差异表达基因预测其靶miRNA，选用miRTarBase(https://mirtarbase.cuhk.edu.cn/)数据库 ......