一种筛选种母鸡贮精能力性状候选位点的方法和应用

摘要

本发明公开了一种筛选种母鸡贮精能力性状候选位点的方法和应用，属于家禽分子标记筛选鉴定技术领域，本发明结合重测序技术和极端群体表型信息，通过计算每个位点的群体分化统计量并与全基因组关联分析相结合，对种母鸡贮精能力性状候选位点和候选基因进行筛选。通过人工构建表型梯度差异群体的方法，使实验样本的数目大幅度降低，有效缩减了研究成本；该筛选方法具有数据量丰富，检测精度高的优点。本发明基于该方法筛选得到了5个种母鸡贮精能力性状候选位点，从而有效提高了种母鸡遗传改良的效率，降低育种工作的成本，为实际的蛋鸡种禽育种工作提供理论依据。

说明书

技术领域

本发明属于家禽分子标记筛选鉴定技术领域，具体涉及一种筛选种母鸡贮精能力性状候选位点的方法和应用。

背景技术

鸡蛋是人们日常生活中必备的食物和消费品，我国是禽蛋消费和生产大国，2019年禽蛋产量为3308万吨且保持着持续增长的趋势，所以如何快速培育产蛋性能好、产受精蛋强的优质蛋鸡成为目前国内家禽种禽行业育种工作者亟待解决的问题。其中种母鸡贮精能力与高产受精蛋高度正相关，因此贮精能力性状的筛选是种禽(种母鸡)育种的思路之一。

“京红1号”蛋鸡作为我国自主培育的蛋鸡配套系，具有生产性能优越，繁殖性能突出等多项优势。“京红1号”蛋鸡的种母鸡繁殖性能受到高强度人工选择，其贮精能力在此过程中是否受到选择以及哪些基因控制着贮精能力的高低，这些问题都没有准确的答案。所以一些常用的筛选候选基因的方法(例如选择信号的方法)不适用于这类性状，而用大群体的重测序数据做全基因组关联分析(GWAS)的方法成本又比较高，基因芯片的方法也不够精确，因此通过人工构建具有表型梯度差异来计算不同群体的分化程度来与目标性状进行关联，并利用GWAS作为辅助方法进一步确定候选位点的方法是一个很好的突破点。虽然群体规模比较小会对GWAS结果造成一定影响，但是我们利用梯度的分化来与目标性状关联的方法在一定程度上弥补了这种影响，并且得到的结果相较于单一的GWAS更具准确性，因此这种挖掘蛋鸡种母鸡贮精能力候选位点和候选基因的方法对于蛋鸡种母鸡培育具有重要的理论意义和经济价值，并对肉鸡和地方品种鸡种母鸡的贮精能力性状筛选具有一定的参考和指导意义。

发明内容

本发明的目的在于提供一种筛选种母鸡贮精能力性状候选位点的方法和应用，该方法运用群体全基因组SNP位点标记的群体分化统计量与GWAS联合分析来评估蛋鸡种母鸡贮精能力性状候选位点和候选基因，具体为：结合重测序技术和极端群体表型信息，通过计算每个位点的群体分化统计量并与GWAS联合分析，用于挖掘影响蛋鸡种母鸡贮精能力性状的候选位点和基因，采用该方法可以有效提高蛋鸡种母鸡遗传改良的效率，降低育种工作的成本，为实际的蛋鸡种禽育种工作提供理论依据。

为实现上述目的，本发明采用的技术方案是：

一种筛选种母鸡贮精能力性状候选位点的方法，包括：

步骤1、饲养同一品种且日龄相同的产蛋种母鸡，从初产期开始每月统计每只种母鸡的最长受精天数表型数据，获得不少于250羽产蛋种母鸡的有效表型数据，并于统计的最后一天提取所述产蛋种母鸡的全部血样；

步骤2、对步骤1统计的最长受精天数表型数据进行分析，根据表型数据差异划分3对梯度亚群，包括最高差异梯度对、中等差异梯度对、普通差异梯度对群体，具体为：按最长受精天数进行秩排序，选取前50个样本与后50个样本组成普通差异梯度对，在此基础上取最长受精天数值最大的35个样本与最长受精天数值最小的35个样本组成中等差异梯度对，再在中等差异梯度对的基础上取最长受精天数值最大的20个样本与最长受精天数值最小的20个样本组成最高差异梯度对；

步骤3、提取步骤2所述的普通差异梯度对的100个产蛋种母鸡的血液样品DNA并进行重测序，去除缺失率大于20％的位点及最小等位基因频率小于0.05的位点，获得SNP位点数据，并基于所述SNP位点数据分别计算步骤2中三个梯度对的F

步骤4、对步骤3计算的F

步骤5、对步骤3获得的100个样品的重测序数据进行全基因组关联分析，取P值最小的10个位点作为候选位点；

步骤6、对步骤4获得的目标性状关联位点和步骤5获得的候选位点进行注释，取所述目标性状关联位点和所述候选位点的的重叠位点作为种母鸡贮精能力性状候选位点，所述贮精能力性状候选位点注释到的基因为贮精能力性状候选基因。

进一步的，步骤1中收集至少5个月的表型数据，并去除饲养过程中死亡或者表型缺失的产蛋种母鸡。

进一步的，所述方法还包括：将筛选得到的贮精能力性状的候选位点的基因型与表型数据进行关联分析，进行基因型效应估计。

本发明还提供了上述方法在筛选种母鸡贮精能力性状候选位点和/或种母鸡贮精能力性状候选基因中的应用。

本发明还提供了采用上述方法筛选得到的种母鸡贮精能力性状候选位点，包括：

SNP1位点：1号染色体的第65489499位点存在一个G＞A碱基突变；

SNP2位点：1号染色体的第87998885位点存在一个A＞G碱基突变；

SNP3位点：27号染色体的第4889482位点存在一个G＞A碱基突变；

SNP4位点：27号染色体的第4898241位点存在一个T＞A碱基突变；

SNP5位点：27号染色体的第4898269位点存在一个C＞T碱基突变。

进一步的，所述种母鸡贮精能力性状候选位点注释的贮精能力性状候选基因包括：ENSGALG00000052327基因、ENSGALG00000052718基因、ENSGALG00000052321基因、MRC2基因和TANC2基因。

进一步的，所述SNP1位点和SNP2位点的碱基突变与种母鸡贮精能力性状呈正相关。

进一步的，所述SNP3位点、SNP4位点和SNP5位点的碱基突变与种母鸡贮精能力性状呈负相关。

本发明还提供了上述种母鸡贮精能力性状候选位点在种母鸡贮精能力性状辅助选择中的应用。

与现有技术相比，本发明的有益效果是：

(1)本发明通过人工构建表型梯度差异群体的方法，使实验样本的数目相比常规的全基因组关联分析的方法大幅度降低，从而有效缩减了研究成本。

(2)本发明通过将群体分化与全基因组关联分析相结合的方法对种母鸡贮精能力性状候选位点和候选基因进行筛选，相比于传统的使用基因芯片技术的全基因组关联分析方法，该筛选方法具有数据量丰富，检测精度高的优点，且该方法适用范围广泛，包括但不限于已明确受到人工选择的表型性状。

(3)本发明基于群体分化与全基因组关联分析相结合的方法，筛选得到了5个种母鸡贮精能力性状候选位点及其注释的候选基因，可将其用于种母鸡贮精能力性状的辅助选择。

附图说明

图1为本发明实施例3中“京红1号”蛋鸡种母鸡贮精能力性状关联位点的分布图，其中等位基因频率差方法具有选择指示作用，大于0代表参考等位基因在贮精能力高的亚群受到选择，小于0代表突变的等位基因在贮精能力高的亚群受到选择；

图2为本发明实施例4中“京红1号”蛋鸡种母鸡贮精能力性状极端梯度群体连锁不平衡衰减图，其中high_3/low_3代表最高差异梯度对，high_2/low_2代表中等差异梯度对，high_1/low_1代表普通差异梯度对；

图3为本发明实施例4中“京红1号”种母鸡贮精能力与SNP关联的曼哈顿图及QQ图；

图4为本发明实施例5中位于“京红1号”种母鸡基因组27号染色体上三个候选位点间的连锁不平衡图。

具体实施方式

下面将结合本发明中的实施例，对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动条件下所获得的所有其它实施例，都属于本发明保护的范围。

以下实施例以“京红1号”蛋鸡的贮精能力性状为例用于说明本发明，但不用来限制本发明的范围。

实施例1表型测定

饲养同一品种的300只日龄相同的“京红1号”产蛋种母鸡，全部鸡在相同条件下，饲养于北京市华都峪口禽业有限责任公司祖代厂。从初产期每个月收集“京红1号”产蛋种母鸡的最长受精天数的表型数据，收集5个月的表型数据，去除饲养过程中死亡或者表型缺失的蛋鸡种母鸡个体，共获得251只种母鸡的有效表型数据，并在最后一次统计的同一天提取试验群京红1号”产蛋种母鸡的全部血样。

按照最长受精天数(DN)值对上述251只种母鸡进行秩排序，选取前50个样本与后50个样本组成普通差异梯度对，在此基础上取最长受精天数值最大的35个样本与最长受精天数值最小的35个样本组成中等差异梯度对，再在中等差异梯度对的基础上取最长受精天数值最大的20个样本与最长受精天数值最小的20个样本组成最高差异梯度对。在本实施例中根据表型变化划分梯度对，相邻梯度对的样本数目差值是15，即公差为15，表型描述性统计量如表1所示，其中H.1和L.1为普通差异梯度对，H.2和L.2为中等差异梯度对，H.3和L.3为最高差异梯度对。

表1最长受精天数(DN)分梯度表型数据

实施例2基因分型检测

(1)根据实施例1的结果，利用SDS-蛋白酶K消化法提取100只差异梯度群体样本的血液样品DNA，具体为：

a)室温解冻血样，取100μl血样于洁净的1.5ml离心管中，加入500μl的灭菌水，静置10min，8000r/min离心5min，弃上清；

b)加入500μl SET溶液重悬沉淀，用枪头将沉淀吹打破碎，8000r/min离心10min，弃上清；

c)重复上一步；

d)配制消化液(500μl SET+30μl SDS+10μl蛋白酶K)，每管加入540μl消化液，将沉淀吹打重悬；

e)将含带消化样品的离心管插入浮漂中，于55℃水浴锅中消化过夜(直到样品完全被消化)；

f)消化好的样品加5-20μl 10mg/ml的RNA酶，37℃水浴或烘箱中1.5-2h；

g)加入等体积的苯酚/氯仿/异戊醇(25：24：1)混合液560μl，缓慢颠倒离心管10次，静置15min，10000rpm室温离心15min，用黄枪头少量多次的吸上清至另一洁净离心管；

h)加入等体积的苯酚/氯仿/异戊醇(25：24：1)混合液及上述上清液，重复上述抽提步骤，将上清转移至1.5ml的EP管中；

i)在上清中加入1/10体积的3M NaAC和2-2.5倍体积预冷的无水乙醇，轻轻翻转至DNA聚集成团，呈白色丝状悬浮沉淀，用枪头直接挑出沉淀至新的EP管中；

j)加入500μl 75％乙醇洗涤沉淀，6000rpm离心5min，离心后倒掉残留液体，再重复洗涤一次，洗涤时彻底悬浮沉淀(悬浮时间与沉淀大小有关)，短暂离心，吸出残留液体，于室温下让酒精挥发干净；

k)加入20-100μl TE溶解后测定DNA的浓度，根据测定浓度稀释分装完毕后于-20℃或-80℃长期保存，原液避免反复冻融。

(2)使用全基因组重测序技术检测SNP位点及其基因型，并去除缺失值大于20％的位点，去除最小等位基因频率小于0.05的位点，共获得了1744595个高质量的SNP位点。

实施例3种母鸡贮精能力性状SNP位点群体分化统计量的计算

(1)根据实施例2获得的SNP位点，分别计算实施例1中所述的3个差异梯度对群体的F

F

其中p

AFD统计量是一种最直观体现群体分化程度的统计量，其计算公式为：

其中n代表在多态位点处观察到的不同等位基因的总数，f

(2)在上述步骤的基础上计算实施例1所述的3个梯度对中每个SNP位点的F

表2“京红1号”蛋鸡种母鸡贮精能力性状关联位点及注释基因

实施例4贮精能力性状与SNP位点之间的GWAS分析

(1)对实施例2得到的100个重测序数据进行全基因组关联分析(GWAS)，其作为辅助方法寻找与目标性状相关联的位点。

GWAS是在全基因组层面上通过对大规模的群体DNA样本进行全基因组高密度遗传标记(如SNP或CNV等)分型，从而寻找与目标性状相关的遗传因素的研究方法。使用GEMMA软件进行GWAS分析，常用混合线性模型，该模型如下：

y＝wα+xβ+Zu+ε

u～MVN

ε～MVN

其中n是个体数量，m是家系或群体的数量，y是n×1个数量性状向量，W＝(w1，w2...wc)是n×c矩阵包括1的列向量的协变量(固定效应)，α是包括截距在内的相应系数的c×1向量，x是标记基因型的n×1向量，β是标记的效应大小，Z是n×m荷载矩阵，u是随机效应的m×1向量，ε是误差的n×1向量，τ

(2)在上述步骤的基础上对100个个体的SNP重测序数据进行GWAS分析，在数据分析过程中，采用GEMMA软件进行无协变量的混合线性模型分析，应用simpleM方法确定全基因组显著临界值，当p＜0.05/M

表3基于GWAS鉴定的“京红1号”种母鸡贮精能力性状候选位点

(3)确定所有的候选SNP位点后，其中同时被F

表4“京红1号”种母鸡贮精能力性状显著关联位点梯度统计量结果

表5“京红1号”种母鸡贮精能力性状显著关联位点基因分型结果

实施例5基因型效应估计

通过将基因型和表型数据联系起来的方法，评估每个SNP性状关联位点的基因型效应(如表6所示)，结果显示，1号染色体的两个突变碱基的出现会使个体的表型值升高，表明突变碱基在这两个位点出现对种母鸡贮精能力是有益的影响(呈正相关)，而27号染色体的三个突变碱基的出现会使个体的表型值降低，表明突变碱基在这三个位点的出现会对贮精能力是负面的影响(呈负相关)。因为27号染色体的4889482位点是与种母鸡贮精能力相关性最显著的位点，该位点注释到的基因MRC2在文献中被报道与子宫内膜异位症的调节有关(Wei,C.,et al.,1-Methyl-tryptophan attenuates regulatory T cellsdifferentiation due to the inhibition of estrogen-IDO1-MRC2 axis inendometriosis.Cell Death Dis,2016.7(12):p.e2489.)，表明该结果作为候选位点具有可信度。我们通过haploview软件绘制了它与27号染色体的另外两个位点之间的连锁不平衡关系(如图4所示)，结果表明27号染色体上的这三个位点是强连锁关系，因此在选育种母鸡贮精能力高的个体的时候，应尽可能选择27号染色体4889482bp位点是GG、4898241bp位点是TT、4898269bp位点是CC的基因型，这为加快优质蛋鸡的育种进程提供帮助。

表6“京红1号”种母鸡最长受精天数的关联位点基因型效应估计

上述研究仅以“京红1号”蛋鸡种母鸡的贮精能力性状为实施例，并不用以限制本发明，对于二倍体生物的数量性状均适用，实施的方法均可按照实施例进行。

以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。