CN110760594B - 肉用西门塔尔牛6号染色体上与后腿长相关的snp位点及应用 - Google Patents

Abstract

本发明提供了肉用西门塔尔牛6号染色体上与后腿长相关的SNP位点及应用，所述SNP标记的位点为国际牛参考基因组UMD3.1版本6号染色体上第38696766个核苷酸位点，该位点的碱基是T或G。本发明通过优选该SNP的优势等位基因，能够逐代增加优势等位基因频率，提高肉用西门塔尔牛的后腿长，加快牛遗传改良进展从而有效提高肉牛育种的经济效益。

Description

肉用西门塔尔牛6号染色体上与后腿长相关的SNP位点及应用

技术领域

本发明涉及肉用西门塔尔牛6号染色体上与后腿长相关的SNP位点及其用途。

背景技术

肉用西门塔尔牛是中国最流行的牛杂交群体配套品种，具有生长速度快、产肉率高的特点，因此在育种生产中常作为终端父本使用。后腿长属于牛的胴体性状，一般而言，后腿长的肉牛体型更大，因此后腿部含肉量更多。牛基因组上控制后腿长度的基因有很多，常规育种手段难以快速加快其遗传进展。而分子标记技术的快速发展则为此类性状的精准选择提供了可能。全基因组关联分析是基于SNP间的连锁不平衡，通过统计分析策略来鉴别影响表型和基因型之间关系的分析方法，在鉴别影响牛重要经济性状的分子标记方面发挥了重要作用。本发明通过全基因组关联分析策略鉴别到了影响肉用西门塔尔牛后腿长的显著SNP，将其用于分子标记辅助选择和基因组选择中，选择对提高后腿长有利的基因型进行留种，从而逐代提高优势等位基因的基因频率，则能加快种牛育种改良的进程，为肉牛养殖带来更多经济效益。

发明内容

为了实现上述目标，本发明的首要目的在于提供牛6号染色体上与肉用西门塔尔牛后腿长相关的SNP位点，所述分子标记的核苷酸序列如SEQ ID NO：1所示，其中序列中的M是G或T，导致牛后腿长的不同。

所述分子标记位于肉用西门塔尔牛的6号染色体上的核苷酸序列上，所述分子标记的SNP位点为SEQ ID NO：1序列标注位置为214位的T214-G214的核苷酸突变；所述的分子标记的SNP位点对应于国际牛基因组UMD3.1版本参考序列6号染色体上第38696766位G＞T突变。

本发明的另一目的在于提供上述的分子标记在筛选高后腿长的牛个体的方法，具体为，检测牛6号染色体上权利要求1所述的分子标记，所述分子标记的5’端第214位单核苷酸是T还是G，淘汰T保留G。具体体现在所述牛选自内蒙古锡林郭勒盟乌拉盖管理区牧场的肉用西门塔尔牛资源群体。

本发明的另一目的在于提供一种用于鉴定上述影响牛后腿长的分子标记的引物对，所述引物对的核酸序列如下所示：

正向引物序列如SEQ ID NO：2所示；

反向引物序列如SEQ ID NO：3所示。

上述的引物对在鉴定影响牛后腿长中的应用。

上述的引物对在牛基因组选择中的应用。

上述的引物对在提高牛后腿长中的应用。

本发明的目的在于提供一种牛的遗传改良的方法，所述方法包括：确定种牛资源群体中的种牛的上述的影响牛后腿长的分子标记的位点，并根据所述分子标记做出相应的选择：在所述种牛资源群体中选择国际牛参考基因组UMD3.1版本6号染色体上第38696766位点为GG和TG、TT基因型的种牛个体，淘汰在第38696766该位点为TT基因型的牛个体，以逐代提高该位点的等位基因G的频率，从而提高后代牛的后腿长。

本发明相对于现有技术具有如下的优点及效果：

本发明研究并确定与牛后腿长相关的分子标记，验证其对后腿长的影响效应，最终建立高效准确的基因组选择育种技术，将其应用于种牛提高后腿长的遗传改良中，从而提高后代牛的后腿长，进而增加养殖企业市场竞争力。

附图说明

图1为肉用西门塔尔牛在6号染色体上关于后腿长的全基因组关联分析(GWAS)曼哈顿图；其中：横坐标表示牛的染色体编号；纵坐标表示-logP值。

图2为不同基因型肉用西门塔尔牛的后腿长。

具体实施方式

下面结合实施例及附图对本发明作进一步详细的描述，但本发明的实施方式不限于此。

本发明的上述发明目的具体是这样实现的：

实施例1

1、实验动物

本发明所使用的实验牛群群体均来自内蒙古锡林郭勒盟乌拉盖管理区牧场中国肉用西门塔尔牛1228头，为中国农业科学院北京畜牧兽医研究所牛遗传育种创新团队组建的肉用西门塔尔牛资源群体。

本实验共选取该资源群体中1228头肉用西门塔尔牛。肉用西门塔尔牛资源群体每年进行扩群，新增个体一般经历出生、育肥和屠宰3个阶段。每年3-5月出生的犊牛经过一段时间的放养管理后，在同年7月份由牛遗传育种创新团队进行统一初生重和体尺测量，并同时对基础母牛进行测量。同年10月份将5～9月龄的青年牛进行统一集中育肥，并收集生长发育性状的表型数据，同时左右用于Illumina BovineHD芯片基因分型获得基因型数据。当所有个体育肥期达10-12个月时，即来年的11月份左右，对所有肉用西门塔尔牛进行分批屠宰。屠宰过程严格按照《肉类采购规范》执行，屠宰数据、胴体数据和肉质数据严格按照GB/T27643-2011《屠宰后胴体性状和肉质性状测定指南》的要求进行测定。

2、样本采集

用采血管收集上述牛群体所有个体静脉血50ml，置于-80℃冰箱保存备用。

3、牛全基因组770K高密度芯片SNP判型

从上述资源群体中选取的1228头肉用西门塔尔牛中的每个个体采集静脉血50ml，用标准苯酚-氯仿法提取全基因组DNA，经Nanodrop2000/2000C核酸蛋白检测仪准确测定每一份样品的DNA的浓度和OD比值(OD260/280、OD260/230)。经Nanodrop2000/2000C核酸蛋白检测仪检测合格的DNA样品，按照检测的浓度将DNA稀释至50ng/μL左右。再将6μl已提取好的待测DNA样品与2μl Loading Buffer混合，上样到1％的琼脂糖凝胶中，150V电压下电泳25min，在紫外分光光度仪和凝胶成像设备下观察并拍照，观察DNA的完整性。

DNA样品送纽勤生物科技(上海)有限公司，根据公司标准流程进行牛全基因组Illumina BovineHD芯片770K SNP芯片(Illumina，美国)基因型判定。利用PLINK v1.90软件对所有样本770K芯片扫描分型数据进行质量控制，剔除检出个体率低于90％、家系孟德尔错误率高于0.1、最小等位基因频率小于0.05且哈代-温伯格平衡显著性水平高于10^-6的SNP，最终得到671,204个SNP的有效基因型数据。

4、全基因组关联(GWAS)分析

为了消除群体层化效应，本发明采用线性混合模型单点回归分析并结合R语言GenABEL软件包进行GWAS分析，分析模型中利用个体间基因组的相似度校正层化效应。采用Bonferrini法确定SNP与后腿长关联程度的显著性阈值，基因组水平显著阈值为0.05除以有效SNP位点数量，即基因组显著水平阈值为7.45e-8，即0.05/671,204(有效SNP数量)；染色体水平显著阈值为1除以有效SNP位点数量，即染色体显著水平阈值为1.49e-6，即1/671,204(有效SNP数量)。

GWAS分析结果如图1所示。从图1可知，在肉用西门塔尔牛6号染色体中存在显著影响后腿长的位点，最强关联的SNP为g.214T＞G(P＝9.45E-10)。

5、不同基因型与后腿长表型的关联性分析

根据表1可知，分子标记的SNP位点g.214T＞G与后腿长极显著相关(P＜0.001)，说明此分子标记显著影响牛的后腿长，可以通过对牛的此SNP位点的辅助选择，从而提高该群体的后腿长，进而加快目标性状的育种进程。

另外根据表1可知，GG型、TG型比TT型的牛后腿更长，说明TT型牛个体对筛选后腿长不利，故优先保留GG、TG型的种牛。因此，在育种过程中需要逐步淘汰TT型的种牛，优先保留GT、TG型的种牛，以逐代提高该位点的等位基因G的频率。

表1分子标记的SNP位点g.214T＞G与后腿长的相关性

6、目的DNG序列扩增及测序

(1)引物设计

通过Ensembl网站(http：//GsiG.ensembl.org/index.html)下载牛的6号染色体上SEQ ID NO：1的DNA序列。并利用引物设计软件primer premier 6.0设计引物。

设计的引物的DNA序列如下所示：

P001正向：5’-GCTGCTGCTAAGTCGCTTCA-3’，

P002反向：5’-TCATCAGACCTTTCCCACTCCA-3’；

(2)PCR扩增

10uL的反应体系中加入DNA模板1uL，双蒸水3.4uL，2×TGg PCR StanMix withLoading Dye 5uL，引物P001和P002各0.3ul。PCR反应条件为：94℃预变性5min后，94℃变性30s、56℃退火30s、72℃延伸45s，35个循环，最后72℃延伸5min。

(3)DNA序列测定

DNA序列测序鉴定：在北京生工科技有限公司进行，基因片段测正反两个反应。将所测得的序列与NCBI基因组序列对比，得出对应SNP位点的突变。测序结果如下所示：

注：序列表中标注的M为突变位点，用标有下划线显示(括号中为突变碱基，为等位基因突变)，在该序列的首尾加粗显示为设计引物序列位置。

7、分子标记的SNP位点g.214T＞G效应分析

通过对分子标记辅助选择，对群体内基因型为TT的牛进行淘汰，可显著提高群体的后腿长，为企业带来更多经济效益。

本发明通过对SEQ ID NO：1序列中的第214个碱基突变位点进行检测，初步进行其基因型与牛的后腿长之间的关联分析的应用，为牛的分子标记辅助选择及基因组选择提供了一个新的分子标记。

上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。

序列表

<110> 中国农业科学院北京畜牧兽医研究所

<120> 肉用西门塔尔牛6号染色体上与后腿长相关的SNP位点及应用

<160> 3

<170> SIPOSequenceListing 1.0

<210> 1

<211> 255

<212> DNA

<213> 人工序列(Artificial Sequence)

<400> 1

gctgctgcta agtcgcttca gtcgtgtccg actctgtgcg accccataga tggcagccca 60

ccaggccccg ccatccctgg gattctccag gcaagaacac tggagtgggt tgtcatttcc 120

tcctccaatg catgaaagtg aaaagtgaaa gtgaagtcgc tcagtcgtgt ccgactctta 180

gcaaccacat ggaccaggct cctccatcca tgggattttc caggcaagag tactggagtg 240

ggaaaggtct gatga 255

<210> 2

<211> 20

<212> DNA

<213> 人工序列(Artificial Sequence)

<400> 2

gctgctgcta agtcgcttca 20

<210> 3

<211> 22

<212> DNA

<213> 人工序列(Artificial Sequence)

<400> 3

tcatcagacc tttcccactc ca 22

Claims (4)

1.一种提高肉用西门塔尔牛后腿长的方法，其特征在于，所述方法包括以下步骤：

检测牛的国际参考基因组UMD3 .1版本6号染色体上第38696766个核苷酸位点的基因型，选择第38696766个核苷酸位点的GG、TG型个体作为种牛。

2.根据权利要求1所述的方法，其特征在于，所述检测牛的国际牛参考基因组UMD3 .1版本6号染色体上第38696766个核苷酸位点的基因型的方法包括以下步骤：

(1)提取待测牛的基因组DNA；

(2)采用引物对，将所述待测牛的基因组DNA进行PCR扩增，以便获得PCR扩增产物；

(3)对所述PCR扩增产物进行测序，以便获得测序结果；

(4)基于所述测序结果，确定所述待测牛的SNP标记的基因型；

所述引物对的核酸序列如SEQ ID NO：2和SEQ ID NO：3所示；

所述SNP标记的序列如SEQ ID NO：1所示，所述SEQ ID NO：1所示序列自5’端起第214位碱基是T或G。

3.一种SNP标记在提高肉用西门塔尔牛后腿长中的应用，其特征在于，所述的SNP标记的序列如SEQ ID NO：1所示，所述SEQ ID NO：1所示序列自5’端起第214位碱基是T或G。

4.一种引物对在提高肉用西门塔尔牛后腿长中的用途，其特征在于，所述引物对的核酸序列如SEQ ID NO：2和SEQ ID NO：3所示。

Images