基于转录组数据高通量发掘扶桑绵粉蚧微卫星引物

【目的】扶桑绵粉蚧Phenacoccus solenopsis Tinsley是我国重要的检疫性害虫。简单重复序列(simple sequence repeat,SSR)研究在遗传图谱和物理图谱的构建、分子标记辅助育种、品种鉴定、基因定位、遗传多样性、动植物分类和进化等方面具有重要意义。筛选的SSR引物将为扶桑绵粉蚧遗传多样性分析、进化分析及入侵生物学等奠定基础。【方法】利用高通量搜索的方法对扶桑绵粉蚧转录组中28 120条unigenes的数据进行搜索。【结果】共找到1 781个SSR位点。扶桑绵粉蚧转录组中SSRs的主要重复类型是单核苷酸重复,占SSR总数的89.44%;其次是三核苷酸重复,占SSR总数的7.52%。单核苷酸重复里主要是A/T基序,占了总量的87.42%。基于筛选的SSRs,运用Primer 3软件进行引物的批量设计,共有481个unigenes成功设计引物,共设计出1 228对引物。【结论】研究表明利用扶桑绵粉蚧转录组数据开发SSR标记是可行的,本研究开发的引物将为扶桑绵粉蚧遗传多样性分析、进化分析及入侵生物学等奠定基础。     

基于转录组数据的齿缘刺猎蝽微卫星分子标记开发

【目的】齿缘刺猎蝽Sclomina erinacea是我国猎蝽科天敌昆虫常见种类,其不同地理种群存在明显形态差异。本研究旨在利用已经获得的齿缘刺猎蝽转录组数据筛选微卫星位点,为齿缘刺猎蝽种群遗传多样性和遗传分化研究开发可靠的微卫星分子标记。【方法】基于高通量测序技术平台Illumina HiSeq^TM 2000获得齿缘刺猎蝽转录组数据(42 215条unigenes),利用MISA软件进行搜索发掘SSR标记;利用Primer Premier 3软件设计SSR引物,从中随机选取54对SSR引物,利用PCR技术在中国齿缘刺猎蝽9个地理种群上进行验证。【结果】利用MISA软件搜索到微卫星位点2 395个,它们分布在2 107条unigenes上,其主要重复类型是三核苷酸重复(占总SSR的44.43%),其次是二核苷酸重复(占总SSR的40.08%),再次是四核苷酸重复(占总SSR的12.94%)。利用Primer Premier 3 软件成功设计出2 000余对SSR引物。随机选取的54对引物对9个齿缘刺猎蝽不同地理种群进行的SSR位点PCR扩增结果表明,共有16对引物较好地扩增出目的片段。【结论】研究表明利用齿缘刺猎蝽转录组数据可以大量发掘微卫星分子标记。本研究为进一步开展齿缘刺猎蝽的种群遗传学研究奠定了基础。     

基于RNA-seq数据的窄足真蚋SSR分子标记开发

【目的】窄足真蚋Simulium(Eusimulium)angustipes属嗜血蚋种,并传播禽类疾病。本研究旨在通过RNA-seq技术获得窄足真蚋的转录组数据,进而开发其微卫星标记,为窄足真蚋的种群遗传学研究提供可靠的分子标记。【方法】以采自陕西宝鸡和新疆阿勒泰地区的窄足真蚋幼虫为材料,构建转录组数据库。使用软件MISA(microsatellite identification tool)搜索窄足真蚋unigenes库中的所有SSR位点,并随机挑选SSR位点,利用软件Primer Premier 5.0和Oligo 6.7设计微卫星引物。经PCR扩增和电泳检测,用Cervus 3.0.7等生物信息学软件进行统计分析,筛选多态性微卫星引物。利用Blast X比对Nr和Swiss-Prot蛋白质数据库,对含多态性SSR位点的unigenes进行功能注释。【结果】共鉴定出29 471个SSRs。最丰富的重复基序是单核苷酸重复类型,占总SSR数的87.05%;其次是三核苷酸重复类型,占7.95%。共发现34种碱基重复基序,主要是A/T基序重复,占86.93%。15对微卫星引物均成功扩增出特异性产物,含12对具多态性的微卫星分子标记。比对结果显示,7个SSR位点来自注释基因序列。【结论】基于RNA-seq二代测序技术可实现窄足真蚋SSR分子标记的高效率开发,对窄足真蚋的种群遗传学研究具有重要意义。     

基于转录组数据的意大利蝗微卫星位点分析与分子标记开发

【目的】意大利蝗Calliptamus italicus (L.)是新疆荒漠半荒漠草原重要害虫。本研究利用已获得的意大利蝗转录组数据,鉴定其微卫星位点。【方法】使用MISA筛选SSR位点,利用Primer Premier 5设计引物,通过PCR扩增对引物进行验证。【结果】在意大利蝗转录组数据库中,共检测出156500个SSR位点,分布在126369条unigene中。其中,单核苷酸重复为60.88%,二核苷酸重复为23.58%,三核苷酸重复和四核苷酸重复分别为12.99%和2.04%。单核苷酸重复主要为A/T(44.38%),二核苷酸重复主要为AC/GT(12.99%)和AG/CT(8.05%)。基于筛选的SSR位点设计引物,随机挑选24对引物,对10个不同地理种群意大利蝗成虫DNA样品进行PCR扩增,共有6对引物扩增成功。【结论】本研究利用转录组数据发掘意大利蝗SSR位点,为意大利蝗分子标记、种群遗传及功能基因等研究奠定基础。     

基于转录组数据高通量发掘沙葱萤叶甲微卫星引物

【目的】沙葱萤叶甲Galeruca daurica(Joannis)是近年来在内蒙古草原上猖獗发生为害的一种新害虫。本研究从其高通量转录组数据中查找微卫星序列,并进行微卫星位点的信息分析和微卫星引物的发掘,将为进一步研究沙葱萤叶甲遗传多样性及遗传分化奠定必要的基础。【方法】应用软件MISA对沙葱萤叶甲转录组中72 352条unigenes的数据进行搜索。【结果】共找到3 880个微卫星位点,分布于3 277条unigenes中,其主要重复类型为单核苷酸重复(80.85%),其次为三核苷酸重复(11.08%),再次为二核苷酸重复(7.37%)。单核苷酸重复中主要是A/T基序,占总量的77.76%。从1 814条unigenes中成功设计出2 160对引物,从中随机选取10对引物对沙葱萤叶甲DNA进行扩增,结果全部扩增出目的 DNA片段。【结论】利用沙葱萤叶甲转录组数据开发微卫星引物是可行的,本研究开发的微卫星引物为研究沙葱萤叶甲种群遗传学和功能基因组学奠定了必要的基础。     

基于RNA-seq数据大规模开发中华蜜蜂幼虫的SSR分子标记

中华蜜蜂是东方蜜蜂的一个亚种,也是我国的特有蜂种。本研究基于已获得的中华蜜蜂幼虫肠道转录组数据预测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引物为研究中华蜜蜂分子遗传学奠定了基础。     

中华大仰蝽转录组学研究及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新标记的开发为中华大仰蝽遗传多样性分析、隐存种鉴定及基因图谱构建提供了更丰富的候选分子标记。     

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

梨小食心虫幼虫中肠转录组及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位点。【结论】本研究成功组装梨小食心虫幼虫中肠参考转录组,为以中肠为靶标的害虫防治提供理论依据。     

闽楠转录组分析及基因功能注释

采用第二代Illumina HiSeq测序技术对闽楠的木质部、韧皮部、叶片进行转录组测序，分别获得Clean Reads片段41 383 707条、43 343 922条、44 191 586条，经转录本拼接后得到序列总长度达120 535 288 bp的383 331条Conting片段，进一步组装得到平均长度为542 bp的151 729条Unigenes。将闽楠转录组Unigenes进行基因功能注释，与NR数据库比对发现，其与葡萄的相似序列最多（34%），与黄瓜、野草莓、大豆的同源性较低（各占3%）；进行GO功能注释，可将其划分为生物过程、细胞成分、分子功能3大类共计52个分支，与eggNOG数据库比对可分为25类，通过KEGG功能注释可知转录组中涉及的基因共参与了176条代谢通路，其中核糖体和碳代谢获得的注释较多。另外通过MISA软件分析，共获得35 972个SSR位点。其中，单核苷酸、二核苷酸和三核苷酸为优势重复类型，SSR位点数分别为21 762（60.50%），8 931（24.83%），4 924（13.69%）。闽楠转录组分析及基因功能注释为深入开展闽楠遗传育种及分子生物学相关研究奠定基础。     

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

[目的]中华大仰蝽Notonecta chinensis为中国和日本冲绳分布的重要水生天敌昆虫,可用于蚊虫的生物防治.本研究旨在建立中华大仰蝽转录组数据库,挖掘其基因信息.[方法]采用高通量测序平台Illumina NextSeq500对中华大仰蝽进行转录组测序、de novo组装及生物信息学分析;利用MISA软件基于...     

基于枇杷转录组序列的SSR分子标记引物开发

为获得更多的枇杷SSR引物,对枇杷转录组测序得到的1 kb以上的11 798条Unigenes进行SSR位点搜索。结果在3515 条Unigenes（6.77%）中共获得4438个SSR位点,其中主要重复类型为双碱基重复和三碱基重复,二者占SSR总数的68.27%,而四、五、六碱基重复类型较少,仅占1.42%。对选出的SSR标记采用Primer3进行引物设计,得到7911 对SSR位点特异引物,可用于枇杷遗传多样性分析、分子标记辅助育种、育种群体的建立等研究。