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

【目的】本研究旨在探讨调控伞裙追寄蝇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通路富集分析结果表明,滞育关联基因主要集中在信号转导、内分泌系统、碳水化合物代谢等途径中。【结论】本研究获得的伞裙追寄蝇转录组数据揭示了其滞育调控的重要代谢途径与关联基因。     

伞裙追寄蝇滞育关联基因和滞育关联代谢物的联合分析

为探索伞裙追寄蝇蛹滞育调控的分子机制,本文对伞裙追寄蝇Exorista civilis滞育蛹和非滞育蛹进行转录组测序以及代谢组检测,鉴定其关键的滞育关联基因(diapause-associated genes, DEGs)和滞育关联代谢物(diapasuse-associated metabolites, DEMs)。本研究基于高通量测序以及液质联用技术,通过筛选,在转录组中获得了差异表达基因7 513个,在代谢组中,获得差异代谢物501个,其中氨基酸占比最多。将所有差异表达基因与差异表达代谢物同时向KEGG映射,获得两者共同的pathway信息,明确差异表达基因与差异代谢物共同参与的主要生化途径和信号转导途径。本研究通过转录组和代谢组的联合分析,在正、负离子模式下,差异表达基因与差异代谢物共同富集到70条通路。在负离子模式下,滞育关联因子主要参与氨基酸代谢和神经系统;在正离子模式下,滞育关联因子主要参与消化系统和信号转导途径。本研究重点分析了柠檬酸循环、cAMP信号通路、氨酰-tRNA的生物合成,将为进一步深入研究伞裙追寄蝇滞育调控的分子机制奠定理论基础。     

沙葱萤叶甲卵滞育的转录组学分析

【目的】建立沙葱萤叶甲Galeruca daurica滞育卵转录组数据库,挖掘卵滞育相关的基因以及代谢和信号通路,在转录组水平探讨卵滞育的分子机制。【方法】采用Illumina NovaSeq6000高通量测序平台对沙葱萤叶甲滞育卵与解除滞育卵进行转录组测序,并进行生物信息学分析;利用DESeq软件分析沙葱萤叶甲滞育卵与解除滞育卵中的差异表达基因,对差异表达基因进行KEGG通路富集分析;利用qRT PCR技术对10个差异表达基因的表达模式进行验证。【结果】基于沙葱萤叶甲滞育卵与解除滞育卵转录组测序结果,共获得53 389个unigene,其中差异表达基因2 145个,24个差异表达基因与保幼激素信号及脂肪酸生物合成和降解相关。与解除滞育卵相比,滞育卵转录组中1 297个基因上调表达,富集于124条KEGG通路,其中核糖体通路显著富集;848个基因下调表达,富集于73条KEGG通路,其中MAPK信号通路和糖胺聚糖生物合成通路显著富集。qRT-PCR结果表明,随机选取的10个差异表达基因的表达趋势与RNA-Seq转录组测序结果完全一致。【结论】保幼激素,脂肪酸生物合成和降解,核糖体,MAPK信号及糖胺聚糖生物合成等通路可能在沙葱萤叶甲卵滞育调控中起着重要的作用。     

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

【目的】通过筛选西伯利亚蝗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个滞育关联基因,为后续深入研究该物种的适应机理奠定了基础。     

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

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

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

【目的】本研究旨在揭示内蒙古草原新害虫沙葱萤叶甲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可能通过影响多种代谢通路调控沙葱萤叶甲的生殖滞育。     

RNA-Seq揭示Foc4在外源氧化胁迫(H_2O_2)下的基因表达及细胞代谢变化

尖孢镰刀菌古巴专化型4号小种(Fusarium oxysporum f. sp. cubense race 4,Foc4)是香蕉枯萎病的强致病性病原菌。Foc4在侵染香蕉植株早期必须面对寄主的活性氧迸发。【目的】了解Foc4应对外源氧化胁迫的分子机制。【方法】利用Illumina 2500 RNA-Seq测序平台分析了经外源氧化胁迫(H_2O_2)处理的Foc4与对照在转录组水平的基因表达差异。【结果】在外源氧化胁迫条件下,Foc4的生长受到抑制。转录组测序获得了超过2千万条clean reads。进一步的差异基因表达分析以差异倍数FC (fold change)≥2且FDA值≤0.001为选择标准,发现496个基因表达上调,298个基因表达下调。GO功能富集分析显示,429个基因比对到GO功能分析数据库,在这些差异表达基因中,许多与代谢过程、生物调节、细胞过程和刺激应答有关。KEGG通路富集分析显示,有141个表达差异显著基因比对到KEGG中的50条代谢途径。其中,主要是各类氨基酸代谢途径、脂肪酸代谢途径。同时也包括与抗氧化胁迫直接相关的代谢途径,包括DNA的损伤修复、类胡萝卜素的生物合成、过氧化物酶体、谷胱甘肽代谢等。【结论】这些结果暗示,为了在强氧化胁迫环境下生存,Foc4细胞从包括直接应对氧化胁迫的信号调控途径在内的物质代谢和能量代谢均发生改变以应对环境变化的胁迫。     

茶足柄瘤蚜茧蜂滞育和非滞育蛹中与能量代谢相关的差异表达蛋白

【目的】本研究旨在从蛋白质组整体层面阐明茶足柄瘤蚜茧蜂Lysiphlebus testaceipes蛹滞育背后的多蛋白调控,重点筛选与能量代谢相关的滞育关联蛋白并分析其功能,有助于更好地理解茶足柄瘤蚜茧蜂蛹滞育的代谢机制。【方法】利用同位素标记相对和绝对定量(isobaric tags for relative and absolute quantification, iTRAQ)技术比较了茶足柄瘤蚜茧蜂滞育蛹与非滞育蛹的蛋白含量;利用GO, KEGG网络数据库等生物信息学方法分析鉴定茶足柄瘤蚜茧蜂滞育蛹与非滞育蛹中差异表达蛋白(differentially expressed proteins, DEPs)。【结果】分析得到茶足柄瘤蚜茧蜂滞蛹与非滞育蛹DEPs有135个,包括滞育蛹中上调表达蛋白有38个,下调表达蛋白有97个。GO和KEGG富集分析表明,与天冬氨酸转运、L-谷氨酸转运、胆碱脱氢酶活性和胆碱生物合成甘氨酸甜菜碱条目以及氧化磷酸化通路相关的蛋白上调表达。【结论】氧化磷酸化通路相关蛋白在茶足柄瘤蚜茧蜂滞育过程中呈显著上调表达,说明能量代谢与该蜂滞育密切相关,并推测氧化磷酸化...     

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

茶足柄瘤蚜茧蜂蛹滞育过程中胰岛素信号通路及其相关途径的初探

本实验通过探索胰岛素信号通路及其相关途径对茶足柄瘤蚜茧蜂蛹滞育的影响,从而方便寻找胰岛素替代物,为害虫防治提供新思路。利用RNA-Seq,对滞育组与非滞育组的茶足柄瘤蚜茧蜂进行转录组测序,结合生物信息学方法对转录组中胰岛素信号通路及其相关途径的差异表达基因进行了分析。与胰岛素信号通路相关差异表达基因共31个,重点分析的PI3K-Akt, FoxO, MAPK三条途径,差异表达基因分别为55, 21和28个。这些滞育关联基因呈现不同程度的上调或下调表达,发现Sos, FASN, TSC1, PRKAB等基因与茶足柄瘤蚜茧蜂滞育密切相关,共同影响茶足柄瘤蚜茧蜂的滞育。胰岛素信号通路及其相关途径对茶足柄瘤蚜茧蜂的滞育起着非常重要的作用,主要体现在影响虫体能量代谢、脂质积累、细胞增殖等方面。     

茶足柄瘤蚜茧蜂碳水化合物代谢相关的滞育关联基因差异表达分析

为明确糖代谢相关途径在茶足柄瘤蚜茧蜂Lysiphlebustestaceipes蛹滞育过程中的作用,揭示滞育调控的分子机制,本试验利用转录组测序技术,对滞育组与非滞育组的茶足柄瘤蚜茧蜂蛹进行转录组测序,并结合生物信息学方法对糖代谢相关途径中的差异表达基因进行了筛选与分析.GO注释到的与碳水化合物代谢条目相关的差异基因共...     

苗期水稻响应褐飞虱取食的基因差异表达分析

【目的】褐飞虱Nilaparvata lugens (St?l)是水稻的重要害虫之一,主要刺吸水稻韧皮部汁液为害,对水稻生产造成严重的产量和经济损失。为研究水稻响应褐飞虱取食的分子机制,对褐飞虱取食6 h后的苗期水稻进行转录组测序及基因差异表达分析。【方法】采用illumina二代测序技术获得褐飞虱取食前后水稻组织的转录组数据,利用RSEM软件进行基因表达定量和DEseq2进行差异表达分析;从差异表达基因中随机选取20个基因采用荧光定量PCR技术进行验证;采用GeneMerge软件对差异表达基因进行KEGG和GO富集分析。【结果】褐飞虱取食后,水稻转录组中的1 104个基因出现了差异表达,其中435个基因表达上调,669个基因表达下调。荧光定量PCR结果显示,20个差异表达基因中18个基因的表达变化趋势和测序结果一致,证明了转录组分析结果可靠。GO和KEGG富集分析表明,表达上调基因主要与水稻氧化应激、海藻糖合成及次生化合物代谢有关,显著富集在14个KEGG通路和30个GO功能分类中;而表达下调基因主要参与水稻纤维素、蛋白质及脂肪酸合成过程,显著富集在29个KEGG通路和26个GO功能分类。在差异表达基因中,分别有61个转录因子和13个水杨酸和茉莉酸信号通路相关基因。【结论】褐飞虱取食激发了水稻的应激反应和保护机制,同时还降低了营养合成的过程,是飞虱为害造成水稻减产的原因之一。本研究初步揭示了苗期水稻响应褐飞虱取食的差异表达基因,为研究水稻-褐飞虱互作机制以及褐飞虱抗性水稻品种培育提供了参考和依据。     

草地贪夜蛾雄性成虫和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中涉及解毒相关的通路显著...     

准噶尔小胸鳖甲短时低温胁迫响应的转录组分析

【目的】拟步甲科昆虫小胸鳖甲Microdera punctipennis是分布于中国新疆古尔班通古特沙漠的特有物种,具有很强的耐寒性,但其低温响应的分子机制尚不明确。本研究旨在利用转录组测序技术丰富小胸鳖甲的已知基因信息,分析在短时低温胁迫下,上调表达基因的主要类群及参与的主要代谢通路。【方法】利用Trinity软件对分别在4℃低温和25℃(对照)下处理3 h的小胸鳖甲成虫RNA-Seq数据进行从头组装,利用Blast2go软件进行基因注释。通过DESeq和GFOLD软件分析差异表达基因,并对上调表达基因进行GO(Gene Ontology)和KEGG(Kyoto Encyclopedia of Genes and Genomes)代谢途径富集分析。【结果】测序过滤后得到68 165 001个有效读序,51 712个平均长度为984 bp的非冗余基因,其中29 925个基因(约57.87%)有同源序列(E-value10-5),与拟步甲科的赤拟谷盗Tribolium castaneum相应基因相似性最高(核苷酸序列一致性为72.75%)。低温响应上调表达基因有514个,富集到35个GO类群,18个KEGG途径。其中胁迫响应类群、生物调节类群和免疫系统过程类群都占有重要比例,嘌呤代谢、硫胺素代谢、糖酵解/糖异生等代谢通路显著富集。在胁迫响应类群中,热激蛋白Hsp90基因、超氧化物歧化酶(SOD)基因和钙联接蛋白Calnexin基因等8个非生物胁迫相关基因显著上调表达。【结论】荒漠昆虫小胸鳖甲低温转录组揭示,在4℃冷驯化过程中,代谢过程、胁迫响应、生物调节和免疫系统等生物学过程相关基因表达显著上调,其中几个非生物胁迫响应基因提示小胸鳖甲可能从分子伴侣、细胞周期阻滞、线粒体稳定性、抗氧化胁迫等多方面应对低温胁迫。KEGG富集分析所揭示的小胸鳖甲在低温下的转录组代谢通路与其他物种有较大差异。研究结果为后续生物信息学分析及小胸鳖甲耐寒关键基因的发掘及耐寒机制的揭示提供了基础数据。     

球孢白僵菌诱导的柑橘木虱免疫应答转录组分析

【目的】筛选柑橘木虱Diaphorina citri Kuwayama应对球孢白僵菌侵染的免疫应答及其网络调控基因,以进一步探讨柑橘木虱对球孢白僵菌的免疫防御机制。【方法】采用Illumina高通量测序平台对感染球孢白僵菌24、48、72 h与健康的柑橘木虱转录组进行了测序,利用生物信息学软件对筛选到的差异表达基因进行了功能注释、分类以及参与的信号通路分析。【结果】组装得到138 313条不可延长的非冗余Unigene,其N50和N90分别为2 532 bp和413 bp,平均长度为1 191.26 bp。在CK vs.S24h、CK vs.S48h和CK vs.S72h 3个转录组里分别获得了971、1 671和752个显著差异表达基因(DEGs),GO富集分析表明分别有405、614和542个DEGs显著富集到56、83和60个GO terms中,KEGG pathway分析结果显示分别有98、333和247个DEGs显著富集到10、27和25条代谢通路里。【结论】差异表达基因主要富集在能量代谢、离子运输、转录和翻译调控、生殖和发育调控以及免疫防御反应等相关通路,大部分编码潜在的与免疫识别及调控相关的基因,并从中筛选到5个显著上调表达的免疫相关基因,为开展柑橘木虱免疫应答虫生真菌的研究奠定理论基础。下一步将进行免疫相关基因的q-PCR验证及表达量分析,研究其在柑橘木虱应对球孢白僵菌侵染过程中的作用。     

感染球孢白僵菌的烟粉虱若虫免疫应答转录组分析

【目的】筛选烟粉虱Bemisia tabaci(Gennadius)若虫应对白僵菌Beauveria bassiana侵染的应答基因,以进一步研究烟粉虱免疫反应的分子机制。【方法】采用新一代高通量测序技术对感染和非感染白僵菌的烟粉虱4龄若虫进行了测序分析,并筛选了差异表达基因;利用生物信息学工具对转录组测序得到的基因进行了功能注释、分类以及参与的信号通路展示。【结果】组装得到非冗余Unigene 232 554个,其N50和N90分别为1 153 bp和260 bp,平均长度为67 424 bp。对所有的基因进行差异表达分析发现:以P<0.05为标准筛选得到了1 166个差异表达基因,其中上调表达基因有474个,下调表达基因有692个,其中,与免疫相关的基因有405个;GO富集分析发现:有416个GO term有富集现象,包括156个生物学过程(66 402个Unigenes),89个细胞组分(27 645个Unigenes)和154个分子功能(73 417个Unigenes);KEGG代谢通路分析将烟粉虱转录组1 145个DEG匹配到309个通路上,其中76个通路得到了富集。【结论】对405个可能参与到烟粉虱若虫对白僵菌侵染的免疫识别和防御的基因进行了测序。研究结果为从分子水平上开展应用虫生真菌防治烟粉虱的研究奠定信息学基础。     

NaCl胁迫下哈茨木霉ACCC32524的转录组和代谢组分析

【目的】通过分析NaCl胁迫下哈茨木霉(Trichoderma harzianum)ACCC32524转录组和代谢组数据,研究差异表达基因及次级代谢产物的变化情况,初步探索响应NaCl胁迫的分子机制。【方法】利用Illumina HiSeq XTen高通量测序平台完成0、0.4、0.6 mol/L NaCl浓度胁迫培养下哈茨木霉ACCC32524的转录组测序,GC-TOF-MS技术完成对0mol/L和0.6mol/LNaCl胁迫培养下的差异次级代谢产物检测,利用相关软件及数据库对差异表达基因(DEGs)和次级代谢产物的注释、筛选和分类,并进行RT-qPCR验证。【结果】本研究分别得到0.4 mol/L和0.6 mol/L NaCl胁迫下417和733条差异表达基因;GO富集分析显示,分别有318和582条差异表达基因注释到生物学过程、分子功能和细胞组分3个一级分类和40个二级分类;COG分类结果表明分别有232和414条转录本为20个类别,涉及差异表达基因最多的分别为氨基酸的转运和代谢、一般功能预测、碳水化合物的转运和代谢;KEGG代谢途径分析结果表明,分别有75和96条基因归到25个代谢通路中(P≤0.05),其中涉及差异基因最多的是氨基酸的生物合成和2-氧代羧酸代谢通路。从转录组数据中共筛选出与渗透调节、离子转运、活性氧清除等22个耐盐相关基因。0 mol/L和0.6 mol/L NaCl胁迫下的代谢组数据中共筛选出101个差异次级代谢产物,包括8种积累量上调和93种下调物质,其中36个得到定性,分属于糖类、有机酸和氨基酸等9个分类中。RT-qPCR验证挑选的差异表达基因的表达量变化,均与RNA-seq分析结果一致。【结论】NaCl胁迫下引起哈茨木霉ACCC32524基因及次级代谢产物发生明显变化,细胞代谢途径发生明显偏移,这些进程共同作用减少NaCl对细胞的毒害作用,为木霉菌的耐盐机理研究提供重要信息。     

miRNA响应印度梨形孢定殖调控大麦生长发育的作用机制

【背景】内生真菌印度梨形孢(Piriformospora indica)定殖植物可以显著促进植物生长发育。miRNA已被证实在植物体的生长发育中具有调控作用。【目的】揭示印度梨形孢定殖大麦促进大麦生长发育过程中miRNA对印度梨形孢定殖的响应及对大麦生长发育的调控作用。【方法】提取大麦总RNA,实施转录组测序并进行序列比对与数据挖掘;使用高效液相色谱检测大麦生长素等激素水平变化。【结果】印度梨形孢对大麦有显著促生作用;全转录组测序结果显示：印度梨形孢侵染3 d较空白对照有18个差异表达的miRNA,其中11个miRNA上调、7个miRNA下调;侵染7 d与空白对照相比24个差异表达的miRNA,其中11个miRNA上调、13个miRNA下调;侵染3 d与侵染7 d相比有3个miRNA上调、6个miRNA下调。GO功能富集分析与KEGG通路分析显示,差异表达miRNA的靶基因主要参与转录、细胞分裂、生长素信号的感知和转导、光合作用和激素刺激响应。靶基因所参与的途径与大麦生长发育密切相关,暗示miRNA对印度梨形孢定殖过程做出了积极响应。代谢产物分析表明miRNA参与的调控路径的代谢产物发生改变。【结论】本研究以miRNA为入手点,探究了miRNA对大麦生长发育的调控机制,为揭示印度梨形孢的促生机制提供了新的研究方向。     

叶面喷水影响芒果叶片午休的转录组分析

【目的】通过组学技术对持续喷水处理的芒果叶片产生的影响进行研究,揭示喷水处理打破“午休”的分子机制,为进一步优化喷水处理条件以实现最佳效果提供科学依据,并对提高芒果产量和品质以及响应国家减肥增效政策具有理论意义。【方法】以元江“台农一号”芒果品种为试验材料,对坐果期芒果在“午休”时段(12:30—14:10)喷水3次,每次喷水20 min,间隔20 min,以不喷水为对照,并以2组3个不同时期叶片进行Illumina Hi SeqTM 62000高通量转录组测序分析,利用基因本体(GO)和京都基因与基因组百科全书(KEGG)数据库分析“午休”时段喷水与不喷水处理的差异表达基因(DEGs)。采用DESeq对基因表达进行差异分析,以|log2αFC|>1,P<0.05为筛选差异表达基因的条件。【结果】3个时期分别获得3 789,2 885,1 667个差异表达基因。GO分析发现差异基因与细胞组分分类中的膜、膜的固有成分、膜的组成成分等条目密切相关,KEGG分析差异表达基因主要富集在植物激素信号转导、植物MAPK信号通路、植物昼夜节律、氨基糖和核苷酸糖的代谢、类黄酮生物合成、光合生物碳固定、角质和亚硫酸盐及蜡的生物合成、卟啉和叶绿素的代谢、淀粉和蔗糖代谢以及甘油酯质代谢等代谢途径。根据富集结果筛选8个差异表达基因表达分析,经实时荧光定量验证,基因表达趋势与转录组测序结果相符。【结论】CAO和POR基因表达延后,推测喷水处理可以延长叶绿素的合成时间;FBPase和SBPase的基因表达略微降低,推测叶面喷水调控基因的表达量降至缓解芒果午休的最适范围内,有利于树体积累碳同化物,从而提高光合效率。