一种筛选奶牛酮病抗性分子标记的方法及其应用技术

本发明专利技术提供一种筛选奶牛酮病抗性分子标记的方法及其应用，属于分子遗传生物技术领域。本发明专利技术以健康和酮病奶牛为样本，进行芯片测序、数据填充、选择信号分析和全基因组关联分析等，筛选出影响中国荷斯坦牛酮病抗性的关键基因APOA1。在APOA1基因启动子区域筛选获得单核苷酸多态性位点g.‑572 A>G，通过关联分析、启动子活性分析等方法，得到GG基因型奶牛血液中β‑羟丁酸浓度显著低于AA基因型，GG基因型个体的启动子活性显著高于AA基因型个体。采用本发明专利技术技术方案，对APOA1基因启动子区域位点g.‑572 A>G进行基因分型，可筛选出酮病抗性强的奶牛个体，达到节省养殖成本的目的。

【技术实现步骤摘要】
一种筛选奶牛酮病抗性分子标记的方法及其应用
本专利技术属于分子遗传生物
，具体涉及一种筛选奶牛酮病抗性分子标记的方法及其应用。
技术介绍
公开该
技术介绍
部分的信息仅仅旨在增加对本专利技术的总体背景的理解，而不必然被视为承认或以任何形式暗示该信息构成已经成为本领域一般技术人员所公知的现有技术。近几十年来，奶业集约化饲养程度不断加强，并一味追求高产，奶牛产后疾病日趋突出。奶牛酮病是由糖类等碳水化合物及挥发性脂肪酸供应不足引发产后奶牛功能失调的代谢性疾病，多发生于泌乳性能良好的高产奶牛(Hibbitt1979)。酮病是奶牛场最常见和最昂贵的代谢疾病之一，患病率大约为30％-40％(ZhangandAmetaj,2017)。患有酮病的奶牛其产奶量、生殖性能和免疫能力降低，罹患其他围产期疾病的风险也更高(McArtetal.,2013)。Duffield指出一例酮病牛的综合各项损失在50到100美元之间(Duffield,2000)。酮病的特征是奶牛乳、尿、血液中的酮体含量增高，临床上乙酰乙酸、β-羟丁酸(BHBA)、丙酮这三种物质总称为酮体。β-羟丁酸是诊断奶牛酮病的最常见酮体(Oetzel,2004)。血浆β-羟丁酸浓度≥1.2mmol/L是判断奶牛是否患酮病的金标准(McArtetal.,2011)。奶牛对酮病具有抗性，且个体间存在较大差异。专利技术人发现，酮病抗性的遗传力低，约为0.02-0.16，若依赖表型的常规育种方法进行选择，难度极大，而分子育种策略对于低遗传力性状更具优势，可以很好地解决这个问题(Koecketal.,2012；Koecketal.,2014)。国外研究者也将酮病抗性作为奶牛育种的热点性状，并开始投入关注(Kroezenetal.,2018；ParkerGaddisetal.,2018)。因此，进一步挖掘影响奶牛酮病抗性形成的功能基因和标记并阐明其分子调控机制，对于高酮病抗性牛的分子育种具有重要意义和价值。
技术实现思路
针对上述现有技术，本专利技术提供一种筛选牛酮病抗性分子标记的方法及其应用，本专利技术针对中国荷斯坦奶牛酮病发生过程的遗传特性，从基因组进化和选择以及疾病易感性角度考虑，采用牛150K基因芯片对患酮病和健康中国荷斯坦奶牛血液样本进行测序，同时整合Fst、XPEHH、XPCLR三种基因组选择信号分析方法以及GEMMA、PLINK全基因组关联分析方法，成功筛选出奶牛酮病抗性基因及其分子标记，并建立相应检测方法，因此具有良好的实际应用之价值。为实现上述目的，本专利技术涉及以下技术方案：本专利技术的第一个方面，提供一种筛选奶牛酮病抗性分子标记的方法，所述方法至少包括：基于SNP芯片检测分析健康和患酮病奶牛DNA样品；基于Fst、XPEHH和XPCLR基因组选择信号分析方法和GEMMA、PLINK全基因组关联分析方法进行分析；筛选受到显著正选择和与酮病显著相关的单核苷酸多态SNPs及候选基因，整合基因功能注释筛选酮病抗性的SNP分子标记。进一步的，基于上述方法，本专利技术筛选鉴定出APOA1基因为奶牛酮病抗性的关键基因，因此可作为奶牛酮病抗性分子标记；更进一步的，奶牛酮病抗性分子标记还包括位于该基因上1个单核苷酸多态位点：g.-572A>G。经试验证明，APOA1基因GG型奶牛血液中β-羟丁酸浓度显著低于AA基因型。并且荧光素酶活性检测，APOA1基因GG型的启动子活性显著高于AA型个体的启动子活性。因此APAO1基因核心启动子多态性位点g.-572A>G可实现对牛的酮病抗性进行评估。其中，所述奶牛为中国荷斯坦奶牛。本专利技术的第二个方面，提供上述方法在如下1)-9)中任一种中的应用：1)鉴定或辅助鉴定待测奶牛的酮病抗性；2)制备鉴定或辅助鉴定待测奶牛的酮病抗性的产品；3)鉴定或辅助鉴定待测奶牛为酮病抗性奶牛；4)制备鉴定或辅助鉴定待测奶牛为酮病抗性奶牛的产品；5)奶牛选育；6)选育奶牛酮病抗性品种；7)制备选育奶牛酮病抗性品种的产品；8)鉴定或辅助鉴定待测奶牛的酮病抗性性状；9)制备鉴定或辅助鉴定奶牛的酮病抗性性状的产品。其中，所述奶牛为中国荷斯坦奶牛。本专利技术的又一具体实施方式中，所述应用方式包括：提取不同奶牛个体的血液DNA；鉴定带有酮病抗性分子标记的个体。所述分子标记为g.-572A>G。与现有技术相比，上述技术方案的有益效果是：1)上述技术方案的检测方法具有独创性，针对牛酮病发生过程的遗传特性，从基因组进化和选择以及疾病易感性角度，以患酮病牛群组作为参考组，健康牛群组作为实验组，整合多种基因组正向选择信号分析技术，分别以奶牛患酮病或健康二元性状为表型以及以奶牛血液中β-羟丁酸含量作为表型，采用全基因组关联分析方法和策略，高效准确筛选与牛酮病发生相关的关键基因和分子标记，方法设计合理，依据关键基因和标记设计的检测方法，具有准确性高、应用操作简便的特点；2)上述技术方案提供了APOA1基因中与酮病抗性相关的一个单核苷酸多态位点，通过检测所述位点的基因型，可以方便的对酮病抗性进行评价，为奶牛的酮病抗性提供有效的评价依据，为家畜养殖企业节省经济成本。3)基于上述基因及位点表达差异的确认，用于检测APOA1基因的试剂、用于检测APOA1基因单核苷酸多态位点的试剂及引物可用于制备奶牛酮病抗性检测的试剂盒，为奶牛酮病抗性提供使用便利的产品，因此具有良好的实际应用之价值。附图说明构成本专利技术的一部分的说明书附图用来提供对本专利技术的进一步理解，本专利技术的示意性实施例及其说明用于解释本公开，并不构成对本公开的不当限定。图1为本专利技术实施例1中筛选奶牛酮病抗性分子标记的方法流程图；图2为本专利技术实施例1中采用XPCLR正向选择信号分析和PLINK全基因组关联分析筛选到奶牛酮病抗性分子标记APOA1基因；其中，A为XPCLR正向选择信号分析；B为PLINK全基因组关联分析；图3为本专利技术实施例2中APOA1基因及其基因组结构示意图，以及SNPg.-572A>G在基因组中位置信息。图4为本专利技术实施例3中APOA1基因中不同基因型的启动子的活性。具体实施方式应该指出，以下详细说明都是例示性的，旨在对本申请提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本申请所属
的普通技术人员通常理解的相同含义。需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本申请的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。正如
技术介绍
所介绍的，专利技术人在研究中发现牛虽然具有相似的遗传背景，其血液中β-羟基丁酸含量的表型以及患酮病与否却存在显著的差本文档来自技高网...

【技术保护点】
1.一种筛选奶牛酮病抗性基因的方法，其特征在于，所述方法至少包括：/n基于SNP芯片检测分析健康奶牛和患酮病奶牛DNA样品；/n基于Fst、XPCLR和XPEHH基因组选择信号分析方法和GEMMA、PLINK全基因组关联分析方法进行分析；筛选受到显著正选择和与酮病显著相关的单核苷酸多态SNPs以及候选基因，整合基因功能注释筛选酮病抗性的SNP分子标记。/n

【技术特征摘要】
1.一种筛选奶牛酮病抗性基因的方法，其特征在于，所述方法至少包括：
基于SNP芯片检测分析健康奶牛和患酮病奶牛DNA样品；
基于Fst、XPCLR和XPEHH基因组选择信号分析方法和GEMMA、PLINK全基因组关联分析方法进行分析；筛选受到显著正选择和与酮病显著相关的单核苷酸多态SNPs以及候选基因，整合基因功能注释筛选酮病抗性的SNP分子标记。


2.如权利要求1所述的方法，其特征在于，所述方法包括：
S1.健康奶牛和患酮病奶牛样品的采集及DNA提取；
S2.SNP芯片检测分析；
S3.基因型填充；
S4.Fst、XPCLR、XPEHH基因组选择信号分析；
S5.GEMMA、PLINK全基因组关联分析；
S6.基于候选基因的筛选策略筛选候选基因；
S7.鉴定位于候选基因内受选择的SNPs；
其中，步骤S4与S5没有先后顺序之分。


3.如权利要求2所述的方法，其特征在于，步骤S1具体方法包括：
S1.1对奶牛血液中β-羟丁酸浓度进行测定；
S1.2采集健康奶牛和患酮病奶牛的血液，提取血液组织中的DNA。


4.如权利要求2所述的方法，其特征在于，步骤S2具体方法包括：
S2.1使用SNP芯片对DNA样品进行分析，并进行基因分型；
S2.2对SNP数据进行过滤，对剩余符合要求的SNP做进一步分析；
S2.3为每条染色体构建单倍型；
S2.4获取数据库中50K基因芯片数据；
S2.5以150K基因芯片为参考，对50K基因芯片数据进行基因型填充。


5.如权利要求2所述的方法，其特征在于，步骤S3具体方法包括：
S3.1选择信号Fst和hapFLK基因组扫描和局部进化树的构建：对奶牛品种获得的所有数据进行Fst分析；
S3.2估算患酮病奶牛和健康奶牛品种之间的XPCLR值；
S3.3估算患酮病奶牛和健康奶牛品种之间的XPEHH值。


6.如权利要求2所述的方法，其特征在于，步骤S4具体方法包括：
S4.1采用GEMMA单变量线性混合模型，以奶牛患酮病或健康二元性状为表型进行全基因组关联分析；
S4.2以奶牛血液中β-羟丁酸含量作为表型，利用多元回归模型，以奶牛的胎次、年龄为协变量进行全基因组关联分析；<...

【专利技术属性】
技术研发人员：黄金明，白佳琛，鞠志花，姜强，王金鹏，张亚冉，王秀革，魏晓超，高亚平，杨春红，刘文浩，李荣岭，高运东，
申请(专利权)人：山东省农业科学院奶牛研究中心，山东奥克斯畜牧种业有限公司，
类型：发明
国别省市：山东;37