牦牛全基因组SNP位点的应用及检测用引物组和试剂盒

摘要

本发明公开了提供牦牛全基因组的16个SNP位点在牦牛表型性状或分子育种分析中的应用，解决了现有技术中尚未有采用SNP位点对牦牛表型性状或分子育种分析的问题。本发明的16个SNP位点的核苷酸序列如SEQ ID No.1‑16所示。本发明还公开了用于检测该16个SNP位点的引物组及包括该引物组的试剂盒。本发明又公开了牦牛表型性状或分子育种分析的方法。本发明的16个SNP位点可为育种上下游过程提供支持，具有测定准确率高、易实现标准化和自动化检测等诸多优势。

说明书

技术领域

本发明属于畜牧技术领域，具体地涉及牦牛全基因组SNP位点的应用及检测用引物组和试剂盒。

背景技术

牦牛是分布于以青藏高原为中心，及其毗邻的高山、亚高山地区的特有牛种，能充分利用高寒草地的牧草资源，对高寒草地的生态环境条件具有极强的适应性。在空气稀薄、牧草生长期短、寒冷、枯草期长的恶劣环境条件下生活自如，繁衍后代，为当地牧民提供奶、肉、毛、役力、燃料等生产生活必需品，是当地畜牧业经济中不可或缺的，可称之为“全能”的家畜。牦牛在遗传资源上是一个极为宝贵的基因库。然而，中国牦牛数量和品种类群多，分布广，各牦牛品种产肉产奶，体型特征等表型各不相同。而中国地方牦牛品种基因组变异信息的不完整，是导致我国牦牛品种的育种及相关研究局限的重要原因之一。

单核苷酸多态性(single nucleotide polymorphism，SNP),主要是指在基因组水平上由单个核苷酸的变异所引起的DNA序列多态性。它是生物可遗传的变异中最常见的一种。SNP是公认的第三代遗传标志，许多表型差异、对疾病的易感性等等都可能与SNP有关。且从实验操作来看，通过SNP发现表型相关基因突变要比通过家系来得容易；有些SNP并不直接导致相关性状，但由于它与某些相关基因相邻，可成为重要的标记。目前并未见关于跟牦牛体重相关的SNP标记位点的报道。

因此，提供牦牛表型相关SNP的位点信息，用于在牦牛表型性状或牦牛分子育种分析中，成为了本领域技术人员亟待解决的问题。

发明内容

本发明的目的之一在于，提供牦牛全基因组的16个SNP位点在牦牛表型性状或分子育种分析中的应用，以解决现有技术中尚未有采用SNP位点对牦牛表型性状或分子育种分析的问题。

本发明的目的之二在于用于检测该16个SNP位点的引物组。

本发明的目的之三在于提供包括该引物组的试剂盒。

本发明的目的之四在于提供牦牛表型性状或分子育种分析的方法。

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

本发明所述的牦牛全基因组的16个SNP位点在牦牛表型性状或分子育种分析中的应用，所述16个SNP位点如下：

AX-174702570，核苷酸序列如SEQ ID No.1所示；

AX-174961896，核苷酸序列如SEQ ID No.2所示；

AX-174407967，核苷酸序列如SEQ ID No.3所示；

AX-174402854，核苷酸序列如SEQ ID No.4所示；

AX-174929694，核苷酸序列如SEQ ID No.5所示；

AX-174547362，核苷酸序列如SEQ ID No.6所示；

AX-174734142，核苷酸序列如SEQ ID No.7所示；

AX-174706158，核苷酸序列如SEQ ID No.8所示；

AX-174783962，核苷酸序列如SEQ ID No.9所示；

AX-174627015，核苷酸序列如SEQ ID No.10所示；

AX-174928167，核苷酸序列如SEQ ID No.11所示；

AX-174555047，核苷酸序列如SEQ ID No.12所示；

AX-174845027，核苷酸序列如SEQ ID No.13所示；

AX-174891371，核苷酸序列如SEQ ID No.14所示；

AX-174570649，核苷酸序列如SEQ ID No.15所示；

AX-174620133，核苷酸序列如SEQ ID No.16所示。

本发明所述的16个SNP位点的引物组，该引物组包括上引物和下引物，其核苷酸序列如Seq ID.No.17-48所示。

包括如上所述的引物组的试剂盒。

本发明所述的牦牛表型性状或分子育种分析的方法，包括以下步骤：

步骤1：合成如上所述的16个SNP位点对应的探针序列；

步骤2：提取牦牛的DNA，采用如上所述的引物组进行PCR扩增，从扩增条带判断检测牦牛是否符合高体重的育种需求。

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

本发明的16个SNP位点来自于32个中国地方牦牛品种特有遗传变异信息，可以应用于性状关联分析和分子辅助育种等方面，可为育种上下游过程提供支持。具有测定准确率高、易实现标准化和自动化检测等诸多优势。

具体实施方式

下面结合实施例对本发明作进一步说明，本发明的方式包括但不仅限于以下实施例。

本发明的思路在于，先对全国范围内32个不同牦牛品种共96个牦牛个体进行全基因组重测序数据。该32个牦牛品种涵盖了我国境内所有的牦牛种类。从96个牦牛的SNP突变中挑选高可信度的位点来制作芯片。然后再扩大样本，用268个牦牛个体，利用制备出的芯片挑选出了16个跟产肉相关的SNP位点。

实施例1

本实施例提供了从96个牦牛的重测序数据中挑选SNP突变可信度最高位点的方法，制备芯片。

首先利用Illumina X Ten测序平台对来自全中国范围内的32个不同地方牦牛群体(每个群体3头个体共计96个牦牛个体)进行全基因组重测序数据(每个个体10X覆盖度)，得到2592G的测序数据。在此基础上，首先进行一个质量控制，过滤掉测序质量较低的数据。接下来我们用BWA比对软件将各个总计的reads比对到参考基因组上，生成sam格式的比对结果文件，然后利用samtools将sam格式的结果文件转换成bam格式并进行排序，最后我们利用鉴定SNP最普遍的GATK来鉴定出96个牦牛中的SNP突变。最后我们从96个牦牛的SNP突变中挑选了一些高可信度的位点来制作芯片。

通过与牦牛现有数据库的比较，鉴定、验证和筛选出与牦牛遗传性状最相关的位点，最终通过Affymetrix公司设计定做出了包含630209个SNP位点的高密度牦牛SNP基因分型芯片。

96个牦牛重测序数据中SNP突变筛选及芯片制备的具体步骤如下：

1.提取DNA。

2.DNA样本检测。

对DNA样品的检测主要包括3种方法：(1)琼脂糖凝胶电泳分析DNA降解程度以及是否有RNA及蛋白等污染。(2)Nanodrop法检测DNA的纯度(OD260/280比值)。(3)Qubit法对DNA浓度进行精确定量。根据上述检测结果，采用OD值在1.8～2.0之间，含量在1.5μg以上的DNA样品用来建库。上述琼脂糖凝胶电泳、Nanodrop法、Qubit法均为现有技术。

3.数据分析

获得原始测序序列(Sequenced Reads)后参考基因组(BosGru_v2.0)进行信息分析流程，大致包括以下两个部分：

1)测序数据质量评估：主要对数据量、碱基质量、比对率、覆盖率、捕获率、均一性等指标进行统计，评估建库测序是否达到了标准，符合标准则进行后续分析。

2)变异检测：将高质量的序列比对到牦牛参考基因组上，检测样本中的变异信息，并对检出的变异进行统计和注释。

4.芯片制备

将重测序挑选出来的SNP信息提供给Affymetrix公司，通过Affymetrix公司平台的位点筛选标准最终选择了630209个SNP位点进行芯片设计及制作。

实施例2

本实施例提供了268个牦牛产肉相关性状SNP位点的筛选方法。

选取了3个国家认证的牦牛群体，包括嘉黎牦牛、帕里牦牛和斯布牦牛总计268个个体，通过实施例1制备的芯片进行SNP位点的筛选。样本严格筛选各群体中4-9岁处于青壮年期的牦牛个体，每个群体都筛选体重偏低和体重偏高的健康个体入组作为芯片筛选的样本。对每个个体的体重数据做详细的记录，根据得到的芯片筛选数据进行GWAS分析，最终找到16个跟牦牛产肉最相关的位点。

具体步骤为：

样品gDNA利用NanoDrop ND-2000(Thermo Scientific)定量并经凝胶电泳检测DNA完整性。DNA质检合格后，样本的扩增、片段化、沉淀、重悬、芯片的杂交以及洗染参照芯片标准流程。gDNA重悬后加入杂交MIX，然后进行质检。质检合格后，利用GeneTitan MCInstrument进行芯片杂交，洗染和扫描。

数据分析部分

将GeneTitan系统扫描得到的原始数据导入到软件Axiom Analysis Suite中，利用该软件对原始数据进行聚类和基因分型。最终通过Axiom AnalysisSuite软件导出PLINK格式数据用于后续的数据分析。

筛选出来的16个SNP位点信息如下表所示：

实施例3

本实施例提供了用于检测实施例1筛选出的16个SNP位点的引物组，所述引物组的核苷酸序列与SNP位点的对应关系见下表：

实施例4

本实施例提供了16个SNP位点在牦牛表型性状和分子辅助育种中的应用方法，具体为：

合成这16个位点对应的探针序列，采集牦牛耳样，提取DNA，质控后采用本发明的引物组进行PCR扩增，从扩增条带判断检测牦牛是否符合高体重的育种需求。其中探针序列的合成，牦牛DNA的提取，以及PCR扩增均为现有技术。如果根据引物扩增出来的条带突变位点为“+”标记，则证明该个体为体重偏高个体，如果根据引物扩增出来的条带突变位点为“-”标记，则证明该个体为体重偏轻个体。据此，可通过对公牛和母牛进行这16个位点的分子检测，如果母牛和公牛都携带高体重的分析标记位点，则可优先配种，对其后代加以选育，经过几代的繁殖则可极大的提升牦牛的体重，增加产肉能力，给牧区人民带来实实在在的经济收入。

上述实施例仅为本发明的优选实施方式之一，不应当用于限制本发明的保护范围，但凡在本发明的主体设计思想和精神上作出的毫无实质意义的改动或润色，其所解决的技术问题仍然与本发明一致的，均应当包含在本发明的保护范围之内。

SEQUENCE LISTING

<110> 西藏自治区农牧科学院畜牧兽医研究所

西南民族大学

成都百因生物科技有限公司

<120> 牦牛全基因组SNP位点的应用及检测用引物组和试剂盒

<130> 20180205

<160> 48

<170> PatentIn version 3.3

<210> 1

<211> 71

<212> DNA

<213> 人工序列

<220>

<221> misc_feature

<222> (36)..(36)

<223> n is a, or g

<400> 1

aggggtcctg ccccatgctc caccaggcag agcttnttct cccagcatct gagctctagg 60

cttttactca t 71

<210> 2

<211> 71

<212> DNA

<213> 人工序列

<220>

<221> misc_feature

<222> (36)..(36)

<223> n is t, or g

<400> 2

gggcatgata ggctaggttt ttaccctcag ggaacntaac atctaatggg agagacaatc 60

aagtaattct a 71

<210> 3

<211> 71

<212> DNA

<213> 人工序列

<220>

<221> misc_feature

<222> (36)..(36)

<223> n is a, or g

<400> 3

ctgctggtcc cgtcaaaact ttctgcctta cagtcngccc tgctgacttg tggccgcaaa 60

gagtcaagat t 71

<210> 4

<211> 71

<212> DNA

<213> 人工序列

<220>

<221> misc_feature

<222> (36)..(36)

<223> n is a, or g

<400> 4

cagtcggccc tgctgacttg tggccgcaaa gagtcnagat tcactcacgc gacttcacag 60

agtagacaag g 71

<210> 5

<211> 71

<212> DNA

<213> 人工序列

<220>

<221> misc_feature

<222> (36)..(36)

<223> n is a, or g

<400> 5

gctggactaa attcagagag agaggggtgg tttcantttc cagccgggac tattcattcg 60

gcaccatgca g 71

<210> 6

<211> 71

<212> DNA

<213> 人工序列

<220>

<221> misc_feature

<222> (36)..(36)

<223> n is a, or g

<400> 6

cctgggcaag acataaaaga aaatcaagaa caaacngtta gcagtatcta aaattcttgc 60

agatcatgag g 71

<210> 7

<211> 71

<212> DNA

<213> 人工序列

<220>

<221> misc_feature

<222> (36)..(36)

<223> n is a, or g

<400> 7

tggattatct gacattgagt ctaagccaat taagancttc ctcttgtttc tcctaaagga 60

aaatacccaa g 71

<210> 8

<211> 71

<212> DNA

<213> 人工序列

<220>

<221> misc_feature

<222> (36)..(36)

<223> n is a, or g

<400> 8

caatgccagt accagagagc atccttcatg ccttcnttca gtaccactat tgagcaacag 60

ctgtgtgcca g 71

<210> 9

<211> 71

<212> DNA

<213> 人工序列

<220>

<221> misc_feature

<222> (36)..(36)

<223> n is a, or g

<400> 9

gagaaaagtc caaactatca gttaagaata acatcnatcc actaggattg tctgcacttt 60

ttaaaaaaaa t 71

<210> 10

<211> 71

<212> DNA

<213> 人工序列

<220>

<221> misc_feature

<222> (36)..(36)

<223> n is a, or g

<400> 10

ctccaaggct attacaacac ctttaatagt gcctcntgat tatctgatta catcacatct 60

ccaagaggtc a 71

<210> 11

<211> 71

<212> DNA

<213> 人工序列

<220>

<221> misc_feature

<222> (36)..(36)

<223> n is t, or c

<400> 11

tcaaacctac tttagccaat agagtctgcc tttcangctg aactggccta ccctctgctc 60

gcatgaatgt g 71

<210> 12

<211> 71

<212> DNA

<213> 人工序列

<220>

<221> misc_feature

<222> (36)..(36)

<223> n is t, or g

<400> 12

cggaagaagt gccccgtaaa agtgactcag ctctcntgag ggcctgactg tctctttgaa 60

cctacccgtg t 71

<210> 13

<211> 71

<212> DNA

<213> 人工序列

<220>

<221> misc_feature

<222> (36)..(36)

<223> n is a, or g

<400> 13

gataactagt caggtttctt tctttgtagc actttntaag ttgagcaggg tcttagcata 60

cttttgtttt g 71

<210> 14

<211> 71

<212> DNA

<213> 人工序列

<220>

<221> misc_feature

<222> (36)..(36)

<223> n is t, or c

<400> 14

ggcccaccac ggcctttgtg ctttgcttag gacttnagga gagagggtct ctccccaaag 60

tccagcacag g 71

<210> 15

<211> 71

<212> DNA

<213> 人工序列

<220>

<221> misc_feature

<222> (36)..(36)

<223> n is a, or g

<400> 15

agaaagcact ttacatctct ttatggactg agaaangaaa ataccggctg gttgggaaag 60

gatagacaga c 71

<210> 16

<211> 71

<212> DNA

<213> 人工序列

<220>

<221> misc_feature

<222> (36)..(36)

<223> n is t, or c

<400> 16

cttgcaattg attccaaagg tccttccacc tatacnctcc tcccgcccag ggcaattatt 60

cacatggcta a 71

<210> 17

<211> 19

<212> DNA

<213> 人工序列

<400> 17

gggtcctgcc ccatgctcc 19

<210> 18

<211> 21

<212> DNA

<213> 人工序列

<400> 18

gcctagagct cagatgctgg g 21

<210> 19

<211> 25

<212> DNA

<213> 人工序列

<400> 19

gggcatgata ggctaggttt tt 22

<210> 20

<211> 25

<212> DNA

<213> 人工序列

<400> 20

tagaattact tgattgtctc tccca 25

<210> 21

<211> 21

<212> DNA

<213> 人工序列

<400> 21

ctgctggtcc cgtcaaaact t 21

<210> 22

<211> 21

<212> DNA

<213> 人工序列

<400> 22

tgactctttg cggccacaag t 21

<210> 23

<211> 21

<212> DNA

<213> 人工序列

<400> 23

gacttgtggc cgcaaagagt c 21

<210> 24

<211> 21

<212> DNA

<213> 人工序列

<400> 24

ctgtgaagtc gcgtgagtga a 21

<210> 25

<211> 24

<212> DNA

<213> 人工序列

<400> 25

gctggactaa attcagagag agag 24

<210> 26

<211> 21

<212> DNA

<213> 人工序列

<400> 26

tggtgccgaa tgaatagtcc c 21

<210> 27

<211> 21

<212> DNA

<213> 人工序列

<400> 27

ggcaagacat aaaagaaaat c 21

<210> 28

<211> 21

<212> DNA

<213> 人工序列

<400> 28

gatactgcta accgtttgtt c 21

<210> 29

<211> 21

<212> DNA

<213> 人工序列

<400> 29

tggattatct gacattgagt c 21

<210> 30

<211> 21

<212> DNA

<213> 人工序列

<400> 30

ttaggagaaa caagaggaag t 21

<210> 31

<211> 21

<212> DNA

<213> 人工序列

<400> 31

tgccagtacc agagagcatc c 21

<210> 32

<211> 22

<212> DNA

<213> 人工序列

<400> 32

gcacacagct gttgctcaat ag 22

<210> 33

<211> 20

<212> DNA

<213> 人工序列

<400> 33

gagaaaagtc caaactatca 20

<210> 34

<211> 20

<212> DNA

<213> 人工序列

<400> 34

aaaagtgcag acaatcctag 20

<210> 35

<211> 20

<212> DNA

<213> 人工序列

<400> 35

caacaccttt aatagtgcct 20

<210> 36

<211> 20

<212> DNA

<213> 人工序列

<400> 36

acctcttgga gatgtgatgt 20

<210> 37

<211> 20

<212> DNA

<213> 人工序列

<400> 37

tcaaacctac tttagccaat 20

<210> 38

<211> 20

<212> DNA

<213> 人工序列

<400> 38

ttcatgcgag cagagggtag 20

<210> 39

<211> 20

<212> DNA

<213> 人工序列

<400> 39

gaagaagtgc cccgtaaaag 20

<210> 40

<211> 20

<212> DNA

<213> 人工序列

<400> 40

acacgggtag gttcaaagag 20

<210> 41

<211> 20

<212> DNA

<213> 人工序列

<400> 41

tagtcaggtt tctttctttg 20

<210> 42

<211> 20

<212> DNA

<213> 人工序列

<400> 42

caaaacaaaa gtatgctaag 20

<210> 43

<211> 20

<212> DNA

<213> 人工序列

<400> 43

acggcctttg tgctttgctt 20

<210> 44

<211> 20

<212> DNA

<213> 人工序列

<400> 44

tgtgctggac tttggggaga 20

<210> 45

<211> 25

<212> DNA

<213> 人工序列

<400> 45

agcactttac atctctttat ggact 25

<210> 46

<211> 23

<212> DNA

<213> 人工序列

<400> 46

gtctatcctt tcccaaccag ccg 23

<210> 47

<211> 20

<212> DNA

<213> 人工序列

<400> 47

attgattcca aaggtccttc 20

<210> 48

<211> 20

<212> DNA

<213> 人工序列

<400> 48

tagccatgtg aataattgcc 20