大垫尖翅蝗转录组分析

【目的】大垫尖翅蝗Epacromius coerulipes(Ivanov)是广泛分布的草原蝗虫之一,但其基因资源缺乏。为了获得大垫尖翅蝗的基因数据,对其进行了转录组测序和分析。【方法】利用Illumina公司paired-end转录组的测序技术进行从头组装。【结果】总计获得了63 033条unigenes,平均长度为772 bp,N50为1 589 bp。通过BLAST搜索,确定有25 132条(39.87%)unigenes与NCBI数据库已知的蛋白质相匹配,其中有24 841,16 490,11 558和8 013条unigenes成功注释到Nr,Swiss-Prot,GO和COG数据库中。KEGG数据库中,7 218条unigenes形成218条代谢或信息通路。其中,189条unigenes参与外源性物质或药物的代谢通路。进一步分析显示,213条unigenes被确认为可能参与外源性物质的解毒作用,29条unigenes被确定为编码杀虫剂的目标蛋白。此外,检测到5 696条简单重复序列。【结论】该转录组测序分析将为进一步研究大垫尖翅蝗的基因功能分析及杀虫剂的抗药性机制分析奠定分子基础。     

意大利蝗转录组及性腺发育相关基因分析

【目的】意大利蝗Calliptamus italicus是新疆荒漠、半荒漠草原优势危害种类之一,其生殖相关基因信息缺乏。本研究对意大利蝗精巢和卵巢组织进行高通量转录组测序,旨在揭示意大利蝗转录组特征,并获得与性腺发育相关的基因数据。【方法】利用Illumina高通量测序平台(Hi Seq/Mi Seq)对意大利蝗进行了转录组测序,并进行序列分析。【结果】获得718 872条unigenes,平均长度631 bp,N50为1 186 bp。通过同源性比对,成功注释254 597条unigenes,其中,Nr数据库注释占28.87%,比例最高。Nr数据库注释的unigenes与黑褐草蚁Lasius niger的相似基因序列最多,为10.80%。与GO数据库比对发现,意大利蝗转录组unigene可分为3大类:涉及分子功能的unigenes有31 392条,涉及细胞组分的unigenes有20 586条,与生物学过程有关的unigenes有57 014条。KEGG pathways分析表明,23 666条unigenes归属于251个通路。涉及生殖系统发育的通路有卵母细胞减数分裂、胰岛素信号通路、Wnt信号通路、促性腺激素释放激素信号通路、转化生长因子β信号通路以及泛素介导的蛋白水解、黄体酮介导的卵母细胞成熟、昆虫激素生物合成等。进一步筛选得到部分与性腺发育相关的基因,包括vg,vgr,sox 21b,testis-specific serine/threonine-protein kinase,spermatogenesis-associated protein和vasa-like gene等。【结论】本研究获得了意大利蝗转录组数据及性腺发育相关基因,为进一步研究意大利蝗生殖调控机理奠定了分子基础。     

疣蝗转录组分析

疣蝗分布广泛,是一种常见的害虫,主要危害农作物,但分子信息缺乏。我们采用高通量Illumina测序技术和从头组装的方法获得了疣蝗转录组unigenes,并鉴定了天然免疫相关信号通路基因。本研究共获得疣蝗转录组unigenes 86 756个,平均长度为670 bp,N50为1 033 bp。使用六大数据库Nr、Nt、Swiss-Prot、GO、COG和KEGG注释疣蝗转录组unigenes,分别对应有29 627,13 866,23 398,13 474,11 968和21 027个unigenes获得注释。本研究共鉴定319个天然免疫信号通路相关的unigenes,其中,参与Toll通路的unigenes 49个,参与免疫缺失(IMD)通路的unigenes 13个,参与JAK-STAT信号通路的unigenes 117个,参与MAPK-fly通路相关的unigenes 58个,参与黑化(Melanization)通路的unigenes 82个。另外,我们还预测到了9 287个微卫星(SSRs)。本研究为疣蝗基因功能分析、免疫通路相关基因和分子标记的进一步研究提供了分子理论基础。     

高丛蓝莓品种“瑞蓝”的幼叶转录组分析

蓝莓营养丰富,富含花青素等抗氧化性活性物质,被誉为"小浆果之王"。为了加深对蓝莓优良品种"瑞蓝"生长发育分子机制的了解,利用第二代测序技术对其2年生苗的幼叶转录组进行测序。获得总共约1 G的有效测序数据,共拼装成41 987条unigenes,其中23 555条unigenes能被GO数据库成功注释,分别涉及到生物学过程、细胞成分及分子功能相关的共35种生理代谢功能;20 379条unigenes能被KOG数据库注视到,共涉及25条代谢通路;9 771条unigenes能被KEGG数据库成功注释,共涉及到5个大类、20个中类、124条代谢通路;27 293条unigenes被Nr数据库成功注释,这些unigenes共对应32 949条蛋白质的编码区。同时,生物信息学分析还显示:全部unigenes中,有5 221条unigenes编码转录因子,涉及到54个家族;2 487条unigenes编码抗性基因,涉及19个家族;共检测到6 261个SSR序列,其中2碱基重复的SSR占68.2%。转录组测序分析结果阐明了"瑞蓝"的生长发育的部分分子信息,对利用其进一步开展分子标记辅助育种等具有较为重要的意义。     

中华大仰蝽转录组学研究及SSR新标记开发

【目的】中华大仰蝽Notonecta chinensis为中国和日本冲绳分布的重要水生天敌昆虫,可用于蚊虫的生物防治。本研究旨在建立中华大仰蝽转录组数据库,挖掘其基因信息。【方法】采用高通量测序平台Illumina NextSeq500对中华大仰蝽进行转录组测序、de novo组装及生物信息学分析;利用MISA软件基于转录组unigenes数据进行SSR新分子标记筛选。毛细管电泳检测SSR多态性。【结果】总计获得34782282条clean reads(NCBI SRA数据库登录号:SRR13259254),组装成37801条unigenes,N50为913 bp。将unigenes与已知数据库比对进行基因功能注释,分别有36474,32470,27781,35079和5638条序列注释到nr,Swiss-Prot,GO,eggNOG和KEGG数据库。通过GO数据库注释,unigenes的功能可分为生物学过程、细胞组分和分子功能三大类,其中参与细胞、细胞部分及结合功能的unigenes比例较大。eggNOG数据库注释结果显示,37801条unigenes归到25个基因家族,注释到未知功能的最多。KEGG代谢通路富集分析显示,5638条unigenes注释到245个代谢通路,注释到核糖体的数目最多。此外,用MISA软件在转录组测序数据中的37801条unigenes中搜索到3124个SSR位点(占总unigenes的8.26%),发生频率为7.07%。通过PCR筛选出16个SSR位点。7个中华大仰蝽地理种群3个位点NcCF/NcCR,NcKF/NcKR和NcLF/NcLR的多态信息含量(PIC)分别为0.870,0.902和0.857,具高度多态性。【结论】本研究成功获得了中华大仰蝽转录组数据,为其基因功能分析提供了分子理论基础;SSR新标记的开发为中华大仰蝽遗传多样性分析、隐存种鉴定及基因图谱构建提供了更丰富的候选分子标记。     

泽兰实蝇雌雄成虫转录组及差异分析

[目的]泽兰实蝇Procecidochares utilis Stone是入侵杂草紫茎泽兰Eupatorium adenophorum Spreng的重要的专食性天敌,已成为控制紫茎泽兰的重要因子.本研究旨在获得泽兰实蝇转录组,并深入研究泽兰实蝇雌雄成虫差异表达基因.[方法]利用Illumina高通量测序技术对泽兰实蝇雌雄成虫进行转录组测序和生物信息分析.[结果]总共获得29 147条unigenes,平均长度为457 bp,同时将所得序列注释到七大数据库进行比对,共获得19 384条unigenes注释结果.分析发现11 331条差异表达基因;与雄成虫相比,雌成虫有2 640条unigenes表达量上调,8 691条unigenes表达量下调.[结论]本研究获得了泽兰实蝇转录组,为今后泽兰实蝇性别相关研究提供了序列支持.     

绿豆象触角转录组及嗅觉相关基因的分析

【目的】建立绿豆象Callosobruchus chinensis触角转录组数据库,挖掘绿豆象的基因数据。【方法】采用高通量测序平台(Illumina Hiseq)对绿豆象成虫触角进行转录组测序、序列组装及生物信息学分析。【结果】共获得51.10 Gb的clean data(NCBI SRA数据库登录号:SRP119884)及83 535条unigenes,长度在1 kb以上的unigenes有15 075条,unigenes的N50为1 492 bp。使用BLAST软件将绿豆象触角Unigene序列与公共数据库比对,共注释到22 148条unigenes,其中NR数据库注释的unigenes最多,为18 744条,且与赤拟谷盗Tribolium castaneum的相似基因最多,达39.57%。通过GO数据库注释,将unigenes功能分为细胞组分、分子功能与生物学过程三大类54个分支,其中参与代谢进程以及结合与催化活性的unigenes比例较大。KEGG代谢途径分析表明,9 084条unigenes形成289条代谢通路,其中核糖体代谢通路所含unigene最多,为516条。进一步基因注释分析筛选到140个嗅觉相关基因,包括12个气味结合蛋白(odorant binding protein,OBP)基因,4个化学感受蛋白(chemosensory protein,CSP)基因,116个气味受体(odorant receptor,Or)基因,1个离子型受体(ionotropic receptor,IR)基因和7个感觉神经元膜蛋白(sensory neuron membrane protein,SNMP)基因,并用FPKM值对基因表达量进行评估。【结论】本研究获得了绿豆象触角转录组数据,为进一步研究绿豆象的基因功能分析及嗅觉感受机制奠定分子基础。     

梨小食心虫幼虫中肠转录组及SSR分子标记分析

【目的】梨小食心虫Grapholitha molesta是一种重要的世界性果树害虫。昆虫中肠在食物消化、免疫反应、生长发育等方面发挥着重要作用。本研究旨在建立梨小食心虫幼虫中肠转录组数据库,挖掘其基因信息。【方法】利用高通量测序平台(Illumina HiSeq X Ten)对梨小食心虫幼虫中肠进行转录组测序、组装及生物信息学分析;进而利用转录组数据进行梨小食心虫幼虫中肠SSR分子标记鉴定。【结果】总计获得梨小食心虫幼虫中肠转录组96 419条unigenes,与公共数据库比对,共注释到57 300条unigenes。通过GO数据库注释,将unigenes的功能分为生物学进程、细胞组分和分子功能三大类55个功能区,其中参与细胞进程、细胞和细胞组分及结合功能的unigenes比例较大。KOG结果显示,10 090条unigenes归到25个基因家族,注释到一般功能预测的最多。KEGG数据库中,10 250条序列注释到232个代谢通路,注释到核糖体的基因数最多。此外,用MISH软件,在10 600条unigenes中搜索到12 690个SSR位点,通过PCR开发出6个SSR位点。【结论】本研究成功组装梨小食心虫幼虫中肠参考转录组,为以中肠为靶标的害虫防治提供理论依据。     

松树蜂毒腺转录组分析

【目的】构建入侵种松树蜂Sirex noctilio毒腺转录组数据库,筛选并分析松树蜂毒腺基因数据。【方法】采用新一代高通量测序平台Illumina HiSeq^TM 4000对松树蜂雌成虫毒腺进行转录组测序、数据组装和生物信息学分析。【结果】共获得12.7 Gb松树蜂雌成虫毒腺有效转录组数据,并组装到37 098条unigenes,平均长度968 bp,N50长度为2 364 bp。将所得的unigenes数据使用BlastX与各大数据库比对,共注释到13 515条unigenes,并且在NR数据库中注释的unigenes最多,共11 108条(占总数的29.94%),其中相似基因占比最高的物种为丽蝇蛹集金小蜂Nasonia vitripennis,达815条(占总数的7.29%)。在GO数据库中注释到5 726条unigenes,根据功能被分为生物学进程、细胞组分和分子功能3大类63个亚类。KEGG代谢通路分析表明,7 602条unigenes注释到357个代谢通路。根据基因注释信息进一步筛选到43条嗅觉相关基因,包括嗅觉受体(odorant receptor, Or)基因25条、化学感受蛋白(chemosensory protein, CSP)基因10条、离子型受体(ionotropic receptor, IR)基因5条和气味结合蛋白(odorant binding protein, OBP)基因3条。此外,还筛选出11条漆酶基因,包括漆酶1(laccase1, LAC1)基因5条、漆酶2(laccase2, LAC2)基因4条、漆酶4(laccase4, LAC4)基因1条和漆酶9(laccase9, LAC9)基因1条,且其中1条LAC2基因在所有被注释的基因中表达量最高(FPKM值=21 126)。【结论】本研究获得的松树蜂毒腺转录组数据为松树蜂毒液组分的鉴定和生物学功能的研究奠定了一定的理论基础。     

叉角厉蝽触角转录组及嗅觉相关基因分析

叉角厉蝽Eocanthecona furcellata是一种在生物防治方面有重要作用和潜力的捕食性天敌昆虫,触角是昆虫进行信息交换的重要器官,而嗅觉相关基因则是调控天敌昆虫捕食行为的重要分子基础。为获得叉角厉蝽触角转录组数据库,挖掘叉角厉蝽嗅觉相关基因,本研究利用Illumina高通量测序平台对叉角厉蝽雌雄成虫与5龄若虫触角进行转录组测序。成功构建了叉角厉蝽触角转录组,获得了67 843条unigenes, N50长度为2 300 bp。与七大公共数据库比对注释到27 686条unigenes,其中NR数据库注释最多(33.33%),且与茶翅蝽Halyomorpha halys相似度最高(64.20%)。14 258条注释到GO数据库中,分为生物过程、细胞组分和分子功能3个大类42个亚类;KEGG代谢途径分析表明,7 703条形成282条代谢通路,其中被注释在信号传导通路中的unigenes最多(11.50%)。进一步基因注释分析,鉴定得到134个候选嗅觉相关基因,32个化学感应蛋白基因(Chemosensory protein genes, CSP),10个气味结合蛋白基因(Odor...     

基于高通量测序的玉米象转录组分析

玉米象Sitophilus zeamais是一种世界性的储粮害虫,但其基因信息尚不完善.本研究利用高通量测序平台Illumina HiSeq TM2000对玉米象成虫进行转录组测序,总计获得64358条unigenes,总长度39481752 bp,最短201 bp,最长29046 bp,平均长度613 bp.将unigenes序列在7大数据库Nr、GO、Swissprot、KOG、Pfam、KEGG、TrEMBL中进行注释,分别有25271、20747、16477、12060、9471、3370、25500条unigenes获得注释.通过GO功能分类,共有20747条unigenes在GO数据库中生物学过程、细胞组分、分子功能3大类68个亚类功能组中找到对应.KOG结果显示,12606条unigenes归到25个基因家族.通过KEGG代谢通路分析,共有3370条unigenes被注释,分别归属于细胞进程、环境信息进程、遗传信息进程、新陈代谢和有机体系统5大类代谢途径,共33个亚类274个功能通路.进一步的数据分析,共鉴定出146081个SNP位点和5002个简单序列重复(SSR)位点.本研究获得的玉米象转录组信息,为玉米象的功能基因挖掘提供了重要的信息资源.     

锤胁跷蝽触角转录组分析

为了获得锤胁跷蝽(Yemma signatus Hsiao)触角基因信息,本研究使用Illumina HiSeq TM 4 000高通量测序技术对锤胁跷蝽成虫触角进行转录组测序和生物学信息分析。共获得61 080 938条clean reads,包括9.16 G (GenBank登录号:SRR3348966)。拼接115 491条unigenes,平均长度为578 bp,N50为863 bp。在7大数据库中注释到46 224条unigenes,其中在NR数据库中注释的基因最多,为32 842 (28.43%)条,与内华达古白蚁(Zootermopsis nevadensis)相似度最高(9.0%)。转录组数据分析鉴定出125个与嗅觉相关的基因,其中66个嗅觉受体、11个味觉受体、1个离子型受体、3个感觉神经元膜蛋白、30个气味结合蛋白和14个化学感受蛋白基因。本研究首次获得锤胁跷蝽触角转录组数据,并鉴定出嗅觉相关基因,本研究为今后利用嗅觉基因靶标防治害虫提供了分子生物学信息数据。     

花椒窄吉丁转录组及化学感受相关基因的分析

【目的】建立花椒窄吉丁 Agrilus zanthoxylumi 转录组数据库,挖掘其基因组数据。【方法】采用高通量测序平台Illumina NovaSeq6000对花椒窄吉丁成虫不同组织(触角、头、胸、腹、足、翅)进行转录组测序、序列组装,使用BLAST软件将花椒窄吉丁unigenes序列与公共数据库进行比对,BLAST同源性搜索的方法从中筛选出与其化学感受相关的基因,并与其他已研究发表的鞘翅目昆虫化学感受相关基因的核酸序列比对,进行系统发育分析。利用RPKM值对这些基因在花椒窄吉丁成虫不同组织中的表达量进行分析。【结果】经测序及序列拼接后共得到80 320条unigenes和 169 398 条contigs,其G+C比例分别是37.63%和39.18%,平均长度分别为828.75和1 084.33 bp;unigenes长度主要分布在200~400 bp的共有37 374条,占全部unigenes的46.53%。在NR数据库中注释的与花椒窄吉丁转录组相似基因所属物种分布为赤拟谷盗 Tribolium castaneum 所占比例26.80%,其后依次是小家鼠 Mus musculus (13.87%),蛀犀金龟 Oryctes borbonicus (5.26%),山松甲虫 Dendroctonus ponderosae (4.10%),黑蚁 Lasius niger (2.36%),其他物种所占比例是47.61%。在GO数据库中比对到27 488个unigenes,可分为分子功能、细胞组分和生物学进程三大类共59个过程;利用KOG数据库功能注释分为25类,注释到unigenes 最多功能类别的是普通功能,共 4 666个,最少的是与细胞移动相关,共40个;将80 320个unigenes映射到KEGG数据库中,21 104条unigens注释到代谢通路35个,占26.27%,注释最多的代谢途径是转录途径,共涉及到2 037条unigenes。从NR数据库中共注释到8个气味受体(odorant receptor, OR)基因、7个离子型受体(ionotropic receptor, IR)基因和6个味觉受体(gustatory receptor, GR)基因,各组织中的基因表达量分析发现,这3类化学感受相关基因在成虫触角中的表达量显著高于其他组织中的表达量,雌、雄触角中表达量差异最显著的为 AzanOR 1,而在雌、雄成虫足中这3类化学感受相关基因表达量没有差异。【结论】本研究首次获得了花椒窄吉丁转录组数据,为进一步研究花椒窄吉丁的基因功能奠定了分子基础。     

李白盾蚧转录组分析及SSR位点开发

李白盾蚧Pseudaulacaspis prunicola (Maskell)寄主范围广泛,是一种重要的入侵害虫。本研究利用高通量测序平台（Illumina NovaSeq 6000）对李白盾蚧进行转录组测序、de novo从头组装及功能注释,在此基础上筛选其微卫星（SSR）位点,并挖掘微卫星引物。研究共获得李白盾蚧转录组60 296条转录本,24 967条单基因（unigenes）序列。通过GO数据库注释,将所有unigenes的功能分为生物学进程、细胞组分和分子功能三大类41个亚类功能区。KOG数据库注释结果显示,5 085条unigenes归到25个基因家族,注释到一般功能预测的数目最多。KEGG代谢通路富集分析显示6 668条unigenes注释到280个代谢通路,其中注释到内质网中蛋白质加工的数目最多。利用MISA软件共搜索到微卫星位点18 193个,分布在9 043条unigenes中,占总unigenes数量的36.22%,平均每2.29 kb出现一个SSR位点。其中主要重复类型为单核苷酸重复,占SSR位点总数的72.03%,其次为三核苷酸重复（15.90%）和二核苷酸重复（8.48%）。单核苷酸重复主要为A/T（71.16%）,二核苷酸重复主要为AG/CT（5.20%）。基于Primer Primer 3软件设计出12 538对李白盾蚧SSR引物,从中随机挑选50对引物进行PCR验证,共29对引物可以稳定扩增出目的片段。本研究成功组装了李白盾蚧转录组数据,并基于转录组数据成功筛选出其微卫星位点,为未来该虫的种群遗传学以及入侵生物学研究提供了数据支撑。     

大蜡螟触角转录组及嗅觉相关基因分析

【目的】大蜡螟Galleria mellonella Linnaeus是世界范围内蜂群中普遍存在的害虫,给世界养蜂业带来了巨大的危害。本研究旨在建立大蜡螟触角转录组数据库,挖掘大蜡螟嗅觉相关基因,为今后利用嗅觉基因靶标防治大蜡螟提供分子生物学信息数据。【方法】利用高通量测序平台(Illumina HiSeqTM 2500)对大蜡螟触角进行转录组测序和组装,通过筛选转录组数据库,鉴定出嗅觉相关基因。【结果】成功构建了大蜡螟触角转录组,经序列拼接获得188 278条unigenes,平均长度为781 bp,N50为1 161 bp。使用BLAST软件将大蜡螟触角unigenes序列与7大公共数据库比对,共注释到108 047条unigenes,其中NR数据库注释的unigenes最多,为78 746条,占41.82%,且NR注释的大蜡螟unigenes中与家蚕Bombyx mori类似基因最多,为18.4%。将unigenes与GO数据库比对发现,共有69 794条unigenes注释到GO数据库中,按照其所参与的生物过程、细胞组分和分子功能对其进行GO功能分类,共含有55个亚类,其中参与结合和催化功能及代谢过程的unigenes比例较大;KEGG代谢途径分析表明,14699条unigenes形成231条代谢通路,其中被注释在信号传导通路中的unigenes最多,为2 401条,占16.33%。进一步基因注释分析,鉴定得到226个候选嗅觉相关基因,包括52个气味结合蛋白(Odorant binding proteins,OBPs)基因,36个化学感受蛋白(Chemosensory proteins,CSPs)基因,80个气味受体(Odorant receptors,ORs)基因,56个离子型受体(Ionotropic receptors,IRs)基因和2个感觉神经元膜蛋白(Sensory neuron membrane proteins,SNMPs)基因;通过比较雌、雄大蜡螟转录组鉴定出114个差异表达基因,包括5个嗅觉相关基因(OBP8、OBP17、CSP3、CSP34和OR15);114个差异表达基因中,66个基因在雌蛾触角中高表达,48个基因在雄蛾触角中高表达。【结论】本研究获得了大蜡螟触角转录组数据,并鉴定出候选嗅觉相关基因,结果为进一步研究大蜡螟的基因功能及嗅觉感受机制奠定分子基础。     

扶桑绵粉蚧转录组分析

【目的】建立扶桑绵粉蚧Phenacoccus solenopsis Tinsley的转录组数据库,揭示扶桑绵粉蚧转录组的整体表达特征。【方法】采用Illumina Hi SeqTM40`00测序平台开展对扶桑绵粉蚧雌成虫转录组测定,对原始数据进行过滤和组装。【结果】共获得58 322 258条序列读取片段(reads),共9.48 Gb(Gen Bank登录号:SAMN06130426)57 422 032条有效转录组数据。进一步组装拼接后,共获得94 475个单基因簇(unigene),平均长度为700 bp。将unigene与数据库中的序列进行BLASTX比对,成功注释20 949个unigenes,其中,Nr注释的unigenes与豌豆蚜Acyrthosiphon pisum unigenes同源性最高,达18.69%。扶桑绵粉蚧转录组unigenes根据GO功能注释大致可分为细胞组分、分子功能和生物过程三大类55个分支,与结合活性、催化活性、细胞进程和代谢进程相关的unigenes较多。此外,本研究还筛选到20条与脂类代谢相关的途径和与性信息素代谢相关的序列。并通过与Nr和Swiss-prot蛋白质数据库比对,获得了15 037条编码序列(CDS)片段。【结论】本研究初步阐明扶桑绵粉蚧雌成虫转录组的整体表达模式,为进一步研究扶桑绵粉蚧的基因功能及性信息素的代谢途径奠定了基础。     

蜜蜂球囊菌的参考转录组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分子标记可推动菌株鉴定、基因图谱构建及基因定位等研究。     

荔枝蒂蛀虫转录组及嗅觉相关基因分析

【目的】荔枝蒂蛀虫Conopomorpha sinensis Bradley是专一性危害我国华南地区荔枝和龙眼的重要害虫,隐蔽性强,防治困难,基因组信息缺乏。本研究的目的是获得荔枝蒂蛀虫的基因数据,寻求有效控制害虫的分子靶标。【方法】采用新一代高通量测序技术Illumina Hi SeqTM4000对荔枝蒂蛀虫进行转录组测序和生物信息学分析。【结果】经序列拼接获得68 996条unigenes。进一步利用七大公共数据库进行同源比对,注释了22 348 unigenes。注释到Nr数据库的unigenes数量最多,达27.01%,其中Nr注释的荔枝蒂蛀虫unigenes中与小菜蛾Plutella xylostella unigenes同源性最高,达34.1%。将unigenes与GO数据库比对发现,15 585条unigenes根据其功能大致可分为3类47亚类。KEGG pathways分析表明,7 272条unigenes定位为267个代谢通路。基因注释进一步筛选鉴定获得100个荔枝蒂蛀虫嗅觉相关基因;与鳞翅目相关气味结合蛋白基因联合分析发现,与荔枝蒂蛀虫气味结合蛋白基因直系同源的有18组,部分基因形成独立一簇。【结论】本研究首次获得了荔枝蒂蛀虫的转录组数据,研究结果为生物控制荔枝蒂蛀虫提供了重要的基础数据和候选分子靶标。荔枝蒂蛀虫独有的气味结合蛋白基因可能与其生境中特有的化学物质相关。     

黑须污蝇不同发育阶段转录组及嗅觉相关基因分析

【目的】建立黑须污蝇Wohlfahrtia magnifica转录组数据库,挖掘黑须污蝇基因数据,为更深入地研究双峰驼阴道蝇蛆病提供理论依据。【方法】采用高通量测序平台Illumina HiseqTM4000对黑须污蝇幼虫、蛹和成虫进行转录组测序,并进行生物信息学分析。利用黑须污蝇幼虫、蛹和成虫转录组数据,通过基因差异表达分析以及GO功能显著性富集分析的方法筛选出嗅觉相关基因,并通过荧光定量PCR技术对黑须污蝇幼虫、蛹和成虫中OBP99b,OBP56a和OBP99a 3个嗅觉相关基因表达水平进行验证。【结果】结果显示,每个黑须污蝇样本的转录组数据量在4.96 Gb以上,G+C含量在35.35%以上,Q20含量在97%以上。与NCBInr, GO, KEGG, Pfam, Swiss-Prot和eggNOG数据库进行对比,共注释到73 303条unigenes。其中eggNOG数据库注释到的unigenes最多,为35 066条,分布在23个蛋白功能中,其中的翻译修饰、蛋白质周转的unigenes所占比例较大;GO数据库注释到的unigenes数为29 193条,功能包括分子功能、细胞组分、生物学过程三大类50个分支,其中参与生物学过程和氧化还原过程的unigenes比例较大;KEGG数据库注释到的unigenes数为15 068条,其中参与信号转导过程的unigenes比例较大。上述数据表明该转录组测序结果较好,为后续研究奠定了基础。依据黑须污蝇幼虫、蛹、成虫转录组数据筛选到30个嗅觉相关基因,其中包含9个气味结合蛋白(OBP)基因,进一步基因注释发现OBP99b,OBP56a和OBP99a基因在黑须污蝇不同发育阶段表达差异显著(|log2FC|>1,P<0.05),这在荧光定量PCR分析结果中得到了验证。【结论】本研究获得了黑须污蝇幼虫、蛹和成虫转录组数据,筛选出黑须污蝇各发育阶段嗅觉相关基因,并验证了OBP99b,OBP56a和OBP99a基因在黑须污蝇幼虫、蛹和成虫3个发育阶段差异表达,为防治双峰驼阴道蝇蛆病提供了新思路。     

绿眼赛茧蜂转录组分析及嗅觉相关蛋白基因的鉴定

[目的]绿眼赛茧蜂Zele chlorophthalmus是草地螟Loxoste gesticticalis幼虫重要天敌之一.本研究通过构建绿眼赛茧蜂触角转录组数据库,挖掘其与嗅觉相关的蛋白基因,为更好的利用绿眼赛茧蜂防治草地螟发挥其生防潜能提供理论依据.[方法]以Illumina Novaseq 6000高通量测序平台为基础,将绿眼赛茧蜂触角的基因进行转录组测序、组装序列,以及完成其生物信息学的研究分析,并对绿眼茧赛蜂触角的相关嗅觉基因做鉴定分析.[结果]成功构建绿眼赛茧蜂转录组数据库,数据库中的unigenes序列为65228条,N50为3882 bp.使用BLAST软件将绿眼赛茧蜂触角unigenes序列各自和Pfam、Swiss-Prot、NR、COG、KEGG、GO权威数据库进行对比,并完成基因功能相关注释,共注释基因数为18662条,占总数的28.61％.其中,NR数据库获得的注释最多,占总数的24.61％,为15863条,KEGG数据库获得的注释最少,为9612条(14.91％),其他依次为Pfam数据库注释数据库12164条(18.86％)、COG数据库注释15584条(24.17％)、GO数据库注释11634条(18.05％),Swiss-Prot数据库注释达到11634条,为总数的18.86％.借助GO数据库对unigenes的注释,其功能供分为三大类,可以细分49个分支,主要包括分子功能、细胞组分以及生物学过程.通过注释基因功能对嗅觉相关基因进行筛选,共发现151条和嗅觉有关的蛋白基因,包括3个感觉神经元膜蛋白基因、22个离子型受体基因、23个味觉受体基因、83个气味受体基因、6个化学感受蛋白基因、14个气味结合蛋白基因.[结论]成功收集了绿眼赛茧蜂触角转录组相关数据,并对与嗅觉相关的蛋白进行鉴定分析,为深入研究基因功能及嗅觉分子机制提供理论依据.