论文部分内容阅读
研究背景及目的:氧化应激损伤是心血管疾病的一个重要的致病因素。缺血缺氧导致心肌细胞产生大量的活性氧自由基（ROS）。而活性氧在控制血管内皮功能、血管张力以及诱导心肌细胞肥大、增殖、凋亡、迁移等病理生理过程中有着重要作用。鼠李素是具有较高的生物活性的黄酮类化合物之一,其对心肌氧化应激损伤具有保护作用,但是关于其发挥作用的分子生物学机制的研究较少,因此本研究的主要目的是采用RNA-Seq分析技术和多种生物信息学方法来揭示鼠李素保护心肌氧化应激损伤的分子机制,为临床治疗提供新的作用靶点。方法:通过H2O2诱导大鼠心肌细胞H9C2,建立损伤模型。H2O2处理H9C2细胞作为对照组,H2O2处理+鼠李素给药为处理组,提取RNA并进行转录组测序,并根据测序结果进行以下分析:（1）筛选差异表达基因与构建相互作用网络。采用NOISeq数据分析包计算差异表达基因,设置筛选条件为:P value<0.05及|log2FC|>2。基于Gene ontology数据库以及KEGG通路数据库对差异表达基因进行功能以及通路富集分析。基于String数据库信息,建立蛋白互作（PPI）网络,设置条件:combined score不低于0.9。基于Cytosape插件MCODE进行子网络模块挖掘,阈值设置:Degree Cut off=2,Node Score Cutoff=0.2,K-Core=2,Max.Depth=100。（2）差异表达基因数据的深度挖掘。基于miRWalk2.0工具检索差异表达基因的调控miRNA,并利用R软件包clusterprofiler对miRNA进行KEGG功能富集分析。基于miRBase,miranda工具筛选score>180的miRNA-lncRNA关系对。ceRNA以及miRNA-TF-gene网络中差异基因的GO-BP功能和KEGG富集运用DAVID工具来分析。运用WebGestalt网站对差异表达基因转录因子进行预测。基于DGIdb药物预测受miRNA和TF调控的关键差异基因的Drug-gene互作关系对。基于软件Cytoscape构建miRNA-gene,miRNA-lncRNA,ceRNA,TF-gene,miRNA-TF-gene以及Drug-gene互作关系网络。（3）差异表达基因的qPCR验证。基于生物信息学预测出的差异表达基因,筛选上调表达基因（NOTCH3,NUMBL,BCL2L11）和下调表达基因（PPP6C,MAPK1）,通过qPCR方法检测其在大鼠心肌细胞H9C2对照组（H2O2处理）和实验组（H2O2处理+鼠李素给药）中的表达差异情况。结果:（1）通过比较对照组（H2O2处理）和实验组（H2O2处理+鼠李素给药）的表达数据,共获得1452个差异表达基因,其中对应上调和下调基因分别为703个和749个。GO功能富集分析结果显示,差异表达基因涉及的生物进程主要包括细胞或大分子代谢、细胞周期或细胞分裂、代谢或信号通路的调控、信号通路等。KEGG通路富集分析主要涉及的关键信号通路包括细胞粘着连接,Notch信号通路、氧化磷酸化、蛋白水解等。PPI网络共包括609个节点,并得到了2个score>6的子网络模块,其中模块1涵盖了17个节点,并富集在mRNA降解、翻译、蛋白质定位相关的功能上;而模块2涵盖15个节点,主要富集于细胞分裂过程相关功能上。（2）通过对差异表达基因做聚类热图分析显示,实验组和对照组之间差异明显。差异表达基因调控miRNA的预测共获得189个miRNA以及577对miRNA-gene调控关系对。其中较为关键的miRNA主要包括miRNA-349,miR-211-5p以及miR-17-5p等。KEGG通路富集显示预测得到的miRNA显著富集到75个通路中,其中AMPK及PI3K-Akt信号通路富集的miRNA最多。通过对miRNA-gene调控关系作图,共获得120个节点和216个互作关系对。基于已有的miRNA对其调控的lncRNA进行预测共获得475个miRNA-lncRNA关系对。ceRNA网络中共有74个miRNA,90个mRNA以及60个lncRNA,其GO功能和KEGG通路富集分析结果显示,主要涉及的生物学过程包括RNA聚合酶II启动子转录的正调控（GO:0045944）,凋亡负调节过程（GO:0043066）以及转录正调控（GO:004589）等过程;主要涉及的关键信号通路包括肾上腺素以及Notch等信号通路。此外,对差异表达基因相关的调控转录因子进行预测,共获得9个转录因子（PAX4,ELK1,ETS2,NRF1,GABP,E4F1,CACBINDINGPROTEIN,AP2ALPHA以及AP2GAMMA）。进一步地,通过构建miRNA-TF-gene调控网络整合得到241个调控关系对,且功能富集和通路富集分别富集到145个GO-BP结果和52个KEGG通路,富集结果与ceRNA网络基本一致。最后,我们对drug-gene互作关系进行预测并构建网络,结合富集生物学过程分析的与药物反应相关的基因分析发现基因BCL2L1,DNMT3A,HDAC2,STAT3等基因为重要的药物小分子基因。（3）qPCR结果显示,与对照组相比鼠李素处理组中的NOTCH3、MAPK1基因表达水平显著性下调（P<0.05）,BCL2L11基因表达显著性上调（P<0.05）,而上调基因BCL2L11与下调基因MAPK1表达趋势与前期预测结果相符,但差异不显著（P>0.05）。结论:总的来说,本研究采用转录组测序技术对鼠李素对心肌氧化应激损伤保护作用的分子机制进行了研究。通过对心肌细胞鼠李素处理组以及对照组分析获得关键差异表达基因包括COX6A2,NDUFA2,UBA52,CDK1,MAPK1,Klf9以及NAA15;重要的miRNA包括miR-211-5p以及miR-17-5p;重要的转录因子包括ELK1及PAX4等;重要的信号通路包括Notch,AMPK以及PI3K/Akt信号通路;重要的药物小分子基因包括BCL2L1,DNMT3A,HDAC2以及STAT3。同时,该研究构建了一系列调控网络关系图,系统阐述了鼠李素对心肌氧化应激损伤保护作用过程中不同作用因子之间的相互作用关系,这也将为今后的研究提供重要的依据。