牛SCD1基因的生物信息學分析

周秀敏，李強子，畢英杰，任衛(wèi)合，張 麗

（西北民族大學生命科學與工程學院，甘肅 蘭州 730030）

牛SCD1基因的生物信息學分析

周秀敏，李強子，畢英杰，任衛(wèi)合，張 麗?

（西北民族大學生命科學與工程學院，甘肅 蘭州 730030）

本研究利用生物信息學方法分析和預測牛SCD1基因的理化特性和編碼產物的結構。結果表明牛SCD1基因編碼359個氨基酸，平均分子質量為41.71 kDa；組成蛋白為不穩(wěn)定的水溶性蛋白，三級結構由μ-螺旋和無規(guī)則卷曲組成，并且牛SCD1基因氨基酸序列和山羊親緣關系最近。牛SCD1基因的生物信息學分析為進一步研究牛SCD1基因的遺傳特性和生理機制奠定了基礎。

牛；SCD1基因；生物信息學

硬脂酰輔酶A去飽和酶（Stearoyl-CoA desaturase-1，SCD）是一種合成單不飽和脂肪酸的限速酶，參與細胞的新陳代謝和前體脂肪細胞的分化。SCD基因的表達及其產物的活性決定了脂肪細胞單不飽和脂肪酸的合成及細胞膜憐脂和甘油三酯的組成，不僅與人類相關疾病如肥胖癥、糖尿病及高血壓等密切相關，也是改善牛奶和牛肉品質的主要候選基因[1]。牛SCD1基因位于26號染色體，由6個外顯子和5個內含子組成，在脂肪組織和乳腺組織都有表達[2]。王小龍等[3-4]表明SCD1基因與牛奶中飽和脂肪酸和單不飽和脂肪酸的比例及總不飽和指數(shù)有關，并且?guī)Щ驅εＨ橹械娜橹屎腿榈鞍茁视袠O大影響。

為深入探討牛SCD1基因的性質和功能，避免人工選育和該基因演化造成的差異，本研究利用生物信息學方法從分子水平對牛SCD1基因序列結構及相關生物信息學特征進行預測和分析，同時構建系統(tǒng)發(fā)育進化樹，為揭示該基因理化性質、功能信息、遺傳特性及相關生理機制提供理論基礎。

1 序列來源與方法

1.1 牛SSCCDD 11基因序列來源 從Gene Bank中查找得牛的SCD1基因的全長編碼區(qū)序列（Coding sequenc-es,CDS）（Genebank no.AY241932）。

1.2 方法 對牛SCD1基因CDS編碼產物的氨基酸序列、其蛋白質特性及功能結構域進行分析與預測所用各種工具軟件與網(wǎng)址見表1。

表1 牛SCD1基因進行生物信息學的主要軟件Table 1 The main software ofSCD1gene bioinformatics in cattle

2 結果

2.1 牛SSCCDD 11基因編碼蛋白質理化特性預測與分析

蛋白質的基本性質包括其相對分子質量、氨基酸組成和等電點等，其預測結果見表2。牛SCD1基因編碼359個氨基酸，20種氨基酸所占比例如圖1，亮氨酸（Leu）數(shù)目最多，為41個，占整個氨基酸組成的11.4%，半胱氨酸（Cys）目最少，為3個，占0.8%。基因編碼產物不穩(wěn)定指數(shù)大于40，表明該基因編碼產物不穩(wěn)定[5]。

表2 牛SCD1基因編碼蛋白質理化特性Table2 The encoding a protein physicochemical properties ofSCD1gene in cattle

2.2 牛SSCCDD 11基因編碼蛋白質疏水性//親水性預測和分析 牛SCD1基因編碼氨基酸序列親水性/疏水性預測結果如圖2所示，其中第64位蛋氨酸（Met）親水性最強（-2.867）；第225位的丙氨酸（Ala）疏水性最強（+3.200）。整個編碼產物中，親水氨基酸占63.23%，疏水氨基酸占 36.74%，總平均數(shù)為-0.223，表現(xiàn)為親水性，由此可推斷編牛SCD1基因編碼的蛋白質是一種可溶性蛋白。

2.3 牛SSCCDD 11基因編碼蛋白質三級結構預測與分析本研究采用SWISS-MODEL同源建模方法對牛SCD1基因編碼產物進行分析，并用PyMOL軟件進行編輯（圖3），結果顯示，牛SCD1主要由α-螺旋、無規(guī)則卷曲-組成。用Came中的ProSa對SCD1氨基酸序列3D模型中總能量進行評估[6-7]，確定此模型的立體構象合理，符合立體化學二面角分布和能量的要求。

2.4 牛SSCCDD 11基因編碼蛋白系統(tǒng)發(fā)育分析 利用Mega 6.0軟件對牛、山羊、小家鼠、大白鼠、綠頭鴨、大黃魚及三角帆蚌SCD1基因編碼的氨基酸序列構建系統(tǒng)發(fā)育樹（圖4）。SCD1基因幾乎在所有的動物中都有表達，從進化樹可看出同為偶蹄目牛和山羊聚為一類，親緣關系最近，與為軟體動物的三角帆蚌親緣關系最遠，嚙齒目中鼠類聚為一類，禽類、魚類以及軟體動物各自聚為一類，這與動物分類學基本一致。

3 討論

硬脂酰輔酶A去飽和酶I是單不飽和脂肪酸生成過程中的關鍵酶酶，單不飽和脂肪酸能影響牛肉大理石花紋等級、深加工性能、脂肪柔軟度以及牛肉的風味等質量[8]。因此對牛SCD1基因的深入研究能對牛的肉質性狀提供重要理論基礎。

圖1 牛SCD1基因編碼蛋白氨基酸組成Fig.1 The encoding a protein amino acid composition ofSCD1gene in cattle

圖2 牛SCD1基因編碼蛋白質的疏水性預測結果Fig.2 The encoding a protein hydrophobic predictions ofSCD1gene in cattle

基因決定生物的性狀，但最終是由蛋白質來體現(xiàn)的。目前對SCD基因的生物信息學分析僅見于羊[9-10]、豬[11]、烏賊[12]和鴨[13]。本研究通過SCD1蛋白質的理化性質、疏水性/親水性、蛋白質三級結構預測牛SCD1蛋白功能。從預測結果來看，牛SCD1蛋白是由359個氨基酸折疊而成的不穩(wěn)定親水蛋白，其氨基酸數(shù)目與羊[10]和豬[11]相同，但比烏賊[12]的稍長，而豬[11]和烏賊[12]的SCD蛋白是疏水蛋白這可能是物種差異造成的。有研究表明蛋白質半衰期越長則蛋白質穩(wěn)定性越高[14]，而牛SCD1編碼蛋白質半衰期為30 h，但其不穩(wěn)定系數(shù)值大于40，與豬[11]的SCD蛋白相同，均是不穩(wěn)定蛋白，可能是不同基因發(fā)揮不同生理功能而出現(xiàn)與其的研究結果相悖。牛SCD1蛋白三級結構主要以α-螺旋和無規(guī)則卷曲為主，并未發(fā)現(xiàn)豬[11]SCD蛋白所報道的β-疊片。雖然豬和牛同為偶蹄目，但屬不同科，因此出現(xiàn)了該差異。從進化樹可看出同為偶蹄目牛和山羊聚為一類，親緣關系最近，與為軟體動物的三角帆蚌親緣關系最遠，嚙齒目中鼠類聚為一類，禽類、魚類以及軟體動物各自聚為一類，這與動物分類學基本一致。

圖3 牛SCD1基因編碼蛋白質三級級結構模型預測Fig.3 Putative result ofSCD1protein tertiary structure in cattle

圖4 牛SCD1基因編碼蛋白質的氨基酸系統(tǒng)進化樹（節(jié)點處數(shù)字表示置信度）Fig.4 Phylogenetic tree with amino acid sequences of SCD1 gene（bootstrap values are in dicated at nodes）

本研究系統(tǒng)地對牛SCD1基因的理化性質和蛋白的結構功能進行了預測，在此基礎上將進一步對SCD1基因的表達差異進行探究，以期能更好地為牛的肉質性狀提供理論基礎，更好地將分子生物學技術應用到生產實踐中去。

[1]盧建雄，劉翊中，臧榮鑫，等.鼠脂肪細胞分化過程中脂代謝相關基因轉錄表達時序的研究[J].中國獸醫(yī)科學，2006，36（4）：306-310.

[2]Myungil Chung,Seckho HA,Soonseog JEONG,et al. Cloning and Characterization of Bovine Stearoyl CoA Deasturasel cDNA from Adipose Tissues[J].Journal of Biosciences,2000,64 (7)∶1526-1530.

[3]王小龍，陳仁金，楊章平，等.中國荷斯坦牛SCD1基因多態(tài)性與泌乳性狀的關聯(lián)分析[J].畜牧獸醫(yī)學報，2011，42（7）：1022-1026.

[4]王小龍，常玲玲，陳瑩，等.中國荷斯坦牛SCD1外顯子5基因多態(tài)性及其對乳中脂肪酸組成的影響[J].中國農業(yè)科學，2014，47（9）：1858-1864.

[5]Guruprasad K,Reddy B V B,Pandit M W.Corre-lation between stability of a protein and its dipeptide composition∶a novel approach for predicting in vivo stability of a protein from its primary sequence[J]. Protein Engineering,1990,4(2)∶155-161.

[6]Wiederstein M,Sippl M J.ProSA-web∶interactive web service for the recognition of errors in three-dimensional structures of proteins[J]. Nucleic acids research, 2007, 35 (suppl 2)∶W407-W410.

[7]Sippl M J.Recognition of errors in three-dimensional structures of proteins[J]. Proteins∶Structure,Function,and Bioinformatics,1993,17(4)∶355-362.

[8]Smith S B,Yang A,Larsen T W,et al.Positional analysis of triacylglycerols from bovine adipose tissue lipids varying in degree of unsaturation[J]. Lipids, 1998, 33(2)∶197-207.

[9]石恒波，羅軍，朱越，等.奶山羊SCD基因CDS區(qū)的克隆、序列分析及過表達[J].中國農業(yè)科學，2012，45（24）：5091-5101.

[10]周恒，曹艷紅，朱江江，等.隆林山羊SCD基因的克隆及序列分析[J].黑龍江畜牧獸醫(yī)，2016（1）：80-83.

[11]曹威榮.豬硬脂酰輔酶A去飽和酶蛋白生物信息學分析[J].畜牧與飼料科學，2014，35（11）：5-7.

[12]夏靈敏，馬晟杰，吳常文，等.曼氏無針烏賊SCD蛋白結構和功能預測[J].浙江海洋學院學報（自然科學版），2014（2）：161-164.

[13]王健，張軍，張蕊，等.番鴨SCD1基因克隆、序列分析及組織表達研究[A].中國畜牧獸醫(yī)學會家禽學分會第九次代表會議暨第十六次全國家禽學術討論會論文集[C].2013.

[14]劉麗霞，張麗，趙生國，等.八眉豬DQA基因編碼區(qū)多態(tài)性及生物信息學分析[J].華北農學報，2015，30（1）：103-108.

Bioinformatics Analysis of CattleSCD1Gene

Zhou Xiumin,Li Qiangzi,Bi Yingjie,Ren Weihe,Zhang Li*
（College of Life Science and Engineering,Northwest University for Nationalities,Gansu Lanzhou 730030）

This study analyzed the physical and chemical characteristics and the prediction cattleSCD1gene structure and function of the encoded product between homologous bioinformatics methods.The results showed thatSCD1gene encoded 359 amino acids,and had a predicted molecular weight of 41.71 kDa,unstable protein composition soluble proteins,the tertiary structure of the alpha helix,and dandom coil.The amino acid sequence of cattle SCD1 gene was closely related to goat.The result lay a foundation for further studying on the genetic characteristics and physiological mechanisms of SCD1 gene.

Cattle;SCD1 Gene;Bioinformatic

S823

B

1672-9692(2016)12-0030-05

2016-10-15

周秀敏（1996-），女，江蘇鹽城人，在讀本科生，研究方向為動物科學。

張麗（1980-），女，寧夏石嘴山人，博士，主要從事動物遺傳育種研究。

西北民族大學引進人才科研項目（xbmuyjrc201316），西北民族大學生命科學與工程學院眾創(chuàng)空間扶持項目。