一种基于转录组测序开发SSR引物的方法技术

本发明专利技术公开了一种基于转录组测序开发SSR引物的方法。该方法基于已有的全基因组序列数据及其注释，针对其基因组外显子序列筛选候选SSR位点并设计SSR扩增引物。然后基于多个个体转录组测序的数据进行ePCR验证，最终筛选出多态性好的SSR位点并得到其扩增引物。本发明专利技术提供的方法与过去采用的实验室PCR与聚丙烯酰胺凝胶电泳大大减少了开发周期和人力，而且相比采用实验室PCR与测序分型的筛选方法相比大大减少了开发周期和开发经费。本发明专利技术提出的方法不仅适用于SSR引物开发，还适用于所有具有参考基因组物种的SSR引物开发，具有广泛的应用前景。

A method for developing SSR primers based on transcriptome sequencing

The invention discloses a method for developing SSR primers based on transcriptome sequencing. Based on the existing genome sequence data and its annotation, the method screened candidate SSR sites for its genomic exon sequences and designed SSR amplification primers. Then, based on the data of multiple individual transcriptome sequencing, ePCR was verified. Finally, SSR polymorphic loci were screened out and primers were amplified. The methods provided by this invention greatly reduce the development cycle and manpower of the laboratory PCR and polyacrylamide gel electrophoresis used in the past, and have greatly reduced the development cycle and development funds compared to the screening methods used in laboratory PCR and sequencing. The proposed method is not only suitable for the development of SSR primers, but also suitable for the development of all SSR primers with reference genome species, and has a wide application prospect.

【技术实现步骤摘要】
一种基于转录组测序开发SSR引物的方法
本专利技术属于SSR引物开发
，具体涉及一种开发SSR引物的新方法。
技术介绍
微卫星标记(microsatellite)，又被称为短串联重复序列(shorttandemrepeats,STRs)或简单重复序列(simplesequencerepeats)，是均匀分布于真核生物基因组中的简单重复序列，由2～6个核苷酸的串联重复片段构成，由于重复单位的重复次数在个体间呈高度变异性并且数量丰富，因此微卫星标记的应用非常广泛。微卫星标记共显性的特点使它能够用于研究等位基因，区分二倍体(或多倍体)的纯合体或杂合体，这是AFLP、RAPD等显性标记所无法做到的。另外，微卫星的特异性引物扩增具有良好的重复性和保真性，方便各实验室间的交流。目前，微卫星标记已被广泛应用到基因连锁与遗传图谱构建、遗传多样性研究、谱系和发育研究、疾病检测以及品种鉴定、亲本分析与个体、纯系检验上。随着微卫星标记在图谱上的丰富,将显现出更大的优势。但微卫星分析中引物的获得需要预先获知核酸序列,其广泛应用受限于从特定的物种中分离微卫星位点的难度和花费。微卫星标记分离最大的困难是引物的获得。同RAPD、AFLP等遗传标记不同，微卫星研究首先需要获知位点的序列信息，以便从重复序列两端侧翼的保守序列中设计引物。因而微卫星引物的开发是应用该技术的关键。目前已知微卫星位点的物种很有限，这是因为微卫星位点的获得需经过克隆、杂交筛选、测序等步骤，因此需花费一定的人力、金钱，这限制了微卫星标记的大量使用。近年来发展起来的各种微卫星位点富集技术已在一定程度上解决了这个问题。已有的最常用的获得微卫星位点的方法包括：1.从公共数据库中查找微卫星位点；2.遗传距离相近物种间引物转移扩增；3.从基因组DNA中筛选微卫星位点。其中，最方便、最经济的方法就是从已发表的文献中获得微卫星引物，但该法在使用对象上，只限于已有引物开发出的物种，而且引物的数量不会有所增加，只能停留在原有基础上。对于近缘种引物的应用，其适用范围究竟有多大；原物种的微卫星基因座在其他物种间转移扩增时，其多态性如何。这些问题使得在不同物种间共用微卫星引物时，盲目性较大，可指导操作的理论依据不足。最好的途径就是从基因组DNA中筛选微卫星位点并进行引物的开发。目前微卫星位点通常通过PCR扩增，扩增产物通过电泳分析或者测序分析并根据大小分离等位基因进行检测。事实上，从庞大的基因组数据中开发得到的SSR位点往往是数量庞大且多态性不佳的。这给我们带来的问题就是，SSR位点识别之后依然工作量巨大。研究者通常动辄需要设计几百对引物，扩增上千个样本，并对他们一一进行电泳或测序验证，不仅耗费的人力物力财力任然不可小觑，最后得到的有效SSR位点也很有限(十几对到几十对)。因此，开发一套高效快速的SSR引物开发方法对开展SSR位点的研究具有十分重大的现实意义。此外，已有SSR标记开发方法还分为基于基因组SSR开发和基于表达序列标签的SSR开发。基于表达序列开发的SSR标记在近缘种属中保守性更好，效率更高。而且基于表达数据的SSR标记是基于某一时期的基因表达序列，与基因功能和性状表型直接相关，在珍稀濒危动植物种质资源和遗传多样性评估等方面有着重要的研究价值。随着测序数据的大爆发，从转录组数据中挖掘开发新的SSR标记成为一种简单高效的方式。然而目前基于转录组测序开发SSR标记的方法依旧存在需要大量验证试验的状况。
技术实现思路
本专利技术的目的是针对现有技术存在的不足，提供一种基于转录组测序开发SSR标记引物的新方法，该方法从基因组的外显子序列中搜索SSR位点，并设计对应PCR扩增引物，先通过外显子进行ePCR验证，最后利用个体无参组装得到的转录本进行ePCR验证，最终得到多态性好的SSR标记及其引物。该方法为研究物种的遗传多样性和种质资源的鉴定提供了一种高效快速的新途径。为了实现上述技术目的，一种基于转录组测序开发SSR引物的方法，包括以下步骤：(1)至少采集2个不同个体的样本，分别提取总RNA并进行转录组测序文库构建；(2)转录组测序以及数据过滤；(3)每个样本单独进行转录本无参组装，得到每个个体的所有转录本；(4)基于外显子序列的SSR筛选与引物设计：从数据库下载参考基因组序列，提取所有外显子序列，并利用GMATA2.0软件识别SSR位点，并设计引物；将得到的引物用GMATA2.0的ePCR功能对外显子序列进行ePCR，排除专一性差的引物对；(5)SSR引物的ePCR验证：第(4)步所筛选的SSR引物用GMATA2.0的ePCR功能对第(3)步得到每个个体的转录本分别进行ePCR，从结果中筛选引物。步骤(1)提取总RNA并进行转录组测序文库构建的具体过程如下：总RNA采用QIAGENRNeasyProtectAnimalBloodKit(73224)进行提取，并用琼脂糖凝胶电泳检查RNA的质量；转录组测序文库构建采用VAHTSTMmRNA-seqv2LibraryPrepKitfor(NR602-01)试剂盒，即每个样本取1μg总RNA，然后用多聚腺苷酸磁珠分离纯化得到mRNA，用二价阳离子在VAHTSFrag/PrimeBuffer中98℃8分钟处理进行片段化；纯化之后进行cDNA第一链和第二链合成，随后进行末端补平，dATP加端和接头连接；琼脂糖凝胶电泳选择150-200bp片段进行磁珠纯化，之后进行RCR文库扩增：98℃10s变性，60℃30s结合，72℃30s延伸，循环扩增12次；最后，对文库质量检测确认之后进行测序。步骤(2)转录组测序以及数据过滤的具体过程如下：mRNA文库采用IlluminaHiseq2500进行测序；测序数据在过滤掉低质量序列和接头污染序列，低质量序列为N含量超过30％或低质量碱基含量超过10％的序列，之后用fastqc软件默认参数进行质量检测，确保数据有效。步骤(3)转录本无参组装的具体过程如下：采用Trinity无参考基因组的组装流程，默认参数对每个样本单独进行转录本组装，每个个体得到自己的转录组序列文件Trinity.fasta。步骤(4)所鉴定的SSR的搜索标准为：最小重复motif为二核苷酸，最大重复motif为十核苷酸；在外显子序列中最少重复5次。步骤(4)设计引物的条件为：扩增产物长度为120bp～400bp之间；产物侧翼序列长度为400bp；最佳退火温度Tm为60℃；最大模板长度为2000bp。步骤(4)和步骤(5)所述ePCR过程参数为：字长为12，连续字长为1，最大缺失为1，最大错配为0，扩增产物长度为100-1000bp。步骤(5)所述的筛选引物的标准为：至少在两个不同的样本中扩增产物片段差异为SSR重复单元motif长度的整数倍。相对现有技术，本专利技术的技术方案带来的有益技术效果有：基于转录组测序利用两次ePCR进行SSR引物的多态性检测验证，与现有方法相比，能够大大节省人力物力财力，是一种高效开发SSR标记的新方法。附图说明图1为本专利技术筛选得到的两对引物的电泳结果。具体实施方式以下通过实施例对本
技术实现思路
作进一步阐述，但不作为对本专利技术权利要求保护范围的限定。所有试剂盒操作按照说明书进行，转录组建库流程未作说明处按照标准流程本文档来自技高网...

【技术保护点】
一种基于转录组测序开发SSR引物的方法，其特征在于：包括以下步骤：(1)至少采集2个不同个体的样本，分别提取总RNA并进行转录组测序文库构建；(2)转录组测序以及数据过滤；(3)每个样本单独进行转录本无参组装，得到每个个体的所有转录本；(4)基于外显子序列的SSR筛选与引物设计：从数据库下载参考基因组序列，提取所有外显子序列，并利用GMATA 2.0软件识别SSR位点，并设计引物；将得到的引物用GMATA 2.0的ePCR功能对外显子序列进行ePCR，排除专一性差的引物对；(5)SSR引物的ePCR验证：第(4)步所筛选的SSR引物用GMATA 2.0的ePCR功能对第(3)步得到每个个体的转录本分别进行ePCR，从结果中筛选引物。

【技术特征摘要】
1.一种基于转录组测序开发SSR引物的方法，其特征在于：包括以下步骤：(1)至少采集2个不同个体的样本，分别提取总RNA并进行转录组测序文库构建；(2)转录组测序以及数据过滤；(3)每个样本单独进行转录本无参组装，得到每个个体的所有转录本；(4)基于外显子序列的SSR筛选与引物设计：从数据库下载参考基因组序列，提取所有外显子序列，并利用GMATA2.0软件识别SSR位点，并设计引物；将得到的引物用GMATA2.0的ePCR功能对外显子序列进行ePCR，排除专一性差的引物对；(5)SSR引物的ePCR验证：第(4)步所筛选的SSR引物用GMATA2.0的ePCR功能对第(3)步得到每个个体的转录本分别进行ePCR，从结果中筛选引物。2.根据权利要求1所述的基于转录组测序开发SSR引物的方法，其特征在于：步骤(1)提取总RNA并进行转录组测序文库构建的具体过程如下：总RNA采用QIAGENRNeasyProtectAnimalBloodKit(73224)进行提取，并用琼脂糖凝胶电泳检查RNA的质量；转录组测序文库构建采用VAHTSTMmRNA-seqv2LibraryPrepKitfor(NR602-01)试剂盒，即每个样本取1μg总RNA，然后用多聚腺苷酸磁珠分离纯化得到mRNA，用二价阳离子在VAHTSFrag/PrimeBuffer中98℃8分钟处理进行片段化；纯化之后进行cDNA第一链和第二链合成，随后进行末端补平，dATP加端和接头连接；琼脂糖凝胶电泳选择150-200bp片段进行磁珠纯化，之后进行RCR文库扩增：98℃10s变性，60℃30s结合，72℃30s延伸，循环扩增12次；最后，对文库质量检测确认之后进行测序。3.根据权利要求1所述的基于转...

【专利技术属性】
技术研发人员：刘学端，张杜，梁伊丽，胡琪，刘宏伟，郭雪，尹华群，高飞，
申请(专利权)人：中南大学，
类型：发明
国别省市：湖南,43