肛瘘患者特征性微小RNA表达的检测及分析

目的：检测并分析肛瘘患者的特征性微小RNA（miRNA）表达谱，探讨其可能的靶基因和临床意义。方法收集2014年3—12月间杭州市第三人民医院肛肠外科收治的3例肛瘘联合混合痔切除手术患者的痔核近端肛管黏膜标本（肛瘘组），另选取年龄、性别和体质量与之相匹配的行单纯混合痔手术的3例患者的肛管黏膜标本作为对照（对照组）。利用miRNA微阵列芯片检测比较两组肛管黏膜的1285个miRNAs的表达差异，并通过聚类分析差异miRNAs的聚集性（差异基因的筛选标准是：P ≤0．05且变化倍数≥2或≤0．5）；采用DIANAmT、miRanda、 miRDB、 miRWalk等10种预测软件对差异miRNAs的靶基因进行预测，并进行综合打分，确定出预测分值最高的基因；采用基因本体（GO）富集法分析与靶基因相关的生物学过程；利用免疫组织化学方法检测两组肛管黏膜中分值最高基因的蛋白表达。结果肛瘘组患者的1285个miRNA中，发现13个miRNA与对照组存在显著差异，其中在肛瘘组上调2个，下调11个。经配对t检验显示，肛瘘组中miR-3609上调是对照组的5．98倍（P＝0．0231），miR-181a-2-3p下调0．13倍（P＝0．0067），是差异表达最大的miRNA。进一步聚类分析显示，上调的hsa-miR-3609和hsa-miR-6086相对表达量在两组间的变化趋势一致；下调的11个miRNAs中，miR-125b-1-3p和miR-548q表达类似，其余9个表达相似（包括：hsa-miR-1185-1-3p、hsa-miR-532-3p、hsa-miR-1233-5p、hsa-miR-769-5p、hsa-miR-149-5p、hsa-miR-99b-3p、hsa-miR-138-5p、hsa-miR-141-3p及hsa-miR-181a-2-3p）。对13个差异miRNAs进行靶基因预测分析显示，其中7个（53．8％）miRNAs能预测到靶基因，共计预测到靶基因104个；其余6个（46．2％）miRNAs未能预测到靶基因。预测分值最高的为几丁质酶基因1（CHIT1），其对应的差异miRNA为hsa-miR-769-5p （r ＝-0．94286，P ＝0．0167）。基因本体（GO）富集分析结果显示，与这104个靶基因最为相关的生物学过程包括角质化、免疫反应和信号传导等。免疫组织化学检测显示，肛瘘组肛管黏膜组织中CHIT1表达明显高于对照组（P ＜0.01）。结论肛瘘患者存在特征性的miRNA表达，其可能在肛瘘发病过程中起了一定的作用。