一种湖羊SLC27A6基因SNP的检测方法及其在肉质性状早期筛选中的应用与流程

一种湖羊slc27a6基因snp的检测方法及其在肉质性状早期筛选中的应用
技术领域
1.本发明属于分子遗传学检测及育种领域，具体涉及利用dna混池测序检测出的湖羊slc27a6基因单核苷酸多态性(single nucleotide polymorphism,snp)位点，并结合荧光多重酶连接反应(improved multiplex ligation detection reaction,imldr)技术对该位点进行基因分型及利用分型结果与脂肪酸百分比进行关联性分析。


背景技术：

2.snp广泛存在于基因组中，大约每300bp的序列长度就会出现一个snp。snp根据突变发生位置可分为编码区snp和非编码区snp，前者主要通过影响蛋白质的结构，进而对表型产生影响。目前筛选重要性状相关的snp在遗传改良等方面具有重要作用。imldr技术是一种新型高通量snp基因分型技术，该技术具有分型准确率高、成本低、速度快、可同时检测多个snp位点等优势。
3.脂肪酸是影响肉品质重要的属性之一，其含量和组成决定了脂肪的质地特性以及肌肉的氧化稳定性，影响肉的风味、嫩度、保质期等，同时脂肪酸是机体的主要能量来源。脂肪酸按照碳链饱和程度分为饱和脂肪酸(saturated fatty acid,sfa)、单不饱和脂肪酸(monounsaturated fatty acids,mufa)、多不饱和脂肪酸(polyunsaturated fatty acids,pufa)。sfa可致使羊肉产生膻味，并能增加人体患心血管疾病的风险，sfa中的己酸(c6:0)、辛酸(c8:0)、癸酸(c10:0)是导致羊肉产生膻味的主要脂肪酸，其中c10:0占主导作用；mufa与肉香味和整体可接受程度呈正相关关系，mufa还具有降低血糖、血脂、胆固醇水平等功能，油酸(c18:1n9c)是机体内含量最高的一种mufa，可降低低密度脂蛋白(low density lipoprotein,ldl)和胆固醇水平；pufa中的n-6pufa和n-3pufa是人体的必须脂肪酸，与人体健康紧密相关，且人体自身合成，只能从食物中获取，pufa对炎症、高血压、缺血性中风、外周动脉疾病和动脉粥样硬化具有有利影响。
4.脂肪酸转运蛋白6(solute carrier family 27 member 6,slc27a6)是一种膜相关的脂肪酸结合蛋白，促进长链脂肪酸的吸收，同时还参与脂肪酸的跨膜转运和脂质代谢。massimo等人研究发现，slc27a6作为牛乳腺组织中的一种跨膜蛋白，可将游离脂肪酸(c16:0、c18:1n9c、c20:4n6、c24:0)从血液转运到牛乳腺细胞。沈子亮等对奶牛乳腺上皮细胞中的slc27a6干扰后，slc27a6的mrna和蛋白质表达量均减少，同时脂肪酸转运及脂肪氧化相关基因(acsl4、cd36、cpt1a)的mrna表达量也减少，pparg、fabp3、fads2基因的mrna表达量升高，此外，乳腺上皮细胞中c16:0、c18:0的含量降低，c18:1n9c、c20:4n6的含量升高。nafikov等人研究发现，slc27a6基因中的snps(g.70916t/c、g.67847t/c、g.16048g/c、g.15975t/c、g.15740a/c、g.390c/t、g.242a/t)与牛奶中的sfa、ufa、mufa及sfa:ufa的值显著相关，同时slc27a6基因中的单倍型h3(g1c2g3t4c5c6t7)与较低的sfa及较高的ufa显著相关。
5.尽管slc27a6在脂肪酸运输及脂质代谢中发挥重要作用，但目前尚无绵羊slc27a6
基因中的snp影响脂肪酸表型的相关报道，缺乏能够用于对绵羊脂肪酸表型进行选育，并全面提升肉品质的分子标记。


技术实现要素：

6.本发明的目的在于提供一种湖羊slc27a6基因snp的检测方法及其在肉质性状早期筛选中的应用。
7.本发明的目的是通过以下技术方案来实现的：
8.一种检测湖羊slc27a6基因单核苷酸多态性的方法，包括以下步骤：
9.以提取的湖羊基因组dna为模板进行pcr扩增，对扩增产物包含的单核苷酸多态性位点进行基因分型，所述单核苷酸多态性位点为位于slc27a6基因上的突变位点g.504g》c(参考序列为nc_056058.1)。
10.优选的，所述pcr扩增的引物为：
11.上游引物：5
’‑
attgagttgtaagacagggcaca-3’12.下游引物：5
’‑
tacgctgcttaggaggaatgag-3’。
13.优选的，所述pcr扩增采用的反应体系包括30ng/μl模板1μl以及10pmol/μl上游引物和10pmol/μl下游引物各0.5μl；pcr扩增采用的反应程序为：95℃预变性4min；94℃变性30s，63.8℃退火60s，72℃延伸90s，共34个循环；72℃再延伸9min。
14.优选的，所述基因分型采用imldr技术。
15.优选的，所述基因组dna采用酚-氯仿法提取自睾丸组织。
16.优选的，所述突变位点是通过dna混池测序确定的。
17.一种检测湖羊slc27a6基因单核苷酸多态性的试剂盒，该试剂盒包括用于扩增slc27a6基因单核苷酸多态性位点的引物(例如上述上游引物和下游引物)，所述单核苷酸多态性位点为位于slc27a6基因上的突变位点g.504g》c(参考序列为nc_056058.1)。
18.上述检测湖羊slc27a6基因单核苷酸多态性的方法在绵羊标记辅助选择育种中的应用。
19.绵羊slc27a6基因突变位点g.504g》c(参考序列为nc_056058.1)及其检测试剂在绵羊标记辅助选择育种中的应用。
20.优选的，具有gg基因型的个体在肉质性状上较优(例如gg基因型的个体c20:3n6、c20:3n3、c22:6n3的百分比在三种基因型中最高，c22:0百分比在三种基因型中最低)。
21.本发明的有益效果体现在：
22.本发明通过对湖羊slc27a6基因上一个新的snp位点(位于基因突变位点g.504g》c)进行基因分型，以及根据关联性分析所揭示的该位点与脂肪酸百分比的显著相关性，发现了slc27a6基因上这一多态性位点中可以作为影响脂肪酸组成及含量的分子标记(snp标记)，为提高绵羊(例如湖羊)肉品质的分子育种和肉质改良提供重要的科学依据，从而加快优异肉品质的绵羊新品系的选育。
23.进一步的，本发明通过pcr高温(94℃)变性、低温(63.8℃)退火，使提取自湖羊睾丸组织的基因组dna双链解开，同时有利于延伸；而采取34次循环，能扩增出足够数量的dna，以便用于dna混池测序。
附图说明
24.图1为湖羊slc27a6基因突变位点g.504g》c附近区域的pcr扩增产物的电泳图。
25.图2为湖羊slc27a6基因突变位点g.504g》c的测序峰图。
具体实施方式
26.下面结合附图和实施例对本发明作进一步详细说明。所述实施例仅用于解释本发明，而非对本发明的保护范围的限制。
27.1、利用dna混池测序检测获得湖羊slc27a6基因snp位点
28.s1样品选取
29.于2018年8月至2021年8月共采集来自民勤德福农业科技有限公司(甘肃省武威市民勤县)相同饲养条件下的1037只(5批次)6月龄健康湖羊的睾丸组织及背最长肌组织，并测定背最长肌组织的脂肪酸组成及含量。
30.s2构建dna混池
31.a)使用酚-氯仿法提取1037只湖羊个体睾丸组织基因组dna，具体步骤如下：
32.i)取10mg研磨好的睾丸组织粉末于2ml离心管中；
33.ii)分别加入750μl组织裂解液和10μl的10mg/ml蛋白酶k，颠倒混匀，在55℃水浴锅中孵育15个小时；
34.iii)孵育完成后冷却至室温，加入800μl的tris饱和酚，冰上摇动15min，12000r/min离心10min，转移800μl上清液至新的2ml离心管中；
35.iv)分别加入500μl的tris饱和酚和500μl的氯仿，冰上摇动15min，12000r/min离心10min，将上清液转移至新的2ml离心管中；
36.v)加入1ml的氯仿，冰上摇动15min，12000r/min离心10min，将上清液转移至新的1.5ml离心管中；
37.vi)加入1ml在-20℃冰箱预冷的无水乙醇，轻轻晃动至沉淀析出，在-20℃冰箱放置30min，12000r/min离心10min，弃上清；
38.vii)加入1ml的70％乙醇，轻轻晃动10min，12000r/min离心10min，弃上清；
39.viii)重复操作步骤ⅶ)，将离心管打开，放在超净工作台干燥至乙醇挥发完，加入200μl的灭菌水，在4℃冰箱过夜溶解dna。
40.ix)通过分光光度计检测dna的浓度和纯度，od260/280的值在1.8-2.0之间表明dna质量良好，将dna在-80℃冰箱保存。
41.b)1037个dna样本经质检合格后，随机选择165个dna样本，统一稀释至30ng/μl，每15个dna样本构建1个混池，共构建11个dna混池。
42.s3 pcr扩增
43.根据ncbi公布的绵羊slc27a6基因组序列(nc_056058.1)作为参考，利用primer-blast在线软件设计能够扩增slc27a6基因外显子1区域的pcr特异性引物，引物序列如下：
44.上游引物slc27a6-f：5
’‑
attgagttgtaagacagggcaca-3’45.下游引物slc27a6-r：5
’‑
tacgctgcttaggaggaatgag-3’46.反应体系(25μl)：12.5μl mix、10.5μl ddh2o、上游引物(10pmol/μl)和下游引物(10pmol/μl)各0.5μl，以及1μl模板(30ng/μl dna混池)。
47.反应程序：95℃预变性4min；94℃变性30s，63.8℃退火60s，72℃延伸90s，共34个循环；72℃再延伸9min。
48.将pcr产物进行电泳，结果如图1所示，可见除marker泳道，其右侧各dna混池泳道均扩增出目标片段(大小为302bp)。
49.s4 pcr产物测序与突变位点基因分型
50.将s3中由dna混池扩增出的pcr产物进行双向测序，根据如图2所示的测序峰图发现了一个位于绵羊slc27a6基因上的突变位点g.504g》c(参考序列为nc_056058.1)；
51.将质量检测合格后的1037个dna样本统一稀释至30ng/μl，分装至200μl的离心管中，并提供待定snp位点(即g.504g》c)在基因组序列上的位置信息，由公司用imldr技术对1037个dna样本进行snp分型。
52.s5 slc27a6基因分型位点的频率统计
53.等位基因频率计算公式如下：
54.pd＝(2ndd+ndd1+ndd2+ndd3+ndd4+
……
+nddn)/2n
55.式中，pd表示等位基因d频率，ndd表示群体中具有dd基因型的个体数量，nddi表示群体中具有ddi基因型个体数量，d1、d2、
…
、dn分别为等位基因d的n个互不相同的复等位基因；
56.基因型频率计算公式如下：
57.pdd＝ndd/n
58.式中，pdd代表某一位点的dd基因型频率，ndd表示群体中具有dd基因型的个体数，n为检测群体的总个体数量；
59.1037个体的等位基因和基因型频率计算结果如表1所示：
60.表1.湖羊slc27a6基因突变位点g.504g》c频率统计结果
[0061][0062]
注：hwb为哈代温伯格平衡；he为期望杂合度；ho为观察杂合度；pic为多态信息含量。
[0063]
从表1可以看出，g为g.504g》c位点的优势等位基因，等位基因频率为90.16％。且结合遗传学指标可以确定该突变位点为snp位点。
[0064]
2、湖羊slc27a6基因snp位点与脂肪酸百分比的关联分析
[0065]
基因型数据：以上采用imldr技术获得的分型结果
[0066]
肉质性状数据：湖羊背最长肌脂肪酸百分比(脂肪酸百分比是指某种脂肪酸占总脂肪酸含量的百分比)
[0067]
采用r-4.0.5软件分析基因位点与脂肪酸百分比的相关性，分析模型如下：
[0068]
y＝μ+g+t+e
[0069]
式中：y是脂肪酸百分比；μ是总体平均值；g是对应于snp的固定效应；t是批次效应；e是随机误差。
[0070]
通过分析数据确定slc27a6基因的snp位点对脂肪酸百分比的影响，寻找影响脂肪
酸百分比的潜在分子标记，结果如表2所示：
[0071]
表2.湖羊slc27a6基因g.504g》c与背最长肌脂肪酸百分比(mean
±
sd)的关联分析
[0072]
[0073]
注：同一行中不同的大写字母表示差异显著(p《0.05)，同一行中不同的小写字母表示差异极显著(p《0.01)。
[0074]
结果表明，g.504g》c位点对c6:0、c17:0、c18:1c6、c18:3n3、n-6/n-3pufa具有极显著影响，对c16:1、c22:0、c20:3n6、c20:3n3、c22:6n3、mufa具有显著影响。其中gg、cg、cc基因型个体的c6:0、c17:0、c18:1c6百分比均具有极显著差异，且在gg、cg、cc基因型中呈极显著降低趋势；gg、cg、cc基因型个体的c22:6n3百分比具有显著差异，且在gg、cg、cc基因型中呈显著降低趋势；cg基因型个体的c16:1、c22:0百分比显著高于gg基因型；cc基因型个体的mufa百分比显著高于gg、cg基因型个体；cc、cg、gg基因型个体的c18:3n3均具有极显著差异，且在cc、cg、gg基因型个体中呈显著降低趋势；cc基因型个体的c20:3n6、c20:3n3百分比显著低于cg和gg基因型，且在cc、cg、gg基因型个体中呈增加趋势；cc和cg基因型个体的n-6/n-3pufa极显著高于gg基因型。
[0075]
由于gg基因型个体中的pufa(具体为c20:3n6、c20:3n3、c22:6n3)百分比最高，gg基因型个体中的sfa(具体为c22:0)百分比最低，而pufa是有益脂肪酸，sfa可致使羊肉产生膻味，并能增加人体患心血管疾病的风险，因此可以通过选择建立基因型为gg(作为分子标记)的群体，对湖羊肉质性状(脂肪酸组成及含量)进行早期筛选，加快选育优异肉品质的绵羊新品系。