监测家禽保种效果的SNP标记筛选方法及其在鸡保种上的应用、以及SNP标记的鉴定方法

申请号 CN201610601296.1 申请日 2016-07-28 公开(公告)号 CN106434867A 公开(公告)日 2017-02-22
申请人 江苏省家禽科学研究所; 发明人 韩威; 苏一军; 李国辉; 王洪志; 朱云芬; 盛中伟; 张会永; 殷建玫; 邹剑敏;
摘要 本 发明 涉及一种监测 家禽 保种效果的高 密度 SNP分子标记筛选方法及其在鸡种保种上的应用、以及SNP分子标记的鉴定方法。选取若干同类家禽品种,每个家禽品种选取若干个体,按照以下方法筛选SNP:1)所有家禽品种共有,各SNP位点在任一品种中至少存在两个等位基因;2)在基因组上分布相对集中,SNP位点数量在基因组1M 区间内 >90个;3)遗传多样性较为丰富,各SNP位点在任一品种中的最小等位基因 频率  >0.1;筛选得到满足上述1)、2)、3)要求的SNP位点。本发明提供一种利用简化基因组RAD-seq测序技术对不同地方鸡种进行测序,通过序列比对,鉴定出在不同鸡种基因组上集中分布、遗传多样性较为丰富的一致性SNP位点。
权利要求

1.监测家禽保种效果的高密度SNP分子标记筛选方法,其特征在于,选取若干同类家禽品种,每个家禽品种选取若干个体,按照以下方法筛选SNP:
1)所有家禽品种共有,各SNP位点在任一品种中至少存在两个等位基因;
2)在基因组上分布相对集中,SNP位点数量在基因组1M 区间内 >90个;
3)遗传多样性较为丰富,各SNP位点在任一品种中的最小等位基因频率 >0.1;
筛选得到满足上述1)、2)、3)要求的SNP位点。
2.根据权利要求1所述的监测家禽保种效果的高密度SNP分子标记筛选方法,其特征在于,应用于鸡品种保种效果监测中。
3.根据权利要求1所述的监测家禽保种效果的高密度SNP分子标记筛选方法,其特征在于,至少选用10个SNP位点进行检测。
4.根据权利要求2所述的监测家禽保种效果的高密度SNP分子标记筛选方法,其特征在于,所述鸡品种:选取鹿苑鸡、藏鸡、丝羽乌骨鸡、河南斗鸡、清远麻鸡、安义瓦灰鸡、大骨鸡、边鸡、寿光鸡、东乡绿壳蛋鸡、汶上芦花鸡、崇仁麻鸡、文昌鸡、金湖乌凤鸡、大围山微型鸡、惠阳胡须鸡、仙居鸡、茶花鸡、固始鸡、狼山鸡、瓢鸡、白黄鸡、北京油鸡、萧山鸡、琅琊鸡、淮南麻黄鸡、麻城绿壳蛋鸡、如皋黄鸡和隐性白羽肉鸡等29种不同鸡种。
5.根据权利要求4所述的监测鸡保种效果的高密度SNP分子标记筛选方法,其特征在于,每个鸡品种随机选择30个个体。
6.根据权利要求4或5所述的监测鸡保种效果的高密度SNP分子标记筛选方法,其特征在于,具体分布区间为 Chr4:88800000-89800000:SNP位点数量为108个;
Chr1:31100000-32100000:SNP位点数量为97个。
7.如上述任一权利要求所述的SNP分子标记的鉴定方法,其特征在于,按照以下步骤进行:
1)保种群个体血液采集:
使用医用一次性注射器采集保种群随机选留公、母鸡翅静脉全血0.5 1.0ml,采集后迅~
速注入含有2 µL 0. 5 mol/L Na2 EDTA 抗凝剂的无酶管中,然后将无酶管置于4℃环境下保存备用;
2)DNA提取和质量检测:
常温下吸出无酶管中保存备用的公、母鸡个体血液0.2 0.3ml,采用常规动物外周血~
苯-酚抽提法提取个体血液DNA;琼脂糖凝胶电泳分析DNA完整度,分光光度计检测DNA的纯度;
3)文库构建和简化基因组测序:
 提取后质检合格的基因组DNA 500 ng,加入0.6 U EcoRI、T4 DNA连接酶、ATP和EcoRI接头在37℃下反应3小时,65℃退火1小时,然后加限制性内切酶NlaIII和NlaIII接头在37℃下反应3小时,反应结束后在65℃ PCR仪中放置30分钟失活内切酶;
 使用琼脂糖凝胶电泳对连接产物进行片段选择,选择400-600bp回收酶切产物;
 使用Qubit3.0对回收产物进行DNA定量,等量混合样品;
 使用Illumina TruSeq试剂盒对混合产物进行DNA文库构建;
4)原始序列数据统计和过滤
为了保证信息分析质量,对测序得到的原始序列raw reads 过滤,得到clean reads用于后续分析;数据过滤的步骤如下: 去除带接头的paired-end reads; 当单端测序read 中含有的N 的含量超过该条read 长度比例的10%时,需要去除此对paired-end reads; 当单端测序read 中含有的低质量基数超过该条read 长度比例的 50% 时,需要去除此对paired-end reads;
5)序列数据比对和SNP鉴定
过滤后得到的clean reads有效测序数据比对到Gallus_gallus,WASHUC鸡参考基因组,根据比对结果,选取样品平均深度在5以上的group,定义ddRAD标记,进行基于样品的SNP鉴定。

说明书全文

监测家禽保种效果的SNP标记筛选方法及其在鸡保种上的应

用、以及SNP标记的鉴定方法

技术领域

[0001] 本发明涉及一种监测家禽保种效果的高密度SNP分子标记筛选方法及其在鸡种保种上的应用、以及SNP分子标记的鉴定方法,属于群体遗传学和分子生物学领域。

背景技术

[0002] 我国是世界上的鸡品种资源最为丰富的国家之一,现有地方鸡品种资源107个,保护好这些宝贵的遗传资源,对我国家禽产业乃至世界家禽产业可持续发展有重大意义。由于家禽特殊的卵生繁殖方式,活体保种是当前鸡品种资源保存的最主要形式。维持一定的保种群数量,通过采取科学合理的保种方法和继代繁育方式,有效减缓群体近交系数增量,保持其遗传多样性。近二十年来,我国地方鸡品种资源保护取得了显著成就,一是大部分地方鸡品种资源在原产地建有保种场或被国家级鸡种基因库引入保存,二是研究提出了多种适合不同条件的保种技术方法,包括家系等量留种随机选配法、家系等量留种轮回选配法、随机交配法等。然而,这些保种技术方法是基于群体遗传学理论,假定保种群处于理想状态下设定的,其具体效如何,需要经过保种实践的检验。因此,保种效果监测是鸡保种实施过程中的一项重要内容,依靠科学、可靠的保种效果监测与评价技术,能够及时发现保种实践中出现的各种问题,从而为进一步优化改进保种技术方案提供依据。
目前,家禽遗传资源保种效果监测主要依赖于:(1)常规表型性状记录比较分析。主要是选择一些常规易于检测的性状列入监测范围,并不能全面地反映出保种群遗传多样性在保种进程中的真实变化;(2)低密度DNA分子标记检测。发掘DNA分子遗传标记用于遗传多样性检测是近二十年来资源保护领域的研究热点,主要包括扩增片段长度多态性(AFLP)、限制性酶切片段长度多态性(RFLP)和微卫星多态(SSR)标记。这些分子标记数量非常有限(如普遍采用的适用于鸡遗传多样性分析的微卫星标记仅有约30个,且这些标记并不完全适用于遗传多样性丰富的中国地方鸡种),其所代表的基因组遗传变异信息量严重匮乏,导致监测效力不足。因此,发掘鉴定出适用于我国地方鸡种保种效果监测的高密度DNA分子标记显得尤为重要和迫切。高通量测序技术的发展为鸡遗传多样性保护提供了全新机遇。基于新一代测序技术(NGS),不仅将遗传变异信息的反映提升到单核苷酸多态(SNP)平,而且能在全基因组水平上高密度地检测出SNP标记。这使得鸡保种效果的全面精细监测成为可能,已经成为目前资源保护领域的一个研究热点。

发明内容

[0003] 本发明针对上述缺陷,目的在于提供一种利用简化基因组RAD-seq测序技术对不同地方鸡种进行测序,通过序列比对,鉴定出在不同鸡种基因组上集中分布、遗传多样性较为丰富的一致性SNP位点。为此本发明采用的技术方案是:本发明选取若干同类家禽品种,每个家禽品种选取若干个体,按照以下方法筛选SNP:
1)所有家禽品种共有,各SNP位点在任一品种中至少存在两个等位基因;
2)在基因组上分布相对集中,SNP位点数量在基因组1M 区间内 >90个;
3)遗传多样性较为丰富,各SNP位点在任一品种中的最小等位基因频率 >0.1;
筛选得到满足上述1)、2)、3)要求的SNP位点。
1)所有家禽品种共有,各SNP位点在任一品种中至少存在两个等位基因;
满足条件1的目的是确保筛选出的位点具有“通用性”,不受品种限制。要求各SNP位点在任一品种中至少存在两个等位基因(即等位基因频率不能为0或1),是确保根据每个SNP位点计算出的杂合度H、多态信息含量PIC不是0或1,若SNP位点的杂合度H、多态信息含量PIC为1或0就失去了保种效果监测的意义。
2)在基因组上分布相对集中,SNP位点数量在基因组1M 区间内 >90个;
满足条件2的目的是确保筛选出的染色体片段为“遗传多样性热区”,能更真实地反映出保种群世代传递过程中遗传多样性改变情况,进而反映出保种效果。
3)遗传多样性较为丰富,各SNP位点在任一品种中的最小等位基因频率 >0.1;
满足条件3的目的是确保实际应用中检测到的SNP位点信息真实、可靠,另外,减少不同世代计算结果间差异显著性分析误差。
表1 基于SNP位点等位基因频率的杂合度和多态信息含量计算
注:Pi、Pj为SNP位点第i和j个等位基因的频率。
(1)杂合度计算公式

(2)多态信息含量计算公式

上述方法应用于鸡品种保种效果监测中。
所述鸡品种:选取鹿苑鸡、藏鸡、丝羽乌骨鸡、河南斗鸡、清远麻鸡、安义瓦灰鸡、大骨鸡、边鸡、寿光鸡、东乡绿壳蛋鸡、汶上芦花鸡、崇仁麻鸡、文昌鸡、金湖乌凤鸡、大围山微型鸡、惠阳胡须鸡、仙居鸡、茶花鸡、固始鸡、狼山鸡、瓢鸡、白黄鸡、北京油鸡、萧山鸡、琅琊鸡、淮南麻黄鸡、麻城绿壳蛋鸡、如皋黄鸡和隐性白羽肉鸡等29种不同鸡种。
每个鸡品种随机选择30个个体。具体分布区间为Chr4:88800000-89800000:SNP位点数量为108个;Chr1:31100000-32100000:SNP位点数量为97个。
SNP分子标记的鉴定方法,其特征在于,按照以下步骤进行:
1)保种群个体血液采集:
使用医用一次性注射器采集保种群随机选留公、母鸡翅静脉全血0.5 1.0ml,采集后迅~
速注入含有2 µL 0. 5 mol/L Na2 EDTA 抗凝剂的无酶管中,然后将无酶管置于4℃环境下保存备用;
2)DNA提取和质量检测:
常温下吸出无酶管中保存备用的公、母鸡个体血液0.2 0.3ml,采用常规动物外周血~
苯-酚抽提法提取个体血液DNA;琼脂糖凝胶电泳分析DNA完整度,分光光度计检测DNA的纯度;
3)文库构建和简化基因组测序:
 提取后质检合格的基因组DNA 500 ng,加入0.6 U EcoRI、T4 DNA连接酶、ATP和EcoRI接头在37℃下反应3小时,65℃退火1小时,然后加限制性内切酶NlaIII和NlaIII接头在37℃下反应3小时,反应结束后在65℃ PCR仪中放置30分钟失活内切酶;
 使用琼脂糖凝胶电泳对连接产物进行片段选择,选择400-600bp回收酶切产物;
 使用Qubit3.0对回收产物进行DNA定量,等量混合样品;
 使用Illumina TruSeq试剂盒对混合产物进行DNA文库构建;
4)原始序列数据统计和过滤
为了保证信息分析质量,对测序得到的原始序列raw reads 过滤,得到clean reads用于后续分析;数据过滤的步骤如下: 去除带接头的paired-end reads; 当单端测序read 中含有的N 的含量超过该条read 长度比例的10%时,需要去除此对paired-end reads;
当单端测序read 中含有的低质量基数超过该条read 长度比例的 50% 时,需要去除此对paired-end reads;
5)序列数据比对和SNP鉴定
过滤后得到的clean reads有效测序数据比对到Gallus_gallus. WASHUC鸡参考基因组,根据比对结果,选取样品平均深度在5以上的group,定义ddRAD标记,进行基于样品的SNP鉴定。
本发明的优点是:本发明通过简化基因组RAD-seq高通量测序技术, 鉴定出在不同鸡种基因组上集中分布、遗传多样性较为丰富的一致性SNP位点,能显著提升鸡保种效果监测的科学性和准确性,也为进一步优化改进保种技术方案提供理论依据。
本发明利用高通量测序技术鉴定出在不同鸡种基因组上集中分布、遗传多样性较为丰富的一致性SNP位点,这些标记位点更适于地方鸡种保种效果监测;根据这些SNP位点信息可以建立保种群不同世代分子标记档案,能够显著提升鸡保种效果监测的科学性和准确性,也为进一步优化改进保种技术方案提供理论依据;在家禽资源保护领域有很好的应用前景。
附图说明
[0004] 图1为4个鸡品种不同SNP位点数量组合的杂合度变化图。图2至图6为本发明监测鸡保种效果的高密度SNP分子标记信息图。

具体实施方式

[0005] 本发明的目的是鉴定出适于监测鸡保种效果的高密度SNP分子标记及其应用。本发明的技术方案是:利用简化基因组RAD-seq测序技术对不同地方鸡种进行测序,通过序列比对,鉴定出在不同鸡种基因组上集中分布、遗传多样性较为丰富的一致性SNP位点。通过统计SNP位点等位基因频率,计算群体杂合度H、多态信息含量PIC等遗传参数,根据保种群世代传递过程中遗传参数变化情况监测保种效果。所采用的材料与方法以及所获得的结果如下所示:1 材料与方法
1.1实验材料
选取鹿苑鸡、藏鸡、丝羽乌骨鸡、河南斗鸡、清远麻鸡、安义瓦灰鸡、大骨鸡、边鸡、寿光鸡、东乡绿壳蛋鸡、汶上芦花鸡、崇仁麻鸡、文昌鸡、金湖乌凤鸡、大围山微型鸡、惠阳胡须鸡、仙居鸡、茶花鸡、固始鸡、狼山鸡、瓢鸡、白耳黄鸡、北京油鸡、萧山鸡、琅琊鸡、淮南麻黄鸡、麻城绿壳蛋鸡、如皋黄鸡和隐性白羽肉鸡等29种不同鸡种,所有鸡种均来自于江苏省家禽科学研究所国家级地方鸡种基因库保种群,每个鸡品种随机选择30个个体。
1.2保种群个体血液采集
使用医用一次性注射器采集保种群随机选留公、母鸡翅静脉全血0.5 1.0ml,采集后迅~
速注入含有2 µL 0. 5 mol/L Na2 EDTA 抗凝剂的无酶管中,然后将无酶管置于4℃环境下保存备用;
1.3 DNA提取和质量检测
常温下吸出无酶管中保存备用的公、母鸡个体血液0.2 0.3ml,采用常规动物外周血~
苯-酚抽提法提取个体血液DNA;琼脂糖凝胶电泳分析DNA完整度,分光光度计检测DNA的纯度;
1.4文库构建和简化基因组测序
 提取后质检合格的基因组DNA 500 ng,加入0.6 U EcoRI、T4 DNA连接酶、ATP和EcoRI接头在37℃下反应3小时,65℃退火1小时,然后加限制性内切酶NlaIII和NlaIII接头在37℃下反应3小时,反应结束后在65℃ PCR仪中放置30分钟失活内切酶;
 使用琼脂糖凝胶电泳对连接产物进行片段选择,选择400-600bp回收酶切产物;
 使用Qubit3.0对回收产物进行DNA定量,等量混合样品;
 使用Illumina TruSeq试剂盒对混合产物进行DNA文库构建;
1.5原始序列数据统计和过滤
为了保证信息分析质量,对测序得到的原始序列raw reads 过滤,得到clean reads用于后续分析;数据过滤的步骤如下: 去除带接头的paired-end reads; 当单端测序read 中含有的N 的含量超过该条read 长度比例的10%时,需要去除此对paired-end reads;
当单端测序read 中含有的低质量碱基数超过该条read 长度比例的 50% 时,需要去除此对paired-end reads;
1.6序列数据比对和SNP鉴定
过滤后得到的clean reads有效测序数据比对到Gallus_gallus. WASHUC鸡参考基因组,根据比对结果,选取样品平均深度在5以上的group,定义ddRAD标记,进行基于样品的SNP鉴定。
1.7 适于保种效果监测的SNP筛选
针对获得的存在于不同鸡种上的SNP位点,根据下述三个原则进行筛选: 所有鸡品种共有,各SNP位点在任一品种中至少存在两个等位基因; 在基因组上分布相对集中,SNP位点数量在基因组1M 区间内 >90个; 遗传多样性较为丰富,各SNP位点在任一品种中的最小等位基因频率 >0.1。筛选得到满足这些条件的SNP位点。
1.8 保种效果监测
通过统计筛选后的SNP位点等位基因频率,计算常规统计量(杂合度、多态信息含量等),根据保种群世代传递过程中常规统计量改变的差异显著性分析,监测保种效果。
2 结果
通过简化基因组测序,鉴定出了在不同鸡种基因组上集中分布、遗传多样性较为丰富的一致性SNP位点,具体分布区间为Chr4:88800000-89800000(SNP位点数量为108个),Chr1:31100000-32100000(SNP位点数量为97个),以上SNP位点分布及突变类型等详细信息见表1。通过统计上述SNP位点等位基因频率,计算群体遗传多样性度量参数,根据遗传参数在保种群世代传递过程中的变化情况监测保种效果。本专利结果能够显著提升鸡保种效果监测的科学性和准确性,也为进一步优化改进保种技术方案提供了科学依据,在家禽资源保护领域有很好的应用前景。
本发明筛选出了符合条件的200多个位点,理论上,在实际应用中使用的位点数量越多越好。可以是全部列出位点,也可以是其中的任意组合。为确定最小适用SNP位点数量,选择
4个鸡品种,每个鸡种随机选择30个个体(10个公鸡,20个母鸡),验证不同SNP位点数量对遗传多样性指标(杂合度H)度量精确度的影响。4个地方鸡品种不同SNP位点数量组合的杂合度见表2,结果表明不同的SNP数量会导致杂合度指标出现一定的波动,但SNP位点数量超过
10个,所测得的杂合度即趋于稳定,无显著变化(P>0.05)。
表2 4个鸡品种不同SNP位点数量组合的杂合度
QQ群二维码
意见反馈