张璐 作品数:14 被引量:69 H指数:5 供职机构: 福建农林大学蜂学学院 更多>> 发文基金: 国家现代农业产业技术体系建设项目 国家自然科学基金 福建省自然科学基金 更多>> 相关领域: 农业科学 医药卫生 生物学 更多>>
基于RNA-seq数据大规模挖掘蜜蜂球囊菌的SSR分子标记 被引量:5 2017年 基于已获得的球囊菌转录组数据预测微卫星标记(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引物是可行的. 李汶东 熊翠玲 王鸿权 侯志贤 童新宇 张璐 付中民 郑燕珍 陈大福 郭睿关键词:蜜蜂球囊菌 微卫星标记 转录组 RNA-SEQ 西方蜜蜂的基因注释信息优化及新基因鉴定 被引量:5 2018年 西方蜜蜂(Apis mellifera)全基因组序列早在2006年就已公布,但其基因注释信息并不完善。本研究利用前期获得的意蜂幼虫肠道RNA-seq数据对西方蜜蜂基因组的已知基因进行结构优化,对新基因进行预测、鉴定和分析。利用Cufflink软件进行转录本重构,进而与参考转录本序列进行比较,共对4 177个已公布西方蜜蜂基因的测序结果进行了优化,其中5′端延长、3′端延长、5′和3′端都延长的基因分别有2 383、2 392和2 384个。利用Cuffcompare软件将重构转录本与参考基因组进行比对,共预测出197个西方蜜蜂新基因,随机挑选20个新基因进行RT-PCR验证,有18对引物成功扩增出符合预期的目的片段。GO分类结果显示这些新基因涉及运输和代谢、能量代谢和信号分子和交互等21个GO terms。KEGG代谢通路(pathway)富集分析结果显示这些新基因富集于13个pathways,其中基因富集数最多的是细胞进程、单一有机体和结合。本研究为西方蜜蜂基因组信息提供了有益的补充和完善,也为新基因的功能研究打下了基础。 郭睿 张璐 熊翠玲 郑燕珍 付中民 陈大福关键词:西方蜜蜂 意大利蜜蜂 基于RNA-seq数据大规模开发中华蜜蜂幼虫的SSR分子标记 被引量:17 2017年 中华蜜蜂是东方蜜蜂的一个亚种,也是我国的特有蜂种。本研究基于已获得的中华蜜蜂幼虫肠道转录组数据预测SSR分子标记,并进行SSR位点的信息分析和SSR引物的发掘。利用MISA软件对幼虫肠道转录组中数据组装得到的43557条unigenes进行搜索,共预测出13448个SSR位点,它们分布于7763条unigene中,其中最主要的重复类型为二核苷酸重复(58.03%),其次为三核苷酸重复(28.23%)和四核苷酸重复(9.72%)。二核苷酸重复中的基序主要是AT/AT(占总量的30.4%)。对于所有的SSR位点,利用Primer Premier 5软件成功设计出21627对引物,随机选取48对引物对5个不同来源的中华蜜蜂幼虫肠道样品进行SSR位点扩增,共有15对成功扩增出符合预期的目的片段。研究结果表明利用转录组数据大规模开发中华蜜蜂幼虫的SSR引物是可行的,本研究开发出的SSR引物为研究中华蜜蜂分子遗传学奠定了基础。 熊翠玲 张璐 付中民 王鸿权 侯志贤 童新宇 李汶东 郑燕珍 陈大福 郭睿关键词:RNA-SEQ 中华蜜蜂 SSR分子标记 胁迫意大利蜜蜂幼虫肠道的球囊菌的转录组分析 被引量:26 2017年 【目的】本研究旨在通过趋势分析对胁迫意大利蜜蜂Apis mellifera ligustica(简称意蜂)幼虫肠道的球囊菌Ascosphaera apis的差异表达基因(DEGs)进行转录组分析。【方法】将纯化的球囊菌孢子配制为1×107孢子/m L的饲料饲喂意蜂3日龄幼虫,利用Illumina Hi Seq 2500平台对球囊菌胁迫的意蜂幼虫肠道c DNA进行测序,将过滤得到的有效读段(clean reads)映射(mapping)到核糖体数据库及意蜂参考基因组,最后将未映射上的reads映射到本课题组组装并注释的球囊菌参考转录组。利用STEM软件对DEGs进行趋势分析。利用WEGO软件对显著表达模式中的DEGs进行GO富集分析。利用Blastall对显著表达模式中的DEGs进行KEGG代谢通路富集分析。最后,通过对随机选取的6个DEGs进行RT-q PCR分析,对RNA-seq数据进行验证。【结果】球囊菌胁迫意蜂幼虫肠道样品的Illumina测序共得到球囊菌的25 454 076条原始读段(raw reads),经过滤得到24 909 820条clean reads,Q30均在93.46%以上。趋势分析结果显示,19 893个DEGs聚类为8个表达模式,其中有12 151个DEGs聚类为3个表现为显著上调趋势的表达模式。GO富集分析结果显示,表现上调趋势的DEGs富集在40个GO term,富集基因数最多的为细胞进程(cellular process)(2 601 unigenes),其次为代谢进程(metabolic process)(2 553 unigenes)和细胞(cell)(2 522unigenes)。KEGG代谢通路富集分析结果显示,上调趋势中的DEGs富集在119个代谢通路上,其中富集基因数最多的是核糖体(ribosome)(213 unigenes),其次为氨基酸生物合成(biosynthesis of amino acids)(154 unigenes)和内质网蛋白加工(protein processing in endoplasmic reticulum)(130unigenes)。共有48个DEGs富集在MAPK信号通路上,聚类分析结果显示,这些DEGs随着胁迫时间的延长表达水平逐渐增强。RT-q PCR结果显示,6个DEGs的表达水平变化趋势与RNA-seq数据一致,证明了本研究中的转录组数据真实可靠。【结论】本研� 陈大福 郭睿 熊翠玲 梁勤 郑燕珍 徐细建 黄枳腱 张曌楠 张璐 李汶东 童新宇 席伟军关键词:意大利蜜蜂 RNA-SEQ 差异表达基因 中华蜜蜂幼虫肠道参考转录组的de novo组装及其响应球囊菌胁迫的免疫应答研究 蜜蜂易受到多种疾病侵袭,其中,白垩病是最具代表性的真菌病,给养蜂业造成巨大损失.目前,尚无一种杀真菌剂被批准用于治疗白垩病.养蜂生产中,意大利蜜蜂幼虫极易受到蜜蜂球囊菌的侵染而罹患白垩病,而中华蜜蜂幼虫具有较强的球囊菌抗... 陈大福 郭睿 熊翠玲 徐细建 Dhiraj Kumar 骆群 郑燕珍 张曌南 张璐 李汶东 席伟军 李龙 杜宇 张琦 余敏 马子韬 侯志贤 曹雪珂 王博 梁勤关键词:中华蜜蜂 SSR 免疫应答 基于比较转录组学分析揭示中华蜜蜂及意大利蜜蜂幼虫的球囊菌抗性差异机制 被引量:4 2019年 为探究中华蜜蜂(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结果显示中蜂的单拷贝基因仅富集在氧化磷酸化。上述结果表明免疫相关基因数量的差异是二者的球囊菌抗性差异产生的重要原因之一,中蜂在与球囊菌的协同进化过程中可能通过提高能量利用从而限制球囊菌的增殖。本研究结果可为阐明中蜂及意蜂幼虫的球囊菌抗性差异产生的分子机制提供参考。 熊翠玲 杜宇 王鸿权 郑燕珍 付中民 王海朋 张璐 陈大福 郭睿关键词:中华蜜蜂 意大利蜜蜂 幼虫 意大利蜜蜂幼虫肠道响应球囊菌胁迫的可变剪接基因分析 被引量:4 2018年 为探讨可变剪接基因在宿主响应球囊菌(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运输和碳代谢。上述分析结果揭示了可变剪接基因在意蜂幼虫肠道响应球囊菌胁迫过程中的重要作用,为深入研究可变剪接基因的功能打下了基础。 郭睿 陈恒 张璐 赵浩宇 熊翠玲 郑燕珍 付中民 梁勤 陈大福关键词:意大利蜜蜂 胁迫 可变剪接 意大利蜜蜂幼虫肠道响应球囊菌胁迫的可变剪接基因分析 本研究基于前期已获得的意大利蜜蜂幼虫肠道转录组数据进行可变剪接基因分析,旨在揭示可变剪接基因在宿主响应胁迫过程中的作用.本研究中,在正常(AmCK)及球囊菌胁迫的意蜂幼虫肠道样品(AmT1、AmT2、AmT3)中共鉴定出... 郭睿 陈恒 张璐 赵浩宇 熊翠玲 郑燕珍 付中民 梁勤 陈大福关键词:意大利蜜蜂 胁迫 可变剪接 中华蜜蜂幼虫肠道参考转录组的de novo组装及SSR分子标记鉴定 被引量:13 2017年 【目的】利用RNA seq技术对中华蜜蜂(Apis cerana cerana,简称中蜂)幼虫肠道参考转录组进行de novo组装,并进行功能及代谢通路注释,进而利用该转录组数据进行中蜂幼虫的SSR分子标记鉴定。【方法】实验室条件下饲养中蜂幼虫,将纯化的蜜蜂球囊菌(Ascosphaera apis,简称球囊菌)孢子饲喂3日龄幼虫,剖取4、5和6日龄幼虫肠道,液氮速冻。将健康幼虫肠道与感染球囊菌的幼虫肠道同时进行Illumina测序。通过对raw reads的过滤得到clean reads,利用Trinity软件组装得到unigenes。通过BLASTx(E-value<10-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分子标记。【结论】成功组装并注释了中蜂幼虫肠道参考转录组,可为中蜂及其幼虫的分子生物学及组学 徐细建 郭睿 骆群 熊翠玲 梁勤 张串联 郑燕珍 张曌楠 黄枳腱 张璐 李汶东 陈大福关键词:中华蜜蜂 UNIGENE SSR 中华蜜蜂幼虫肠道的高表达基因分析 2017年 中华蜜蜂(中蜂)是中国特有的本土蜜蜂资源,也是养蜂生产中的常用蜂种之一,具有重要的生态和经济价值。本研究利用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个。研究结果不仅可提供中蜂幼虫肠道发育过程中的基因表达谱数据,也为分子水平上深入研究中蜂幼虫肠道发育提供重要信息。 解彦玲 王鸿权 张璐 侯志贤 刀晨 江亮亮 熊翠玲 郑燕珍 徐细建 黄枳腱 郭睿 陈大福关键词:RNA-SEQ 中华蜜蜂 高表达基因