鉴别或辅助鉴别哺乳动物物种的探针组合物及其试剂盒与应用的制作方法

1.本发明涉及生物技术领域中鉴别或辅助鉴别哺乳动物物种的探针组合物及其试剂盒与应用。


背景技术：

2.细胞色素c氧化酶i基因(coi基因)是线粒体基因编码的3种细胞色素氧化酶亚基之一，它是其中分子量最大、功能结构最为保守的基因。coi基因具有多变异、易被通用引物扩增、序列本身又很少存在插入和缺失等特点，因此coi基因被选为dna分类的标记基因(dna条形码)，该段编码基因的长度一般大约在658bp左右，其除可以用于dna分类外，还可以用于物种的系统发育关系和分子进化的研究。
3.以往获得coi基因序列，一般采用一代测序的方法，该方法有灵活、方便、可以快速得到结果的优点；但同时也存在批量样本操作繁琐，一代测序部分结果不稳定造成无法得到准确coi基因序列的风险。
4.发明公开
5.本发明所要解决的一个技术问题是如何批量鉴别或辅助鉴别哺乳动物物种。
6.为了解决以上技术问题，本发明提供了鉴别或辅助鉴别哺乳动物物种的探针组合物。
7.本发明所提供的鉴别或辅助鉴别哺乳动物物种的探针组合物为对哺乳动物coi基因序列进行序列聚类得到代表性序列，针对代表性序列，每个snp位点设计探针覆盖，得到的探针组合物。
8.上述探针组合物，所述序列聚类为使用angiosperms353方法进行序列聚类，设置遗传距离为0.05(即序列相似度为95％)，覆盖深度2x。
9.上述探针组合物中，所述设计探针，按照gc含量＞30％的标准进行。
10.上述探针组合物中，所述每个snp位点设计探针覆盖，为每个snp位点设计2条探针覆盖。
11.本发明所提供的鉴别或辅助鉴别哺乳动物物种的探针组合物具体为由序列表序列1-序列3590所示的3590个单链dna组成的组合。
12.本发明还提供一种利用上述探针组合物鉴别或辅助鉴别哺乳动物物种的方法：包括以上述探针组合物捕获待测哺乳动物的coi基因，建库，通过二代测序高通量检测获得coi基因的dna序列，根据得到的coi基因序列确定待测哺乳动物物种。
13.上述方法中，所述根据得到的coi基因序列确定待测哺乳动物物种，为与已知物种coi基因比对，例如与线粒体全基因数据中的coi基因比对。
14.上述方法中，可批量鉴别哺乳动物物种。
15.本发明还保护鉴别或辅助鉴别哺乳动物物种的试剂或试剂盒，所述试剂或试剂盒包括所述的探针组合物。
16.上述探针组合物、上述方法和/或上述试剂或试剂盒在鉴别或辅助鉴别哺乳动物物种中的应用，或在制备鉴别或辅助鉴别哺乳动物物种产品中的应用属于本发明的保护范围。
17.本发明应用靶向测序基因型分型技术，根据评估分析选择的序列设计合成液相探针，并对探针进行捕获效率测试，最终形成coi基因捕获试剂盒，实现批量进行哺乳动物鉴别的目的。
附图说明
18.图1为本发明实施例1中探针组合开发的流程图。
具体实施方式
19.下面结合具体实施方式对本发明进行进一步的详细描述，给出的实施例仅为了阐明本发明，而不是为了限制本发明的范围。以下提供的实施例可作为本技术领域普通技术人员进行进一步改进的指南，并不以任何方式构成对本发明的限制。
20.下述实施例中的实验方法，如无特殊说明，均为常规方法。下述实施例中所用的材料、试剂等，如无特殊说明，均可从商业途径得到。
21.实施例1
22.发明人根据不同物种间coi序列的多态性，设计捕获探针，用于捕获目标物种的coi基因，然后建库用二代测序高通量检测的方法获得coi基因的dna序列。该方案可以同时兼顾一代测序灵活、方便、可以快速得到结果的优点，同时解决一代测序中存在的批量样本操作繁琐的问题，以及避免一代测序部分结果不稳定造成无法得到准确coi基因序列的风险；将捕获探针开发成哺乳动物coi基因捕获试剂盒，用于大量哺乳动物物种dna条形码研究，可实现对哺乳动物物种的鉴别。
23.具体开发过程按照图1中的流程图进行：
24.1、技术路线及方案的确定
25.两种二代测序技术genobaits和genoplexs的对比见表1:
26.表1 genobaits和genoplexs技术的对比
27.[0028][0029]
发明人应用靶向测序基因型分型(genotyping by target sequencing,gbts，又称genobaits)技术，根据评估分析选择的序列设计合成液相探针，并对探针进行捕获效率测试，最终形成coi基因捕获试剂盒，实现批量进行哺乳动物鉴别的目的。
[0030]
根据300多种哺乳动物的coi基因序列，以及部分从线粒体全基因数据中提取的coi基因，共413条coi基因进行序列聚类，目的是找到每一个分类中最具有代表性的coi序列，并以最终聚类得到的序列进行探针设计。因设计时使用不同参数有不同效果，给出两套评估方案，具体如下：
[0031]
方案一：使用angiosperms353方法进行序列聚类，设置遗传距离为0.1即序列相似度为90％，覆盖深度1x；以此参数进行聚类，最终得到的代表性序列为378条，开发及检测成本较低，存在风险是当目标序列再次发生变异的时候即检测时捕获到的序列和提供的序列都存在差异的时候，可能会出现无法捕获相应序列，从而影响物种鉴结果的情况；
[0032]
方案二：使用angiosperms353方法进行序列聚类，设置遗传距离为0.05即序列相似度为95％，覆盖深度2x；以此参数进行聚类，最终得到的代表性序列较多，为479条，开发及检测成本较高，对目标序列捕获的概率会大幅提升，最终鉴定结果的准确性更高。
[0033]
最终确认以方案二进行试剂盒的开发，进行探针设计。
[0034]
2、探针设计和挑选原则
[0035]
2.1探针设计原则：
[0036]
用genobaits probe designer软件对方案二中获得的479条代表性序列进行液相捕获探针设计，探针长度设置为110bp，gc含量＞30％，每个snp位点设计由2条捕获探针覆盖。
[0037]
2.2探针挑选原则：
[0038]
1)选取探针含量在30％-80％之间的；
[0039]
2)选取同源性区域个数＜10的；
[0040]
3)选取探针区域不包含ssr、n区域的。
[0041]
探针设计结果见表2:
[0042]
表2探针设计结果
[0043][0044]
注：该评估结果为探针设计结果，并非只要设计出探针就一定能捕获目标位点。
[0045]
挑选出针对哺乳动物的探针总数3590个，每个长度110bp。利用芯片原位合成技术合成带有生物素(biotin，b)修饰的核苷酸探针序列(生物素位于探针5’端)。
[0046]
3、引物合成及开发测试
[0047]
合成步骤2挑选的3590个探针，对7个测试物种共计10个样品进行测试。7个测试物种分别为牦牛(样本号qh421、qh565)、盘羊(样本号m-15、m-11)、北山羊(样本号661、660)、白颊长臂猿(样本号b12)、赤麂(样本号62)、黔金丝猴(样本号125)、黑麂(样本号188)，所用样本为血液样本。所有样品在其体检时获得，样本获取经北京动物园学术委员会的审核及许可。
[0048]
3.1样品dna质检
[0049]
将每个测试样品dna用qubit fluorometric quantitation(thermo fisher公司)对dna浓度进行测定，用1％琼脂糖凝胶电泳检测dna的完整性。检测合格的样品放入4℃冰箱，保存、备用。
[0050]
3.2样品dna测序文库构建
[0051]
针对每个测试样，取12μl步骤3.1质检合格的dna(200ng)放置于0.2μl pcr管中，将管置于超声波破碎仪中对dna进行随机物理破碎，片段破碎至200-400bp。然后向管中加入4μl genobaits endrepairbuffer(博瑞迪公司)和2.7μl genobaits end repair enzyme(博瑞迪公司)，补水至20μl，放入abi 9700pcr仪中37℃温育20分钟，完成破碎片段的末端修复和加a过程。
[0052]
从pcr仪中取出小管加入2μl genobaits ultra dna ligase(博瑞迪公司)、8μlgenobaits ultra dna ligase buffer(博瑞迪公司)和2μl genobaits adapter(博瑞迪公司)，补水至40μl，然后放置于abi9700 pcr仪上22℃反应30分钟，完成测序接头的连接。向连接产物中加入48μl的beackman ampure xp beads(beackman公司)对连接产物进行纯化，纯化后按照0.65 0.2倍磁珠进行片段筛选，保留插入片段在200-300bp的连接产物。
[0053]
向上一步的pcr管中加入5μl带有barcode序列(所用barcode序列选自序列表序列3591-序列3686中的barcode序列，不同的barcode用于区分不同的样品)的测序接头、1μl p5接头、10μl genobaits pcr master mix(博瑞迪公司)，并用纯水补至20μl；用abi 9700pcr仪进行扩增，扩增程序为：95℃预变性5min，95℃变性30s，60℃退火30s，72℃延伸30s；重复2-4步，共8个循环；72℃延伸5min。
[0054]
向第二轮pcr产物中加入24μl beckmen ampure xp beads(beckmen公司)，用移液器上下吸打均匀后，将0.2μl的pcr管置于磁力架上至溶液澄清，弃去上清并用75％乙醇洗涤磁珠一次，用ph值为8.0的tris-hcl将文库dna洗脱下来。
[0055]
3.3样品dna杂交捕获
[0056]
取500ng已完成构建的样品dna测序文库，加入5μl genobaits block i(博瑞迪公司)和2μl genobaits block ii(博瑞迪公司)，置于eppendorfconcentratorplus(eppendorf公司)真空浓缩仪上，在≤70℃的温度下蒸干至干粉。向干粉管中加入8.5μl genobaits 2x hyb buffer(博瑞迪公司)、2.7μl genobaits hyb buffer enhancer(博瑞迪公司)、2.8μl nuclease-free water，用移液器吸打混匀后放置于abi 9700pcr仪上95℃温育10分钟，然后取出pcr管加入3μl已经合成的探针(60ng/ul)，旋涡震荡混匀后放置于abi 9700pcr仪上65℃温育2小时，完成探针杂交反应。
[0057]
向上一步杂交完成的反应体系中加入100μl genobaits dna probe beads(博瑞迪公司)，上下吸打10次，放入abi 9700pcr仪上65℃温育45分钟，使磁珠与探针结合。用100μl genobaits wash buffer i(博瑞迪公司)、150μl genobaits wash bufferii(博瑞迪公司)分别对结合探针后的磁珠进行65℃热洗，然后再用100μl genobaits wash buffer i(博瑞迪公司)、150μl genobaits wash buffer ii(博瑞迪公司)和150μlgenobaits wash buffer iii(博瑞迪公司)分别对磁珠进行常温洗涤。洗涤完成的磁珠用20μl nuclease-free water进行重悬。
[0058]
取13μl重悬后的dna(带磁珠)加入到新的0.2ml pcr管中，然后加入15μlgenobaits pcr master mix(博瑞迪公司)、2μl genobaits primer mix(博瑞迪公司)配置post-pcr体系，用abi 9700pcr仪进行文库扩增，扩增程序为：95℃预变性5min，95℃变性30s，60℃退火30s，72℃延伸30s；重复2-4步，共15个循环；72℃延伸5min。
[0059]
向post-pcr产物中加入45μl beckmen ampure xp beads(beckmen公司)并用移液器上下吸打均匀，然后将0.2ml pcr管置于磁力架上至溶液澄清，弃去上清并用75％乙醇洗涤磁珠两次，用ph为8.0的tris-hcl将文库dna洗脱下来。完成测试样品的杂交捕获工作。
[0060]
3.4样品dna杂交捕获文库的质检与测序
[0061]
用qubit fluorometric quantitation(thermo fisher公司)对文库的dna浓度进行测定，然后用琼脂糖凝胶电泳检测文库dna的片段大小是否在300-400bp之间。构建好的dna文库用illumina hiseq x ten测序仪进行测序。
[0062]
3.5测试数据分析
[0063]
测序得到的原始测序序列(sequenced reads)或者raw reads，里面含有带接头的、低质量的reads。对于二代测序技术，测序错误率分布具有下面两个主要原因：1)由于测序过程中化学试剂的消耗，导致测序错误率会随着测序序列(sequenced reads)长度的增加而升高。2)pcr过程中随机引物和dna模版的不完全结合可能导致前几个碱基测序错误率
较高。
[0064]
为了保证信息分析质量，必须对raw reads过滤，得到clean reads，使用clean reads进行后续分析。使用软件fastp(version 0.20.0,参数：-n 10-q 20-u 40)对raw reads进行过滤，数据处理的步骤如下:
[0065]
3.5.1去除接头序列(adapter)；
[0066]
3.5.2当测序read中含有的n的含量超过该条read长度比例的10％时，需要去除此对paired reads；
[0067]
3.5.3当测序read中含有的低质量(q≤20)碱基数超过该条read长度比例的40％时，需要去除此对paired reads。
[0068]
得到检测数据后，将测序结果进行全长组装后，再次进行分析，可以采用聚类的方式，也可采用blast的方式，最终找到与测序结果亲缘关系最近的序列进而确定目标样本所属物种情况。
[0069]
上述7个测试物种共计10个样品利用针对哺乳动物的整套探针组合物(共3590条探针，每条长度110bp，具体探针序列见序列表序列1-序列3590)进行测试，测试结果中覆盖率为1(即100％)的表明待测样品为对应的物种，覆盖率不为1的表明待测样品并不是对应的物种，测试结果具体如下表1：
[0070]
表1
[0071]
[0072][0073][0074]
结果表明，10个样品检测出的物种和其实际所述物种完全相符合，准确率为100％。以序列1-序列3590的单链dna组成的探针组合物生产coi基因捕获试剂盒，可用于批
量样品的检测。
[0075]
以上对本发明进行了详述。对于本领域技术人员来说，在不脱离本发明的宗旨和范围，以及无需进行不必要的实验情况下，可在等同参数、浓度和条件下，在较宽范围内实施本发明。虽然本发明给出了特殊的实施例，应该理解为，可以对本发明作进一步的改进。总之，按本发明的原理，本技术欲包括任何变更、用途或对本发明的改进，包括脱离了本技术中已公开范围，而用本领域已知的常规技术进行的改变。按以下附带的权利要求的范围，可以进行一些基本特征的应用。
[0076]
工业应用
[0077]
本发明公开了鉴别或辅助鉴别哺乳动物物种的探针组合物及其试剂盒与应用。所述探针为序列表序列1-序列3590所示的核苷酸探针。利用本发明的探针捕获目标物种的coi基因，以此建库用二代测序高通量检测的方法获得coi基因的dna序列，可以同时兼顾一代测序灵活、方便、可以快速得到结果的优点，同时解决一代测序中存在的批量样本操作繁琐的问题，以及避免一代测序部分结果不稳定造成无法得到准确coi基因序列的风险。将探针开发成哺乳动物coi基因捕获试剂盒，用于大量哺乳动物物种dna条形码研究，可实现对未知哺乳动物物种的快速、批量鉴别。