侧孢短芽孢杆菌S2-31拮抗下细辛叶枯病菌转录组差异表达分析

通过分析侧孢短芽孢杆菌拮抗作用下叶枯病菌的转录组学特征,研究差异表达基因(DEGs)和代谢通路的富集情况,初步探索侧孢短芽孢杆菌拮抗叶枯病菌的分子机制。首先利用S2-31与叶枯病菌的对峙培养观察其拮抗作用,然后利用转录组测序探究侧孢短芽孢杆菌拮抗下和正常生长下的叶枯病菌的基因表达水平差异,并进行RT-qPCR验证,最后,利用相关数据库对DEGs进行注释和富集分析。对峙培养后,叶枯病菌菌丝出现皱缩、扭曲等现象。测序后共得到3681个DEGs,其中2224个相对上调,1457个相对下调。GO功能注释分析表明,上调和下调表达的DEGs均注释到生物过程8个亚群、细胞组分8个亚群和分子功能4个亚群。KEGG富集分析中,共有732个unigenes定位到115条生物学通路中,其中富集程度相对较高的通路有氨基糖和核苷糖代谢、糖基磷脂酰肌醇(GPI)锚定生物合成和过氧化物酶体;涉及差异基因较多的通路有淀粉和蔗糖代谢、剪接体和胞吞作用;上调表达的DEGs主要富集在代谢途径,下调表达的DEGs主要富集在遗传信息处理和细胞过程途径。从显著富集的通路中筛选出与菌丝生长发育相关的差异表达基因,结果表明侧孢短芽孢杆菌主要抑制病原菌菌丝过氧化物酶体的形成、胞吞过程和GPI锚定的合成等过程。差异表达基因的RT-qPCR验证结果与转录组测序的结果一致。在侧孢短芽孢杆菌S2-31的拮抗下抑制了叶枯病菌的生长,其转录组特征也发生显著变化,差异基因主要涉及代谢、细胞过程以及遗传信息处理等途径的相关通路。     

忍冬属植物中环烯醚萜苷类化合物的研究进展

本文主要综述了忍冬属植物中环烯醚萜类化合物的结构类型、分离方法、测试方法及其药理活性研究;由报道可知,忍冬属植物中分离得到的环烯醚萜类化合物有70多个,并且分属于三种不同的结构类型;无论属于哪种结构类型,环烯醚萜的分离方法均可概括为溶剂浸泡、萃取、柱层析细分,柱层析细分主要有反相柱层析、硅胶柱层析、凝胶柱层析。对于环烯醚萜的定性或定量测试,目前采用的主要方法是HPLC;环烯醚萜类化合物的药理活性研究目前主要集中在抗炎、止痛、解热、保肝及抗病毒方面;除此之外,文中对部分环烯醚萜化合物的生物合成途径也进行了探讨研究。     

药用植物裂环烯醚萜苷类化合物生物合成途径关键酶基因研究进展

裂环烯醚萜苷类化合物是植物体内产生的具有保肝、消炎、降血糖血脂等多重生物活性的次级代谢产物,在临床上应用广泛。该文依据近年来国内外有关裂环烯醚萜苷类化合物生物合成途径及关键酶基因的挖掘与调控机理研究进展,主要对药用植物裂环烯醚萜苷类化合物生物合成途径、关键酶(GPPS、GES、G10H、8HGO、IS、IO、7DLGT、DL7H、LAMT、SLS)与编码基因的研究进展进行综述,为进一步阐明其生物合成途径机制与关键酶调控作用、提高有效活性成分积累、减缓药用植物野生资源紧张等问题提供参考。     

滇龙胆与青叶胆环烯醚萜类物质计量特征分析

龙胆科龙胆属滇龙胆(Gentiana rigescens)与獐牙菜属青叶胆(Swertia mileensis)为我国特有种。采用高效液相色谱法测定其中马钱苷酸、獐牙菜苦苷、龙胆苦苷和当药苷含量,结合多变量分析研究环烯醚萜类成分含量及其构成比例在龙胆科植物种内和种间的差异。研究结果显示,4种环烯醚萜类成分的计量特征呈现种间差异,依据环烯醚萜类成分含量及其构成比例,所有滇龙胆样品可分为4类。马钱苷酸含量与纬度呈极显著负相关(R=-0.348,P0.05),与海拔呈极显著正相关(R=0.307,P0.01);獐牙菜苦苷含量与经度呈极显著负相关(R=-0.592,P0.01),高海拔有利于植株中当药苷构成比例的增加(R=0.245,P0.05);地理因子对植株中龙胆苦苷含量变化影响不显著(P0.05)。     

不同鲜食品质橄榄果实转录组测序及酚类代谢途径相关调控基因挖掘

【目的】果实酚类物质类别、含量是与橄榄营养及风味密切相关的重要品质性状,探究橄榄酚类物质生物合成的分子调控机制。【方法】以总酚含量差异显著的橄榄（低酚/高酚）为试材,分别取花后80-160 d的果实进行转录组分析,对莽草酸-水解单宁/苯丙烷-类黄酮生物合成通路差异基因进行表征,并对差异表达基因进行WGCNA分析,从中挖掘与酚类代谢途径相关的转录因子。【结果】转录组测序共获得296 314条Unigene,其中73%的Unigene被注释到数据库;4个品种（系）橄榄成熟果共鉴定到1 628个差异表达基因（DEGs）,KEGG分析显示,DEGs在酚类代谢途径的“类黄酮生物合成”通路上显著富集;进一步对果实成熟过程莽草酸-水解单宁/苯丙烷-类黄酮生物合成途径的DEGs进行表征,结合WGCNA分析中每个模块与性状的R2和P值,筛选到4个关键模块,根据模块内基因的调控关系利用MCC拓扑分析法以度值≥1挖掘模块内关键转录因子,挖掘到137个转录因子Unigene与30个酚类合成结构基因Unigene共表达,转录因子Unigene功能注释来自35个基因家族,最多的为锌指蛋白（C2C2、C3H、C2H...     

滇杨正/倒扦插苗叶片和树皮组织的转录组比较分析

该研究以滇杨正、倒扦插苗为对象,于速生期(8月)收集插穗主枝萌发处下侧树皮组织和主枝的叶片组织进行转录组比较分析。结果显示:(1)正、倒扦插滇杨苗于2种组织中的转录调控变化存在差异;测序得到358 587条Unigenes,并于叶片组织和树皮组织中分别筛选得到6 512和12 020条DEGs。(2)KEGG在叶片组织中成功注释了552个DEGs,被分配至122个通路中;显著富集通路主要为植物病原菌互作、氨基糖和核苷酸糖代谢、戊糖和葡糖醛酸互换、半乳糖代谢、氨基酸生物合成、抗坏血酸和新陈代谢,其中39个DEGs编码了USPase,影响叶片正常的细胞壁代谢。(3)树皮组织中,分布于128个通路中的1 623个DEGs被KEGG成功注释,主要的富集通路为碳代谢、核糖体、氨基酸生物合成、苯丙素生物合成。(4)大量DEGs富集到与ROS代谢相关的作用酶中,包括ROS产生机制相关的卡尔文循环相关酶、ROS清除机制相关的氧化还原酶(POD)和半胱氨酸合成酶(如cysK和MET17)等。研究认为,上述关键酶类的基因变化,表明倒插扰乱了树皮组织的ROS浓度平衡,激发了倒插苗的生理保护机制。     

基于转录组测序的宁夏枸杞不同品种果实活性成分合成差异表达基因分析

为探究不同品种宁夏枸杞果实活性成分生物合成相关基因的表达水平,筛选关键差异表达基因(differentially expressed genes,DEGs),揭示宁夏枸杞品种间活性成分含量差异的分子机制,本研究采用Illumina NovaSeq 6000高通量测序技术,对宁夏枸杞‘宁杞1号’和‘宁杞7号’青果期、转色期及成熟期果实进行转录组测序,比较2个品种果实不同发育期相关基因表达谱的变化。结果显示:转录组测序共获得811818178条clean reads,有121.76 Gb有效数据。‘宁杞1号’和‘宁杞7号’在青果期、转色期和成熟期差异表达基因分别有2827、2552和2311个;分别有2153、2050和1825个差异基因在基因本体论(gene ontology,GO)、京都基因与基因组百科全书(Kyoto encyclopedia of genes and genomes,KEGG)富集分析和同源蛋白簇(clusters of orthologous groups of proteins,KOG)分析等6个数据库中被成功注释。青果期、转色期和成熟期果实的差异表达基因,在GO数据库分别有1307、865和624个被富集到生物学过程、细胞组分及分子功能3个部分中;KEGG通路富集结果均集中在代谢途径、次生代谢物生物合成和植物-病原互作过程;在KOG数据库,3个发育期分别注释了1775、1751和1541个差异表达基因。对注释的基因进行PubMed数据库检索,在青果期、转色期和成熟期分别筛选到与枸杞活性成分合成相关的差异表达基因18、26和24个,这些基因主要参与类胡萝卜素、类黄酮、萜类、生物碱和维生素等代谢途径。选取7个差异表达基因进行RT-qPCR验证,结果与转录组测序数据表达趋势一致。本研究从转录水平为不同品种宁夏枸杞活性成分含量差异提供了初步证据,为进一步挖掘枸杞活性成分生物合成的关键基因及解析其表达调控机制提供了研究基础。     

基于整合方法分析茶树响应病原真菌胁迫的共有模式

为探讨茶树(Camellia sinensis)对病菌胁迫的共有响应模式和抗病机制,运用生物信息学方法对多组RNA-seq数据进行提取、整合及功能富集,结合多种工具和数据库资源对主要调控分子及蛋白互作模块加以分析。结果表明,病原真菌胁迫下,茶树有较多细胞色素P450家族成员表达显著上调;类固醇和激素的代谢过程、苯丙烷合成途径被激活,有丝分裂细胞周期调控、DNA甲基化等生物过程及光合作用途径受到抑制;主要调控分子如转录因子WRKY和NAC、激酶RLK-Pelle和CAMK等以上调为主。差异表达的蛋白互作模块分析表明,有丝分裂周期调控、基于微管运动、淀粉和蔗糖代谢、细胞壁多糖合成、光合作用、类黄酮代谢模块明显下调,木质素合成和萜类生物合成模块上调;且模块之间可能存在互作。病菌胁迫激活的木质素和萜类合成途径的关键基因包括阿魏酸-5-羟基化酶基因F5H、过氧化物酶基因POD和萜类合成酶基因HMGR等。细胞色素P450基因可能在病菌胁迫中起关键作用,增强木质素和萜类物质的合成、削弱光合作用可能是茶树响应真菌胁迫的核心模式。     

滇龙胆与青叶胆环烯醚萜类物质计量特征分析?

龙胆科龙胆属滇龙胆( Gentiana rigescens)与獐牙菜属青叶胆( Swertia mileensis)为我国特有种。采用高效液相色谱法测定其中马钱苷酸、獐牙菜苦苷、龙胆苦苷和当药苷含量,结合多变量分析研究环烯醚萜类成分含量及其构成比例在龙胆科植物种内和种间的差异。研究结果显示,4种环烯醚萜类成分的计量特征呈现种间差异,依据环烯醚萜类成分含量及其构成比例,所有滇龙胆样品可分为4类。马钱苷酸含量与纬度呈极显著负相关( R＝－0?348, P<0?05),与海拔呈极显著正相关( R＝0?307, P<0?01);獐牙菜苦苷含量与经度呈极显著负相关( R＝－0?592, P<0?01),高海拔有利于植株中当药苷构成比例的增加( R＝0?245, P<0?05);地理因子对植株中龙胆苦苷含量变化影响不显著( P>0?05)。     

苦荞转录因子FtMYBF的克隆、亚细胞定位及表达分析

R2R3-MYB转录因子SG7亚家族成员在植物黄酮醇生物合成中具有极其重要的调控作用。探究苦荞中SG7 R2R3-MYB转录因子在黄酮醇生物合成中的功能,为苦荞黄酮醇生物合成的分子调控机制研究奠定基础。采用RT-PCR技术从苦荞中克隆到一个前期经转录组与代谢组联合分析鉴定到的SG7 R2R3-MYB转录因子,命名为FtMYBF。利用生物信息学、亚细胞定位、基因共表达、基因表达量与总黄酮含量相关性对该基因进行了分析。结果表明,FtMYBF全长CDS序列为1119 bp,编码372个氨基酸。蛋白多序列比对和系统进化树分析表明,FtMYBF属于R2R3-MYB SG7亚家族成员。亚细胞定位显示,FtMYBF蛋白定位于细胞核。基因共表达分析表明,FtMYBF与多个苦荞黄酮醇生物合成结构基因共表达。基因表达量与总黄酮含量间相关性分析表明,FtMYBF在不同组织部位表达量与总黄酮含量高度正相关。FtMYBF属于R2R3-MYB SG7亚家族成员,是一个核定位转录因子,其可能通过正调控苦荞黄酮醇生物合成结构基因的表达来正调控苦荞黄酮和黄酮醇的生物合成。     

三叶木通转录组分析及三萜皂苷生物合成途径关键酶基因的发掘

为了解三叶木通(Akebia trifoliata(Thunb.)Koidz.)的三萜皂苷合成途径及其关键酶,本研究对其花、叶、根、茎进行转录组测序,组装获得了57.25 Gb数据,含140 859个unigenes,序列平均长度为1350 bp。KEGG代谢通路富集结果显示,517个unigenes参与三萜皂苷合成相关的3条代谢途径,其中415个unigenes编码三萜皂苷生物合成途径的19个关键酶。对三萜皂苷生物合成过程中的关键酶角鲨烯环氧酶(SE)进行序列分析和同源建模,发现其具有保守的底物结合结构域。将三叶木通茎与花、叶、根的基因表达水平进行比较,发现茎与根相比较其上调基因数目最多,其中295个差异表达基因(DEGs)与三萜皂苷生物合成途径相关。     

基于转录组的瑶药半枫荷根和叶的差异分析

为了探究半枫荷(Semiliquidambar cathayensis)根和叶的基因表达差异和关键活性成分合成通路中关键基因的表达规律,本研究对半枫荷的根和叶进行转录组测序和生物信息分析。59378个差异表达基因归类到在GO分类的3个大类中,主要与生物学过程有关(50.59%)。626个差异表达基因注释KOG数据库的24个分类中。81个差异表达基因参与苯丙烷类化合物的生物合成,110个差异表达基因参与黄酮类化合物的生物合成,211个差异表达基因参与萜类化合物的生物合成。本研究获得了半枫荷根和叶的转录组信息特征,为今后半枫荷基因功能鉴定、次生代谢途径解析及调控机制的研究提供依据。     

正交法优选调控秦艽次生代谢产物含量变化的条件

为优选可抑制秦艽次生代谢产物含量变化的条件，该文采用3因素4水平正交试验设计方法，共设计16组处理，研究洛伐他汀(MVA途径抑制剂)、膦胺霉素(MEP途径抑制剂)和取样天数对秦艽中马钱苷酸、獐牙菜苷、獐牙菜苦苷、龙胆苦苷4种主要环烯醚萜类化合物含量的影响。结果表明：(1)4种环烯醚萜类化合物含量变化受取样天数影响最大，其次为膦胺霉素浓度，次之为洛伐他汀浓度。(2)以最佳抑制条件处理后，马钱苷酸、獐芽菜苦苷、龙胆苦苷和獐芽菜苷含量分别下降了69%、36%、33%和4%。基于正交法优选的抑制条件，对4种化合物均可抑制。综上所述，可确定调控秦艽次生代谢产物含量变化的最佳抑制条件为膦胺霉素400μmol·L^-1,洛伐他汀50μmol·L^-1,取样天数6 d,该条件为进一步研究MEP和MVA途径在环烯醚萜类化合物代谢合成中的调控机制奠定基础。     

草地贪夜蛾雄性成虫和5龄幼虫的转录组比较分析

【目的】草地贪夜蛾Spodoptera frugiperda是一种新近入侵我国的重要害虫。本研究旨在对草地贪夜蛾雄性成虫和5龄幼虫两个不同发育阶段的转录组进行比较分析。【方法】利用高通量测序技术对草地贪夜蛾雄性成虫和5龄幼虫进行转录组测序和数据组装,并对转录组数据进行功能注释和比较分析。【结果】经de novo组装共获得209 002条转录本,平均长度为687.55 bp,N50为982 bp。共有46 198条(57.43%) unigene在至少一个数据库中获得功能注释,其中1 713条(2.13%) unigene在所有数据库中均能获得注释。在GO数据库中获得205 269条unigene的注释,主要包括68个功能分类;在KEGG数据库中共有3 408条unigene得到注释,涉及277个代谢通路。共鉴定到424个嗅觉相关的基因,并且在雄性成虫和5龄幼虫之间的表达存在差异。通过比较转录组分析,在雄性成虫中鉴定到9 162个上调和6 399个下调差异表达基因(DEGs);功能富集分析发现在上调DEGs中涉及信息素以及信号转导的代谢通路显著富集,而下调DEGs中涉及解毒相关的通路显著富集。【结论】这些转录组数据为探究草地贪夜蛾的生长发育、嗅觉相关功能基因以及候选分子靶标提供资源信息。     

干旱胁迫下丹参数字基因表达谱分析

丹参是中国传统大宗药材,具有重要的药用价值,但其品质受干旱胁迫影响较大,而对丹参响应干旱胁迫的应答机制尚不明确。本研究基于数字基因表达谱技术,对干旱胁迫处理后的丹参植株cDNA文库进行了差异基因表达谱分析。分析结果表明,共有5 740条基因受干旱胁迫诱导差异表达,其中1 143条表达上调,4 597条表达下调。Gene Ontology (GO)功能显著性富集分析表明,干旱胁迫应答基因的生物学功能主要集中在催化活性和结合活性中。Pathway显著性富集分析结果表明,上调基因显著富集在黄酮生物合成途径、核糖体途径、苯基丙酸类生物合成途径、半胱氨酸蛋氨酸代谢途径和植物激素生物合成等途径,下调基因则主要富集在鞘糖脂的合成途径、糖胺聚糖降解途径、卟啉代谢和叶绿素代谢途径上。随机抽取9个差异表达基因进行实时荧光定量PCR验证,其结果与EDG检测结果一致。本研究的分析结果有助于进一步了解丹参响应干旱胁迫的分子机制,为干旱应答功能基因的筛选奠定基础。     

川西獐牙菜SmSLS2基因的克隆、生物信息学及表达分析

裂环马钱子苷合成酶（SLS）是环烯醚萜类化合物合成途径的限速酶。该研究根据川西獐牙菜转录组信息获得SmSLS2基因序列,对该基因进行克隆、生物信息分析、原核表达载体构建及蛋白体外诱导。结果表明SmSLS2 cDNA全长1 566 bp（GenBank编号MH535904）,编码521个氨基酸,推测其蛋白相对分子质量为59.79 kDa,等电点为8.92。结构域分析表明SmSLS2含一个跨膜结构域,SmSLS2的gDNA全长为2 576 bp,含5个外显子和4个内含子。SmSLS2蛋白与其他植物中SLS蛋白具有高度的相似性。进一步将SmSLS2构建到原核表达载体pET-28a-sumo上,转入大肠杆菌BL21（DE3）中,在37℃、0.5 mmol·L^-1 IPTG条件下进行原核表达,诱导出与预测蛋白大小一致的目的蛋白。该研究成功从川西獐牙菜中克隆了SmSLS2,为进一步阐明该基因在环烯醚萜合成途径中的作用和通过基因工程手段改良环烯醚萜类化合物产量提供了理论基础。     

胁迫意大利蜜蜂幼虫肠道的球囊菌的转录组分析

【目的】本研究旨在通过趋势分析对胁迫意大利蜜蜂Apis mellifera ligustica(简称意蜂)幼虫肠道的球囊菌Ascosphaera apis的差异表达基因(DEGs)进行转录组分析。【方法】将纯化的球囊菌孢子配制为1×107孢子/m L的饲料饲喂意蜂3日龄幼虫,利用Illumina Hi Seq 2500平台对球囊菌胁迫的意蜂幼虫肠道c DNA进行测序,将过滤得到的有效读段(clean reads)映射(mapping)到核糖体数据库及意蜂参考基因组,最后将未映射上的reads映射到本课题组组装并注释的球囊菌参考转录组。利用STEM软件对DEGs进行趋势分析。利用WEGO软件对显著表达模式中的DEGs进行GO富集分析。利用Blastall对显著表达模式中的DEGs进行KEGG代谢通路富集分析。最后,通过对随机选取的6个DEGs进行RT-q PCR分析,对RNA-seq数据进行验证。【结果】球囊菌胁迫意蜂幼虫肠道样品的Illumina测序共得到球囊菌的25 454 076条原始读段(raw reads),经过滤得到24 909 820条clean reads,Q30均在93.46%以上。趋势分析结果显示,19 893个DEGs聚类为8个表达模式,其中有12 151个DEGs聚类为3个表现为显著上调趋势的表达模式。GO富集分析结果显示,表现上调趋势的DEGs富集在40个GO term,富集基因数最多的为细胞进程(cellular process)(2 601 unigenes),其次为代谢进程(metabolic process)(2 553 unigenes)和细胞(cell)(2 522unigenes)。KEGG代谢通路富集分析结果显示,上调趋势中的DEGs富集在119个代谢通路上,其中富集基因数最多的是核糖体(ribosome)(213 unigenes),其次为氨基酸生物合成(biosynthesis of amino acids)(154 unigenes)和内质网蛋白加工(protein processing in endoplasmic reticulum)(130unigenes)。共有48个DEGs富集在MAPK信号通路上,聚类分析结果显示,这些DEGs随着胁迫时间的延长表达水平逐渐增强。RT-q PCR结果显示,6个DEGs的表达水平变化趋势与RNA-seq数据一致,证明了本研究中的转录组数据真实可靠。【结论】本研究为在分子水平揭示球囊菌的致病机理提供了重要信息,也为阐释球囊菌胁迫意蜂幼虫过程中的病原-宿主互作机制奠定了基础。     

功能雄性不育茄子差异基因及代谢途径分析

以茄子功能雄性不育系 (S13) 和可育系 (F142) 为材料,对花蕾不同发育时期的细胞学特征进行观察,发现S13花药中环形细胞簇的解体时期比F142延迟2 d,开花当天F142的裂口组织及内绒毡层的细胞解体,而S13中无此现象发生。转录组测序分析表明,F142和S13开花前8天、5天和开花当天的花药共有1 436个差异表达基因 (DEGs) (651个表达上调,785个表达下调)。GO显著性分析表明差异基因主要富集在生物学过程中涉及单细胞生物过程、代谢过程和细胞过程的单基因簇 (Unigene) 较多,分子功能中涉及较多的有催化活性和结合功能;KEGG注释发现,差异基因主要富集在次生代谢产物的生物合成、代谢途径、内质网蛋白质加工、氨基酸的生物合成、碳代谢和植物激素信号转导。选取16 465个基因进行加权基因共表达网络构建分析 (Weighted gene co-expression network analysis,WGCNA),共鉴定到15个基因共表达模块,其中3个为S13在3个花蕾发育时期高度相关特异性模块 (Plum2、Royalblue和Bisque4模块);KEGG富集表明,特异性模块可以富集到苯丙烷类生物合成、光合作用、卟啉与叶绿素代谢、α-亚麻酸代谢、多糖的生物合成和代谢、脂肪酸降解以及戊糖与葡萄糖醛酸的相互转化等相关通路,这些基因可能在S13花蕾发育过程中发挥重要作用。研究结果为进一步探索茄子花药开裂机制提供了一定的参考。     

重度抑郁症发病机制相关基因的生物信息学分析

本文旨在采用生物信息学方法筛选重度抑郁症发病机制相关基因。以美国国家生物技术信息中心(NCBI)网站GEO公共数据库中GSE98793芯片数据为研究对象,用R语言limma包筛选出116个差异表达基因(differentially expressed genes,DEGs),其中上调基因66个,下调基因50个。Gene Ontology(GO)基因功能注释分析结果显示,DEGs主要分布于线粒体内膜、线粒体内,参与铜离子结合、半胱氨酸肽链内切酶活性、白细胞介素-1的细胞应答、蛋白质加工等生物过程。KEGG通路富集分析结果显示,DEGs主要富集于氧化磷酸化反应、帕金森病、非酒精性脂肪肝病、阿尔茨海默病和亨廷顿氏舞蹈病等发病机制。蛋白质相互作用网络分析结果显示,54个DEGs编码的蛋白之间存在相互作用。结合上述多种方法的分析结果,UQCRC1、GZMB、NDUFB9、NSF、SLC17A5、CTSH、NDUFB10、UQCR10、ATOX1、CST7和CTSW共11个关键基因被筛选出来,可作为诊断和治疗重度抑郁症的候选基因。综上,本研究通过对重度抑郁症芯片及其DEGs进行深入分析得到关键基因,为揭示抑郁症的分子机制和临床靶向治疗提供了重要的线索。     

基于生物信息学分析的肝癌核心基因的筛选及验证

本研究基于GEO数据库,选取由慢性乙型肝炎诱导的肝细胞癌芯片数据GSE121248为研究对象,利用GEO2R软件分析数据,筛选出差异表达基因,利用DAVID数据库进行GO分析和KEGG pathway富集分析.利用STRING数据库构建PPI网络,分析筛选核心基因.利用GEPIA对核心基因的表达进行验证,Kaplan Meier Plotter在线分析工具对核心基因与患者生存情况的相关性进行验证.通过上述方法筛选出309个DEGs,其中上调基因94个,下调基因215个.差异基因功能分析显示上调的DEGs主要参与细胞周期和卵母细胞减数分裂通路等途径,下调的DEGs则在补体和凝血级联、代谢途径以及咖啡因代谢途径富集.筛选出15个具有高度关联性的核心基因(BUB1,BUB1B,BIRC5,CCNB1,CCNB2,CDC20,CDK1,KIF-20A,MAD2L1,NCAPG,ZWINT,PBK,BTL,TTK和NUSAP1),它们与肝癌患者的总体生存率具有明显相关性,并为其构建了miRNA调控网络.本研究通过生物信息学方法有效分析了肝细胞癌发生、发展相关的差异表达基因,筛选出15个核心基因,分析其生物学相关功能,以期探索肝细胞癌发病机制,并为临床诊断标志物的改进以及筛选提供一定的理论基础.