肺动脉高压(pulmonary hypertension，PH)主要特征为肺动脉阻塞引起肺血管阻力及肺动脉压力渐进性升高，伴随不可逆的肺血管重构，最终导致右心衰竭而死亡^[1]。PH时肺血管重构主要表现为内膜增生、中膜肥厚、外膜增生、原位血栓、不同程度的炎症以及丛状动脉样改变^[2-3]。PH时重构的肺血管内膜中α-SMA标记的细胞明显增加^[1]。传统的观念认为静止的肺动脉平滑肌细胞(pulmonary arterial smooth muscle cell，PASMCs)自身增殖是这些细胞的唯一来源。近期有学者在PH动物模型和患者病变肺血管中发现了内皮细胞间质转化(endothelial-mesenchymal transition, EndMT)的存在，并通过体外实验证实雷帕霉素可抑制各种因素诱导的肺血管EndMT^[4-5]。提示，肺血管EndMT是肺血管重构时α-SMA样细胞的另一重要来源。进一步探索PH时肺血管EndMT的病理机制，为寻找防治PH的新靶点和药物具有重要意义。

微小非编码RNA(microRNA，miRNA)是一类存在于真核生物中具有调控功能的非编码小分子单链RNA，大约长度为21~25 nt。大量研究表明其参与生命过程中一系列的重要进程，包括个体发育、器官形成以及细胞增殖、死亡与分化等^[6]。众所周知，一个miRNA有多个靶点，一个基因又可被多个miRNAs调控，他们之间往往形成复杂的网络调控机制参与疾病发生发展。因此，通过网络药理学技术探索相互作用网络，而不是单一的研究某个miRNA或基因的功能，能更好地解释miRNA或基因在疾病中的整体作用。

目前研究虽然提示miRNAs在PH的病理生理过程中扮演着重要的角色。但仅发现了少部分miRNAs在PH中起重要作用，且对于这类疾病分子机制的整体功能研究也尚在起步阶段。鉴于PH时肺血管EndMT的机制尚不明确，本文拟运用网络药理学方法探索与PH时肺血管EndMT相关miRNAs及其相关下游靶点的相互作用机制，为后续研究PH时肺血管EndMT的分子机制提供明确的方向及策略。

通过Medline(PubMed)搜索“pulmonary hypertension”，通过此我们寻找到与PH疾病相关的基因，我们称为PH模块。通过进一步切换关键词，包含“pulmonary hypertension”、“microRNA”、“miRNA”、“endothelial to mesenchymal transition”、“epithelial-mesenchymal transition”。搜索得到与PH时EMT/EndMT相关的miRNAs及EMT/EndMT相关的miRNAs。

从Medline(PubMed)基因数据库中随机选取与PH模块个数相同的蛋白编码基因，作为随机网络模块。随机网络模块重复选取3次。同时将PH模块及随机网络模块中的基因分别通过蛋白质-蛋白质相互作用软件数据库BioGRID(Biological General Repository for Interaction Datasets)数据库得到各自模块基因相互作用关系，并计算各自网络的最大连接组件的“节点数”(largest connected component，LCC)和直接相互作用的“边”(edges)。用Cytocape 3软件做相互作用网络图，孤立的节点不在图上画出。

通过生物信息数据库(DIANA3、Miranda4、PicTar5、TargetScan6、miRDB7和microT-CDS8) 进行搜索得到相应miRNAs靶基因关系。

通过miRNAtap9工具包超几何分析列表中每个miRNA的靶基因在PH模块中的富集效应。具体方法为选择PH模块中的230个基因作为基因集合，每个miRNAs预测得到的所有miRNAs为基因列表，以在2个或以上不同靶基因预测数据源中出现过的基因为全部基因背景，四联表卡方分析，采用蒙特卡洛模拟方法计算确切显著性概率值P值。对于在基因集合中出现基因列表中的基因的比例小于随机概率的，P值赋予0.95，即直接认定富集不显著。校正P值计算用到BH(Benjamini Hochberg)方法。通过设定筛选条件[校正P值＜0.05, 靶基因在基因集合中富集的基因个数(setsize)≥5]得出与PH时肺血管EndMT相关miRNAs。将参与PH疾病的功能通路主要分为TGF-β、低氧和炎症，通过miRNAs基因靶点所属功能类别进一步对miRNAs进行功能分类评分，筛选出评分较高的miRNAs，并用Cytocape 3软件构建miRNAs-靶点相互作用网络图。

清洁级Sprague-Dawley(SD)♂大鼠20只，合格号为scxk(湘)2009-0004，体质量180~200 g，适应性喂养1周后称量体质量、标号；按体质量随机分为对照组(n=10) 和低氧模型组(n=10) 。对照组常氧下(21% O₂)饲养；低氧模型组于低氧仓中(10% O₂，放有无水氯化钙和钠石灰分别用来吸收水分和CO₂)饲养，两组大鼠自由饮水、进食。3周后颈静脉插管分别检测右心室收缩压(right ventricular systolic pressure, RVSP)及肺动脉平均压(mean pulmonary artery pressure, mPAP)以确定造模是否成功。

肺动脉用液氮磨碎后加TRIzol。根据逆转录试剂盒(TaKaRa)说明书逆转录为cDNA。cDNA经适当稀释后，用SYBR Premix Ex TaqⅡ试剂盒(TaKaRa)进行实时荧光定量PCR。PCR反应条件：95℃ 30 s，(95℃ 5 s，60℃ 31 s)40个循环。内参选用U6，内参及目的miRNAs引物由广州锐博生物有限公司提供。

所有数据均用x±s表示。统计分析采用SPSS 19.0软件。两组均数比较采用独立样本的Student′s t检验。

通过Medline(PubMed)搜索pulmonary hypertension(2016年2月14日)，共50 307篇文献，寻找到与PH疾病相关的基因共230个，我们将其称为PH模块。将这230个基因通过22个不同功能进行分类并做相互作用关系图(Fig 1A)。图中每个节点表示一个基因，颜色为该基因对应的功能；字体大小表示该节点的度的大小，度越大节点字体越大；每条边的粗细表示该两两基因作用的研究报道支撑多寡，线条越粗则该两两作用在不同研究中报道的次数越多。分别计算PH模块及随机网络模块的网络属性得到PH模块共有195个节点，LCC为76，Edges为727；随机模块LCC为3.7±3，Edges为36.7±3.5。PH模块LCC及Edges显著高于随机网络模块(Fig 1B~C，P＜0.01) 。表明PH模块具有巨大而稠密的相互关联，是作为筛选相关miRNAs的理想模块。

通过Medline(PubMed)搜索得到PH时EndMT相关的文献23篇；PH时EndMT相关miRNAs文献0篇；EndMT/EMT有关的miRNAs文献56篇(2016年2月14日)。通过阅读摘要，必要时阅读全文，得到相关miRNAs189个，对应成熟体322个。通过富集分析发现可能与PH时EndMT相关的miRNAs共98个，将其称为“目标miRNAs”(Tab 1)。

Fig 1 The PH-network

图选项

由于“目标miRNAs”数目较多，为了筛选评分最高的miRNAs，将参与PH疾病的功能通路主要分为TGF-β、低氧和炎症并对得到的目标miRNAs进行功能分类。其中有22个miRNAs同时参与了TGF-β、低氧和炎症3个功能通路，评分最高，它们分别是miR-let-7家族、miR-124、miR-130家族、miR-135、miR-144、miR-149、miR-155、miR-16-1、miR-17、miR-181家族、miR-182、miR-200家族、miR-204、miR-205、miR-21、miR-224、miR-27、miR-29家族、miR-301a、miR-31、miR-361和miR-375(Fig 2)。

Fig 2 A subset of miRNAs previously associated with hypoxia, inflammation and TGF-β is predicted to target PH network

图选项

为进一步验证生物信息学分析所得miRNAs，我们构建了低氧诱导PH大鼠模型，探讨miRNAs在肺动脉中的表达变化。结果显示，低氧模型组大鼠的RVSP及mPAP均显著高于对照组(Fig 3，P＜0.01) 。提示：低氧PH模型大鼠构建成功。

随机挑选在EMT和PH中作用均肯定的miRNAs 2个(miR-21、miR-124) 、只在EMT中作用肯定的miRNAs 1个(let-7g)及只在PH中作用肯定的miRNAs 1个(miR-130家族)进行实时荧光定量PCR初步验证。结果发现与对照组相比，低氧诱导PH大鼠肺动脉中let-7g和miR-124表达明显下调，miR-21、miR-130家族表达明显上调(Fig 4，P＜0.05) 。进一步提示通过生物信息学预测得到的let-7g、miR-21、miR-124、miR-130家族可能与低氧PH时EndMT有关。

Fig 3 Homodynamic index for hypoxia-induced PH rats(n=10) **P＜0.01 vs Control

图选项

Fig 4 Change of miRNAs in pulmonary artery in hypoxia induced PH rats(n=8) *P＜0.05, **P＜0.01 vs Control

图选项

对这22个miRNAs进行miRNAs-靶点(PH网络中的)网络构建，以直观观察该22个miRNAs在PH内皮间质转化中可能发挥的整体作用(Fig 5)。

Fig 5 Predicted target network of 22 miRNA groups

图选项

EndMT即在外界刺激及体内微环境的影响下，单层内皮细胞间固有的联系遭到破坏，内皮细胞失去原有的表型，获得狭长的间质样及肌成纤维细胞样表型^[7]。在EndMT过程中，内皮细胞标志蛋白(VE-cadherin, Tie-1/2, VEGFR1/2, PECAM/CD31) 表达明显下调，肌成纤维细胞和平滑肌细胞标志蛋白波形蛋白(vimentin)和α-平滑肌肌动蛋白(α-smooth muscle actin, α-SMA)表达增加。实验证据提示^[8]，EndMT在心血管发育(包括肺动脉发育)、各种血管疾病(如动脉粥样硬化和血管再狭窄时内膜增生)以及伤口愈合中发挥重要作用。有学者发现在PH患者的肺动脉内膜中，同时存在着内皮细胞和多功能间质祖细胞，其中肌成纤维细胞在肺动脉内膜增生中发挥重要作用^[9]；PH患者肺动脉内皮细胞形态发生改变，并且通过标志蛋白的表达检测发现有平滑肌样细胞(即α-SMA标记细胞)和小部分移行细胞(即内皮样细胞标志蛋白和平滑肌样细胞标志蛋白同时出现的细胞)的存在^[10]。提示EndMT可促进α-SMA样细胞在病变血管部位聚集从而参与PH时肺血管重构。但PH时肺血管EndMT发生的机制尚不明确。

自2007年英国药理学家Hopkins^[11]率先提出“网络药理学”概念后，由于该方法成本较低，并可有效地帮助科研工作者从大数据中缩小范围寻找靶标，显示出重要的理论和实际应用价值。目前被迅速应用于众多领域的研究中，以探索miRNAs和基因在人类疾病中的整体作用。早在2012年便有学者运用网络药理学的方法探索了与PH相关的miRNAs，并最终通过实验证实了miR-21在多种PH模型中表达上调并可通过抑制其靶点RhoB改善PH肺血管重构，延缓病情的发展^[12]。本文首次运用网络药理学方法探索与PH时肺血管EndMT相关miRNAs及其相关下游靶点的相互作用机制，并通过实验对部分miRNAs进行实验验证。其中let-7家族属于肿瘤抑制因子，可通过HMGA2抑制多种肿瘤细胞的EMT^[13]；miR-130可促进多种类型的PH肺血管重构^[14]；miR-124及miR-21已确证参与了PH的发生发展，且参与了多种细胞的EndMT/EMT。PH时，miRNA124在肺动脉外膜及肺组织中均表达下调，通过多个靶点(NFATc1、CAMTA1、PTBP1) 抑制PASMCs的增殖和表型转化及肺动脉成纤维细胞的增殖、迁移和炎症^[15-16]。miR-124还可通过多个靶点(CDH2、RHOG)抑制TGF-β或低氧诱导的肺癌细胞^[17]及视网膜色素上皮细胞^[18]EMT。与其他学者结果一致，我们验证时发现低氧PH时let-7g下调而miR-130上调，这从侧面反映了低氧PH时可能发生了EndMT。而miR-124在低氧诱导的肺动脉中表达下调，提示其可能在PH时具有抑制肺血管EndMT的作用。总之，网络药理学方法不仅可通过生物信息学预测从而为后续的研究提供方向和策略，还可通过探究信号分子在疾病中的整体作用而解释表面看似矛盾的结果。

对于本研究的实验验证，尚存在以下不足。首先，我们通过分子生物学实验验证了部分miRNAs在低氧PH大鼠肺动脉中的表达变化，对于其它miRNAs如何改变未做探讨；其次，我们仅在低氧诱导的PH模型中进行了miRNAs验证而未涉及其他动物模型如野百合碱诱导PH模型，而两种模型代表的PH病理过程有差异；再次，肺动脉由肺动脉内皮细胞、平滑肌细胞及成纤维细胞等多种细胞组成，我们仅检测了miRNAs在肺动脉中的表达变化，故这不足以说明其参与了PH时肺血管EndMT。因此，后续还需在多种动物模型及细胞水平进一步验证这些miRNAs在PH时肺血管EndMT的整体作用及机制。

( 致谢: 本文实验在中南大学药学院药理学系心血管药理研究室完成。湖南医药学院药理学教研室吴卫华博士以及中南大学药学院药理学系博士生祝田田、邹小舟、刘汀和硕士生葛晓月参与完成本文实验，并协助数据分析。 )