东方蜜蜂微孢子虫侵染意大利蜜蜂工蜂过程中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在病原侵染中的功能及调控机制提供了理论和实验依据。     

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

【目的】微孢子虫(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蛋白在东方蜜蜂微孢子虫成熟孢子中有表达;该蛋白定位于东方蜜蜂微孢子虫孢壁上,为东方蜜蜂微孢子虫新的孢壁蛋白。本研究为深入研究该蛋白的生物学功能奠定了基础。     

东方蜜蜂微孢子虫侵染意大利蜜蜂工蜂过程中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-23928及其靶基因表达谱

【背景】东方蜜蜂微孢子虫(Nosema ceranae)专性侵染成年蜜蜂中肠上皮细胞而导致的微孢子虫病给养蜂业造成严重损失。【目的】检测东方蜜蜂微孢子虫nce-miR-23928及其靶基因在侵染意大利蜜蜂(Apis mellifera ligustica)工蜂过程的表达谱,为深入探究nce-miR-23928在东方蜜蜂微孢子虫侵染中的功能及调控机制提供依据。【方法】通过RNAhybrid、miRanda和TargetScan软件预测nce-miR-23928的靶基因。使用BLAST工具将上述靶基因比对到基因本体论(geneontology,GO)、京都基因和基因组百科全书(Kyoto encyclopedia of genes and genomes, KEGG)、Nr和Swiss-Prot数据库以获得相应注释。采用实时荧光定量PCR(realtimequantitativePCR,RT-qPCR)技术检测nce-miR-23928及其靶基因在东方蜜蜂微孢子虫侵染意蜂工蜂过程中的相对表达量。【结果】相较于接种后1 d (1 day post infection, 1 dpi),nce-...     

东方蜜蜂微孢子虫类甲基转移酶样蛋白5基因NcMettl5的克隆、分子特性与表达模式

N⁶-甲基腺苷(N⁶-methyladenosine,m⁶A)是真核生物中分布广泛、含量丰富且作用关键的一类RNA修饰。东方蜜蜂微孢子虫Nosemaceranae专性侵染成年蜜蜂导致微孢子病,给养蜂业造成很大损失。本研究对东方蜜蜂微孢子虫的甲基转移酶样蛋白5(Methyltransferase-likeprotein5,Mettl5)基因NcMettl5进行克隆,通过生物信息学预测NcMettl5蛋白的理化性质和分子特性,并通过RT-qPCR检测NcMettl5在东方蜜蜂微孢子虫侵染意大利蜜蜂Apis mellifera ligustica工蜂过程的相对表达量,以期丰富NcMettl5的信息,并为探究东方蜜蜂微孢子侵染中NcMettl5的功能及调控机制提供基础。结果表明,扩增出约600bp的目的片段,与GenBank数据库收录的预测出的NcMettl5序列一致;NcMettl5的分子量约为22.64kDa,分子式为C₁₀₃₉H₁₆₄₅N₂₆₉O_2...     

东方蜜蜂微孢子虫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)内动态表达且均呈上升-下降的表达模式。     

东方蜜蜂微孢子虫侵染过程中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₂₂₆     

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

【目的】东方蜜蜂微孢子虫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基因在相应的日龄可对病原物应答,然而血淋巴中的葡萄糖和海藻糖的浓度不因侵染与否变化。     

东方蜜蜂微孢子虫LCFAL 2基因的分子特性及系统进化分析

本研究旨在解析东方蜜蜂微孢子虫Nosema ceranae的长链脂肪酸辅酶A连接酶(long chain fatty acid COA ligase)基因LCFAL 2的分子特性,鉴定东方蜜蜂微孢子虫和其他物种的LCFAL 2蛋白的保守基序和结构域并进行系统进化分析,以期丰富东方蜜蜂微孢子虫LCFAL 2的信息,为深入开展相关功能研究提供基础。通过Expasy网站上的相关软件预测和分析LCFAL 2的理化性质、跨膜螺旋域、信号肽、磷酸化位点、二级结构和三级结构。使用MEME软件预测东方蜜蜂微孢子虫和其他物种LCFAL 2蛋白的保守基序。采用TBtools软件预测东方蜜蜂微孢子虫和其他物种LCFAL 2的结构域。通过Mega 11.0软件构建东方蜜蜂微孢子虫和其他物种的LCFAL 2的系统进化树。结果表明,东方蜜蜂微孢子虫LCFAL 2含有1 836个核苷酸,可编码611个氨基酸,分子式为C₃₁₄₅H₄₉₈₅N₈₁₃O₉₀₄S₂₃,分子量约为69.39 kDa,脂溶系数为96.8...     

参与调控意大利蜜蜂工蜂中肠基因表达的东方蜜蜂微孢子虫miRNA的组学解析及其调控网络

[目的]本研究结合前期已获得的miRNA和mRNA组学数据对东方蜜蜂微孢子虫Nosema ceranae 的差异表达 miRNA(differentially expressed miRNA,DEmiRNA)靶向意大利蜜蜂 Apis mellifera ligustica 工蜂中肠的 mRNA 和差异表达 mRNA(d...     

亚致死剂量吡虫啉胁迫对意大利蜜蜂工蜂嗅觉学习行为的影响及其脑部转录组分析

【目的】分析吡虫啉处理对意大利蜜蜂Apis mellifera ligustica嗅觉学习行为及脑部基因转录的影响,为新烟碱类杀虫剂对蜜蜂的负面影响提供依据。【方法】实验室条件下一次性饲喂意大利蜜蜂成年工蜂含有4 ng吡虫啉的50%蔗糖溶液,以饲喂不含吡虫啉的50%蔗糖溶液为对照,通过伸吻反应(proboscis extension response, PER)行为实验测定其对意大利蜜蜂成年工蜂嗅觉学习行为的影响;收集上述测试的意大利蜜蜂工蜂脑部提取RNA,进行RNA-Seq测序和生物信息学分析,利用实时荧光定量PCR挑选6个差异表达基因(DEGs)检测其在脑部的表达量以验证RNA-Seq测序结果。【结果】一次性饲喂意大利蜜蜂成年工蜂含4 ng吡虫啉的50%蔗糖溶液后,意大利蜜蜂的嗅觉学习能力与对照组(饲喂50%蔗糖溶液)相比显著降低。RNA-Seq测序结果显示饲喂意大利蜜蜂含4 ng吡虫啉蔗糖溶液后,处理组与对照组之间共有123个DEGs［校正后的P值(padj)<0.05)］,包括82个下调DEGs和41个上调DEGs。GO聚类分析发现下调DEGs主要富集在S-腺苷甲硫氨酸依赖性甲基转移酶活性、酸性磷酸酶活性、氧化还原酶活性、蛋白质异二聚化活性等,上调DEGs主要富集在跨膜受体活性、分子传感器活性、神经学系统过程等功能条目。KEGG富集分析显示下调DEGs主要富集在核糖体和溶酶体等细胞器、碳代谢和色氨酸代谢等代谢通路及Toll和Imd信号通路,上调DEGs未富集在KEGG代谢通路。荧光定量PCR结果显示测试的6个DEGs相对表达量变化趋势与RNA-Seq测序结果FPKM(fragments per kilobase million)值变化趋势一致,证明RNA-Seq测序结果可靠。【结论】亚致死剂量吡虫啉胁迫意大利蜜蜂成年工蜂显著降低意大利蜜蜂的嗅觉学习能力,影响意大利蜜蜂脑部免疫解毒等相关基因的表达和酶活性及氧化还原等生物代谢过程;亚致死剂量吡虫啉短时胁迫会刺激意大利蜜蜂的嗅觉感觉过程及神经信号传导过程。     

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

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

意大利蜜蜂哺育蜂学习记忆相关基因的转录组学分析

【目的】筛选与意大利蜜蜂Apis mellifera ligustica哺育蜂学习记忆密切相关的基因。【方法】在人工组建意大利蜜蜂蜂群中收集10日龄哺育蜂、21日龄哺育蜂、21日龄采集蜂,通过喙伸反应(proboscis extension reflex, PER)实验测定这3组样本之间学习记忆能力的差异。利用RNA-seq技术对具有学习能力和不具有学习能力的10日龄哺育蜂和21日龄哺育蜂脑中基因表达量进行全面分析,筛选出与哺育蜂学习记忆密切相关的差异表达基因(DEGs),并对这些差异表达基因进行GO和KEGG富集分析。qPCR检测随机选取的3个DEGs(上调基因TpnCⅢa和MED23以及下调基因Pkc)在具有学习能力和不具有学习能力的10日龄哺育蜂和21日龄哺育蜂脑中的表达量。【结果】PER实验结果显示,经过5次训练后,意大利蜜蜂21日龄哺育蜂的学习能力显著高于10日龄哺育蜂的,而21日龄哺育蜂和21日龄采集蜂的学习能力无显著差异;同样地,21日龄哺育蜂的记忆能力显著高于10日龄哺育蜂的,而21日龄哺育蜂和21日龄采集蜂的记忆能力无显著差异。RNA-seq分析筛选到88个与哺育蜂学习记忆密切相关的DEGs,其中18个上调表达,70个下调表达。GO富集结果显示,上调DEGs在生物学进程分类中富集的基因数最多,主要富集在信号转导、蛋白质加工修饰相关;下调DEGs也是在生物学进程分类中富集的基因数最多,主要富集在转录、信号转导、蛋白质生物合成相关,其中显著性富集在转录相关。KEGG富集结果显示,下调DEGs显著性富集在吞噬、光转导、AGE-RAGE信号通路。qPCR结果显示差异表达基因TpnCⅢa,MED23和Pkc在具有学习能力和不具有学习能力的10日龄哺育蜂和21日龄哺育蜂脑中的表达水平的变化趋势与RNA-seq数据中的变化趋势一致。【结论】PER实验表明日龄是影响意大利蜜蜂哺育蜂的学习和记忆能力的重要因素;同时本研究获得学习后哺育蜂脑部的基因表达变化趋势和富集分析,为深入研究哺育蜂学习记忆的分子机制奠定了重要的理论基础。     

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

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

球囊菌侵染对意大利蜜蜂幼虫肠道免疫与解毒相关基因表达及蛋白质、脂质和糖类含量的影响

[目的]球囊菌Ascosphaera apis是一种真菌性病原体,可以侵染意大利蜜蜂Apis mellifera ligustica幼虫,导致蜂群患白垩病.本研究通过球囊菌侵染意大利蜜蜂幼虫,探究对意大利蜜蜂幼虫肠道免疫解毒相关基因表达及蛋白质、脂质和糖类含量的影响,分析球囊菌侵染胁迫下意大利蜜蜂免疫应答反应与营养物质...     

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

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

【目的】本研究旨在揭示内蒙古草原新害虫沙葱萤叶甲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结果高度一致。【结论】糖孝解/糖异生、脂肪酸生物合成及钙离子信号通路可能在沙葱萤叶甲滞育调节中起着重要的作用。本研究为进一步研究沙葱萤叶甲成虫专性夏滞育的分子机理奠定了基础。     

ame-miR-79负调控意大利蜜蜂工蜂幼虫肠道中靶基因CYP450和FG的表达

【目的】本研究旨在通过饲喂ame-miR-79的模拟物(mimic)和抑制物(inhibitor)对意大利蜜蜂Apis mellifera ligustica工蜂幼虫肠道内的ame-miR-79分别进行过表达和敲减,探究ame-miR-79对幼虫肠道内靶基因表达的调控作用。【方法】通过分子克隆与Sanger测序验证意大利蜜蜂工蜂幼虫肠道内ame-miR-79的序列真实性。通过饲喂mimic-miR-79和inhibitor-miR-79对意大利蜜蜂工蜂4-6日龄幼虫肠道内的ame-miR-79分别进行过表达和敲减。采用相关生物信息学软件进行ame-miR-79的靶基因预测与分析。通过RT-qPCR检测ame-miR-79的过表达和敲减效果及过表达和敲减ame-miR-79后靶基因的相对表达量。【结果】ame-miR-79在意大利蜜蜂工蜂幼虫肠道内真实存在。与无义模拟物(nonsense mimic, mimic-NC)组相比,mimic-miR-79组的4-6日龄幼虫肠道内ame-miR-79的表达量皆极显著上调;与无义抑制物(nonsense inhibitor, inhibitor-NC)组相比,inhibitor-miR-79组的4和5日龄幼虫肠道内ame-miR-79的表达量显著下调,6日龄幼虫肠道内ame-miR-79的表达量极显著下调。ame-miR-79共靶向303个基因,涉及27个GO条目和179条KEGG通路。相较于mimic-NC组,靶基因细胞色素P450基因CYP450在mimic-miR-79组的4和5日龄幼虫肠道内均极显著下调表达,而在6日龄幼虫肠道内上调表达;靶基因fringe糖基化转移酶(fringe glycosyltransferase, FG)基因在4日龄幼虫肠道内下调表达,在5日龄幼虫肠道内显著下调表达,而在6日龄幼虫肠道内上调表达。相较于inhibitor-NC组,CYP450在inhibitor-miR-79组的4日龄幼虫肠道内极显著上调表达,在6日龄幼虫肠道内上调表达,而在5日龄幼虫肠道内下调表达;FG的表达水平在4日龄幼虫肠道内显著上调,在5和6日龄幼虫肠道内极显著上调。【结论】ame-miR-79在意大利蜜蜂工蜂幼虫肠道内真实存在;通过饲喂模拟物和抑制物能分别实现意大利蜜蜂工蜂幼虫肠道内ame-miR-79的有效过表达和敲减;ame-miR-79负调控意大利蜜蜂工蜂幼虫肠道内CYP450和FG的表达。