CN113151501A - 一种黄牛wbp1l基因cnv标记辅助检测生长性状的方法及其应用 - Google Patents

Abstract

本发明公开了一种黄牛WBP1L基因CNV标记辅助检测生长性状的方法及其应用：基于实时荧光定量PCR技术，以黄牛个体基因组DNA为模板，以一对特异性引物扩增WBP1L基因拷贝数变异区域的部分片段，以另外一对特异性引物扩增BTF3基因的部分片段作为内参，最后利用‑ΔΔCt判定个体的拷贝数变异类型。本发明提供的方法基于WBP1L基因拷贝数变异与生长性状之间的关联，可用于快速建立遗传资源优势种群，有利于加快黄牛分子标记辅助选择育种工作，该方法简单、快速，便于推广应用。

Description

一种黄牛WBP1L基因CNV标记辅助检测生长性状的方法及其 应用

技术领域

本发明属于分子遗传育种领域，具体涉及一种检测牛WBP1L基因拷贝数变异的方法，该方法利用基因组DNA进行实时荧光定量PCR，以BTF3基因为参照、根据-ΔΔCt确定个体的拷贝数变异类型。

背景技术

基因组DNA携带生物体的所有遗传信息，研究证明，遗传信息的改变是相关重要经济性状变异的主要来源，随着高通量测序技术的日渐成熟，筛选及鉴定经济性状相关基因，以及发掘和精准定位重要表型相关位点成为推动分子育种发展的一项重要任务。

目前分子育种主要的研究方向集中在标记辅助选择(molecular mark-assistselection,MAS)技术上，该技术是通过DNA分子标记对遗传资源或育种材料进行选择，从而实现对畜禽品种的经济性状进行改良。在畜禽育种中，通过对与生长性状密切相关，并且与数量性状紧密关联的DNA分子标记的选择，达到早期选种和提高育种值准确性的目的，从而在畜禽育种中获得更大的遗传进展。

拷贝数变异(Copy Number Variations，CNVs)为基因结构变异，是指1kb到数MbDNA片段的拷贝数重复、缺失等类型的突变。CNVs是由基因重排导致的，是基因组变异的重要组成部分，可导致由基因剂量效应、基因断裂、基因融合和位置效应等引起的表型多态性。研究表明，拷贝数变异可能影响基因网络并调节相关基因的表达，促成个体表型，因此，检测与黄牛生长性状相关基因的CNV标记有助于加速黄牛遗传选育的进程。

在检测已知CNV的各种方法中，实时荧光定量PCR(quantitive Real-Time PCR,qPCR)是使用比较广泛的一种技术。该方法操作简单、敏感性高、速度快，通过对目的基因(有拷贝数变异)和内参基因(无拷贝数变异)进行相对定量，然后利用-ΔΔCt判定个体的拷贝数变异类型。

WBP1L(WW域结合蛋白1Like)是一种跨膜适配蛋白编码基因。近年来国内外对WBP1L基因的研究较少，对WBP1L基因的功能及作用机制的认识都还处在初步阶段。目前为止，尚未见到关于WBP1L基因拷贝数变异与黄牛生长性状关联的报道。

发明内容

本发明的目的在于提供一种黄牛WBP1L基因CNV标记辅助检测生长性状的方法及其应用，从而加快黄牛良种选育进程。

为达到上述目的，本发明采用了以下技术方案：

一种检测牛WBP1L基因拷贝数变异的方法，包括以下步骤：以牛个体基因组DNA为模板，以引物对P1及引物对P2为引物，分别通过实时荧光定量PCR扩增WBP1L基因的拷贝数变异区域及作为内参的BTF3基因的部分片段，然后根据定量结果鉴定个体WBP1L基因的拷贝数变异类型；

所述引物对P1为：

上游引物F1：5′-CCAGAGCTGGACTTCGTGGG-3′

下游引物R1：5′-TCGTGACATATCTCAGACGCAG-3′；

所述引物对P2为：

上游引物F2：5′-AACCAGGAGAAACTCGCCAA-3′

下游引物R2：5′-TTCGGTGAAATGCCCTCTCG-3′。

优选的，所述WBP1L基因的拷贝数变异区域位于牛WBP1L基因候选区域Chr26:23641732-23643331(参考序列为AC_000183.1)。

优选的，所述个体WBP1L基因的拷贝数变异类型是根据-ΔΔCt将定量结果分为的三类：重复型(Duplication类型)，-ΔΔCt>0.5，并可按拷贝数(CN)再进行细分，例如CN＝3或CN≥4的多种拷贝数不同的重复型；缺失型(Deletion类型)，-ΔΔCt<-0.5，并可按拷贝数(CN)再进行细分，例如CN＝0或CN＝1的拷贝数不同的缺失型；正常型(Normal类型)，-0.5≤-ΔΔCt≤0.5，按拷贝数(CN)对应为CN＝2。

优选的，所述实时荧光定量PCR所用的扩增反应体系包括10～50ng/μL模板DNA 1μL以及10μmol/L的引物对P1或引物对P2所对应的上、下游引物各0.5μL。

优选的，所述实时荧光定量PCR所用的反应程序包括以下步骤：95℃预变性10min；95℃变性15s，60℃退火1min，共39个循环。

优选的，基于引物对P1扩增的PCR产物片段大小为145bp，基于引物对P2扩增的PCR产物片段大小为166bp。

上述检测牛WBP1L基因拷贝数变异的方法在黄牛分子标记辅助选择育种中的应用。

优选的，在皮南牛群体中，Deletion类型(CN＝0)的个体在胸围、体重、尻长性状方面的表现均优于其他类型；在郏县红牛群体中，Duplication类型(CN＝3)的个体在体高、体长、尻长、胸深、体重性状方面的表现均优于其他类型；在固原牛群体中，Deletion类型(CN＝1)的个体在体高性状方面的表现均优于其他类型；在吉安牛群体中，Normal类型(CN＝2)的个体在体高性状方面的表现均优于其他类型。

一种检测牛WBP1L基因拷贝数变异的试剂盒，该试剂盒包括以上引物对P1及引物对P2。

本发明的有益效果体现在：

本发明以BTF3基因作为内参，采用实时荧光定量PCR对牛WBP1L基因的拷贝数变异进行检测，可以对个体的WBP1L基因拷贝数变异类型进行鉴定，从而依据所发现的定位于WBP1L基因的DNA分子标记(CNV标记)辅助检测具有生长性状优势的个体，快速建立牛遗传资源优势种群，从而加快牛(例如，地方黄牛品种)的育种进程。

本发明具有以下优点：

1.本发明提供的WBP1L基因拷贝数变异检测方法，不受个体年龄和性别的限制，可用于黄牛早期选育。

2.与高通量测序方法、基因芯片方法相比，本发明快速、简单、成本低，能够准确的鉴定个体的拷贝数变异类型。

3.本发明在一定程度上为黄牛生长性状的分子标记辅助选择提供科学依据。

附图说明

图1为本发明实施例中进行qPCR(WBP1L基因)绘制的扩增曲线。

图2为本发明实施例中进行qPCR(WBP1L基因)绘制的溶解曲线。

具体实施方式

下面结合附图和实施例对发明做进一步详细说明，所述实施例是对本发明的解释，而不是对本发明保护范围的限制。

本发明对WBP1L基因的拷贝数变异进行检测并用于黄牛分子育种，主要包括以下步骤：参考NCBI数据库的牛WBP1L基因序列，利用Primer-BLAST网站进行引物设计，并用普通PCR对引物进行检验；采用实时荧光定量PCR技术检测候选位点(Chr26:23641732-23643331)在群体中的拷贝数变异情况；将拷贝数变异类型与生长性状进行关联分析，筛选到与生长性状相关的CNV标记；根据拷贝数变异类型筛选生长性状优异的个体，并建立用于选育的牛种群。具体说明如下。

1.黄牛样本采集

本发明以秦川牛、皮南牛、云岭牛、夏南牛、郏县红牛、固原牛、吉安牛为检测对象，采集和使用了732头黄牛个体的血液样本(见表1)，并记录个体生长性状数据，如体高、体重、体长、胸围、胸深、尻长等，以用于后续的关联分析。

表1.样本采集信息

2.血样基因组DNA的提取

(1)冷冻血样(主要为血细胞)室温解冻，吸取500μL血液于1.5mL离心管中，加入等体积的磷酸缓冲液(PBS)混匀，温和摇动，4℃、12000r/min离心5min，弃去上清液；重复本步骤至上清液透明、沉淀呈透明色。

(2)在离心管中加入DNA抽提缓冲液500μL，轻轻吹打，使血细胞沉淀脱离离心管壁，37℃水浴1h。

(3)加蛋白酶K至5μL(20mg/mL)并混匀，55℃水浴中消化过夜(16h左右)至絮状沉淀不见，尚未澄清的，可补加10μL蛋白酶K混匀后继续消化直至澄清，将反应液冷却至室温。

(4)加入500μL Tris饱和酚，温和摇动15min，使其充分混匀，4℃、12000r/min离心10min，将上层水相转入另一灭菌离心管；重复本步骤1次。

(5)加入氯仿500μL，温和摇动20min，使其充分混匀，4℃、12000r/min离心15min，将上层水相转入另一灭菌的1.5mL离心管。

(6)加入氯仿、异戊醇混合液(24:1)500μL，充分混合20min，4℃、12000r/min离心10min，将上清液转入另一1.5mL离心管中。

(7)加入0.1倍体积的NaAc缓冲液及2倍体积的冰冷乙醇，转动离心管直至白色絮状沉淀析出。

(8)4℃、12000r/min离心10min，弃去上清液，用70％的冰冷乙醇漂洗DNA沉淀2次。

(9)4℃、12000r/min离心10min，弃去上清液，室温下使乙醇挥发干净。

(10)干燥后的DNA加入80～100μL的TE溶解，4℃保存直至DNA完全溶解，利用紫外分光光度仪检测其质量，-80℃保存。

3.目的基因及内参基因特异性引物设计

以NCBI数据库(http://www.ncbi.nlm.nih.gov/)公布的牛WBP1L基因(AC_000183.1)为参考序列，查找重测序中筛选出来的WBP1L基因(目的基因)拷贝数变异候选区域的序列(Chr26:23641732-23643331)，利用Primer-BLAST网站设计用于扩增该区域内靠近中间位置的一段145bp序列的引物(引物对P1)，同时，以NCBI所公布的牛BTF3基因序列(AC_000177.1)为参考序列，采用相同的方法设计用于扩增内参基因(BTF3基因)中一段166bp序列的引物(引物对P2)。引物对序列信息如表2所示(引物合成时间2020年10月)。

表2.实时荧光定量PCR的引物信息

注：F1或F2为上游引物，R1或R2为下游引物

4.实时荧光定量PCR

qPCR反应体系如表3所示。

表3.qPCR的反应体系

qPCR的扩增反应程序为：

(1)95℃预变性10min，然后按照(2)进行扩增反应；

(2)95℃变性15s，60℃退火1min，共39个循环。

参见图1及图2，通过绘制扩增曲线和溶解峰确定引物适用于实时荧光定量PCR分析。

5.个体CNV类型判定

根据目的基因和内参基因的qPCR实验结果计算样本的-ΔΔCt的值：

ΔΔCt＝ΔCt_(实验组)-ΔCt_(参照组)

ΔCt_(实验组)＝Ct_{(实验组目的基因)}-Ct_{(实验组内参基因)}，ΔCt_(参照组)＝Ct_{(参照组目的基因)}-Ct_{(参照组内参基因)}

其中，实验组即为待检测有无拷贝数变异的个体样本，参照组即为已知无拷贝数变异的个体样本，可以采用重测序试验中所选择的参照组各品种黄牛个体。Ct即Cyclethreshold，为PCR扩增过程中扩增产物的荧光信号达到设定的阈值时所经过的扩增循环次数。

根据以上公式计算得出每个个体样本的-ΔΔCt，并根据以下判定标准鉴定个体的CNV类型：-ΔΔCt>0.5为Duplication类型(重复型)；-0.5≤-ΔΔCt≤0.5为Normal类型(正常型)，-ΔΔCt<-0.5为Deletion类型(缺失型)。

在划分出三种不同CNV类型的同时，分析每个个体的拷贝数：因利用qPCR已得出不同个体的-ΔΔCt的值，所以接下来首先计算出2^-ΔΔCt的值，然后计算2×2^-ΔΔCt的值，最后将所得结果取整处理，便可分别得到不同个体的拷贝数。

6.WBP1L基因CNV位点与生长性状的关联分析

生长性状：体高、体长、十字部高、胸围、腰角宽、尻长、体重、坐骨端宽、荐高、胸深、胸宽、管围。

CNV类型：Deletion、Normal和Duplication。

关联分析模型：利用SPSS(23.0)进行关联性分析。在数据处理中，根据影响生长性状指标的因素的不同，考虑到环境效应、年龄、品种、遗传效应及其互作效应，采用固定模型进行分析，同时根据实际情况进行简化。最终的模型如下：

Z_ijk＝μ+M_i+N_j+CNV_k+e_ijk

其中，Z_ijk为个体表型记录；μ为群体均值；M_i为第i个年龄效应；N_j为第j个品种效应；CNV_k为第k个WBP1L基因拷贝数变异类型效应；e_ijk为随机误差。

关联分析结果如表4、表5、表6及表7所示。

表4.皮南牛WBP1L基因拷贝数变异与生长性状之间的关联分析

注：CN表示拷贝数(Copy Number)；平均值肩标上具有相同字母表示差异不显著(P>0.05)，平均值肩标上字母不同表示差异显著(P<0.05)，^*P<0.05，^**P<0.01；n表示具有相同拷贝数的个体样本数。

从表4可以看出，皮南牛的WBP1L基因存在与胸围、体重、尻长这三个生长性状的优势密切相关的CNV标记，即Deletion类型(CN＝0)。

表5.郏县红牛WBP1L基因拷贝数变异与生长性状之间的关联分析

注：CN表示拷贝数(Copy Number)；平均值肩标上具有相同字母表示差异不显著(P>0.05)，平均值肩标上字母不同表示差异显著(P<0.05)，^*P<0.05；n表示具有相同拷贝数的个体样本数。

从表5可以看出，郏县红牛的WBP1L基因存在与体高、体长、尻长、胸深、体重这五个生长性状的优势密切相关的CNV标记，即Duplication类型(CN＝3)。

表6.固原牛WBP1L基因拷贝数变异与生长性状之间的关联分析

注：CN表示拷贝数(Copy Number)；平均值肩标上具有相同字母表示差异不显著(P>0.05)，平均值肩标上字母不同表示差异显著(P<0.05)，^*P<0.05；n表示具有相同拷贝数的个体样本数。

从表6可以看出，固原牛的WBP1L基因存在与体高这个生长性状的优势密切相关的CNV标记，即Deletion类型(CN＝1)。

表7.吉安牛WBP1L基因拷贝数变异与生长性状之间的关联分析

注：CN表示拷贝数(Copy Number)；平均值肩标上具有相同字母表示差异不显著(P>0.05)，平均值肩标上字母不同表示差异显著(P<0.05)，^*P<0.05；n表示具有相同拷贝数的个体样本数。

从表7可以看出，吉安牛的WBP1L基因存在与体高这个生长性状的优势密切相关的CNV标记，即Normal类型(CN＝2)。

总之，关联分析结果表明(见表4至表7)，WBP1L基因的拷贝数变异位点与牛体高、胸围、体重等重要生长性状具有显著的关联性。因此，WBP1L基因的拷贝数变异位点(Chr26:23641732-23643331)可以作为一个有效提高牛生长性状的候选分子遗传标记(CNV标记)位点，从而加快黄牛选育进程。

<110> 西北农林科技大学

<120> 一种黄牛WBP1L基因CNV标记辅助检测生长性状的方法及其应用

<160> 4

<210> 1

<211> 20

<212> DNA

<213> F1

<400> 1

ccagagctgg acttcgtggg 20

<210> 2

<211> 22

<212> DNA

<213> R1

<400> 2

tcgtgacata tctcagacgc ag 22

<210> 3

<211> 20

<212> DNA

<213> F2

<400> 3

aaccaggaga aactcgccaa 20

<210> 4

<211> 20

<212> DNA

<213> R2

<400> 4

ttcggtgaaa tgccctctcg 20

Claims (10)

1.一种检测牛WBP1L基因拷贝数变异的方法，其特征在于：包括以下步骤：

以牛基因组DNA为模板，以引物对P1及引物对P2为引物，分别通过实时荧光定量PCR扩增WBP1L基因的拷贝数变异区域及作为内参的BTF3基因的部分片段，然后根据定量结果鉴定WBP1L基因的拷贝数变异类型；

所述引物对P1为：

上游引物F1：5′-CCAGAGCTGGACTTCGTGGG-3′

下游引物R1：5′-TCGTGACATATCTCAGACGCAG-3′；

所述引物对P2为：

上游引物F2：5′-AACCAGGAGAAACTCGCCAA-3′

下游引物R2：5′-TTCGGTGAAATGCCCTCTCG-3′。

2.根据权利要求1所述一种检测牛WBP1L基因拷贝数变异的方法，其特征在于：所述WBP1L基因的拷贝数变异区域位于牛WBP1L基因候选区域Chr26:23641732-23643331。

3.根据权利要求1所述一种检测牛WBP1L基因拷贝数变异的方法，其特征在于：所述拷贝数变异类型是根据-ΔΔCt将定量结果分为的三类：重复型，-ΔΔCt>0.5；缺失型，-ΔΔCt<-0.5；正常型，-0.5≤-ΔΔCt≤0.5。

4.根据权利要求1所述一种检测牛WBP1L基因拷贝数变异的方法，其特征在于：所述实时荧光定量PCR的扩增反应体系包括10～50ng/μL模板DNA 1μL以及10μmol/L的引物对P1或引物对P2所对应的上、下游引物各0.5μL。

5.根据权利要求1所述一种检测牛WBP1L基因拷贝数变异的方法，其特征在于：所述实时荧光定量PCR的反应程序包括以下步骤：95℃预变性10min；95℃变性15s，60℃退火1min，共39个循环。

6.根据权利要求1所述一种检测牛WBP1L基因拷贝数变异的方法，其特征在于：基于引物对P1扩增的PCR产物片段大小为145bp，基于引物对P2扩增的PCR产物片段大小为166bp。

7.如权利要求1-6中任意一项权利要求所述的方法在黄牛分子标记辅助选择育种中的应用。

8.根据权利要求7所述的应用，其特征在于：在皮南牛中，拷贝数变异类型为缺失型的个体在生长性状上较优；在郏县红牛中，拷贝数变异类型为重复型、且拷贝数接近正常型的个体在生长性状上较优；在固原牛中，拷贝数变异类型为缺失型、且拷贝数接近正常型的个体在生长性状上较优；在吉安牛中，拷贝数变异类型为正常型的个体在生长性状上较优。

9.根据权利要求8所述的应用，其特征在于：所述生长性状选自胸围、体重、尻长、体高、体长、胸深中的一种或多种。

10.一种检测牛WBP1L基因拷贝数变异的试剂盒，其特征在于：该试剂盒包括用于分别通过实时荧光定量PCR扩增牛WBP1L基因拷贝数变异区域及作为内参的BTF3基因的部分片段的引物对P1及引物对P2；

所述引物对P1为：

上游引物F1：5′-CCAGAGCTGGACTTCGTGGG-3′

下游引物R1：5′-TCGTGACATATCTCAGACGCAG-3′；

所述引物对P2为：

上游引物F2：5′-AACCAGGAGAAACTCGCCAA-3′

下游引物R2：5′-TTCGGTGAAATGCCCTCTCG-3′。

Images