CN114107518A

CN114107518A - 一种猪eepd1基因第八外显子变异作为肌内脂肪含量性状的遗传标记及应用

Info

Publication number: CN114107518A
Application number: CN202111443176.0A
Authority: CN
Inventors: 吴俊静; 孙华; 梅书棋; 彭先文; 乔木; 周佳伟; 董斌科; 李梓芃; 张宇; 徐忠; 宋忠旭; 赵海忠; 李良华; 朱昌友
Original assignee: Institute of Animal Science and Veterinary of Hubei Academy of Agricultural Sciences
Current assignee: Institute of Animal Science and Veterinary of Hubei Academy of Agricultural Sciences
Priority date: 2021-11-30
Filing date: 2021-11-30
Publication date: 2022-03-01
Anticipated expiration: 2041-11-30
Also published as: CN114107518B

Abstract

一种猪EEPD1基因第八外显子变异作为肌内脂肪含量性状的遗传标记及应用。本发明提供了一种猪EEPD1基因核苷酸序列第123857bp处的一个C>G碱基突变和第124200处的一个A>G碱基突变作为猪肌内脂肪含量性状相关的分子遗传标记。本发明通过制备得到一种与肌内脂肪含量性状相关的遗传标记，其核苷酸序列如序列表SEQ ID NO：1和SEQ ID NO：2所示，在它们第241bp处存在一个C>G碱基突变，在584bp处存在一个A>G碱基突变，两个突变位点高度连锁，可显著影响猪肌内脂肪含量性状。本发明还提供了该遗传标记的制备方法及其在猪肌内脂肪含量性状关联分析中的应用。本发明为猪肌内脂肪含量性状的分子标记辅助育种提供了一个新的遗传标记。

Description

一种猪EEPD1基因第八外显子变异作为肌内脂肪含量性状的遗传标记及应用

技术领域

本发明属于家畜遗传标记筛选技术领域，具体涉及一种猪EEPD1基因第八外显子变异作为肌内脂肪含量性状的遗传标记及应用。

背景技术

随着社会经济发展和人民生活水平的提高，高品质、多元化、个性化消费需求的优质猪肉市场需求激增，但是当前综合性状优良、特色鲜明、满足中国人口味的猪肉产品不多，猪肉产品市场供需不匹配问题逐渐显现，因此，猪肉品质遗传改良已越来越被重视。肌内脂肪含量(Intramuscular Fat,IMF)是影响肉品质性状的一项重要指标，它与猪肉的嫩度、风味、多汁性存在显著相关(Lu etal.2008)。然而，猪肉质性状较难活体测量，很难用常规育种技术进行遗传改良，因此寻找控制猪肉质性状(尤其是IMF)的关键分子遗传标记，应用于该性状的分子辅助育种，对提高猪肉品质性能具有重要意义。

EEPD1(Endonuclease/Exonuclease/Phosphatase Family Domain Containing1)编码核酸内切酶/核酸外切酶/磷酸酶家族结构域包含蛋白1，主要在人类DNA损伤过程中修复应激复制叉中起作用(Kim etal.,2017)。Ji等(2012)在对鸡脂肪组织的转录组分析也发现，鸡EEPD1基因表达量的降低与禁食和胰岛素中和能力均相关，表明其对脂肪酸合成具有负调控作用。Li等(2020)利用高通量转录组测序技术发现EEPD1基因在荷斯坦公牛和阉公牛的腰最长肌中差异表达，且牛EEPD1基因与牛肉品质性状相关。而有关猪EEPD1基因对肉质性状的影响还未见报道。

申请人前期利用RNA-seq技术对比分析了IMF极端个体背最长肌组织基因表达差异，发现EEPD1基因显著差异表达。因此，申请人扩增了猪EEPD1基因的部分核苷酸序列，利用该序列进行了多态性变异位点的筛选鉴定及与IMF性状的关联分析研究，以期获得与猪IMF性状相关的新分子标记。

发明内容

本发明的技术目的在于扩增猪EEPD1基因的部分第七内含子和第八外显子DNA序列，以寻找该基因序列突变位点多态性，获得一种与猪肌内脂肪含量性状相关的遗传标记，并利用该遗传标记作为猪标记辅助选择中。

为实现上述技术目的，本发明提供了一种猪EEPD1基因第八外显子变异作为肌内脂肪含量性状的遗传标记，它是猪EEPD1基因的部分DNA序列，其核苷酸序列如SEQ ID NO：1和SEQ ID NO：2所示，在该序列第241碱基处存在一个C>G碱基突变，在该序列第584碱基处存在一个A>G碱基突变。

本发明还提供了所述遗传标记在猪肌内脂肪含量性状标记辅助选择中的应用。

优选地，所述的应用是所制备的遗传变异标记在猪IMF性状的关联分析中的应用或者是所制备的单倍型标记在猪IMF性状的关联分析中的应用。

本发明的有益效果是：本发明为猪IMF性状的分子标记辅助育种提供了一个新的遗传标记。

附图说明

图1是本发明的总体技术流程图；

图2是本发明中猪EEPD1基因PCR扩增产物琼脂糖凝胶电泳图；

图3是本发明中猪EEPD1基因第八外显子g.123857C>G突变位点SEQ ID NO：4引物反向测序的测序图谱；

图4是本发明中猪EEPD1基因第八外显子g.124200A>G突变位点SEQ ID NO：4引物反向测序的测序图谱；

图5是猪EEPD1基因片段核苷酸序列直观图(与SEQ ID NO：1和SEQ ID NO：2所示的序列表相对应)。

具体实施方式

下面结合附图1-5，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。

如未特别说明，以下实施例中所采用的方法和参数均为本领域常规的方法和参数。

需要说明的是：

序列表SEQ ID NO：1是扩增的一个“硒都黑猪”EEPD1基因部分第七内含子和第八外显子的核苷酸序列，序列长度为674bp。在该序列第241碱基处存在一个C>G碱基的替换，在该序列第584碱基处存在一个A>G碱基的替换(在序列表SEQ ID NO：1中显示的结果是已经替换的碱基)；

序列表SEQ ID NO：2是扩增的一个“硒都黑猪”EEPD1基因部分第七内含子和第八外显子的核苷酸序列，序列长度为674bp，在该序列第241碱基处存在一个G>C碱基的替换，在该序列第584碱基处存在一个G>A碱基的替换(在序列表SEQ ID NO：1中显示的结果是已经替换的碱基)；

序列表SEQ ID NO：3是扩增SEQ ID NO：1和SEQ ID NO：2特异基因片段所用的正向引物序列；

序列表SEQ ID NO：4是扩增SEQ ID NO：1和SEQ ID NO：2特异基因片段所用的反向引物序列；

图2中泳道M是DNA Marker，泳道1、2是扩增的SEQ ID NO：1和SEQ ID NO：2特异基因片段；

图3中从上往下依次是CC纯合型、GG纯合型和CG杂合型；

图4中从上往下依次是CC纯合型、GG纯合型和CG杂合型；

图5中加粗带框字母代表存在突变处，突变位置分为位于该序列的第241碱基和第584碱基处，带下划线序列代表引物位置，带灰色背景的表示EEPD1基因第八外显子区域。

本发明通过以下技术方案实现：

本发明获得了一种猪肌内脂肪含量性状的遗传标记，它们分别是猪EEPD1基因第八外显子的部分DNA片段，其核苷酸序列如序列表SEQ ID NO：1(长度为674bp)和SEQ IDNO：2(长度为674bp)所示。在SEQ ID NO：1和SEQ ID NO：2第241碱基处存在一个C>G碱基突变(如图3所示)，在该序列第584碱基处存在一个A>G碱基突变(如图4所示)。

申请人提供了一种扩增猪EEPD1基因部分DNA片段的引物对，该引物对的DNA序列如SEQ ID NO：3和SEQ ID NO：4所示。

申请人提供了一种制备猪肌内脂肪含量性状相关基因EEPD1遗传标记的方法，其步骤如下：

参考GenBank数据库中猪EEPD1基因(登录号NC_010460.4)的DNA序列设计特异引物(该引物的序列如序列表SEQ ID NO：3和SEQ ID NO：4所示)，以硒都黑猪基因组DNA混合池为模板进行PCR扩增，对PCR产物回收并进行DNA测序分析，得到如序列表SEQ ID NO：1和SEQ ID NO：2所示的基因片段。根据测序结果筛查DNA序列变异，在241bp处存在一个C>G碱基突变(即等位基因突变)，在584bp处存在一个A>G碱基突变(即等位基因突变)，所述的遗传标记如图3和图4所示，利用该遗传标记在对硒都黑猪群中基因型与IMF性状进行了关联分析检测。

更详细技术细节如下列实施例所述。

实施例1猪EEPD1基因片段的获得及多态性检测方法的建立

1.猪基因组DNA的提取

本发明的试验猪品种为硒都黑猪，样本来源于湖北天之力优质猪育种有限公司。硒都黑猪是以恩施黑猪、梅山猪、湖北白猪等为基础，由湖北省农业科学院畜牧兽医研究所、湖北天之力优质猪育种有限公司等单位历经12年培育而成的黑猪新品种。猪基因组DNA的提取采用北京百泰克生物技术有限公司生产的基因组DNA试剂盒(按该试剂盒说明书进行操作)提取，具体步骤如下所述：

(1)采取猪背最长肌组织，放入2mL的离心管中，加入200μL裂解液TL，用枪头吹打均匀；

(2)加入20μL蛋白酶K(20mg/ml)，剧烈颠倒充分混匀，55℃水浴锅中消化过夜；

(3)加入200μL结合液CB(试剂盒自带)，充分颠倒混匀，70℃放置10min；

(4)冷却后加入100μL异丙醇，剧烈颠倒充分混匀；

(5)用1mL的枪头吸取上述混合物，加入吸附柱AC中，10000rpm离心30s，倒掉收集管中的废液；

(6)加入500μL抑制物去除液IR(该试剂盒自带)，12000rpm离心30s，弃废液；

(7)加入700μL漂洗液WB(试剂盒自带)，12000rpm离心30s，倒掉废液；

(8)重复操作步骤7；

(9)将吸附柱AC放回收集管中，12000rpm离心2min，尽量去除漂洗液，以免残留的乙醇抑制下游反应；

(10)取出吸附柱AC，放入一个干净的离心管中，向吸附膜的中间部位加50-100μL洗脱缓冲液EB(该试剂盒自带)，室温放置3-5min，12000rpm离心1min，将溶液收集到离心管中；

(11)对提取出的DNA的浓度及质量进行检测后置于-20℃下保存备用。

剩下肌肉样本封袋保存于干冰中，送至华中农业大学农业部种猪质量监督检验测试中心(武汉)用于IMF的测定。

2.猪EEPD1基因片段的获得

(1)PCR扩增

根据猪EEPD1基因的基因组序列(GenBank登陆号：NC_010460.4)设计以下引物对：

正向引物EEPD1-E8-F：5'ACAGGCAGACAGGATGGC 3'，

反向引物EEPD1-E8-R：5'CGGGTGGTCAGAAACGAG 3'。

利用上述引物在37个硒都黑猪混合基因组DNA池中进行PCR扩增，PCR反应体系50μL，体系中各组分的浓度为100ng模板DNA、10×buffer(含Mg2+)4μL、上述上下游引物各0.5μM、2.5μM dNTPs、1U TaqDNA聚合酶。

PCR的运行程序为：95℃预热2min；94℃变性20s，57℃退火40s，72℃延伸60s，共35个循环；72℃延伸10min；4℃保存。PCR产物用1％琼脂糖凝胶电泳检测。

(2)PCR产物纯化

上述PCR产物用上海生工生物工程有限公司的Gel Extraction Kit试剂盒进行纯化(按照该试剂盒的说明书操作)，具体步骤如下：首先从琼脂糖凝胶上切下含目的片段的凝胶，放入1.5mL离心管，加入400μL溶胶液，50-60℃水浴至胶彻底融化，加热融胶时，每2min混匀一次，冷却至室温；将离心柱放入收集管中，把混合液移至离心柱，室温放置2min；12000r/min离心1min，此时DNA被吸附到柱上；倒掉收集管中废液，将离心柱放入同一个收集管中，加入700μL洗脱液，12000r/min离心1min；倒掉收集管中的废液，12000r/min离心1min；将离心柱放入一预先准备好的灭菌1.5mL离心管中，加入40μL洗脱液或双蒸水(pH＞7.0)，室温或37℃放置2-3min；12000r/min离心1min，离心管中的液体即为回收的DNA片段。

3.猪EEPD1基因片段变异位点的获得

将上述回收获得的DNA片段送到武汉奥科鼎盛生物科技有限公司采用ABI3730XL测序仪测序，发现了2个单碱基突变位点(图3和图4)，分别位于EEPD1基因组核苷酸序列(即全序列)(GenBank NC_010460.4)的第123857bp处的C>G突变，以及第124200bp处的A>G突变(对应于本发明的EEPD1基因片段的突变位点分别是：在SEQ ID NO：1是241bp处存在一个C>G碱基突变(即等位基因突变)，在584bp处存在一个A>G碱基突变(即等位基因突变)；在SEQID NO：2是241bp处存在一个G>C碱基突变(即等位基因突变)，在584bp处存在一个G>A碱基突变(即等位基因突变)。

4.分子标记基因分型

将待检测个体的DNA样本作为模板，按照以上步骤2所述方法进行猪EEPD1基因序列片段的扩增，将获得的PCR纯化产物直接送武汉奥科鼎盛生物科技有限公司进行测序，直接从测序结果中读取基因分型结果。

实施例2本发明制备的遗传变异标记在猪IMF性状的关联分析中的应用

申请人在218头硒都黑猪群(来自湖北天之力优质猪育种有限公司)检测本发明制备的两个分子遗传标记与猪IMF性状的关联性。

采用实施例1建立的PCR直接测序方法来进行基因分型检测，采用SPSS统计软件(Statistical Package for the Social Sciences,Version 17.0)一般线性模型GLM进行统计分析。所采用模型为：Y_ijklm＝μ+G_i+A_j+X_k+S_l+e_ijklm，其中：Y_ijklm表示IMF值；μ表示群体均值；G_i表示基因型效应；A_j表示年季效应；X_k表示性别效应；S_l表示父本效应；e_ijklm表示随机残差效应。结果以最小二乘均值±标准误表示，P<0.05判定为差异显著；P<0.01判定为差异极显著。

关联分析结果如表1所示，发现g.123857C>G和g.124200A>G均极显著影响IMF性状(P<0.01)。g.123857C>G位点CG基因型个体的IMF极显著高于CC基因型个体(P<0.01)，且GG基因型个体的IMF也显著高于CC基因型个体(P<0.05)，G等位基因是IMF性状的一个优势等位基因。g.124200A>G位点AG型个体的IMF极显著高于AA型个体(P<0.01)，GG型个体的IMF也存在高于AA型个体的趋势，因此G等位基因是IMF性状的一个优势等位基因。

表1猪EEPD1基因g.123857C>G和g.124200A>G与IMF性状关联分析结果

表1注：肩标不同小写字母表示同一列的数据之间差异显著，P<0.05；*基因型：g.123857C>G位点基因型0表示CC型，1表示CG型，2表示GG型；g.124200A>G位点基因型0表示AA型，1表示AG型，2表示GG型。

实施例3本发明制备的单倍型标记在猪IMF性状的关联分析中的应用

通过计算g.123857C>G和g.124200A>G的连锁不平衡关系，发现两突变位点存在一定程度的连锁(D＝0.83)，因此进一步分析了两个位点的组合单倍型分析。发现在硒都黑猪群体中，存在CA、GA、CG、GG四种单倍型，其基因频率分别为79.72％、1.89％、5.66％、12.74％。其中，CA单倍型比例最高，其次是GG单倍型，GA和CG型比较稀少。

在此基础上，进一步分析了单倍型所致组合基因型对IMF性状的影响，统计分析模型参考应用实例2。结果见表2所示，发现GA/GG组合基因型个体的IMF最高，显著高于其他六种基因型个体(P<0.05)。因此，在IMF性状的选育改良中，建议优先选择GA/GG组合基因型个体留种。

表2猪EEPD1基因g.123857C>G和g.124200A>G组合基因型与IMF性状关联分析结果

表2注：肩标不同小写字母表示同一列的数据之间差异显著，P<0.05；不同大写字母表示同一列的数据之间差异极显著，P<0.01。

综上，EEPD1基因全序列第123857bp处的G等位基因以及第124200bp处的G等位基因是IMF性状的优势等位基因，但在联合筛选中，应优先保留携带CA/GG组合基因型的个体，有利于快速提高群体的IMF性状。

序列表

<110> 湖北省农业科学院畜牧兽医研究所

<120> 一种猪EEPD1基因第八外显子变异作为肌内脂肪含量性状的遗传标记及应用

<141> 2021-11-22

<160> 4

<170> SIPOSequenceListing 1.0

<210> 1

<211> 674

<212> DNA

<213> Artificial Sequence

<400> 1

acaggcagac aggatggctg ggggggtctt gtcggggagc ccaccccacc gctgccggct 60

ggcctatgac ctcccagccc ccagcctgac ctctgggtcc ctctcccctc tgcctccagg 120

tcactgggct gtcgtgagag aaggcctcac gaacccttgg attccagaca actggtcctg 180

gggcggagtg gcttcggaac actgccccgt gctcgccgag ttttacacgg agaaggacgg 240

cagcaggaag gagggccctc ggaacagcag cggggtcacc ctagagccaa gtgaggccaa 300

cattaagcac gagcgatgat gacgccaagt cggttttcca cccgcaaccc cggaggcgag 360

gagctcgccc tccgggtgaa gatgcagagc ccgagctccg gcctcagccc cggcctcgcc 420

cgcctgcgcc ctgtgtggac ctctctggac tctcaaggaa ggagcccgtc cggccctcga 480

ggagccggat gctcacatca gaagcccagg ggagacccgg aatatcttct ctctgccgat 540

gccgccaacc ggggtgcccc cggcggcccc cccgctcctg gacagcagcg ccgttggagg 600

ctcaaggcgg agagagacgc tccgccagca ggacctgccc ggcagtgccc ggagccctcg 660

tttctgacca cccg 674

<210> 2

<211> 674

<212> DNA

<213> Artificial Sequence

<400> 2

acaggcagac aggatggctg ggggggtctt gtcggggagc ccaccccacc gctgccggct 60

ggcctatgac ctcccagccc ccagcctgac ctctgggtcc ctctcccctc tgcctccagg 120

tcactgggct gtcgtgagag aaggcctcac gaacccttgg attccagaca actggtcctg 180

gggcggagtg gcttcggaac actgccccgt gctcgccgag ttttacacgg agaaggacgg 240

gagcaggaag gagggccctc ggaacagcag cggggtcacc ctagagccaa gtgaggccaa 300

cattaagcac gagcgatgat gacgccaagt cggttttcca cccgcaaccc cggaggcgag 360

gagctcgccc tccgggtgaa gatgcagagc ccgagctccg gcctcagccc cggcctcgcc 420

cgcctgcgcc ctgtgtggac ctctctggac tctcaaggaa ggagcccgtc cggccctcga 480

ggagccggat gctcacatca gaagcccagg ggagacccgg aatatcttct ctctgccgat 540

gccgccaacc ggggtgcccc cggcggcccc cccgctcctg gacggcagcg ccgttggagg 600

ctcaaggcgg agagagacgc tccgccagca ggacctgccc ggcagtgccc ggagccctcg 660

tttctgacca cccg 674

<210> 3

<211> 18

<212> DNA

<213> Artificial Sequence

<400> 3

acaggcagac aggatggc 18

<210> 4

<211> 18

<212> DNA

<213> Artificial Sequence

<400> 4

cgggtggtca gaaacgag 18

Claims

1.一种猪EEPD1基因第八外显子变异作为肌内脂肪含量性状的遗传标记，它是猪EEPD1基因的部分DNA序列，其核苷酸序列如SEQ ID NO：1和SEQ ID NO：2所示，在该序列第241碱基处存在一个C>G碱基突变，在该序列第584碱基处存在一个A>G碱基突变。

2.如权利要求1所述遗传标记在猪肌内脂肪含量性状标记辅助选择中的应用。

3.如权利要求2所述的应用，其特征在于是所制备的遗传变异标记在猪IMF性状的关联分析中的应用或者是所制备的单倍型标记在猪IMF性状的关联分析中的应用。