空间与加权基因共表达分析揭示跨种属的缺血性心力衰竭机制

[目的]从空间与基因共表达网络的角度揭示跨种属的缺血性心力衰竭机制。[方法]从美国国立生物技术信息中心基因表达数据库(NCBI-GEO)中检索获得GSE210374和GSE57338高通量测序数据,利用R语言软件程序包分析筛选在心肌梗死大鼠不同心肌区域中差异表达基因(DEG)以及缺血性心力衰竭患者与健康对照者的心肌样本间的DEG,分析共同基因的分区域表达情况。利用加权基因共表达网络分析(WGCNA)筛选心肌梗死相关的基因并进行富集分析,构建蛋白质-蛋白质相互作用网络(PPI)筛选核心基因(HG)。[结果]在心肌梗死大鼠和正常对照组中共筛选出4835个差异基因,在缺血性心力衰竭患者和正常对照样本共筛选出51个差异基因,揭示心肌梗死后左心室心肌梗死区(I区)、边缘区(BZ区)及远端区(R区)的代表性基因集。空间表达分析发现两个种属样本:在每个心肌分区均有20个共表达基因,其中16个在3个分区均有表达,在I区、BZ区和R区特异表达的基因数分别为2、0和2个。富集分析显示:各分区的共表达基因功能各异,I、BZ区与胶原纤维组装、应激诱导的心肌细胞肥大、下调c-Jun氨基末端激酶(JNK信号)与细胞增殖等功能以及补体信号通路相关;而I、R区则富集于Wnt和胶原结合;作为非缺血的远端R区,共表达基因显著富集于细胞外基质的抗压性、细胞溶解以及抑制T细胞增殖等功能。此外值得注意的是:3个分区的共表达基因的产物大部分位于细胞外空间和细胞外基质当中,提示可能存在活化的细胞分泌与互作调控。进一步PPI分析提示无孢蛋白(ASPN)、骨诱导因子(OGN)和ⅩⅣ型胶原蛋白α链(COL14A1)基因可能是前述机制的核心基因。[结论]大鼠和人类缺血性心力衰竭的共性机制涉及补体与凝血级联信号、Wnt等多种信号通路;可能与细胞外基质、外泌体介导的细胞凋亡密切相关;ASPN、OGN和COL14A1可能是其中的核心基因。本工作有望为缺血性心力衰竭相关转化研究中的干预靶点与路径选择提供空间与通路参考。

基于血清药物化学与加权基因共表达网络分析探究葛兰心宁胶囊治疗冠心病的作用

目的探究葛兰心宁胶囊治疗冠状动脉粥样硬化性心脏病(CAD)可能调控的关键靶点基因。方法HPLC法对入血成分进行鉴定,同时借助GSE20680数据集,基于P<0.05标准,筛选出15549个差异表达基因。利用加权基因共表达网络分析,鉴定与疾病相关的模块基因。确定入血成分可能调控的关键靶点基因,进一步取交集,对其靶点进行富集分析,探索关键靶点的功能和参与的通路,并完成其与相关成分的分子对接。通过动物实验验证加权基因共表达网络分析结果。结果鉴定出21种入血成分。GO分析显示,交集基因主要参与凋亡、炎症反应的调控。KEGG分析显示,入血成分治疗CAD主要参与NF-κB信号通路、糖尿病心肌病等相关通路。分子对接实验验证了同时调控关键靶点PTGS2、PLAU 4种入血成分的相互作用。动物实验显示,葛兰心宁胶囊能够降低小鼠主动脉窦的脂质蓄积,缩小斑块面积,减少PTGS2、PLAU表达。结论本研究揭示了葛兰心宁胶囊入血成分及其治疗CAD可能调控的关键靶点基因,为其进一步临床应用提供了重要依据。

基于加权基因共表达网络和癌症基因组图谱临床数据分析并鉴定肝细胞癌的Hub基因研究

背景 肝细胞癌(HCC)是全球常见的癌症相关死亡的第三大原因,约占所有原发性肝癌病例的90%,其复发率和死亡率较高,目前发生的分子机制仍不清楚。目的 探索HCC潜在的分子机制,发掘新的生物标志物。方法 从TCGA数据库下载RNA-seq表达数据和临床相关信息,通过差异基因表达分析正常肝脏组织与HCC组织的差异基因;对差异表达基因进行富集分析;基于TCGA中HCC的基因表达数据概况,使用WGCNA R包建立共表达网络,进行加权基因共表达网络分析(WGCNA),选择具有临床意义的模块,并筛选候选Hub基因;进一步分析候选Hub基因在HCC组织和正常肝脏组织显著差异表达、与HCC患者总体生存期和无病生存期是否显著相关,最终确定Hub基因;通过人类蛋白质图谱数据库对Hub基因蛋白表达进行验证。结果 本研究的基因表达数据来自50个正常肝脏组织样本和373个HCC组织样本。通过差异基因表达分析发现7 230个在HCC和正常肝脏组织之间差异表达的基因(HCC中3 691个上调基因和3 539个下调基因)。富集分析表明,上调的差异表达基因主要参与细胞周期调控和有丝分裂过程;下调的差异表达基因主要参与小分子代谢和有机酸代谢等过程。WGCNA确定了19个与HCC患者临床特征相关基因模块,通过分析模块与临床特征之间的关系,筛选出青色模块和紫色模块。青色模块基因中同时与患者总生存期和无病生存期强烈相关的前两个基因为VPS45和FAM189B;紫色模块基因中同时与患者总生存期和无病生存期强烈相关的前两个基因分别为CLEC1B和FCN3,因此将VPS45、FAM189B、CLEC1B和FCN3确定为最终的Hub基因。人类蛋白质图谱数据库免疫组织化学染色显示:VPS45和FAM189B在HCC组织中的表达高于正常肝脏组织,FCN3在HCC组织中的表达低于正常肝脏组织,CLEC1B在HCC组织和正常肝脏组织中表达差异不明显。结论 初步确定VPS45、FAM189B、CLEC1B和FCN3可能是HCC的新型潜在生物标志物,这些Hub基因可能为HCC的靶向治疗提供理论基础。

应用加权基因共表达网络分析识别肝细胞癌发生和进展过程中关键通路和基因作用研究

目的探讨肝细胞癌(HCC)发生进展过程中功能富集通路和关键基因表达。方法从基因表达数据库(GEO)下载HBV感染不同阶段和正常对照肝脏转录组数据,构建基因网络并使用加权基因共表达网络分析(WGCNA)将基因分类为不同的模块,对相关模块中的基因进行富集分析。应用GEO数据集进一步验证重要基因水平。结果共6145个组间差异基因参与构建加权基因共表达网络,被分为9个模块,进一步描绘了从肝脏早期病变到肿瘤的进化轨迹,从正常组织到癌组织过程中的细胞增殖、DNA损伤修复和细胞衰老相关通路线性激活;随疾病进展,脂质代谢和凝血等肝脏功能相关通路逐渐受到抑制;在肿瘤发生前,慢性炎症期免疫相关通路被激活,后期逐渐趋向于抑制状态;共鉴定出3个重要衰老相关基因,即CCNA2、UBE2C和ANAPC1,并在外部数据集验证了这3个基因水平变化;进一步分析显示上述3个基因水平与肝癌患者不良预后密切相关(P<0.05)。结论通过生物信息学分析,我们初步确定了肝癌发生和进展过程中潜在途径和重要参与基因,为诊断和治疗干预提供了潜在的靶标。

加权基因共表达网络分析联合差异表达基因分析法鉴定重度抑郁症关键基因

目的 应用加权基因共表达网络分析(WGCNA)联合差异表达基因(DEG)分析法鉴定重度抑郁症(MDD)的关键基因。方法 从基因表达数据库中获取MDD患者和健康志愿者死后脑组织mRNA微阵列数据。应用WGCNA和DEG分析法筛选与MDD最显著相关的模块和DEG,对WGCNA模块基因和DEG进行基因本体论富集分析和京都基因与基因组百科全书通路富集分析,取WGCNA模块基因与DEG的交集,绘制受试者工作特征曲线评估交集基因诊断MDD的能力。结果 最终获得10个基因(FLJ20021、FLJ30058、NNAT、MRPS12、KCNK4、TROVE2、MESP1、MIF、TMEM93、SYNGR1),均具有良好诊断MDD的潜力。结论 FLJ20021、FLJ30058、NNAT、MRPS12、KCNK4、TROVE2、MESP1、MIF、TMEM93、SYNGR1可作为辅助诊断MDD的生物标志物。

基于加权基因共表达网络分析鉴定PE患者的潜在生物标志物

目的:本研究旨在通过对子痫前期(PE)患者基因数据集进行加权基因共表达网络分析(WGCNA),寻找PE的分子标志物。方法:从GEO数据库下载GSE10588表达数据集,该数据集包括17例PE患者和26例对照者的胎盘组织转录组数据。使用差异表达分析及WGCNA筛选与PE相关的重点差异表达基因(DEGs),利用基因本体论(GO)和京都基因与基因组百科全书(KEGG)通路分析来探索重点基因的生物学作用。构建蛋白-蛋白相互作用(PPI)网络筛选关键枢纽基因,并通过风险预测模型和受试者工作特征(ROC)曲线评价这些枢纽基因对PE的预测效果。结果:共获得138个DEGs,其中包括100个上调和38个下调DEGs。通过WGCNA确定了23个与PE发病最相关的重点DEGs。GO和KEGG通路分析显示,这些重点DEGs在调节促性腺激素分泌、白细胞分化、造血功能、内分泌功能及B细胞的激活等生物学作用中发挥重要作用。在PPI网络中,卵泡抑素相关基因3(FSTL3)、酪氨酸激酶受体1(FLT1)、抑制素亚单位A(INHBA)、N-myc下游调节基因1(NDRG1)和瘦素(LEP)等成为节点最多的基因。风险预测和ROC曲线表明,FSTL3和LEP具有较高的风险得分和诊断价值。结论:FSTL3和LEP可以作为预测和诊断PE的生物标志物。

基于加权基因共表达网络分析筛选肝细胞癌预后相关生物标志物及潜在治疗药物

目的 筛选肝细胞癌(hepatocellular carcinoma,HCC)预后相关生物标志物,从分子水平探讨HCC可能的发病机制,并预测对HCC具有潜在治疗作用的候选药物。方法 下载TCGA数据库和GEO数据库中HCC的转录组数据及其临床记录信息。分别对其基因表达谱数据进行加权基因共表达网络分析(weighted gene co-expression network analysis,WGCNA)和差异表达基因分析,选取两个数据集中与疾病正相关性最高的模块基因与差异表达基因取交集作为关键基因。利用GO和KEGG对关键基因进行功能富集分析。利用STRING数据库构建PPI网络,使用Cytoscape软件对关键基因进行相关性分析,筛选核心基因。使用R语言程序包K-M进行生存分析明确核心基因与HCC患者预后关系。利用在线数据库DGIdb、DREMIT联合进行HCC潜在治疗药物的筛选,依据药物靶点匹配数量、优选得分、特异性得分对预测结果进行排序,选择排名靠前的药物作为可能的候选治疗药物。结果 最终共得到64个关键基因,主要富集于细胞周期与分裂、DNA复制和损伤修复、病毒感染、P53信号通路等。PPI分析发现,CDC20、KIF2C、CCNB2、KIF20A、CCNA2、TOP2A、UBE2C、NUSAP1、AURKA和TRX2为核心基因。K-M生存分析显示,CDC20、KIF20A、CCNA2、TOP2A与HCC的预后显著相关。DGIdb、DREMI联合筛选对HCC有潜在疗效的候选药物,其中索拉非尼、多维替尼、氟尿嘧啶、米托蒽醌等在DGIdb、DREMI中皆排名靠前。结论 CDC20、KIF20A、CCNA2、TOP2A可能是HCC预后的潜在生物标志物,但仍需进一步临床研究加以验证。HCC的发生发展可能与病毒感染、细胞周期与分裂,DNA复制和损伤修复、P53信号通路有关。索拉非尼、多韦替尼、氟尿嘧啶、米托蒽醌等可能作为HCC治疗的潜在临床候选药物。

通过加权基因共表达网络和LASSO分析识别儿童呼吸道合胞病毒的诊断基因

目的:运用生物信息学方法挖掘儿童呼吸道合胞病毒(RSV)的诊断标志物,为RSV的靶向治疗提供依据。方法:从GEO中获得RSV的基因表达数据集。运用差异表达分析、免疫浸润分析、WGCNA和LASSO确定RSV的关键生物标志物。同时,绘制ROC曲线对关键基因的诊断性能进行评估。结果:共识别出592个差异基因和10种差异浸润的免疫细胞,例如中性粒细胞。WGCNA和LASSO分析筛选出7个关键的生物标志物,包括ADM、LTF、UBE2J1、PDK4、S100A4、BST1、CAMP。此外,ROC曲线分析发现这些基因具有较高的诊断价值。结论:筛选的7个基因可作为RSV的诊断生物标志物。

基于加权基因共表达网络分析探索射血分数保留的心力衰竭中微RNA功能模块和分子调控网络

目的通过生物信息学方法探索射血分数保留的心力衰竭(HFpEF)潜在的生物标志物和治疗靶点。方法使用NCBI-GEO数据库检索获得GSE53437高通量测序数据进行分析。通过加权基因共表达网络分析(WGCNA)获取与HFpEF密切相关的基因模块和核心微RNA(miRNA),利用miRNet数据库预测核心miRNA的靶向mRNA。应用R语言“clusterProfiler”程序包对靶基因进行KEGG通路和基因本体论(GO)富集分析。最后采用STRING数据库和Cytoscape软件构建蛋白质-蛋白质相互作用(PPI)网络和miRNA-mRNA基因网络。结果在基因模块与疾病相关性分析中,棕色模块与HFpEF的相关系数(r)绝对值为0.37,P值为3e-4,具有最高的正相关性。各基因模块的基因显著性(GS)分布显示,棕色模块具有最高的GS值。模块内分析显示,棕色模块的模块相关性(MM)与GS之间密切相关(r=0.61,P值为7e-52)。以MM>0.85且GS>0.40为条件,共筛选出17个与HFpEF高度相关的核心miRNA。应用miRNet数据库预测到1578个靶向mRNA。GO注释分析结果显示,该组基因参与的生物学过程包括造血调控、T细胞激活和细胞组分大小调控等,涉及的细胞成分包括转录调控复合物、运输囊泡、早期内体等的构成,分子功能主要富集在DNA结合转录抑制因子活性/RNA聚合酶Ⅱ特异性、泛素蛋白连接酶结合和生长因子活性等方面。KEGG通路分析结果显示,该组基因主要富集于Hippo信号通路、ErbB信号通路、TNF-α信号通路、Ras信号通路等。miRNA-mRNA基因网络分析显示,miRNA-3188与UBA52、HDAC2、PSMF1和SUFU相互作用,miRNA-3909与HDAC2、PSMF1、SUFU和RELA相互作用,miRNA-762与PSMF1、SUFU和RELA相互作用,miRNA548b-3p与UBA52、HDAC2和PSMF1相互作用,miRNA-3198与UBA52、PSMF1、SUFU和RELA相互作用。结论本研究可能有助于阐释HFpEF的分子机制,为识别HFpEF的生物标志物及潜在治疗靶点奠定新的理论依据。

基于加权基因共表达网络分析和机器学习筛选脓毒症免疫抑制特征基因及其靶向中药活性成分预测

目的本研究通过加权基因共表达网络分析(weighted gene co-expression network analysis,WGCNA)和多种机器学习方法识别脓毒症免疫抑制特征基因,并预测其靶向中药活性成分。方法从美国国家生物技术信息中心基因表达综合数据库中下载GSE182522数据集,使用Bioconductor的“limma”软件包提取出差异表达基因(differentially expressed genes,DEGs),对DEGs进行京都基因与基因组百科全书(Kyoto encyclopediaofgenesandgenomes,KEGG)通路和基因本体论(geneontology,GO)富集分析,同时采用WGCNA,以及最小绝对收缩和选择算子(least absolute shrinkage and selection operator,LASSO)回归、支持向量机递归特征消除(support vector machine-recursive feature elimination,SVM-RFE)、弹性网络回归分析3种机器学习法筛选出脓毒症免疫抑制特征基因。随后,构建受试者操作特征曲线(receiver operating characteristic curve,ROC curve,简称ROC曲线)验证特征基因的诊断效能,并绘制箱线图展示特征基因的表达模式。选取SPF级BALB/c小鼠39只,采用随机数字表法将小鼠随机分为脓毒症免疫正常组(n=20)和脓毒症免疫抑制组(n=19)。盲肠结扎穿孔术构建脓毒症小鼠模型,在术后12h取脓毒症免疫正常组小鼠的外周血,术后24 h取脓毒症免疫抑制组小鼠的外周血。采用逆转录定量聚合酶链反应(reverse transcription quantitative polymerase chain reaction,RT-qPCR)检测两组小鼠TRBV7-2的表达情况,验证TRBV7-2的诊断效能。最后,通过分子对接技术预测特征基因潜在靶点中药活性成分。结果共筛选出445个DEGs,其中上调基因173个,下调基因272个。DEGs主要富集于氮代谢、胆汁分泌、胃酸分泌3条信号通路,分子功能的正调控、细胞对含氧化合物的反应、细胞对含氮化合物的反应、对肽的响应、葡萄糖输入的调控等生物过程,质膜蛋白复合体、T细胞受体复合物、细胞侧膜等细胞组分,核苷三磷酸酶调节活性、GTPase激活剂活性、外源蛋白结合等分子功能。通过WGCNA共筛选出56个枢纽基因,3种机器学习交集得到1个特征基因:TRBV7-2。ROC曲线分析得出TRBV7-2的AUG=0.72,具有较高的临床诊断价值,其表达量在脓毒症免疫抑制患者样本中显著下调。RT-qPCR结果显示,脓毒症免疫抑制组小鼠血液中TRBV7-2基因表达低于脓毒症免疫正常组。通过分子对接技术得出小檗碱、黄芩苷、苍术素等10个潜在靶向中药活性成分。结论TRBV7-2可能作为脓毒症免疫抑制的诊断生物标志物,小檗碱等10味靶向中药活性成分可能成为降低脓毒症致死率的立足点。

基于加权基因共表达网络分析及机器学习构建铁死亡相关基因在多囊卵巢综合征的疾病模型及相关中药预测

目的:利用机器学习分析铁死亡基因与多囊卵巢综合征(PCOS)的相关性,构建诊断及预后风险模型,探讨PCOS的潜在铁死亡机制,并预测相关中药。方法:多芯片标准化合并GEO数据库的卵巢组织芯片,采用加权基因共表达网络分析(WGCNA)算法筛选与PCOS相关的铁死亡特征基因,对特征基因进行功能富集分析,利用多种机器学习进行构建并筛选疾病模型,提取最优的疾病特征基因,并构建风险模型,分析其与免疫浸润细胞及功能之间的相关性,对铁死亡基因进行中药及小分子化合物预测,并对小分子化合物与铁死亡靶点进行分子对接验证。结果:共提取26个铁死亡与PCOS的强相关交集靶基因,主要在铁死亡通路中的Xc-系统和铁代谢通路上发挥作用。通过WGCNA结合机器学习筛选以HMOX1、IFNA1、LPCAT3、FOXO4和PARP3为特征基因的RF模型有较优的模拟特性,并通过绘制列线图来预测各基因表达的患病风险性。上述模型基因与PCOS的发病呈正相关性。免疫浸润分析提示多种免疫细胞在两组之间存在显著浸润差异,预测的中药中丹参、三七可能是其潜在的分子药物来源,分子对接中丹酚酸B、葛根素与黄芪甲苷A的匹配程度最高。结论:以HMOX1、IFNA1、LPCAT3、FOXO4和PARP3为特征基因的铁死亡相关基因模型可为早期识别和治疗POCS提供新思路。铁死亡相关基因可通过干预多种途径在PCOS中发挥重要作用,中医药在干预PCOS中可能通过调节铁稳态代谢来介导铁死亡途径。

用权重基因共表达网络分析识别心脏重构关键节点基因

目的采用权重基因共表达网络分析方法(WGCNA)挖掘心肌梗死后心脏重构的关键节点基因,并研究其与ACE1和ACE2的关系。方法从NCBI的GEO下载2个心肌梗死后心脏重构的全基因组表达数据GSE7487和GSE738;数据初处理后,用WGCNA构建基因共表达网络,识别与心脏重构相关的模块与关键节点基因,分析关键节点基因与ACE1和ACE2的关联性;并在心肌梗死后心脏重构大鼠模型中验证它们的关系。结果分析发现在GSE7487,17个模块中有6个模块与心脏重构相关,模块基因富集于16条KEGG信号通路。在GSE738,5个模块与心脏重构相关,模块基因富集于15条KEGG信号通路,其中有10条信号通路与第一组数据结果相同,这些信号通路涉及心肌肥厚病理、氧化磷酸化、代谢等。进一步利用模块内连通性和基因重要性找到了一些心脏重构的关键调控基因,如钙依赖磷酸酶调节子(RCAN1)。RCAN1表达与ACE1表达高度相关,但与ACE2不相关。在动物模型中验证结果与上述结果一致。结论权重基因共表达网络分析方法是一个高效的系统生物学方法,应用本方法发现了心脏重构的关键节点基因,其中RCAN1可能影响ACE1-ACE2在肾素血管紧张素系统中的平衡。

骨关节炎中枢纽基因及在免疫浸润中作用的生物信息学分析与鉴定

背景:低度慢性炎症被认为在骨关节炎的发病机制中起着核心作用,但是具体的分子机制目前仍不清楚。目的:筛选和探讨骨关节炎中的潜在枢纽基因及免疫细胞的浸润情况。方法:将来自GPL570平台的GSE206848数据集与来自GPL96平台的GSE55235和GSE55457数据集合形成原始数据集。使用加权基因共表达网络分析去除离群样本,随后鉴定差异表达基因,并对差异表达基因进行功能富集分析。此外,构建差异表达基因的蛋白质-蛋白质相互作用网络,并使用Cytoscape软件中的两种不同算法筛选枢纽基因。利用CIBERSORT算法评估骨关节炎样本与正常对照之间免疫细胞浸润比例的差异。通过对收集到的骨关节炎患者滑膜组织样本进行定量反转录聚合酶链反应(RT-qPCR)实验,并结合来自GPL96测序平台的GSE12021数据集作为独立数据集,验证枢纽基因对骨关节炎的诊断能效。结果与结论:①在去除5个离群样本后,共鉴定出340个差异表达基因,其中包括159个上调基因和181个下调基因。通过加权基因共表达网络分析和Cytoscape共获得了6个枢纽基因。②CIBERSORT分析显示,骨关节炎组织中多种类型免疫细胞浸润比例与正常组织相比存在差异,6个枢纽基因的表达水平与骨关节炎中多种免疫细胞的相对比例密切相关。③RT-qPCR结果表明,6个基因的相对表达水平相对于正常组织呈下调趋势,但是NFKBIA和PTGS2的表达差异不显著(P>0.05);原始和外部数据集中的ROC曲线表明,6个枢纽基因对骨关节炎具有较强的诊断能力(AUC>0.8)。④结论:最终确定了4个枢纽基因,分别为CDKN1A、MYC、CXC趋化因子配体2和血管内皮生长因子A,可能通过介导免疫应答和炎症反应成为未来治疗骨关节炎的分子靶点。

缺乏运动降低成年人胰岛素敏感性的加权基因共表达网络分析

目的通过加权基因共表达网络分析(weighted gene co-expression network analysis,WGCNA)分析久坐影响成年人胰岛素敏感性的关键基因及其相关功能和机制。方法针对下载于NCBI的人类基因表达汇编(gene expression omnibus,GEO)数据库的数据集GSE9103,利用WGCNA进行模块化分析,并通过相关性分析和网络可视化确定模块内关键hub基因及其调控网络,再对关键hub基因进行基因本体GO分析和基因集富集GSEA分析。结果通过构建共表达网络,最终得到5个基因模块,其中yellow模块与久坐呈正相关,而NDUFA9为其枢纽基因;tan模块与久坐呈负相关,而MYBPC2为其枢纽基因。进一步对GO和GSEA分析发现,yellow模块基因包括NDUFA9主要发挥线粒体形成、有氧呼吸、脂肪酸降解等生物学功能并参与物质代谢、三羧酸循环、阿尔茨海默氏症、帕金森病相关通路;而tan模块基因包括MYBPC2则与骨骼肌收缩、糖原代谢、心脏发育相关。结论本研究通过构建基因共表达网络筛选出与久坐生活方式显著相关的2个关键基因模块和2个枢纽基因;这2个模块基因分别通过线粒体功能、物质代谢和骨骼肌收缩、糖原代谢而影响成年人对胰岛素的敏感性。

基于加权基因共表达网络分析的口腔黏膜下纤维化与其他器官纤维化的比较分子分析

目的:目前尚无成熟且经济的动物模型模拟口腔黏膜下纤维化(oral submucous fibrosis,OSF),这是制约OSF的机制和药物研究的主要障碍之一。本研究参考泛癌分析的底层逻辑,从分子、信号通路、生物学过程等方面对OSF及其他4种器官纤维化进行系统比较,这有利于学者们发现不同器官纤维化基因组之间的异同,寻找到某些普遍规律,也为OSF的研究提供新思路。方法:从基因表达综合数据库(Gene Expression Omnibus,GEO)网站下载GSE64216、GSE76882、GSE171294、GSE92592和GSE90051数据集的芯片数据。用“Limma”软件包检测各类型纤维化的差异表达mRNAs(differentially expressed mRNAs,DEmRNAs)。采用加权基因共表达网络分析(weighted gene co-expression network analysis,WGCNA)鉴定各类型纤维化相关模块。通过功能富集和通路富集分析各纤维化相关模块基因的异同。结果:OSF、肾间质纤维化(kidney intestinal fibrosis,KIF)、肝纤维化(liver fibrosis,LF)、特发性肺纤维化(idiopathic pulmonary fibrosis,IPF)、皮肤纤维化(skin fibrosis,SF)中分别检测到6057、10910、27990、10480和4801个DEmRNAs。通过WGCNA对各器官纤维化相关模块进行识别,分别构建各类型纤维化的共表达网络。除了KIF和LF有5个共同的hub基因外,其他纤维化疾病之间没有共同的hub基因。OSF、KIF、LF、IPF和SF的共同通路主要集中于免疫相关通路。结论:OSF和其他4种器官纤维化在分子水平上具有组织和器官特异性,但它们有许多共同的信号通路和生物学过程,主要是在炎症和免疫方面。

应用基因集富集分析和加权基因共表达网络分析方法分析类风湿关节炎的枢纽基因及意义

背景:用目前的免疫方法检测类风湿关节炎患者的传统血清标志物,仍有30%为阴性;近年发现的血清生物标志物对类风湿关节炎诊断的敏感性和特异性较低,不适于临床推广。目的:应用基因集富集分析和加权基因共表达网络分析方法分析类风湿关节炎枢纽基因及其意义。方法:从GEO数据库下载GSE1919、GSE55235数据集,使用R软件“limma”包初筛差异表达基因、“clusterProfiler”包进行加权基因共表达网络分析,取交集基因作为候选枢纽基因。将结果导入String数据库构建蛋白质相互作用网络,Cytoscape软件cytoHubba插件筛选枢纽基因,并进行基因本体分析、京都基因与基因组百科全书通路富集分析。同时对2个数据集用基因集富集分析方法进行京都基因与基因组百科全书通路分析。结果与结论:①共筛选出10个枢纽基因,其基因本体分析集中在淋巴细胞分化、T细胞分化、质膜信号受体复合物、主要组织相容性复合体蛋白结合等生物学功能;②京都基因与基因组百科全书通路富集于T细胞受体信号通路、Th1辅助细胞和Th2辅助细胞分化、Th17辅助细胞分化、肿瘤细胞程序性死亡-配体1表达和程序性死亡受体1检查点通路及原发性免疫缺陷5条通路;③基因集富集分析2个数据集类风湿关节炎样本发现3条共同通路:哮喘、自身免疫性甲状腺疾病和细胞黏附分子通路显著上调;④受试者工作特征曲线结果显示,5个基因(CD27、IL2RG、CD8A、LCK、NKG7)对诊断类风湿关节炎具有良好的特异性和敏感性(曲线下面积>0.85);⑤CD4可能是多能的基因网络协调员,在免疫反应中发挥重要功能。

利用加权基因共表达网络分析构建食管腺癌预后枢纽基因网络

目的 利用加权基因共表达网络分析(WGCNA)构建食管腺癌预后枢纽基因网络,筛选与食管腺癌预后相关的枢纽基因.方法 提取癌症和肿瘤基因图谱计划(TCGA)数据库中食管腺癌的标本数据,利用WGCNA匹配基因表达数据与临床预后数据,构建食管腺癌的基因共表达网络模块,并筛选出与预后最相关的模块和枢纽基因,再根据共表达权重系数大小在Cytoscape软件中进行可视化.结果 共获得9例正常食管组织和78例食管腺癌的标本,其中男67例,女11例;年龄28.0~86.6[68.4(58.0,77.1)]岁;TNM分期:Ⅰ期10例,ⅡA期9例,ⅡB期16例,Ⅲ期33例,Ⅳ期10例.共筛选得到的多个模块中深蓝色模块与预后最相关,在深蓝色模块中识别出19个枢纽基因(PAQR8、MAMDC4、CEP44、SCRG1、UGT2B15、NUDT9、FOLH1、SFTPB、FBXO8、TENM1、CASR、NEB、SMIM19、SLC20A2、RNF170、SCN2A、GOLIM4、ICK、DNAH2),并构建了枢纽基因的共表达网络.其中基因间共表达权重系数最大的3对基因分别是FOLH1和SCRG1、FOLH1和UGT2B15、FOLH1和SFTB.结论 利用WGCNA可识别与食管腺癌预后相关的枢纽基因,为食管腺癌的治疗提供新靶点.

基于Copula函数的权重基因共表达网络分析法

利用Frank-Copula函数构建相关系数,并应用于权重基因共表达网络分析模型的改进。同时,利用该改进模型分析小鼠基因数据,发现在相关度最大的模块中有67个基因与小鼠肥胖导致的体质量问题相关。其中,模型的筛选结果有效率为62.6%,说明其在功能基因筛选应用前景中的科学性、有效性和可行性。

加权基因共表达网络分析识别强迫症中的关键基因

目的:利用加权基因共表达网络分析(WGCNA)及套索算法(LASSO)筛选可能影响强迫症发生发展的生物学通路及关键基因,为进一步发现强迫症的分子标志物,确立新的诊断及治疗提供理论依据。方法:从GEO数据库中下载强迫症基因芯片数据GSE78104,包括30例强迫症患者的外周血样本和30例正常人的外周血样本。用WGCNA构建一个基因共表达网络,选择出与强迫症发生高度相关的枢纽模块,分析枢纽模块的生物学功能。通过构建LASSO回归模型筛选与强迫症发生相关的关键基因并检验模型的准确性,最后对关键基因进行差异表达分析。结果:分析发现GSE78104的13模块中,magenta模块与强迫症最为相关(r=0.32,P=0.01)。将magenta作为枢纽模块进一步进行GO和KEGG分析,表明magenta模块主要涉及中性粒细胞活化、调控白细胞凋亡、白细胞介素-1介导的信号通路和TNF信号通路等功能和信号通路(P<0.05)。采用LASSO回归模型筛选得出WDFY3、ZCCHC6、HAL3个关键基因并发现模型准确性较高(AUC=0.78),最后发现3个关键基因在强迫症中均为高表达。结论:WDFY3、ZCCHC6、HAL 3个关键基因可能参与强迫症的发生发展过程。

基于加权基因共表达网络分析及生物信息学分析鉴定肝细胞癌临床特征关键基因

目的:通过加权基因共表达网络分析鉴定肝细胞癌(hepatocellular carcinoma,HCC)临床特征密切相关模块基因,为临床早期诊断及治疗提供参考。方法:从GEO数据库下载GSE84598芯片数据,综合加权基因共表达网络分析提取与HCC临床特征密切相关的模块基因。通过最大集团中心性(maximal clique centrality,MCC)算法对蛋白相互作用网络分析,鉴定出枢纽基因;最后通过TCGA数据库验证枢纽基因的表达情况、Kaplan-Meier Plotter在线数据库评估枢纽基因与HCC患者的预后关系。结果:通过对比HCC组织样本与正常肝组织样本的基因表达数据,共获得6262个差异表达基因,其中上调2207个,下调4055个。运用加权基因共表达网络分析鉴定出关键模块的120个基因;通过与差异表达基因取交集,得到候选枢纽基因115个。富集分析结果显示,候选枢纽基因与细胞有丝分裂、p53信号通路等密切相关。进一步运用MCC算法对115个候选枢纽基因的蛋白相互作用网络进行分析,鉴定出5个枢纽基因,即NUF2、RRM2、UBE2C、CDC20和MAD2L1。通过TCGA数据库对枢纽基因的验证,结果发现,与正常肝组织相比,在HCC组织中5个枢纽基因均明显上调;而且生存分析显示枢纽基因的高表达与HCC患者的不良预后密切相关。结论:本研究通过结合多个数据库鉴定出了5个枢纽基因,为HCC的临床诊断及治疗提供方向。