豬OLR1基因多態性與肉質性狀的關聯分析

郭咪咪 梅書棋 喬木 吳俊靜 趙康 程堯 李助南

摘要：為鑒定與豬的肉質性狀有重要影響的分子標記，經克隆測序，鑒定出OLR1基因第5內含子存在C52T突變，通過PCR-PstⅠ-RFLP方法進行了基因型檢測，分別計算了該SNP在大白豬、長白豬、梅山豬、大白×梅山F2代豬中的基因型頻率和基因頻率，并且對大白×梅山F2代豬的基因型與肉質相關性狀進行了關聯分析。結果表明，CC基因型的肩部膘厚、胸腰椎間膘厚、平均背膘極顯著高于CT基因型個體（P<0.01）；CC基因型的皮率、瘦肉率、6-7腰椎間膘厚、花油重、臀部膘厚、背最長肌大理石紋、股二頭肌大理石紋顯著高于CT基因型（P<0.05），C等位基因對豬的脂肪沉積具有正效應。

關鍵詞：豬；OLR1；PCR-PstⅠ-RFLP；肉質性狀；關聯分析

中圖分類號：S828 文獻標識碼：A 文章編號：0439-8114（2015）11-2686-04

DOI：10.14088/j.cnki.issn0439-8114.2015.11.034

Association Analysis on Meat Quality Traits and Polymorphism in Porcine OLR1 Gene

GUO Mi-mi1，MEI Shu-qi2，QIAO-Mu2，WU Jun-jing2，ZHAO Kang2，CHENG Yao1，LI Zhu-nan1

（1. College of Animal Science， Yangtze University， Jingzhou 434025， Hubei， China；

2.Institute of Animal Science and Veterinary Medicine， Hubei Academy of Agricultural Science， Wuhan 430064， China）

Abstract： To identify DNA markers which have significant effect on pig meat quality traits， the single nucleotide polymorphism （SNP） and genotype were tested by sequencing and PCR-PstⅠ-RFLP method，respectively. The results showed that a SNP （C52T） was found in intron 5，then the gene and genotype frequency in the Large White， Landrace， Meishan，Large White×Meishan F2s populations and the association analysis was performed in Large White×Meishan F2s pigs were analyzed. The results showed that the fat thickness， thorax-waist fat thickness and average backfat thickness of CC genotype pigs were extremely significantly higher than that of CT genotype（P<0.01）； CC genotype pigs had significant higher skin percentage， lean meat percentage，6-7 rib fat thickness，dressing percentage， hip fat thickness， meat marbling of longissimus dorsi and meat marbling of biceps femoris than that of CT genotype（P<0.05）. The results also showed that C allele had positive effect on fat deposition of pigs.

Key words： pig； OLR1； PCR-PstⅠ-RFLP； meat quality traits； association analysis

OLR1基因是動脈粥樣硬化的候選基因，主要氧化低密度脂蛋白（OxLDL）內皮細胞，調節動脈粥樣硬化病變和參與一些細胞過程，影響動脈粥樣硬化的發病機制[1]。OLR1基因對血管內皮細胞的損傷、參與動脈粥樣硬化的形成和發展過程中發揮了重要作用[2，3]。

目前關于OLR1基因的研究主要集中在對人的心血管疾病和奶牛的產奶性狀方面[4-6]。Liu等[7]對3′-UTR C188T和G501C多態性進行研究，結果表明，G501C突變位點的SNP對腦梗塞有顯著影響，3′-UTR處的SNP C188T對腦梗塞無顯著影響。Khatib等[8]在OLR1基因上發現8個單核苷酸多態性，其中第4外顯子上有2個SNP；第4內含子有5個SNP；3′-UTR有1個SNP，3′-UTR的g.8232的SNP多態性可能控制奶牛的乳脂量和乳脂率，該突變位點也影響OLR1 mRNA的翻譯及穩定性。Schennink等[9]對OLR1基因g.8232 C/A突變位點與荷斯坦奶牛產奶性狀進行關聯分析，發現該位點對奶牛乳脂率有顯著影響（P<0.05）。李托等[10]對秦川牛OLR1基因g.10497 C/A多態性與肉質性狀進行關聯分析，研究發現該位點對秦川牛眼肌面積和眼肌深度有顯著性影響。豬的OLR1基因位于豬5號染色體上。目前，關于豬的OLR1基因的報道較少，國內外還未對該基因多態性與豬的生產性狀進行關聯分析。本試驗對4種不同品種豬的OLR1基因基因型進行檢測，并對大白×梅山F2代豬的胴體性狀和肉質性狀進行了關聯分析，尋找影響豬肉性狀的分子標記，為豬的優良育種提供理論依據，以推動豬肉品質改良的育種進程。

1 材料與方法

1.1 材料

大白豬、長白豬、梅山豬、大白×梅山F2代豬均來自湖北省農業科學院優良育種豬場。

1.2 試劑與儀器

主要試劑：組織細胞基因組DNA提取試劑盒購自艾德萊生物公司、2×Taq Master Mix購自Yeastern Biotech公司、DL 2000 DNA Marker 購自TaKaRa生物公司、PstⅠ限制性內切酶購自湖北晶茂生物技術有限公司。

主要儀器：超微量分光光度計Photometer（北京天翔飛域儀器設備有限公司）、DNA Engine PCR儀、DYY-6B電泳儀（北京市六一儀器廠）、LRH-250F生化培養箱（上海一恒科學儀器有限公司）、凝膠成像系統（美國Bio-rad伯樂公司）。

1.3 方法

1.3.1 DNA提取與豬群性狀測定 采集試驗豬的耳組織約0.5 g，剪碎后按照DNA提取試劑盒（艾德萊生物公司）說明提取基因組DNA，用超微量分光光度計測定其濃度及OD值，將提取后的DNA置于-20 ℃保存。

試驗豬群的性狀測定包括活體性狀測定，體重達100 kg左右時，集中屠宰測定胴體性狀、肉質性狀和其他性狀。測定方法根據《種豬測定原理與方法》進行[11]。

1.3.2 引物設計 根據GenBank數據庫中的豬DNA序列（GenBank登錄號為CU856657.3），篩查SNP位點，設計引物P1擴增第4、5內含子序列；設計引物P2，進行PCR-RFLP分型，引物均由北京奧科鼎盛生物科技有限公司合成，引物序列及PCR條件見表1。

1.3.3 PCR擴增測序，篩查SNP 用引物P1分別對5頭梅山豬和5頭大白豬的DNA進行PCR擴增，PCR反應體系如下：2.0 μL DNA混合液模板，上下游引物各1.5 μL （10 μmol/L），25.0 μL 2×Taq Master Mix，加去離子水至50 μL。

PCR反應程序：94 ℃預變性4 min；94 ℃變性40 s，57 ℃退火40 s，72 ℃延伸50 s，共35個循環；72 ℃延伸10 min，4 ℃保存。將PCR產物送北京奧科生物技術有限責任公司進行雙向序列測定，利用Chromas軟件對測序結果進行分析。

1.3.4 PCR-PstⅠ-RFLP分析 用引物P2，PCR擴增體系如下：1.0 μL DNA模板，上下游引物各1.0 μL （10 μmol/L），12.5 μL 2×Taq Master Mix，加去離子水至總體積25 μL。PCR反應程序：94 ℃預變性 4 min；94 ℃變性40 s，53 ℃退火40 s，72 ℃延伸 45 s，共35個循環；72 ℃延伸10 min，最后4 ℃保存。

酶切體系：取5.5 μL PCR產物，加入PstⅠ內切酶0.5 μL （10 U/μL），1 μL 10×Buffer（含10×BSA），加去離子水至總體積10 μL，在37 ℃溫育過夜，用瓊脂糖凝膠電泳檢測酶切產物，成像系統觀察并記錄酶切結果。

1.3.5 統計分析 OLR1基因及基因型頻率計算方法參考文獻[12]進行，設等位基因C、T的頻率分別為p和q，則：

p=（2n1+n2）/2N

q=（2n3+n2）/2N

式中，n1為CC型個體數，n2為CT型個體數，n3為TT型個體數，N為總個體數。

基因型頻率=某一基因型個體數/所測個體總數

應用模型對242頭大白×梅山F2代豬OLR1基因PCR-PstⅠ-RFLP不同基因型與豬生產性狀進行關聯分析，所分析的性狀主要為胴體性狀和肉質性狀。

根據Liu[13]所述方法建立單標記回歸統計模型，采用SAS統計軟件（SAS Institute Inc，Version 8.0）的GLM過程進行單標記方差的統計分析，同時采用REG程序對基因加性效應和顯性效應進行計算，并進行顯著性檢驗，所采用模型為：

Yijkl=？滋+Gi+Fj+Sk+Yl+bijklXijkl+eijkl

式中，Yijkl為性狀表型值，？滋為平均值，Gi為基因型效應；Sk、Yl、Fj為固定效應，分別為性別、年度、家系效應，bijkl為屠宰體重或屠宰日齡的回歸系數，胴體性狀以屠宰體重為協變量，肉質性狀以屠宰日齡為協變量，eijkl為殘差效應。

2 結果與分析

2.1 豬OLR1基因的SNP研究

將獲得的5頭梅山豬和5頭大白豬OLR1基因測序結果，用軟件Chromas進行序列分析，分析表明，位點第5內含子的C52T突變，測序圖中雙峰處即為C/T突變，測序圖譜如圖1所示。

2.2 豬OLR1基因在不同群體中PCR-PstⅠ-RFLP基因型分析

分析發現，該SNP位點會被PstⅠ（CTGCA^G）限制性內切酶識別，針對該處突變位點設計了表 1所示的引物 P2進行PCR-PstⅠ-RFLP多態基因型分析。PCR產物長度為583 bp，當第5內含子的 C52T位點為C時，不能被PstⅠ酶切，酶切后只有一個片段，長度為583 bp（C等位基因）；當該位點為T時，PstⅠ酶切后產生兩個片段，長度分別為172 bp和 411 bp（T等位基因）。豬OLR1基因的PCR-PstⅠ-RFLP酶切分型結果如圖2所示。

利用PCR-PstⅠ-RFLP方法，對大白豬、長白豬、梅山豬和大白×梅山F2代豬群體進行基因分型。如表2所示，在長白豬群體中C等位基因的頻率等于T等位基因的頻率，大白、梅山和大白×梅山豬群體中C等位基因的頻率高于T等位基因的頻率，尤其是梅山豬群體中C等位基因的頻率高達1.00，等位基因C在群體中明顯占優勢。

2.3 豬OLR1基因的PCR-PstⅠ-RFLP基因型與生產性狀的統計分析

對豬OLR1基因PCR-PstⅠ-RFLP不同基因型與豬生產性狀進行了關聯分析，用于統計分析的試驗材料為242頭大白×梅山F2代豬，所分析的性狀主要為胴體性狀和肉質性狀。PCR-PstⅠ-RFLP基因型檢測結果表明，在所檢測的242頭試驗豬群體中，CC基因型159頭，CT基因型83頭，TT基因型0頭。不同基因型與生產性狀間的統計分析結果總結于表3。用SAS統計軟件的GLM和REG程序，對242頭大白×梅山F2代豬個體基因型與生產性狀進行了關聯分析，由表3可知，豬OLR1基因多態性與肩部膘厚、胸腰椎間膘厚、平均背膘極顯著相關（P<0.01）；與皮率、瘦肉率、6-7腰椎間膘厚、花油重、臀部膘厚、背最長肌大理石紋、股二頭肌大理石紋顯著相關（P<0.05）。基因效應方面花油重顯性效應達到極顯著水平（P<0.01）；6-7腰椎間膘厚、臀部膘厚的加性效應和內脂率、背最長肌肌肉色值、股二頭肌肌肉色值的顯性效應達到顯著水平（P<0.05）。

3 討論

豬的肉質性狀屬多基因控制的數量性狀，由一系列主基因或數量性狀位點（Quantitative trait locus，QTLs）決定[14]。分子標記輔助選擇在豬的選育方面起著重要的作用，促進了研究者利用相關的分子標記對豬的胴體及肉質性狀進行有效的選育，推動了豬肉質改良的育種進程。

對于該位點在不同品種豬群中的多態性分析，在長白豬群體中C等位基因頻率與T等位基因頻率相等，在大白豬、梅山豬和大白×梅山F2代豬群體中均是C等位基因高于T等位基因，尤其是脂肪型梅山豬群中C等位基因頻率高達1.00，該結果表明C等位基因為優勢等位基因，且C等位基因對豬的脂肪沉積具有正效應。在對4種不同豬群體的PCR-PstⅠ-RFLP分析中，均未檢測到TT基因型的存在，這可能與T等位基因在4種豬群體中含量少或者試驗的樣本量較少有關。

通過對豬OLR1基因PCR-PstⅠ-RFLP基因型與生產性狀關聯分析，發現CT基因型比CC基因型具有更高的瘦肉率，差異顯著（P<0.05）；胴體性狀的肩部膘厚、6-7腰椎間膘厚、胸腰椎間膘厚、花油重、臀部膘厚、平均背膘CC基因型比CT基因型高；肉質性狀中的背最長肌大理石紋、股二頭肌大理石紋CC基因型比CT基因型高，以上結果表明OLR1基因在豬的脂肪沉積方面具有積極的調控作用，其中CC基因型的個體具有更強的脂肪沉積能力。本研究發現該SNP位點是與豬脂肪沉積相關的分子標記，為標記輔助選擇奠定了基礎，由于豬脂肪生長發育是一個復雜的生理生化過程，受多種基因調控[15]，其具體的作用機理及遺傳機制有待進一步的研究。

參考文獻：

[1] PIRILLO A，CATAPANO A L. Soluble lectin-like oxidized low density lipoprotein receptor-1 as a biochemical marker for atherosclerosis-related diseases[J].Diseas Markers，2013，35（5）： 413-418.

[2] SAWAMURA T，KUME N，AOYAMA T，et al.An endothelial receptor for oxidized low-density lipoprotein[J].Nature，1997，386（6620）：73-77.

[3] CRUCET M，WUST S J，SPIELMANN P.Hypoxia enhances lipid uptake in macrophages：Role of the scavenger receptors Lox1，SRA and CD36[J].Atherosclerosis，2013，229（1）：110-117.

[4] 邊云飛，李茂蓮，楊慧宇，等.LOX-1特異性小發夾RNA表達載體的構建及其效應[J].中國藥物與臨床，2010，10（1）：42-44.

[5] WANG X，PENAGARICANO F，TALSTEIN R，et al.Association of an OLR1 polymorphism with milk production traits in the Israeli Holstein population[J]. Journal of Dairy Science，2012， 95（3）：1565-1567.

[6] JAWALE C V， SHABIR N， BHONG C D， et al. Characterization of SNPs in oxidised low density lipoprotein receptor 1 （OLR1） gene in Mehsana and Banni breeds of buffalo[J]. Biomedical Research，2014，25（3）：303-306.

[7] LIU X， ZHU R X，LI L，et al.Association of LOX-1 gene polymorphisms with cerebral infarction in northern Chinese Han population[J]. Lipids in Health and Disease，2014，13（55）：1-6.

[8] KHATIB H， LEONARD S D， SCHUTZKUS V. Association of the OLR1 gene with milk composition in holstein dairy cattle[J]. American Dairy Science Association，2006，89：1753-1760.

[9] SCHENNINK A，BOVENHUIS H， VISKER M H， et al. Effect of polymorphisms in the FASN，OLR1，PPARGC1A，PRL and STAT5A genes on bovine milk-fat composition[J]. Animal Genetics， 2009， 40（6）：909-916.

[10] 李 托，趙志東，王 婧，等.秦川牛OLR1基因多態性與肉質性狀的關聯分析[J].西北農林科技大學學報，2012，21（4）：17-20.

[11] 熊遠著，鄧昌彥.種豬測定原理與方法[M].北京：中國農業出版社，1999.

[12] 李 寧.動物遺傳學[M].北京：中國農業出版社，2003.

[13] LIU B H.Statistical genomics：linkage，mapping，and QTL analysis[M].Florida：CRC Press，1998.

[14] 李 鳳，黃業傳，吳照明.影響豬肉質的主效基因和候選基因[J].食品工業科技，2012，33（24）：442-444.

[15] 官紅平.SMPX基因3′UTR區中一個SNP的發現、檢測及與性狀的關聯分析[A].李輝.中國動物遺傳育種研究進展[M].北京：中國農業科學技術出版社，2005.