意大利蜜蜂幼虫肠道的miRNAs的生物信息学预测及分析

【目的】蜜蜂肠道是食物消化和营养吸收的主要部位,同时也是抵御病原侵染的重要场所。本研究旨在对意大利蜜蜂Apismelliferaligustica(简称意蜂)幼虫肠道的microRNAs(miRNAs)及其靶基因进行深入分析,进而解析miRNAs在幼虫肠道发育和生长过程中的作用。【方法】利用small RNA-seq(sRNA-seq)技术对意蜂4、5和6日龄幼虫肠道进行测序,通过相关生物信息学软件对意蜂幼虫肠道的miRNAs进行预测及分析。利用TargetFinder软件预测miRNAs的靶基因,然后将靶基因通过BLAST软件进行GO和KEGG数据库的功能注释。利用Cytoscape软件构建miRNAs与其靶向结合的mRNAs的调控网络。【结果】意蜂幼虫肠道样品的测序共获得96 329 456条有效标签序列(Clean tags),预测出560个miRNAs,包括45个novelmiRNAs。上述miRNAs的长度主要分布在17-27nt之间,且不同长度的miRNAs的首位碱基偏向性具有明显差异。表达量聚类分析结果显示,331个miRNAs为4、5和6日龄幼虫肠道所共有且表达量较为稳定,随着发育时间延长,特有miRNAs的表达水平呈总体增加的趋势。利用TargetFinder软件共预测出16479个意蜂幼虫肠道的靶基因,其中有8132和3361个可分别注释到GO和KEGG数据库,进一步分析发现有224个靶基因注释在Wnt信号通路,分别有2个和27个靶基因与保幼激素和蜕皮激素的调控密切相关。【结论】本研究在全基因组水平对意蜂幼虫肠道的miRNAs进行预测和分析,研究结果不仅丰富了意蜂的miRNAs信息,也为深入研究意蜂幼虫肠道的生长发育机理打下了初步基础。     

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

中华蜜蜂幼虫肠道响应球囊菌胁迫的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及其靶基因分析

【目的】蜜蜂球囊菌Ascosphaera apis(简称球囊菌)是一种特异性侵染蜜蜂幼虫肠道的致死性真菌病原。微小RNA(microRNA,mi RNA)是一种在转录后水平对m RNA进行负调控的关键调控因子。本研究旨在对意大利蜜蜂Apis mellifera ligustica(简称意蜂)幼虫肠道在球囊菌胁迫后期的差异表达mi RNAs(DEmiRNAs)及其靶基因进行深入分析,为揭示DEmiRNAs在胁迫应答中的作用提供重要信息。【方法】利用small RNA-seq(s RNA-seq)技术对正常及球囊菌侵染的意蜂6日龄幼虫肠道(分别表示为Am CK和Am T)进行测序,通过相关生物信息学软件对DEmiRNAs进行预测和分析。利用TargetFinder软件预测DEmiRNAs的靶基因,然后利用Blast软件对靶基因进行GO和KEGG数据库的功能注释。通过Cytoscape软件构建DEmiRNAs与其靶m RNAs的调控网络。通过茎环实时荧光定量PCR(stem-loop RT-qPCR)验证测序数据的可靠性。【结果】Am CK和Am T的测序分别得到9 230 496和10 823 667条有效序列标签; Am CK vs Am T比较组中包含15个上调和6个下调mi RNAs,分别结合3 503和3 252个靶基因,它们可分别注释到40和39个GO terms以及104和99条代谢通路;进一步分析发现调控网络中的17个DEmiRNAs靶向结合116个与丝氨酸蛋白酶相关的m RNAs,14个DEmiRNAs靶向结合54个与泛素介导的蛋白水解相关的m RNAs。Stem-loop RT-qPCR验证结果显示随机选择的4个DEmiRNAs (mi R-251-x,mi R-9277-y,mi R-1672-x和mi R-4968-y)的表达量变化趋势与测序结果一致,表明测序数据真实可信。【结论】本研究率先对意蜂幼虫肠道在球囊菌胁迫后期的DEmiRNAs及其靶基因进行预测与分析,提供了mi RNAs的表达谱和差异表达信息。ame-miR-927b,mi R-429-y和mi R-8440-y等DEmiRNAs可能参与对丝氨酸蛋白酶的调控和对泛素介导的蛋白水解的调控,DEmiRNAs与靶基因之间存在复杂的调控网络关系。DEmiRNAs可参与到意蜂幼虫肠道与球囊菌间的互作。     

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

蜜蜂球囊菌的参考转录组de novo组装及SSR分子标记开发

【目的】通过RNA seq技术对纯培养的蜜蜂球囊菌Ascosphaera apis孢子和球囊菌侵染的蜜蜂幼虫肠道组织进行测序,de novo组装球囊菌的参考转录组,并对其进行功能与代谢通路注释,进而基于该转录组数据开发球囊菌的SSR分子标记。【方法】首先通过差速离心获得活化的球囊菌孢子,配制含1×107孢子/m L的饲料饲喂4,5和6日龄的意大利蜜蜂Apis mellifera ligustica幼虫和中华蜜蜂Apis cerana cerana幼虫,通过Illumina Hi SeqTM2500平台同时对上述蜜蜂幼虫肠道及纯化球囊菌孢子进行深度测序,原始数据过滤后通过Trinity软件组装得到unigenes,进而通过BLASTX比对NCBI Nr,Swiss-Prot,KOG和KEGG数据库对unigenes进行功能与代谢通路注释。利用MISA软件对所有unigenes进行SSR搜索,并利用Primer Premier 5软件设计SSR特异性引物,通过PCR对不同来源的球囊菌SSR位点进行扩增。【结果】本研究共获得146 135 308条高质量reads,de novo组装得到42 609个unigenes。BLASTX比对结果显示,29 316个unigenes在上述公共数据库中具有功能和代谢通路注释。注释到法夫酵母Xanthophyllomyces dendrorhous上的unigenes最多,达6 050个。KEGG注释结果显示,unigenes可注释到117个代谢通路,其中富集在核糖体(ribosome)上的unigenes数量最多(529)。所有unigenes中共预测到7 968个SSRs,通过PCR开发出5个球囊菌SSR分子标记。【结论】本研究成功组装球囊菌的参考转录组,并进行了功能与代谢通路注释,可为在分子水平深入研究球囊菌提供重要的参考信息。基于该转录组信息开发的5个SSR分子标记可推动菌株鉴定、基因图谱构建及基因定位等研究。     

侵染中华蜜蜂6日龄幼虫的蜜蜂球囊菌的微小RNA差异表达谱及调控网络

【目的】蜜蜂球囊菌(Ascosphaeraapis,简称球囊菌)是专性侵染蜜蜂幼虫的致死性真菌病原。MicroRNA(miRNA)作为一类重要的基因表达调控因子,能够广泛参与真菌及其宿主的相互作用过程。本研究通过比较分析球囊菌孢子(AaCK)和侵染中华蜜蜂(Apis cerana cerana,简称中蜂) 6日龄幼虫肠道内的球囊菌(AaT)的smallRNA(sRNA)组学数据对球囊菌的差异表达miRNA(differentiallyexpressed miRNA,DEmiRNA)、靶mRNA及二者间的调控网络进行全面解析,旨在揭示miRNA介导的球囊菌对中蜂幼虫的侵染机制。【方法】对于球囊菌侵染的中蜂6日龄幼虫肠道的small RNA-seq (sRNA-seq)数据,利用BLAST工具连续比对东方蜜蜂(Apiscerana)和球囊菌的参考基因组筛滤得到AaT的sRNA组学数据。分别将AaCK和AaT的sRNA组学数据比对miRBase数据库,对球囊菌侵染宿主前后miRNA的数量和结构特征进行分析。联用RNAhybrid+svm_light、Miranda和TargetScan软件预测AaCK vs AaT比较组中DEmiRNA的靶mRNA,进而利用相关生物信息学软件对上述靶mRNA进行GO分类和KEGG代谢通路富集分析。通过Cytoscape软件对DEmiRNA-mRNA调控网络进行可视化。利用Stem-loop RT-PCR、RT-qPCR和分子克隆验证测序结果的可靠性。【结果】在AaCK和AaT中分别鉴定到380和387个miRNA。结构特征分析结果显示,AaCK和AaT的mi RNA皆集中分布在18–25 nt,且首位碱基主要偏向于U。AaCKvsAaT比较组共有270个DEmiRNA,包含155个上调miRNA和115个下调miRNA,分别靶向结合6091和6145个mRNA。GO分类结果显示,上述靶mRNA主要涉及代谢进程、细胞进程、应激反应等15个生物学进程;细胞、细胞组分、细胞器等12个细胞组分;催化活性、结合、转运子活性等11个分子功能。KEGG代谢通路富集分析结果显示,上述靶mRNA富集在123条代谢通路,参与对氨基酸代谢、碳水化合物代谢以及核苷酸代谢等物质代谢,氧化磷酸化、硫代谢、氮代谢等能量代谢,以及MAPK和Hippo等信号通路的调控。球囊菌DEmiRNA与靶mRNA之间存在复杂的调控关系,其中miR-29-x、miR-250-x、miR-4968-y、miR-11200-x、novel-m0023-5p、novel-m0130-5p和novel-m0135-5p等DEmiRNA可靶向结合与球囊菌的半胱氨酸蛋白酶、DNA甲基化转移酶以及几丁质酶相关的mRNA;此外,miR-7-x、miR-9-z、miR-319-y和miR-5951-y等同时参与调控MAPK信号通路;进一步分析发现,miR-250-x同时参与对DNA甲基化转移酶、MAPK信号通路及其他酶类合成与代谢途径的调控,并可能参与球囊菌与中蜂6日龄幼虫之间的跨界调控。通过Stem-loopRT-PCR和RT-qPCR验证了4个DEmiRNA的差异表达,并利用分子克隆和Sanger测序证实miR-7-x的序列与测序结果一致。【结论】本研究解析了侵染中蜂6日龄幼虫的球囊菌的miRNA差异表达谱及DEmiRNA的调控网络,揭示了球囊菌DEmiRNA可能通过调控病原的物质和能量代谢、增殖、毒力、信号通路及相关mRNA参与对中蜂幼虫的侵染过程。miR-7-x、miR-250-x、novel-m0023-5p等关键DEmiRNA有望作为白垩病治疗的新型分子靶点。     

猪伪狂犬病毒Fa株侵染对PK-15细胞microRNAs表达谱的影响

【目的】分析猪伪狂犬病毒Fa株(PRV-Fa)侵染对猪肾传代细胞PK-15 microRNAs(miRNAs)表达谱的影响。【方法】利用Illumina高通量测序技术,鉴定感染和非感染PRV-Fa的PK-15细胞的miRNAs;筛选并利用实时荧光定量RT-PCR(RT-q PCR)验证差异表达miRNAs;对差异miRNAs进行靶基因预测和Gene ontology(GO)分析。【结果】在感染和未感染PK-15细胞中分别检测到384个和405个miRNAs,其中感染PRV-Fa后差异表达的miRNAs共127个(60个上调,67个下调)。荧光定量结果显示差异miRNAs的表达趋势与高通量测序结果一致。GO分析显示,miRNAs广泛参与信号传导、细胞代谢、免疫反应、基因表达等生物学进程,其中miR-10b、miR-16、miR-18a、miR-19b、miR-20a、miR-145-5p、miR-146a、miR-181a、miR-499-5p等miRNAs与免疫相关。在靶基因调控网络图中,ssc-miR-30a-5p与ssc-miR-30d处于关键位置。研究鉴定出5个新的病毒编码miRNAs,其中PRV-miR-LLT2与PRV-miR-LLT4靶向PRV早期蛋白基因EPO。【结论】伪狂犬病毒Fa株感染对PK-15细胞编码miRNAs有显著影响。     

结合三代测序与二代测序技术揭示蜜蜂球囊菌孢子转录组的复杂性

【目的】蜜蜂球囊菌(Ascosphaera apis)是一种专性侵染蜜蜂幼虫的致死性真菌病原。本研究旨在利用PacBio单分子实时(singlemoleculereal-time,SMRT)测序技术对蜜蜂球囊菌孢子(AaS)中基因的可变剪切(alternative splicing,AS)和可变多聚腺苷酸化(alternative polyadenylation,APA)以及长链非编码RNA (long non-coding RNA,lncRNA)进行鉴定和分析,进而揭示蜜蜂球囊菌孢子中转录组的复杂性。【方法】采用Suppa软件对蜜蜂球囊菌孢子中基因的AS事件进行鉴定。通过RT-PCR对不同类型的AS事件进行验证。采用TAPIS pipeline对蜜蜂球囊菌孢子基因的APA位点进行鉴定。利用MEME软件分析孢子全长转录本的poly(A)剪接位点上游50bp的序列特征并鉴定motif。联用CPC和CNCI软件和比对Swiss-prot数据库的方法预测lncRNA,取三者的交集作为高可信度的lncRNA集合。进一步比较lncRNA和mRNA的转录本长度,外显子数量与长度,内含子长度,GC含...     

意大利蜜蜂工蜂中肠piRNA的鉴定与分析

【目的】Piwi蛋白互作RNA(PIWI-interacting RNA, piRNA)在昆虫的发育和免疫等重要生物学过程发挥重要的调控作用。本研究旨在丰富西方蜜蜂Apis mellifera的piRNA信息,并为进一步探究差异表达piRNA(differentially expressed piRNA, DEpiRNA)调控意大利蜜蜂Apis mellifera ligustica工蜂中肠发育的分子机理提供基础。【方法】基于已获得的意大利蜜蜂7日龄(Am7)和10日龄(Am10)工蜂中肠的small RNA(sRNA)组学数据,对piRNA进行预测和分析。将质控后的sRNA数据比对西方蜜蜂参考基因组,再将比对上的序列标签(tags)进一步比对数据库以滤除rRNA和tRNA等小非编码RNA(non-coding RNA,ncRNA),进而根据piRNA的长度特征鉴定piRNA。采用TPM(tags per million)算法对piRNA表达量进行计算和归一化处理。根据|log₂ fold change|≥1且P≤0.05的标准筛选Am7 vs Am10 比较组的DEpiRNA。通过相关软件预测DEpiRNA的靶mRNA并进行GO和KEGG数据库注释。利用Cytoscape软件对DEpiRNA-mRNA调控网络进行可视化。通过Stem-loop RT-PCR对随机选取的Am7与Am10两组共有的6个piRNA的表达进行验证。利用RT-qPCR对随机挑选的6个DEpiRNA的表达趋势进行验证。【结果】从意大利蜜蜂工蜂中肠中共鉴定到596个piRNA,长度介于24~33 nt,且Am7和Am10组中不同长度分布范围的piRNA的首位碱基偏向性具有明显差异。在Am7 vs Am10比较组共筛选出41个DEpiRNA,其中piR-ame-11093, piR-ame-1111451,piR-ame-190949和piR-ame-932156可分别靶向1 195, 1 018, 4 040和1 063条mRNA。上述靶mRNA可分别注释到45个功能条目和45条通路。Stem-loop RT-PCR验证结果显示6个piRNA(piR-ame-1084826, piR-ame-11093, piR-ame-14476, piR-ame-24995, piR-ame-39500和piR-ame-774987)均真实表达。RT-qPCR结果显示共有6个DEpiRNA (piR-ame-1084826,piR-ame-11093, piR-ame-14476, piR-ame-24995, piR-ame-39500和piR-ame-774987)的表达趋势与测序数据中的表达趋势一致,证实了本研究中sRNA-seq数据的可靠性和piRNA差异表达趋势的真实性。【结论】本研究从意大利蜜蜂工蜂中肠中鉴定到596个piRN,长度介于24~33 nt。不同长度分布范围的piRNA的首位碱基偏向性具有明显差异;piR-ame-11093, piR-ame-1111451,piR-ame-190949和piR-ame-932156具有靶向调控基因表达参与工蜂中肠发育过程的潜力。     

中华蜜蜂6日龄幼虫响应蜜蜂球囊菌侵染的长链非编码RNA应答研究

[目的] 本研究旨在探究长链非编码RNA（long non-coding RNA,lncRNA）在中华蜜蜂（Apis cerana cerana,简称中蜂）6日龄幼虫应答蜜蜂球囊菌（Ascosphaera apis,简称球囊菌）侵染过程中的差异表达谱及调控作用。[方法] 利用链特异性cDNA建库的RNA-seq技术对未被侵染及球囊菌侵染的中蜂6日龄幼虫肠道（AcCK和AcT）进行深度测序。通过相关生物信息学软件分析lncRNA的结构特征和表达谱。筛选并分析差异表达lncRNA（differentially expressed lncRNA,DElncRNA）的顺式（cis）作用及竞争性内源RNA（competing endogenous RNA,ceRNA）调控网络。采用RT-qPCR验证测序数据及DElncRNA差异变化趋势的可靠性。[结果] AcCK和AcT共预测出642个已知lncRNA和487个新lncRNA。与蛋白编码基因相比,上述中蜂lncRNA外显子数更少、长度更短且表达量更低。43个antisense lncRNA与40个正义链mRNA之间互补配对。AcCK和AcT比较组包含367个上调lncRNA和268个下调lncRNA。有194个DElncRNA潜在调控461个上下游基因,并涉及细胞进程、代谢进程和催化活性等38个功能条目以及氨基酸代谢、内吞作用和MAPK等191条通路。此外,有180个DElncRNA可靶向结合50个DEmiRNA,进而调控6365个mRNA;三者之间形成较为复杂的ceRNA调控网络。[结论] 中蜂6日龄幼虫肠道的部分lncRNA可作为antisense lncRNA参与应答球囊菌侵染;部分DElncRNA可通过cis作用调节物质代谢和免疫途径相关的上下游基因,从而介导宿主的侵染应答;TCONS_00010661和TCONS_00003104等DElncRNA可通过ceRNA网络调控Jak-STAT和氧化磷酸化等通路及富集基因,进而参与宿主的侵染应答。     

基于第三代长读段测序数据解析蜜蜂球囊菌基因的可变剪切与可变腺苷酸化

【目的】蜜蜂球囊菌(Ascosphaera apis,简称球囊菌)是一种专性侵染蜜蜂幼虫的真菌病原,导致的白垩病是严重影响养蜂生产的顽疾,每年给养蜂业造成较大损失。本研究旨在基于已获得的第三代长读段测序数据对球囊菌菌丝(Aam)和孢子(Aas)中基因的可变剪切(alternative splicing,AS)和可变多聚腺苷酸化(alternative polyadenylation,APA)进行深入分析。【方法】利用Astalavista软件鉴定Aam和Aas中基因的AS事件类型。利用IGV浏览器对部分剪切异构体(isoform)的结构进行可视化。通过TAPIS pipeline对Aam和Aas中基因的APA位点进行鉴定。通过MEME软件对Aam和Aas中全长转录本的APA位点上游50bp的序列特征进行分析并对motif进行鉴定。【结果】在Aam中共鉴定到286次AS事件,包括162次RI(Retained intron),87次A3(Alternative 3'splice-site)、32次A5(Alternative 5'splice-site)和5次SE (Skipping exon);在Aas中共鉴定到559次AS事件,包括305次RI、155次A3、85次A5、13次SE和1次MEE (Mutually exclusive exon)。进一步分析发现,现有参考基因组上的多数注释基因结构并不完整;部分注释基因在菌丝和孢子中转录形成的isoform在数量和结构方面均存在差异;部分isoform在参考基因组中没有对应的注释基因。Aam中共鉴定到2748个基因含有1个及以上的APA位点,其中含有1个APA位点的基因数量最多(726,26.42%);Aas中共鉴定到2768个基因含有1个及以上的APA位点,其中含有5个以上APA位点的基因数量最多(1180,42.63%)。部分基因在菌丝和孢子中含有不同的APA位点数。序列特征分析结果显示,球囊菌全长转录本的3'UTR的上下游表现出明显的碱基倾向性,U和A分别富集在3'UTR的上游和下游。此外,在球囊菌全长转录本的APA位点上游鉴定到4个motif,分别是UCUCCU、UCUUCU、CCCACC和CCCCCU。【结论】本研究通过对球囊菌菌丝和孢子中基因的AS和APA进行深入分析,揭示了球囊菌转录组的复杂性,为完善现有的基因组和转录组注释提供了宝贵信息,也为探究AS和APA在球囊菌的基因表达调控中的作用提供了关键基础。     

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...     

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

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

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

miRNA响应印度梨形孢定殖调控大麦生长发育的作用机制

【背景】内生真菌印度梨形孢(Piriformospora indica)定殖植物可以显著促进植物生长发育。miRNA已被证实在植物体的生长发育中具有调控作用。【目的】揭示印度梨形孢定殖大麦促进大麦生长发育过程中miRNA对印度梨形孢定殖的响应及对大麦生长发育的调控作用。【方法】提取大麦总RNA,实施转录组测序并进行序列比对与数据挖掘;使用高效液相色谱检测大麦生长素等激素水平变化。【结果】印度梨形孢对大麦有显著促生作用;全转录组测序结果显示：印度梨形孢侵染3 d较空白对照有18个差异表达的miRNA,其中11个miRNA上调、7个miRNA下调;侵染7 d与空白对照相比24个差异表达的miRNA,其中11个miRNA上调、13个miRNA下调;侵染3 d与侵染7 d相比有3个miRNA上调、6个miRNA下调。GO功能富集分析与KEGG通路分析显示,差异表达miRNA的靶基因主要参与转录、细胞分裂、生长素信号的感知和转导、光合作用和激素刺激响应。靶基因所参与的途径与大麦生长发育密切相关,暗示miRNA对印度梨形孢定殖过程做出了积极响应。代谢产物分析表明miRNA参与的调控路径的代谢产物发生改变。【结论】本研究以miRNA为入手点,探究了miRNA对大麦生长发育的调控机制,为揭示印度梨形孢的促生机制提供了新的研究方向。     

利用基因组学和MALDI-TOF MS技术鉴定放线菌纲细菌的核糖体蛋白质标志物

【目的】基质辅助激光解吸电离飞行时间质谱(MALDI-TOF MS)法基于微生物的特征蛋白指纹图谱鉴定菌种,本研究利用基因组学和MALDI-TOFMS技术鉴定放线菌纲细菌的核糖体蛋白质标志物。【方法】从MALDI-TOF MS图谱数据库选取放线菌纲代表菌种,在基因组数据库检索目标菌种,获取目标菌株或其参比菌株的核糖体蛋白质序列,计算获得分子质量理论值,用于注释目标菌株MALDI-TOFMS指纹图谱中的核糖体蛋白质信号。【结果】从8目,24科,53属,114种,142株放线菌的MALDI-TOFMS图谱中总共注释出31种核糖体蛋白质。各菌株的指纹图谱中核糖体蛋白质信号数量差异显著。各种核糖体蛋白质信号的注释次数差异显著。总共15种核糖体蛋白质在超过半数图谱中得到注释,注释次数最高的是核糖体大亚基蛋白质L36。【结论】本研究找到了放线菌纲细菌MALDI-TOF MS图谱中常见的15种核糖体蛋白质信号,可为通过识别核糖体蛋白质的质谱特征峰鉴定放线菌的方法建立提供依据。