CN103757002A

CN103757002A - 猪6号染色体上一个与产仔数相关的分子标记

Info

Publication number: CN103757002A
Application number: CN201310727885.0A
Authority: CN
Inventors: 张亚平; 黄路生; 谢海兵; 任军; 艾华水; 徐丹
Original assignee: Kunming Institute of Zoology of CAS; Jiangxi Agricultural University
Current assignee: Kunming Institute of Zoology of CAS; Jiangxi Agricultural University
Priority date: 2013-12-26
Filing date: 2013-12-26
Publication date: 2014-04-30
Anticipated expiration: 2033-12-26
Also published as: CN103757002B

Abstract

本发明涉及与猪产仔数性状相关的分子标记，具体涉及猪6号染色体上一个与产仔数相关的分子标记及其引物，属生物技术领域。该分子标记位点为基因组版本NCBI Build Sscrofa10.2的Chr6:103465181；该位点的等位基因为T和C，有C/C、T/C和T/T三种基因型，该位点及其侧翼序列如SEQ ID NO:1所示。扩增该分子标记位点的引物及延伸引物为：Chr6-3-PCRU：5’-AAA AAC AGT CGT GTG CCC T-3’；Chr6-3-PCRL：5’-ATA CTA TCT GCA TGT ATA GAG GAA CTT AAA-3’；Chr6-3-SNPU：5’-GAC TGA CTG ACT GAC TGA CTT TCA GTG TCT CAT GTA GAA TCA TCA-3’。带有C/C基因型的猪产仔数高于T/C和T/T基因型。本发明的有益效果在于：运用全基因组重测序方法，能够快速、批量地筛选与性状相关的分子标记或候选基因，Chr6:103465181这个位点能够作为遗传标记用于猪的遗传育种提高母猪的产仔数。

Description

猪6号染色体上一个与产仔数相关的分子标记

技术领域：

本发明涉及与猪产仔数性状相关的分子标记，具体涉及猪6号染色体上一个与产仔数相关的分子标记及其引物，属生物技术领域。

背景技术：

猪的产仔数性状是猪繁殖性状的主要组成部分，提高猪产仔数对提高养猪业的总体经济效益意义重大。但产仔数属于遗传力很低的数量性状，难以用常规育种技术进行遗传改良。我国的太湖猪以高产而闻名，窝产仔数比大多数其它中国地方猪种以及国外猪种的窝产仔数多4-5头（《中国猪品种志》张仲葛主编，上海科技出版社，1986）。二花脸猪是太湖猪的一种，是世界上繁殖力最高的猪种。太湖猪中与高产仔数相关的多态位点，能够作为遗传标记用于提高母猪的产仔数。

国内外研究人员采用全基因组扫描和候选基因的方法对产仔数性状进行QTL定位及候选基因的研究。Li等人(Li,K.,J.Ren,et al.Anim Genet.2009.40(6):963-966)利用覆盖了猪19条染色体的183个微卫星标记对与猪产仔数相关的总产仔数，产活仔数，产死仔数进行QTL（Quantitative Trait Locus,数量性状基因座）定位。他们研究了299个杜洛克猪与中国二花脸猪杂交的F₂代，在6，7，8和15号染色体上定位了与总产仔数，产活仔数和产死仔数相关的QTL。

Rothschild等(Rothischild,M.,Jacobson,D.et al.Proc Natl Acad Sci U SA.1996.93(1):201-205)通过候选基因法，利用包括中国梅山系太湖猪在内的两个极端品种杂交，发现1号染色体上雌激素受体座位的一个基因(ESR)与高产仔数主效基因连锁。通过雌激素受体基因座位的RFLPs分析，认为具有3.7kb条带的纯合型母猪(BB型)比具有4.3kb条带的纯合型母猪(AA型)初产时要多2.3头(P<0.01)，每胎平均多产1.5头(P<0.01)。目前已报道与产仔数相关的基因还有OPN(US6410227B1)，PRLR(Tomas,A.,J.Casellas,et al._J Anim Sci.2006.84(8):1991-1998)，FSHβ(Zhao,Y.Li N,et al.Sci China C Life Sci.1998.41(6):664-668)和RBP4(Spotter,A.,S.Muller,et al.Reprod Domest Anim.2009.44(1):100-105)等。

应用QTL定位方法能够定位与性状紧密连锁的区域，但是需要研究几代家系的很多个体，周期长，费时费力，在一般的实验室无法完成，而且有些QTL定位很宽泛。利用候选基因法对性状进行研究，相对来说简单直接，但是利用候选基因法主要研究的是外显子区域的突变，而很多发挥重要功能的突变很可能在内含子或者其它功能区域；有些复杂性状是受到多基因调控的；不能检出所有的QTL，而且有些用候选基因法研究的基因并没有找到QTL的定位。

随着技术的发展和成熟，基因组及基因注释越来越完善。由于测序技术的快速发展，测序成本的降低，能够利用全基因组测序或者转录组测序对性状相关基因或突变位点进行批量筛选。本发明选择二花脸猪（太湖猪的一种），藏猪，巴马香猪，莱芜猪，滇南小耳猪等5个中国地方猪种以及野猪进行了全基因组重测序。希望从中筛选出与高产仔数相关的SNPs（single nucleotide polymorphisms,单核苷酸多态性）位点作为遗传标记用于猪的遗传育种，以提高母猪的产仔数。

发明内容：

本发明的目的是提供猪6号染色体上一个与产仔数相关的分子标记及其引物。

本发明选择了二花脸猪（太湖猪的一种），藏猪，巴马香猪，莱芜猪，及滇南小耳猪5个中国地方猪种以及野猪进行了全基因组重测序。6个猪种的重测序个体数分别为11个，11个，9个，10个，10个和10个。重测序数据跟参考基因组（NCBI Build Sscrofa 10.2）进行比对，利用SAMTools软件进行SNPs位点的提取。发明人比较了二花脸猪与其它5个猪种，在6号染色体上发现一个SNP位点（基因组版本NCBI Build Sscrofa 10.2的Chr6:103465181）可能与猪的产仔数相关。该位点的等位基因为T和C，有C/C、T/C和T/T三种基因型。该位点及其侧翼序列（SEQ ID NO:1）为：

筛选到该位点后，发明人以猪的DNA为模板扩增包含该位点在内的核苷酸序列片段并进行测序，对该位点进行验证和相关性分析。其中发明人选择的个体是白色杜洛克×二花脸猪（太湖猪的一种）的F2后代。结果发现该位点与猪产仔数有显著相关性。

本发明的有益效果在于：运用全基因组重测序方法，能够快速、批量地筛选与性状相关的分子标记或候选基因，基因组版本NCBI Build Sscrofa 10.2的Chr6:103465181这个位点能够作为遗传标记用于猪的遗传育种提高母猪的产仔数。

附图说明：

图1：白色杜洛克×二花脸F₂代母猪资源群中Chr6:103465181位点3种基因型的产仔数比较（*表示两种基因型产仔数比较时，0.01≤P≤0.05；**表示P＜0.01）。

具体实施方案：

（1）SNP位点的群体分布特征

本发明选择了二花脸猪（太湖猪的一种），藏猪，巴马香猪，莱芜猪，及滇南小耳猪5个中国地方猪种以及野猪进行了全基因组重测序。6个猪种的重测序个体数分别为11个，11个，9个，10个，10个和10个。重测序数据跟参考基因组（NCBI Build Sscrofa 10.2）进行比对，利用SAMTools软件进行SNPs位点的提取。发明人比较了二花脸猪与其它5个猪种，在6号染色体上发现一个SNP位点（基因组版本NCBI Build Sscrofa 10.2的Chr6:103465181）可能与猪的产仔数相关。从重测序的分析结果来看，其它4个中国地方猪种和野猪中该位点等位基因T的频率平均为92%以上，而在二花脸猪中等位基因C的频率为100%。

（2）SNP位点的鉴定

在引物设计网站设计SNAPSHOT引物，包括一对PCR扩增引物（Chr6-3-PCRU：5’-AAA AAC AGT CGT GTG CCC T-3’；Chr6-3-PCRL：5’-ATACTA TCT GCA TGT ATA GAG GAA CTT AAA-3’）和延伸引物（Chr6-3-SNPU：5’-GAC TGA CTG ACT GAC TGA CTT TCA GTG TCT CAT GTA GAA TCATCA-3’）。该位点与其它11个位点共24条（12对）引物各取1μl混合为新的多重PCR引物池。首先进行12重PCR反应，反应体系所用试剂为：多重PCR引物池(6.25μM each)8μl、dNTP(10mM each)8μl、10×PCR Buffer II 116μl、MgCl₂(25mM Stock)234μl、Fast Start Taq(5U/μl)24μl和ddH₂O 460μl，将以上试剂混合液8.5μl和1.5μl DNA(50ng/μl)混合成10μl体系进行12重PCR扩增反应。扩增反应条件为：94℃ 4min；94℃ 30s，55℃ 30s，72℃ 1min，共40个循环。下一步进行PCR反应纯化，所需试剂为：Exo I(10U/μl)60μl、SAP(1U/μl)130μl、10×SAPBuffer 35μl和ddH₂O 125μl，在每个样品的PCR产物中加入3.5μl上述纯化混合液离心后按以下程序进行反应：37℃ 40min，96℃ 10min。然后进行单碱基延伸反应，所需试剂为：SNAPSHOT Multiplex ready reaction Mix 100μl、延伸引物池(6.25μM each，12条延伸引物各取1μl混合)8μl、5×Sequencing Buffer 200μl和ddH₂O 300μl，将以上试剂混合液6μl和4μl的纯化产物混合成10μl体系进行延伸反应。延伸反应条件为：96℃ 10s，50℃ 5s，60℃ 30s，共25个循环。接着进行延伸反应产物的纯化，所需试剂为：SAP(1U/μl)120μl、10×SAP Buffer 70μl和ddH₂O 140μl，在每个样品的延伸反应产物中加入上述3.3μl纯化混合液离心后按以下程序进行反应：37℃ 40min，75℃ 20min，4℃保存备用。每个样品取5-4μl上述纯化后的SNAPSHOT反应产物和5-6μl的高纯甲酰胺(Hi-Di Formamide)混合成10μl体系离心后上样于ABI 3730XL测序仪进行基因分型，使用GeneMarker软件分析分型结果。

(注：以上4个步骤的反应试剂混合液用量均为100个样品所需试剂量，从单碱基延伸反应开始，需要避光进行。)

（3）白色杜洛克×二花脸繁殖性能杂交资源群构建

后续的验证实验所用资源群个体样品及其性状数据由江西农业大学动物生物技术国家重点实验室培育基地提供。

1998年12月，在江西农业大学科研猪场将从江苏省3个二花脸猪保种场购买的二花脸母猪进行场间血缘的交叉配种。2001年1月根据亲本产仔性能，从3头公猪和11头母猪的纯繁后代中选择了17头二花脸母猪作为资源群体的祖代母本。把Sygen PIC公司惠赠的2头优质白色杜洛克公猪作为资源家系的祖代父本。将这2头白色杜洛克和17头二花脸猪杂交产生的F₁代个体进行配种，避免全同胞交配，且每头公猪通常与同一头母猪配种，以获得大样本的F₂半同胞家系。2003年1月、3月和6月，项目组分别将第1批和第2批F₂母猪送往江西省宜春市畜牧良种场（59头）、江西省国营红星种猪场（52头）和江西上犹县四元杂交猪场（118头）用于母猪繁殖性状的测定。

（4）SNP位点在杂交资源群F₂代母猪群体的基因分型

按照上述步骤(2)所述的方法，以192个白色杜洛克×二花脸猪杂交资源群F₂代母猪个体DNA为模板检测该SNP位点的基因型。F₂代资源群个体中187个个体成功分型，3种基因型C/C，T/C和T/T的个体数分别是39，88和60。

（5）SNP位点对产仔数性状的调控效应

在F₂代资源群个体中，带有C/C基因型个体的平均产仔数为11.5541，标准差为0.9051；带有T/C基因型个体的平均产仔数为10.8831，标准差为0.8526；带有T/T基因型个体的平均产仔数为10.3249，标准差为0.8692。

本发明采用一般线性模型（GLM）分析了SNP位点对猪总产仔数、产活仔数的影响。全部统计分析采用SAS软件进行。模型如下：

Y_ijklm=u+A_i+G_j+B_k+P_l+e_ijklm

其中Y_ijklm为猪总产仔数或产活仔数，u为群体均值，A_i为第i个个体的加性效应，G_j为第j个SNP位点基因型的固定效应，B_k为批次效应（k=1,2,3,4），P_l为随机效应（l=1,2,3）。

分析结果发现Chr6:103465181这个SNP位点跟产仔数有显著相关性（如图1所示）。带有C/C基因型的F₂代母猪的产仔数比带有T/T基因型的母猪平均多1.2292（P=0.0277）；虽然带有C/C基因型的F₂代母猪的产仔数与带有T/C基因型的母猪没有显著性差异，但是带有C/C基因型的F₂代母猪的产仔数比带有T/C基因型的母猪平均多0.671。从图1看出，带有T/T基因型的母猪产仔数最低。C为提高猪产仔数的有益等位基因。

从SNP位点与产仔数相关性分析看出，Chr6:103465181这个SNP位点能够作为分子标记应用于分子标记辅助选择（MAS）来提高母猪的产仔数。

SEQUENCE LISTING

<110> 中国科学院昆明动物研究所

江西农业大学

<120> 猪6号染色体上一个与猪产仔数相关的分子标记

<130> 7

<160> 1

<170> PatentIn version 3.5

<210> 1

<211> 601

<212> DNA

<213> Sus scrofa

<400> 1

agctgggtca tcagtcacaa agttatctgc atgtagagga tacttcatta ttgttctatg 60

gggctccaat taagtaggac atgataaaca gcgatcaagt tgacgatcca tccaaatatt 120

cagcaattac cacaacccat gaccatgaaa ggaacatggg cctctctgct tcctccacag 180

gtgaaggttg acaggaatga aagatttact tcaaacagcc tctgagcatc ttgctcctta 240

tcaagtctcc tcgaaaaaca gtcgtgtgcc ctcaattcag tgtctcatgt agaatcatca 300

tgagttaatg ctgattttaa gttcctctat acatgcagat agtatgacac taatgcaatg 360

atatcccatg aaaaaaagga tcatgatagg tttttattca ctttcagagc tgctaggaag 420

tcttcgtcat tgcatcagat tttaggccca agctcagaga gataaccacg tctaagtaca 480

ggacttagga gaagtacaaa tcattgcatt gattttcata ttgatgcaga catgcttagc 540

ctctaccaca acgagggaga tgcagcactc cacaggaatg ataggtgttt tcagatcttc 600

a 601

SEQUENCE LISTING

<110> 中国科学院昆明动物研究所

江西农业大学

<120> 猪6号染色体上一个与产仔数相关的分子标记

<130> 7

<160> 3

<170> PatentIn version 3.5

<210> 1

<211> 19

<212> DNA

<213> Sus scrofa

<400> 1

aaaaacagtc gtgtgccct 19

<210> 2

<211> 30

<212> DNA

<213> Sus scrofa

<400> 2

atactatctg catgtataga ggaacttaaa 30

<210> 3

<211> 45

<212> DNA

<213> Sus scrofa

<400> 3

gactgactga ctgactgact ttcagtgtct catgtagaat catca 45

Claims

1.猪6号染色体上一个与产仔数相关的分子标记，其特征在于：该分子标记位点为基因组版本NCBI Build Sscrofa 10.2的Chr6:103465181，该位点的等位基因为T和C，有C/C、T/C和T/T三种基因型。

2.权利要求1的分子标记，其中所述三种基因型特征在于：带有C/C基因型的猪产仔数高于带有T/C和T/T基因型的猪。

3.权利要求1的分子标记，其中所述等位基因C为提高猪产仔数的有益等位基因。

4.权利要求1的分子标记，其特征在于：该分子标记位点及其侧翼序列如SEQ ID NO:1所示。

5.权利要求1-4的分子标记，其中所述猪为太湖猪的后代。

6.以猪的DNA为模板扩增包含权利要求1-7所述分子标记在内的核苷酸序列片段。

7.以猪的DNA为模板对包含权利要求1-7所述分子标记在内的核苷酸序列片段进行测序。

8.扩增权利要求1所述分子标记位点的引物及延伸引物为：

Chr6-3-PCRU：5’-AAA AAC AGT CGT GTG CCC T-3’；

Chr6-3-PCRL：5’-ATA CTA TCT GCA TGT ATA GAG GAA CTT AAA-3’；

Chr6-3-SNPU：5’-GAC TGA CTG ACT GAC TGA CTT TCA GTG TCT CAT

GTA GAA TCA TCA-3’。