意大利蜜蜂幼虫肠道响应球囊菌胁迫的可变剪接基因分析

为探讨可变剪接基因在宿主响应球囊菌(Ascosphaera apis)胁迫过程中的作用,基于前期已获得的意大利蜜蜂(Apis mellifera ligustica)幼虫肠道转录组数据进行可变剪接基因分析。在正常的意蜂4日龄幼虫肠道(AmCK),球囊菌胁迫的意蜂4、5和6日龄幼虫肠道(AmT1、AmT2、AmT3)中共鉴定出55 895个可变剪接事件,以外显子跨越和可变3′端剪接为主。各样品的共有可变剪接基因数为8 157个,AmT1、AmT2、AmT3的特有可变剪接基因数分别为207、334和542个。GO分类结果显示共有可变剪接基因分布在50个GO条目。KEGG代谢通路富集分析结果显示,共有可变剪接基因富集在290个pathway,基因富集数最多的是内质网蛋白加工、RNA运输和碳代谢。上述分析结果揭示了可变剪接基因在意蜂幼虫肠道响应球囊菌胁迫过程中的重要作用,为深入研究可变剪接基因的功能打下了基础。     

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

【目的】白垩病是困扰养蜂生产的顽疾。目前,尚无利用二代测序技术研究中华蜜蜂(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),进一步分析表明中蜂幼虫肠道的部分细胞免疫对球囊菌的胁迫产生应答,而体液免疫不产生应答,宿主的代谢相关基因受到球囊菌的显著抑制。【结论】揭示了中蜂幼虫在球囊菌入侵早期的胁迫应答,为深入解析中蜂幼虫的胁迫应答机制提供了重要信息,也为在分子水平研究关键应答基因打下了基础。     

意大利蜜蜂6日龄幼虫肠道内球囊菌的差异表达基因分析

为检测蜜蜂球囊菌(Ascosphaera apis)在胁迫意大利蜜蜂(Apis mellifera ligustica,简称意蜂)幼虫后期的基因表达谱,利用RNA-seq技术对纯培养的球囊菌孢子及球囊菌胁迫的意蜂6日龄幼虫肠道进行深度测序。通过比较处理组和对照组得到21 361个差异表达基因(DEGs),包括15 306个上调基因和6 055个下调基因。GO富集分析结果显示,上调基因富集于39个GO条目(term),基因富集数最多的为细胞进程(2 416unigenes)、细胞(2 408 unigenes)和细胞组件(2 408 unigenes);下调基因富集于38个GO term,基因富集数最多的为代谢进程(1 227 unigenes)、细胞进程(1 185 unigenes)和细胞(1 131 unigenes)。KEGG代谢通路富集分析结果显示,上调基因富集在119个通路,其中基因富集数最多的是核糖体(221 unigenes),其次为内质网中蛋白加工(190 unigenes)和内吞作用(177 unigenes);下调基因富集在113个通路上,基因富集数最多的是核糖体(144 unigenes),其次为RNA转运(113 unigenes)和氨基酸生物合成(108 unigenes)。进一步分析发现,富集在MAPK信号通路上的64个基因上调表达,说明该通路在球囊菌胁迫后期被激活。研究结果不仅为揭示球囊菌在胁迫意蜂幼虫后期的病原-宿主互作提供了重要信息,也为阐明球囊菌致病的分子机制奠定了基础。     

意大利蜜蜂幼虫响应东方蜜蜂微孢子虫胁迫的高表达基因分析

为探究意大利蜜蜂(Apis mellifera ligustica,简称意蜂)幼虫响应东方蜜蜂微孢子虫(Nosema ceranae)胁迫的高表达基因(HEG)的作用,利用RNA-seq技术对正常的意蜂4、5和6日龄幼虫(Am4CK、Am5CK和Am6CK)及N.ceranae胁迫的意蜂4、5和6日龄幼虫(Am4T、Am5T和Am6T)进行测序。共得到268 765 544条原始读段,经过滤得到的有效读段数为266 318 442条,各样品的Q30均在92.42%及以上。根据FPKM值≥15的标准从上述6个样品中分别筛选出2 460、2 692、2 908、2 367、2 549和2 270个HEG。Venn分析结果显示,去除与对照组共有的HEG,Am4T、Am5T和Am6T的共有HEG数为3个,特有HEG数分别为373、210和79个。GO分类结果显示,Am4T、Am5T和Am6T的特有HEG分别能够注释26、27和16个GO条目,注释基因数最多的分别是代谢进程、结合和结合。KEGG分析结果显示这些特有HEG可分别富集在145、138和68条代谢通路,富集基因数最多的分别是过氧化物酶体、碳代谢和内质网蛋白加工。进一步分析结果表明N.ceranae胁迫可导致宿主的部分物质代谢、能量代谢和免疫代谢通路水平提高。从6个样品的共有HEG中随机选取10个进行RT-PCR验证,均能成功扩增出目的条带,说明本研究中的HEG真实存在。研究结果揭示了意蜂幼虫响应N.ceranae胁迫的HEG表达谱及作用,为深入解析宿主的胁迫应答机制提供了重要信息。     

意大利蜜蜂4、5和6日龄幼虫肠道发育过程中差异表达基因的趋势分析

为探讨意大利蜜蜂(简称意蜂)幼虫肠道发育过程中的基因表达谱和基因差异表达信息,研究基于前期已获得的高质量的意蜂幼虫肠道转录组数据,利用STEM软件对意蜂4、5和6日龄幼虫肠道的差异表达基因(DEGs)进行趋势分析,进一步通过相关生物学软件对显著上调趋势包含的DEGs进行GO与KEGG代谢通路富集分析及表达量聚类分析。趋势分析结果显示,5 920个DEGs共聚类在8个基因表达模式,其中935个DEGs聚类在1个显著上调趋势。GO富集分析结果显示,显著上调趋势包含DEGs富集于39个GO条目,富集基因数最多的为结合,其次是细胞进程和单组织进程。KEGG代谢通路富集分析结果显示,显著上调趋势包含DEGs富集于94个代谢通路,富集基因数最多的是嘌呤代谢、内质网蛋白质加工和内吞作用。表达量聚类分析结果显示,随着意蜂幼虫肠道发育时间的延长,代谢和免疫相关通路富集基因的表达水平逐渐增强。研究提供了意蜂幼虫肠道发育过程的基因表达谱和基因差异表达信息,揭示了基因表达规律,为深入研究意蜂幼虫肠道发育的分子机制打下了一定基础。     

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

利用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技术对中蜂5和6日龄幼虫肠道进行深度测序,得到的有效读段(clean reads)数分别为29290630与26636038,Q20分别为98.45%和98.38%。GO富集分析结果显示,Ac5中的高表达基因(HEGs)富集在44个GO term,其中基因富集数最多的为代谢进程(753 unigenes),其次是催化活性(679 unigenes)与结合(657 unigenes);Ac6中的HEGs富集在44个GO term,其中基因富集数最多的为代谢进程(792 unigenes),其次是结合(760 unigenes)和催化活性(744 unigenes)。KEGG代谢通路(pathway)富集分析结果显示,Ac5的HEGs富集在129个pathway,其中基因富集数最多的是核糖体(103 unigenes)、内质网蛋白质加工(68 unigenes)和碳代谢(67 unigenes);Ac6的HEGs富集在129个代谢通路,其中基因富集数最多的是核糖体(103 unigenes)、氧化磷酸化(88 unigenes)及内质网蛋白质加工(68 unigenes)。进一步对HEGs进行Venn分析,结果显示,二者共有的HEGs为2 015个,特有的HEGs分别为219和491个。研究结果不仅可提供中蜂幼虫肠道发育过程中的基因表达谱数据,也为分子水平上深入研究中蜂幼虫肠道发育提供重要信息。     

意大利蜜蜂响应东方蜜蜂微孢子虫胁迫的免疫应答

【目的】通过对意大利蜜蜂（Apis mellifera ligustica,简称意蜂）工蜂中肠响应东方蜜蜂微孢子虫（Nosema ceranae）胁迫的差异表达基因（differentially expressed gene,DEG）及免疫通路进行深入细致的分析,揭示宿主的细胞和体液免疫应答,为关键免疫应答基因的筛选及功能研究打下基础。【方法】基于前期获得的正常及东方蜜蜂微孢子虫胁迫的意蜂7日龄和10日龄工蜂中肠（Am7CK、Am7T、Am10CK、Am10T）转录组数据,按照FDR≤1,P≤0.05和|log₂ fold change|≥1的标准筛选出各组的DEG,利用相关生物信息学软件对DEG进行Pearson相关性分析、Venn分析、GO分类和KEGG代谢通路富集分析,进而对免疫通路富集基因进行统计和分析,利用实时荧光定量PCR（real-time quantitative PCR,RT-qPCR）验证转录组数据的可靠性。【结果】 差异分析结果显示,Am7CK vs Am7T比较组包含472个上调基因和385个下调基因;Am10CK vs Am10T比较组包含611个上调基因和360个下调基因。Venn分析结果显示,两个比较组特有的DEG分别为739和853个,共有的DEG为118个。GO分类结果显示,Am7CK vs Am7T比较组中上调和下调基因分别涉及23和29个功能条目,其中富集上调基因数最多的前5位分别为结合、催化活性、代谢进程、细胞进程和单组织进程;富集下调基因数最多的前5位分别为代谢进程、单组织进程、催化活性、细胞进程和结合。Am10CK vs Am10T比较组中上调和下调基因分别涉及36和26个功能条目,其中富集上调基因数最多的前5位分别为单组织进程、结合、细胞进程、催化活性和代谢活性;富集下调基因数最多的前5位分别为结合、细胞进程、催化活性、代谢进程和单组织进程。KEGG代谢通路富集分析结果显示,Am7CK vs Am7T比较组中上调和下调基因分别富集在38和33条代谢通路,富集上调基因数最多的前5条分别为胆汁分泌、内质网蛋白加工、泛素介导的蛋白水解、PI3K-Akt信号通路和神经营养因子信号通路;富集下调基因数最多的前5条分别为胞质DNA传感通路、嘌呤代谢、嘧啶代谢、RNA聚合酶和核糖体;涉及泛素介导的蛋白水解等3条细胞免疫通路,以及PI3K-Akt信号通路等7条体液免疫通路。Am10CK vs Am10T比较组中上调和下调基因分别富集在54和43条代谢通路,富集上调基因数最多的前5条分别为Hippo信号通路、药物代谢-细胞色素P450、P450对外源物质的代谢、泛素介导的蛋白水解和鞘脂类代谢;富集下调基因数最多的前5条分别为mRNA监测、鞘脂类信号通路、果糖和甘露糖代谢、半乳糖代谢和鞘脂类代谢;涉及泛素介导的蛋白水解等7条细胞免疫通路,以及NF-κB信号通路等2条体液免疫通路。RT-qPCR验证结果显示6个随机挑选的DEG的表达量变化趋势与测序结果一致,证实了本研究中测序数据的可靠性。进一步分析发现意蜂7日龄和10日龄工蜂中肠的NF-κB信号通路均被东方蜜蜂微孢子虫激活,随即启动了3种抗菌肽apidaecin、defensin-1和hymenoptaecin的合成,体现了它们在宿主抵御东方蜜蜂微孢子虫入侵中的重要性。【结论】在转录组水平解析了意蜂工蜂中肠对东方蜜蜂微孢子虫的免疫应答,揭示宿主在胁迫前期同时作出细胞和体液免疫应答,前者可能在抵御病原入侵方面发挥主要作用;宿主的细胞免疫在胁迫后期持续增强,但体液免疫较大程度地减弱;泛素介导的蛋白水解通路及富集DEG,NF-κB信号通路及富集DEG,以及抗菌肽编码基因apidaecin、defensin-1和hymenoptaecin可能在意蜂工蜂对东方蜜蜂微孢子虫的免疫应答中起到关键作用。     

中华蜜蜂和意大利蜜蜂耐寒性能差异比较

【目的】研究中华蜜蜂(Apis cerana cerana,简称中蜂)和意大利蜜蜂(Apis mellifera ligustica,简称意蜂)体内抗寒系统的组成,为探究中蜂和意蜂耐寒性能的差异提供理论依据。【方法】试验选取室外相同环境条件越冬的中蜂和意蜂各5群,于2015年12月20日每群随机选取蜜蜂立即测定两者的过冷却点、冰点,比较其耐寒性能差异。然后分别测定蜂体游离水、结合水、脂肪、糖原和蛋白质含量以及海藻糖酶的活性。利用高效液相色谱测定血淋巴中山梨醇、海藻糖、甘露醇的含量。在海藻糖的生物合成通路中确定4个重要调控基因,选择β-action和β-肌动蛋白分别作为中蜂和意蜂的内参基因,采用实时荧光定量PCR检测中蜂和意蜂海藻糖4个相关调控基因在越冬期时m RNA相对表达量的差异。【结果】越冬期中蜂和意蜂过冷却点差异显著,中蜂的过冷却点显著低于意蜂(P0.05),但二者冰点没有显著性差异(P0.05);越冬期中蜂和意蜂蜂体含水量存在显著差异,中蜂体内游离水含量显著低于意蜂(P0.05),而结合水含量显著高于意蜂(P0.05);越冬中期蜂体脂肪含量在两蜂种间差异不显著(P0.05);中蜂和意蜂蜂体和血淋巴蛋白质含量差异不显著(P0.05);意蜂在越冬期时体内糖原含量明显低于中蜂,并且差异显著(P0.05);中蜂血淋巴中甘露醇、山梨醇含量接近于意蜂的两倍,海藻糖含量也相对高于意蜂,差异显著(P0.05);中蜂蜂体和血淋巴的海藻糖酶活性均显著低于意蜂(P0.05);以中蜂为对照组,实时荧光定量PCR检测意蜂的海藻糖酶基因(Tre-2)和尿苷二磷酸葡萄糖焦磷酸化酶基因(UGP)表达量明显高于中蜂(P0.05);而基因Tre-1表达情况相反;海藻糖合成酶基因(TPS)表达量在中蜂和意蜂体内无显著差异(P0.05)。【结论】与意蜂相比,中蜂能够通过降低体内游离水含量,提高糖原和小分子糖含量,调节海藻糖酶活性和相关基因表达等方式降低体内过冷却点,提高自身耐寒性能。     

感染意大利蜜蜂工蜂的东方蜜蜂微孢子虫及其纯化孢子的高表达基因分析

为探究东方蜜蜂微孢子虫(Nosema ceranae)在侵染意大利蜜蜂(Aips mellifera ligustica)过程的基因表达谱及高表基因(highly expressed gene,HEG)的功能,本研究利用RNA-seq技术对N.ceranae感染的意蜂工蜂7和10 d工蜂中肠(NcT1和NcT2)及其纯化孢子(NcCK)进行高通量测序。根据FPKM值≧15的标准筛选出各组的HEG并进行Venn分析,进而通过GO分类和KEGG代谢通路富集分析对各组的共有HEG进行功能注释。测序得到的原始读段经严格质控后分别得到5 816 465、28 501 225和210 824 312条高质量的有效读段,分别从NcT1、NcT2和NcCK中筛选出1 263、1 233和1 511个HEG,其中有1 079个基因在各组均高量表达,分别有15和7个基因仅在NcT1和NcT2中高量表达。GO分类结果显示,共有HEG涉及9个生物学进程相关条目,10个细胞组分相关条目,以及5个分子功能相关条目。KEGG代谢通路富集分析结果表明,共有HEG富集在218条代谢通路,富集基因数最多的是核糖体、Epstein-Barr病毒感染、细胞周期、内质网蛋白加工和真菌核糖体生物合成。深入分析结果显示分别有5和5个共有HEG富集在MAPK和cAMP信号通路,表明此2条信号通路及富集的HEG可能在N.ceranae的环境应激和侵染过程中发挥重要作用。研究结果提供了N.ceranae感染意蜂工蜂过程中的基因表达谱,揭示了HEG的作用,为深入解析N.ceranae的侵染机制打下了初步基础。     

基于中华蜜蜂幼虫肠道转录组数据对东方蜜蜂基因组的基因结构优化及新基因鉴定

利用前期获得的中华蜜蜂(Apis cerana cerana)幼虫肠道转录组数据对东方蜜蜂(Apis cerana)基因组的已注释基因进行结构优化、对未注释的新基因进行预测、分析和鉴定:利用Cufflink软件对reads进行转录本重构,通过与参考转录本序列进行比较,共对3 366个东方蜜蜂的基因结构进行了优化;利用Cuffcompare软件将重构转录本与参考基因组进行比对,共预测出527个东方蜜蜂的新基因,其中234个新基因在Nr数据库中存在注释信息;随机选取12个新基因进行PCR鉴定,有11个能够扩增出符合预期的目的片段,表明多数预测出的新基因真实存在;GO分类和KEGG代谢通路分析结果显示上述新基因可能在东方蜜蜂的生长发育、物质代谢、遗传信息传递等方面具有重要功能。本研究结果丰富和完善了东方蜜蜂基因组的基因结构与功能注释信息,也为新基因的功能研究打下基础。     

中华蜜蜂酚氧化酶和丝氨酸蛋白酶相关基因 及其全长转录本发掘与分析

为深入开展中华蜜蜂相关研究提供参考信息和基础,利用已获得的纳米孔(Nanopore)测序数据对中华蜜蜂Apis cerana cerana的酚氧化酶(phenoloxidase, PO)和丝氨酸蛋白酶(serine protease, SP)相关基因和全长转录本进行发掘和分析。通过Blast工具将中华蜜蜂所有全长转录本比对KEGG和Nr数据库以筛选出PO和SP相关基因和全长转录本。采用gffcompare软件将PO和SP相关全长转录本与东方蜜蜂参考基因组(ACSNU-2. 0)已注释基因进行比较以优化结构。使用Astalavista软件鉴定PO和SP相关基因的可变剪接(alternative splicing, AS)事件,再利用IGV浏览器进行结构可视化。通过PCR验证AS事件的真实性。利用TAPIS pipeline预测和分析可变多聚腺苷酸化(alternative polyadenylation, APA)位点,并通过TBtool软件鉴定APA位点上游的基序(motif)。鉴定到中华蜜蜂PO和SP相关的42个基因与146条转录本。优化了16个参考基因组已注释PO和SP相关基因的结构,其中5′端延长和3′端延长的基因均有7个,5′端和3′端同时延长的基因有2个。鉴定到PO和SP相关的10个基因的389次AS事件,其中最丰富的AS类型是5′端可变剪接。PCR结果证实了2次AS事件的真实性。共鉴定到34个PO和SP相关基因含有1个及以上的APA位点,其中含有3个APA位点的基因数量最多;在APA位点上游鉴定到多个motif,一致性序列为：GGHKSYWSHHTRATWTCNBHDMRRYWYRTNYTVACNGCKGCDCAYTGYR。鉴定到中华蜜蜂PO和SP相关的42个基因和146条全长转录本,优化了东方蜜蜂参考基因组已注释的PO和SP相关基因结构,并发掘出PO和SP相关基因的389次AS事件和237个APA位点。     

利用转录组数据开发意大利蜜蜂的SSR分子标记

意大利蜜蜂(简称意蜂)是自然界重要的授粉昆虫,广泛用于世界各国的养蜂生产,具有很高的生态和经济价值。为了开发意蜂的SSR分子标记,利用MISA软件对基于前期获得的意蜂幼虫肠道转录组数据组装得到的12 115条unigenes进行搜索,共预测出分布于2 149条unigenes的6 312个SSR位点,其中主要的重复类型为二核苷酸重复(54.42%)和三核苷酸重复(32.49%),主要基元为AT/AT(31.2%)和AG/CT(18.6%)。进一步利用软件设计出18 444对特异性引物,随机选取24对引物对国内3个不同来源的意蜂幼虫样品进行SSR位点扩增,有22对成功扩增出目的片段。研究开发出的SSR分子标记可用于意蜂的种群遗传和分子进化等研究,结果表明利用转录组数据鉴定非模式生物SSR位点的方法可行且高效。     

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

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

捻转血矛线虫丙硫咪唑敏感虫株和耐药虫株miRNA-mRNA关联比较分析

旨在比较分析捻转血矛线虫敏感虫株和耐药虫株miRNA表达谱,探讨miRNA与捻转血矛线虫丙硫咪唑耐药性相关基因之间的调控机制。运用Illumina Hiseq2000平台进行测序并进行cDNA文库的构建,利用RNAhybrid、Miranda和TargetScan对测序得到的数据进行靶基因预测,并取交集作为miRNA的靶基因预测结果,使用Bowtie、RepeatMasker、MIREAP、Rfam、miRBase和DAVID等生物信息学分析软件和权威数据库系统,筛选出差异的miRNA以及与捻转血矛线虫耐药相关的miRNA并进行靶基因预测,同时对KEGG pathway和GO功能富集到的靶基因和可能参与调控的通路进行预测分析。结果表明,敏感虫株和耐药虫株共筛选显著差异表达的miRNA共294个,其中113个上调,181个下调。对miRNA和其靶向调节的mRNA进行关联分析,共关联到1 770个miRNA-mRNA对为负调控关系,其中涉及到274个差异的miRNA和603个差异的mRNA。显著差异表达的miRNA靶基因被注释到了1 430条GO terms和327条KEGG pathway,富集到FoxO信号通路(FoxO signaling pathway),mTOR信号通路(mTOR signaling pathway),PI3K-Akt信号通路(PI3K-Akt signaling pathway)等与耐药相关的一些抗性通路。综上,通过对miRNA表达谱分析以及miRNA-mRNA靶向分析,发现捻转血矛线虫敏感虫株和耐药虫株显著差异表达的miRNA和其对应的靶基因,并对部分耐药相关通路中存在的靶基因及显著差异的靶基因所在的通路进行筛选和分析,为后续进一步通过转录组学深层次挖掘捻转血矛线虫产生耐药性的关键基因及相关调控分子提供科学依据。     

西方蜜蜂Sirtuin蛋白家族基因鉴定及表达分析

Sirtuin蛋白是一类进化保守的NAD⁺依赖性脱酰基酶家族，可调控生物的新陈代谢、基因组稳定性和寿命。为探究西方蜜蜂Sirtuin蛋白家族基因种类、分布和分子特征，以及不同级型、不同发育时期的表达模式，基于西方蜜蜂全基因组数据，采用生物信息学方法对西方蜜蜂Sirtuin蛋白家族成员进行鉴定，预测家族基因的结构和染色体位置，分析家族基因的结构域和系统发育关系；采用qRT-PCR分析Sirtuin蛋白家族基因在西方蜜蜂雌性蜂不同发育时期的表达情况。结果表明，在西方蜜蜂全基因组中鉴定出10个Sirtuin蛋白基因，属于6个家族成员（ AmSirt1、 AmSirt2、 AmSirt4~ AmSirt7）， AmSirt1和 AmSirt2家族成员各包含3个基因。6个家族成员基因分别定位于不同的染色体上，其中 AmSirt1和 AmSirt2呈现典型的串联重复分布， AmSirt1、 AmSirt6和 AmSirt7主要定位于胞质/细胞核， AmSirt2和 AmSirt5主要定位于胞质， AmSirt4定位于线粒体。西方蜜蜂Sirtuin蛋白家族基因的编码氨基酸数为265~899个，分子质量为 29.45~102.91 ku，等电点为4.58~9.09，外显子数为4~9个，内含子长56~2 719 bp，Sirtuin蛋白家族基因包含12个保守结构域，并发现3个西方蜜蜂所特有的保守氨基酸位点。Sirtuin蛋白家族基因鉴定和系统发育关系表明，西方蜜蜂缺少 Sirt3基因，10个Sirtuin蛋白基因聚合为4类6个分支，与东方蜜蜂、大蜜蜂的同源性较高。6个家族基因在蜂王和工蜂不同发育时期均有表达，蜂王整体基因表达量高于工蜂，各家族基因表达量均在2日龄蛹达到峰值。综上表明，西方蜜蜂有6个Sirtuin蛋白家族成员，系统进化中分为4类，缺少Sirt3家族成员。推测在西方蜜蜂雌性蜂不同发育时期，Sirtuin蛋白基因表达受禁食行为、热量限制和器官分化的影响。     

ETV5基因敲除对小鼠肌肉基因表达谱的影响

【目的】与野生型小鼠相比较,采用CRISPR/Cas9技术制备的ETV5基因纯合敲除小鼠在表现出内源性精原干细胞消融的同时伴随着体型和体质的明显弱势,本研究的目的是探究ETV5基因敲除对小鼠肌肉表达谱的影响。【方法】采集6周龄的3只野生型公鼠和3只ETV5纯合敲除公鼠的肌肉组织并抽提RNA,进行转录组测序,并对测序结果进行生物信息学分析,对2组小鼠肌肉样品的差异表达基因进行聚类分析、GO和KEGG富集分析。【结果】2组小鼠的肌肉组织共筛选出了574个差异表达基因,其中,上调基因292个,下调基因282个。发现了多个基因影响ETV5敲除小鼠的生长发育,其中,包括影响小鼠肌肉发育的基因Amd1和影响脂肪累积的基因Chrna2。GO和KEGG分析富集到的通路大多与脂肪代谢和生长发育相关。【结论】本研究结果为ETV5敲除小鼠生长发育缓慢和延迟的分子机制提供了解释,为研究ETV5在生理和病理上的相关功能提供了参考。