东方蜜蜂微孢子虫侵染意大利蜜蜂工蜂过程中nce-miR-10660及其靶基因表达谱

【目的】本研究旨在为探究nce-miR-10660调控东方蜜蜂微孢子虫Nosema ceranae侵染的作用机制提供理论和实验依据。【方法】采用Stem-loop RT-PCR对前期鉴定到的东方蜜蜂微孢子虫nce-miR-10660进行表达验证，再通过Sanger测序验证nce-miR-10660的序列。利用相关生物信息学软件预测和分析nce-miR-10660的靶基因。通过RT-qPCR检测nce-miR-10660及其靶基因在东方蜜蜂微孢子虫侵染意大利蜜蜂Apis mellifera ligustica工蜂过程中中肠中的表达谱。【结果】Stem-loop RT-PCR和Sanger测序结果分别证实了nce-miR-10660在东方蜜蜂微孢子虫孢子中的表达和真实存在。靶向预测结果显示nce-miR-10660共靶向RRDRP和RCDP 42等9个基因；分别有2和6个靶基因可被分别注释到KEGG数据库中的4条通路和GO数据库中的23个条目。RT-qPCR结果显示，相较于东方蜜蜂微孢子虫侵染后1 d时意大利蜜蜂工蜂中肠中nce-miR-10660的表达量，东方蜜蜂微孢子虫侵染后2 d时意...     

东方蜜蜂微孢子虫侵染过程中nce-miR-34537和靶基因PIP5KI的表达谱

【目的】本研究旨在对前期鉴定到的nce-miR-34537进行表达和序列验证,预测nce-miR-34537的靶基因并明确其分子特性,进而检测nce-miR-34537及其靶基因在东方蜜蜂微孢子虫(Nosema ceranae)侵染意大利蜜蜂(Apis mellifera ligustica)工蜂过程的表达谱,为进一步探究nce-miR-34537调控东方蜜蜂微孢子虫侵染的功能和作用机制提供基础。【方法】通过Stem-loop-RT-PCR和Sanger测序验证nce-miR-34537的表达和序列。通过生物信息学软件预测nce-miR-34537的靶基因PIP5KI(I型磷脂酰肌醇4-磷酸-5-激酶基因)的理化性质等分子特性和保守基序,并构建基于氨基酸序列的系统进化树。采用RT-qPCR检测nce-miR-34537及其靶基因的表达谱。【结果】nce-miR-34537在东方蜜蜂微孢子虫孢子中真实存在和表达。nce-miR-34537共靶向PIP5KI等151个基因。PIP5KI蛋白的分子式为C₈₈₂H_{1 364}N₂₂₆     

东方蜜蜂微孢子虫侵染意大利蜜蜂工蜂过程中nce-miR-12220及其靶基因的表达谱

东方蜜蜂微孢子虫Nosema ceranae侵染成年蜜蜂导致蜜蜂微孢子虫病。本研究旨在验证东方蜜蜂微孢子虫nce-miR-12220的存在和表达,并检测nce-miR-12220及其靶基因在病原侵染意大利蜜蜂Apis mellifera ligustica (简称意蜂)工蜂过程的表达谱。Stem-loop RT-PCR和Sanger测序结果显示nce-miR-12220真实存在和表达。靶向预测结果显示nce-miR-12220共靶向KRAB-A和γ tubulin等15个基因。上述靶基因可注释到19个GO条目和3条KEGG通路。RT-qPCR结果显示,相较于接种后1 d (1 day post infection,1 dpi),nce-miR-12220在2 dpi上调表达,而在3-12 dpi阶段总体表现出显著下调表达的趋势。类似地,与1 dpi相比,靶基因KRAB-A在2 dpi上调表达,而在3-12 dpi阶段总体呈下调表达的趋势。另外,与1 dpi相比,靶基因γ tubulin在2-12 dpi阶段总体表现出显著下调表达的趋势。上述结果表明nce-miR-12220与KRAB-A和γ tubulin之间存在潜在的靶向结合和正向调控关系;东方蜜蜂微孢子虫通过下调表达nce-miR-12220抑制KRAB-A的表达进而促进增殖;意蜂工蜂可能通过抑制东方蜜蜂微孢子虫的γ tubulin表达抵御病原侵染。研究结果明确了nce-miR-12220及其靶基因KRAB-A和γ tubulin在东方蜜蜂微孢子虫侵染意蜂工蜂过程中的动态表达规律,为深入探究nce-miR-12220在病原侵染中的功能及调控机制提供了理论和实验依据。     

转录组分析揭示东方蜜蜂微孢子虫侵染意大利蜜蜂的分子机制

【目的】本研究旨在通过差异表达基因(differentially expressed gene, DEG)分析以及毒力因子和其他侵染相关因子分析,在转录组水平揭示东方蜜蜂微孢子虫Nosema ceranae侵染意大利蜜蜂Apis mellifera ligustica的分子机制。【方法】基于前期已获得高质量的东方蜜蜂微孢子虫纯化孢子(NcCK)及侵染意大利蜜蜂工蜂7和10 d的东方蜜蜂微孢子虫(分别为NcT1和NcT2)转录组数据,根据P≤0.05且|log₂(Fold change)|≥1的标准,通过比较分析筛选出NcCK vs NcT1, NcCK vs NcT2和NcT1 vs NcT2比较组的DEG。通过相关生物信息学软件对上述DEG进行Venn分析、GO分类和KEGG代谢通路富集分析。根据Nr和KEGG数据库注释信息和相关文献进行对东方蜜蜂微孢子虫的毒力因子和侵染相关因子的统计和分析。通过RT-qPCR验证转录组数据及DEG表达趋势。【结果】从NcCK vs NcT1, NcCK vs NcT2和NcT1 vs NcT2比较组分别鉴定出1 397, 1 ...     

东方蜜蜂微孢子虫孢子中微小RNA的鉴定与分析

【目的】丰富东方蜜蜂微孢子虫Nosema ceranae的微小RNA(microRNA, miRNA)信息,并为深入探究miRNA在病原孢子和病原侵染中的功能提供理论和实验依据。【方法】基于已获得的small RNA-seq数据,利用生物信息学软件对东方蜜蜂微孢子虫的纯净孢子中的miRNA进行鉴定和分析。采用茎环反转录PCR(stem-loop RT-PCR)检测已鉴定的miRNA的表达;通过分子克隆与Sanger测序验证miRNA的序列。使用TargetFinder软件预测这些miRNA的靶基因,并对靶基因进行数据库注释。根据miRNA与靶基因的靶向结合关系构建调控网络,再利用Cytoscape软件进行可视化。【结果】在东方蜜蜂微孢子虫孢子中共鉴定到10个miRNA;这些miRNA的长度分布介于21~25 nt,首位碱基表现出U偏向性,每一位碱基的偏向性差异明显。Stem-loop RT-PCR检测结果表明这10个miRNA均真实表达;Sanger测序结果证实了随机选取的其中2个miRNA的序列真实性。共预测出249个靶基因,其中分别有249, 118, 136和3个靶基因可注释到Nr,Swiss-Prot, KOG和eggNOG数据库。此外,分别有134和71个靶基因可分别注释到GO数据库的30个功能条目和KEGG数据库的54条通路。【结论】本研究揭示了东方蜜蜂微孢子虫孢子中miRNA的存在和表达;这些miRNA通过调控潜在靶基因的表达参与孢子的生命活动。     

东方蜜蜂微孢子虫N6-腺苷特异性甲基化转移酶基因的克隆、分子特性及表达模式

【目的】东方蜜蜂微孢子虫(Nosme ceranae)专性侵染成年蜜蜂导致微孢子虫病,给养蜂生产造成很大损失。目前,东方蜜蜂微孢子虫的N6-腺苷特异性甲基化转移酶(N6-adenine-specific methyltransferase,N6AMT)基因NcN6AMT的研究仍然缺失。本研究对NcN6AMT的编码序列(coding sequence,CDS)区进行克隆,并解析NcN6AMT蛋白的理化性质和分子特性,进而测定东方蜜蜂微孢子虫侵染意大利蜜蜂(Apis mellifera ligustica)和中华蜜蜂(Apis cerana cerana)工蜂过程中NcN6AMT的相对表达量,以期丰富NcN6AMT的信息,并为探究东方蜜蜂微孢子虫侵染过程NcN6AMT的功能及表观调控机制提供基础。【方法】采用Protparam和ProtScale软件对NcN6AMT进行等电点和亲水性分析。通过SignalP 5.0、NetPhos 3.1、TMHMM-2.0、SOPMA和SWISS-MODEL等软件分别预测NcN6AMT的信号肽、磷酸化位点、跨膜结构域、二级结构和三级结构。使用WoLF PSORT II软件预测NcN6AMT的亚细胞定位。根据N6AMT氨基酸序列,通过TBtools软件对智人(Homo sapiens)、小鼠(Mus musculus)、褐飞虱(Nilaparvata lugens)、兔脑炎微孢子虫(Encephalitozoon cuniculi)、肠脑炎微孢子虫(Encephalitozoon intestinalis ATCC 50506)、蚱蜢脑炎微孢子虫(Encephalitozoon romaleae SJ-2008)、美洲思普雷格孢虫(Spraguea lophii 42_110)、家蚕微孢子虫(Nosema bombycis CQ1)、隐生菱形藻(Nitzschia inconspicua)和东方蜜蜂微孢子虫(Nosema ceranae)的N6AMT进行结构域预测和分析。利用MEME软件和MEGA 11.0软件进行东方蜜蜂微孢子虫和其他物种N6AMT的保守基序预测及进化树构建。采用实时荧光定量聚合酶链式反应(real-time fluorescence quantitative polymerase chain reaction,RT-qPCR)检测NcN6AMT在东方蜜蜂微孢子虫侵染意大利蜜蜂和中华蜜蜂工蜂过程的相对表达量。【结果】通过PCR扩增出大小约500 bp的目的片段,克隆测序结果显示其与GenBank数据库收录的预测序列一致;NcN6AMT蛋白的分子量约为18.7 kDa,分子式为C845H1374N214O249S6,理论等电点为5.88,脂溶系数是119.76,不稳定系数为37.47,平均亲水系数为0.025,含166个氨基酸和15个磷酸化位点,不含典型的跨膜结构域和信号肽,可同时定位于细胞质、线粒体、细胞核和液泡膜;NcN6AMT含1个甲基转移酶小结构域(methyltransferase small domain,MTS),该结构域同样存在于家蚕微孢子虫和兔脑炎微孢子虫等8个其他物种的N6AMT;在东方蜜蜂微孢子虫、兔脑炎微孢子虫、肠脑炎微孢子虫和蚱蜢脑炎微孢子虫的N6AMT中均预测到5个相同的保守基序;NcN6AMT与家蚕微孢子虫、肠脑炎微孢子虫、兔脑炎微孢子虫和蚱蜢脑炎微孢子虫的N6AMT序列一致性达到70.92%;东方蜜蜂微孢子虫和家蚕微孢子虫的N6AMT在系统进化树上聚为一支;东方蜜蜂微孢子虫接种后1–4 d,NcN6AMT在意大利蜜蜂和中华蜜蜂工蜂中肠内均呈现先上升后下降的表达趋势。【结论】成功克隆到NcN6AMT基因的CDS区,明确了NcN6AMT蛋白的理化性质和分子特性,并揭示东方蜜蜂微孢子虫和家蚕微孢子虫的N6AMT蛋白具有较高的保守性,NcN6AMT在东方蜜蜂微孢子虫侵染意大利蜜蜂和中华蜜蜂工蜂的第一个增殖周期(1–4 dpi)内动态表达且均呈上升-下降的表达模式。     

东方蜜蜂微孢子虫N6-腺苷特异性甲基化转移酶基因的克隆、分子特性及表达模式北大核心

【目的】东方蜜蜂微孢子虫(Nosme ceranae)专性侵染成年蜜蜂导致微孢子虫病,给养蜂生产造成很大损失。目前,东方蜜蜂微孢子虫的N6-腺苷特异性甲基化转移酶(N6-adenine-specific methyltransferase,N6AMT)基因NcN6AMT的研究仍然缺失。本研究对NcN6AMT的编码序列(coding sequence,CDS)区进行克隆,并解析NcN6AMT蛋白的理化性质和分子特性,进而测定东方蜜蜂微孢子虫侵染意大利蜜蜂(Apis mellifera ligustica)和中华蜜蜂(Apis cerana cerana)工蜂过程中NcN6AMT的相对表达量,以期丰富NcN6AMT的信息,并为探究东方蜜蜂微孢子虫侵染过程NcN6AMT的功能及表观调控机制提供基础。【方法】采用Protparam和ProtScale软件对NcN6AMT进行等电点和亲水性分析。通过SignalP 5.0、NetPhos 3.1、TMHMM-2.0、SOPMA和SWISS-MODEL等软件分别预测NcN6AMT的信号肽、磷酸化位点、跨膜结构域、二级结构和三级结构。使用WoLF PSORT II软件预测NcN6AMT的亚细胞定位。根据N6AMT氨基酸序列,通过TBtools软件对智人(Homo sapiens)、小鼠(Mus musculus)、褐飞虱(Nilaparvata lugens)、兔脑炎微孢子虫(Encephalitozoon cuniculi)、肠脑炎微孢子虫(Encephalitozoon intestinalis ATCC 50506)、蚱蜢脑炎微孢子虫(Encephalitozoon romaleae SJ-2008)、美洲思普雷格孢虫(Spraguea lophii 42_110)、家蚕微孢子虫(Nosema bombycis CQ1)、隐生菱形藻(Nitzschia inconspicua)和东方蜜蜂微孢子虫(Nosema ceranae)的N6AMT进行结构域预测和分析。利用MEME软件和MEGA 11.0软件进行东方蜜蜂微孢子虫和其他物种N6AMT的保守基序预测及进化树构建。采用实时荧光定量聚合酶链式反应(real-time fluorescence quantitative polymerase chain reaction,RT-qPCR)检测NcN6AMT在东方蜜蜂微孢子虫侵染意大利蜜蜂和中华蜜蜂工蜂过程的相对表达量。【结果】通过PCR扩增出大小约500 bp的目的片段,克隆测序结果显示其与GenBank数据库收录的预测序列一致;NcN6AMT蛋白的分子量约为18.7 kDa,分子式为C845H1374N214O249S6,理论等电点为5.88,脂溶系数是119.76,不稳定系数为37.47,平均亲水系数为0.025,含166个氨基酸和15个磷酸化位点,不含典型的跨膜结构域和信号肽,可同时定位于细胞质、线粒体、细胞核和液泡膜;NcN6AMT含1个甲基转移酶小结构域(methyltransferase small domain,MTS),该结构域同样存在于家蚕微孢子虫和兔脑炎微孢子虫等8个其他物种的N6AMT;在东方蜜蜂微孢子虫、兔脑炎微孢子虫、肠脑炎微孢子虫和蚱蜢脑炎微孢子虫的N6AMT中均预测到5个相同的保守基序;NcN6AMT与家蚕微孢子虫、肠脑炎微孢子虫、兔脑炎微孢子虫和蚱蜢脑炎微孢子虫的N6AMT序列一致性达到70.92%;东方蜜蜂微孢子虫和家蚕微孢子虫的N6AMT在系统进化树上聚为一支;东方蜜蜂微孢子虫接种后1–4 d,NcN6AMT在意大利蜜蜂和中华蜜蜂工蜂中肠内均呈现先上升后下降的表达趋势。【结论】成功克隆到NcN6AMT基因的CDS区,明确了NcN6AMT蛋白的理化性质和分子特性,并揭示东方蜜蜂微孢子虫和家蚕微孢子虫的N6AMT蛋白具有较高的保守性,NcN6AMT在东方蜜蜂微孢子虫侵染意大利蜜蜂和中华蜜蜂工蜂的第一个增殖周期(1–4 dpi)内动态表达且均呈上升-下降的表达模式。     

东方蜜蜂微孢子虫感染对中华蜜蜂免疫基因表达和血淋巴中糖水平的影响

【目的】东方蜜蜂微孢子虫Nosemaceranae主要感染蜜蜂的消化系统,它被认为是引起蜂群崩溃综合症(Colony collapse disorder,CCD)的主要原因之一。东方蜜蜂微孢子虫对新寄主西方蜜蜂Apis mellifera的影响已得到广泛的研究,但它对原宿主中华蜜蜂Apis cerana cerana影响的报到则相对较少。本实验以意大利蜜蜂Apis mellifera ligustica来源的东方蜜蜂微孢子虫感染中华蜜蜂,分析感病中华蜜蜂的免疫基因表达和血淋巴中糖含量变化,评价东方蜜蜂微孢子虫对中华蜜蜂健康的影响。【方法】利用荧光定量检测样本的免疫基因(vitellogenin、abaecin、apidaecin、hymenoptaecin、defensinl、defensin2和eater)表达量,分析东方蜜蜂微孢子虫的感染能否引起中华蜜蜂免疫应答。通过高效液相色谱法对样本的血淋巴进行糖浓度定量分析,检测其血淋巴中葡萄糖、海藻糖含量。【结果】所研究基因中,hymenoptaecin基因表达量在14日龄时感染组和对照组差异显著、vitellogenin基因表达量在7日龄时感染组和对照组间差异显著、defensin2基因表达量在7日龄和14日龄时感染组和对照组均差异显著,其他基因的表达量在均未达到显著水平;血淋巴中的葡萄糖与海藻糖的浓度在处理组和对照组间差异不显著。【结论】东方蜜蜂微孢子虫侵染中华蜜蜂后,hymenoptaecin、vitellogenin,defensin2基因在相应的日龄可对病原物应答,然而血淋巴中的葡萄糖和海藻糖的浓度不因侵染与否变化。     

东方蜜蜂微孢子虫一新型孢壁蛋白的基因克隆、 表达及亚细胞定位

【目的】微孢子虫(Microsporidia)孢壁在孢子构成及孢子侵染宿主过程中扮演重要的角色。本研究旨在鉴定获得的东方蜜蜂微孢子虫Nosema ceranae新型孢壁蛋白,并进行基因克隆和原核表达,明确其亚细胞定位。【方法】通过在线软件对东方蜜蜂微孢子虫新型孢壁蛋白AAJ76_1400036761序列进行生物信息学分析。利用PCR法获取目的片段并将其克隆至原核表达载体pCold II中,利用IPTG诱导表达重组蛋白并通过镍柱亲和层析法纯化目的蛋白。以获得的重组蛋白为抗原免疫小鼠制备多克隆抗体,通过间接免疫荧光技术和免疫胶体金定位技术对该蛋白进行亚细胞定位分析;利用蛋白质免疫印迹法检测该蛋白与东方蜜蜂微孢子虫几丁质壳的互作。【结果】在MicrosporidiaDB数据库中获得AAJ76_1400036761基因序列,基因全长681 bp,编码226个氨基酸;预测等电点为6.84,分子量为26.19 kD。SDS-PAGE电泳和Western blot结果表明AAJ76_1400036761重组蛋白能够在大肠杆菌Eescherichia coli Rosetta中高量表达。Western blot结果表明,制备的多克隆抗体能够特异地识别东方蜜蜂微孢子虫总蛋白中的AAJ76_1400036761,说明其在成熟东方蜜蜂微孢子虫中有表达。亚细胞定位结果显示,AAJ76_1400036761定位于东方蜜蜂微孢子虫孢壁上。重组蛋白AAJ76_1400036761能够与蜜蜂微孢子虫的几丁质壳结合。【结论】AAJ76_1400036761蛋白在东方蜜蜂微孢子虫成熟孢子中有表达;该蛋白定位于东方蜜蜂微孢子虫孢壁上,为东方蜜蜂微孢子虫新的孢壁蛋白。本研究为深入研究该蛋白的生物学功能奠定了基础。     

ame-miR-14在意大利蜜蜂工蜂幼虫肠道发育过程的调控作用

【目的】本研究旨在揭示ame-miR-14在意大利蜜蜂Apis mellifera ligustica工蜂幼虫肠道发育过程的调控作用。【方法】通过Stem-loop RT-PCR和Sanger测序分别验证ame-miR-14在意大利蜜蜂工蜂6日龄幼虫肠道中的表达和序列真实性。饲喂ame-miR-14的模拟物(mimic-ame-miR-14)和抑制物(inhibitor-ame-miR-14)及其相应的阴性对照mimic-NC和inhibitor-NC对ame-miR-14分别进行过表达和敲降,利用RT-qPCR检测意大利蜜蜂工蜂4-6日龄幼虫肠道中ame-miR-14的表达量。利用生物信息学软件预测ame-miR-14的靶基因并进行相关分析。利用RT-qPCR检测ame-miR-14的过表达和敲降后其靶基因FoxO和Hedgehog在4-6日龄幼虫肠道中的相对表达量。【结果】ame-miR-14在意大利蜜蜂工蜂6日龄幼虫肠道中真实存在和表达。相较于饲喂mimic-NC,饲喂mimic-ame-miR-14后ame-miR-14表达量在意大利蜜蜂工蜂4-6日龄幼虫肠道中均为显著上调;相...     

中华蜜蜂幼虫肠道响应球囊菌胁迫的microRNA应答分析

【目的】蜜蜂球囊菌(Ascosphaera apis,简称球囊菌)是一种能够侵染中华蜜蜂(Apis cerana cerana,简称中蜂)幼虫的致死性真菌病原。微小RNA(microRNA,miRNA)可通过在转录后水平靶向抑制或降解mRNA而参与宿主与病原互作过程。本研究旨在对球囊菌胁迫的中蜂6日龄幼虫肠道的差异表达miRNA(DEmiRNA)及其靶基因进行深入分析,进而揭示DEmiRNA在中蜂响应球囊菌胁迫应答过程中的作用。【方法】利用Illumina MiSeq平台对正常及球囊菌胁迫的中蜂6日龄幼虫肠道(AcCK和AcT)进行测序,通过相关生物信息学软件预测DEmiRNA及其靶基因。通过Blast将靶基因注释到GO和KEGG数据库。利用Cytoscape软件构建DEmiRNA与其靶mRNA的调控网络。通过Stem-loop RT-PCR和qPCR验证测序数据的可靠性。【结果】本研究共预测出537个miRNA,其长度分布介于16–35 nt之间,且不同长度的miRNA首位碱基偏向性差异明显。通过Stem-loop RT-PCR证实了10个novel miRNA的表达。AcCK vs AcT比较组共有54个DEmiRNA,包含31个上调和23个下调miRNA,可分别靶向结合6170和8199个靶基因。GO分类结果显示上调和下调miRNA的靶基因分别涉及47和47个条目,富集基因数最多的皆为结合细胞进程和催化活性。KEGG代谢通路(pathway)富集分析结果表明上调和下调miRNA的靶基因分别富集在134和126条pathway,富集基因数最多的均为内吞作用和内质网中的蛋白质加工。调控网络分析结果表明,DEmiRNA及其靶mRNA形成十分复杂的调控关系;31个DEmiRNA可靶向结合51个与泛素介导的蛋白水解相关的mRNA,18个DEmiRNA可靶向结合14个与Jak-STAT信号通路相关的mRNA;miR-1277-x、miR-26-x、miR-27-y、miR-30-x、miR-6052-x等16个miRNA共同参与了上述两条免疫通路的调控。最后,随机挑选3个DEmiRNA进行qPCR验证,结果证明了测序数据的可靠性。【结论】本研究提供了中蜂幼虫肠道在球囊菌胁迫后期的miRNA的表达谱和差异表达信息,揭示了球囊菌与宿主之间在miRNA组学水平存在复杂的互作。miR-6052-x和miR-1277-x作为调控网络的核心可能通过影响细胞凋亡参与宿主的免疫防御,miR-26-x和miR-30-x可能通过调控Jak-STAT信号通路参与宿主的胁迫应答。本研究筛选出的关键DEmiRNA有望作为治疗白垩病的分子靶标。     

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

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

意大利蜜蜂工蜂幼虫nkd基因的生物信息学分析及功能研究

【目的】意大利蜜蜂(Apis mellifera ligustica,简称意蜂)是西方蜜蜂(Apis mellifera)的亚种之一。蜜蜂球囊菌(Ascosphaera apis)侵染意蜂幼虫导致白垩病。本研究对西方蜜蜂裸表皮蛋白(naked cuticle, Nkd)进行保守基序预测和系统进化分析,并通过RNAi明确nkd基因对意蜂工蜂幼虫体重及宿主响应蜜蜂球囊菌胁迫的免疫应答的影响,以期丰富西方蜜蜂基因nkd的信息,并揭示意蜂幼虫nkd的功能。【方法】通过MEME软件预测西方蜜蜂和其他9个物种Nkd蛋白的保守基序。采用MEGA X软件对西方蜜蜂及其他9个物种的Nkd蛋白进行系统进化分析。通过饲喂dsRNA对意蜂幼虫肠道内的nkd进行RNAi。使用电子天平对幼虫进行称重。利用RT-qPCR检测nkd基因的干扰效率及免疫基因的相对表达量。【结果】西方蜜蜂与东方蜜蜂、柑橘凤蝶、家蚕和金凤蝶的Nkd蛋白均含有3个保守基序(motif 1、motif 2 和motif 3),说明上述5个昆虫物种的Nkd具有较高的保守性。西方蜜蜂与东方蜜蜂的Nkd蛋白聚为一支,说明二者的亲缘关系近。与dsRNA-egfp组相比,dsRNA-nkd组5日龄和6日龄幼虫肠道内nkd的表达量均极显著下调(P<0.001),干扰效率分别为49.60%和56.40%。另外,dsRNA-nkd组幼虫体重较dsRNA-egfp组显著下降,说明nkd显著影响幼虫体重。RT-qPCR结果显示,4日龄幼虫肠道内abaecin、apidaecin、birc5、defensin-1和PGRP-S2均被激活表达;5日龄幼虫肠道内abaecin、apidaecin、birc5和defensin-1均被激活表达,PGRP-S2的表达受到抑制;6日龄幼虫肠道内abaecin被激活表达,而apidaecin、birc5、defensin-1和PGRP-S2的表达均受到抑制,说明上述5个免疫基因在宿主响应胁迫的过程中呈不同的表达趋势,均参与宿主的免疫应答,nkd与abaecin和apidaecin的表达存在负向调控关系。【结论】西方蜜蜂的Nkd蛋白含有3个保守基序(motif 1、motif 2和motif 3),西方蜜蜂与东方蜜蜂的Nkd蛋白亲缘关系最近,通过饲喂dsRNA能有效干扰意蜂工蜂幼虫肠道内nkd表达,nkd影响意蜂工蜂幼虫体重及宿主对蜜蜂球囊菌胁迫的免疫应答。     

基于转录组与翻译组的海洋食烷菌(Alcanivorax)烷烃代谢调控研究

【目的】食烷菌是海洋烃类降解优势菌,其烷烃代谢调控机制有待深入研究。本研究拟从食烷菌转录和翻译水平上认识烷烃降解的调控过程。【方法】分别以乙酸和正十六烷(C16)为唯一碳源与能源,获取柴油食烷菌(Alcanivorax dieselolei) B5菌株的转录组和翻译组数据,并整合数据计算得到该菌在2种碳源培养条件下基因的翻译效率。采用基因本体论(gene ontology, GO)和京都基因和基因组百科全书(Kyoto encyclopedia of genes and genomes, KEGG)对差异翻译和翻译效率基因进行功能和代谢通路注释。【结果】当以C16为唯一碳源与能源时,B5菌株烷烃代谢途径的关键基因在转录与翻译水平均大量提升,包括烷烃单加氧酶、细胞色素P450氧化酶、醇脱氢酶和醛脱氢酶等。KEGG富集结果表明,翻译水平显著上调基因参与了肽聚糖生物合成、脂肪酸降解、氯代烷烃降解、氧化磷酸化和生物膜形成等通路;翻译效率差异基因主要富集在铁载体非核糖体肽的生物合成、氧化磷酸化和不饱和脂肪酸的生物合成等途径。通过转录组和翻译组学的联合分析显示,为了适应烷烃氧化,B5有效地协调了转...     

Transcriptome analysis of Verticillium lecanii based on RNA-Seq technology

[目的]蜡蚧轮枝菌是一种防治植物病原物和害虫的生防菌,弄清其致病的分子基础具有重要的意义。[方法]利用Illumina HiSeq技术测序蜡蚧轮枝菌的cDNA文库,并进行RNA-Seq序列拼接和功能注释。[结果]共获得1634670条高品质reads,碱基GC含量48.50%。14856条Unigene经组装拼接后,在Nr、Swiss-Prot、COG和KEGG数据库中分别比对得到1467、9379、6518和4188条Unigenes。采用GO分类工具将21403条Unigenes分成24个功能组,主要包含在细胞组成（1369）、分子功能（1679）和生物学过程（1928）3个大类里,在COG数据库注释的基础上,把6518条Unigenes分成38个功能组,通过GO和COG丰度识别分类,聚类的主要是一般功能预测、次生代谢产物的合成、运输和分解代谢以及信号转导机制起作用的独立基因。能与KEGG路径匹配的108条Unigenes中大部分和新陈代谢途径及次生代谢产物的合成有关。[结论]这些结果将有助于找寻蜡蚧轮枝菌的致病因子,从而丰富其致病机理。     

LncRNA13164通过ace-miR-4968-y调节中华蜜蜂幼虫对蜜蜂球囊菌侵染的免疫应答

【目的】长链非编码RNA (long non-coding RNA,lncRNA)可通过竞争性内源RNA(competing endogenous RNA,ce RNA)等多种方式发挥重要的调控作用,但蜜蜂lncRNA的功能研究至今依然缺失,lncRNA在蜜蜂免疫应答中的作用尚不清楚。本研究旨在揭示中华蜜蜂(Apis cerana cerana,简称中蜂)幼虫响应蜜蜂球囊菌(Ascosphaera apis)侵染的免疫应答中lncRNA的调控功能及作用机制。笔者团队前期通过深度测序和生物信息学分析发现,lncRNA13164与ace-miR-4968存在靶向关系且潜在参与在中蜂幼虫对蜜蜂球囊菌侵染的应答。【方法】利用实时荧光定量PCR (real-time fluorescence quantitative polymerase chain reaction,RT-qPCR)检测蜜蜂球囊菌接种后中蜂幼虫肠道内lncRNA13164的表达谱。采用ILoc-Lnc RNA软件预测lncRNA13164的亚细胞定位。联用RNAhybrid、Miranda和Target Scan软件预测lnc...     

环链虫草Cordyceps cateniannulata对番茄生长和抗氧化酶活性的影响

【目的】评估环链虫草Cordyceps cateniannulata对植物促生和植物抗氧化酶活性的影响。【方法】本研究利用浸种法将环链虫草接种于番茄植物体,在接种后的第30天和60天,通过番茄株高、根长、地上和地下部分的干鲜重指标评价其对番茄生长的影响;在接种后第10、20、30、60和90天,通过选择性培养基分析其在番茄不同组织中的生存情况,使用形态学及DNA序列比对的方法检验所分离菌株与原有菌株的一致性。在处理后的第30天,检测番茄叶片中的过氧化物酶(POD)、过氧化氢酶(CAT)、超氧化物歧化酶(SOD)及丙二醛(MDA)含量,观察环链虫草对番茄的抗氧化酶活性影响。【结果】环链虫草可定殖于番茄幼苗且对番茄生长有显著促进作用,菌株对植株的定殖偏好性分别为根部>茎部>叶部。酶活检测结果表明,处理组番茄叶片防御酶活性均呈显著升高的趋势,其中POD、CAT、SOD活性分别比对照增加了52.21%、75.31%和158.59%,MDA含量下降了35.15%。【结论】环链虫草可以通过浸种的方法感染并定殖番茄幼苗的根、茎、叶,促进番茄幼苗的生长并提高番茄抗氧化酶活性,具有较好的田间...     

蜜蜂球囊菌的microRNA鉴定及其调控网络分析

【目的】本研究利用small RNA-seq技术对球囊菌的纯培养进行测序,对球囊菌的micro RNAs miRNAs)进行预测、鉴定和分析,进而构建miRNAs-mRNAs的调控网络。【方法】利用Illumina Hiseq Xten平台对球囊菌菌丝与孢子进行测序,通过相关生物信息学软件对球囊菌的miRNAs进行预测和分析,通过茎环(Stem-loop)PCR对部分miRNAs进行鉴定,利用Cytoskype软件构建miRNAs-mRNAs的调控网络。【结果】本研究共获得48268696条clean reads,预测出118个球囊菌的miRNAs,它们的长度分布介于18–25 nt之间,不同长度的mi RNA的首位碱基偏好性差异明显。Stem-loop PCR验证结果显示共有10个miRNAs能够扩增出符合预期的目的片段,说明多数miRNAs可能真实存在。共预测出6529个球囊菌miRNAs的靶基因,其中5725个能够注释到Nr、Swissprot、KOG、GO和KEGG数据库。进一步分析结果显示有24个靶基因注释在MAPK信号通路。Cytoskype软件分析结果显示球囊菌的miRNAs与mRNAs之间存在复杂的调控网络,绝大多数的miRNAs处于调控网络的内部且同时结合多个mRNAs。【结论】本研究率先对球囊菌的miRNAs及miRNAs-mRNAs调控网络进行全面分析,研究结果丰富了对球囊菌miRNAs的认识,为其基础生物学信息提供了有益补充,也为阐明球囊菌致病的分子机理打下了一定基础。