降糖甲治疗2型糖尿病的机制研究

基于生物信息学、网络药理学和分子对接技术探讨降糖甲颗粒(JTJ)治疗2型糖尿病(T2DM)的作用机制。通过TCMSP、Uniprot数据库、swisstarget和SEA数据库筛选降糖甲颗粒的活性化合物和作用靶点,从GEO数据库获取T2DM患者组胰岛样本和非T2DM患者组胰岛样本数据,利用R语言limma包筛选差异表达基因,二者取交集得到JTJ-T2DM共有靶点。通过Cytoscape3.9.1绘制JTJ-T2DM-靶点网络图,利用STRING数据库、Cytoscape3.9.1的cytoHubb插件制作JTJ-T2DM共有靶点的PPI网络图并获取关键靶点(Hubs)。利用DAVID平台对JTJ-T2DM共有靶点进行GO功能富集分析及KEGG通路富集分析,对共有靶点有差异表达并进行ROC分析,最后通过DS BIOVIA软件进行分子对接验证。结果显示,最终获得41个JTJ活性化合物,909个药物作用靶点,2952个T2DM靶点,JTJ-T2DM交集靶点160个,Hubs 12个;GO及KEGG富集分析结果显示,JTJ主要参与细胞凋亡等过程,调控MAPK及AGE-RAGE等信号通路而治疗T...     

不同浓度IAA影响平菇菌丝生长的调控机制研究

采用添加不同浓度(10~(-3)、0和10~(-8 )mol/L)IAA的PDA培养基分别培养P99平菇(Pleurotus ostreatus)菌丝7 d和15 d.收集菌丝体,提取总RNA并进行转录组测序.对测序信息进行生物信息学分析,包括差异表达基因(differentially expressed genes, DEGs)的筛选、聚类分析、gene onotology(GO)和kyoto encyclopedia of genes and genomes(KEGG)富集分析.最后采用实时荧光定量PCR验证色氨酸代谢通路中关键基因的表达.结果显示:与对照相比,高浓度(10~(-3 )mol/L)和低浓度(10~(-8)mol/L)IAA处理后,在菌丝生长的第7 d,分别获得806个和412个DEGs;在菌丝生长后期,即第15 d,分别获得145个和26个DEGs.GO富集分析表明,两种浓度IAA处理条件下,DEGs主要富集在次生代谢、木质素代谢、几丁质代谢和氧化还原酶代谢等生物学过程中.KEGG富集分析结果表明,两种浓度IAA处理后,生长前期和后期的DEGs富集通路均包括了色氨酸代谢通路和MAPK代谢通路.实时荧光定量PCR结果显示,色氨酸代谢通路中关键基因DN3209和DN3495的表达水平和测序数据一致.以上结果表明,几丁质代谢、色氨酸代谢和MAPK信号通路可能参与了不同浓度的IAA对平菇菌丝生长的促进和抑制作用.     

大肠杆菌在葡萄糖和甘油碳源培养条件下的基因表达谱分析

【目的】探讨葡萄糖和甘油碳源对大肠杆菌(Escherichia coli)基因表达谱的影响。【方法】通过基因表达数据库GEO下载基因芯片数据集GSE2037和高通量测序数据集GSE156143,分别对两个数据集中的葡萄糖和甘油碳源样本进行差异表达基因筛选,对筛选到的差异表达基因进行并集后,进行GO功能富集分析、KEGG通路富集分析以及蛋白质相互作用网络构建。【结果】筛选出701个差异表达基因,差异表达基因主要富集于小分子分解代谢过程、碳水化合物运输、趋化性等生物学过程,细胞器、细菌型鞭毛、甲基接受趋化蛋白复合物等细胞组分,离子跨膜转运蛋白活性、碳水化合物结合、肌动活动等分子功能以及不同环境中的微生物代谢和ABC转运体相关通路。此外,筛选出的14个核心基因在鞭毛合成和趋化中发挥作用。【结论】大肠杆菌在葡萄糖和甘油碳源培养条件下具有不同的基因表达模式,所确定的14种与趋化和鞭毛合成有关的基因有助于进一步揭示大肠杆菌在应对不同营养环境变化时采取的分子调控机制。     

基于生物信息学方法分析GABRD基因在结肠癌中的表达及预后情况

采用生物信息学方法探讨GABRD基因在结肠癌样本中的表达及预后情况。通过UCSC XENA下载33种肿瘤类型和正常组织的RNA序列数据和相关临床数据,使用R软件分析GABRD基因在结肠癌样本中的表达,并筛选共表达基因,对其进行富集分析;分析GABRD基因对结肠癌患者生存及预后的影响,并建立预后列线图;构建GABRD基因的蛋白质-蛋白质相互作用(protein-proteininteraction, PPI)网络并筛选关键模块及枢纽基因,验证枢纽基因的生存及临床诊断价值。结果表明：GABRD基因在结肠癌样本中高表达并影响患者生存,筛选得到369个共表达基因,基因本体论(gene ontology, GO)功能富集发现其主要参与G蛋白偶联等生物学过程,京都基因与基因组百科全书(Kyoto encyclopedia of genes and genomes, KEGG)通路富集显示其主要参与AMPK等信号通路;构建出由51个节点和523个连接组成的PPI网络,筛选枢纽基因5个,其中2个显著影响生存,5个具有临床诊断价值。综上,GABRD基因在结肠癌样本中高表达,影响结肠癌患者生存及预后,可能...     

具有预后价值的乳腺癌发病关键基因鉴别研究

目的:筛选与乳腺癌发病相关的关键基因,为研究乳腺癌的诊疗提供新的潜在分子靶标.方法:使用GEO2R在线工具比较乳腺癌与正常乳腺组织基因表达情况,进行差异显著性分析,利用基因功能注释工具DAVID对差异基因进行GO和KEGG富集分析.通过STRING数据库构建差异基因的蛋白互作网络,基于最大团中心性(maximal clique centrality,MCC)算法鉴别关键基因,利用Kaplan Meier生存分析验证关键基因对患者生存时间的影响.结果:乳腺癌基因表达差异分析共产生491个差异基因,其中有254个上调,237个下调.GO富集分析表明差异基因显著富集于肿瘤相关的生物学过程、细胞组分和分子功能,KEGG通路富集分析主要集中于PI3K-Akt信号通路、黏附斑和癌症通路.基于MCC算法共选取10个关键基因,其中的4个基因(ISG15,IFIT1,GBP1和IFI27)过表达与患者生存时间显著降低密切相关(P0.05).结论:共鉴别了4个与乳腺癌发病密切相关的关键基因,相关研究结果可为研究乳腺癌的诊疗提供新的潜在分子靶标.     

转录组分析sgf73基因缺失对粟酒裂殖酵母生长代谢影响的分子机制

为研究粟酒裂殖酵母sgf73基因缺失后的关键基因和关键代谢通路,采用RNA-Seq测序技术对野生型酵母菌株及sgf73基因缺陷型菌株进行测序及生物信息学分析,并进行GO和KEGG功能富集分析.结果表明：sgf73基因缺陷型菌株中高表达基因有1 834个,极高表达基因有6个,且有1 714个差异表达基因,包括934个上调基因,780个下调基因;sgf73基因敲除导致细胞代谢和跨膜转运过程异常变化;MAPK信号通路中参与周期调控cki1,cki2和cdc25基因的下调导致sgf73Δ菌株有丝分裂时间延长;调控细胞骨架rgf2,rho1和stt4基因的下调导致sgf73Δ菌株肌动蛋白环收缩异常.     

肝细胞癌关键基因的生物信息学分析

肝细胞癌是一种高死亡率的原发性肝癌,其有限治疗和较低的化疗敏感性,使得迫切需要寻找潜在的临床治疗靶点和预后判断的生物标志物.为此,采用生物信息学方法对肝细胞癌发生发展的关键基因进行挖掘.从基因表达数据库(GEO)中下载数据集,并筛选差异表达基因,使用在线数据库DAVID 6.8对筛选到的DEGs进行基因本体(GO)富集...     

大气压低温等离子体抑制甲状腺未分化癌CAL-62细胞生长的机制研究

为了探究大气压低温等离子体(CAP)抑制间变性甲状腺癌的作用机制，应用转录组测序技术，借助生物信息学分析手段研究CAP对CAL-62的抑制作用，筛选出差异基因进行KEGG富集分析和GO分类，最后选取部分基因进行Western blot (WB)验证.结果表明，与对照组相比CAP处理组共有674个差异表达基因，以BIK和CDKN1C等为代表的287个基因上调表达，以TANC2、ESM1和CBL等为代表的387个基因下调表达.KEGG富集分析表明差异基因主要富集在MAPK、Rap1等信号通路及癌症转录失调等.GO分类表明差异基因主要富集在免疫反应、细胞凋亡过程等方面.WB检测部分差异基因的变化趋势，得到与转录组相一致的结果.该研究从转录组学角度去解释CAP抑制CAL-62细胞的机制，以期为CAP应用于甲状腺未分化癌治疗提供一定的理论基础.     

基于网络药理学和分子对接技术探讨五味子治疗慢性咳嗽的机制研究

目的 基于网络药理学和分子对接技术探讨五味子治疗慢性咳嗽的关键成分和潜在作用机制.方法 通过TCMSP数据库筛选五味子的活性成分及相应靶点,运用TTD、DisGeNET和DrugBank数据库收集咳嗽靶点集,再通过Venny平台筛选五味子活性成分和咳嗽的相关交集靶点.利用String网络平台和Cytoscape 3.8.2软件构建相关网络图及直方图,并得到关键靶点,通过Webgestalt网站对关键靶点进行KEGG和GO的富集分析.利用AutoDock Tools和Autodock Vina软件将关键靶点和核心成分进行分子对接,再通过Pymol软件进行可视化分析.结果 共筛选出8种有效活性成分和877个与活性成分相关的作用靶点;与咳嗽交集的相关靶点有88个,其中Akt1和Tnf为关键靶点.KEGG通路富集分析筛选得到27条通路,包括FcεRI信号通路、脂肪细胞因子信号通路、Toll样受体信号通路等;GO功能富集分析得到205个GO条目,主要涉及阴离子跨膜转运的负调控、细胞表面蛋白质定位的正向调节等生物过程,还涉及吞噬杯、囊泡等细胞组成,一氧化氮合酶调节活性、GTP酶激活蛋白结合等分子功...     

基于网络药理学探究百合抗抑郁的作用机制

通过网络药理学的方法探究百合抗抑郁的作用机制.利用中药系统药理学技术平台(TCMSP)等数据库,以生物利用度(OB)≥30%、类药性(DL)≥0.18为条件筛选百合活性成分及抗抑郁的靶点,构建成分-靶点、关键靶点蛋白网互作图、靶点基因GO富集及KEGG通路分析.结果表明,通过筛选获得7个符合标准的活性成分、85个与抑郁相关的靶点,GO富集分析得到17个生物学功能,KEGG通路分析得到11条显著富集的通路.百合通过多靶点、多通路、多成分发挥抗抑郁的作用,可能是通过胆碱能突触、神经活性配体-受体、p53信号通路、钙信号通路等信号通路发挥复杂的调节作用.     

有氧运动对伴随衰老代谢综合征基因表达的调控

目的：探索有氧运动对代谢综合征患者血液基因表达的影响，分析运动前后血液转录表达谱的差异，研究差异表达基因涉及的功能与内在联系.方法：从美国国立生物技术信息中心(NCBI)公共基因芯片数据平台(GEO)下载GSE10540数据集mRNA芯片数据，对5名代谢综合征患者运动前后的mRNA表达谱数据进行基因集富集分析(GSEA),应用metascape软件进行功能富集(GO)和信号通路(KEGG)富集分析，用STRING数据库构建蛋白质相互作用网络.结果：运动后血样比较运动前共筛选出152个差异表达基因(differentially expressed genes, DEGs),其中有54个基因上调，98个显著下调.GESA分析发现12个基因或通路在运动后表达谱中显著富集负相关，蛋白质相互作用网络筛选出9个关键蛋白.结论：代谢综合征患者运动前后血样相关DEGs揭示了代谢综合征的分子特征，这些基因可能与癌症发生之间呈负相关，这为进一步探讨有氧运动对代谢综合征患者健康改善的作用机制及关键调控分子提供了研究思路.     

基于BATMAN-TCM在线分析平台的浙贝母-瓜蒌配伍治疗慢性阻塞性肺病的网络药理学研究

基于中药分子机理的生物信息学分析工具（bioinformatics analysis tool for molecular mechanism of traditional Chinese medicine,BATMAN-TCM）挖掘浙贝母-瓜蒌治疗慢性阻塞性肺病的分子作用机制。应用BATMAN-TCM预测浙贝母-瓜蒌的潜在靶点,GeneCards数据库检索慢性阻塞性肺病的靶点;应用Cytoscape 3.7.1软件构建化合物-疾病-靶点相互作用网络、蛋白质-蛋白质相互作用( protein-protein interaction,PPI) 网络,进行基因本体( gene ontology,GO) 功能富集分析和京都基因与基因组百科全书( Kyoto encyclopedia of genes and genomes,KEGG) 的生物通路富集分析。结果筛选得到浙贝母-瓜蒌中21种活性成分,基于STRING构建的浙贝母-瓜蒌靶点PPI网络包含131个节点和236条相互作用关系,关键节点包括ANXA1、EDN1、HDAC1、RXRA、ADRA2A、OXT、RXRB等。GO富集分析途径主要涉及离子门控通道活性、底物特异性通道活性、被动跨膜转运活性、细胞外配体门控离子通道活性、神经递质受体活性、配体门控离子通道活性等。KEGG通道分析主要涉及神经活动配体-受体相互作用、cAMP信号通路、钙离子信使通路、类胆碱能突触、逆行内源性大麻素信号传导、5-羟色胺能突触、多巴胺能突触等途径。研究结果提示浙贝母-瓜蒌关键的作用靶点包括ANXA1、HDAC1等,信号通路方面涉及类胆碱能、5-羟色胺、多巴胺等神经受体相关通路,体现了中药复方的整体性特点,为进一步研究浙贝母-瓜蒌的分子作用机制提供了新思路。     

乳腺癌-冠心病共享生物标志物筛选

运用生物信息学方法筛选乳腺癌-冠心病标志物,为乳腺癌诱发的冠心病治疗提供潜在的作用靶点.从基因表达数据库(gene expression omnibus,GEO)中下载乳腺癌和冠心病相关表达谱芯片数据,使用GEO2R筛选差异表达基因,依据Venn图交集获取差异共表达基因,通过DAVID网站进行基因功能注释(gene ontology,GO)及京都基因与基因组百科全书(Kyoto Encyclopedia of Genes and Genomes,KEGG)生物功能富集分析,STRING网站和Cytoscape 3.7.2软件进行蛋白互作分析.GE-PIA和Kaplan-Meier Plotter进行对乳腺癌患者hub基因mRNA表达水平和预后分析.结果表明:2个数据集筛选得到差异表达基因286个,基于在乳腺癌中mRNA显著性表达水平筛选出45个基因.GO功能富集分析发现差异表达基因主要在泛素蛋白转移酶活性、糖蛋白结合、泛素蛋白连接酶结合等生物学过程发挥作用.KEGG分析显示差异基因主要参与缝隙连接、肾素分泌、5-羟色胺能突触、谷氨酸能突触、血管平滑肌收缩、血小板活化、癌细胞蛋白多糖等多条信号通路.基因mRNA表达水平和预后分析显示NLN、POSTN、MAPT、MYO6、MAP1B、FBXO31、KIT、PIK3R1等8个与冠心病相关的hub基因参与乳腺癌的发生、发展过程.NLN、POSTN、MAPT、MYO6、MAP1B、FBXO31、KIT、PIK3R1可作为检测乳腺癌诱导冠心病的潜在标志物.     

膜性肾病差异性表达microRNA的靶基因功能分析

目的:探讨膜性肾病(MN)差异性表达microRNA的靶基因主要参与的Gene Ontology(GO)富集与Kyoto Encyclopedia Genes And Genomes(KEGG)通路.方法:100例MN患者与100例健康对照者(NC)作为实验对象.100例NC随机分成10组(NC1~10),每组10人.同样100例MN随机分成10组,每组10人(MN1~10).分别从实验参与者静脉采取全血10 mL,来源于静脉血的单个核细胞用于提取RNA.总RNA以组为单位等量混合用于高通量技术测序(high-throughput sequencing),对RNA测序结果进行长度分布,基因比对,RNA分类注释,RNA特有序列及其公共序列分析后,获得小核糖核酸(microRNA)进行MN与NC差异性表达分析来寻找显著表达microRNA.应用相关的软件对差异性表达的microRNA进行靶基因预测,靶基因借助功能注释软件对其参与的生物过程进行GO富集及KEGG通路分析.结果:靶基因在GO富集分析中主要集中在细胞器官组成,细胞膜组成,离子结合,生物代谢过程,生物调节过程等.靶基因在KEGG通路分析中主要参与嘌呤代谢通路,代谢产物生物合成,癌症通路,蛋白激酶信号通路等.结论:MN与NC之间存在着显著差异性表达microRNA.差异性表达的microRNA的靶基因参与的GO富集与KEGG通路可能与MN的发病机制与临床症状有着密切的联系.     

亚麻木酚素干预下ApoE基因敲除小鼠肝脏转录组分析

利用转录组测序技术研究了亚麻木酚素干预下ApoE-/-(Apolipoprotein E,ApoE)小鼠肝脏组织中差异表达基因及相关信号通路，为明确亚麻木酚素缓解高脂血症的分子机制奠定基础。结果表明，亚麻木酚素组和模型组之间，共有372个差异表达基因，其中包括234个显著上调基因，138个显著下调基因。这些关键的差异表达基因包括Hsd3b5、Acot1、Elovl3等与脂质代谢相关的基因。GO(Gene ontology, GO)富集分析发现差异表达基因主要富集在脂肪酸代谢等生物过程、细胞外间隙等细胞组分以及单加氧酶活性等分子功能。KEGG(Kyoto encyclopedia of genes and genomes, KEGG)信号通路富集分析表明差异表达基因在PPAR(Peroxisome proliferation activated receptors, PPAR)信号通路被富集。PPAR信号通路的多个基因的表达被激活，如Fabp、Cyp、Acsl、Ehhadh、Fads、Acaa、Scp、Pltp等，推测亚麻木酚素是通过PPAR信号通路调节ApoE-/-小鼠血脂水平的。     

水分胁迫下马铃薯萌芽出苗期 根系转录组差异表达分析

以马铃薯‘克新1号’为材料,通过Illumina HiSeqTM 4000高通量测序技术,对不同水分胁迫后的马铃薯萌芽出苗期根系cDNA文库进行了RNA-Seq分析,拟探明马铃薯感应水分胁迫的分子机制,挖掘出与抗旱密切相关的功能基因.将比对到基因组上的基因利用RPKM衡量各样本间的基因表达丰度,以P0.05筛选出差异表达基因,通过GO (基因本体,gene ontology)数据库、KEGG代谢途径数据库,对不同水分胁迫样本差异表达基因的功能和参与的调控路径进行分析;选择6个差异表达基因,利用实时荧光定量PCR验证RNASeq结果的可靠性.结果表明:(1)DLWR(重度水分胁迫)样本中,比对到基因组的表达基因30 114 133个,差异表达基因5 241个;LWR(中度水分胁迫)样本中,比对到基因组的表达基因32 858 890个,差异表达基因1 488个.2个样本中同时存在的差异基因369个,其中上调表达基因78个,下调表达基因291个;显著差异表达基因主要与蛋白激酶、水解酶、氧化还原酶、水通道蛋白、转运蛋白、转录因子等相关.(2)在GO富集分析中,2个水分胁迫处理样本在生物学功能中富集的类别并不完全相同,但显著富集的主要是与氧化还原、转运、膜结合、转录调节子等相关的基因.(3)KEGG代谢分析显示,差异表达基因显著富集在谷胱甘肽代谢、植物激素信号转导等途径,与植物抗逆能力密切相关.(4)利用qRT-PCR验证的6个基因表达模式与RNA-Seq分析结果一致,证实了RNA-Seq结果的可靠性.通过研究得到了不同水分胁迫下马铃薯根系转录组信息,获得了差异基因及其功能信息,阐明了差异基因所参与的代谢通路及与植物水分胁迫的关系.     

基于生物信息学对肝细胞癌潜在枢纽基因的筛选及验证

目的:筛选网络数据库中肝细胞癌组织与非癌性肝组织中差异表达的基因,并探讨其中适合作为早期筛查肝癌及治疗晚期肝细胞癌(HCC)的分子标记.方法:从GEO数据库中下载3个数据集GSE121248、GSE101685和GSE14520,筛选差异表达显著的基因(DEGs),并对其进行基因功能富集和信号通路富集分析,利用STRING网站构建对应蛋白的相互作用网络,导入Cytoscape软件后用CytoHubb插件筛选处于关键位置前10位的枢纽基因.利用Kaplan-Meier生存分析比较枢纽基因表达差异对肝癌患者平均生存率的影响,并通过GEPIA验证枢纽基因在肝癌组织中的表达变化情况,进一步分析各枢纽基因所涉及的信号通路.结果:数据集GSE121248、GSE101685和GSE14520分别筛选出为187、536和253个DEGs(|log_2FC|2, adj.P0.01).Venn图显示3个数据集中共同包含的DEGs有87个,其所涉及的细胞功能为类固醇代谢和氧化还原酶活性,所涉及的通路为视黄醇代谢及p53信号通路.87个DEGs对应的69个蛋白的相互作用网络中,筛选最大团中心性(maximal clique centrality, MCC)前10位的枢纽基因为CDK1、CCNB1、TOP2A、PRC1、ECT2、RACGAP1、DTL、RRM2、ASPM和ESR1.Kaplan-Meier生存分析提示9个基因(CDK1、CCNB1、TOP2A、PRC1、ECT2、RACGAP1、DTL、RRM2和ASPM)高表达可能导致HCC患者的总生存期(OS)较差,而ESR1基因高表达时患者OS延长(P0.05).GEPIA分析369例HCC组织和160例正常肝组织的结果显示,CDK1、CCNB1、TOP2A、PRC1、ECT2、RACGAP1、DTL、RRM2和ASPM表达升高,而ESR1表达下调(|log_2FC|1,P0.05).这10个枢纽基因所涉及的信号通路主要为p53信号通路、细胞周期和孕酮介导的卵母细胞成熟.结论:9个枢纽基因CDK1、CCNB1、TOP2A、PRC1、ECT2、RACGAP1、DTL、RRM2和ASPM的高表达以及ESR1的低表达与肝癌患者的不良生存率密切相关,提示它们作为检测对象预测肝癌预后及作为干预靶点改善肝癌预后的可能.     

中华按蚊雌雄蚊转录组比较分析

【目的】对中华按蚊(Anopheles sinensis)雌、雄蚊的转录组测序数据进行差异表达分析,获得两者间差异表达基因,对它们进行功能富集分析,挖据这些基因的功能表现和参与的通路。【方法】基于已知中华按蚊雌、雄蚊转录组数据进行了RNA-Seq数据相关性检查,利用R语言DESeq包比较雌蚊和雄蚊中基因的差异表达情况,并用GOseq和KOBAS对筛选出来的差异表达基因进行GO和KEGG富集分析。【结果】Pearson相关系数的平方(R2)大于0.92,说明样本在3个生物学重复间具有高度相关性,可用于进一步分析;中华按蚊雌、雄蚊差异表达基因共3 298个,其中1 589个基因在雌蚊中表达上调,1 709个基因在雌蚊中表达下调;差异表达基因明显富集在酰胺生物合成、有机酸代谢、肽代谢过程等功能类别上,参与核糖体在真核生物中的生物发生、RNA转运等KEGG通路。【结论】上述结果对深入研究不同性别的中华按蚊基因表达及差异具有潜在价值。     

三七-白及药对作用机制的网络药理学分析

通过网络药理学的方法,研究三七-白及药对在疾病治疗过程中潜在的药理作用机制.采用TCMSP数据库筛选三七-白及药对的活性成分和作用靶点,并用Uniprot数据库校正靶点信息得到靶点基因,利用CTD数据库获得靶点基因相关的疾病类型,通过STRING数据库构建靶点蛋白相互作用网络,分析得到核心蛋白,运用DAVID数据库富集分析靶点基因参与的基因本体论(GO)生物学过程及京都基因与基因百科全书(KEGG)通路.共筛选得到17个活性成分,涉及193个作用靶点.结果表明:槲皮素、β-谷甾醇和豆甾醇的作用靶点数目最多;靶点基因与35类疾病相关,主要包括癌症、神经系统疾病、心血管疾病等;蛋白相互作用网络分析得到核心蛋白AKT1,MAPK1,c-Jun,P53,TNF等;靶点基因主要涉及药物反应,RNA聚合酶Ⅱ启动子转录的正调控等91条GO生物学过程,KEGG通路显著富集到癌症通路、乙型肝炎等66条通路,与前述疾病分析结果相符.     

乳腺癌脑转移差异表达基因的生物学功能分析

乳腺癌脑转移(breast cancer brain metastasis,BCBM)的发病机制尚未明确。为了探究BCBM的发病机制,对BCBM差异表达基因的生物学功能进行研究并筛选关键调控基因。从基因表达综合数据库(gene expression omnibus,GEO)下载4个BCBM基因表达谱数据(GSE12237、GSE100534、GSE125989以及GSE43837),采用R语言筛选差异表达基因,采用富集分析包括基因本体分析(gene ontology,GO)和京都基因与基因组百科全书分析(Kyoto encyclopedia of genes and genomes,KEGG)进行生物学功能分析,采用STRING和Cytoscape分析蛋白质相互作用网络,采用Kaplan-Meier进行生存分析。结果表明,同时存在于2个及以上基因表达谱数据中的差异表达基因261个,GO分析主要涉及细胞外基质组织、细胞外结构组织等生物过程,细胞外基质结构组成、胶原结合等分子功能,含有胶原的细胞外基质、胶原蛋白三聚物等细胞组分;KEGG分析主要涉及蛋白质消化和吸收、局部黏附等通路。蛋白质相互作用网络分析得到9个关键调控基因,其中,DCN、COL6A1与BCBM的生存率显著相关,可作为潜在的BCBM关键调控基因,并为BCBM分子机制的研究提供思路。