微小RNA介导意大利蜜蜂工蜂对东方蜜蜂微孢子虫的跨界调控

【目的】东方蜜蜂微孢子虫（Nosema ceranae）感染意大利蜜蜂（Apis mellifera ligustica,简称意蜂）导致蜜蜂微孢子虫病。本研究结合前期已获得的miRNA和mRNA组学数据,通过生物信息学方法对意蜂工蜂中肠的差异表达miRNA（differentially expressed miRNA,DEmiRNA）靶向结合的东方蜜蜂微孢子虫的mRNA和差异表达mRNA（DEmRNA）进行预测、数据库注释和调控网络分析,以期在组学水平解析miRNA介导意蜂工蜂对东方蜜蜂微孢子虫的跨界调控机制。【方法】通过比较东方蜜蜂微孢子虫侵染7 d和10 d的意蜂工蜂中肠（AmT1、AmT2）和未受侵染的工蜂中肠（AmCK1、AmCK2）的miRNA组学数据筛选出宿主的显著性DEmiRNA,通过比较侵染意蜂工蜂中肠的东方蜜蜂微孢子虫（NcT1、NcT2）和东方蜜蜂微孢子虫纯净孢子（NcCK）的mRNA数据筛选出病原的DEmRNA。利用TargetFinder软件预测宿主显著性DEmiRNA靶向结合的病原mRNA和DEmRNA。利用相关生物信息学工具对上述靶DEmRNA进行GO和KEGG数据库注释。结合前期研究结果筛选出孢壁蛋白、极管蛋白、蓖麻毒素B凝集素、ABC转运蛋白、ATP/ADP移位酶和糖酵解/糖异生途径等毒力因子和能量代谢通路相关的病原DEmRNA及与其存在靶向结合关系的宿主显著性DEmiRNA,并构建和分析二者的调控网络。【结果】AmCK1 vs AmT1比较组中宿主的48条显著上调miRNA和36条显著下调miRNA分别靶向病原的1 345和1 046条mRNA;进一步分析发现,宿主的47条显著上调miRNA和34条显著下调miRNA可分别靶向NcCK vs NcT1比较组中病原的584条显著下调mRNA和265条显著上调mRNA,它们可分别注释到19和22个功能条目以及66和64条通路。AmCK2 vs AmT2比较组中宿主的56条显著上调miRNA和51条显著下调miRNA分别靶向病原的1 260和1 317条mRNA;进一步分析发现,宿主的52条显著上调miRNA和49条显著下调miRNA可分别靶向NcCK vs NcT2比较组中病原的587条显著下调mRNA和336条显著上调mRNA,它们可分别注释到20和23个功能条目以及64和65条通路。AmCK1 vs AmT1和AmCK2 vs AmT2比较组的8条共同显著上调miRNA和1条共同显著下调miRNA分别靶向NcCK vs NcT1和NcCK vs NcT2比较组中的144条共同显著下调和10条共同显著上调mRNA,可分别注释到18和13个功能条目以及38和7条通路。此外,AmCK1 vs AmT1和AmCK2 vs AmT2比较组中宿主的显著上调miRNA可靶向结合NcCK vs NcT1和NcCK vs NcT2比较组中与RNAi途径,孢壁蛋白和蓖麻毒素B凝集素等毒力因子,糖酵解/糖异生途径以及MAPK信号通路相关的病原下调表达mRNA。【结论】在东方蜜蜂微孢子虫的侵染过程中,意蜂工蜂中肠的DEmiRNA与病原的DEmRNA之间存在复杂的靶向结合关系以及潜在的跨界调控关系;宿主的DEmiRNA可能通过抑制或降解病原的RNAi途径、毒力因子、糖酵解/糖异生通路、ATP/ADP移位酶、ABC转运蛋白及MAPK信号通路相关靶DEmRNA影响病原的侵染和增殖。     

福州地区中华蜜蜂(Apis cerana cerana)工蜂春季活动规律初探

为探究中华蜜蜂工蜂在春季的活动规律,采用无线射频识别技术,对福州地区中华蜜蜂工蜂春季出巢活动进行跟踪研究.结果表明,中华蜜蜂工蜂在3日龄内不出巢活动;4至16日龄每天中午前后短暂出巢飞行2-3 min;17日龄之后,工蜂转变为采集蜂,出巢活动扩展至整个白天,出勤频率由开始的4-6次逐渐增加至30日龄的15次,单次出巢采集时间持续(32.61±10.04)min,每次巢内停留(28.75±10.79)min;此外,研究揭示试验期间成年工蜂的寿命为(26.89±5.51)d.     

鳖的人工养殖（下）

鳖的人工养殖（下）陈大福广西农业大学５３０００５四、鳖的人工繁殖亲鳖的选择亲鳖要选择健壮、无病、无损伤、反应灵活、５～Ｂ年龄、体重在１０００～２０００克以上的鳖作种用。鳖的雌雄鉴别（图３）。如下表；鳖的雌雄鉴别表图３鳖的雌雄鉴别（ａ）华鳖腹面观；（）...     

胁迫中华蜜蜂幼虫肠道的球囊菌及其体外培养的高表达基因分析

利用RNA-seq技术对胁迫中华蜜蜂(简称中蜂)6日龄幼虫肠道的球囊菌及其体外培养进行深度测序,根据基因的FPKM值筛选得到高表达基因(HEGs),进而对HEGs进行GO及KEGG代谢通路(pathway)富集分析.Illumina测序数据经质控和过滤得到100814558条有效读段(clean reads),平均Q30均在93.81%以上.GO富集分析结果显示处理组(Aac T)的HEGs富集于39个GO term,基因富集数最多的是细胞(1872 unigene)、细胞组件(1872 unigene)和细胞进程(1748 unigenes);对照组(Aa CK)的HEGs富集于37个GO term,基因富集数最多的是细胞(785 unigenes),其次是细胞组件(785 unigenes)和代谢进程(776 unigenes).KEGG pathway富集分析显示,Aac T的HEGs富集在119个代谢通路上,基因富集数最多的是核糖体(176 unigenes)、碳代谢(147 unigenes)及氨基酸的生物合成(134 unigenes);Aa CK的HEGs富集在110个代谢通路上,基因富集数最多的是核糖体(179 unigenes)、氨基酸的生物合成(70 unigenes)以及碳代谢(62 unigenes).深入分析发现,对于富集在MAPK信号通路上的高表达基因数量,胁迫中蜂幼虫肠道的球囊菌远多于体外培养的球囊菌,说明病原的该通路在胁迫后期被显著激活.     

基于RNA-seq数据大规模挖掘蜜蜂球囊菌的SSR分子标记

基于已获得的球囊菌转录组数据预测微卫星标记(SSR),并进行SSR位点的信息分析和SSR引物的挖掘.利用软件MISA对球囊菌转录组中42610条unigenes进行搜索,共预测出7968个SSRs,分布于5233条unigenes中,其中最主要的重复类型为三核苷酸重复(53.15%),其次为二核苷酸重复(32.32%)和四核苷酸重复(8.46%).二核苷酸重复中的基序主要是AG/CT(占总量的15.8%).针对所有的SSR位点,利用Primer Premier 5软件设计出6956对SSR特异性引物,随机选取20对引物对两个不同来源的球囊菌样品进行SSR位点扩增,共有6对引物成功扩增出符合预期的目的片段.研究结果表明,利用球囊菌转录组数据开发SSR引物是可行的.     

中华蜜蜂幼虫肠道响应球囊菌早期胁迫的转录组学

【目的】白垩病是困扰养蜂生产的顽疾。目前,尚无利用二代测序技术研究中华蜜蜂(Apis cerana cerana,简称中蜂)幼虫白垩病的报道。本研究利用RNA-seq技术对健康(Ac CK)及球囊菌(Ascosphaera apis)胁迫的中蜂4日龄幼虫肠道(Ac T)进行深度测序,在转录组水平研究中蜂幼虫在球囊菌胁迫早期的胁迫应答。【方法】通过Illumina Hi Seq 2500平台对Ac CK和Ac T进行双端(PE125)测序,首先对测序数据进行质控和评估,利用edge R软件进行差异表达基因(DEG)分析,进而对DEGs进行GO富集分析及KEGG代谢通路富集分析,最后,利用实时荧光定量PCR(q RT-PCR)验证测序数据的可靠性。【结果】Ac CK和Ac T的转录组测序共得到188 457 338条原始读段(raw reads),经过滤得到182 088 448条有效读段(clean reads),两端Q20与Q30均在97.96%和94.97%以上,说明测序数据质量良好;主成分分析(PCA)结果显示第一与第二主成分可分别解释基因表达总体差异的75.8%和10.7%;DEG分析结果显示上调基因与下调基因的数量分别为344和239个;GO富集分析结果显示DEGs共富集在36个GO分类(term)上,其中基因富集数最多的细胞(106 unigenes)、细胞组件(106 unigenes)和代谢进程(104 unigenes);KEGG代谢通路富集分析结果显示上调与下调基因分别富集在72个和45个代谢通路上,其中上调基因富集数最多的是核糖体(72 unigenes)、碳代谢(16 unigenes)和糖酵解(14 unigenes),而下调基因富集数最多的是碳代谢(9 unigenes)、二羧酸代谢(8 unigenes)和氨基酸生物合成(7 unigenes),进一步分析表明中蜂幼虫肠道的部分细胞免疫对球囊菌的胁迫产生应答,而体液免疫不产生应答,宿主的代谢相关基因受到球囊菌的显著抑制。【结论】揭示了中蜂幼虫在球囊菌入侵早期的胁迫应答,为深入解析中蜂幼虫的胁迫应答机制提供了重要信息,也为在分子水平研究关键应答基因打下了基础。     

蜜蜂球囊菌基因结构优化及新基因鉴定

为补充和完善蜜蜂球囊菌Ascosphaera apis基因组序列和功能注释信息,本研究利用RNA-seq技术对球囊菌菌丝和孢子进行深度测序,并基于高质量的转录组数据对已注释基因进行结构优化,对未注释基因进行预测、鉴定和分析。将测序得到的有效读段进行参考基因组比对和转录本重构;利用Cuffcompare软件将重构转录本与参考基因组进行比对。结果表明:共对101个已注释基因的5′端或3′端进行了延长;共鉴定出373个新基因,随机挑选10个新基因进行RT-PCR验证,其中8个能扩增出符合预期的目的片段,表明预测出的多数新基因真实存在;共计147个球囊菌新基因可注释到Nr和eggNOG数据库中,85个新基因注释到29个GO条目,66个新基因注释到33条代谢通路。本研究结果补充了球囊菌的基因结构和功能注释信息,也为新基因的功能研究打下了初步基础。     

基于比较转录组学分析揭示中华蜜蜂及意大利蜜蜂幼虫的球囊菌抗性差异机制

为探究中华蜜蜂(Apis cerana cerana,以下简称中蜂)与意大利蜜蜂(Apis mellifera ligustica,以下简称意蜂)幼虫的球囊菌(Ascospahaera apis)抗性差异产生的原因,利用Venn分析、GO分类分析、KEGG代谢通路分析以及Ka/Ks分析对二者的转录组进行比较研究。Venn分析结果显示,中蜂与意蜂的同源基因有17 656个,归属于8 111个基因家族,二者的特有基因分别有1 158和241个,分别归属于468和86个基因家族。GO分类结果显示,中蜂和意蜂的特有基因分别富集在38和28个GO条目。KEGG代谢通路富集分析结果显示,中蜂和意蜂的特有基因分别富集于96和21个代谢通路;进一步分析发现中蜂的特有基因富集在内吞作用、溶酶体、泛素介导的蛋白水解和黑化作用等4个细胞免疫通路,以及MAPK信号通路和Toll-like受体信号通路等2个体液免疫通路,而意蜂仅有1个特有基因富集在MAPK信号通路。Ka/Ks分析结果显示受到强烈正向选择、弱正向选择、中性选择和负向选择的单拷贝同源基因分别有37、281、2 630和3 269个。对受到强烈正向选择的基因进行GO分类和KEGG代谢通路富集分析,GO结果显示中蜂和意蜂的单拷贝同源基因富集的GO条目类型相同;KEGG结果显示中蜂的单拷贝基因仅富集在氧化磷酸化。上述结果表明免疫相关基因数量的差异是二者的球囊菌抗性差异产生的重要原因之一,中蜂在与球囊菌的协同进化过程中可能通过提高能量利用从而限制球囊菌的增殖。本研究结果可为阐明中蜂及意蜂幼虫的球囊菌抗性差异产生的分子机制提供参考。     

中华蜜蜂幼虫肠道中微小RNA的鉴定及分析

[目的]利用small RNA-seq(sRNA-seq)技术和生物信息学方法对中华蜜蜂(Apis cerana cerana,简称中蜂)幼虫肠道的微小RNA(microRNA,miRNA)进行全转录组鉴定和分析,旨在丰富中蜂的miRNA信息,并为深入研究miRNA调控中蜂幼虫肠道发育的分子机理提供依据.[方法]利用s...     

东方蜜蜂微孢子虫RCDP42基因的分子克隆与生物信息学分析

对东方蜜蜂微孢子虫RCDP 42基因进行分子克隆,并通过生物信息学解析RCDP 42蛋白的分子特性,进而对东方蜜蜂微孢子虫和其他物种的RCDP 42蛋白进行保守基序和结构域预测及分析。结果表明,成功克隆到东方蜜蜂微孢子虫RCDP 42基因的CDS区,RCDP 42基因在孢子中真实表达,含有540个核苷酸,可编码179个氨基酸;RCDP 42蛋白的分子质量约为20.3 ku,分子式为C₉₁₅H₁₄₁₆N₂₃₀O₂₇₄S₁₀,脂溶系数为86.03,等电点为6.57;包含的亲水氨基酸数量多于疏水氨基酸,但不含信号肽与跨膜结构域,说明RCDP 42可能为亲水性蛋白、胞内蛋白和非跨膜蛋白,可能包含16个磷酸化位点、63个α-螺旋、43条延长链、11个β-转角及62个无规则卷曲。在东方蜜蜂微孢子虫、蚱蜢脑炎微孢子虫和肠脑炎微孢子虫的RCDP 42蛋白中均鉴定到3个保守基序(motif 1、motif 2和motif 3)与1个结构域(Rho)。