CN110218799A - 猪剩余采食量性状的分子遗传标记及应用 - Google Patents
猪剩余采食量性状的分子遗传标记及应用 Download PDFInfo
- Publication number
- CN110218799A CN110218799A CN201910489835.0A CN201910489835A CN110218799A CN 110218799 A CN110218799 A CN 110218799A CN 201910489835 A CN201910489835 A CN 201910489835A CN 110218799 A CN110218799 A CN 110218799A
- Authority
- CN
- China
- Prior art keywords
- pig
- feed intake
- genetic marker
- molecular genetic
- boar
- 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
- 235000021050 feed intake Nutrition 0.000 title claims abstract description 76
- 239000003550 marker Substances 0.000 title claims abstract description 73
- 230000002068 genetic effect Effects 0.000 title claims abstract description 42
- 238000009395 breeding Methods 0.000 claims abstract description 34
- 230000001488 breeding effect Effects 0.000 claims abstract description 33
- 238000009360 aquaculture Methods 0.000 claims abstract description 6
- 244000144974 aquaculture Species 0.000 claims abstract description 6
- 238000000034 method Methods 0.000 claims description 25
- 108090000623 proteins and genes Proteins 0.000 claims description 22
- 241000282887 Suidae Species 0.000 claims description 12
- 230000000869 mutational effect Effects 0.000 claims description 12
- 238000001514 detection method Methods 0.000 claims description 7
- 230000035772 mutation Effects 0.000 claims description 5
- 210000000349 chromosome Anatomy 0.000 claims description 4
- 238000011144 upstream manufacturing Methods 0.000 claims description 4
- 101000685271 Homo sapiens Protein sel-1 homolog 2 Proteins 0.000 claims description 3
- 238000002372 labelling Methods 0.000 abstract description 9
- 238000004519 manufacturing process Methods 0.000 abstract description 9
- 238000012098 association analyses Methods 0.000 abstract description 7
- 235000013305 food Nutrition 0.000 abstract description 6
- 230000036528 appetite Effects 0.000 abstract description 2
- 235000019789 appetite Nutrition 0.000 abstract description 2
- 239000012636 effector Substances 0.000 abstract description 2
- 241000282898 Sus scrofa Species 0.000 description 62
- 230000000694 effects Effects 0.000 description 25
- 239000011159 matrix material Substances 0.000 description 16
- 238000004458 analytical method Methods 0.000 description 12
- 239000002773 nucleotide Substances 0.000 description 11
- 125000003729 nucleotide group Chemical group 0.000 description 11
- 235000019786 weight gain Nutrition 0.000 description 10
- 230000004584 weight gain Effects 0.000 description 9
- 108700028369 Alleles Proteins 0.000 description 6
- 241001465754 Metazoa Species 0.000 description 5
- 238000012216 screening Methods 0.000 description 5
- 235000015097 nutrients Nutrition 0.000 description 4
- 238000003908 quality control method Methods 0.000 description 4
- 230000000875 corresponding effect Effects 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000007613 environmental effect Effects 0.000 description 3
- 238000003205 genotyping method Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 241000283690 Bos taurus Species 0.000 description 2
- 108010078791 Carrier Proteins Proteins 0.000 description 2
- 108090000790 Enzymes Proteins 0.000 description 2
- 102000004190 Enzymes Human genes 0.000 description 2
- 238000007476 Maximum Likelihood Methods 0.000 description 2
- 235000015278 beef Nutrition 0.000 description 2
- 239000008280 blood Substances 0.000 description 2
- 210000004369 blood Anatomy 0.000 description 2
- 230000037396 body weight Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000002596 correlated effect Effects 0.000 description 2
- 238000012217 deletion Methods 0.000 description 2
- 230000037430 deletion Effects 0.000 description 2
- 235000019621 digestibility Nutrition 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 230000014509 gene expression Effects 0.000 description 2
- 230000000968 intestinal effect Effects 0.000 description 2
- 244000144972 livestock Species 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 230000000813 microbial effect Effects 0.000 description 2
- 238000011056 performance test Methods 0.000 description 2
- 244000144977 poultry Species 0.000 description 2
- 235000013594 poultry meat Nutrition 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000010187 selection method Methods 0.000 description 2
- 238000013179 statistical model Methods 0.000 description 2
- 210000001519 tissue Anatomy 0.000 description 2
- 241000252095 Congridae Species 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- 241000287828 Gallus gallus Species 0.000 description 1
- 239000002253 acid Substances 0.000 description 1
- 238000009412 basement excavation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000012214 genetic breeding Methods 0.000 description 1
- 101150044508 key gene Proteins 0.000 description 1
- 230000004060 metabolic process Effects 0.000 description 1
- 239000002777 nucleoside Substances 0.000 description 1
- 125000003835 nucleoside group Chemical group 0.000 description 1
- 238000012372 quality testing Methods 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
- 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)
- Analytical Chemistry (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Organic Chemistry (AREA)
- Zoology (AREA)
- Wood Science & Technology (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Microbiology (AREA)
- Immunology (AREA)
- Molecular Biology (AREA)
- Biotechnology (AREA)
- Biophysics (AREA)
- Physics & Mathematics (AREA)
- Biochemistry (AREA)
- Bioinformatics & Cheminformatics (AREA)
- General Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Genetics & Genomics (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本公开提供猪剩余采食量性状的分子遗传标记及应用，通过鉴定一个影响公猪剩余采食量的分子标记ALGA0093815，该标记不同基因型公猪的剩余采食量有极显著差异，并通过公猪剩余采食量与全基因组分子遗传标记的关联分析，成功筛选影响猪剩余采食量的大效应分子遗传标记，并通过杜长大商品猪群体验证所筛选的遗传标记，可有效降低生产过程中饲料消耗量，选留低剩余采食量的纯合猪，可有效降低生产过程中饲料消耗量，降低养猪生产成本，提高企业经济效益和竞争力，加快高饲料利用效率品系的选育进度，通过选留TT纯合种猪进入核心群，降低剩余采食量，有效减少饲料消耗量和养殖成本。
Description
技术领域
本公开涉及猪遗传基因技术领域，特别涉及猪剩余采食量性状的分子遗传标记及应用。
背景技术
饲料效率作为一个重要的经济性状，一直受到国内外养猪企业和种猪改良公司的重点关注。20世纪中期，国内外沿用饲料增重比(Feed/Gain ratio,F/G)研究饲料效率，由于利用饲料增重比研究饲料效率具有2头具有相同的饲料增重比的猪，其采食量和增重量可能各不相同的局限性(杨运清，王忠新，尹炳北.比值性状的指数选择方法探讨.畜牧兽医学报，31(3)：198-202(2000))；1963年，Koch(Koch R M,Swiger L A,Chambers D,etal.Efficiency of feed use in beef cattle.J Anim Sci,22(2):486-494(1963))首先提出了剩余采食量(Average Day Feed Intake,DFI)的概念，反映的是动物本身由遗传背景决定的代谢差异。研究表明，剩余采食量作为评价饲料利用效率的表型性状具有中等遗传力，有很好的选择反应。
利用SNP标记辅助选育饲料利用效率相关性状对养猪生产管理和企业经济效益具有重要影响。①低RFI猪的饲料利用效率高(Vigors S,Sweeney T,Oshea C J,et al.Pigsthat are divergent in feed efficiency,differ in intestinal enzyme andnutrient transporter gene expression,nutrient digestibility and microbialactivity.Animal,10(11):1848-1855(2016))，能降低生产中的饲料使用量和生产成本，进而节约了饲料资源，亦能一定程度减少猪的排污量，从而减缓猪与人类竞争粮食资源和养猪业环保问题的压力。②开发有效的分子标记用于饲料效率相关性状选育工作，大大缩短培育周期，降低培育成本，提高了选种准确度，加快遗传进展，能避免引种-退化-再引种的现象。
因此，挖掘和利用新的与剩余采食量相关基因对于猪的遗传育种有着重大意义。基于覆盖全基因组的高密度SNP数据和大群体的性状表型记录，可通过全基因组关联分析技术(GWAS)(Hirschhorn,J.N.&Daly,M.J.Genome-wide association studies forcommon diseases and complex traits.Nat.Rev.Genet.6,95–108(2005))准确定位控制性状的候选基因。尽管该技术仍然存在一些缺陷(De,R.,Bush,W.S.&Moore,J.H.Bioinformatics challenges in genome-wide association studies(GWAS).Methods Mol.Biol.1168,63–81(2014))，其已被广泛应用于人类复杂疾病候选基因挖掘和畜禽重要经济性状关键基因的定位。经典的GWAS一般基于Plink(Purcell,S.etal.PLINK:A Tool Set for Whole-Genome Association and Population-Based LinkageAnalyses.Am.J.Hum.Genet.813,559–575(2007))等软件对所有标记逐个进行单标记回归分析，继而设定一个显著阈值来筛选显著位点。这类方法往往面临计算强度大、过高估计标记效应、显著性阈值设定不合理等问题。为了进一步提高GWAS的效率，新方法和软件不断被提出。其中，一步法全基因组关联分析(wssGWAS(WANG,H.,MISZTAL,I.,AGUILAR,I.,LEGARRA,A.&MUIR,W.M.Genome-wide association mapping including phenotypes fromrelatives without genotypes.Genet Res 94,73–83(2012))；(Wang,H.et al.Genome-wide association mapping including phenotypes from relatives withoutgenotypes in a single-step(ssGWAS)for 6-week body weight in broilerchickens.Front.Genet.5,1–10(2014)))同时利用系谱、历史个体表型记录和基因型数据进行关联分析，适用于大量个体拥有表型记录而只有少量个体拥有基因型数据的情况，尤其适用于畜禽重要经济性状的全基因组关联分析。基于GBLUPf90软件(Misztal,I.etal.BLUPF90and related programs(BGF90).in Proc.7th WorldCongr.Genet.Appl.Livest.Prod.21–22(2002).doi:9782738010520)，可轻易实现wssGWAS。利用wssGWAS筛选出与猪剩余采食量相关的SNP分子标记，为家猪饲料效率性状的遗传选量提供了可行途径，对养猪业具有重要意义。
发明内容
针对上述技术问题，本公开提供猪剩余采食量性状的分子遗传标记及应用，通过检测该分子标记，可以应用于种猪的选育，加快高饲料利用效率品系的选育进度，通过鉴定一个影响公猪剩余采食量的分子标记ALGA0093815，该标记不同基因型公猪的剩余采食量有极显著差异，并通过公猪剩余采食量与全基因组分子遗传标记的关联分析，成功筛选影响猪剩余采食量的大效应分子遗传标记，应用于种猪的选育，选留低剩余采食量的纯合猪，并通过杜长大商品猪群体验证所筛选的遗传标记，可有效降低生产过程中饲料消耗量，提高企业经济效益和竞争力，其中涉及的ALGA0093815遗传标记即SNP号为ALGA0093815的突变位点，见NCBI中猪基因组数据库(Sscrofa11.1)。
本公开的SNP分子遗传标记(ALGA0093815)，参阅Ensembl数据库(http://asia.ensembl.org/Sus_scrofa/Search/New？db＝core)，得到登录号为ALGA0093815的基因片段(RS号为rs80818852),SNP分子遗传标记位于猪17号染色体22568436bp位置，且属于SEL1L2基因的内含子上，该位置为一个T>C突变(突变位点)，T>C即T为大频率的等位基因，C为小频率的等位基因，符号>为等位基因频率大小，该SNP分子遗传标记在突变位点上下游100bp序列如下：
5’-TTCAGTCAGGTCAGGAACAAGGCAGGAGCCAGGGGCAAATTCTTCTAGCTCATGAGTGACTTTAAAGATTAGAAGAGAGCTGCTCTCTAATCCACCTTCCR(T/C)ATCATCCACGTCCATGTGCATGTTTTGAAATAATATGAGTTAAATTTTATATAGGGGGATTTCCCATCATTGCTCAGTGGTAAAGAACCTGACTAGTGTC-3’；R为突变位点，上述序列101位核苷酸处的R是C或T时，即R(T/C)时，导致上述序列多态性；当上述核苷酸序列的第101位核苷酸为T时，猪显著降低剩余采食量，5’-和-3’分别表示核苷酸序列的5’端和3’端。
上述ALGA0093815标记TT基因型比CC基因型纯种杜洛克公猪个体间剩余采食量显著降低了63.52g，所以，T是有利于显著降低剩余采食量的等位基因。
在杜长大商品猪群体中，TT基因型猪的剩余采食量比CC基因型显著降低53.19g；所以，该位点多态性有利于降低杜长大商品猪剩余采食量，T是有利于降低剩余采食量的等位基因，通过检测ALGA0093815标记基因型辅助种猪选育，可通过选留TT基因型纯合种猪进入核心群，降低剩余采食量，有效减少饲料消耗量和养殖成本(由于猪的DNA是反向螺旋双链结构，两条链的突变位点核苷酸均为T时为TT基因型纯合猪，其中，每一条链有一个核苷酸序列，T表示一个突变位点为T，TT基因型是双链的突变位点都是T的纯合猪，同理，CC基因型为双链的突变位点都是C的纯合猪；CT基因型为一条链的突变位点是T另一条的突变位点是C的猪)。
筛选猪剩余采食量性状的分子遗传标记及应用的方法具体包括以下步骤：
1、分子标记的获得的流程步骤
1.1、采集公猪的耳组织样品或者血样，提取总DNA，并对DNA进行质量检测。采用GGP 50k SNP(GeneSeek，US)芯片进行基因分型，获得覆盖全基因组的SNP标记基因型。
1.2、根据最新版的猪参考基因组(Sscrofa11.1)，采用NCBI基因组比对程序(https://www.ncbi.nlm.nih.gov/)对所有SNP标记的物理位置进行更新。基因组位置未知的SNP不用于关联分析。
1.3、对于所有常染色体上的SNP标记，利用Plink软件进行质量控制，标准为：个体检出率≥90％；SNP检出率≥90％；小等位基因频率≥0.01；哈迪-温伯格平衡p值≥10-6。对于缺失基因型，采用Beagle软件(version 4.1)进行填充。
2、分子标记的验证的流程步骤
2.1、整理种猪系谱，主要包括公猪个体号、父亲、母亲和初生日期等信息。采用公式对全自动种猪生产性能测定***(FIRE，美国)记录的生长数据进行分析获得饲料增重比表型数据，用于表型-基因型关联分析。其中，FCR为饲料增重比；Wa为活体增重量；Wf为饲料消耗量。
2.2、统计模型，采用加权的一步法全基因组关联分析法(weighted single stepgenome-wide association study，wssGWAS)进行全基因组关联分析。该方法首先基于混合模型方程组来估计个体育种值，继而基于育种值模型与标记效应模型的等价关系将育种值转换为标记效应。本发明采用的全基因组关联分析模型如下：
y＝Xb+Za+Wp+e，
其中，y为饲料增重比观测值向量；X，Z和W为设计矩阵；b为固定效应向量(环境、日龄)；为育种值向量；为个体永久环境效应；为残差。H为同时整合系谱和SNP标记的亲缘关系矩阵，其逆矩阵计算公式如下：
其中，A为基于系谱的亲缘关系矩阵；A22为A中有基因型个体对应的分块矩阵；Gω＝0.9G+0.1A22，为基于全基因组SNP标记的亲缘关系矩，Z为小等位基因频率(minor allele frequency，MAF)校正后的基因型矩阵，其中0-2p，1-2p和2-2p分别代表AA，Aa和aa三种基因型，p为小等位基因频率；D为对角线矩阵，表示SNP的权重；pi为第i个标记的小等位基因频率；m为标记数量。
对于上述混合模型，采用AI-REML(average information restricted maximumlikelihood)法估计方差组分，并通过求解混合模型方程组获得育种值。通过迭代的方式获得标记权重，主要步骤如下：
第1步：初始化(t＝1)，D(t)＝I，G(t)＝λZD(t)Z′，
第2步：通过ssGBLUP计算个体育种值；
第3步：通过公式将个体育种值转换为SNP效应，其中为有基因型个体的育种值；
第4步：利用公式计算SNP权重用于下一轮迭代；
第5步：利用公式对SNP权重进行标准化，以保证方差一致；
第6步：利用公式G(t+1)＝λZD(t+1)Z′计算亲缘关系矩阵用于下一轮迭代；
第7步：令t＝t+1，并从第2步开始下一轮迭代。
上述步骤迭代三次，最终获得SNP标记效应。将第三轮迭代输出的标记效应作为最终的结果。计算过程主要通过在R统计分析平台编程调用BLUPF90软件来实现，其中AIREMLF90程序用于方差组分估计，BLUPF90程序用于计算育种值，postGSf90用于计算标记效应。
3、标记筛选
对于所有标记的效应值，取其绝对值画曼哈顿图，展示和筛选大效应的SNP标记。并采用方差分析和多重比较(R统计分析平台)，分析ALGA0093815标记不同基因型群体公猪剩余采食量差异情况。
本公开的有益效果为：本公开提供猪剩余采食量性状的分子遗传标记及应用，通过检测该分子标记，可以应用于种猪的选育，选留低剩余采食量的纯合猪，可有效降低生产过程中饲料消耗量，降低养猪生产成本，提高企业经济效益和竞争力，加快高饲料利用效率品系的选育进度，影响公猪剩余采食量的分子标记ALGA0093815，该标记不同基因型公猪的剩余采食量有极显著差异，通过检测ALGA0093815标记基因型辅助种猪选育，可通过选留TT纯合种猪进入核心群，降低剩余采食量，有效减少饲料消耗量和养殖成本。
附图说明
通过对结合附图所示出的实施方式进行详细说明，本公开的上述以及其他特征将更加明显，本公开附图中相同的参考标号表示相同或相似的元素，显而易见地，下面描述中的附图仅仅是本公开的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图，在附图中：
图1所示为本公开的筛选猪剩余采食量性状的分子遗传标记及应用的方法工作流程图；
图2所示为本公开的ALGA0093815标记基因组位置及剩余采食量全基因组SNP效应分布。
具体实施方式
以下将结合实施例和附图对本公开的构思、具体结构及产生的技术效果进行清楚、完整的描述，以充分地理解本公开的目的、方案和效果。需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。
如图1所示为根据本公开的筛选猪剩余采食量性状的分子遗传标记及应用的方法工作流程图，下面结合图1来阐述根据本公开的筛选猪剩余采食量性状的分子遗传标记及应用的方法。
本公开筛选猪剩余采食量性状的分子遗传标记及应用的方法，具体包括以下步骤：
(1)表型-系谱数据采集
本公开的基础研究群体为杜洛克公猪，全部来自广西某核心种猪场。完整系谱中包含4个世代735头种猪，其中2015-2018年间记录了370头杜洛克公猪的饲料增重比性状表型数据。验证群体为杜长大商品猪，记录了1157头商品猪(577头母猪，582头公猪)的饲料增重比等性状表型数据。饲料增重比采用公式对全自动种猪生产性能测定***(FIRE，美国)记录的生长数据进行分析获得，用于表型-基因型关联分析。其中，FCR为饲料增重比；Wa为活体增重量；Wf为饲料消耗量。
(2)基因分型与质量控制
采集1733头公猪的耳组织样品或者血样，提取总DNA，并采用GGP 50k SNP(GeneSeek，US)芯片进行基因分型，获得覆盖全基因组的50705个SNP标记。根据最新版的猪参考基因组(Sscrofall.1)，采用NCBI基因组比对程序(https：//www.ncbi.nlm.nih.gov/)对所有SNP标记的物理位置进行更新。基因组位置未知的SNP不用于关联分析。对于所有常染色体上的SNP标记，利用Plink软件进行质量控制，标准为：个体检出率≥90％；SNP检出率≥90％；小等位基因频率≥0.01；哈迪-温伯格平衡p值≥10-6。对于缺失基因型，采用Beagle软件(version 4.1)进行填充。基于以上质量控制标准，剩余1623头公猪和28289个SNP标记用于关联分析。
(3)统计模型
为了充分利用所有表型数据和基因型数据，本发明公开采用加权的一步法全基因组关联分析法(weighted single step genome-wide association study，wssGWAS)进行全基因组关联分析。该方法首先基于混合模型方程组来估计个体育种值，继而基于育种值模型与标记效应模型的等价关系将育种值转换为标记效应。本发明采用的全基因组关联分析模型如下：
y＝Xb+Za+Wp+e，
其中，y为饲料增重比观测值向量；X，Z和W为设计矩阵；b为固定效应向量(环境、日龄)；为育种值向量；为个体永久环境效应；为残差。H为同时整合系谱和SNP标记的亲缘关系矩阵，其逆矩阵计算公式如下：
其中，A为基于系谱的亲缘关系矩阵；A22为A中有基因型个体对应的分块矩阵；Gω＝0.9G+0.1A22，为基于全基因组SNP标记的亲缘关系矩，Z为小等位基因频率(minor allele frequency，MAF)校正后的基因型矩阵，其中0-2p，1-2p和2-2p分别代表AA，Aa和aa三种基因型，p为小等位基因频率；D为对角线矩阵，表示SNP的权重；pi为第i个标记的小等位基因频率；m为标记数量。
对应上述混合模型，采用AI-REML(average information restricted maximumlikelihood)法估计方差组分，并通过求解混合模型方程组获得育种值。通过迭代的方式获得标记权重，主要步骤如下：
第1步：初始化(t＝1)，D(t)＝I，G(t)＝λZD(t)Z′，
第2步：通过ssGBLUP计算个体育种值；
第3步：通过公式将个体育种值转换为SNP效应，其中为有基因型个体的育种值；
第4步：利用公式计算SNP权重用于下一轮迭代；
第5步：利用公式对SNP权重进行标准化，以保证方差一致；
第6步：利用公式G(t+1)＝λZD(t+1)Z′计算亲缘关系矩阵用于下一轮迭代；
第7步：令t＝t+1，并从第2步开始下一轮迭代。
上述步骤迭代三次，最终获得SNP标记效应，即得到SNP标记效应。将第三轮迭代输出的标记效应作为最终的结果。计算过程主要通过在R统计分析平台编程调用BLUPF90软件来实现，其中AIREMLF90程序用于方差组分估计，BLUPF90程序用于计算育种值，postGSf90用于计算标记效应。
(4)标记筛选
对于所有标记的效应值，取其绝对值画曼哈顿图，展示和筛选大效应的SNP标记。并采用方差分析和多重比较(R统计分析平台)，分析ALGA0093815标记不同基因型群体公猪剩余采食量差异情况。
分析不同基因型公猪剩余采食量
对于所有标记的效应值，取其绝对值画曼哈顿图，展示和筛选大效应的SNP标记(如图2所示)。并采用方差分析和多重比较(R统计分析平台)，分析不同基因型群体公猪剩余采食量差异情况(表1)。
关于ALGA0093815标记基因序列。
本公开鉴定了一个影响公猪剩余采食量的分子标记ALGA0093815，由表1可知，该标记不同基因型公猪的剩余采食量有极显著差异(表1为ALGA0093815标记不同基因型公猪剩余采食量)。
表1 ALGA0093815标记不同基因型公猪剩余采食量
通过检测ALGA0093815标记基因型辅助种猪选育，可通过选留TT纯合种猪进入核心群，降低剩余采食量，有效减少饲料消耗量和养殖成本
由表2可知(表2为ALGA0093815标记不同基因型杜长大商品猪剩余采食量)，在杜长大商品猪群体中，TT基因型猪的剩余采食量比CC基因型显著降低了53.19g(P＝0.05)；所以，该位点多态性有利于降低杜长大商品猪剩余采食量，T是有利于降低剩余采食量的等位基因。
表2 ALGA0093815标记不同基因型杜长大商品猪剩余采食量
本公开的SNP分子遗传标记(ALGA0093815)，参阅Ensembl数据库(http://asia.ensembl.org/Sus_scrofa/Search/New？db＝core)，得到登录号为ALGA0093815的基因片段(RS号为rs80818852),SNP分子遗传标记位于猪17号染色体22568436bp位置，且属于SEL1L2基因的内含子上，该位置为一个T>C突变(突变位点)，该SNP分子遗传标记在突变位点上下游100bp序列如下：
5’-TTCAGTCAGGTCAGGAACAAGGCAGGAGCCAGGGGCAAATTCTTCTAGCTCATGAGTGACTTTAAAGATTAGAAGAGAGCTGCTCTCTAATCCACCTTCCR(T/C)ATCATCCACGTCCATGTGCATGTTTTGAAATAATATGAGTTAAATTTTATATAGGGGGATTTCCCATCATTGCTCAGTGGTAAAGAACCTGACTAGTGTC-3’；R为突变位点，上述序列101位核苷酸处的R是C或T时，即R(T/C)时，导致上述序列多态性；当上述核苷酸序列的第101位核苷酸为T时，猪显著降低剩余采食量，5’-和-3’分别表示核苷酸序列的5’端和3’端。
(上述序列在突变点位为T时的序列如序列表SEQ IDNo.1所示的核苷酸序列),序列表SEQ IDNo.1是本发明筛选得遗传标记(即SNP号为ALGA0093815，RS号为rs80818852)的突变位点上下游100bp的核苷酸序列。
通过检测ALGA0093815标记基因型辅助种猪选育，通过选留TT纯合种猪进入核心群，SNP分子遗传标记的SNP号为ALGA0093815，ALGA0093815标记TT基因型比CC纯种杜洛克公猪个体间剩余采食量显著降低了63.52g(由表1可知)，所以，T是有利于显著降低剩余采食量的等位基因。
通过选留TT纯合种猪进入核心群，SNP分子遗传标记的SNP号为ALGA0093815，ALGA0093815标记TT基因型比CC纯种杜洛克公猪个体间剩余采食量显著降低了63.52g。
6.根据权利要求5所述的猪剩余采食量性状的分子遗传标记在猪剩余采食量性状辅助选择中的应用，其特征在于，所述猪剩余采食量性状辅助选择中的应用方法为通过选留ALGA0093815位点多态性的种猪个体，TT基因型猪的剩余采食量比CC基因型显著降低53.19g
主要参考文献：
1.杨运清，王忠新，尹炳北.比值性状的指数选择方法探讨.畜牧兽医学报，31(3)：198-202(2000)。
2.Koch R M,Swiger L A,Chambers D,et al.Efficiency of feed use in beefcattle.J Anim Sci,22(2):486-494(1963)。
3.Vigors S,Sweeney T,Oshea C J,et al.Pigs that are divergent in feedefficiency,differ in intestinal enzyme and nutrient transporter geneexpression,nutrient digestibility and microbial activity.Animal,10(11):1848-1855(2016)。
4.Hirschhorn,J.N.&Daly,M.J.Genome-wide association studies for commondiseases and complex traits.Nat.Rev.Genet.6,95–108(2005)。
5.De,R.,Bush,W.S.&Moore,J.H.Bioinformatics challenges in genome-wideassociation studies(GWAS).Methods Mol.Biol.1168,63–81(2014)。
6.Purcell,S.et al.PLINK:A Tool Set for Whole-Genome Association andPopulation-Based Linkage Analyses.Am.J.Hum.Genet.813,559–575(2007)。
7.WANG,H.,MISZTAL,I.,AGUILAR,I.,LEGARRA,A.&MUIR,W.M.Genome-wideassociation mapping including phenotypes from relatives withoutgenotypes.Genet Res 94,73–83(2012)。
8.Wang,H.et al.Genome-wide association mapping including phenotypesfrom relatives without genotypes in a single-step(ssGWAS)for 6-week bodyweight in broiler chickens.Front.Genet.5,1–10(2014)。
9.Misztal,I.et al.BLUPF90 and related programs(BGF90).in Proc.7thWorld Congr.Genet.Appl.Livest.Prod.21–22(2002).doi:9782738010520。
序列表
<110> 佛山科学技术学院
<120> 猪剩余采食量性状的分子遗传标记及应用
<141> 2019-06-05
<160> 1
<170> SIPOSequenceListing 1.0
<210> 1
<211> 201
<212> DNA
<213> Sscrofa11.1
<220>
<221> gene
<222> (1)..(201)
<220>
<221> mutation
<222> (101)..(101)
<400> 1
ttcagtcagg tcaggaacaa ggcaggagcc aggggcaaat tcttctagct catgagtgac 60
tttaaagatt agaagagagc tgctctctaa tccaccttcc tatcatccac gtccatgtgc 120
atgttttgaa ataatatgag ttaaatttta tataggggga tttcccatca ttgctcagtg 180
gtaaagaacc tgactagtgt c 201
Claims (6)
1.猪剩余采食量性状的分子遗传标记，其特征在于，所述分子遗传标记为SNP分子遗传标记，SNP分子遗传标记位于猪17号染色体22568436bp位置，且属于SEL1L2基因的内含子上，该位置为一个T>C突变，猪参考基因组为Sscrofa11.1。
2.根据权利要求1所述的猪剩余采食量性状的分子遗传标记，其特征在于，所述SNP分子遗传标记的序列为突变位点的上下游100bp序列。
3.根据权利要求2所述的猪剩余采食量性状的分子遗传标记，其特征在于，所述SNP分子遗传标记的序列为以下所示：
5’-TTCAGTCAGGTCAGGAACAAGGCAGGAGCCAGGGGCAAATTCTTCTAGCTCATGAGTGACTTTAAAGATTAGAAGAGAGCTGCTCTCTAATCCACCTTCCRATCATCCACGTCCATGTGCATGTTTTGAAATAATATGAGTTAAATTTTATATAGGGGGATTTCCCATCATTGCTCAGTGGTAAAGAACCTGACTAGTGTC-3’；R为突变位点，当R为T时，猪显著降低剩余采食量。
4.权利要求3所述的猪剩余采食量性状的分子遗传标记在猪剩余采食量性状辅助选择中的应用。
5.根据权利要求4所述的猪剩余采食量性状的分子遗传标记在猪剩余采食量性状辅助选择中的应用，其特征在于，所述猪剩余采食量性状辅助选择中的应用方法为通过检测ALGA0093815标记基因型辅助种猪选育，通过选留TT纯合猪进入核心群，SNP分子遗传标记的SNP号为ALGA0093815，ALGA0093815标记TT基因型比CC纯种杜洛克公猪个体间剩余采食量显著降低了63.52g，通过检测ALGA0093815标记基因型辅助种猪选育，可通过选留TT基因型纯合猪，降低剩余采食量，有效减少饲料消耗量和养殖成本。
6.根据权利要求4所述的猪剩余采食量性状的分子遗传标记在猪剩余采食量性状辅助选择中的应用，其特征在于，所述猪剩余采食量性状辅助选择中的应用方法为通过选留ALGA0093815位点多态性的猪个体，TT基因型猪的剩余采食量比CC基因型显著降低53.19g，通过检测ALGA0093815标记基因型辅助种猪选育，可通过选留TT基因型纯合种猪，降低剩余采食量，有效减少饲料消耗量和养殖成本。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910489835.0A CN110218799B (zh) | 2019-06-06 | 2019-06-06 | 猪剩余采食量性状的分子遗传标记及应用 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910489835.0A CN110218799B (zh) | 2019-06-06 | 2019-06-06 | 猪剩余采食量性状的分子遗传标记及应用 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110218799A true CN110218799A (zh) | 2019-09-10 |
| CN110218799B CN110218799B (zh) | 2024-03-26 |
Family
ID=67819631
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910489835.0A Active CN110218799B (zh) | 2019-06-06 | 2019-06-06 | 猪剩余采食量性状的分子遗传标记及应用 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110218799B (zh) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111500746A (zh) * | 2020-05-22 | 2020-08-07 | 华中农业大学 | 一种与猪的饲料转化效率相关的snp分子标记 |
| CN112002371A (zh) * | 2020-07-31 | 2020-11-27 | 中国农业科学院北京畜牧兽医研究所 | 一种白羽肉鸡剩余采食量的基因组选择方法 |
| CN113215285A (zh) * | 2021-05-25 | 2021-08-06 | 江西农业大学 | 影响猪生长性状的snp标记 |
| CN113584181A (zh) * | 2021-07-26 | 2021-11-02 | 温氏食品集团股份有限公司 | 一种与猪剩余采食量相关的snp分子标记及其用途 |
| CN113699247A (zh) * | 2021-07-26 | 2021-11-26 | 华南农业大学 | 一种猪1号染色体上与猪剩余采食量相关的snp分子标记及其用途 |
| CN114292927A (zh) * | 2022-03-11 | 2022-04-08 | 佛山科学技术学院 | 与母猪产仔均匀度相关的分子标记及其获得方法和应用 |
| CN114568388A (zh) * | 2022-03-30 | 2022-06-03 | 上杭傲农槐猪产业发展有限公司 | 槐猪的保种和育种方法 |
| CN116042849A (zh) * | 2022-11-09 | 2023-05-02 | 佛山科学技术学院 | 一种用于评估猪采食量的遗传标记及其筛选方法和应用 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104250646A (zh) * | 2013-06-27 | 2014-12-31 | 华中农业大学 | 一种与猪饲料转化效率性状相关的分子标记及检测方法和应用 |
| WO2018098671A1 (en) * | 2016-11-30 | 2018-06-07 | China Agricultural University | A method for crispr library screening |
| CN108559781A (zh) * | 2018-03-28 | 2018-09-21 | 中国农业科学院北京畜牧兽医研究所 | 一种培育高饲料利用效率猪的方法 |
| KR20190034125A (ko) * | 2017-09-22 | 2019-04-01 | 고려대학교 산학협력단 | 모돈의 생애총생산성 형질 예측용 snp 마커 및 이의 용도 |
-
2019
- 2019-06-06 CN CN201910489835.0A patent/CN110218799B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104250646A (zh) * | 2013-06-27 | 2014-12-31 | 华中农业大学 | 一种与猪饲料转化效率性状相关的分子标记及检测方法和应用 |
| WO2018098671A1 (en) * | 2016-11-30 | 2018-06-07 | China Agricultural University | A method for crispr library screening |
| KR20190034125A (ko) * | 2017-09-22 | 2019-04-01 | 고려대학교 산학협력단 | 모돈의 생애총생산성 형질 예측용 snp 마커 및 이의 용도 |
| CN108559781A (zh) * | 2018-03-28 | 2018-09-21 | 中国农业科学院北京畜牧兽医研究所 | 一种培育高饲料利用效率猪的方法 |
Non-Patent Citations (3)
| Title |
|---|
| JEREMY T HOWARD 等: "Genome-wide association study on legendre random regression coefficients for the growth and feed intake trajectory on Duroc Boars", BMC GENET ., vol. 16, no. 59, 30 May 2015 (2015-05-30), pages 1 - 11 * |
| 张跃博等: "杜洛克猪HLCS基因组织表达分析及其编码区多态性与剩余采食量的关联", 《畜牧兽医学报》, no. 06, pages 1 - 3 * |
| 蔡飞翔等: "猪剩余采食量的研究进展", 《中国畜牧杂志》, no. 01, pages 77 - 79 * |
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111500746B (zh) * | 2020-05-22 | 2021-09-17 | 华中农业大学 | 一种与猪的饲料转化效率相关的snp分子标记 |
| CN111500746A (zh) * | 2020-05-22 | 2020-08-07 | 华中农业大学 | 一种与猪的饲料转化效率相关的snp分子标记 |
| CN112002371A (zh) * | 2020-07-31 | 2020-11-27 | 中国农业科学院北京畜牧兽医研究所 | 一种白羽肉鸡剩余采食量的基因组选择方法 |
| CN112002371B (zh) * | 2020-07-31 | 2023-09-26 | 中国农业科学院北京畜牧兽医研究所 | 一种白羽肉鸡剩余采食量的基因组选择方法 |
| CN113215285A (zh) * | 2021-05-25 | 2021-08-06 | 江西农业大学 | 影响猪生长性状的snp标记 |
| CN113215285B (zh) * | 2021-05-25 | 2022-04-26 | 江西农业大学 | 影响猪生长性状的snp标记 |
| CN113699247B (zh) * | 2021-07-26 | 2023-07-14 | 华南农业大学 | 一种猪1号染色体上与猪剩余采食量相关的snp分子标记及其用途 |
| CN113584181A (zh) * | 2021-07-26 | 2021-11-02 | 温氏食品集团股份有限公司 | 一种与猪剩余采食量相关的snp分子标记及其用途 |
| CN113699247A (zh) * | 2021-07-26 | 2021-11-26 | 华南农业大学 | 一种猪1号染色体上与猪剩余采食量相关的snp分子标记及其用途 |
| CN114292927A (zh) * | 2022-03-11 | 2022-04-08 | 佛山科学技术学院 | 与母猪产仔均匀度相关的分子标记及其获得方法和应用 |
| CN114292927B (zh) * | 2022-03-11 | 2022-05-13 | 佛山科学技术学院 | 与母猪产仔均匀度相关的分子标记及其获得方法和应用 |
| CN114568388A (zh) * | 2022-03-30 | 2022-06-03 | 上杭傲农槐猪产业发展有限公司 | 槐猪的保种和育种方法 |
| CN116042849A (zh) * | 2022-11-09 | 2023-05-02 | 佛山科学技术学院 | 一种用于评估猪采食量的遗传标记及其筛选方法和应用 |
| CN116042849B (zh) * | 2022-11-09 | 2023-11-17 | 佛山科学技术学院 | 一种用于评估猪采食量的遗传标记及其筛选方法和应用 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110218799B (zh) | 2024-03-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110218799A (zh) | 猪剩余采食量性状的分子遗传标记及应用 | |
| Lee et al. | Hanwoo cattle: origin, domestication, breeding strategies and genomic selection | |
| WO2023165128A1 (zh) | 一种基于靶向捕获测序的鸡低密度snp液相芯片及其应用 | |
| CN110358840A (zh) | 与剩余采食量相关的tpp2基因的snp分子遗传标记 | |
| CN112002371B (zh) | 一种白羽肉鸡剩余采食量的基因组选择方法 | |
| CN110358838A (zh) | Fa2h基因片段中与猪饲料转化相关的snp遗传标记 | |
| CN110295236A (zh) | 猪饲料转化率的snp分子遗传标记 | |
| CN110358839A (zh) | 与猪饲料转化率相关的gckr基因的snp分子遗传标记 | |
| CN109371143A (zh) | 与猪生长性状相关联的snp分子标记 | |
| CN114941033A (zh) | 一种基于snp位点辅助培育地方优质白羽鸡高产蛋品系的方法 | |
| CN109371144A (zh) | 一种与猪生长性状相关联的snp分子标记 | |
| CN114921561B (zh) | 杜洛克猪全基因组低密度snp芯片及其制备方法和应用 | |
| CN110273006A (zh) | 一种公猪有效***数相关的分子遗传标记 | |
| Qiao et al. | Population genetic structure analysis and identification of backfat thickness loci of Chinese synthetic Yunan pigs | |
| CN111370058B (zh) | 一种基于全基因组snp信息追溯水牛血统来源以及进行基因组选配的方法 | |
| US20240043912A1 (en) | Genomic selection (gs) breeding chip of huaxi cattle and use thereof | |
| CN113549699A (zh) | 一种白羽肉鸡产蛋数的基因组选择方法 | |
| CN114752678B (zh) | 与猪达115公斤体重背膘厚相关的snp分子标记及其应用 | |
| CN113322333B (zh) | 基于全基因组测序筛选的与广西麻鸡体尺和屠宰性状相关的cnv分子标记组合及应用 | |
| CN110714085A (zh) | 一种猪***畸形率相关的分子遗传标记及应用 | |
| CN112176073B (zh) | 鸡屠体性状相关的pros1基因分子标记及应用 | |
| CN111910009B (zh) | 一种影响鸡法氏囊指数的分子标记及其应用 | |
| CN110878363B (zh) | Vipr1基因与鸡种睾丸性状相关性的检测方法及应用 | |
| Masters et al. | Phylogenomic analysis reveals five independently evolved African forage grass clades in the genus Urochloa | |
| CN108707678A (zh) | 与猪种繁殖性状相关的pacsin2基因分子标记及应用 |
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 |