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

【目的】白垩病是困扰养蜂生产的顽疾。目前,尚无利用二代测序技术研究中华蜜蜂(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条代谢通路。本研究结果补充了球囊菌的基因结构和功能注释信息,也为新基因的功能研究打下了初步基础。     

基于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,以下简称中蜂)与意大利蜜蜂(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-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信号通路上的高表达基因数量,胁迫中蜂幼虫肠道的球囊菌远多于体外培养的球囊菌,说明病原的该通路在胁迫后期被显著激活.     

蚕豆自噬基因鉴定及响应干旱胁迫分析

[目的]根据转录组测序结果鉴定蚕豆自噬(Augophagy,ATG)基因,分析其对干旱胁迫的响应,为深入研究该类基因在干旱胁迫下的作用机制提供理论参考.[方法]利用蚕豆干旱胁迫转录组测序(RNA-seq)数据,根据swiss-port、NCBI、Ensemble等数据库的注释,筛选和鉴定蚕豆ATG基因;通过Smart、...     

蜜蜂球囊菌菌丝和孢子中微小RNA及其靶mRNA的比较分析

【目的】蜜蜂球囊菌(Ascosphaera apis,球囊菌)专性侵染蜜蜂幼虫而导致白垩病。本研究旨在通过small RNA-seq（sRNA-seq）技术和生物信息学方法对球囊菌纯化菌丝（AaM）和纯化孢子（AaS）进行深度测序和比较分析,明确球囊菌菌丝miRNA和孢子miRNA的数量、结构和表达谱差异,并揭示菌丝和孢子共有miRNA、特有miRNA和差异表达miRNA（differentially expressed miRNA,DEmiRNA）及其靶mRNA与球囊菌菌丝和孢子生长、发育和病原致病性的潜在关系。【方法】实验室条件下获得纯培养的球囊菌,利用sRNA-seq技术对AaM和AaS分别进行测序,通过对原始读段（raw reads）进行过滤和质控获得有效标签序列（clean tags）。通过Venn分析筛选菌丝和孢子共有miRNA和特有miRNA。根据P≤0.05且|log₂ fold change|≥1的标准筛选AaM vs AaS的DEmiRNA。对上述共有miRNA、特有miRNA和DEmiRNA的靶mRNA进行预测,并对靶mRNA进行GO及KEGG数据库注释。根据靶向结合关系构建DEmiRNA和靶mRNA的调控网络。利用RT-qPCR验证测序数据的可靠性。【结果】AaM和AaS中分别得到12 982 320和12 708 832条raw reads,经过滤和质控分别得到10 800 101和9 888 848条clean tags。AaM中miRNA的长度介于18—26 nt,AaS中miRNA的长度介于18—24 nt,分布miRNA数量最多的长度均为18 nt,AaM和AaS中首位碱基为U的miRNA数量最多。AaM和AaS中表达量最高的miRNA均为miR6478-x、miR10516-x和miR482-x。菌丝和孢子共有miRNA靶向结合5 946个mRNA,二者特有miRNA分别靶向结合6 141和6 346个mRNA。共有miRNA的靶mRNA主要参与代谢进程、细胞进程和催化活性等42个功能条目,以及翻译、碳水化合物代谢和能量代谢等120条通路。AaM vs AaS比较组包含93个DEmiRNA,可靶向结合6 090个mRNA,这些靶mRNA可注释到38个功能条目和120条通路。DEmiRNA与靶mRNA之间形成较为复杂的调控网络,miR-4968-y位于调控网络的中心且能够靶向结合多达118个mRNA。RT-qPCR结果显示5个DEmiRNA的表达趋势与测序数据一致,证实了本研究中测序数据的可靠性。【结论】球囊菌菌丝和孢子中的miRNA具有类似的结构特征,但表达谱表现出明显差异;菌丝和孢子可能通过特异性表达和差异表达部分miRNA对其生长、发育和生殖进行调控。     

中华蜜蜂幼虫肠道参考转录组的de novo组装及SSR分子标记鉴定

【目的】利用RNA seq技术对中华蜜蜂(Apis cerana cerana,简称中蜂)幼虫肠道参考转录组进行de novo组装,并进行功能及代谢通路注释,进而利用该转录组数据进行中蜂幼虫的SSR分子标记鉴定。【方法】实验室条件下饲养中蜂幼虫,将纯化的蜜蜂球囊菌(Ascosphaera apis,简称球囊菌)孢子饲喂3日龄幼虫,剖取4、5和6日龄幼虫肠道,液氮速冻。将健康幼虫肠道与感染球囊菌的幼虫肠道同时进行Illumina测序。通过对raw reads的过滤得到clean reads,利用Trinity软件组装得到unigenes。通过BLASTx(E-value10-5)比对NCBI Nr、Swiss-Prot、KOG和KEGG数据库,对unigenes进行功能和代谢通路注释。利用MISA软件对所有unigenes进行SSR搜索,并利用Primer Premier 5软件设计特异性SSR引物,通过常规PCR对来源于北京、辽宁兴城和四川成都的中蜂幼虫肠道样本进行SSR位点鉴定。【结果】中蜂幼虫肠道的RNA seq共得到35 670 000条reads,de novo组装得到43 557个unigenes,平均长度为898 nt。共有18 225个unigenes可被注释到上述公共蛋白数据库,单独注释到NCBI Nr、Swiss-Prot、KOG和KEGG数据库的unigenes数分别为3 899、443、37和10个。KOG注释结果显示,11 442条unigenes分布于25个基因家族,其中注释到RNA加工和修饰家族的基因数最多,达1 249个。9 679个unigenes的GO分类结果显示,在生物学进程分类中,注释到细胞进程的基因最多,达4 201个,在分子功能和细胞组分类中,注释到结合与细胞的基因数最多,分别为4 935和2 900个。4 517个unigenes可注释到KEGG数据库中的216个代谢通路,注释到核糖体的基因数最多,达385个。利用MISA软件,可在7 763个unigenes搜索到13 448个SSR位点,随机选取20对SSR引物对国内3个不同来源的中蜂幼虫肠道样本的SSR位点进行扩增,有6对引物可鉴定出SSR分子标记。【结论】成功组装并注释了中蜂幼虫肠道参考转录组,可为中蜂及其幼虫的分子生物学及组学研究提供重要的参考信息,也可用于补充、丰富和检验东方蜜蜂的参考基因组,基于此转录组数据开发出6个中蜂的SSR分子标记,可应用于中蜂的基因图谱构建、基因多样性分析、基因定位等研究,也说明利用转录组数据开发非模式生物SSRs的方法可行。     

球囊菌侵染对意大利蜜蜂成年工蜂的影响

为了探究球囊菌侵染对意大利蜜蜂成年工蜂存活率、免疫基因、发育基因及肠道菌群的影响，将纯化培养的球囊菌孢子接种6日龄意大利蜜蜂成年工蜂，利用荧光定量PCR(qRT-PCR)检测球囊菌侵染对意大利蜜蜂成年工蜂免疫基因、发育基因及肠道菌群的影响。结果表明，连续3 d接种球囊菌孢子对意大利蜜蜂成年工蜂存活率无显著影响(χ²=0.133,P=0.715)。意大利蜜蜂工蜂中肠与后肠球囊菌孢子含量在接种球囊菌孢子4 d时分别显著高于接种8 d时(P<0.05)。与未接种组相比，接种球囊菌孢子4 d，意大利蜜蜂成年工蜂免疫相关基因abaecin、apidaecin、defensin-1、defensin-2、hymenoptaecin、GNBP-1和IMD均显著上调表达(P<0.05)；接种球囊菌孢子8 d，意大利蜜蜂成年工蜂免疫相关基因相对表达量无显著差异(P>0.05)；接种球囊菌孢子4 d和8 d，发育相关基因VGMC、usp、kr-h1、AMHex10869相对表达量无显著差异(P>0.05)。与未接种组相比，接种球囊菌孢子4 d，意大利蜜蜂成年工...     

蜜蜂患白垩病虫体内一株球囊菌拮抗细菌的分离与鉴定

【目的】寻找抑制蜜蜂球囊菌的拮抗菌,用于白垩病的生物防治。【方法】利用细菌纯化技术从患白垩病的蜜蜂幼虫体内分离纯化了一株对蜜蜂球囊菌(Ascosphaera apis)有抑制作用的细菌,并结合形态学、革兰氏染色以及16S rDNA序列分析技术进行鉴定,同时通过蜜蜂体内、体外接种试验研究其对蜜蜂球囊菌的抑制作用。【结果】该细菌被初步鉴定为蜡样芽胞杆菌(Bacillus cereus),它在培养基中能够完全抑制蜜蜂球囊菌的生长,但其发酵液对蜜蜂球囊菌的生长没有影响。将含菌饲料喂蜜蜂幼虫,蜜蜂幼虫的死亡率和生长发育没有受到影响,说明该细菌对蜜蜂幼虫没有致病性。蜜蜂幼虫体内同时接种细菌和球囊菌孢子,发现细菌对球囊菌在蜜蜂体内的萌发和初期菌丝的生长没有影响,蜜蜂依然能够患病,但生长后期该细菌能够降解球囊菌菌丝,对球囊菌子代孢子的产生表现出一定的抑制作用。【结论】 该研究结果为蜜蜂白垩病的生物防治提供了理论依据。     

基于纳米孔全长转录组数据完善东方蜜蜂微孢子虫的基因组注释

【目的】利用已获得的纳米孔全长转录组数据对现有的东方蜜蜂微孢子虫（Nosema ceranae）参考基因组的基因序列和功能注释进行完善。【方法】采用TransDecoder软件预测东方蜜蜂微孢子虫基因的开放阅读框（open reading frame,ORF）及相应的氨基酸。利用gffcompare软件将全长转录本与参考基因组注释的转录本进行比较,对基因组注释基因的非编码区向上游或下游延伸,修正基因的边界。利用MISA软件鉴定长度在500 bp以上的全长转录本的简单重复序列（simple sequence repeat,SSR）位点,包括单核苷酸重复、双核苷酸重复、三核苷酸重复、四核苷酸重复、五核苷酸重复、六核苷酸重复、混合SSR等类型。通过Blast工具将鉴定到的新基因和新转录本比对Nr、KOG、eggNOG、GO和KEGG数据库,从而获得功能注释。【结果】共预测出2 353个完整ORF,其中长度分布在0—100个氨基酸的ORF最多,占总ORF数的72.12%。共对东方蜜蜂微孢子虫的2 340个基因进行了结构优化,其中5′端延长的基因有1 182个,3′端延长的基因有1 158个。共鉴定到1 658个SSR,其中单核苷酸重复、双核苷酸重复、三核苷酸重复、四核苷酸重复的数量分别为1 622、23、7和6个;单核苷酸重复类型的SSR密度最大,达到182.32个/Mb,其次为混合SSR、双核苷酸重复和三核苷酸重复,分别达到6.90、2.78和0.73个/Mb。共鉴定出954个新基因,其中分别有951、333、371、422和321个新基因可注释到Nr、KOG、eggNOG、GO和KEGG数据库。此外,还鉴定出6 164条新转录本,其中分别有6 141、2 808、2 932、3 196和2 585条新转录本可注释到Nr、KOG、eggNOG、GO和KEGG数据库。新基因和新转录本注释数量最多的物种均为东方蜜蜂微孢子虫,其次是蜜蜂微孢子虫（Nosema apis）。【结论】研究结果较好地完善了现有的东方蜜蜂微孢子虫参考基因组已注释基因的序列和功能注释,并补充和注释了大量参考基因组未注释的新基因和新转录本。     

蜜蜂球囊菌中长链非编码RNA的调控作用

【目的】长链非编码RNA（long non-coding RNA,lncRNA）是一类二、三级结构高度保守,长度>200 nt且不具蛋白编码能力的RNA,在转录和转录后水平广泛参与调控剂量补偿、细胞分化和生长发育等生命活动。本研究基于前期获得的蜜蜂球囊菌（Ascosphaera apis,简称球囊菌）菌丝和孢子混合样品的高质量lncRNA组学数据进行球囊菌lncRNA的顺式（cis）作用、反义lncRNA（antisense lncRNA）作用和竞争性内源RNA（competing endogenous RNA,ceRNA）作用的分析和探讨,以期揭示lncRNA在球囊菌中的潜在功能。【方法】基于lncRNA基因在染色体上的位置,预测lncRNA上下游100 kb以内的蛋白编码基因;使用Blast软件将上下游基因比对到GO和KEGG数据库,以获得功能和通路注释。利用LncTar软件对反义lncRNA的靶mRNA进行预测,并使用Blast软件将上述靶mRNA比对到KEGG和eggNOG数据库。利用TargetFinder软件预测lncRNA靶向结合的miRNA及miRNA靶向结合的mRNA,根据靶向结合关系建立lncRNA-miRNA和lncRNA-miRNA-mRNA调控网络,进而通过Cytoscape v3.7.1软件进行调控网络的可视化。利用RT-PCR对调控网络中的lncRNA、靶miRNA和靶mRNA进行表达验证。【结果】共预测出371个lncRNA的5 852个上下游基因,可注释到细胞进程、代谢进程和应激反应等48个功能条目,以及新陈代谢途径、次生代谢产物的生物合成和抗生素的生物合成等121条通路,表明部分lncRNA可通过顺式作用调节上下游基因的表达,从而参与调控球囊菌的生长发育及物质能量代谢等基础生命活动。球囊菌的7个lncRNA与7个靶mRNA存在序列互补关系,其中5个mRNA在eggNOG数据库中仅注释为假定蛋白,仅gene3444在KEGG数据库注释为核孔复合体蛋白An-Nup120和假定蛋白,表明上述1个反义lncRNA可能参与调控核孔复合体蛋白的生物合成等生物学过程。此外,共预测出227个lncRNA与73个miRNA之间存在靶向结合关系,其中多数lncRNA(79.02%)仅能结合1—2个miRNA,部分miRNA可被多个lncRNA靶向结合;进一步构建氧化磷酸化通路和MAPK信号通路相关的lncRNA-miRNA-mRNA调控网络,分析结果显示氧化磷酸化通路相关的222个lncRNA靶向78个miRNA及50个mRNA,MAPK信号通路相关的222个lncRNA靶向76个miRNA及46个mRNA,表明部分lncRNA通过ceRNA作用调控此两条通路,从而影响球囊菌的能量合成、环境适应以及生长发育等过程。【结论】部分lncRNA可能通过顺式作用和ceRNA作用调节球囊菌的生长发育、物质能量代谢以及环境适应等生物学过程;MSTRG.5393.1可能作为反义lncRNA调控球囊菌的核孔复合体的蛋白合成。     

卵形鲳鲹卵巢不同发育时期的转录组测序分析

【目的】鉴定筛选出与卵形鲳鲹卵巢发育相关的候选基因及信号通路,为揭示其卵巢性成熟过程的分子机制打下基础。【方法】挑选卵巢发育处于?期和Ш期的雌性卵形鲳鲹,分别构建卵形鲳鲹卵巢?期和Ш期的cDNA文库,采用Illumina HiSeqTM 2500进行转录组测序,经过滤、质量控制及拼接组装后获得的Unigenes在七大数据库（Nr、Nt、Pfam、KOG/COG、Swiss-Prot、KEGG和GO）中进行比对;通过FPKM及DEGseq筛选出差异表达基因,以GOseq和KOBAS对差异表达基因分别进行功能注释及信号通路富集分析,并采用MISA和GATK3进行SSR鉴定及SNP分析。【结果】卵形鲳鲹卵巢组织转录组测序获得的325156432条Raw reads,经过滤筛选得到317206752条Clean reads,拼接组装后得到59554条Unigenes;69.65%的Unigenes在Nr、Nt、Pfam、KOG/COG、Swiss-Prot、KEGG和GO等七大数据库中注释成功,其中有24599条Unigenes被注释到GO数据库,15997条Unigenes被注释到KEGG数据库。在卵形鲳鲹卵巢组织的2个发育时期共鉴定获得56115个基因,经差异表达分析后获得17737个差异基因,其中8169个基因在卵巢Ш期上调表达、9568个基因在卵巢Ш期下调表达。GO功能注释分析发现,卵形鲳鲹卵巢差异表达基因主要注释在细胞过程、氮化合物代谢过程、初级代谢过程、核、核部分、离子结合及水解酶活性等条目上;而KEGG信号通路富集分析结果显示,17737个差异表达基因显著富集在318条代谢途径上,其中前20条KEGG信号通路包括2-氧代羧酸代谢、PI3K-Akt信号通路、甲状腺激素信号通路、磷脂酶D信号通路、Fc εRI信号通路和细胞周期等。卵形鲳鲹卵巢转录组（59554条Unigenes）中共存在30133个SSRs和82490个SNPs。【结论】GnRHR、FSHR、FSHβ、CYP11A、SIRT3和PEG3等差异表达基因及PI3K-Akt信号通路和VEGF信号通路等与卵形鲳鲹卵巢的发育密切相关,共同调节卵巢的发育与成熟,在卵巢性成熟过程中发挥重要作用。     

青花菜花蕾转录组测序分析及蜡粉合成相关基因挖掘

【目的】对青花菜花蕾进行转录组测序分析，并挖掘与蜡粉合成相关基因，为探明青花菜花球表面蜡粉形成的分子机制提供理论参考。【方法】分别提取野生型和蜡粉缺失型青花菜花球总RNA，采用Illumina HiSeqTM 2500平台进行转录组测序，获得高质量Clean reads，采用Trinity进行序列组装后获得青花菜Unigene库，将获得的Unigene序列与Nr、Nt、KEGG、Pfam、KOG/COG、Swiss-Prot和GO数据库比对，获得基因功能注释信息；使用DESeq2进行差异表达分析。【结果】共获得44.68 Gb Clean data，De novo组装得到41244条Unigenes，N50长度为1847 bp。从所获得的Unigenes中筛选出8685个差异表达基因（DEGs）（上调基因5747个，下调基因2938个），共有8038个基因被注释到不同数据库，其中，5220个基因注释到Pfam数据库； 2066个基因注释到COG数据库，3866个基因注释到KOG数据库； 2580个差异表达基因被注释到75个转录因子家族中，注释最多的是MYB家族（235个）；GO数据库中6095个差异表达基因注释到细胞组分、分子功能和生物学过程三大类的52个功能分类； KEGG数据库中，1671个差异表达基因富集到138条代谢通路，其中13个差异表达基因与脂肪酸合成有关，7个差异表达基因与蜡粉生物合成途径有关。【结论】转录因子MYB家族在调控青花菜蜡粉合成中发挥重要作用。蜡粉合成过程中相关酶基因的差异表达是调控青花菜蜡粉合成的关键，尤其是野生型和蜡粉缺失突变体中特异性表达的差异表达基因，可作为后续研究青花菜花球表面蜡粉形成分子机制的对象。     

凡纳滨对虾转录组测序分析及肌肉生长发育相关基因的筛选

【目的】筛选出凡纳滨对虾生长发育的功能基因及代谢调控网络,揭示其生长发育的分子机制,为后续开展凡纳滨对虾分子生物学研究提供宝贵的基因数据来源。【方法】以快速生长群体和慢速生长群体的凡纳滨对虾肌肉组织为研究材料,通过Illumina HiSeqTM2500平台对构建的cDNA文库进行高通量测序分析,以StringTie进行拼接组装后利用DESeq2筛选差异表达基因,并基于KOG、GO、nr、COG、Swiss-Pro、KEGG和Pfam等数据库进行差异表达基因功能注释分析。【结果】经拼接组装共获得53458844条Clean reads,各样品Clean reads的Q30均在93.00%以上;Illumina测序获得的Clean reads与参考基因组的比对效率在87.75%~87.80%,说明转录组测序数据真实可靠。采用SnpEff进行SNP/InDel变异注释分析,结果显示,SNP位点中以A>G、G>A、C>T和T>C等4种类型的数量较多（14950~21562个）,InDel位点则以SYNONYMOUS_CODING的数量最多（33630个）。使用StringTie对Mapped reads（比对到参考基因组的Reads）进行拼接,共发掘到4607个新基因,分别输入COG、GO、KEGG、KOG、Pfam、Swiss-Prot、eggNOG和nr数据库中进行序列比对,最终发现共有1098个新基因被注释,以被nr数据库注释的新基因数量最多（1077个）,而被COG数据库注释的新基因数量最少（416个）。基于Q-value<0.05且Fold Change>2的筛选条件,共获得1408个差异表达基因（661个为显著上调表达基因,747个为显著下调表达基因）;1408个差异表达基因被注释到53个GO功能条目中,其中,22条被注释到生物学过程（Biological process）,16条被注释到细胞组分（Cellular component）,15条被注释到分子功能（Molecular function）;KEGG信号通路富集分析发现3条重要的信号通路,分别是溶酶体通路（Lysosome）、氨基糖和核苷酸糖新陈代谢（Amino sugar and nucleotide sugar metabolism）及鞘脂类代谢（Sphingolipid metabolism）。在溶酶体通路中,CTSL基因、Nramp基因、MyoG基因和Myf5基因在凡纳滨对虾快速生长群体肌肉组织中呈上调表达,而Trypsin基因呈下调表达。【结论】通过转录组测序分析从快速生长群体和慢速生长群体的凡纳滨对虾肌肉组织中筛选出1408个差异表达基因（661个为显著上调表达基因,747个为显著下调表达基因）,主要富集在溶酶体、氨基糖和核苷酸糖新陈代谢及鞘脂类代谢等通路上,在对虾肌肉生长发育过程中发挥重要作用。     

意大利蜜蜂幼虫肠道在球囊菌侵染前期的 差异表达microRNA及其调控网络

【目的】 微小RNA（microRNA,miRNA）是一类重要的基因表达调控因子,可影响宿主与病原间的互作过程。蜜蜂球囊菌（Ascosphaera apis）是一种特异性侵染蜜蜂幼虫的致死性真菌病原。本研究旨在对意大利蜜蜂（Apis mellifera ligustica,简称意蜂）幼虫肠道在球囊菌胁迫前期的差异表达miRNA（differentially expressed miRNA,DEmiRNA）及其靶基因进行深入分析,在miRNA组学水平探究意蜂幼虫在球囊菌侵染前期的胁迫应答,并通过构建显著DEmiRNA的调控网络筛选出与宿主应答相关的关键miRNA。【方法】 利用small RNA-seq（sRNA-seq）技术对正常及球囊菌侵染的意蜂4日龄幼虫肠道（AmCK和AmT）进行高通量测序,首先对原始数据进行质控和评估,随后将过滤后的数据与西方蜜蜂（Apis mellifera）参考基因组进行比对;将比对上的序列标签（tags）注释到miRBase数据库,得出已知miRNA的表达量;通过TPM（tags per million）算法对各样本中miRNA的表达量进行归一化处理,以|log₂ fold change|≥1且P≤0.05作为标准筛选得到显著DEmiRNA;利用TargetFinder 软件预测显著DEmiRNA的靶基因,并对其进行GO和KEGG代谢通路（pathway）富集分析。利用Cytoscape软件对miRNA-mRNA调控网络进行可视化。最后,利用茎环反转录PCR（Stem-loop RT-PCR）和荧光定量PCR（qPCR）验证测序数据的可靠性。【结果】 AmCK和AmT样品的测序分别得到13 553 302和10 777 534条原始读段（raw reads）,经严格过滤后得到的有效读段（clean reads）数分别为13 186 921和10 480 913条。各样品的生物学重复间的Pearson相关性系数分别在0.9822和0.9508以上。共有10个显著DEmiRNA,包括4个上调miRNA和6个下调miRNA。显著DEmiRNA在AmT的整体表达水平低于AmCK。10个显著DEmiRNA可靶向结合3 788个靶基因,其中上调miRNA的1 240个靶基因可注释到GO数据库中的39个GO条目,主要富集在结合、细胞进程、代谢进程和应激反应等;下调miRNA的749个靶基因可注释到34个GO条目,主要富集在细胞进程、结合、代谢进程和应激反应等。KEGG数据库注释结果显示,上调miRNA和下调miRNA的靶基因分别注释到95和66条代谢通路,富集基因数最多的分别是Wnt信号通路、Hippo信号通路、光传导以及内吞作用、磷脂酰肌醇信号系统、嘌呤代谢。对于上调和下调miRNA,分别有31和52个靶基因注释到内吞作用,15和7个靶基因注释到泛素介导的蛋白水解,11和5个靶基因注释到Jak-STAT信号通路,1和3个靶基因注释到MAPK信号通路。显著DEmiRNA与靶mRNA之间形成复杂的调控网络,7个显著DEmiRNA靶向结合96个与Wnt信号通路相关的mRNA,8个显著DEmiRNA靶向结合55个与内吞作用相关的mRNA。Stem-loop RT-PCR和qPCR结果验证了测序数据的可靠性。【结论】 对意蜂幼虫肠道在球囊菌侵染前期的DEmiRNA及其靶基因进行预测和分析,并构建和分析了DEmiRNA-mRNA调控网络,研究结果提供了宿主miRNA的表达谱和差异表达信息,揭示了DEmiRNA通过调控细胞生命活动、新陈代谢以及部分细胞和体液免疫等生物学过程参与宿主的胁迫应答。miR-4331-y、miR-4968-y、miR-8440-y、novel-m0023-5p和novel-m0025-3p共同参与了宿主的Wnt信号通路和内吞作用的调控,可作为白垩病治疗的潜在分子靶标。     

基于Iso-Seq技术的野生马齿苋叶片转录组分析

基于全长转录组测序(isoform sequencing, Iso-Seq)技术对野生马齿苋叶片(对照、干旱和低温胁迫3种处理的叶片)转录组进行测序及其转录组信息分析的结果显示，总共得到641 732条原始片段，质控后得到合格的插入片段26 758 071条，总碱基量为44 Gb,利用CCS软件获得442 760条一致性序列，根据smrtlink软件得到全长的插入片段和全长的去除嵌合体的插入片段分别为391 258条和301 412条，去除冗余得到103 298条转录本，最终得到39 717条unigene,预测出40 952条LncRNA和38 419条CDS序列；并将unigene与Uniprot、Pfam、GO、KEGG、eggNOG、Pathway和Nr等数据库进行同源比对，其中36 381条unigene被注释，从这些unigene里查找到25 898个简单重复序列标记(simple sequence repeats, SSR)位点。利用差异表达分析软件DESeq2筛选不同处理组间的差异表达基因，对照和干旱处理、对照和低温胁迫处理及干旱和低温胁迫处理间分别鉴定出16 194(...     

蜜蜂球囊菌菌丝和孢子中长链非编码RNA的比较及潜在功能分析

[背景]蜜蜂球囊菌(Ascosphaera apis,简称球囊菌)是一种专性侵染蜜蜂幼虫的真菌病原,可引起成年蜜蜂数量和蜂群群势的急剧下降.长链非编码RNA(long non-coding RNA,lncRNA)是一类新近发现的非编码RNA,在表观遗传、细胞周期、剂量补偿等众多生命活动中发挥重要生物学功能.[目的]明确...     

2种不同生长规格墨瑞鳕肌肉组织转录组测序分析

【目的】掌握不同生长规格墨瑞鳕个体间差异表达基因的表达特点,为其功能相关基因深度挖掘及分子遗传育种提供科学依据。【方法】挑选同一养殖条件下极大个体和极小个体的墨瑞鳕,构建肌肉组织cDNA文库后,采用Illumina HiSeqTM 4000测序平台对存在生长差异的墨瑞鳕肌肉组织进行转录组测序分析,获得的Unigenes在Nr、Nt、Pfam、KOG/COG、Swiss-Prot、KEGG和GO等数据库中进行比对;通过FPKM及DEGseq筛选出差异表达基因,以GOseq和KOBAS对差异表达基因分别进行GO功能注释及KEGG信号通路富集分析,并采用MISA进行SSR鉴定分析。【结果】从墨瑞鳕肌肉组织中共测序获得39749条Unigenes,其长度范围在301~55230 bp,平均长度为1705 bp。注释到Nt、Nr、Swiss-Prot、Pfam数据库的Unigenes分别有27046、20824、18268和17772条,在7个数据库中均得到注释的Unigenes共计6742条,占Unigenes总数的16.96%。根据差异表达基因筛选条件P<0.05且|log₂ Fold Change|>1,共筛选出722个差异表达基因,其中上调基因308个、下调基因414个。差异表达基因GO功能注释分析结果表明,注释基因数目较多的GO功能条目包括细胞过程、代谢过程、膜、细胞器及结合等;KEGG信号通路富集分析发现,差异表达基因被成功富集到234条信号通路上,主要涉及磷脂酰肌醇3激酶/蛋白激酶信号通路、MAPK信号通路、胰岛素信号通路及FoxO信号通路等。在39749条Unigenes中鉴定筛选出22120个SSRs,占Unigene总数的55.65%,SSR的平均间距为3063 bp。【结论】基于转录组测序分析获得的墨瑞鳕肌肉组织差异表达基因以发挥结合、细胞过程及代谢过程等功能为主,且主要富集在PI3K-Akt信号通路、核糖体信号通路、FoxO信号通路及细胞凋亡等能量代谢相关通路上,通过共同协调而对墨瑞鳕的生长发育起调控作用。