CN118086517A - 用于检测摩拉水牛繁殖性状的snp分子标记组合及应用 - Google Patents
用于检测摩拉水牛繁殖性状的snp分子标记组合及应用 Download PDFInfo
- Publication number
- CN118086517A CN118086517A CN202410105921.8A CN202410105921A CN118086517A CN 118086517 A CN118086517 A CN 118086517A CN 202410105921 A CN202410105921 A CN 202410105921A CN 118086517 A CN118086517 A CN 118086517A
- Authority
- CN
- China
- Prior art keywords
- buffalo
- murrah
- type
- breeding
- molecular marker
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000003147 molecular marker Substances 0.000 title claims abstract description 39
- 241000357437 Mola Species 0.000 title abstract 2
- 238000009395 breeding Methods 0.000 claims abstract description 69
- 230000001488 breeding effect Effects 0.000 claims abstract description 67
- 230000001850 reproductive effect Effects 0.000 claims description 37
- 108090000623 proteins and genes Proteins 0.000 claims description 34
- 108020004414 DNA Proteins 0.000 claims description 28
- 210000000349 chromosome Anatomy 0.000 claims description 19
- 230000013011 mating Effects 0.000 claims description 19
- 238000000034 method Methods 0.000 claims description 18
- 238000012163 sequencing technique Methods 0.000 claims description 17
- 108091028043 Nucleic acid sequence Proteins 0.000 claims description 13
- 239000002773 nucleotide Substances 0.000 claims description 11
- 125000003729 nucleotide group Chemical group 0.000 claims description 11
- 238000005516 engineering process Methods 0.000 claims description 6
- 230000003321 amplification Effects 0.000 claims description 4
- 238000003205 genotyping method Methods 0.000 claims description 4
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 4
- 238000012098 association analyses Methods 0.000 description 12
- 108700028369 Alleles Proteins 0.000 description 11
- 230000000694 effects Effects 0.000 description 8
- 239000008280 blood Substances 0.000 description 7
- 210000004369 blood Anatomy 0.000 description 7
- 230000002068 genetic effect Effects 0.000 description 6
- 239000003550 marker Substances 0.000 description 5
- 238000007400 DNA extraction Methods 0.000 description 4
- 238000001914 filtration Methods 0.000 description 4
- 230000006872 improvement Effects 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 241000283690 Bos taurus Species 0.000 description 3
- 238000012408 PCR amplification Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 210000004731 jugular vein Anatomy 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 235000013372 meat Nutrition 0.000 description 3
- 239000008267 milk Substances 0.000 description 3
- 210000004080 milk Anatomy 0.000 description 3
- 235000013336 milk Nutrition 0.000 description 3
- 238000003908 quality control method Methods 0.000 description 3
- 238000007619 statistical method Methods 0.000 description 3
- 238000013517 stratification Methods 0.000 description 3
- 238000000246 agarose gel electrophoresis Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 108020004707 nucleic acids Proteins 0.000 description 2
- 102000039446 nucleic acids Human genes 0.000 description 2
- 150000007523 nucleic acids Chemical class 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000000717 retained effect Effects 0.000 description 2
- 238000012216 screening Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000011144 upstream manufacturing Methods 0.000 description 2
- 238000010200 validation analysis Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 238000001712 DNA sequencing Methods 0.000 description 1
- 241001465754 Metazoa Species 0.000 description 1
- 238000000137 annealing Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 244000309466 calf Species 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000004925 denaturation Methods 0.000 description 1
- 230000036425 denaturation Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 235000013325 dietary fiber Nutrition 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 239000012154 double-distilled water Substances 0.000 description 1
- 239000000499 gel Substances 0.000 description 1
- 230000007614 genetic variation Effects 0.000 description 1
- 229920000669 heparin Polymers 0.000 description 1
- ZFGMDIBRIDKWMY-PASTXAENSA-N heparin Chemical compound CC(O)=N[C@@H]1[C@@H](O)[C@H](O)[C@@H](COS(O)(=O)=O)O[C@@H]1O[C@@H]1[C@@H](C(O)=O)O[C@@H](O[C@H]2[C@@H]([C@@H](OS(O)(=O)=O)[C@@H](O[C@@H]3[C@@H](OC(O)[C@H](OS(O)(=O)=O)[C@H]3O)C(O)=O)O[C@@H]2O)CS(O)(=O)=O)[C@H](O)[C@H]1O ZFGMDIBRIDKWMY-PASTXAENSA-N 0.000 description 1
- 238000009396 hybridization Methods 0.000 description 1
- 244000144972 livestock Species 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 230000032696 parturition Effects 0.000 description 1
- 102000054765 polymorphisms of proteins Human genes 0.000 description 1
- 238000012257 pre-denaturation Methods 0.000 description 1
- 238000000746 purification Methods 0.000 description 1
- 238000005215 recombination Methods 0.000 description 1
- 230000006798 recombination Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
Classifications
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6876—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
- C12Q1/6888—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for detection or identification of organisms
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6844—Nucleic acid amplification reactions
- C12Q1/6858—Allele-specific amplification
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q2600/00—Oligonucleotides characterized by their use
- C12Q2600/124—Animal traits, i.e. production traits, including athletic performance or the like
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q2600/00—Oligonucleotides characterized by their use
- C12Q2600/156—Polymorphic or mutational markers
Landscapes
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Organic Chemistry (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Zoology (AREA)
- Wood Science & Technology (AREA)
- Analytical Chemistry (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Biophysics (AREA)
- Immunology (AREA)
- Microbiology (AREA)
- Molecular Biology (AREA)
- Biotechnology (AREA)
- Physics & Mathematics (AREA)
- Biochemistry (AREA)
- Bioinformatics & Cheminformatics (AREA)
- General Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Genetics & Genomics (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明公开了一种用于检测摩拉水牛繁殖性状的SNP分子标记组合及其应用。SNP分子标记组合由5个SNP分子标记组成，包括chr7_783387、chr19_1328487、chr4_86313、chr11_503605和chr5_174492，该SNP分子标记可用于摩拉水牛不同的繁殖性状鉴定，具有准确性高、稳定性好、成本低、操作简单、方便快速等优点，可用于摩拉水牛分子标记辅助育种，加快选育良种摩拉水牛，可显著提高选择效率，缩短育种年限，降低饲养和育种成本，提高育种效率，具有广泛的应用前景。
Description
技术领域
本发明属于分子生物学技术领域，涉及用于检测摩拉水牛繁殖性状的SNP分子标记组合及应用。
背景技术
繁殖性状是水牛经济性状的重要组成部分，提高水牛繁殖率是增加养殖业经济效益的关键环节。水牛是单胎动物，通常一胎产一犊，繁殖效率低、世代间隔长，导致种牛的培育和群体遗传改良需要很长的时间跨度。因此，提高水牛的繁殖力是新品种培育的主要目标之一。而水牛的繁殖性状表现出低到中等的遗传力，利用传统育种方法通常很难在短时间内在一个品种内进行遗传改良。因此，如何利用现代分子育种技术培育出具有优良经济性状的水牛品种具有深远意义。
摩拉水牛自印度引入我国经过60多年的繁育，该品种在繁殖性能、生长发育、产肉和产奶性能方面均与原产地相当。在我国南方亚热带气候环境下，具有耐热、耐粗饲、抗病力强、生长正常、适应性强等特点，与本地沼泽型水牛相比具有较好的乳肉兼用性能。用摩拉水牛与本地水牛进行杂交改良，可大幅度提高杂交后代的产肉产奶性能。
目前，纯种摩拉水牛饲养数量少，规模化程度较低，种质资源匮乏，大部分牛场也只是通过传统的育种方法进行选种选育，造成水牛遗传改良进展缓慢。近些年来，随着科技的飞速发展，分子标记辅助育种已成为改良遗传性状的新方法。全基因组关联分析(GWAS)是一种对全基因组范围内的常见遗传变异(单核苷酸多态性和拷贝数)总体关联分析的方法，该方法以自然群体为研究对象，以长期重组后保留下来的基因(位点)间连锁不平衡(LD)为基础，将目标性状表型的多样性与基因(或标记位点)的多态性结合起来分析，可直接鉴定出与表型变异密切相关且具有特定功能的基因位点或标记位点。采用GWAS技术在全基因组范围内进行研究，能够一次性对多个性状进行定位，适用于定位性状关联区间、功能基因研究、开发性状选育和功能标记等方面的研究，其得到的结果更具可靠性。该技术作为一种新的方法必将在家畜育种领域得到广泛应用。
发明内容
本发明的发明目的是，针对上述问题，提供了一种用于检测摩拉水牛繁殖性状的SNP分子标记组合及应用。
为达到上述目的，本发明所采用的技术方案是：
一种用于检测摩拉水牛繁殖性状的SNP分子标记组合，所述SNP分子标记组合包括chr7_783387、chr19_1328487、chr4_86313、chr11_503605和chr5_174492，其中：
chr7_783387位于如SEQ ID NO.1所示的核苷酸序列的第51位，多态性为A或G；
chr19_1328487位于如SEQ ID NO.2所示的核苷酸序列的第51位，多态性为T或G；
chr4_86313位于如SEQ ID NO.3所示的核苷酸序列的第51位，多态性为A或G；
chr11_503605位于如SEQ ID NO.4所示的核苷酸序列的第51位，多态性为C或T；
chr5_174492位于如SEQ ID NO.5所示的核苷酸序列的第51位，多态性为A或G；
所述繁殖性状为初配月龄和/或产犊间隔繁殖性状。
本发明的另一个目的在于提供用于扩增如如上所述的SNP分子标记组合的引物组。
进一步说明，所述的引物组，所述引物组包含Forward primer和Reverse primer；其中，
所述chr7_783387分子标记的引物组中Forward primer和Reverse primer引物的序列分别如核苷酸序列SEQ ID NO.6、SEQ ID NO.7所示；
所述chr19_1328487分子标记的引物组中Forward primer和Reverse primer引物的序列分别如核苷酸序列SEQ ID NO.8、SEQ ID NO.9所示；
所述chr4_86313分子标记的引物组中Forward primer和Reverse primer引物的序列分别如核苷酸序列SEQ ID NO.10、SEQ ID NO.11所示；
所述chr11_503605分子标记的引物组中Forward primer和Reverse primer引物的序列分别如核苷酸序列SEQ ID NO.12、SEQ ID NO.13所示；
所述chr5_174492分子标记的引物组中Forward primer和Reverse primer引物的序列分别如核苷酸序列SEQ ID NO.14、SEQ ID NO.15所示。
本发明的另一个目的在于提供一种试剂盒，包括如上述所述引物组。
本发明的另一个目的在于提供所述的SNP分子标记组合或上述所述的引物组在摩拉水牛分子标记辅助育种中的应用。
本发明的另一个目的在于提供所述的SNP分子标记组合或上述所述的引物组在摩拉水牛品种繁殖性状鉴定中的应用。
本发明的另一个目的在于提供利用分子生物学技术检测摩拉水牛基因型的方法，包括以下步骤：根据权利要求1所述的5个SNP分子标记位点两侧翼的核苷酸序列分别设计扩增引物，以摩拉水牛个体DNA为模板扩增，将扩增产物进行一代测序，根据测序结果对摩拉水牛个体进行基因分型，从而鉴定待测摩拉水牛个体的基因型。
进一步说明，所述的方法在摩拉水牛分子标记辅助育种中的应用。
本发明还提供一种应用上述的SNP分子标记选育/辅助选育与摩拉水牛繁殖性状相关的摩拉水牛品种或品系的方法，所述方法是：
提取摩拉水牛的基因组DNA，检测7号染色体的第783387位点的脱氧核苷酸，测出第783387位点的脱氧核苷酸为A或G，确定待测摩拉水牛的基因型是AA型、GA或GG型，据育种需求，选择AA型、GA或GG型基因的摩拉水牛进行下一步的选种和/或育种，，所述GG型基因的摩拉水牛初配月龄远低于GA和AA型；
提取摩拉水牛的基因组DNA，检测19号染色体的第1328487位点的脱氧核苷酸，测出第1328487位点的脱氧核苷酸为T或G，确定待测摩拉水牛的基因型是TT型、GT或GG型，据育种需求，选择TT型、GT或GG型基因的摩拉水牛进行下一步的选种和/或育种，所述TT型基因的摩拉水牛产犊间隔远低于GT和GG型；
提取摩拉水牛的基因组DNA，检测4号染色体的第86313位点的脱氧核苷酸，测出第86313位点的脱氧核苷酸为A或G，确定待测摩拉水牛的基因型是AA型、GA或GG型，据育种需求，选择AA型、GA或GG型基因的摩拉水牛进行下一步的选种和/或育种，所述AA型基因的摩拉水牛产犊间隔远低于GA和GG型；
提取摩拉水牛的基因组DNA，检测11号染色体的第503605位点的脱氧核苷酸，测出第503605位点的脱氧核苷酸为T或C，确定待测摩拉水牛的基因型是TT型、CT或CC型，据育种需求，选择TT型、CT或CC型基因的摩拉水牛进行下一步的选种和/或育种，所述TT型基因的摩拉水牛产犊间隔远低于CT和CC型；
提取摩拉水牛的基因组DNA，检测5号染色体的第174492位点的脱氧核苷酸，测出第174492位点的脱氧核苷酸为A或G，确定待测摩拉水牛的基因型是AA型、GA或GG型，据育种需求，选择AA型、GA或GG型基因的摩拉水牛进行下一步的选种和/或育种，所述AA型基因的摩拉水牛产犊间隔远低于GA和GG型。
由于采用上述技术方案，本发明具有以下有益效果：
本发明提供摩拉水牛繁殖性状表型相关的SNP分子标记组合，可用于摩拉水牛的分子标记辅助育种，运用全基因组关联分析筛选与摩拉水牛繁殖性状显著相关的SNP分子标记，可用于摩拉水牛分子标记辅助育种，加快选育良种摩拉水牛。本发明利用SNP分子标记对摩拉水牛繁殖性状进行鉴定，可显著提高选择效率，缩短育种年限，降低饲养和育种成本，提高育种效率，具有广泛的应用前景。
附图说明
图1是本发明摩拉水牛繁殖性状全基因组关联分析的初配月龄曼哈顿图。附图标记说明：涉及摩拉水牛初配月龄性状，高于虚线上的为本发明筛选的SNP分子标记chr7_783387，该标记位于摩拉水牛7号染色体上。
图2是本发明摩拉水牛繁殖性状全基因组关联分析的初配月龄qq图。
图3是本发明摩拉水牛繁殖性状全基因组关联分析的产犊间隔曼哈顿图。附图标记说明：涉及摩拉水牛产犊间隔性状，高于虚线上的为本发明筛选的4个SNP分子标记chr19_1328487、chr4_86313、chr11_503605、chr5_174492，这些标记分别位于摩拉水牛19、4、11、5号染色体上。
图4是本发明摩拉水牛繁殖性状全基因组关联分析的产犊间隔qq图。
具体实施方式
以下结合附图对发明的具体实施进一步说明。
实施例1用于检测摩拉水牛繁殖性状的SNP分子标记组合的筛选与开发
1、群体选择与繁殖性状表型数据收集
在广西水牛研究所种牛场选取身体健康、发育良好的摩拉水牛79头，同时收集其初配月龄和1～2胎产犊间隔繁殖性状相关数据。
2、繁殖性状表型数据收集整理和统计分析
利用SPSS20软件对收集到的表型数据进行统计分析，包括最小值、最大值、平均值、标准差、变异系数，结果见表1。
表1摩拉水牛繁殖性状表型数据统计分析
| 表型 | 最小值 | 最大值 | 平均值 | 标准差 | 变异系数 |
| 初配月龄(月) | 16.33 | 55.6 | 30.62 | 8.28 | 27.04 |
| 产犊间隔(天) | 331 | 584 | 426.71 | 70.47 | 16.51 |
3、DNA提取和检测
收集每头牛颈静脉血5mL，肝素钠抗凝，轻轻摇晃使其混匀，采用血液基因组DNA提取试剂盒提取水牛基因组DNA，通过核酸浓度测定仪检测DNA浓度，1％的琼脂糖凝胶电泳测定DNA样品质量。将检测合格的DNA样品送往北京诺禾致源生物信息科技有限公司进行全基因组重测序，利用IlluminaHiSeq PE150测序平台，进行双末端(Paired-End)150测序，重测序的深度在5X以上，对测序下机数据进行过滤并获取有效信息，原始数据过滤方法：(1)需要过滤掉含有接头序列的reads；(2)当单端测序read中N的含量超过该条read长度比例的10％时，需要去除此对paired reads；(3)当单端测序read中含有的低质量(<＝5)碱基数超过该条read长度比例的50％时，需要去除此对paired reads；经过对测序数据的严格过滤，得到高质量的clean data。对产出数据进行统计，包括测序数据产量，测序错误率，Q20含量，Q30含量，GC含量如表2所示，最后得到高质量的clean data数据量为56.65Gb。
表2摩拉水牛基因组DNA测序结果统计
4、基因组数据比对
有效的高质量测序数据通过BWA生物信息软件(参数：mem-t 4-k 32-M)比对到参考基因组。参考基因组下载地址：ftp://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/000/ 471/725/GCF_000471725.1_UMD_CASPUR_WB_2.0/GCF_000471725.1_UMD_CASPUR_WB_2.0_ genomic.fna.gz群体样本平均比对率为99.4645％，对基因组的平均测序深度为13.9792X， 平均1X覆盖度(至少有一个碱基的覆盖)3.60％，信息见表3。
表3摩拉水牛基因组数据比对结果统计
5、SNP质控及过滤
采用SAMTOOLS软件进行群体SNP的检测。利用贝叶斯模型检测群体中的多态性位点。通过质控方法：(1)将质量值在20以上(错误率大于1％)的SNPs过滤掉；(2)若检测到两个SNP之间距离在5bp之内，将SNP去除；(3)SNP的覆盖深度需在平均深度的1/3-5倍之间。经过质控初步获得1176310个SNP位点，对获得的SNPs进行过滤，以获得高质量的SNPs，过滤条件为dp2、Miss0.2、Maf0.01，最后共获得122549个SNP位点用于后续的关联分析。
6、全基因组关联分析
采用GEMMA软件结合繁殖性状表型信息及基因组SNP信息进行关联分析，通过关联的显著性(P-value)，当-log10(P)>5时差异显著，筛选出潜在的候选SNPs。结果如图1～4所示，其中曼哈顿图，为遗传标记效应值即经F检验的全基因组P值按染色体上物理位置排序图，横坐标为基因组染色体上物理位置，纵坐标为-log10P，P值越小关联性越强，表现为纵坐标越大。Manhattan图中水平的虚线表示显著性水平，当-log10(P)>5时，认为该SNP与该性状显著关联。qq图(Quantile-quantileplot)，表示实际P值和无关联零假设期望P值的分布，用于检测群体分层和个体亲缘关系对关联分析的影响，如果观测到的P值与期望P值仅在分布最右端出现，表明性状并非由群体分层造成。QQ plot主要是用来估计数量性状观测值与预测值之间的差异。一般我们所取得的数量性状数据都为正态分布数据。在GWAS研究中QQ plot的X和Y轴主要是代表各个SNP的-lgPvalues。预测的线是一条从原点发出的45°角的虚线。GWAS分析过程中，个体亲缘关系和群体分层是造成假关联的主要因素。因此，采用混合线性模型进行性状关联分析，群体遗传结构作为固定效应，个体亲缘关系作为随机效应，以校正群体结构和个体亲缘关系的影响：
y＝Xα+Zβ+Wμ+e
其中y为表型性状，X为固定效应的指示矩阵，α为固定效应的估计参数；Z为SNP的指示矩阵，β为SNP的效应；W为随机效应的指示矩阵，μ为预测的随机个体，e是随机残差，服从e～(0，δe2)。
7、筛选并提取与摩拉水牛初配月龄和1～2胎产犊间隔繁殖性状相关的SNP分子标记
利用PLINK的“--indep-pairwise 2550.2”参数获得基因组上独立的SNP位点数，以1/SNP位点数作为多重检验阈值。通过R软件提取达到显著关联水平的多态位点，优选出5个分别与摩拉水牛的初配月龄和1～2胎产犊间隔相关联的SNP分子标记位点，可以用于摩拉水牛初配月龄和1～2胎产犊间隔繁殖性状的选育。详细信息见表4。
表4筛选得到的5个分子标记的信息
实施例2：5个SNP标记位点引物组在摩拉水牛繁殖性状检测中的应用
挑选摩拉水牛52头。颈静脉采集血液并采用血液基因组DNA提取试剂盒提取基因组DNA，采用核酸浓度测定仪检测DNA浓度，1％的琼脂糖凝胶电泳测定DNA样品质量。针对筛选出来的5个候选SNP位点设计扩增引物。引物序列为：
表55个候选SNP的引物设计
注：表中的F即上游引物Forwardprimer，R即下游引物Reverse primer
利用上述引物，以摩拉水牛血液基因组DNA为模板进行PCR扩增。采用50μLPCR反应体系：ddH2O 19μL，Premix TaqTM25μL，DNA模板2.0μL，引物(上游引物和下游引物均为10μmol/L)各2.0μL。
PCR反应条件：
95℃预变性4min；94℃变性10s，退火30s(温度为60℃)，72℃延伸1min，共35个循环；72℃延伸5min。
用上海生工生物工程有限公司的Gel Extraction Kit试剂盒对上述PCR扩增产物进行纯化，具体步骤见试剂盒说明书。将上述获得的PCR纯化产物回收后直接送至华大基因(深圳)生物科技公司进行一代测序。根据测序结果对个体基因分型。基因分型结果见表5。
表65个候选SNP位点在摩拉水牛基因型频率和等位基因频率
从表6中可以看出突变位点chr7_783387的G等位基因频率要明显大于A等位基因频率；chr19_1328487的G等位基因频率要明显大于T等位基因频率；chr4_86313的A等位基因频率要明显大于G等位基因频率；chr11_503605的T等位基因频率要明显大于C等位基因频率；chr5_174492的G等位基因频率要明显大于A等位基因频率。
表7摩拉水牛初配月龄显著关联SNP验证
初配月龄表型值以“最小二乘均值±标准差”表示，同列数据肩标不同字母表示差异显著(P<0.05)；肩标相同字母或无字母标注表示差异不显著(P>0.05)；SNP位点基因型按照突变型、杂合型、参考型依次排列。
结果表明，通过验证发现候选chr7_783387不同基因型与摩拉水牛初配月龄具有显著差异(表7)，GG型初配月龄显著小于GA和AA型，GA和AA型差异不显著，说明初配月龄小，则繁殖间隔短，能显著缩短繁殖年限，可作为鉴定摩拉水牛繁殖性状初配月龄的分子标记辅助育种，应用到生产中。
表8摩拉水牛产犊间隔显著关联SNP验证
产犊间隔表型值以“最小二乘均值±标准差”表示，同列数据肩标不同字母表示差异显著(P<0.05)；肩标相同字母或无字母标注表示差异不显著(P>0.05)；SNP位点基因型按照突变型、杂合型、参考型依次排列。
结果表明，通过验证发现候选chr19_1328487、chr4_86313、chr11_503605、chr5_174492不同基因型与摩拉水牛产犊间隔具有显著差异(见表8)，chr19_1328487的TT型产犊间隔显著短于GT和GG型，GT和GG型差异不显著；chr11_503605的TT型产犊间隔显著短于CT和CC型，CT和CC型差异不显著；chr4_86313和chr5_174492的AA型产犊间隔显著少于GA和GG型，GA和GG型差异不显著，说明产犊间隔短，则繁殖周期短，能显著缩短繁殖年限，可作为鉴定摩拉水牛繁殖性状产犊间隔的分子标记辅助育种，应用到生产中。
实施例4：
在摩拉水牛育种时，可以利用以上5个SNP分子标记，在其两侧翼的核苷酸序列上分别设计扩增引物，在摩拉水牛3～6月龄时采集水牛颈静脉血液，采用血液基因组DNA提取试剂盒提取基因组DNA，以摩拉水牛个体DNA为模板进行PCR扩增，将扩增产物一代测序，根据测序结果进行基因分型，保留带有目的性状基因型的个体，可以缩短育种年限，降低饲养和育种成本，提高育种效率，加快育种进程。
可应用本申请的SNP分子标记选育/辅助选育与摩拉水牛繁殖性状相关的摩拉水牛品种或品系的方法，所述方法是：
提取摩拉水牛的基因组DNA，检测7号染色体的第783387位点的脱氧核苷酸，测出第783387位点的脱氧核苷酸为A或G，确定待测摩拉水牛的基因型是AA型、GA或GG型，据育种需求，选择AA型、GA或GG型基因的摩拉水牛进行下一步的选种和/或育种，，所述GG)型基因的摩拉水牛初配月龄远低于GA和AA型；
提取摩拉水牛的基因组DNA，检测19号染色体的第1328487位点的脱氧核苷酸，测出第1328487位点的脱氧核苷酸为T或G，确定待测摩拉水牛的基因型是TT型、GT或GG型，据育种需求，选择TT型、GT或GG型基因的摩拉水牛进行下一步的选种和/或育种，所述TT型基因的摩拉水牛产犊间隔远低于GT和GG型；
提取摩拉水牛的基因组DNA，检测4号染色体的第86313位点的脱氧核苷酸，测出第86313位点的脱氧核苷酸为A或G，确定待测摩拉水牛的基因型是AA型、GA或GG型，据育种需求，选择AA型、GA或GG型基因的摩拉水牛进行下一步的选种和/或育种，所述AA型基因的摩拉水牛产犊间隔远低于GA和GG型；
提取摩拉水牛的基因组DNA，检测11号染色体的第503605位点的脱氧核苷酸，测出第503605位点的脱氧核苷酸为T或C，确定待测摩拉水牛的基因型是TT型、CT或CC型，据育种需求，选择TT型、CT或CC型基因的摩拉水牛进行下一步的选种和/或育种，所述TT型基因的摩拉水牛产犊间隔远低于CT和CC型；
提取摩拉水牛的基因组DNA，检测5号染色体的第174492位点的脱氧核苷酸，测出第174492位点的脱氧核苷酸为A或G，确定待测摩拉水牛的基因型是AA型、GA或GG型，据育种需求，选择AA型、GA或GG型基因的摩拉水牛进行下一步的选种和/或育种，所述AA型基因的摩拉水牛产犊间隔远低于GA和GG型。
上述说明是针对本发明较佳可行实施例的详细说明，但实施例并非用以限定本发明的专利申请范围，凡本发明所提示的技术精神下所完成的同等变化或修饰变更，均应属于本发明所涵盖专利范围。
Claims (9)
1.一种用于检测摩拉水牛繁殖性状的SNP分子标记组合，其特征在于：所述SNP分子标记组合包括chr7_783387、chr19_1328487、chr4_86313、chr11_503605和chr5_174492，其中：
chr7_783387位于如SEQ ID NO.1所示的核苷酸序列的第51位，多态性为A或G；
chr19_1328487位于如SEQ ID NO.2所示的核苷酸序列的第51位，多态性为T或G；
chr4_86313位于如SEQ ID NO.3所示的核苷酸序列的第51位，多态性为A或G；
chr11_503605位于如SEQ ID NO.4所示的核苷酸序列的第51位，多态性为C或T；
chr5_174492位于如SEQ ID NO.5所示的核苷酸序列的第51位，多态性为A或G；
所述繁殖性状为初配月龄和/或产犊间隔繁殖性状。
2.用于扩增如权利要求1所述的SNP分子标记组合的引物组。
3.根据权利要求2所述的引物组，其特征在于：所述引物组包含Forwardprimer和Reverse primer；其中，
所述chr7_783387分子标记的引物组中Forwardprimer和Reverse primer引物的序列分别如核苷酸序列SEQ ID NO.6、SEQ ID NO.7所示；
所述chr19_1328487分子标记的引物组中Forwardprimer和Reverse primer引物的序列分别如核苷酸序列SEQ ID NO.8、SEQ ID NO.9所示；
所述chr4_86313分子标记的引物组中Forward primer和Reverse primer引物的序列分别如核苷酸序列SEQ ID NO.10、SEQ ID NO.11所示；
所述chr11_503605分子标记的引物组中Forward primer和Reverse primer引物的序列分别如核苷酸序列SEQ ID NO.12、SEQ ID NO.13所示；
所述chr5_174492分子标记的引物组中Forwardprimer和Reverse primer引物的序列分别如核苷酸序列SEQ ID NO.14、SEQ ID NO.15所示。
4.一种试剂盒，其特征在于，包括如权利要求2或3所述引物组。
5.如权利要求1所述的SNP分子标记组合或权利要求2或3任一项所述的引物组在摩拉水牛分子标记辅助育种中的应用。
6.如权利要求1所述的SNP分子标记组合或权利要求2或3任一项所述的引物组在摩拉水牛繁殖性状鉴定中的应用。
7.利用分子生物学技术检测摩拉水牛基因型的方法，其特征在于，包括以下步骤：根据权利要求1所述的5个SNP分子标记位点两侧翼的核苷酸序列分别设计扩增引物，以摩拉水牛个体DNA为模板扩增，将扩增产物进行一代测序，根据测序结果对摩拉水牛个体进行基因分型，从而鉴定待测摩拉水牛个体的基因型。
8.权利要求7所述的方法在摩拉水牛分子标记辅助育种中的应用。
9.一种应用权利要求1所述的SNP分子标记选育/辅助选育与摩拉水牛繁殖性状相关的摩拉水牛品种或品系的方法，其特征在于，所述方法是：
提取摩拉水牛的基因组DNA，检测7号染色体的第783387位点的脱氧核苷酸，测出第783387位点的脱氧核苷酸为A或G，确定待测摩拉水牛的基因型是AA型、GA或GG型，据育种需求，选择AA型、GA或GG型基因的摩拉水牛进行下一步的选种和/或育种，，所述GG型基因的摩拉水牛初配月龄远低于GA和AA型；
提取摩拉水牛的基因组DNA，检测19号染色体的第1328487位点的脱氧核苷酸，测出第1328487位点的脱氧核苷酸为T或G，确定待测摩拉水牛的基因型是TT型、GT或GG型，据育种需求，选择TT型、GT或GG型基因的摩拉水牛进行下一步的选种和/或育种，所述TT型基因的摩拉水牛产犊间隔远低于GT和GG型；
提取摩拉水牛的基因组DNA，检测4号染色体的第86313位点的脱氧核苷酸，测出第86313位点的脱氧核苷酸为A或G，确定待测摩拉水牛的基因型是AA型、GA或GG型，据育种需求，选择AA型、GA或GG型基因的摩拉水牛进行下一步的选种和/或育种，所述AA型基因的摩拉水牛产犊间隔远低于GA和GG型；
提取摩拉水牛的基因组DNA，检测11号染色体的第503605位点的脱氧核苷酸，测出第503605位点的脱氧核苷酸为T或C，确定待测摩拉水牛的基因型是TT型、CT或CC型，据育种需求，选择TT型、CT或CC型基因的摩拉水牛进行下一步的选种和/或育种，所述TT型基因的摩拉水牛产犊间隔远低于CT和CC型；
提取摩拉水牛的基因组DNA，检测5号染色体的第174492位点的脱氧核苷酸，测出第174492位点的脱氧核苷酸为A或G，确定待测摩拉水牛的基因型是AA型、GA或GG型，据育种需求，选择AA型、GA或GG型基因的摩拉水牛进行下一步的选种和/或育种，所述AA型基因的摩拉水牛产犊间隔远低于GA和GG型。
Priority Applications (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202411130066.2A CN118910277A (zh) | 2024-01-25 | 2024-01-25 | 一种与摩拉水牛繁殖性状相关的分子标记及其应用 |
| CN202411130020.0A CN118957090A (zh) | 2024-01-25 | 2024-01-25 | 用于检测摩拉水牛产犊间隔繁殖性状相关的snp分子标记及其应用 |
| CN202410105921.8A CN118086517B (zh) | 2024-01-25 | 2024-01-25 | 用于检测摩拉水牛繁殖性状的snp分子标记组合及应用 |
| CN202411129970.1A CN118755851A (zh) | 2024-01-25 | 2024-01-25 | 与摩拉水牛产犊间隔繁殖性状相关的分子标记及其应用 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410105921.8A CN118086517B (zh) | 2024-01-25 | 2024-01-25 | 用于检测摩拉水牛繁殖性状的snp分子标记组合及应用 |
Related Child Applications (3)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202411130020.0A Division CN118957090A (zh) | 2024-01-25 | 2024-01-25 | 用于检测摩拉水牛产犊间隔繁殖性状相关的snp分子标记及其应用 |
| CN202411130066.2A Division CN118910277A (zh) | 2024-01-25 | 2024-01-25 | 一种与摩拉水牛繁殖性状相关的分子标记及其应用 |
| CN202411129970.1A Division CN118755851A (zh) | 2024-01-25 | 2024-01-25 | 与摩拉水牛产犊间隔繁殖性状相关的分子标记及其应用 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN118086517A true CN118086517A (zh) | 2024-05-28 |
| CN118086517B CN118086517B (zh) | 2024-10-15 |
Family
ID=91154111
Family Applications (4)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202411130020.0A Pending CN118957090A (zh) | 2024-01-25 | 2024-01-25 | 用于检测摩拉水牛产犊间隔繁殖性状相关的snp分子标记及其应用 |
| CN202410105921.8A Active CN118086517B (zh) | 2024-01-25 | 2024-01-25 | 用于检测摩拉水牛繁殖性状的snp分子标记组合及应用 |
| CN202411129970.1A Pending CN118755851A (zh) | 2024-01-25 | 2024-01-25 | 与摩拉水牛产犊间隔繁殖性状相关的分子标记及其应用 |
| CN202411130066.2A Pending CN118910277A (zh) | 2024-01-25 | 2024-01-25 | 一种与摩拉水牛繁殖性状相关的分子标记及其应用 |
Family Applications Before (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202411130020.0A Pending CN118957090A (zh) | 2024-01-25 | 2024-01-25 | 用于检测摩拉水牛产犊间隔繁殖性状相关的snp分子标记及其应用 |
Family Applications After (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202411129970.1A Pending CN118755851A (zh) | 2024-01-25 | 2024-01-25 | 与摩拉水牛产犊间隔繁殖性状相关的分子标记及其应用 |
| CN202411130066.2A Pending CN118910277A (zh) | 2024-01-25 | 2024-01-25 | 一种与摩拉水牛繁殖性状相关的分子标记及其应用 |
Country Status (1)
| Country | Link |
|---|---|
| CN (4) | CN118957090A (zh) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113308554A (zh) * | 2021-06-28 | 2021-08-27 | 广西壮族自治区水牛研究所 | 一种与牛生长性状相关的snp分子标记及其应用 |
| CN116397032A (zh) * | 2023-01-30 | 2023-07-07 | 广西壮族自治区水牛研究所 | 与水牛体重生长性状相关的snp分子标记及其应用 |
-
2024
- 2024-01-25 CN CN202411130020.0A patent/CN118957090A/zh active Pending
- 2024-01-25 CN CN202410105921.8A patent/CN118086517B/zh active Active
- 2024-01-25 CN CN202411129970.1A patent/CN118755851A/zh active Pending
- 2024-01-25 CN CN202411130066.2A patent/CN118910277A/zh active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113308554A (zh) * | 2021-06-28 | 2021-08-27 | 广西壮族自治区水牛研究所 | 一种与牛生长性状相关的snp分子标记及其应用 |
| CN116397032A (zh) * | 2023-01-30 | 2023-07-07 | 广西壮族自治区水牛研究所 | 与水牛体重生长性状相关的snp分子标记及其应用 |
Non-Patent Citations (2)
| Title |
|---|
| JUN LI: "Novel insights into the genetic basis of buffalo reproductive performance", BMCGENOMICS, vol. 19, no. 814, 12 November 2018 (2018-11-12), pages 1 - 11, XP021262471, DOI: 10.1186/s12864-018-5208-6 * |
| NCBI: "GENBANK Reference: OY997248.1", NCBI, 11 January 2024 (2024-01-11), pages 1 - 2 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN118957090A (zh) | 2024-11-15 |
| CN118910277A (zh) | 2024-11-08 |
| CN118086517B (zh) | 2024-10-15 |
| CN118755851A (zh) | 2024-10-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111705140B (zh) | 体重性状相关的SNPs分子标记及其在湖羊辅助育种中的应用 | |
| CN111996264B (zh) | 一种猪snp分子标记在猪繁殖性状筛选和猪育种中的应用 | |
| CN114941033A (zh) | 一种基于snp位点辅助培育地方优质白羽鸡高产蛋品系的方法 | |
| CN117701738A (zh) | 与水牛初产月龄性状相关的snp分子标记及其应用 | |
| CN118441073A (zh) | 与绒山羊胸宽相关的单倍型分子标记及其应用 | |
| CN108085400A (zh) | 白来航鸡红羽致因突变基因型的鉴定及专用标记 | |
| CN111979332A (zh) | 用于选择母鸡储精能力期间受精率的snp分子标记及其应用 | |
| CN118086517A (zh) | 用于检测摩拉水牛繁殖性状的snp分子标记组合及应用 | |
| CN104694651B (zh) | 一种与二花脸母猪产仔性状相关的snp标记、检测方法及应用 | |
| CN118531133A (zh) | 一种与尼里-拉菲水牛产犊间隔性状相关的snp分子标记及其应用 | |
| CN118109605B (zh) | 与尼里-拉菲水牛生长性状相关的snp分子标记组合及应用 | |
| CN117701727B (zh) | 基于全基因组测序筛选的与摩拉水牛体尺、初生重相关的snp分子标记组合及应用 | |
| CN111979335A (zh) | 用于选择母鸡储精能力期间受精率的snp分子标记及其应用 | |
| CN111154888B (zh) | 一种用于检测比格犬封闭群的微卫星标记及其组合 | |
| CN111979334A (zh) | 与母鸡储精能力期间受精率相关的snp分子标记 | |
| CN119162341A (zh) | 一种与黄羽肉鸡饲料转化率相关st3gal4基因snp分子标记及其应用 | |
| CN119220703A (zh) | 一种与黄羽肉鸡饲料转化率相关cdk5rap3基因snp分子标记及其应用 | |
| CN118345174A (zh) | 一种文昌鸡特异snp位点及应用 | |
| CN118581233A (zh) | 一种商品猪液相芯片及其制备方法和应用 | |
| CN119287029A (zh) | 一种猪出生重相关的snp分子标记及其应用 | |
| CN119162342A (zh) | 一种与黄羽肉鸡饲料转化率相关tmem165基因snp分子标记及其应用 | |
| CN119410788A (zh) | 一种与猪总乳头数性状相关的snp分子标记及应用 | |
| CN118932080A (zh) | 与绵羊生长性状相关的InDel分子标记及其应用 | |
| CN115851962A (zh) | 一种与鸡的腹脂重性状相关的分子标记及其应用 | |
| CN104726578B (zh) | 一种与二花脸母猪产仔性状相关的snp标记、检测方法及应用 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |