沙葱萤叶甲成虫不同夏滞育阶段的转录组分析(英文)

【目的】本研究旨在揭示内蒙古草原新害虫沙葱萤叶甲Galeruca daurica专性夏滞育相关的重要基因以及代谢通路。【方法】应用RNA-Seq技术,对沙葱萤叶甲成虫不同夏滞育阶段［滞育前期(PD)、滞育期(D)及滞育后期(TD)］进行转录组测序、分析及基因功能预测,基于RNA-Seq数据筛选夏滞育不同阶段差异表达基因;利用qPCR对基于RNA-Seq数据筛选的10个差异表达基因的表达水平进行验证。【结果】从9个文库中获得202 770 198 clean reads,将12 078 060条转录本组装获得82 292 条unigene,平均长度为783.59 bp,N50为1 545 bp。沙葱萤叶甲D vs PD和TD vs D 比较组分别有2 395(2 119上调和277下调)和62(59上调和3下调)个差异基因。KEGG分析表明,D vs PD和TD vs D比较组差异表达基因分别显著富集于糖孝解/糖异生通路和脂肪酸生物合成通路;此外,许多与钙离子信号转导相关的基因在滞育期间差异表达。10个差异表达基因的qPCR分析表明,RNA-Seq与qPCR结果高度一致。【结论】糖孝解/糖异生、脂肪酸生物合成及钙离子信号通路可能在沙葱萤叶甲滞育调节中起着重要的作用。本研究为进一步研究沙葱萤叶甲成虫专性夏滞育的分子机理奠定了基础。     

沙葱萤叶甲对蜕皮激素响应的转录组分析

为建立沙葱萤叶甲Gauleruca daurica成虫对蜕皮激素（20-hydroxyecdysone,20E）响应的转录组数据库,挖掘对20E响应的基因以及代谢和信号通路,并在转录组水平探讨20E调控生殖滞育的分子机制,本研究采用Illumina HiSeqTM4000高通量测序平台对20E及二甲基亚砜（dimethyl sulfoxide,DMSO）处理后的沙葱萤叶甲成虫进行了转录组测序,共获得80 313个unigene;与阴性对照DMSO相比,共获得201个差异表达基因,其中106个上调、95个下调。GO和KEGG富集分析表明,多数差异表达功能基因富集于各种代谢通路,其中核黄素代谢（riboflavin metabolism）、溶酶体（lysosome）和泛酸与乙酰辅酶A合成（pantothenate and CoA biosynthesis）通路显著富集（q<0.05）。结果表明,20E可能通过影响多种代谢通路调控沙葱萤叶甲的生殖滞育。     

西伯利亚蝗卵发育过程的比较转录组分析及滞育关联基因筛选

【目的】通过筛选西伯利亚蝗Gomphocerus sibiricus卵滞育基因及代谢通路,初步明确西伯利亚蝗卵滞育分子机理。【方法】利用Illumina NovaSeq 6000测序平台对西伯利亚蝗早期发育(ES)、滞育(DS)和滞育后发育阶段(PS)的卵进行转录组测序;采用GO和KEGG富集分析并结合文献,预测滞育相关通路并聚类热图分析,筛选蝗卵滞育关联基因,并选取其中6个重要基因进行qRT-PCR验证。【结果】DSvs ES和PS vs DS两比较组分别富集到12 419和4 789个差异表达基因,且多上调表达;两比较组共表达差异基因为2 206个,主要与糖代谢、环境胁迫和生长发育相关。DS vs ES组富集最显著的GO条目为蛋白结合,PS vs DS组GO条目主要包含酶活性、细胞骨架构建及蛋白结合等过程。滞育关联基因主要集中在Wnt信号通路、胰岛素信号通路、细胞周期通路以及昆虫激素生物合成等通路。6个重要的滞育关联基因的表达量变化趋势与转录组数据结果基本一致。【结论】本研究初步明确了西伯利亚蝗卵滞育调控的重要代谢途径,并筛选出20个滞育关联基因,为后续深入研究该物种的适应机理奠定了基础。     

沙葱萤叶甲保幼激素结合蛋白基因GdJHBP的克隆及表达分析

【目的】克隆沙葱萤叶甲Galerucadaurica保幼激素结合蛋白基因(Juvenilehormonebinding protein,JHBP)cDNA全长序列,分析其分子特征和表达特性,为进一步明确其在沙葱萤叶甲生长发育及滞育中的作用奠定基础。【方法】基于本实验室组装的沙葱萤叶甲转录组数据库,采用RACE技术,克隆沙葱萤叶甲GdJHBP基因cDNA全长序列;运用ORF Finder、SignaIP、DNAMAN和TMHMM等软件分析其分子特征;利用MEGA6.0软件中的邻接法(Neighbor-joining,NJ)构建系统发育树;应用荧光实时定量PCR(RT-qPCR)技术分析GdJHBP在沙葱萤叶甲不同发育时期、成虫不同组织及高温胁迫下的表达模式。【结果】克隆获得了沙葱萤叶甲保幼激素结合蛋白基因GdJHBP cDNA全长序列(GenBank登录号:MG460309),cDNA全长为826 bp,开放阅读框(ORF)为714 bp,编码237个氨基酸;蛋白质预测分子量为26.58ku,等电点为4.37;包含1条信号肽,无跨膜区,且在第27-189位氨基酸之间存在一个保幼激素结合蛋白家族JHBP保守结构域。序列比对分析表明,不同昆虫JHBP间氨基酸序列一致性较低,沙葱萤叶甲GdJHBP与棕榈象RhynchophorusferrugineusRfJHBP和马铃薯甲虫Leptinotarsa decemlineata JHBP 3p2的氨基酸序列一致性最高也仅为30%。系统发育分析表明,GdJHBP与棕榈象血淋巴JHBP亲缘关系最近。RT-qPCR结果显示,GdJHBP在沙葱萤叶甲不同发育阶段均有表达,在幼虫期表达量最高,在卵和蛹期微量表达;在成虫滞育期间低表达,在滞育前与滞育结束后则有较高表达;在成虫发育过程中,头部的表达量显著低于腹部和胸部;高温(30-40℃)可诱导GdJHBP上调表达,在35℃时表达量达到最高值。【结论】沙葱萤叶甲GdJHBP属于血淋巴JHBP,在沙葱萤叶甲生长发育和成虫夏滞育中可能发挥着重要作用。     

沙葱萤叶甲核糖体蛋白基因GdRpS3a的克隆、分子特性及表达分析

核糖体蛋白(ribosomal protein,Rp)是一类参与蛋白质合成、细胞代谢、机体免疫及信号传导等重要功能的蛋白质.本研究旨在克隆沙葱萤叶甲Galeruca daurica核糖体蛋白S3a基因cDNA全长序列,分析其分子特性和表达模式,以期为进一步研究其在沙葱萤叶甲生长发育及滞育中的作用提供必要的基础.根据现有的沙葱萤叶甲转录组测序数据,应用cDNA末端快速扩增(RACE)技术克隆获得了沙葱萤叶甲RpS3a基因的全长cDNA序列,命名为GdRpS3a(GenBank登录号:MN660144).该基因全长为800 bp,开放阅读框全长为687 bp,编码228个氨基酸,具有核糖体蛋白S3a家族的保守功能域;预测分子量为25.63 kDa,等电点pI为10.53,无信号肽和跨膜结构.同源序列比对和系统发育分析表明,GdRpS3a与玉米根萤叶甲Diabrotica virgifera virgifera RpS3a的亲缘关系最近,其氨基酸序列一致性最高为54.10％.采用qPCR技术检测该基因在沙葱萤叶甲不同发育阶段(卵、1～3龄幼虫、预蛹、蛹及成虫羽化3、7、10、15、25、40、60、80和100 d)和3日龄成虫在不同温度(0、5、10、15、20、25、30、35和40℃)处理1h后的表达谱.结果 表明:GdRpS3a在不同发育阶段的沙葱萤叶甲体内均有表达,在成虫滞育结束后表达量最高,在卵期和成虫滞育期间低表达.温度对GdRpS3a的表达有显著的影响,30℃下最高,0℃下最低.GdRps3a的表达与沙葱萤叶甲生长发育阶段及环境温度有关,可能在沙葱萤叶甲成虫夏滞育中起着重要作用.     

沙葱萤叶甲钙结合蛋白基因的鉴定及表达谱分析

【目的】沙葱萤叶甲Galeruca daurica是一种近年来在内蒙古草原上猖獗成灾的新害虫。本研究旨在克隆沙葱萤叶甲钙结合蛋白(calcium-binding protein, CaBP)基因,分析其在沙葱萤叶甲成虫不同发育阶段及不同温度下的表达谱,为进一步探究其在沙葱萤叶甲生长发育及滞育过程中的作用奠定基础。【方法】根据沙葱萤叶甲转录组和蛋白质组数据,筛选CaBP基因序列信息,应用RT-PCR技术克隆获得CaBP基因的开放阅读框(ORF)全长序列,并对其进行生物信息学分析;通过qPCR检测其在沙葱萤叶甲成虫不同日龄(羽化后3, 7, 10, 15, 20, 30, 40, 60, 90和110 d)及3日龄成虫在不同温度(0, 5, 10,15, 20, 25, 30和35℃)下处理1 h后的表达水平。【结果】克隆得到4条具有完整ORF的沙葱萤叶甲CaBP基因cDNA序列,分别命名为GdCaM, GdCAPSL, GdTnCl和GdCRT(GenBank登录号: MN695412-695415),ORF全长分别为480, 648, 516和1 209 bp,分别编码149, 215, 171和402个氨基酸;只有GdCRT拥有信号肽。同源序列比对和系统发育分析表明,GdCaM, GdCAPSL, GdTnCl和GdCRT分别与玉米根萤叶甲Diabroticavirgifera virgifera的CaM, CAPSL, TnCl及马铃薯甲虫Leptinotarsa decemlineata CRT的氨基酸序列一致性最高,分别为100.0%, 74.0%, 88.2%和92.5%。qPCR结果表明,4个CaBP基因在沙葱萤叶甲不同日龄成虫中均差异表达,且表达模式不同。GdCaM在成虫滞育前(羽化后3 d)表达量较高,进入滞育后(羽化后7 d)表达量降低,在滞育期间(羽化后7-60 d)表达量变化较小,而解除滞育后(羽化后90 d)表达量进一步下调。GdCAPSL在滞育初期(羽化后10 d)维持在最低水平,进入滞育中后期(羽化后40和60 d)开始回升,滞育解除后又突然下调至最低水平,而羽化后110 d急剧上升至最高水平。GdTnCl在滞育初期(羽化后15-20 d)高表达,进入滞育中后期(羽化后30-60 d)急剧下降至最低水平,滞育解除后再次上调,但羽化后110 d突然下调至最低水平。GdCRT在进入滞育后表达量开始逐渐下调,在滞育维持期间(羽化后15-60 d)维持在低水平,滞育解除后又开始上升。温度对3日龄成虫中除GdCRT外的其他3个CaBP基因的表达有显著影响。温度低于20℃时,GdCaM的表达量随温度的降低而升高,但0℃时突然下降;温度高于20℃时,GdCaM的表达量随着温度的升高而上升。GdCAPSL的表达量随着温度的升高而呈现上升的趋势,25℃时达到最高,然后下降。GdTnCl随着温度的升高,表达量呈现下降的趋势。【结论】钙结合蛋白可能在沙葱萤叶甲成虫生长发育及夏滞育调控过程中发挥着重要作用。     

伞裙追寄蝇滞育关联基因的转录组学分析

【目的】本研究旨在探讨调控伞裙追寄蝇Exorista civilis滞育的重要关联基因以及代谢途径,为明确该虫在转录组水平下的滞育分子机制提供理论基础。【方法】采用新一代高通量测序平台Illumina HiSeq^TM2000对伞裙追寄蝇非滞育以及滞育的蛹进行转录组测序分析以及生物信息学分析;应用KAAS在线pathway比对分析工具对筛选出的满足padj＜0.05且fold change≥32或padj＜0.05且fold change≤1/32条件的差异表达基因进行KEGG通路富集分析。【结果】根据测序结果,共获取58 050个基因。差异倍数在32倍以上的滞育关联基因(diapause-associated genes, DAGs)有454个,其中406个表达上调,共涉及134条通路,包括氧化磷酸化、柠檬酸循环及其他重要通路;48个表达下调,涉及32条通路。KEGG通路富集分析结果表明,滞育关联基因主要集中在信号转导、内分泌系统、碳水化合物代谢等途径中。【结论】本研究获得的伞裙追寄蝇转录组数据揭示了其滞育调控的重要代谢途径与关联基因。     

沙葱萤叶甲海藻糖酶基因GdTre1的克隆、分子特性和表达分析

【目的】海藻糖酶(trehalase,Tre)作为昆虫体内海藻糖代谢的关键性酶,在昆虫能量调节和生长发育中具有重要作用。本研究旨在克隆获得沙葱萤叶甲Galeruca daurica海藻糖酶基因,并对其表达模式进行定量分析,以期探究海藻糖酶在沙葱萤叶甲生长发育和成虫夏滞育中的作用。【方法】根据沙葱萤叶甲转录组数据,采用RACE技术,克隆得到Tre基因的c DNA全长,并进行生物信息学分析;应用实时荧光定量PCR(RT-q PCR)技术检测其在沙葱萤叶甲不同发育阶段[卵、1-3龄幼虫、预蛹、蛹、成虫(羽化后3,7,10,15,25,40,60,80和100 d)]、成虫不同组织(头、胸和腹)以及3日龄成虫在不同温度(15,20,25,30,35和40℃)胁迫下的表达水平;采用3,5-二硝基水杨酸法测定不同日龄成虫体内海藻糖酶活性。【结果】克隆获得了沙葱萤叶甲可溶性海藻糖酶基因,并命名为Gd Tre1(Gen Bank登录号:KY697913),该基因全长1 933 bp,开放阅读框1 704bp,编码567个氨基酸;蛋白预测分子量为66.56 k D,等电点为6.62;编码蛋白具有海藻糖酶超基因家族典型的功能结构域,包含1条信号肽,不具有跨膜结构。同源序列比对和系统发育分析表明,Gd Tre1与马铃薯甲虫Leptinotarsa decemlineata Tre1b的同源性最高,氨基酸序列一致性达70.25%。RT-q PCR检测结果表明,Gd Tre1在沙葱萤叶甲各发育阶段均有表达,其中在卵期和成虫滞育期间高表达,而在幼虫、预蛹、蛹及成虫滞育前低表达;在成虫不同组织中,通常在腹部中表达量最高,其次为胸部,最低为头部;Gd Tre1表达量随着温度升高而上升,30℃达最高值,而后随温度升高略有下降。沙葱萤叶甲不同日龄成虫体内Gd Tre1表达量及Tre活性存在显著差异,且Gd Tre1表达量与Tre活性变化趋势一致。【结论】海藻糖酶与沙葱萤叶甲生长发育和成虫夏滞育有密切的关系。该结果为进一步揭示该虫的夏滞育分子机理提供了必要的基础。     

沙葱萤叶甲化学感受蛋白基因的鉴定及表达谱分析（英文）

【目的】沙葱萤叶甲Galeruca daurica是一种在内蒙古草原上爆发成灾的新害虫。昆虫的化学感受蛋白(CSPs)是一类水溶性小分子蛋白,其主要功能是识别和传导环境中的化学刺激至受体,参与众多昆虫行为。本研究旨在鉴定沙葱萤叶甲化学感受蛋白基因,并对其表达谱进行分析。【方法】通过筛选本实验室组装的沙葱萤叶甲转录组鉴定沙葱萤叶甲化学感受蛋白基因;采用实时荧光定量PCR方法分析其在沙葱萤叶甲不同发育阶段(卵、1-3龄幼虫、蛹和成虫)和成虫组织[触角、头(去触角)、胸、腹、足和翅]中的表达水平。【结果】鉴定出10个化学感受蛋白基因,将其命名为Gdau CSP1-10(Gen Bank登录号:KY885471-KY885480)。10条编码蛋白的氨基酸序列一致性范围为17.27%~62.79%,彼此间分化程度较高。通过NCBI的Blast比对结果显示,Gdau CSPs与榆黄毛萤叶甲Pyrrhalta maculicollis的Pmac CSP氨基酸序列一致性最高,为90%。系统发育分析表明,4种Gdau CSPs首先与榆黄毛萤叶甲的Pmac CSPs聚为一支。实时荧光定量PCR结果显示,Gd CSPs在不同发育阶段中的表达量有显著差异,Gdau CSP4-5,Gdau CSP7-8和Gdau CSP10 5个Gd CSPs在成虫中的表达量显著高于其他发育阶段,而Gd CSP2在卵中的表达量显著高于其他发育阶段;Gd CSPs不仅表达于成虫触角中,在头(去触角)、胸、腹、足、翅等其他部位也有表达,且10个Gd CSPs在沙葱萤叶甲成虫各组织中有着不同的表达谱,其中Gdau CSP2,Gdau CSP4,Gdau CSP5,Gdau CSP8和Gdau CSP9 5个Gd CSPs在雌成虫触角中的表达量显著高于其他组织。【结论】结果提示化学感受蛋白在沙葱萤叶甲的生长发育和化学感受过程中可能起着不同的作用。本研究为进一步研究沙葱萤叶甲化学感受蛋白的生理功能及化学通讯的分子机理奠定了必要的基础。     

沙葱萤叶甲丝氨酸蛋白酶基因GdSP的克隆及对温度胁迫的响应

【目的】丝氨酸蛋白酶(serine protease,SP)是一类以丝氨酸为活性中心的重要的蛋白水解酶。本研究旨在克隆获得沙葱萤叶甲Galeruca daurica丝氨酸蛋白酶基因,分析其对温度胁迫的响应,以期为进一步揭示沙葱萤叶甲耐温性的调控机制及其他生理功能奠定基础。【方法】根据沙葱萤叶甲2龄幼虫转录组数据,采用RACE技术克隆得到沙葱萤叶甲丝氨酸蛋白酶基因的cDNA全长序列,并进行生物信息学分析;应用qPCR技术检测其在不同温度(-10,-5,0,5,25和35℃)下处理1 h后及25℃下恢复30 min后在沙葱萤叶甲2龄幼虫中的表达量变化。【结果】自沙葱萤叶甲克隆获得一个丝氨酸蛋白酶基因,命名为GdSP(GenBank登录号:MG797556)。该基因全长1 110 bp,开放阅读框969 bp,编码322个氨基酸;蛋白预测分子量35.41 kD,等电点5.61;编码蛋白具有丝氨酸蛋白酶的典型特征,具有一个跨膜结构,无信号肽。同源序列比对和系统发育分析表明,GdSP与光肩星天牛Anoplophora glabripennis SP的同源性最高,氨基酸序列一致性为30.53%。qPCR测定结果表明,不同温度处理间2龄幼虫中GdSP表达量差异不显著,但对各高低温(-10℃除外)胁迫处理回温后GdSP表达量显著上升。【结论】快速冷驯化对沙葱萤叶甲丝氨酸蛋白酶基因表达无显著影响,而回温可诱导其上调表达。     

七星瓢虫滞育关联基因的转录组学分析

对七星瓢虫正常发育、滞育以及滞育解除的雌成虫进行RNA测序,并对筛选出来的滞育相关基因进行KEGG通路富集分析,从分子水平解析七星瓢虫滞育发机理。本研究以正常发育产卵、滞育30 d以及滞育贮存30 d后解除产卵的七星瓢虫雌成虫为研究对象,分别抽提RNA,合成c DNA,构建c DNA文库,文库检测合格后在Illummina Hiseq 2500测序仪上进行双向测序。根据测序结果,共获取unigene 82820个。采用两两比较法对正常发育组和滞育组、滞育组和滞育解除组进行差异表达分析,分别获得差异表达基因3501个和1427个。深入分析两组比对结果,将在滞育组上调且滞育解除组下调的unigene定义为滞育关联基因,共有443个基因为滞育关联基因。应用KEGG KAAS在线pathway比对分析工具对滞育关联基因进行通路富集分析,结果发现这些基因主要集中在碳水化合物代谢、脂质代谢以及信号转导等途径中。     

感染球孢白僵菌的家蚕免疫应答差异表达基因分析

【目的】筛选家蚕Bombyx mori应对白僵菌Beauveria bassiana侵染的应答基因, 以进一步研究家蚕抵御真菌侵染的分子机制。【方法】采用新一代Solexa高通量测序技术对感染白僵菌及未感染白僵菌的对照组家蚕进行了测序分析, 筛选差异表达基因; 结合生物信息学工具分析差异表达基因的功能注释、 分类及涉及的信号通路等; 应用荧光定量PCR技术验证10个基因的差异表达。【结果】通过测序和生物信息学分析共获得377个差异表达基因, 其中表达上调基因236个, 下调基因141个; KEGG通路分析表明, 各通路中既有表达上调的基因, 也有下调基因; 12个上调基因、 26个下调基因参与3个显著性富集的KEGG通路, 即核糖体、 氨酰tRNA生物合成和剪接体通路。定量PCR与测序结果显示, 溶菌酶、 热激蛋白、 谷胱甘肽S-转移酶、 肽聚糖识别蛋白等与免疫应激相关的蛋白基因均呈现表达上调。【结论】本研究筛选获得的差异表达基因, 特别是上调表达的基因可能与家蚕应对白僵菌侵染的应答机制有关, 其中与免疫应激相关的蛋白基因如溶菌酶、 热激蛋白、 谷胱甘肽S 转移酶、 肽聚糖识别蛋白基因等可能直接参与了家蚕对白僵菌的免疫识别和防御, 研究结果为从分子水平阐明家蚕抵御真菌侵染的防御机制和白僵菌对家蚕的致病机理提供新的依据。     

家蚕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及潜在靶基因,为探究家蚕精巢和卵巢发育差异奠定了基础。     

中华蜜蜂精巢和卵巢差异表达基因分析

[目的]在蜂群中,雄蜂与蜂王都有着发育完全的性腺,但两者达到性成熟的时间却是不同步的.本研究旨在探究中华蜜蜂Apis cerana cerana雄蜂与蜂王生殖腺的基因表达差异.[方法]利用Illumina测序技术对中华蜜蜂雄蜂精巢与蜂王卵巢转录组差异表达基因(differentially expressed genes...     

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

【目的】蜜蜂是典型的具有发育狭温性的全变态昆虫。本研究以对低温最敏感的意大利蜜蜂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可能是低温胁迫抑制蜜蜂化蛹的一个关键基因。本研究为探索低温胁迫影响蜜蜂发育变态的分子调控机制奠定了基础。     

草地贪夜蛾雄性成虫和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中涉及解毒相关的通路显著富集。【结论】这些转录组数据为探究草地贪夜蛾的生长发育、嗅觉相关功能基因以及候选分子靶标提供资源信息。