一种基于全基因组str建立动物亲权鉴定系统的方法
【专利摘要】本发明涉及种质资源保护领域，特别涉及动物亲权鉴定领域，公开了一种基于全基因组STR建立动物亲权鉴定系统的方法，包括：1)全基因组STR位点查找；2)STR位点的初步筛选；3)引物设计；4)动物群体随机抽样并提取基因组DNA；5)STR位点群体多态性验证；6)对具有遗传多态性的STR位点进行哈迪-温伯格平衡检验，筛选出基因型频率分布满足哈迪-温伯格平衡的STR位点；7)STR位点的法医学学应用价值评价；8)建立STR亲权鉴定系统。该方法具有操作简便的优点，所构建亲权鉴定系统的非父排除率可达0.99以上，能准确进行动物亲权鉴定从而可用于指导珍稀物种的繁育及经济类物种的鉴别及育种。
【专利说明】一种基于全基因组STR建立动物亲权鉴定系统的方法
【技术领域】
[0001] 本发明涉及一种亲权鉴定系统的建立方法，特别涉及一种基于全基因组序列中 STR位点建立用于动物亲权鉴定系统的方法。
【背景技术】
[0002] 无论是对珍稀动物种质资源的保护，还是对具有较高经济价值动物的繁育，要求 对动物个体识别与亲权鉴定的事件越来越多。例如，对大熊猫进行亲子鉴定和个体识别，有 助于避免人工饲养繁殖中近亲繁殖的难题，对拯救大熊猫物种有重要意义。又如，在现代育 马业及赛马业中，对一些优秀品种的马匹进行出生来源鉴定也备受重视。
[0003] DNA分子标记由于能直接反映遗传物质本身或基因组的变化差异等，且大多数呈 中性突变，遗传较为稳定，符合孟德尔遗传规律，不受生理期和环境的影响，并在群体中具 有高度多态性和较高的杂合度等，因此在动物个体识别和亲权鉴定方面有重要的应用前 景。DNA分子标记主要包括限制性片段长度多态性、随机扩增DNA多态性、DNA指纹图、小卫 星DNA、微卫星DNA标记和单核苷酸多态性分析等。
[0004] 微卫星DNA，又称短串联重复序列（Short Tandem Repeat，STR)，是以1?6个喊基 为核心而组成的短串联重复序列，广泛随机地分布于真核生物基因组中，在动物基因组DNA 序列中，平均6?IOkb就可能出现一个，其一般大小为100?350碱基。微卫星DNA比其 他标记更具有检测方法简便、快速，分析材料广泛，重复性好，结果可靠等优点，因此被认为 是进行个体识别与亲权鉴定的绝佳分子标记。
[0005] 亲权鉴定的基本原理有两点：
[0006] 1，在肯定子代的某个标记基因是来自生父，而假设父亲并不带有这个基因的情况 下，可以排除它是该子代的父亲。由此可知，被检查的血型系统或DNA分子标记个数越多， 假设父亲被排除的几率越大；
[0007] 2,在肯定子代的某些标记基因是来自生父，而假定父亲也带有这些基因的情况 下，不能排除它是该子代的生父。这时可以计算出它是该子代生父的概率，在使用多个遗 传标记时，通常计算其累积排除概率：假设若干相互独立的标记个数排除概率分别为P n P2……Pn，则η个相互独立的标记个数的累积排除概率为：
[0008] Pc= (ι-p D (I-P2)……（I-Pn)
[0009] 在动物育种及家畜问题引起的民事纠纷中，STR分型技术已经被用于多种动物的 亲子鉴定，但是由于传统方法在开发和利用微卫星DNA标记的同时，必须首先通过克隆STR 位点，测序得到其两侧特异的保守序列，并以此设计引物。该方法整个过程周期长，并且筛 选通量低，从而给STR的应用带来了一定的困难和障碍。目前随着测序技术的发展和测序 通量的不断提高，越来越多的全新物种基因组得到解码，可以通过STR查找软件非常方便 的查找STR位点及其两侧序列。基于此，本发明开发了从基因组参考序列出发，建立动物 STR亲权鉴定系统的方法。
[0010] 经过文献检索，未见与本发明相同的公开报道。
【发明内容】
[0011] 本发明的目的在于提供一种在动物全基因组序列中筛选STR位点进而建立STR亲 权鉴定系统的方法。
[0012] 本发明建立了一套基于全基因组参考序列信息建立该物种STR亲权鉴定系统的 整套方法。依据本发明中的方法，能够快速建立目标物种的STR亲权鉴定系统，该系统能够 有效鉴定动物个体间的亲子关系及亲缘关系，因此能够为动物育种提供有力指导。
[0013] 一种基于全基因组STR建立动物亲权鉴定系统的方法，包括以下步骤：
[0014] 1)全基因组STR位点查找：通过串联重复序列查找软件Tandem R印eat Finder(TRF)在所研宄物种全基因组寻找所有STR位点；
[0015] 2)STR位点的初步筛选：统计出该物种全基因组参考序列中重复单位为4个碱基 的STR位点，参考该物种群体遗传多态性高低情况，并从中随机选取50?300个STR位点；
[0016] 3)引物设计：对于步骤2)中所选取的每一个STR位点，在串联重复序列前后两端 200碱基长度的保守序列内，分别设计用于多核苷酸扩增反应的正向引物和反向引物；
[0017] 4)动物群体随机抽样并提取基因组DNA :在所研宄的动物群体中通过随机抽样的 方式获得样本，并提取样本中所有个体的基因组DNA ;
[0018] 5) STR位点群体多态性验证：
[0019] ①使用步骤3)中的引物分别扩增样本中所有个体的基因组DNA目标基因组区 域；
[0020] ②扩增产物经分离后进行基因分型；
[0021] ③对于具有遗传多态性的STR位点，估算各等位基因频率；
[0022] 6)采用皮尔森卡方检验对步骤5)中各STR位点的基因型分布频率进行哈迪-温 伯格平衡检验，筛选出符合哈迪-温伯格平衡的STR位点；
[0023] ①计算基因型分布频率的观察值；
[0024] ②计算基因型分布频率的理论值；
[0025] ③用皮尔森卡方检验判断每个STR位点基因型频率分布的观察值和理论值是否 具有显著差异，如果二者具有显著差异P〈〇. 05,则该位点不符合哈迪-温伯格平衡，在后 续步骤中将该位点剔除；如果二者不具有显著差异P>〇. 05,则该位点符合哈迪-温伯格平 衡；
[0026] 7) STR位点的法医学学应用价值评价：对步骤6)中筛选出的各STR位点分别计算 其法医学参数，如杂合度、个体识别率、多态信息量、非父排除率等；
[0027] 8)建立STR亲权鉴定系统：综合以上参数，选取法医学应用价值高的若干独立位 点，使得累积排除概率大于〇. 99,由这些位点及对应的引物组成该物种的STR亲权鉴定系 统。
[0028] 本发明的有益效果在于：
[0029] 1)从全基因参考序列中筛选STR位点，与传统分子克隆方法相比，该方法周期短， 筛选通量高；具有方便快速，成本低廉的优点。
[0030] 2)该方法也可用于其他重复单元长度的STR。
[0031] 3)该方法也可用于建立动物的STR个体识别系统。
[0032] 4)该方法可以根据鉴定的精度需求，增加 STR检测位点。
[0033] 该方法具有操作简便的优点，所构建亲权鉴定系统的非父排除率可达0. 99以上， 能准确进行动物亲权鉴定从而可用于指导珍稀物种的繁育及经济类物种的鉴别及育种。
【专利附图】
【附图说明】
[0034] 图1是用STR亲权鉴定系统鉴定的一个四代家系家谱，以及家系中所有个体在各 STR位点的基因型。
【具体实施方式】
[0035] 下面结合具体实施例对本发明做进一步详细说明。
[0036] 实施例：朱鹮STR亲权鉴定系统的建立
[0037] 具体步骤：
[0038] 1)定义重复单元为1?6bp并且总长度不少于15bp的串联重复序列为STR。
[0039] 用串联重复序列查找软件Tandem R印eat Finder(TRF)在朱鹮基因组中寻找STR。 查找程序参数设定如下："Match = 2，Mismatch = 7，Delta = 7，PM = 80，PI = 10，Minscore =30，Maxperiod = 6"。在朱鹤基因组中共发现至少82万个STR位点，其中4碱基重复的 STR位点数约5万个。
[0040] 2)从重复单元为4bp的STR位点中选取300个STR位点。
[0041] 选取原则：i )该位点为非简并性STR位点；ii )各位点独立遗传。
[0042] 3)使用引物设计软件Primer Premier 5. 0在STR两端200bp以内的保守序列中 分别设计上游和下游引物。设计引物参数设定如下：引物长度18?22bp，产物长度160? 400bp〇
[0043] 4)在陕西洋县4个不同的朱鹮亚群中，随机抽取105个个体构成样本，其中包括 48个雄性和57个雌性。采用血液DNA抽提试剂盒E. Z. N. A. ? Blood DNA Kit从朱鹮血液 中提取基因组DNA。所得的DNA样本用试剂盒TIANamp Genemic DNA Kit进行定量。所有 步骤按照试剂盒说明书操作。
[0044] 5)对于步骤2)中选取的位点，在步骤4)中抽取的样本中进行遗传多态性验证。 具体操作方法包括：
[0045] ①通过多核苷酸扩增反应方法对基因组中包含重复序列的目标区域进行扩增。反 应在GeneAmp 9700多核苷酸扩增反应仪上进行。多核苷酸扩增反应程序：
[0046] 95°C /5 分钟
[0047] 95 °C /30 秒-52 °C /30 秒-72 °C /30 秒，30 个循环
[0048] 72°C /10 分钟
[0049] 对上述扩增产物进行聚丙烯酰胺凝胶电泳，根据扩增产物片段大小挑选出等位基 因数目最高的50个STR位点，对于长度重叠的位点使用不同的荧光标记，第一组用FAM(蓝 色）标记，第二组用HEX (绿色）标记，第三组用TMR (黄色）标记。
[0050] ②扩增产物通过 ABI3730DNA Genetic Analyzer48_capillary array system 分 离后，根据扩增产物片段大小用软件Genemapperf. 5进行基因分型；
[0051] ③以片段大小作为等位基因，计算各位点等位基因频率，即计算各等位基因数量 占该基因座全部等位基因总数的比率。
[0052] 6)采用皮尔森卡方检验对各STR位点的基因型分布频率进行哈迪-温伯格平衡检 验。基因型频率是指群体中某特定基因型个体数占全部个体数的比率，哈迪-温伯格平衡 是指基因型频率分布的观察值和理论值无显著差异（P>〇. 05)。
[0053] ①计算基因型分布频率的观察值：观察到的特定基因型个体数除以全部个体数。
[0054] ②计算基因型分布频率的理论值：对于一个STR位点，假设有η个等位基因，用Xi 代表该第i个等位基因的频率，则该位点的基因型分布频率的理论值计算公式为：
[0055]
【权利要求】
1. 一种基于全基因组STR建立动物亲权鉴定系统的方法，其特征在于，包括以下步骤： 1) 全基因组STR位点查找：通过串联重复序列查找软件TandemR印eatFinder(TRF) 在所研宄物种全基因组参考序列中寻找所有STR位点； 2. STR位点的初步筛选：统计出该物种全基因组参考序列中重复单位为4个碱基的STR 位点，并从中随机选取50?300个STR位点； 3) 引物设计：对于步骤2)中所选取的每一个STR位点，在串联重复序列前后两端200 碱基长度的保守序列内，分别设计用于多核苷酸扩增反应的正向引物和反向引物； 4) 动物群体随机抽样并提取基因组DNA:在所研宄的动物群体中通过随机抽样的方式 获得样本，并提取样本中所有个体的基因组DNA; 5. STR位点群体多态性验证： ① 使用步骤3)中的引物分别扩增样本中所有个体的基因组DNA目标基因组区域； ② 扩增产物经分离后进行基因分型； ③ 以片段大小作为等位基因，计算各STR位点等位基因分布频率； 6) 采用皮尔森卡方检验对步骤5)中各STR位点的基因型分布频率进行哈迪-温伯格 平衡检验，筛选出符合哈迪-温伯格平衡的STR位点； 7)STR位点的法医学学应用价值评价：对步骤6)中筛选出的符合哈迪-温伯格平衡的 STR位点分别计算其法医学参数； 8) 建立STR亲权鉴定系统：根据法医学参数，选取若干独立位点，使得累积排除概率大 于0. 99,由这些位点及对应的引物组成该物种的STR亲权鉴定系统。
2. 根据权利要求1所述的基于全基因组STR建立动物亲权鉴定系统的方法，其特征在 于，所述步骤5)-②中，扩增产物经分离后根据片段大小，利用软件GenemappeU. 5进行基 因分型。
3. 根据权利要求1所述的基于全基因组STR建立动物亲权鉴定系统的方法，其特征在 于，所述步骤5)-③中，计算各位点等位基因频率，即计算各等位基因数量占该基因座全部 等位基因总数的比率。
4. 根据权利要求1所述的基于全基因组STR建立动物亲权鉴定系统的方法，其特征在 于，步骤6)中，采用皮尔森卡方检验对各STR位点的基因型分布频率进行哈迪-温伯格平 衡检验，哈迪-温伯格平衡是指基因型频率分布的观察值和理论值无显著差异，具体为： ① 计算基因型分布频率的观察值：观察到的特定基因型个体数除以全部个体数； ② 计算基因型分布频率的理论值：对于一个STR位点，假设有n个等位基因，用Xi代表 该第i个等位基因的频率，则该位点的基因型分布频率的理论值计算公式为：
③ 用皮尔森卡方检验判断每个STR位点基因型频率分布的观察值和理论值是否具有 显著差异，如果二者具有显著差异P〈〇. 05,则该位点不符合哈迪-温伯格平衡，在后续步骤 中将该位点剔除；如果二者不具有显著差异P>〇. 05,则该位点符合哈迪-温伯格平衡。
5. 根据权利要求4所述的基于全基因组STR建立动物亲权鉴定系统的方法，其特征在 于，基因型频率分布的观察值是观察到的特定基因型个体数除以全部个体数。
6.根据权利要求1所述的基于全基因组STR建立动物亲权鉴定系统的方法，其特征 在于，步骤7)中，分别计算其法医学参数包括杂合度、个体识别率、多态信息量和非父排除 率。
【文档编号】C12Q1/68GK104480205SQ201410757128
【公开日】2015年4月1日 申请日期:2014年12月10日 优先权日:2014年12月10日
【发明者】李生斌, 李波, 刘清波, 常辽, 熊子军 申请人:西安交通大学