目的:构建胃癌相关竞争内源性RNA(ceRNA)网络，探索胃癌发生、发展的分子机制。方法:应用生物芯片技术测定胃癌及对应癌旁组织的信使RNA(mRNA)、微小RNA(miRNA)和长链非编码RNA(lncRNA)表达谱，应用R软件中的edgeR包筛选差异表达基因并绘制火山图，基于miRNA－lncRNA在线预测工具lncBase数据库和miRNA靶基因预测数据库预测mRNA、miRNA和lncRNA的相互作用关系，采用Cytoscape软件构建胃癌lncRNA－miRNA－mRNA ceRNA网络，根据cytohubba计算各节点degree得分筛选关键基因(hub基因)。运用注释、可视化和集成发现数据库(DAVID)，对差异表达的mRNA进行基因本体论(GO)功能富集和京都基因与基因组百科全书(KEGG)通路富集分析。在OncoLnc数据库中，选择肿瘤基因组计划(TCGA)数据集的胃癌项目，输入hub基因，以标准化后基因表达量中位数为界值，将样本分为高表达组和低表达组，并进行生存分析。结果:共筛出差异表达的mRNA 766个，miRNA 10个，lncRNA 110个，其中90个mRNA、6个miRNA和4个lncRNA构成ceRNA网络，20个hub基因中有2个与患者预后有关。经TCGA胃癌数据库验证，LncRNA MAP3K20的反义RNA 1(MLK7－AS1)、分泌性磷酸化糖蛋白(SPP1)、硫酸酯酶1(SULF1)hsa－miR－1307－3p等基因在胃癌组织生物芯片中高表达，金属硫蛋白2A(MT2A)、金属硫蛋白1X(MT1X)等基因低表达，与TCGA胃癌数据库一致。生物学功能和通路分析显示，差异表达的mRNA主要在生物学过程(BP)和矿物吸收的通路中显著富集。在TCGA胃癌数据中，碳水化合物磺基转移酶1(CHST1)和miR－183－5p均与患者生存时间有关，CHST1表达量越高患者生存时间越短，而miR－183－5p表达量越高患者生存时间越长。结论:胃癌ceRNA网络的构建有助于深入了解胃癌的分子生物学机制，CHST1和mi－183－5p可作为胃癌预后的预测指标。

目的:研究二甲双胍对2型糖尿病患者微小RNA(miRNA)表达水平的影响，并利用网络药理学进行分析，预测其治疗作用潜在的靶点及途径。方法:筛选本院入院治疗的初发2型糖尿病患者15例，住院期间及出院后均规律服用盐酸二甲双胍片。患者服药6个月后，对比其用药前后血清基质金属蛋白酶(MMP)-9、转化生长因子(TGF)-β1及心肌纤维化相关miRNA的表达水平。对呈现统计学差异表达的miRNA进行基因本体论(GO)和KEGG通路富集分析，并构建差异表达miRNA对应靶基因信使RNA(mRNA)网络图。结果:患者用药后，血清MMP-9水平较用药前明显下降( P<0.05)。二甲双胍可显著下调患者血清人类微小RNA(hsa-miR) 29a-3p、hsa-miR133a-5p、hsa-miR21-5p、hsa-miR30c-5p、hsa-miR1-3p表达水平( P<0.05或 P<0.01)。GO和KEGG分析结果显示，差异表达miRNA主要集中在内质网腔、突触、基膜等细胞组分中，通过胶原分解代谢过程、短时神经元突触可塑性调节、轴突延伸等生物过程来发挥Rho GTP酶(Rho GTPase)结合、参与细胞外基质结构成分等分子功能；其主要富集于糖尿病并发症晚期糖基化终末产物-糖基化终末产物受体(AGE-RAGE)信号通路、肿瘤坏死因子(TNF)信号通路、Wnt信号通路等多种细胞信号转导通路。miRNA-mRNA网络分析结果显示，与差异表达miRNA对应mRNA共230个。 结论:二甲双胍可通过下调相关miRNA表达，参与相关细胞信号通路转导，调控染色质、核酸结合及酶活性过程而发挥其治疗2型糖尿病的作用。

目的:筛选与新生儿支气管肺发育不良(BPD)发病相关的关键转录因子(TF)、微小核糖核酸(miRNA)和信使核糖核酸(mRNA)，并构建了调控网络。方法:从高通量基因表达数据库GEO获取mRNA芯片数据集GSE108756，筛选出的差异表达基因(DEGs)进行京都基因和基因百科全书(KEGG)、基因本体注释(GO)功能富集分析；参照HumanTFDB数据库获取与DEGs对应的TF后利用cytoHubba插件筛选出degree值前10的Hub TF；通过hTFtarget数据库找出Hub TF中7个TF所对应的靶基因；在GSE108756中获取7个TF及其靶基因所对应的探针表达量后做相关性分析，筛选出最终纳入研究的6个Hub TF及所对应的正相关TF-mRNA关系队；利用Targetscan Human数据库获取靶向调控6个Hub TF的miRNAs，并与miRNA芯片数据集GSE166762差异miRNAs取交集后匹配出负向调节的miRNA-TF关系队，最后构建BPD患儿潜在miRNA-TF-mRNA调控网络，进一步筛选出关键miRNA、mRNA。结果:共筛选出201个差异表达基因，KEGG富集结果主要涉及B细胞受体信号通路、造血系统、原发免疫缺陷等作用，GO功能富集分析主要富集在B细胞增殖活化及免疫应答调节细胞表面受体等生物学过程中。与Hub TF(EBF1、TCF4、JUN、BCL11A、SPIB、E2F1)表达谱强关联的正相关下游靶基因有159个，负相关的上游miRNA有120个。BPD患儿发病潜在的关键miRNA-TF-mRNA网络包括has-miR-217-TCF4-Pou2AF1、has-miR-150-5p-E2F1-ADM。结论:本研究构建了新生儿BPD的miRNA-TF-mRNA共调控网络，从转录调控的崭新视角为揭示BPD发病机制奠定理论基础。

目的:通过高通量测序鉴定子痫前期（PE）孕妇胎盘组织中环状RNA（circRNA）的表达谱，构建circRNA-微小RNA（miRNA）-信使RNA（mRNA）相互作用网络，揭示PE发病的相关通路及调控机制。方法:收集2019年11月至2021年6月在四川大学华西第二医院分娩的42例PE孕妇（PE组）和30例正常妊娠孕妇（对照组）的临床资料并采集其胎盘组织。（1）采用高通量测序技术建立其中5对PE组和对照组孕妇胎盘组织中的差异表达circRNA谱。（2）实时荧光定量PCR（qRT-PCR）技术验证其中6个差异表达circRNA在PE组和对照组胎盘组织中的表达水平。（3）使用生物信息学方法预测靶miRNA和分析共表达mRNA，构建竞争性内源RNA（ceRNA）网络；并对差异表达circRNA进行基因本体（GO）功能注释分析和京都基因与基因组百科全书（KEGG）通路分析。（4）采用logistic回归分析、Pearson相关、Kendall tau-b相关分析检验3个差异表达circRNA与PE发生风险、临床特征的相关性。（5）选择circRNA_05393进行后续功能研究，使用小分子干扰RNA、过表达质粒分别敲低或提高滋养细胞系HTR-8/SVneo细胞中circRNA_05393的表达水平，采用穿膜小室（transwell小室）体外实验检测细胞的迁移、侵袭能力，活细胞计数法检测细胞的增殖能力，细胞成管实验检测细胞的成管能力。结果:（1）高通量测序发现了72个差异表达circRNA，其中35个上调，37个下调。（2）qRT-PCR技术检测显示，与对照组比较，PE组孕妇胎盘组织中circRNA_00673（分别为1.306±0.168、2.059±0.242； t=2.356， P=0.021）和circRNA_07796（分别为1.275±0.232、1.954±0.230； t=2.018， P=0.047）的表达水平显著升高，而circRNA_05393（分别为1.846±0.377、0.790±0.094； t=3.138， P=0.002）的表达水平显著降低。（3）circRNA-miRNA-mRNA相互作用网络中包含3个circRNA、8个miRNA和53个mRNA。GO分析显示，生物过程方面主要富集于铁离子稳态、动作电位膜去极化和神经元动作电位等通路，细胞组分方面主要富集于皮质细胞骨架，质膜成分等，分子功能方面主要富集于电压门控钠通道活性、碱性氨基酸跨膜转运蛋白活性等。KEGG信号通路富集分析显示，相互作用网络中的mRNA主要富集于补体和凝血级联反应、甘氨酸、丝氨酸和苏氨酸代谢，p53信号通路，过氧化物酶体增殖物激活受体信号通路。（4）logistic回归分析结果显示，circRNA_05393表达下调是PE发生的危险因素（ OR=0.044，95% CI为0.003~0.596； P=0.019）；相关性分析结果显示，circRNA_05393与PE孕妇的收缩压、舒张压显著相关（ P均<0.05）。（5）敲低或过表达circRNA_05393显著降低或增加HTR-8/SVneo细胞的迁移、侵袭能力（ P均<0.05），但对细胞成管能力、增殖能力无显著影响（ P均>0.05）。 结论:胎盘组织中circRNA表达谱的构建、circRNA-miRNA-mRNA相互作用网络的探索为揭示特定circRNA在PE发生中参与的调控机制提供了可能。抑制circRNA_05393可能通过减少滋养细胞的迁移和侵袭来参与PE的进展。

目的:构建反复种植失败（RIF）相关竞争性内源RNA（ceRNA）调控网络并进行临床样本验证。方法:从高通量基因表达数据库（GEO）得到RIF相关的长链非编码RNA（lncRNA）、微小RNA（miRNA）和信使RNA（mRNA）表达谱数据集，构建lncRNA-miRNA-mRNA的ceRNA调控网络。同时利用加权基因共表达网络分析（WGCNA）探索网络中的枢纽基因。回顾性收集2020—2021年于郑州大学第一附属医院生殖医学中心行体外受精（IVF）/卵胞浆内单精子显微注射（ICSI）助孕的RIF患者和对照组（助孕后至少有1次妊娠史）患者的子宫内膜组织，应用实时荧光定量聚合酶链反应（qRT-PCR）验证RIF相关枢纽基因的mRNA表达水平，并应用Western印迹和免疫组化技术验证血管细胞黏附分子-1（VCAM1）的蛋白表达水平。结果:构建了由32个lncRNA、31个miRNA和88个mRNA组成的RIF相关ceRNA调控网络，并利用WGCNA鉴定出7个RIF相关枢纽基因。将ceRNA网络中的88个mRNA与枢纽基因取交集得到2个RIF相关关键基因：VCAM1和白细胞介素2受体α（IL-2RA）。临床验证中，对照组和RIF组的年龄［ M（ Q1， Q3）］分别为26.50（25.00，34.00）和30.50（25.75，35.25）岁（ P>0.05）。与对照组相比，RIF组子宫内膜中VCAM1的mRNA［ M（ Q1， Q3）］［0.30（0.15，0.42）比0.99（0.69，1.34）， P=0.001］和蛋白表达水平［ M（ Q1， Q3）］［0.44（0.16，1.27）比2.39（1.58，2.58）， P<0.001］均降低。 结论:本研究成功构建了RIF相关ceRNA调控网络，并通过临床样本证实RIF患者子宫内膜组织中VCAM1的表达水平降低。

目的:使用生物信息学方法构建铁死亡相关微小RNA(miRNAs)-信使RNA(mRNAs)的调控网络并探讨其在骨关节炎(OA)中的分子作用机制。方法:从基因表达综合数据库(GEO)下载软骨相关数据集GSE175961和GSE1140007，利用生物信息学方法筛选铁死亡相关miRNAs和mRNAs，并进行富集分析、miRNA-mRNA网络构建以及诊断价值分析。最后采用荧光实时定量聚合酶链反应(qRT-PCR)验证关键基因表达，并采用独立样本 t检验比较其差异。 结果:GSE175961数据集中共筛选获得了15个差异表达miRNAs。使用TargetScan和miRDB数据库分别预测差异表达miRNAs的下游靶基因，共获得了5 115个共存的靶基因。与GSE1140007数据集中19个差异表达铁死亡基因取交集后共获得了4个下游铁死亡基因(PARP15、KLF2、CDKN1A、PDK4)。基于上述研究结果构建了由15个差异表达的miRNAs和4个下游铁死亡基因组成的miRNA-mRNA调控网络。进一步富集分析表明，4个下游铁死亡基因与多种生物学功能和信号途径有关，如脂肪酸代谢、细胞周期调控、FOXO信号途径以及P53信号途径等。受试者工作曲线分析表明，4个下游铁死亡基因具有良好的诊断价值。此外，使用外部数据集证实，CDKN1A和PDK4在OA软骨组织中表达降低。而列线图表明，OA发生率随CDKN1A和PDK4表达降低而升高。qRT-PCR证实，正常组和OA组中，CDKN1A相对表达量分别为1.322±0.890、0.035±0.031；PDK4相对表达量分别为1.085±0.567、0.393±0.180。结论:本研究构建的铁死亡相关miRNA-mRNA调控网络与软骨退变密切相关。

目的:基于miRNA-mRNA调控网络进行生物信息学分析，探讨骨关节炎病变的相关分子机制。方法:从GEO数据库下载人血清样本miRNA表达数据集，通过R语言limma包获取差异表达miRNA,采用miRwalk 2.0版数据库预测其对应靶基因(mRNA)，并构建miRNA-mRNA调控网络。对靶基因进行功能GO分析和KEGG信号通路分析，构建靶基因所编码的蛋白质相互作用网络(PPI)，从中筛选出骨关节炎病变的核心基因。结果:共筛选出7个差异表达的miRNA(表达均为下调)和900个mRNA，这些基因主要涉及蛋白结合、DNA结合、转录等生理过程，参与Cell cycle、p53、Neurotrophin、PI3K-Akt等信号通路。蛋白相互作用分析表明MAPK1、TP53、MAPK14、CCND1、EP300、POLR2E、POLR2F、ABL1、RAC1、SKIV2L2为该调控网络的核心靶基因。结论:OA的发生和发展涉及多个的miRNA、靶基因和作用途径，而通过构建骨关节炎相关miRNA-mRNA调控网络，可为找出骨关节炎病的分子机制和今后临床上诊疗提供新的思路。

目的:探讨环状RNA（circRNA）在RA PBMCs表达谱的差异及临床意义。方法:收集4例RA患者（T组）及4名健康体检者（C组）静脉血，运用Arraystar circRNA微阵列芯片检测2组PBMCs中circRNA表达谱，应用相关数据库进行生物信息学分析。采用实时荧光定量PCR（RT-qPCR）验证表达失调circRNA的表达情况，并对目标circRNA构建circRNA-微小RNA（miRNA）-信使RNA（mRNA）调控网络（即ceRNA网络）。建立受试者工作特征曲线（ROC曲线），对潜在诊断价值进行分析。采用SPSS Statistics 23.0和Graphpad Prism 8.0分析数据，采用两独立样本 t检验分析2组间差异。 结果:①微阵列芯片结果显示，与C组比较，T组PBMCs中异常表达的circRNA共有399个，其中149个上调，250个下调。②生物信息学分析：circRNA与miRNA结合位点预测提示RA中差异表达的circRNA可能通过靶向结合miR-140-5p、miR-338-5p和miR-9-5p等分子影响炎症反应；基因本体论分析发现差异表达的circRNA主要在细胞质、蛋白质结合等方面；京都基因与基因组百科全书通路分析发现涉及的通路多与人体免疫调节相关。③多样本RT-qPCR验证结果显示T组中hsa_circRNA_009012的相对表达量显著高于C组（ t=-4.417， P<0.01），hsa_circRNA_101328的表达水平显著低于C组（ t=-1.042， P<0.01），hsa_cir-cRNA_058230的表达无明显改变（ t=4.691， P>0.05）。④ ROC曲线分析提示hsa_circRNA_009012具有诊断RA的潜在价值（ROC曲线下面积=0.96）。 结论:circRNA在RA患者PBMCs中表达失调，可能参与RA发生发展的调控，其中hsa_circRNA_009012有望成为RA诊疗的新型生物学标志物。

目的:基于生物信息学方法挖掘胆道闭锁(BA)差异表达基因(DEGs)并初步构建微小RNA(miRNA)-信使RNA(mRNA)调控网络，探索其在BA中的分子调控机制。方法:从基因表达综合数据库(GEO)下载数据集GSE46960，利用GEO2R在线工具分析BA和正常肝组织之间的DEGs，利用注释、可视化和富集发现数据库(DAVID6.8)对DEGs进行富集和功能注释。基于蛋白质互作数据库(STRING11.0)和Cytoscape_v3.7.1软件构建蛋白质相互作用(PPI)网络，分析其关键基因。通过人类miRNA疾病数据库(HMDD_V3.2)查询与BA相关且经过验证的miRNA，并利用miRNA目标预测数据库(miRDB)预测其靶基因。将预测的靶基因和GSE46960中的差异表达基因求得交集，获得miRNA、mRNA相互作用关系。利用Cytoscape软件构建miRNA-mRNA调控网络图。结果:从GSE46960数据集中共获得565个DEGs，其中352个为上调基因，213个为下调基因。其中上调的DEGs显著富集于细胞外基质受体相互作用、黏着斑激酶、阿米巴病通路、磷酸酰肌醇3激酶/蛋白激酶B(PI3K/Akt)信号通路中，下调的DEGs显著富集于新陈代谢、抗生素的生物合成、视黄醇代谢通路中。根据PPI网络，筛选出10个关键基因。在DEGs基础上成功构建了BA miRNA-mRNA分子调控网络。结论:通过生物信息学筛选出的DEGs及构建的miRNA-mRNA调控网络有助于全面了解BA发生的分子机制，为探索联合miRNA及DEGs作为临床诊断及治疗的分子靶点提供了新方向。

目的:通过构建内源竞争性RNA(ceRNA)网络为探究脑动脉瘤破裂分子机制提供线索。方法:从基因表达综合数据库(GEO)下载含有21个破裂脑动脉瘤和21个未破裂脑动脉瘤组织样本的数据集，利用加权基因共表达网络和差异基因表达分析方法分别获得与脑动脉瘤破裂相关的信使RNA(mRNA)和长链非编码RNA(lncRNAs)，并基于miRWalk2.0和miRcode数据库及文献证据分别预测和筛选与上述mRNA和lncRNA相互作用的微小RNA(miRNAs)；最后，通过分析mRNA-miRNA和lncRNA-miRNA相互关系构建脑动脉瘤破裂相关的ceRNA网络。结果:通过共表达网络和差异基因分析共获得与脑动脉瘤破裂相关的470个mRNAs和78个lncRNAs，数据库预测到49个miRNAs可与上述mRNAs和lncRNAs结合，结合文献证据共筛选出13个miRNAs、7个lncRNAs和73个mRNAs构建ceRNA网络。结论:通过生物信息学分析构建了一个包含13个miRNAs、7个lncRNA和73个mRNAs的ceRNA网络，为探究脑动脉瘤破裂具体分子机制提供了线索。