东方蜜蜂微孢子虫侵染意大利蜜蜂工蜂过程中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-...     

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

基于第三代纳米孔测序技术的东方蜜蜂微孢子虫全长转录组构建及注释

【目的】本研究旨在利用Oxford Nanopore测序技术组装和注释东方蜜蜂微孢子虫Nosema ceranae的高质量全长转录组。【方法】采用Nanopore PromethION系统对东方蜜蜂微孢子虫的纯净孢子进行转录组测序。通过识别每条clean read两端引物鉴定全长转录本序列。利用Blast工具将全长转录本比对Nr, Swiss-Prot, KOG, eggNOG, Pfam, GO和KEGG数据库,获得相应注释信息。分别利用蛋白结构域分析方法CPC, CNCI, CPAT和Pfam对长链非编码RNA(long noncoding RNA, lncRNA)进行预测,获得高可信度lncRNA。利用CPM(counts per million)法计算每一条全长转录本的表达量。【结果】利用Nanopore PromethION系统对东方蜜蜂微孢子虫转录组测序共测得6 988 795条raw reads,经质控获得6 953 469条clean reads,其中包含5 143 999条全长转录本。共鉴定到10 243条非冗余全长转录本,N50和平均读长分别为1 042 bp和89...     

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

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

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

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

【目的】蜜蜂是典型的具有发育狭温性的全变态昆虫。本研究以对低温最敏感的意大利蜜蜂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个基因;Profil...     

东方蜜蜂微孢子虫孢子中微小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通过调控潜在靶基因的表达参与孢子的生命活动。     

东方蜜蜂微孢子虫的SNP与InDel位点鉴定及分析

【目的】东方蜜蜂微孢子虫Nosema ceranae是一种专性侵染成年蜜蜂中肠上皮细胞的单细胞真菌病原,广泛感染世界各地的蜂群。本研究拟利用已获得的东方蜜蜂微孢子虫纯净孢子的高质量转录组数据进行单核苷酸多态性（Single nucleotide polymorphism,SNP）和插入缺失（Insertion-Deletion,InDel）位点的鉴定和分析,旨在丰富东方蜜蜂微孢子虫的SNP和InDel信息,并为新型分子标记的开发提供基础。【方法】使用GATK软件识别东方蜜蜂微孢子虫的SNP和InDel位点。采用SnpEff软件预测变异位点发生的基因组区域及变异产生的影响。通过相关生物信息学软件分别将SNP和InDel位点所在基因分别比对GO和KEGG数据库,获得相应的功能和通路注释。【结果】共鉴定到28195个SNP位点,其中发生转换和颠换的SNP位点分别有21403和6792个;上述SNP位点的突变类型有12种,其中最丰富的突变类型为C/T;分布在CDS区的SNP位点最多,其次是基因间区、上游区、下游区和内含子区;最丰富的密码子突变类型是同义突变;SNP位点所在基因可注释到代谢进程、细胞组分和催化活性等43个GO条目以及代谢途径、核糖体和等次生代谢产物的生物合成等85条KEGG通路。共鉴定到2831个InDel位点,其中分布在基因间区InDel位点最多,分布在CDS区的InDel位点最少;最丰富的密码子突变类型是移码突变;InDel位点所在基因可注释到细胞进程、细胞和结合等38个GO条目以及代谢途径、次生代谢产物的生物合成及核糖体等73条KEGG通路。【结论】东方蜜蜂微孢子虫中存在大量的SNP和InDel位点,SNP位点的突变类型主要为转换,与其他物种类似;SNP与InDel位点的基因组功能元件分布规律和突变类型具有明显差异;SNP和InDel位点所在基因与东方蜜蜂微孢子虫适应宿主细胞内环境及病原增殖过程具有潜在关系。     

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

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

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

[目的]分析吡虫啉处理对意大利蜜蜂Apis mellifera ligustica嗅觉学习行为及脑部基因转录的影响,为新烟碱类杀虫剂对蜜蜂的负面影响提供依据.[方法]实验室条件下一次性饲喂意大利蜜蜂成年工蜂含有4 ng吡虫啉的50％蔗糖溶液,以饲喂不合吡虫啉的50％蔗糖溶液为对照,通过伸吻反应(proboscis ex...     

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

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

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

【目的】筛选与意大利蜜蜂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个与哺育蜂学...     

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

【目的】筛选家蚕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的表达。