一种基于杜鹃花转录组数据的ssr分子标记引物及其筛选方法与应用
【专利摘要】本发明公开了基于转录组数据开发的杜鹃花SSR分子标记引物及其具体应用，属于生物技术领域。该引物是基于转录组序列开发得到，在转录组测序的基础上，对大量测序数据进行批量筛选，获取SSR序列并进行SSR标记的引物设计，克服了目前杜鹃花SSR分子标记数量少、开发效率低等障碍。利用杜鹃花不同品种验证了SSR引物的有效性，为杜鹃花遗传多样性研究、遗传连锁图谱构、分子标记辅助育种等奠定了基础。
【专利说明】
一种基于杜鹃花转录组数据的SSR分子标记引物及其筛选方 法与应用
技术领域
[0001] 本发明涉及分子标记技术开发及应用领域。更具体涉及基于杜鹃花转录组数据的 SSR分子标记引物开发和应用，该引物适用于杜鹃花遗传多样性鉴定、遗传连锁图谱构建、 QTL定位、分子标记辅助育种等。
【背景技术】
[0002] 杜鹊花（Rhododendron simsii Planch)为杜鹊花科（Ericaceae)杜鹊花属 (Rhododendron)常绿或者落叶灌木，重要的园林绿化和观赏植物。杜鹃花一般春季开花，每 簇花2-6朵，花冠漏斗形，有红、淡红、杏红、雪青、白色等，花色繁茂艳丽，为著名的花卉植 物，具有较高的观赏价值。该物种全株供药用有行气活血、补虚，治疗内伤咳嗽，肾虚耳聋， 月经不调和风湿等。长期以来，众多育种学者对杜鹃花种质资源的收集、分类和鉴定做了大 量工作，为杜鹃花育种研究提供了大量优良性状材料，选育出了多种符合生产需求的杜鹃 花品种。分子标记作为材料品种鉴定、遗传图谱构建和基因定位的基础，已经成为花卉研究 和应用的重要技术手段，但是目前对杜鹃花育种中应用的分子标记十分有限，限制了其在 杜鹃花优质品种资源快速、有效的选育和应用开发。
[0003] 简单重复序列（Simple sequence repeats，SSR)广泛分布在真核生物基因组上， 由于其序列重复类型和次数不同，使其具有高度的多态性。与其他分子标记相比，SSR标记 具有多态信息量高，共显性遗传、重复性好、特异性强和技术简单等特点，是较为可靠的分 子标记类型。但是使用该技术必须要有相应的SSR序列做支撑，目前SSR标记可分为基因组 SSR(gSSR)和转录组SSR(EST-SSR)，由于杜鹃花的基因组尚未测序或者测序数据尚未公布， 而其同源物种基因组开发的SSR尚无确定是否能与其通用，所以目前现有技术中尚未有利 用杜鹃花基因组SSR或转录组数据开发SSR标记的相关报道。
【发明内容】
[0004] 本发明的目的在于提供一种基于杜鹃花转录组数据的SSR分子标记引物筛选方法 和引物。我们利用高通量测序技术对杜鹃花花蕾材料进行转录组测序，转录组数据进行组 装拼接后挖掘带有SSR标记的序列，将符合SSR引物设计要求的序列进行批量SSR引物开发， 并随机选取部分对SSR引物进行适用性验证。该引物的开发为后续杜鹃花的品种鉴定、分子 标记辅助育种、重要性状基因定位和克隆等方面具有重要意义。
[0005] 为了实现上述目的，本发明首先提供了一种基于杜鹃花转录组数据的SSR分子标 记引物的筛选方法，包括以下步骤：
[0006] (1)提取总 RNA;
[0007] (2)构建杜鹃花转录组测序文库，再对所述转录组测序文库进行高通量测序；
[0008] (3)对所述转录组测序文库的数据进行过滤与组装，得到Unigene数据库；
[0009] (4)对所述Unigene数据库的序列进行SSR位点搜索，筛选出符合SSR引物设计的序 列；
[0010] (5)对含有SSR位点的所述Unigene数据库的序列进行批量SSR引物设计，获得SSR 分子标记引物。
[0011] 优选地，本发明的基于杜鹃花转录组数据的SSR分子标记引物的筛选方法中，所述 步骤(1)的RNA选自杜鹃花早期花蕾材料;所述早期花蕾材料为获取后用锡箱纸包裹置于液 氮中速冻，迅速转移至-70°C冰箱保存；采用高盐Trizol试剂和RNA结合柱的方法提取总 RNA〇
[0012] 优选地，本发明的基于杜鹃花转录组数据的SSR分子标记引物的筛选方法中，所述 步骤(2)中，使用特异性转录组文库构建试剂盒构建了所述杜鹃花转录组测序文库，所述特 异性转录组文库构建试剂盒的操作步骤为:采用特异性试剂盒对总RNA先进行mRNA 01 igo (dT)磁珠富集;富集到的mRNA进行双链cDNA合成、末端修复和接头连接;再用USER酶对所述 双链cDNA第二链降解，保留真实转录的mRNA第一链信息；后续进行PCR扩增和2%的琼脂糖 凝胶电泳，回收300_500bp的片段，即转录组测序文库；所述转录组文库采用Illumina Hiseq2500PE125进行高通量测序。
[0013] 优选地，本发明的基于杜鹃花转录组数据的SSR分子标记引物的筛选方法中，所述 步骤(3)中，所述数据的过滤采用FASTX和⑶TADAPT软件，所述数据的组装采用Trinity软 件。
[0014] 优选地，本发明的基于杜鹃花转录组数据的SSR分子标记引物的筛选方法中，所述 步骤(4)中，使用MISA软件对Unigene数据库的序列进行SSR位点搜索，所述搜索的标准为二 核苷酸、三核苷酸、四核苷酸、五核苷酸、六核苷酸最少重复次数分别为10、7、5、5和5次。
[0015] 优选地，本发明的基于杜鹃花转录组数据的SSR分子标记引物的筛选方法中，所述 步骤(5)中，用primer 3 · 0软件批量设计所述含有SSR位点的Unigene数据库的序列设计引 物，并且所述SSR位点离侧翼序列长度大于50bp，所述引物退火温度54-65Γ之间，PCR产物 大小80-300bp，所述引物长度18-25bp之间，所述引物的CG含量40 % -60 %。
[0016]更优选地，本发明杜鹃花SSR引物开发和应用的具体流程为:选取杜鹃花材料提取 总RNA，构建了杜鹃花转录组测序文库，采用Illumina二代测序仪对转录组文库进行高通量 测序。测序数据进过严格的数据过滤后，对转录组数据进行组装，得到Unigene基因库序列， 将其作为后续SSR引物开发的背景数据。对Unigene基因库序列进行SSR位点搜索，筛选出符 合SSR引物设计的序列。利用相关软件对含有SSR位点的Unigene序列进行SSR引物设计，并 将设计好的SSR引物对多个杜鹃花品种进行PCR扩增验证其有效性。
[0017]本发明的另一目的在于，本发明通过上述筛选方法获得了一种基于杜鹃花转录组 数据的SSR分子标记引物，包括以下引物对：
[0020] 优选地，本发明所述的SSR分子标记引物中，所述引物对的退火温度54_65°C之间， 所述引物长度18_25bp之间，所述引物的CG含量40%-60%，PCR产物大小80-300bp。
[0021] 优选地，本发明所述的SSR分子标记引物中，所述引物对退火温度为：
[0024] 本发明的一种基于杜鹃花转录组数据的SSR分子标记引物筛选方法，是大规模SSR 引物开发中效率最高和最为切实可行的策略之一。本发明中开发了40对杜鹃花SSR引物，并 对7个不同杜鹃花品种进行PCR扩增效率和多态性检测。本发明的特点在于：
[0025] 1、本发明通过转录组数据获取杜鹃花的SSR引物开发的背景数据，与基因组上的 gSSR相比具有更高的通用性；
[0026] 2、本发明的筛选方法获得的SSR分子标记引物，具有以下特点：引物退火温度54- 65°C之间，PCR产物大小80-300bp，引物长度18-25bp之间，CG含量40%-60%，通过本发明的 后续实施例表明，本发明设计的引物可以避免产生二聚体结构，更有利于后续的操作；
[0027] 3、本发明的40个EST-SSR分子标记引物在7个杜鹃花品种中有多态性，可将本发明 提供的标记应用于更多的杜鹃花植物进行种质资源遗传多样性、遗传图谱构建、重要性状 定位和分子标记辅助育种等方面。
【附图说明】
[0028]图1是本发明基于转录组SSR引物筛选方法的流程示意图；
[0029 ]图2是实施例中40对引物在7个杜鹃花品种中6 % PAGE胶电泳检测。
【具体实施方式】
[0030] 本发明实施例中基于杜鹃花转录组数据的SSR分子标记引物的筛选方法，包括以 下步骤：
[0031] (1)提取总 RNA;
[0032] (2)构建杜鹃花转录组测序文库，再对所述转录组测序文库进行高通量测序；
[0033] (3)对所述转录组测序文库的数据进行过滤与组装，得到Unigene数据库；
[0034] (4)对所述Unigene数据库的序列进行SSR位点搜索，筛选出符合SSR引物设计的序 列；
[0035] (5)对含有SSR位点的所述Unigene数据库的序列进行批量SSR引物设计，获得SSR 分子标记引物。
[0036] 本发明实施例中基于杜鹃花转录组数据的SSR分子标记引物的筛选方法中，所述 步骤(1)的RNA选自杜鹃花早期花蕾材料;所述早期花蕾材料为获取后用锡箱纸包裹置于液 氮中速冻，迅速转移至-70°C冰箱保存；采用高盐Trizol试剂和RNA结合柱的方法提取总 RNA〇
[0037] 本发明实施例中基于杜鹃花转录组数据的SSR分子标记引物的筛选方法中，所述 步骤（2)中，使用特异性转录组文库构建试剂盒NEB Next Poly(A)mRNAMagnetic Isolation Module构建了所述杜鹊花转录组测序文库，其他可商购的转录组文库试剂盒也 可以用于本发明，所述特异性转录组文库构建试剂盒的操作步骤为:采用特异性试剂盒对 总RNA先进行mRNA 01 igo (dT)磁珠富集；富集到的mRNA进行双链cDNA合成、末端修复和接头 连接;再用USER酶对所述双链cDNA第二链降解，保留真实转录的mRNA第一链信息；后续进行 PCR扩增和2%的琼脂糖凝胶电泳，回收300-500bp的片段，即转录组测序文库;所述转录组 文库采用Illumina Hiseq2500PE125进行高通量测序。
[0038] 本发明实施例中基于杜鹃花转录组数据的SSR分子标记引物的筛选方法中，所述 步骤(3)中，所述数据的过滤采用FASTX和⑶TADAPT软件，所述数据的组装采用Trinity软 件。
[0039] 本发明实施例中基于杜鹃花转录组数据的SSR分子标记引物的筛选方法中，所述 步骤(4)中，使用MISA软件对Unigene数据库的序列进行SSR位点搜索，所述搜索的标准为二 核苷酸、三核苷酸、四核苷酸、五核苷酸、六核苷酸最少重复次数分别为10、7、5、5和5次。
[0040] 本发明实施例中基于杜鹃花转录组数据的SSR分子标记引物的筛选方法中，所述 步骤(5)中，用primer 3 · 0软件批量设计所述含有SSR位点的Unigene数据库的序列设计引 物，并且所述SSR位点离侧翼序列长度大于50bp，所述引物退火温度54-65Γ之间，PCR产物 大小80-300bp，所述引物长度18-25bp之间，所述引物的CG含量40 % -60 %。
[0041] 本发明实施例中杜鹃花SSR引物开发和应用的具体流程为:选取杜鹃花材料提取 总RNA，构建了杜鹃花转录组测序文库，采用Illumina二代测序仪对转录组文库进行高通量 测序。测序数据进过严格的数据过滤后，对转录组数据进行组装，得到Unigene基因库序列， 将其作为后续SSR引物开发的背景数据。对Unigene基因库序列进行SSR位点搜索，筛选出符 合SSR引物设计的序列。利用相关软件对含有SSR位点的Unigene序列进行SSR引物设计，并 将设计好的SSR引物对多个杜鹃花品种进行PCR扩增验证其有效性。
[0042] 实施例基于杜鹃花转录组数据的SSR分子标记引物的筛选和在检测杜鹃花品种多 态性的应用
[0043]本实施例中提供了40对SSR引物是从杜鹃花转录组数据中挖掘得到，使用MISA软 件对Unigene进行SSR位点搜索，搜索标准为:二、三、四、五、六核苷酸最少重复次数分别为 10、7、5、5和5次。用primer 3.0软件设计含有SSR位点的Unigene序列设计引物，并且SSR位 点离侧翼序列长度大于50bp，引物退火温度54-65°C之间，PCR产物大小80-300bp，引物长度 18-25bp之间，CG含量40%-60%，引物避免产生二聚体结构。
[0044]具体包括以下步骤：
[0045] (1 )DNA提取。利用2 X CTAB法提取7种杜鹃花叶片组织的基因组DNA。叶片经液氮研 磨、酚/氯仿(体积比1:1)抽提、异丙醇沉淀、75%乙醇洗涤和RNA酶处理等步骤后，将DNA溶 于50yL TE溶液中。本发明实验方案中所用到的不同杜鹃花品种信息见表1。
[0046]表1本发明所采用的7个杜鹃花品种信息
[0049] (2)转录组SSR序列筛选。选取杜鹃花材料提取总RNA，构建杜鹃花转录组测序文 库，采用Illumina二代测序仪对转录组文库进行高通量测序。测序数据进过严格的数据过 滤后，对转录组数据进行组装，得到Unigene基因库，将其作为后续SSR引物开发的背景数 据。对Unigene序列进行SSR位点搜索，筛选出符合SSR引物设计的序列。其中，拼接得到的 Unigene基因库的59620个序列，含SSR序列的Unigene序列有9930个，从中统计出SSR信息， 见表2。
[0050] 表2 Unigene基因库中SSR重复序列和重新次数统计
[0052]使用MISA软件对获得的Unigene基因库基因进行SSR位点搜索，搜索标准为：二核 苷酸、三核苷酸、四核苷酸、五核苷酸、六核苷酸最少重复次数分别为1〇、7、5、5和5次。用 primer 3.0软件批量设计含有SSR位点的Unigene序列引物，并且SSR位点离侧翼序列长度 大于50bp，引物退火温度54-65°C之间，PCR产物大小80-300bp，引物长度18-25bp之间，CG含 量40 %-60 %，引物避免产生二聚体结构。SSR引物序列见表3。引物合成由南京金斯瑞合成， 采用RPC的方式进行纯化。
[0053] 表3 SSR引物序列
[0056] (3)PCR 扩增。15μ1 反应体系包含：ddH20 11 · 9μ1，10 X Buf f er (含 Mg2+) 1 · 5μ1，dNTPs (lOmM)，0·2μ1，ΙΟμΜ正反引物各0· 2μ1，Taq聚合酶0· 2μ1，DNA模板（50ngAU)ΙμL。反应程序 为：94°C 5min;94°C 30sec，退火温度（见表3)30sec，72°C SOsecJOcyclesJST： 5min;4 °C hold。所选引物序列见表3。
[0057] (4)6%PAGE胶检测。将上述扩增产物取3μ1，加 2μ1 6 XDNA loading buffer，进行 6 % PAGE胶电泳，PAGE胶分别用硝酸银溶液和氢氧化钠溶液进行染色和显色。结果见图2:6 个随机挑选的 EST-SSR 引物（SSRl-6:RhE-l、RhE-12、RhE-14、RhE-15、RhE-19、RhE-22)在7 个 杜鹃花品种中DNA的扩增结果:六个标记在不同杜鹃品种内的扩增条带存在明显的多态性； 其中SSR1、SSR2、SSR4、SSR5、SSR6五个标记还能区分出纯合子和杂合子。
[0058]以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人 员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应 视为本发明的保护范围。
【主权项】
1. 一种基于杜酷花转录组数据的SSR分子标记引物，包括w下引物对：2. 根据权利要求1所述的SSR分子标记引物，其特征在于:所述引物对的退火溫度54- 65°(：之间，所述引物长度18-2化口之间，所述引物的〔6含量40%-60%，口〔1?产物大小80- 300bp。3. 根据权利要求1所述的SSR分子标记引物，其特征在于:所述引物对退火溫度为：4. 一种基于杜酷花转录组数据的SSR分子标记引物的筛选方法，其特征在于，包括W下 步骤： (1) 提取总3臟； (2) 构建杜酷花转录组测序文库，再对所述转录组测序文库进行高通量测序； (3) 对所述转录组测序文库的数据进行过滤与组装，得到Unigene数据库； (4) 对所述化igene数据库的序列进行SSR位点捜索，筛选出符合SSR引物设计的序列； (5) 对含有SSR位点的所述化igene数据库的序列进行批量SSR引物设计，获得SSR分子 标记引物。5. 根据权利要求4所述的筛选方法，其特征在于，所述步骤（1)的RNA选自杜酷花早期花 蕾材料;所述早期花蕾材料为获取后用锡锥纸包裹置于液氮中速冻，迅速转移至-70°C冰 箱保存;采用高盐化izol试剂和RNA结合柱的方法提取总RNA。6. 根据权利要求4所述的筛选方法，其特征在于:所述步骤(2)中，使用特异性转录组文 库构建试剂盒构建了所述杜酷花转录组测序文库，所述特异性转录组文库构建试剂盒的操 作步骤为:采用特异性试剂盒对总RNA先进行mRNA Oligo(dT)磁珠富集;富集到的mRNA进行 双链cDNA合成、末端修复和接头连接;再用US邸酶对所述双链cDNA第二链降解，保留真实转 录的mRNA第一链信息；后续进行PCR扩增和2 %的琼脂糖凝胶电泳，回收300-50化P的片段， 即转录组测序文库;所述转录组文库采用111皿ina化seq2500阳125进行高通量测序。7. 根据权利要求4所述的筛选方法，其特征在于:所述步骤(3)中，所述数据的过滤采用 FASTX和CUTADAPT软件，所述数据的组装采用Trinity软件。8. 根据权利要求4所述的筛选方法，其特征在于：所述步骤（4)中，使用MISA软件对 Unigene数据库的序列进行SSR位点捜索，所述捜索的标准为二核巧酸、Ξ核巧酸、四核巧 酸、五核巧酸、六核巧酸最少重复次数分别为1〇、7、5、5和5次。9. 根据权利要求4所述的筛选方法，其特征在于:所述步骤(5)中，用primer3. ο软件批 量设计所述含有SSR位点的化igene数据库的序列设计引物，并且所述SSR位点离侧翼序列 长度大于50bp，所述引物退火溫度54-65Γ之间，PCR产物大小80-300bp，所述引物长度18- 25bp之间，所述引物的CG含量40 % -60 %。10. 根据权利要求1-3任意一项所述的SSR分子标记引物在杜酷花的品种鉴定中的应 用。
【文档编号】C12N15/11GK105969767SQ201610564377
【公开日】2016年9月28日
【申请日】2016年7月18日
【发明人】王书珍, 方元平, 金卫斌, 向福, 项俊, 程华
【申请人】黄冈师范学院