白化菠萝蜜茎2个蛋白激酶和4个转录因子家族分析

为探究蛋白激酶(PKs)和转录因子(TFs)在白化菠萝蜜(Artocarpus heterophyllus)幼苗茎次生生长中的表达变化,基于转录组数据对其差异表达基因(DEGs)进行预测及分类,并对挑选出的2个PKs和4个TFs家族构建系统进化树。结果表明,胞质类受体激酶(RLCK)-VIII家族的DEGs上下调表达各4个,亮氨酸富集重复类受体激酶(LRR-RLK)-X家族Xa和Xb-2分支中的DEGs均下调表达,Xb-1中的均上调,TCP家族的20个DEGs中有15个上调表达,zf-HD和GRF家族中的大多数DEGs上调表达,Alfin-like家族中的DEGs均下调表达。因此,这表明6个家族可能在菠萝蜜茎的次生生长过程和应对非生物胁迫中发挥重要作用。     

自毒胁迫对甜瓜种子萌发与幼苗保护酶活性和转录组的影响

甜瓜自毒作用是导致其栽培实践中连作障碍严重的关键原因之一。该实验以甜瓜植株水浸提液处理模拟甜瓜自毒胁迫,通过测定自毒条件下甜瓜种子萌发、幼苗根系保护酶活性和MDA含量的变化及转录组分析,以探讨甜瓜自毒作用机理。结果显示:(1)甜瓜自毒胁迫总体上抑制了甜瓜种子萌发和后续生长,胁迫处理的阈值为0.03 g/mL。(2)甜瓜植株水浸提液处理后其幼苗根系保护酶活性和MDA含量变化非常剧烈,SOD活性表现为降-升-降趋势,POD活性先降后升,CAT活性先升后降,MDA含量持续增加。(3)转录组分析结果共鉴定出2 599个差异表达基因(differentially expressed genes, DEGs),胁迫2 d后共产生2 251个DEGs,显著多于胁迫4 d后的329个DEGs;且2 d与0 d对比有923个DEGs上调,1 328个DEGs下调,表明甜瓜幼苗在基因水平对自毒胁迫产生了积极响应。(4)相关生物信息学分析表明,自毒胁迫导致的DEGs主要与苯丙烷代谢、活性氧代谢、光合作用和植物激素信号转导有关,同时也涉及到渗透调节、膜和蛋白保护等过程;且这些DEGs主要富集于2 d与0 d,说明自毒胁迫条件下幼苗光合作用发生了改变,且这种变化主要发生在胁迫早期。(5)对6个甜瓜自毒胁迫密切相关DEGs(AP2-2、bZIP1、bZIP2、AP2-1、bHLH和HIS)的qRT-PCR分析显示,AP2-1、bHLH和HIS等3个基因在自毒胁迫2 d时表达量出现峰值,对照组AP2-2基因在2 d时表达量达到峰值,结果与转录组测序分析结果一致。研究表明,甜瓜自毒胁迫引起了其植物细胞的异常,进而会对生长和器官结构产生不良影响并诱发了大量与刺激或胁迫相关的基因差异表达,且甜瓜自毒胁迫下短期响应的基因或转录因子数量明显多于长期响应的数量,暗示甜瓜幼苗可以对自毒胁迫做出快速响应。     

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

【目的】本研究旨在从转录组水平筛选和分析西方蜜蜂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的功能研究和其作用机制是进一步研究蜜蜂对低温响应机制的重点内容,为探究蜜蜂蛹期响应低温胁迫的分子机制提供了基础数据。     

低温胁迫下意大利蜜蜂预蛹的差异表达基因趋势分析

【目的】蜜蜂是典型的具有发育狭温性的全变态昆虫。本研究以对低温最敏感的意大利蜜蜂Apis mellifera ligustica预蛹为研究对象,通过低温胁迫不同时间的差异表达基因(differentially expressed genes, DEGs)趋势分析,探讨低温胁迫对蜜蜂发育影响的关键基因。【方法】对3日龄意大利蜜蜂封盖子预蛹进行20℃低温胁迫18 h(T18)和36 h(T36),以未经低温胁迫的预蛹为对照(CK),通过Illumina HiSeq~(TM)平台进行转录组学测定。利用Short Time-series Expression Miner(STEM)软件对2个处理组与对照组比较共有的差异表达基因进行趋势分析,再进一步对显著富集趋势模式中富集的差异表达基因进行GO分类和KEGG pathway分析。利用RT-qPCR对随机挑选的5个DEGs的表达模式进行验证。【结果】对检测到1 062个T18 vs CK和T36 vs CK共有的DEGs进行趋势分析,发现3个显著的基因表达模式,包括2个上调表达模式(Profile 6,有539个基因;Profile 7,有271个基因),1个下调表达模式(Profile 1,有183个基因)。对3个显著富集趋势模式DEGs分别进行GO富集分析和KEGG pathway分析,在Profile 6中找到同时富集在FoxO信号通路和寿命调节信号通路上的胰岛素样肽基因ILP和叉头蛋白O基因FoxO持续上调,说明低温胁迫会影响蜜蜂预蛹蜕皮激素信号传递;在Profile 7中CREB结合蛋白基因CBP持续上调,蜜蜂预蛹受到低温胁迫会影响细胞发育;在Profile 1中细胞色素P450基因CYP450 306a1 (phm)和WNT1显著下调,说明化蛹时受到低温胁迫会影响蜕皮激素的合成。通过RT-qPCR分析挑选的5个基因的表达结果和高通量测序结果一致。【结论】通过对蜜蜂响应低温胁迫的差异表达基因趋势分析发现,蜜蜂胰岛素和蜕皮激素共同调控的FoxO可能是低温胁迫抑制蜜蜂化蛹的一个关键基因。本研究为探索低温胁迫影响蜜蜂发育变态的分子调控机制奠定了基础。     

盐度胁迫对幼体大海马基因转录表达的影响

为探究大海马(Hippocampus kuda Bleeker)幼体在高盐、低盐胁迫条件下的基因表达水平的变化规律, 对实验条件下的大海马幼体的肝脏样品进行了转录组测序。对照组(CK, 25‰)、高盐(HS-test, 31‰)和低盐 (LS-test, 17‰)胁迫组共获得71794个单基因簇(Unigenes), N50为1780 bp, 平均长度为820.71 bp。高盐胁迫组与对照组比较, 共获得2740个差异表达基因 (DEGs), 其中495个DEGs上调, 2245个DEGs下调; 与对照组相比, 低盐胁迫组共获得3715个DEGs, 其中1854个DEGs上调, 1861个DEGs下调。高/低盐胁迫组DEGs经 KEGG 数据库富集发现, 高/低盐度胁迫均能导致大海马幼体体内氨基酸代谢、免疫代谢、能量和脂肪酸代谢相关基因受到影响。其中, 低盐胁迫时能量代谢和氨基酸代谢的相关基因显著上调, 高盐胁迫时脂肪酸代谢的相关基因显著下调, 而高/低盐胁迫时免疫代谢的相关基因都显著上调。从经过盐度胁迫的大海马幼体的肝脏转录组中分别筛选到免疫相关基因Gst、Hsp70、Hsp90、Sod、Bcl-2、Gadd45α、Tcrβ、Tap2和Traf3, 脂肪酸代谢相关基因Fadsd6、Fas、Sqle、Cyp51、Elovl6和Slc27a6, 能量代谢相关基因Vlcad、Pdha1、Mdh1、Idh3b、G6pd和Sdhd, 及一些氨基酸代谢相关基因Gldc、Atp6v1e1、Sms、Fadh、Asl、Ass1和Glud1等, 可作为大海马幼体响应环境盐度变化应激的候选基因。研究结果为盐度胁迫下大海马幼体的稳态调控机制的研究奠定了一定基础, 有助于在养殖实践中预防极端的盐度改变对大海马幼体所造成的影响。     

罗氏沼虾急性低温应激响应的转录组分析

为探究罗氏沼虾(Macrobrachium rosenbergii)急性低温应激响应的基因表达模式,研究以罗氏沼虾成虾为研究对象,分别设置实验组(16℃)和对照组(24℃),实验组从24℃急性降温(2℃/1h)至16℃后的0、1h、3h、6h、12h、24h和回温至24℃后共7个时间点采集肝胰腺组织进行转录组测序分析。差异表达基因(DEGs)筛选结果显示,胁迫3h与回温后的DEGs数量(5062个)几乎是1h与回温后的1.5倍(3516个),随着胁迫时间的延长,各时间点与回温后的DEGs数量逐渐降低。KEGG富集发现, DEGs显著富集在溶酶体、淀粉和蔗糖代谢、抗原加工与呈递等途径中。此外,罗氏沼虾在急性低温应激下,黏着斑、ECM-受体互作、细胞色素P450对异生物质的代谢、谷胱甘肽代谢、氧化磷酸化、p53信号通路等相关基因也发生了显著变化。WGCNA分析发现各组间的共有DEGs聚类在与细胞功能及免疫相关模块和与能量物质代谢相关模块上。另外,花生四烯酸代谢信号通路中的Cytochrome P450 2L1-like基因、谷胱甘肽代谢信号通路中的NADP-specific isocitr...     

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

【目的】本研究利用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的表达水平随着胁迫时间的延长而逐渐下降,推测中蜂幼虫通过抑制该通路而阻遏球囊菌增殖。【结论】本研究不仅为揭示白垩病过程中的球囊菌-中蜂幼虫互作提供了重要信息,也为阐明不同抗性蜂种的球囊菌抗性差异奠定了基础。     

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

【目的】本研究旨在通过趋势分析对胁迫意大利蜜蜂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数据一致,证明了本研究中的转录组数据真实可靠。【结论】本研究为在分子水平揭示球囊菌的致病机理提供了重要信息,也为阐释球囊菌胁迫意蜂幼虫过程中的病原-宿主互作机制奠定了基础。     

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

脂肪的过度积累严重危害人类健康。前体脂肪细胞分化是脂肪发育的关键过程,研究前体脂肪细胞分化相关基因的表达有助于认识脂肪沉积的机理。尽管家兔是一种理想的研究脂肪发育的动物模型,但是针对其前体脂肪细胞分化不同时期基因表达谱的研究鲜见报道。本研究通过诱导家兔前体脂肪细胞分化,在分化第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)基因可能通过调控细胞周期而影响家兔前体脂肪细胞分化。     

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

【目的】本研究旨在分析西方蜜蜂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结果一致。【结论】本研究发现低温状态下西方蜜蜂中期蛹主要通过体内激素调控发育进程,使其减速或停止发育以应对低温;低温通过影响其抗氧化酶表达量进而可能积累氧化损伤,从而对羽化后蜜蜂产生影响。本研究结果有助于理解发育阶段温度生态幅狭窄的昆虫对于低温的应对机制及低温对其的不利影响。     

急性氨氮胁迫下大口黑鲈的肝脏转录组特征分析

为研究大口黑鲈(Micropterus salmoides)受到氨氮胁迫时基因表达的变化规律,文章采用Illumina平台的HiSeq测序策略分析了氨氮胁迫12h和24h后大口黑鲈肝脏的基因表达谱。氨氮胁迫后在大口黑鲈肝脏共获得111.66 Gb的有效数据,组装获得133486个单基因簇(Unigenes), N50为1000 bp。通过比较转录组分析在2个时间点共获得2072个差异表达基因,其中氨氮胁迫12h获得1516个差异表达基因, 907个上调, 609个下调;氨氮胁迫24h获得556个差异表达基因,其中330个上调, 226个下调。GO功能注释分析结果表明,在氨氮胁迫12h差异基因富集数目明显多于24h,但在“氧化还原酶活性”条目中富集的差异基因数目则是氨氮胁迫24h多于12h。此外, KEGG功能富集分析发现,氨氮胁迫对大口黑鲈肝脏氧化应激、细胞自噬和凋亡相关等途径产生影响。选取10个差异表达基因进行qPCR验证,其表达水平与转录组差异基因分析结果一致,证明了测序结果的可靠性。其中,与氧化应激相关基因缺氧诱导因子1α (HIF1α)、生长停滞DNA损伤可诱导蛋白β (Gad...     

烟草幼苗响应弱光胁迫的转录组分析

为了研究烟草幼苗对弱光胁迫反应的分子机制,以'云烟87'为材料,构建全光照和弱光处理条件下大十字期烟苗cDNA文库,利用Illumina测序技术进行转录组测序,筛选差异表达基因.结果表明:借助RNA-seq技术共筛选到2 956个DEGs,其中弱光相对于全光照表达上调的DEGs有691个,表达下调的DEGs为2 265个;通过GO功能分类和KEGG富集分析,将差异基因划分为30个功能类别和113个代谢途径,包括代谢进程、碳水化合物结合、核酸结合转录因子活性、水解酶活性、转录因子活性、催化活性、内质网蛋白质加工、植物-病原体互作、苯丙烷生物合成、苯丙氨酸代谢、胁迫反应、激素信号转导以及光合作用等.本研究将为揭示烟草对弱光胁迫反应的分子机制提供理论依据.     

不同油葵品种对盐碱地根际土壤酶活性及微生物群落功能多样性的影响

通过盆栽试验,研究了新葵4号、新葵6号、新葵10号和美国矮大头4个油葵品种对盐碱地土壤理化、酶活性和微生物群落功能多样性的影响,以期筛选出更适宜改善盐碱地土壤质量的油葵品种。结果表明,种植新葵6号对降低盐碱地根际土壤pH值、提高土壤全氮含量和蔗糖酶活性的效果最为显著,新葵4号对提高根际土壤碱解氮、速效磷、速效钾含量以及脲酶和磷酸酶活性的效果最为显著;种植这4个品种的油葵均能显著(72 h,P0.05)提高盐碱地根际土壤微生物对31种碳源的平均利用率(Average well color development,AWCD),并呈现如下规律:新葵4号新葵6号美国矮大头新葵10号CK。种植这4个油葵品种均不同程度地提高了盐碱地根际土壤微生物群落的Shannon多样性指数(H)、Shannon优势度指数(D)和碳源利用丰富度指数(S),并呈现出相似的规律:新葵4号根际土的微生物多样性指数最大,而CK的最小,且显著高于CK。主成分分析表明,种植油葵改善了盐碱地根际土壤微生物的群落组成;碳水化合物、氨基酸、羧酸类化合物和聚合物是盐碱地土壤微生物利用的主要碳源。因此,在盐碱地中种植油葵可提高相关土壤理化性质和酶的活性,改善微生物功能多样性,优化盐碱地微生物的群落结构,尤其是种植新葵4号对盐碱地的改良效果最为显著。     

典型淹水稻田土壤微生物群落的基因转录活性及其主要生理代谢过程

【目的】利用环境转录组技术,研究复杂稻田土壤中微生物群落主要生理代谢过程的基因表达水平及其对长期施氮磷钾肥(Mineral nitrogen,phosphorus,and potassium,NPK)的响应规律。【方法】针对中国科学院常熟农田生态系统长期定位试验的NPK施肥处理和不施肥对照处理(Control check,CK)稻田土壤,淹水培养2周后提取土壤微生物总RNA进行高通量转录组测序,利用MG-RAST网络分析平台(Metagenomics Analysis Server)进行活性微生物组成分析、基因功能注释及基因功能分类。【结果】细菌是CK和NPK处理稻田土壤微生物的优势类群,占比高达95%以上,细菌中的活性基因主要源于变形菌门(Proteobacteria,占细菌的50%以上)。同时也检测到古菌、真核生物和病毒等多种微生物的活性基因,而古菌中的活性基因主要源于奇古菌门(Thaumarchaeota,约占古菌的70%)。酸杆菌门(Acidobacteria)在NPK处理土壤中的转录活性显著高于CK处理土壤,而其他的细菌及古菌类群的转录活性在CK和NPK处理土壤间无显著性差异。CK和NPK处理土壤中表达量最高的基因是ABC transporter编码基因,与物质跨膜运输紧密相关。基于COG(Clusters of Orthologous Genes)、Subsystem、KEGG(Kyoto Encyclopedia of Genes and Genomes)3种基因功能分类数据库,发现CK和NPK处理土壤中微生物的主要代谢活动均为能量产生与转化、碳水化合物代谢、蛋白代谢和氨基酸代谢,而最活跃的代谢路径为氧化磷酸化及氨酰-tRNA合成。【结论】淹水状态下CK和NPK处理稻田土壤中的活性微生物组成较为一致,仅Acidobacteria的转录活性在两者间差异较大;在微生物的主要代谢活动方面,CK和NPK处理土壤之间基本一致,均以能量获取与蛋白代谢为主,长期施用无机化肥对复杂土壤微生物群落水平的主要代谢活动影响较小。     

NaCl胁迫下黑果枸杞转录组测序分析

研究黑果枸杞在不同浓度盐胁迫下基因表达谱变化情况,为进一步研究黑果枸杞抗盐分子机制奠定研究基础。对0(CK)、50、250 mmol/L NaCl溶液胁迫的黑果枸杞组培苗的根和叶在胁迫时间为0、1、12 h时分别取样,采用转录组测序(RNASeq)技术进行测序分析。结果表明,转录组测序共产生222.49 Gb原始数据,拼接出Unigenes 86 037条,注释到7大功能数据库(GO、KEGG、KOG、NR、Pfam、Swiss-Prot和egg NOG)上的Unigenes总数为46 594个,占总Unigenes的54.76%,还有38 929个Unigenes在这些数据库中没有得到注释。通过GO分类和KEGG Pathway富集性分析,分别归于51个GO类别和211条代谢途径。差异表达基因分析显示,黑果枸杞叶片和根的上调基因和下调基因数随着NaCl浓度的增大和处理时间的延长均呈增加趋势,叶片中的上调基因数(7 514)小于下调基因数(9 032),根中的上调基因数(12 347)大于下调基因数(11 559)。在黑果枸杞盐胁迫下转录组中发现28 325个SSR位点,最多的为单核苷酸SSR,占70.47%。综合分析表明,黑果枸杞对盐胁迫的反应是一个多基因参与、多个生物过程协同调控的过程,基因表达量的变化可能是基因调控的主要方式。     

阿尔泰蝠蛾(鳞翅目:蝙蝠蛾科)幼虫低温胁迫转录组分析

[目的]本研究旨在通过转录组测序技术分析低温胁迫引起的阿尔泰蝠蛾Hepialus altaicola Wang幼虫基因转录差异及上调表达基因的主要功能类群和参与的主要代谢通路.[方法]采用Illumina HiSeq TM 2500测序平台对阿尔泰蝠蛾幼虫进行转录组测序、组装,利用Blast软件进行数据库比对和基因功能注释,用DEGSeq R软件包分析4℃低温处理与室内适温饲养试虫的差异表达基因,并对上调表达基因进行GO和KEGG代谢途径富集分析.[结果]经序列拼接后共获得100300个unigenes,总长度81600309 bp,平均长度813 bp,N50长度1719 bp.与7大数据库同源比对,共获得34691(34.59％)条unigenes.低温胁迫转录组分析得到11569个差异表达基因(DEGs),7158条基因上调,4411条基因下调.富集到47个GO类群,217个KEGG途径.其中代谢过程、催化、结合活性类群占有重要比例,核糖体、碳代谢、剪接体等途径显著富集.另外,热激蛋白、昆虫表皮蛋白、海藻糖酶、超氧化物歧化酶等非生物胁迫相关基因显著上调表达.[结论]阿尔泰蝠蛾幼虫低温胁迫转录组分析揭示,代谢过程、细胞过程、生物调节、对刺激的反应等生物学过程相关基因和部分非生物胁迫响应基因显著上调表达,提示蝠蛾幼虫可能从抗氧化防御、分子伴侣、体温调节和维持细胞的渗透平衡等多方面应对低温胁迫.     

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

本文以肺形侧耳栽培菌株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分析结果一致。本文通过全面分析肺形侧耳低温胁迫时期基因表达情况,为进一步研究肺形侧耳低温胁迫相关重要基因的功能鉴定与信号通路调控机理提供基础。     

阿尔泰蝠蛾(鳞翅目:蝙蝠蛾科)幼虫低温胁迫转录组分析

[目的]本研究旨在通过转录组测序技术分析低温胁迫引起的阿尔泰蝠蛾Hepialus altaicola Wang幼虫基因转录差异及上调表达基因的主要功能类群和参与的主要代谢通路.[方法]采用Illumina HiSeq TM 2500测序平台对阿尔泰蝠蛾幼虫进行转录组测序、组装,利用Blast软件进行数据库比对和基因功能注释,用DEGSeq R软件包分析4℃低温处理与室内适温饲养试虫的差异表达基因,并对上调表达基因进行GO和KEGG代谢途径富集分析.[结果]经序列拼接后共获得100300个unigenes,总长度81600309 bp,平均长度813 bp,N50长度1719 bp.与7大数据库同源比对,共获得34691(34.59％)条unigenes.低温胁迫转录组分析得到11569个差异表达基因(DEGs),7158条基因上调,4411条基因下调.富集到47个GO类群,217个KEGG途径.其中代谢过程、催化、结合活性类群占有重要比例,核糖体、碳代谢、剪接体等途径显著富集.另外,热激蛋白、昆虫表皮蛋白、海藻糖酶、超氧化物歧化酶等非生物胁迫相关基因显著上调表达.[结论]阿尔泰蝠蛾幼虫低温胁迫转录组分析揭示,代谢过程、细胞过程、生物调节、对刺激的反应等生物学过程相关基因和部分非生物胁迫响应基因显著上调表达,提示蝠蛾幼虫可能从抗氧化防御、分子伴侣、体温调节和维持细胞的渗透平衡等多方面应对低温胁迫.     

高浓度氧或HCl处理条件下家蚕胚胎早期比较转录组分析

本研究以阐明阻断蚕卵滞育机制为目标,通过获得连续时间点的转录组数据,深入解析滞育过程的分子变化.采用高通量测序技术对高浓度氧、HC1处理的滞育蚕卵(以下简称"处理组")和未处理的对照组,在产后26 h、32 h和44h取样进行转录组测序,在产后26 h到32h阶段,处理组和对照组的共同差异基因(DE-Gs)233个,占整个数量的44.5％;32 h到44 h的阶段,处理组和对照组共同DEGs仅16个,占整体数量的6.7％,但是此时处理组的共同DEGs为93个,占整体数量的38.9％.进一步对筛选的DEGs进行GO和KEGG通路分析,GO分析表明DEGs在处理组中在氧化还原功能上得到了显著富集,但是在对照组中没有得到显著富集;KEGG分析发现在产后26 h到32 h,处理组和对照组的DEGs都在溶酶体等代谢通路上调富集,但在产后32 h到44 h,处理组的蚕卵代谢通路变化相似,都在甘氨酸、丝氨酸和苏氨酸等代谢通路富集,而对照组在精氨酸和脯氨酸等代谢通路富集.综上,在产后初期,发育时间对蚕卵发育的影响大于处理的影响,然而到了后期处理的影响则逐渐突显.     

盐生植物海滨锦葵幼苗盐胁迫下基因差异表达分析

利用cDNA-AFLP技术对海滨锦葵幼苗盐胁迫下叶片和根部的基因差异表达模式进行分析和比较, 并对部分盐胁迫应答的转录衍生片段进行了回收、测序和功能推测, 以从转录水平分析海滨锦葵的耐盐分子机制。结果显示：(1) 盐胁迫下海滨锦葵幼苗叶片和根部的基因差异表达多以量的变化为主, 包括盐胁迫下基因表达上调、下调或随盐处理浓度高低和胁迫时间长短而波动的差异表达模式; 只有少量基因的差异表达表现出质的变化, 如盐胁迫下基因沉默或诱导表达; (2) 仅在盐胁迫处理2 h的海滨锦葵幼苗根部, 基因的差异表达以质的变化为主的类型比例略高于量的变化类型比例; (3) 盐胁迫应答基因在不同组织中上调、下调、诱导或沉默的比例随胁迫处理时段而动态变化, 在刚胁迫时基因表达的差异加剧, 而后随胁迫处理时段的延长而渐趋稳定。结果预示, 从基因表达水平探讨植物的耐盐分子机理, 尽管有一定的规律可循, 但由于不同组织对盐胁迫的应答是动态变化的过程, 海滨锦葵不同组织在盐胁迫不同阶段的基因时、空、序表达特征并没有固定的程式。对部分盐胁迫下上调或诱导表达的转录衍生片段(Trivially distributed file system, TDFs)进行的序列分析和功能推测表明, 苗期海滨锦葵在盐胁迫下应答基因至少涉及3类：(1) 离子平衡重建或减少胁迫损伤相关基因(特别是运转蛋白类); (2) 恢复盐胁迫下植物生长和发育相关基因：如参与能量合成和激素调节途径相关基因等; (3)信号转导相关基因及功能未确定的新基因。文章并对盐胁迫应答基因的差异表达模式与海滨锦葵的耐盐性关系进行了讨论。