重要食药同源植物余甘子转录组微卫星特征分析

为全面了解余甘子转录组SSR位点的分布特征和变异规律,本研究利用Illumina Hiseq 4000平台对余甘子叶片转录组进行测序,通过MISA软件对获得的Unigenes进行SSR位点搜索和统计分析。结果发现9 538条包含SSR位点的Unigenes,共检测到9 991个SSR位点,平均每5.49 kB出现1个SSR。单碱基和二碱基为余甘子转录组SSR主要重复类型,分别占SSR总数的42.3%和30.79%。位于基因编码区的SSR位点共有1 731个,出现频率为0.039 SSRs/kB,优势重复类型为三碱基重复。余甘子转录组SSR中共有169种重复基元,其中所占比例最高的是A/T（42.10%）,其次是AG/CT（22.91%）和AAG/CTT（5.02%）。SSR各基元的重复次数波动于4~75次,且多数集中于4~20次。重复片段长度≥ 20 bp的SSR占21.20%,且SSR发生频率与片段长度呈显著负相关（P<0.01）,相关系数为-0.561。本研究获得的余甘子转录组SSR位点出现频率较高、分布密度较大、低级重复基元较多,重复次数较高、长片段较多,大多数SSR位点的多态性潜能较高,用于余甘子遗传多样性分析的潜力较大,为下一步余甘子转录组SSR标记的大规模开发和群体遗传学研究提供了重要的数据信息,进而为余甘子野生资源的保护和合理开发利用提供了参考依据。     

基于转录组测序的缺须盆唇鱼微卫星标记特征分析

缺须盆唇鱼(Placocheilus cryptonemus)是云南省怒江流域特有的一种高原经济鱼类.通过在Illumi-naHiSeq 2500平台上对缺须盆唇鱼进行转录组测序(RNA-seq)分析,利用MicroSAtellite (MISA)微卫星识别工具对1～6个碱基重复、碱基数在10 bp以上的微卫星进行鉴定.结果 表明,在缺须盆唇鱼转录组中共发现分布在36764条Unigene上的30500个SSR位点,包含序列的SSR数量为18078,发生频率为49.17％.获得的单纯SSR中共有215种重复基元,其中单核苷酸重复基元类型数量最多,为88个,六碱基重复的SSR最少,其他的重复从二、三、四、五、六碱基开始,SSR的形成呈现递减规律.SSR位点的重复次数主要发生在6～10次和11～15次重复之间,发生频率分别为38.59％和36.02％.本研究对缺须盆唇鱼的转录组数据进行测序鉴定,得到了数量较多的SSR标记,同时对得到的缺须盆唇鱼SSR数据特点进行归纳总结.本研究不仅提供了一个有价值的转录组资源,而且建立了一套SSR标记,可以为缺须盆唇鱼的基础研究、应用研究提供有效的分子标记.     

基于转录组数据的杜仲红叶性状SSR分子标记分析

对‘华仲12号’杜仲的幼嫩叶片（绿色）和成熟叶片（红色）及‘华仲11号’杜仲成熟叶片（绿色）进行转录组测序,进行测序数据的拼接和组装,且对转录组获得的基因（Unigenes）进行SSR分析。研究得到54 517条平均长度为806.90 bp的Unigenes,其中25 993条Unigenes在Nr、Swiss-Prot、KEGG和COG蛋白数据库获得功能注释,占所有Unigenes的47.68%。参照KEGG数据库,可将注释到的6 910条Unigenes划分到122个代谢途径分支,其中花色苷代谢途径相关酶基因39个,类黄酮代谢途径38个,类胡萝卜素合成途径34个。54 517条Unigenes中共包含17 010个完整型SSR位点,占总SSR位点的96.28%。完整型SSR位点共包含67种重复基元,其中出现频率最高的重复基元类型为单核苷酸重复中的A/T （7 747个）,其次是AG/CT （5 039个）和AT/AT （850个）从花色苷代谢途径、类黄酮代谢途径及类胡萝卜素代谢途径中共找到13个SSR位点。为今后杜仲遗传多样性分析、遗传图谱构建及杜仲红叶性状分子标记开发等方面奠定了分子基础。     

珍贵树种火力楠转录组SSR特征分析

为了全面了解珍贵乡土树种火力楠转录组SSR位点的分布及序列特征,为火力楠遗传资源的保存和合理开发利用提供遗传学资料,为同属植物及近缘种SSR标记的开发及遗传研究提供便利。利用Illumina Hiseq2000高通量测序平台对火力楠进行转录组测序,再通过MISA软件对测序所得Unigenes进行SSR位点的发掘和分析。分析的结果显示发现含SSR的序列21 218条,共得到27 379个SSR,出现频率为28.08%,平均约每3 kb出现1个SSR。单碱基和二碱基重复为火力楠SSR主要重复单元类型,分别占SSR总数的42.18%和35.66%,所有重复基元共85种,其中(A/T) n所占比例最高41.65%,然后是(AG/CT)n(29.47%)、(AAG/CTT)n (6.83%)和(AC/GT)n (4.11%)。在SSR和CDS的交集基因中,共发现24 668个SSR位点,其中3 805个位于编码区,出现频率为0.099 SSR/kb,而非编码区为0.287 SSR/kb,在基因编码区中出现频率最高的是三碱基重复(2 030, 53.35%)。在SSR序列长度方面,长度变化范围最大的为单碱基重复SSR,其次是二碱基重复。火力楠转录组SSR位点的出现频率高、分布密度大、基元类型丰富、重复次数较高、长片段较多,具有较高的多态性潜能,用于遗传分析的潜力很大,能满足该物种的保护遗传学研究。     

沟眶象转录组微卫星特征分析

本研究基于高通量测序获得的沟眶象Eucryptorrhynchus chinensis(Olivier)转录组数据,采用MISA(Micro Satellite)软件对其简单序列重复(simple sequence repeat,SSR)位点进行了高通量发掘。对筛选得到的1 kb以上的Unigene(17590条,占Unigene总数的26.98%)进行了SSR分析,结果发现含SSR的序列有1665条,共发现1823个SSR,平均32.52 kb出现一个SSR。沟眶象微卫星主要以三碱基重复类型为主(637,34.94%),其次为单碱基重复(576,31.60%)。共发现104种碱基重复基元,所占比例最高的为(A/T)n(30.39%),其次为(AT/AT)n(6.91%)。沟眶象SSR的平均长度为15.74 bp,长度大于20 bp的SSR仅占总数的12.70%。研究还发现沟眶象微卫星出现的频率和其长度呈极显著负相关(P0.01),相关系数为-0.572。研究结果为开发高多态性微卫星引物来进行沟眶象功能基因组学、种群遗传结构、种群遗传多样性等研究奠定了基础。     

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

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

基于转录组数据的荔枝蒂蛀虫SSR位点信息分析

荔枝蒂蛀虫Conopomorpha sinensis Bradley是专一性危害我国荔枝和龙眼的重要害虫。简单重复序列标记(Simple sequence repeat,SSR)为短串联重复序列或微卫星标记,其在荔枝蒂蛀虫偏嗜选择寄主的遗传进化机制研究和荔枝蒂蛀虫综合治理中具有重要意义。本研究基于高通量测序获得的荔枝蒂蛀虫转录组数据,利用MISA软件从68996条转录组unigenes结果中发掘出10521个SSR位点,出现频率为15.25%。荔枝蒂蛀虫转录组中SSR的主要重复类型为单碱基重复,占SSR总数的66.22%。其次是三碱基重复,占SSR总数的24.94%。在发现的33种重复基元中共筛选获得8种优势重复基元,其中A/T在单碱基重复基元中所占的比例达98.55%。基于筛选的SSR设计的9对引物中,有4对引物得到扩增预期大小的条带。荔枝蒂蛀虫SSR位点的信息分析将为探究荔枝蒂蛀虫种群遗传结构、遗传多样性和进化关系、害虫综合治理等研究提供重要科学依据。     

基于转录组的沙棘木蠹蛾简单重复序列特征分析

为开发沙棘木蠹蛾微卫星信息,利用已获得的转录组数据,对其EST-SSR位点进行发掘,进而分析其特征。结果发现含SSR的序列5126条,识别的SSR总数为7499个,SSR出现频率为51.41%。微卫星序列主要以单碱基重复为主,发生频率为39.52%。研究共发现77种碱基重复基元,所占比例最高的为(A/T)n(73.74%),其次是(AT/AT)n(3.37%)。微卫星多为重复次数为10且长度为10 bp的短序列。研究结果为沙棘木蠹蛾的SSR分子标记研究,遗传多样性分析,种群遗传结构以及关键性状基因的发掘等研究奠定基础。     

两种黄连转录组中SSR位点信息分析与多态性引物开发

本研究通过对两种黄连,即黄连(Coptis chinensis Franch.)和云南黄连(C.teeta Wall.)转录组中SSR位点信息的分析,设计SSR引物,为开发新的SSR标记提供理论基础。利用MISA工具筛选黄连及云南黄连转录组测序获得的55 903和49 741条Unigenes,对其SSR位点信息进行了分析;Primer3设计SSR引物,随机选取60对引物对两种黄连,每种20株,共40株进行多态性扩增分析。黄连和云南黄连转录组分别搜索到4 071和4 041个SSR位点。两种黄连中二核苷酸和三核苷酸重复皆是主要的类型,分别占总SSR的87.71%和87.58%。黄连中,二核苷酸重复基元出现最多的为AG/CT,共953个,占总SSR的23.40%,三核苷酸重复基元出现最多的为AAG/CTT,共746个,占总SSR的18.3%;在云黄连中,AG/CT和AAG/CTT分别为934(23.10%)和773(19.10%)。随机选取60对引物进行PCR扩增,其中12对(20.00%)表现出多态性差异。利用Populations软件作图,将40个材料分为2类。两种黄连转录组中SSR位点出现频率高,类型丰富;大量的SSR为其遗传多样性分析和遗传图谱构建提供了丰富的候选分子标记。     

基于转录组测序的茄子SSR标记开发

应用Trinity软件对茄子转录组测序数据组装,得到总长度为47919660 bp的45404条Unigenes,MISA从中检索到8316个SSR位点,发生频率为18.32%,平均5.63 kb一个位点。SSR位点中单碱基重复类型最多,为5372个,占到64.60%;其次为三碱基重复1628个,占到19.58%。三碱基重复中AAG/CTT是优势重复单元,占三碱基重复数的31.6%;二碱基重复中AG/CT是优势重复单元,占二碱基重复数的42.3%。利用Primer 3设计引物,共得到858对SSR引物,随机选取100对引物对17份茄子材料进行扩增,结果表明:有84对可以扩增出条带清晰的片段,有47对引物扩增片段为多态性片段。对47对多态性引物进行分析,多态性信息含量范围为0.10~0.64,平均多态性信息含量为0.32,UPGMA聚类分析可将17份材料分为3类。以上结果表明,基于茄子转录组测序开发的SSR标记可以为茄子的遗传多样性分析和遗传图谱构建提供更加丰富的标记来源。     

山地虎耳草转录组SSR信息分析

利用MISA（MicroSatellite）软件对山地虎耳草转录组拼接序列进行微卫星位点信息分析,为后期SSR标记的开发和物种遗传多样性检测提供候选序列。结果发现,在拼接得到的63 763条Unigene序列中含有4 622个SSR,发生频率为7.25%,有110种重复基元,平均每10.00 kB出现一个SSR位点。山地虎耳草转录组序列的SSR主要集中在三核苷酸重复（55.50%）,其次为二核苷酸重复（30.23%）。二核苷酸重复和三核苷酸重复中的优势重复基元分别为AG/TC和AAG/TTC。二核苷酸重复基元的重复次数类型最多,跨度最大,具有更高的多态性,三核苷酸次之,而四、五、六核苷酸重复类型很少。山地虎耳草转录组SSR以5~9次重复为主,且SSR数量随着重复次数的增加逐渐减少,基序长度主要集中于12~30 bp,多态性均在中等以上。     

山地虎耳草和棒腺虎耳草转录组SSR和SNP分析

利用Illumina HiSeqTM2000对山地虎耳草和棒腺虎耳草进行转录组测序,分析和比较其SSR和SNP特征。结果表明:山地虎耳草63 763条Unigene序列中含有4 622个SSR,发生频率为7.25%,有110种重复基元,平均每10.00kb出现一个SSR位点;棒腺虎耳草60 972条Unigene序列中含有4 542个SSR,发生频率为7.45%,有85种重复基元,平均每10.40kb出现一个SSR位点,略低于山地虎耳草。山地虎耳草和棒腺虎耳草转录组序列的SSR优势基元均为三核苷酸重复。2个物种的转录组SSR以5~10次的较低重复次数为主,长度主要集中于12~30bp,具有较高的多态性。山地虎耳草和棒腺虎耳草中分别获得118 424个和112 006个SNP位点,编码区的SNP位点分别占30.40%和28.59%,且在编码SNP中同义突变所占比例(30.27%、28.48%)远高于非同义突变(0.13%、0.11%)。比较发现,2个物种的各项检索结果基本一致,推测与选取的组织部位、组织的发育阶段以及物种的亲缘关系有关。     

星天牛转录组SSR位点特征分析

【目的】为了获得星天牛Anoplophora chinensis的SSR位点信息并开发其SSR分子标记技术,进一步为其遗传多样性以及综合治理提供理论依据。【方法】利用MISA软件,对星天牛转录组数据进行简单重复序列(SSR)位点筛选与分析;使用Primer3软件设计引物,采用PCR扩增以及电泳检测,筛选SSR引物,开发星天牛SSR分子标记技术。【结果】在9 325条unigene序列中共挖掘到2 360个SSR位点,出现频率为25.31%,涉及SSR位点序列1 758条,发生频率为18.85%。星天牛转录组中SSR的主要重复类型为单碱基重复,其次是三碱基重复,分别占总数的79.03%、12.54%。在核苷酸重复类型中,A/T基元种类数目最多,所占比例高达99.30%。SSR长度为10-11 bp的占比最高,为56.10%;重复次数为10次的数量最多,SSR位点数为1 188(50.34%)。重复次数和长度的分析结果对SSR位点的多态性获得了初步验证。在随机挑选序列设计的60对引物中,53对扩增产物达到预期大小,候选引物可用率高达88%,可在今后的研究中利用。【结论】本文对星天牛SSR位点的信息分析以及引物的设计与验证将有助于星天牛基因挖掘、种群遗传结构、遗传多样性、进化关系和综合治理的研究。     

基于RNA-Seq高通量测序技术的大鲵转录组分析

为发掘大鲵的功能基因,该研究以大鲵组织作为研究对象,利用RNA-seq技术对大鲵进行转录本测序和数据分析,经拼接组装共获得132 912条Unigenes,序列平均长度690 bp,N50为1 263 bp。另外从长度分布、GC含量与表达水平等方面对Unigenes进行评估,数据显示测序质量好,可信度高。此外,本研究也预测出132 416个能编码蛋白的Unigenes。使用9大数据库CDD、KOG、NR、NT、PFAM、Swiss-Prot、Tr EMBL、GO和KEGG注释大鲵转录组Unigenes,分别对应有24 049、18 406、36 711、15 858、20 500、27 515、36 705、28 879和10 958条Unigenes获得注释。其中,6 323条Unigenes在以上所有数据库中同时注释成功,39 672条Unigenes至少被一个数据库注释。KEGG分析结果显示:获得注释的10 958条Unigenes被划分到343个代谢通路中,参与信号转导通路的Unigenes数最多,共有2 644条(24.12%);其次是免疫系统通路的Unigenes有2 450条(15.29%)。最后本研究还检测到了29 790个SSR位点。通过对大鲵转录组测序,获得了大量的转录组信息,有助于进一步研究大鲵功能基因克隆、基因组学、遗传多样性分析和分子标记开发等。     

NaCl胁迫下黑果枸杞转录组测序分析

研究黑果枸杞在不同浓度盐胁迫下基因表达谱变化情况,为进一步研究黑果枸杞抗盐分子机制奠定研究基础。对0(CK)、50、250 mmol/L NaCl溶液胁迫的黑果枸杞组培苗的根和叶在胁迫时间为0、1、12 h时分别取样,采用转录组测序(RNASeq)技术进行测序分析。结果表明,转录组测序共产生222.49 Gb原始数据,拼接出Unigenes 86 037条,注释到7大功能数据库(GO、KEGG、KOG、NR、Pfam、Swiss-Prot和egg NOG)上的Unigenes总数为46 594个,占总Unigenes的54.76%,还有38 929个Unigenes在这些数据库中没有得到注释。通过GO分类和KEGG Pathway富集性分析,分别归于51个GO类别和211条代谢途径。差异表达基因分析显示,黑果枸杞叶片和根的上调基因和下调基因数随着NaCl浓度的增大和处理时间的延长均呈增加趋势,叶片中的上调基因数(7 514)小于下调基因数(9 032),根中的上调基因数(12 347)大于下调基因数(11 559)。在黑果枸杞盐胁迫下转录组中发现28 325个SSR位点,最多的为单核苷酸SSR,占70.47%。综合分析表明,黑果枸杞对盐胁迫的反应是一个多基因参与、多个生物过程协同调控的过程,基因表达量的变化可能是基因调控的主要方式。     

药用资源植物山莨菪的转录组信息分析

为了增强对资源植物山莨菪的深入了解,本研究采用高通量测序技术对山莨菪进行转录组测序分析,经过处理得到71 463个Unigenes。通过与多个数据库进行比对,对基因进行分类和分析注释,最终成功获得注释的基因有47 624条。将Unigenes比对到KOG蛋白质库中,有13 110个基因被注释,共有26个子类;比对到NR库中后有39 621个Unigenes被注释;转录本与Swissprot、Tr EMBL的比对结果得到GO功能注释信息,注释得到的29 309个Unigenes可被分为分子功能、生物学过程和细胞组分3个大类,62个子类;以KEGG数据库为参考,3 679条基因被注释,参与的代谢通路可归为4个大类,分别是代谢相关的通路、遗传信息处理、细胞过程、环境信息处理,其中与代谢相关的通路最多,约占所有代谢通路的一半。对山莨菪的药用活性成分的代谢通路及相关Unigenes数量和类型的统计结果表明,与生物碱相关的代谢通路最多,萜类和苯丙素类所对应的Unigenes数量最多。另外,结果还检测到31 382个SNP位点,6种SSR重复类型,其中单碱基重复类型所占的比例最高,每百万碱基中出现的单碱基重复的SSR个数有56. 52个,占45. 30%。该结果丰富了山莨菪的转录组信息数据,为该物种分子生物学方面的研究奠定了基础,有助于进一步开展对山莨菪的合理保护及开发利用工作。     

中间锦鸡儿转录组EST-SSR标记系统性识别与引物筛选

旨在对中间锦鸡儿转录组数据库EST信息进行SSR系统性识别和初步验证,为进一步SSR分子标记开发提供依据。对Hi Seq2000测序技术获得的中间锦鸡儿转录组Unigenes进行SSR位点搜索,共获得45 706个SSR位点,出现频率为10.38%,平均4.30kb出现一个SSR位点。SSR重复类型以单核苷酸重复序列基元为主,所占比例为56.47%;二核苷酸、三核苷酸重复序列基元的数量所占比例分别是20.56%和21.04%;其他数量的基元所占比例仅为1.9%。多核苷酸重复类型中最多的为2核苷酸重复AG/CT;其次为3核苷酸重复AAG/CTT。针对EST-SSR位点随机挑选了150对引物,通过琼脂糖凝胶电泳进行PCR验证,其中有79对能获得扩增条带,21对引物扩增出单一条带,比例为14.0%。     

红豆越橘幼叶转录组高通量测序及生物信息学分析

红豆越橘是一种新兴水果,具有较高的营养及食疗药用价值。本研究利用高通量测序技术对野生红豆越橘的幼叶转录组进行了测序,并对测序数据进行生物信息学分析。总共获得大约2 G的纯净数据,拼装了48 736条Unigenes。GO数据库注释到的Unigenes分别涉及到生物学过程、细胞成分及分子功能相关的共43种生理代谢功能;18 728条Unigenes能被KOG数据库注释,共涉及25条代谢通路;8 321条Unigenes能被KEGG数据库成功注释,共涉及到5个功能大类、20个功能中类、123条代谢通路;30 467条Unigenes可被Nr数据库注释;21 551条Unigenes可被Swissprot数据库注释;以上4个数据库共注释到30 545条Unigenes,占全部Unigenes的62.67%;被以上4个数据库均注释到的Unigene为6 075条,占全部Unigenes的12.47%。同时,生物信息学分析还显示:全部Unigenes中,有5 197条Unigenes编码转录因子,涉及57个家族;2 747条Unigenes编码抗性基因,涉及19个家族。共检测到8 099个SSR位点,其中2碱基重复的SSR位点达5 791个,占SSR位点总数的71.52%。以上研究结果对了解红果越橘的生长发育、生理生化等分子机制,对进一步进行红豆越橘分子方面的相关研究均具有较强的参考价值。     

基于转录组测序的波纹唇鱼SSR和SNP多态特征分析

波纹唇鱼是一种极度濒危的珍贵观赏鱼和高价值食用鱼,全面了解其转录组中SSR和SNP位点的分布及序列特征,可有助于开展波纹唇鱼遗传资源的保存和合理开发利用。利用二代高通量RNA-seq测序技术对波纹唇鱼进行转录组测序,通过MISA和Samtools对所得Unigene进行SSR与SNP位点的发掘与分析。结果显示,在150 218条Unigene序列中共发现22 428个SSR,出现频率为14.93%,平均约5.35 kb出现1个SSR。波纹唇鱼SSR的主要重复单元类型为单碱基和二碱基重复,分别占SSR总数的61.32%和19.12%,除复合类型以外,所有重复基元共65种,其中(A/T)n所占比例最高(55.64%),(AG/CT)n和(AC/GT)n分别占8.31%和6.80%。22 427个SSR处于CDS中,其中1 773个位于编码区。SSR重复次数集中在5～15次,序列平均长度为13.3 bp。3 438个SSR位点共获得669对候选引物。随机选取50对验证,发现21对可扩增出与预期产物长度大小一致的特异性条带,其中5对在4个波纹唇鱼个体间具有多态性。在Unigene中共发现245 373个SNP位点,发生频率为1/490 bp,其中转换类型发生频率显著高于颠换类型,转换类型中A/G (33.61%)和C/T (32.51%)发生频率最高,颠换类型中则是A/T的最高(11.12%)。波纹唇鱼转录组中SSR和SNP位点非常丰富,可为波纹唇鱼遗传多样性分析、亲缘关系鉴定与遗传资源开发利用等方面提供丰富的基础数据信息。     

陆地棉EST长度多态性与其SSR分布特征相关性分析

目的:分析陆地棉EST长度多态性与其SSR分布特征的相关性。方法:从NCBI公共数据库下载陆地棉EST序列,应用SSRIT搜索SSR,分析20 000条无冗余的EST序列。结果:在剔除低质量和冗余的序列后,得到全长为7 363.878kb的无冗余EST序列7 322条,其中含有SSR位点的EST序列数520条,占被分析EST比例的2.60%。长度在400bp以下的EST序列含SSR的比例为1.46%;长度在400bp以上的EST序列含SSR的比例为8.94%。在1～6bp的重复基元中,二核苷酸重复基元的SSR重复频率最高,占总数的63.46%,其次是三核苷酸,占总数的34.04%。二核苷酸类型(AG)n、(AT)n和三核苷酸类型(AAG)n、(ACC)n、(ACT)n、(AAT)n是SSR的主要重复基元。结论:棉花EST-SSR可用于棉花分子标记,为有针对性设计陆地棉EST-SSR引物奠定基础。