MDCK细胞和MDCK-G1细胞感染H1N1流感病毒后病毒滴度和细胞转录组表达差异分析

目的通过分析细胞表达基因及其相关信号通路的差异等信息,探索MDCK和MDCK-G1细胞株在接种H1N1流感病毒后出现病毒滴度和细胞生长趋势差异的内在生物学特性等原因。方法通过Illumina测序平台对2株细胞进行转录组测序(RNA sequencing, RNA-seq)和测序结果生物信息学分析后,选择了17个差异基因进行实时定量聚合酶链反应(quantitative real-time polymerase chain reaction, qRT-PCR)验证。结果 MDCK和MDCK-G1细胞差异基因(differentially expressed genes, DEGs)|log_2(FoldChange)|0和padj≤0.05总数为2 786个。相比于MDCK细胞,MDCK-G1细胞有967个基因的表达显著上升,1 819个基因的表达显著下降。差异基因通过基因本体(Gene Ontology, GO)数据库富集功能注释和京都基因和基因组数据库(Kyoto Encyclopedia of Genes and Genomes, KEGG)进行通路分析发现,2株细胞在免疫调控和细胞生长周期调控上均存在差异,qRT-PCR验证的17个基因中有16个基因的表达趋势与转录组测序结果一致。结论 2株细胞的转录谱存在显著差异。     

家蚕5龄幼虫精巢和卵巢microRNA芯片及转录组比较分析

【目的】本研究旨在通过对家蚕Bombyx mori 5龄幼虫精巢和卵巢组织微小RNA (microRNA, miRNA)基因芯片及转录组进行分析,找到参与家蚕性腺发育相关的miRNA分子及可能的靶基因。【方法】采用新一代高通量测序平台对家蚕5龄幼虫精巢和卵巢(分别定义为Test和Control)进行miRNA基因芯片检测及转录组测序分析,根据P＜0.05且log₂(fold change, FC)≥2的标准,通过比较筛选出Test vs Control的差异表达miRNA;根据q≤0.05且|log₂(fold change)|≥1的标准,通过比较筛选出Test vs Control的差异表达基因 (differentially expressed genes, DEGs);随机选取8个上调和12个下调差异表达miRNA,对其表达及其预测的5个靶基因进行qRT-PCR验证;对DEGs以及差异表达miRNA的靶基因进行KEGG通路富集分析。【结果】从精巢和卵巢样本中(Test vs Control)分别鉴定出68个差异表达miRNA和3 991个DEGs,其中上调和下调miRNA分别为36和32个,上调和下调DEGs分别为2 033和1 958个。差异表达miRNA的qRT PCR验证结果均与芯片数据一致。KEGG通路富集分析结果显示DEGs在新陈代谢及核糖体的信号通路显著富集。对差异表达miRNA在DEGs中的可能靶基因进行预测,结果找到了4组表达趋势相反的miRNA与靶基因：分别是bmo-miR-2774a与LOC101745556;bmo-miR-92b与LOC101735954以及bmo-miR-3266与LOC733130和LOC778467;1组表达趋势一致的miRNA与靶基因：bmo-miR-3321与LOC101744895。5个靶基因的qRT-PCR验证结果与转录组测序结果一致。【结论】本研究获得了家蚕5龄幼虫精巢和卵巢转录组及miRNA芯片数据,筛选并验证了4组差异表达和1组一致表达miRNA及潜在靶基因,为探究家蚕精巢和卵巢发育差异奠定了基础。     

三峡库区消落带野生狗牙根对水淹适生性的转录组分析

为探究狗牙根对水淹环境的适生性机制,取长势一致的狗牙根分为未水淹对照组(A1)和30 cm沉水处理组(A2),A2水淹15 d后两组同时采样,提取样本总RNA作转录组测序。A1和A2测序数据经de novo组装分别获得128031和83363条高质量Unigene;A1和A2对比发现28844条差异表达基因(differentially expressed genes, DEGs),其中11858条DEGs上调表达,16986条DEGs下调表达;功能注释显示,这些DEGs主要集中在转录、翻译、碳水化合物代谢和环境适应等途径;信号通路富集显示,包括植物对病原体的防御、物质代谢、光合作用、形态发育和抗氧化等过程的103个信号通路与狗牙根对水淹生境的适生性有关;随机挑选4个DEGs作qRT-PCR分析,荧光定量结果与测序数据基本一致,说明基于转录组分析DEGs的表达情况较为可靠。该工作或可为进一步挖掘与狗牙根水淹适生性相关的调控基因和功能基因提供基础。     

轮状病毒感染患者外周血单个核细胞的转录组学分析

轮状病毒(Rotavirus,RV)是引起急性肠胃炎的主要病原体,分析RV感染患者的人外周血单个核细胞(Peripheral blood mononuclear cell,PBMC)中差异表达基因(Differentially expressed genes,DEGs)有利于探讨人PBMC在清除RV中的作用。为此,本研究采集2019年2月-2019年6月长春儿童医院中RV感染患者和健康儿童血液,分离PBMC,通过转录组测序(RNA sequencing,RNA-seq)技术比较RV感染患者与健康儿童之间的RNA表达图谱,借助基因本体论(Gene Ontology,GO)数据库功能富集分析、京都基因与基因组百科全书(Kyoto Encyclopedia of Genes and Genomes,KEGG)、Reactome通路富集分析DEGs,使用实时荧光定量PCR(Real-time quantitative PCR,qPCR)技术进行验证。结果显示,与健康对照组相比,RV感染轻症患者PBMC中有1619个DEGs;重症患者PBMC中有2816个DEGs,主要与干扰素(Interferon,IFN)反应、中性粒细胞、溶酶体、核小体、染色质等相关。qPCR验证轻症患者干扰素刺激基因(IFN-stimulated genes,ISGs)15表达上调,白介素(Interleukin,IL)1β表达下调;重症患者IL15、ISG15表达上调,IL1β表达下调,与转录组结果相一致。本研究提示,RV感染可能激活人I型和II型IFNs反应抵御病毒感染,但也会抑制溶酶体相关基因,对细胞自噬过程产生影响。     

西方蜜蜂工蜂蛹响应低温胁迫的差异表达基因分析

【目的】本研究旨在从转录组水平筛选和分析西方蜜蜂Apis mellifera工蜂不同时期蛹对低温胁迫反应的差异表达基因(DEGs)。【方法】将西方蜜蜂工蜂封盖后3 d预蛹以及封盖后6和9 d蛹分别置于低温环境(20℃)和最适发育温度(35℃)中4 h,分别作为处理组和对照组,通过转录组学技术筛选低温处理组与对应的对照组之间的DEGs,并进行GO功能分类和KEGG通路分析。利用RT qPCR分别对封盖后3 d预蛹以及封盖后6和9 d蛹的8, 6和5个DEGs的表达谱进行验证。【结果】与对照组相比,封盖后3 d预蛹以及封盖后6和9 d蛹受低温胁迫后的DEGs分别有220, 50和26个;GO功能分类发现,DEGs富集数最多的条目为代谢进程、细胞进程、催化活性和结合,封盖后3 d预蛹的DEGs在生物学进程调控、细胞部分和细胞器等有较多富集。KEGG通路分析显示,处理组和对照组间西方蜜蜂各日龄DEGs在整体概述图、氨基酸代谢、信号转导、运输和分解代谢有富集。封盖后3 d预蛹以及封盖后6和9 d蛹受低温胁迫后共有的DEG 3-磷酸肌醇依赖性蛋白激酶1基因PDK1上调表达,同时富集在自噬-动物、mTOR和FoxO信号通路。低温胁迫后,封盖后3 d预蛹的胰岛素受体底物1-B基因IRS1-B和Kruppel同源物1基因Kr-h1上调表达,而核激素受体FTZ-F1基因Ftz-F1与蜕皮启动激素基因Eth显著下调表达,说明封盖后3 d预蛹响应低温细胞自噬和蜕皮受到抑制程度更大。在低温胁迫后封盖后6 d蛹中与昆虫角质层着色与免疫相关的酪氨酸羟化酶基因TyHyd下调表达;与对照组相比,在低温胁迫后封盖后9 d蛹中DEGs最少,说明在3个蛹期中,其受低温的影响较小。【结论】本研究测定了西方蜜蜂3个不同发育阶段蛹响应低温的DEGs,结果显示大部分DEGs为阶段特有,说明西方蜜蜂不同发育阶段响应低温的机制不同。一些共有DEGs以及阶段特有DEGs的功能研究和其作用机制是进一步研究蜜蜂对低温响应机制的重点内容,为探究蜜蜂蛹期响应低温胁迫的分子机制提供了基础数据。     

宁乡猪皮下脂肪与肌内脂肪组织转录组差异分析

为了比较分析宁乡猪皮下脂肪(subcutaneous fat,SAF)和肌内脂肪(intramuscular fat,IMF)组织脂肪沉积的分子机制,本文利用RNA-seq技术鉴定和分析宁乡猪SAF和IMF组织中差异基因表达谱。选取6头250日龄健康状况良好、种内个体体重相近(约85 kg)的雄性宁乡猪为实验材料,采集SAF与IMF组织样品。通过对两个脂肪组织转录组测序并进行GO(Gene Ontology)功能注释及KEGG(Kyoto Encyclopedia of Genes and Genomes)信号通路富集分析,得到与脂肪沉积和脂质代谢相关的差异基因。为验证测序数据结果的可靠性,本研究随机选取6个差异基因进行qRT-PCR验证。结果表明,宁乡猪SAF和IMF组织中有2406个基因差异表达,其中1422个基因表达上调,984个基因表达下调。通过GO功能注释分析发现,差异表达基因主要通过参与类固醇生物合成、不饱和脂肪酸的生物合成、甘油磷脂代谢及自噬途径等与脂质代谢相关途径,调控宁乡猪SAF和IMF的沉积。KEGG通路富集结果显示,差异基因主要富集在脂质结合、脂肪酸代谢过程、甘油...     

侧孢短芽孢杆菌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的拮抗下抑制了叶枯病菌的生长,其转录组特征也发生显著变化,差异基因主要涉及代谢、细胞过程以及遗传信息处理等途径的相关通路。     

基于高通量测序技术分析番鸭呼肠孤病毒感染雏番鸭肠道转录组学的变化

番鸭呼肠孤病毒(Muscovy duck reovirus,MDRV)感染雏番鸭可导致肠道黏膜受损,造成肠黏膜免疫功能低下甚至丧失。为寻找番鸭呼肠孤病毒感染雏番鸭肠道组织免疫相关的差异表达基因,本研究建立番鸭呼肠孤病毒自然发病模型,在感染后3d、6d、21d采集同居感染组和对照组试验雏番鸭的空肠,对其进行高通量测序,对差异表达基因进行GO功能分类和KEGG数据库分析。结果显示,感染后3d、6d和21d分别筛选出2 297、5 860和3 721个差异表达基因;GO分子功能结果显示,免疫相关差异表达基因主要和免疫球蛋白的产生、T细胞活化和单核细胞趋化性等有关;KEGG通路富集结果显示,差异表达基因主要涉及在吞噬体、细胞溶质DNA传感途径、细胞因子-细胞因子受体相互作用、Toll样受体信号通路、Jak-STAT信号通路和RIG-I受体信号通路;荧光定量RT-PCR对差异表达基因进行验证,其结果与转录组学结果基本一致。本研究为进一步阐明MDRV感染导致雏番鸭肠道损伤和黏膜免疫抑制的机制奠定基础。     

西方蜜蜂工蜂蛹中期响应低温胁迫的差异表达基因分析

【目的】本研究旨在分析西方蜜蜂Apis mellifera工蜂蛹中期可能存在的低温应对机制以及低温对中期蛹发育的不利影响。 【方法】对西方蜜蜂工蜂8 d封盖子进行低温(20℃)处理96 h(T),以未经低温胁迫的8 d封盖子为对照(CK),通过转录组测序(RNA-seq)并筛选差异基因(differentiallyexpressed genes, DEGs);对DEGs进行GO分类和KEGG通路分析。利用RT-qPCR对随机选取的5个DEGs的表达量进行验证。【结果】西方蜜蜂8 d封盖子低温(20℃)胁迫96 h的DEGs有1 101个;GO分类发现DEGs富集数最多的条目为代谢过程(142个DEGs)和细胞过程(142个DEGs),其次是结合(131个DEGs),催化活性(120个DEGs)和单有机体过程(118个DEGs);KEGG通路分析结果显示,DEGs显著富集于与氧化损伤密切相关的过氧化物酶体通路、D-谷氨酰胺代谢通路、D-谷氨酸代谢通路和谷胱甘肽代谢通路;此外,与激素调控相关的昆虫激素代谢通路、mTOR信号通路和FOXO信号通路上也有DEGs显著富集。随机挑选的5个DEGs的RNA-Seq与RT-qPCR结果一致。【结论】本研究发现低温状态下西方蜜蜂中期蛹主要通过体内激素调控发育进程,使其减速或停止发育以应对低温;低温通过影响其抗氧化酶表达量进而可能积累氧化损伤,从而对羽化后蜜蜂产生影响。本研究结果有助于理解发育阶段温度生态幅狭窄的昆虫对于低温的应对机制及低温对其的不利影响。     

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

本研究基于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个核心基因,分析其生物学相关功能,以期探索肝细胞癌发病机制,并为临床诊断标志物的改进以及筛选提供一定的理论基础.     

球囊菌胁迫中华蜜蜂幼虫肠道过程中病原的转录组学研究

【目的】本研究利用RNA-seq技术对球囊菌胁迫的中华蜜蜂(中蜂)幼虫肠道进行深度测序,经趋势分析得到差异表达基因(DEGs)的显著表达模式,进而对胁迫过程中的球囊菌进行转录组学分析。【方法】利用Illumina HiSeq 2500平台对球囊菌胁迫的中蜂幼虫肠道进行深度测序,并利用相关软件进行了深入分析。最后,通过RT-qPCR对RNA-seq数据进行了验证。【结果】本研究共得到球囊菌的41133932条高质量clean reads。22865个DEGs共聚类为8个基因表达模式,其中,16769个DEGs聚类为2个显著上调趋势与2个显著下调趋势。GO富集分析结果显示,显著上调与显著下调趋势中的DEGs分别富集于40与37个GO term,基因富集数最多的为细胞进程(2486 unigenes)。KEGG代谢通路(pathway)富集分析结果显示,显著上调与显著下调趋势中的DEGs分别富集于119和112个pathway,基因富集数最多的分别是氨基酸生物合成(127 unigenes)与核糖体(98 unigenes)。进一步分析表明球囊菌在胁迫中蜂幼虫肠道的过程中通过提高物质合成促进其增殖,而宿主通过抑制球囊菌的蛋白合成抵御病原入侵。富集在MAPK信号通路的11个DEGs的表达水平随着胁迫时间的延长而逐渐下降,推测中蜂幼虫通过抑制该通路而阻遏球囊菌增殖。【结论】本研究不仅为揭示白垩病过程中的球囊菌-中蜂幼虫互作提供了重要信息,也为阐明不同抗性蜂种的球囊菌抗性差异奠定了基础。     

肺形侧耳低温胁迫时期的转录组分析

本文以肺形侧耳栽培菌株X57为研究对象,利用高通量RNA测序技术对4℃低温处理0h、6h和12h后的肺形侧耳进行基因表达分析来探讨肺形侧耳低温应答机制。分析结果筛选出低温处理6h差异表达基因742个,其中上调表达基因374个,下调表达基因368个;低温处理12h差异表达基因1 489个,其中上调表达基因占53%,下调表达基因占47%。Gene Ontology（GO）功能聚类分析表明,差异表达基因主要富集在结合、催化分子功能组和代谢过程、细胞过程生物学过程中。KEGG功能富集分析结果显示,差异基因主要富集在氨基酸、核糖体和类固醇生物合成,及氮类物质代谢的通路上,富集到与低温胁迫相关通路MAPK signaling pathway-yeast上的基因表达随低温处理时间的增加呈上调趋势。已报道HOG-MAPK通路是MAPK途径研究较为明确、信号传递单一的真菌低温胁迫中重要的通路,本文运用生物信息学软件构建肺形侧耳的HOG-MAPK通路,并利用荧光定量PCR对相关基因表达进行验证,结果显示以低温处理0h为参照,随着低温胁迫时间增加通路上大部分基因的表达量持续上升,且差异显著,与RNA‐Seq分析结果一致。本文通过全面分析肺形侧耳低温胁迫时期基因表达情况,为进一步研究肺形侧耳低温胁迫相关重要基因的功能鉴定与信号通路调控机理提供基础。     

RNA-seq用于获得共轭亚油酸对贵妃鸡肌内脂肪代谢调控的关键基因

通过转录物组测序获得在贵妃鸡基础日粮中添加共轭亚油酸(CLA)对肌内脂肪代谢的差异表达基因,经生物信息学分析获得相关的信号通路及可能发挥重要作用的候选基因,为CLA对肌内脂肪沉积的分子机制奠定基础。本研究选用55日龄健康的贵妃鸡为试验动物,在基础日粮中添加CLA 0%、1%和2%,预饲期1周,正饲期6周。屠宰采集胸肌组织进行转录物组测序,对测序数据进行差异表达分析,差异表达基因GO功能和差异表达基因KEGG通路富集分析,筛选出与胸肌脂类代谢相关的差异表达基因,利用qRT-PCR对差异表达基因进行验证。结果显示,共获得1 065个差异表达基因,其中上调基因703个,下调基因362个。GO富集结果显示,差异表达基因主要富集在生物过程的细胞过程、单一生物过程、生物调节和代谢过程。KEGG信号通路富集显示,差异表达基因显著富集在黏着斑、不饱和脂肪酸生物合成、脂肪酸生物合成和类固醇生物合成等信号通路中,发现11个主要与肌内脂肪代谢相关的候选基因,分别是FADS1、FADS2、ELOVL5、ACOX2、SLC27A1、FABP5、LPL、LOC107050163、ENSGALG00000030996、ENSGALG00000005043和ENSGALG00000048882。并随机选取6个基因进行qRT-PCR验证,其相对表达量变化趋势与测序结果一致。本研究筛选到CLA影响贵妃鸡胸肌脂类代谢相关的差异表达基因,并对11个主要参与脂肪代谢相关的基因进行分析,为揭示CLA调控肌内脂肪沉积的分子机制奠定基础。     

黄瓜二倍体及其同源四倍体果实转录组分析

果形和糖含量分别是与瓜类蔬菜作物外观、营养及风味密切相关的重要品质性状。为探究黄瓜倍性差异导致果形及糖含量变化的分子机制，以果形和糖含量差异显著的黄瓜二倍体（长果/低糖）和四倍体（短果/高糖）为试材，分别取开花当天和花后9 d的果实进行转录组分析。结果表明，在二、四倍体花后9 d果实中共鉴定到1 874个差异表达基因（DEGs），其中上调基因818个和下调基因1 056个；GO分析显示，DEGs显著富集在“DNA复制”（GO:0006260）、“DNA代谢过程”（GO:0006259）、“蛋白质复合体”（GO:0032993）、“染色体”（GO:0005694）、“蛋白质二聚活性”（GO:0046983）等GO条目上；KEGG分析表明，DEGs富集在“植物激素信号转导”（csv04075）、“淀粉与蔗糖代谢”（csv00500）等代谢通路上；参与“植物激素信号转导”相关基因AUX1、DELLA、B-ARR、TCH4上调表达，GH3、GID1、PP2C、CYCD3下调表达，以及参与“淀粉与蔗糖代谢”相关基因SPS、SUSY上调表达，可能与果形指数减小和果实糖含量增加有关。荧光定量PCR(...     

干酪乳杆菌刺激小鼠肠道的转录组学分析

干酪乳杆菌(Lactobacillus casei)刺激小鼠肠道后，利用高通量测序技术对干酪乳杆菌饲喂组和空白组小鼠的肠道组织进行分析，以期查询和验证干酪乳杆菌对肠道免疫的影响。转录组数据的生物信息学分析发现差异表达基因共751个，通过GO富集分析发现有14个基因与细胞免疫相关，聚焦在T细胞激活、细胞分化、免疫调节负调控等功能上。KEGG通路富集显示聚集在PPAR信号通路、B细胞受体信号通路和趋化因子信号通路。对基因的基本特性分析结果显示，14个基因分别定位在10条染色体上，蛋白的分子量介于11.37～83.45 kDa之间，等电点pI 4.42～11.36,均为不稳定脂溶性蛋白，并具有相关功能结构域。通过保守基序与结构分析发现Motif分布具有保守性，大部分基因含有2个内含子。qRT-PCR验证结果表明，14个基因中有6个基因(Ces1d、Lzts1、Paqr7、Aloxe3、Zbtb16、OX40)在不同时间的处理下整体表达水平较高。验证结果与转录组测序结果一致，且这些基因功能与细胞免疫相关，该结果为研究干酪乳杆菌对机体的免疫作用效果提供一定的理论依据。     

沙棘刺和芽的组织解剖与转录组分析北大核心CSCD

该研究通过组织切片分析比较了沙棘刺与芽的组织染色差异,通过转录组测序分析比较了刺和芽的基因表达谱差异,通过qRT-PCR验证候选基因的表达水平,为进一步通过基因工程手段调控沙棘刺的发育提供候选基因。结果表明:(1)沙棘芽横切面的组织边界层次分明,而沙棘刺横切面的组织边界模糊,且刺尖中心的红色区域较大,木质素自发荧光结果表明刺中木质化程度显著增强。(2)转录组分析组装得到81560个单基因(unigenes),其中包含9385个差异表达基因(DEGs)。(3)GO和KEGG富集分析表明,富集程度最高的基因主要与细胞壁、表皮、气孔发育有关,尤其是与木质素合成途径相关;芽刺相关DEGs分析结果表明,木质素合成相关DEGs共20个,角质、木栓质和蜡质生物合成相关DEGs共29个,植物表皮发育和气孔复合体发育的相关DEGs相同共14个;结合文献分析推测,部分热头蛋白基因参与控制沙棘刺顶端分生组织消失,木质素合成相关基因则与刺硬化过程密切相关。(4)qRT-PCR分析验证显示,与芽相比,挑选的6个木质素合成相关DEGs中有3个分别上调了9.5倍、41.1倍和13.7倍,3个下调至5%、14%和24%;挑选的4个角质、木栓质和蜡质生物合成相关DEGs中有2个分别上调3.6倍和22.1倍,另外2个下调至3%和6%;挑选的4个表皮和气孔发育相关DEGs均下调至6%、13%、3%和4%;关键基因的相对表达水平与转录组测定的TPM值可相互印证,表明验证结果与转录组分析结果一致。研究发现,沙棘刺和芽具有明显的形态和组织差异,二者间的DEGs主要与木质素合成以及细胞壁、表皮和气孔发育有关。     

鲤在低温胁迫下肝胰腺转录组测序分析

检测低温耐受与低温敏感鲤品种低温胁迫下的差异表达基因,探究鲤低温适应的分子应答过程。低温胁迫后对松浦镜鲤及荷包红鲤肝胰腺进行转录组测序,并对表达差异基因进行功能、通路富集分析。通过差异表达分析,共获得10 521个差异表达基因,其中5 246个基因仅在低温耐受品种发生表达变化。功能、通路富集分析将差异基因显著富集在包括糖酵解/糖质新生代谢途径在内的144个代谢途径。选取5个基因进行qRT-PCR验证,定量结果与RNA-seq结果相关性达0.89。研究结果表明低温胁迫会显著改变鲤肝胰腺的转录活动,编码糖酵解/糖质新生途径的多个关键酶基因在低温耐受品种的表达量显著高于低温敏感品种,推测这些基因功能有助于其长期低温适应。     

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

本文旨在采用生物信息学方法筛选重度抑郁症发病机制相关基因。以美国国家生物技术信息中心(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进行深入分析得到关键基因,为揭示抑郁症的分子机制和临床靶向治疗提供了重要的线索。     

家兔前体脂肪细胞分化不同时期基因表达谱分析

脂肪的过度积累严重危害人类健康。前体脂肪细胞分化是脂肪发育的关键过程,研究前体脂肪细胞分化相关基因的表达有助于认识脂肪沉积的机理。尽管家兔是一种理想的研究脂肪发育的动物模型,但是针对其前体脂肪细胞分化不同时期基因表达谱的研究鲜见报道。本研究通过诱导家兔前体脂肪细胞分化,在分化第0 d、3 d和9 d收集脂肪细胞,利用转录组测序(RNA-seq),在分化第3 d样本与第0 d样本的比较中筛选出1352个差异表达基因(differentially expressed genes, DEGs),在分化第9 d样本与第3 d样本的比较中筛选出888个DEGs。GO (gene ontology)功能富集和KEGG (kyoto encyclopedia of genes and genomes)通路分析发现,0~3 d分化期上调的DEGs显著富集在PPAR信号通路和PI3K-Akt信号通路上,3~9d分化期上调的DEGs显著富集到与细胞周期调控有关的GO条目和KEGG信号通路,0~3d和3~9d阶段特异上调的DEGs可能分别作用于细胞质和细胞核。通过DEGs的蛋白-蛋白互作(protein-protein interaction, PPI)网络分析发现,筛选出的核心节点(hub node)基因可能通过调控细胞周期而影响家兔前体脂肪细胞分化。     

生物膜形成中间状态下副溶血弧菌的基因转录谱分析

【目的】探究生物膜形成中间状态下副溶血弧菌的差异基因表达情况,为今后研究生物膜形成调控机制提供基因信息。【方法】以非生物膜形成条件下为参照,采用Illumina HiSeq测序平台进行转录组测序（RNA sequencing,RNA-seq）研究,分析生物膜形成中间状态下副溶血弧菌的基因表达情况,并采用实时定量PCR （quantitative real-time PCR,qPCR）进行验证。【结果】本研究共获得979个差异显著性表达基因（differentially expressed gene,DEG）,其中下调基因379个,上调基因600个。基因本体（gene ontology,GO）分类分析结果显示,共有363个DEGs注释到分子功能、细胞组分和生物学过程3个一级分类和30个二级分类;京都基因和基因组百科全书（Kyoto encyclopedia of genes and genomes,KEGG）代谢途径分析结果显示,共有706个DEGs归到37个代谢通路中（Q value<0.05）,差异表达基因主要集中在细胞代谢和转运通路上;蛋白相邻类的聚簇（clusters of orthologous groups,COG）分类结果显示,有888个DEGs可归为20个类别,涉及DEGs最多的为氨基酸转运与代谢、一般功能预测基因、能量产生与转换以及未知功能基因。qPCR验证挑选的DEGs变化趋势均与RNA-seq的结果一致。此外,从转录组数据中共筛选出10个c-di-GMP代谢相关基因、1个侧生鞭毛蛋白基因（lafA）、13个极生鞭毛合成相关基因、6个荚膜多糖合成相关基因、6个胞外多糖合成相关基因、5个IV型菌毛合成相关基因、膜融合蛋白（membrane fusion protein,Mfp）基因（cpsQ-mfpABC）、14个III型分泌系统1（T3SS1）相关基因、14个Vp-PAI基因（1个tdh2和13个T3SS2基因）、3个VI型分泌系统1（T6SS1）相关基因、6个T6SS2基因、45个推定调控子基因和15个推定的外膜蛋白基因。【结论】生物膜形成引起副溶血弧菌基因表达谱发生明显变化,差异表达基因中包含生物膜形成相关基因、关键毒力基因和许多推定调控子基因等,这为后续研究生物膜形成调控机制提供重要信息。