靜寧雞PPARa基因克隆與生物信息學分析

楊妍梅　李玉　覃圣　蔣劉芽　趙金潭　陸曉東　熊杰　陳偉基　扎西英派

摘要：為闡明靜寧雞PPARa基因的結構及功能，根據NCBI上登錄的原雞PPARc基因的CDS序列設計引物，以靜寧雞腎臟為材料采用PCR的方法，對目的基因進行克隆，并測序。根據測序的結果，采用各類分析軟件，對所獲得的DNA片段進行生物信息學分析。結果成功克隆了靜寧雞PPARa基因的全長CDS序列。生物信息學分析結果表明，該基因CDS全長1407 bp，所編碼的蛋白質包含468個氨基酸。其分子量為52300，等電點為5.85。在二級和三級結構上，a-螺旋和無規則卷曲為蛋白的主要結構形式。PPARc蛋白是一個親水性蛋白質，不存在跨膜區域，不含信號肽，沒有N-糖基化位點，但存在16個0-糖基化位點，存在25個絲氨酸（Ser）磷酸化位點、12個蘇氨酸（Thr）磷酸化位點和6個酪氨酸（Tyr）磷酸化位點。靜寧雞PPARa蛋白在第100～168個氨基酸之間有1個ZnF_C4結構域，在第264～449個氨基酸之間有1個HOLI結構域。同源性分析結果表明，靜寧雞PPARa基因與原雞的親緣關系最近。靜寧雞PPARa基因的成功克隆和功能預測為進一步研究PPARa基因在靜寧雞脂肪代謝中的作用提供了基礎。

關鍵詞：靜寧雞;過氧化物酶體增殖劑激活受體;生物信息學分析

中圖分類號：Q785

文獻標識碼：A

文章編號：1000-4440（2019）02-0370-08

過氧化物酶體增殖劑激活受體（Peroxisome proliferators-activated receptors，PPARs）在1990年被英國科學家Isseman和Green首先發現"。它是一種新型的固醇類激素受體和配體依賴的核轉錄因子，屬于核受體第一亞家族C群（NR1C）[2]。PPARs是一個由3種核受體組成的家族，目前已知有3種亞型：PPARa、PPARγ、PPARβ。這3種亞型由不同的基因編碼，在組織中的表達和功能也不同[34。PPARs基因在心臟，脂肪組織，腦，腸，肌肉，脾臟，肺臟，腎上腺和大鼠脊髓中普遍存在。PPARs基因可被長鏈脂肪酸激活，廣泛存在于體內脂肪代謝旺盛的組織中。

近年來，國內外學者對畜禽PPARs進行了深入研究，如Sundvold等[5]對豬10種組織PPARs的Northern blot檢測發現，PPARa較高表達于肝臟和腎臟，中等程度表達于心肌、骨骼肌和小腸。Diot等[6]對雞9種組織進行Northern blot檢測，結果顯示雞PPARa mRNA在肝臟、心臟和腎臟中高表達，在尾脂腺也有較高的表達。孟和等[7]的研究結果表明，PPARa mRNA只在雞心臟、肝臟、腎臟和胃這4種組織中表達，在肝臟雜交信號最強。田亞東等[8]的研究結果表明，PPARa基因影響雞體脂肪代謝。孟和等[9]的研究結果表明，PPARa基因已在AA肉雞中檢測到單堿基突變位點，在3種中國地方品種中也檢測到，如石岐雜雞、北京油雞和白耳雞。該突變位點顯著影響雞群的腹脂質量和胴體性狀，推測該基因是影響雞體脂肪代謝的主效基因或與控制該性狀的主效基因連鎖，因此，它可作為雞脂肪性狀的標記輔助選擇的分子標記[9]，該研究結果與田亞東等[8]所得結果相似。

靜寧雞主產區分布在甘肅省靜寧縣和寧夏回族自治區固原市，甘肅省莊浪、通渭、華亭、秦安和會寧，寧夏回族自治區隆德、涇源和西吉等也有分布，靜寧雞以靜寧縣的雞源多和品質好而得名。靜寧雞以低脂肪、高蛋白、肉質鮮嫩、雞汁鮮美而聞名[10]。由于在實際生產中，靜寧雞生長速度較緩，在一定程度上影響了該品種的開發利用，以致靜寧雞日趨減少且表現出混雜退化的趨勢。隨著人們消費水平的提高，對低脂肪高蛋白質肉質的需求也逐漸提升。目前，PPARo的研究主要集中在嚙齒動物和人體上，并且在家禽中也有PPARc的報道，但是沒有關于靜寧雞PPARa的報道。因此，本研究以靜寧雞為研究對象，利用PCR獲得PPARa基因序列，然后利用生物信息學手段對該基因的結構及功能進行分析和預測，以期為進一步研究靜寧雞脂質代謝中PPARa的調控作用提供基礎。

1 材料與方法

1.1 試驗動物

6月齡靜寧雞購于靜寧綠洲生態農業科技發展有限公司。

1.2 試劑

Reverse Transcriptase M-MLV（RNase H-）等反轉錄相關試劑、TaKaRaExTaq等PCR反應相關試劑、pMD18-T Vector均購自寶生物工程（大連）有限公司，Trizol購自賽默飛世爾科技（中國）有限公司，氨芐青霉素購自北京索萊寶科技有限公司，Axyen膠回收試劑盒購自Axyen公司，BL2感受態細胞購于上海碧云天生物技術有限公司，質粒小量提取試劑盒購自生工生物工程（上海）股份有限公司。

1.3 方法

1.3.1 組織總RNA的提取和eDNA合成利用Trizol法提取靜寧雞腎臟的總RNA，以ddH2O溶解總RNA，-80°C保存備用;參考RTase M-MLV（RNaseH-）反轉錄操作步驟合成eDNA，-20°C保存備用。1.3.2引物的設計及PPARa基因的克隆根據NCBI公布的原雞的PPARa基因的CDS序列（NM_001001464），用Primer 5.0軟件進行引物設計。所設計的引物序列為：PPARa-F：5'-AGTGACCGCTC-TACTTGACCAA-3';PPARo-R：5'-CTCCGAAC-CGAGTGAACAGC-3 '，引物由寶生物工程（大連）有限公司合成。

以合成的eDNA第一鏈為模板，進行PCR擴增。反應體系為：TaKaRa Ex Taq （5 U/μl）0.5 μl，10xExTaq Buffer（mg2+ Plus）（20 mmol/L）2.5 μl，dNTPMixture（2.5 mmol/L）2.0 μl，引物PPARa-F（10μmol/L）0.5 μl，引物PPARa-R（10 μmol/L）0.5 μl，eDNA2.0μl，添加滅菌水至25.0μl。擴增反應程序為：95°C預變性5 min;95°C變性3 min，62.3°C退火45s，72°C延伸80s，35個循環，最后在72°C條件下進行10 min的延伸，4°C保存。擴增產物經過1.2%瓊脂糖凝膠電泳檢測后，切膠回收目的片段。與pMD18-T Vector載體16°C恒溫連接0.5 h后，利用熱激法轉BL2感受態細胞，挑出陽性克隆進行培養，然后提取質粒進行PCR驗證，并將質粒送至寶生物工程（大連）有限公司進行測序。測序結果用DNAStar、MegAlign和Editseq軟件比對和拼接獲得靜寧雞PPARa基因CDS序列。

1.3.3 PPARa基因的生物信息學分析采用在線軟件Protparam（http：//www.expasy.ch/tools/prot-param.html）預測PPAR蛋白的理化性質，采用Protscal（http：//www.expasy.ch/tools/pr-otscale.ht-ml）預測PPARa蛋白的疏水性，采用SignalP4.1（http：//www.cbs.dtu.dk/services/SignalP/）預測PPARa蛋白的信號肽，采用NetNGlye 1.0 Server（http：//www.cbs.dtu.dk/services/NetNGlyc/）預測PPARa蛋白的N-糖基化位點，采用NetOGlye 4.0Server（http：//www.cbs.dtu.dk/services/NetOGlyc/）預測PPARa蛋白的0-糖基化位點，采用NetPhos-3.1（http：//www.cbs.dtu.dk/services/NetPhos-3.1/）預測PPARo蛋白的磷酸化位點，采用TMpred（http：//ch.embnet.org/software/TMPRED-form.html）預測PPARa蛋白的跨膜結構，采用PSORT II Prediction（http：//psort.hge.jp/form.html）預測PPARo蛋白的亞細胞定位，采用APSSP2（https：//npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl？page=/NPSA/npsa_gor4.html）預測PPAR&蛋白二級結構，采用SWISS-MODEL（http：//swissmodel.Expasy.org/）預測PPARa蛋白三級結構，采用HMMER（https：//www.ebi.ac.uk/Tools/hmmer/search/jackhm-mer）預測PPARa蛋白保守結構，采用NCBI中BLAST進行PPARa蛋白同源性分析，最后，通過CLUSTALX對PPARo氨基酸序列進行比對，然后采用MEGA 5.10軟件基于NJ（Neighbor-joining）法構建PPARa蛋白的系統進化樹，距離校正使用泊松模型，計算時將空位完全刪除，自舉分析通過1000次循環實現。

2 結果與分析

2.1 靜寧雞PPARa基因的克隆

采用TRIzol法對靜寧雞腎臟的總RNA進行提取，并進行瓊脂糖凝膠電泳檢測。結果表明，總RNA條帶清晰，無拖尾，整體質量較高，符合基因克隆的要求。目的基因的PCR擴增產物進行瓊脂糖凝膠電泳檢測，結果顯示，在1500 bp左右有明亮的條帶，如圖1所示。分析測序結果表明PPARa基因CDS區的24 bp處是起始密碼子ATG，1430 bp處是終止密碼子TAA。該序列BLAST的結果顯示目的片段與NCBI上登錄的原雞PPARa基因的CDS序列（NM_001001464）99%相同，說明靜寧雞PPARa基因已被成功克隆出來。

2.2 靜寧雞PPARa蛋白的生物信息學分析

2.2.1 PPARx蛋白的氨基酸序列及理化性質分析通過在線軟件Protparam分析PPARo蛋白的理化性質，結果表明，PPARo的CDS區包含1407個核苷酸，編碼468個氨基酸，氨基酸組成見圖2，分子量為52300，理論等電點為5.85，含強酸性氨基（Asp+Glu）61個，強堿性氨基酸（Arg+Lys）52個，PPARa蛋白的分子式為C2301H3661 N6250gS3z2，組成的原子總數為7317個。預測該蛋白質不穩定系數為44.20，表明該蛋白質為一個不穩定蛋白質。在哺乳動物細胞中的半衰期為30 h，在酵母中的半衰期大于20 h，在細菌中的半衰期大于10 h。

2.2.2 靜寧雞PPARa蛋白的疏水性和跨膜結構域分析利用在線軟件Protscal預測靜寧雞PPARa蛋白氨基酸序列的親水性/疏水性，結果如圖3所示，第381位賴氨酸（Lys）疏水性最強（+3.156），第150位的丙氨酸（Ala）親水性最強（-3.078），Z值為-0.225，表明PPARo蛋白為親水性蛋白質。使用TMpred跨膜結構預測服務器對PPARa蛋白的氨基酸序列進行預測，結果發現，PPARo蛋白不存在跨膜區域，結果如圖4所示。

2.2.3 靜寧雞PPARa蛋白的信號肽和糖基化位點分析采用在線軟件SignalP 4.1 Server進行信號肽預測，結果顯示，PPARa蛋白不包含信號肽。使用NetNGlyc1.0 Server和NetOGlyc 4.0 Server分別預測PPARa蛋白N-糖基化位點和0-糖基化位點，結果發現PPARo蛋白沒有N-糖基化位點，但是在56、66、69、71、73、75、77、79、80、85、89、93、95、167、179、233處發現了O-糖基化位點。

2.2.4 靜寧雞PPARa蛋白磷酸化位點分析利用在線軟件NetPhos-3.1 Server預測PPARa蛋白磷酸化位點，結果顯示，有25個絲氨酸（Ser）磷酸化位點，分別在21、24、38、40、45、50、56、63、66、73、76、77、79、80、89、93、95、110、163、179、198、293、323、346和452位;有12個蘇氨酸（Thr）磷酸化位點，分別在12、71、129、246、253、279、283、285、288、307、438和450位;有6個酪氨酸（Tyr）磷酸化位點，分別在9、194、214、314.464和468位，結果如圖5所示。

2.2.5 靜寧雞PPARa蛋白的亞細胞定位預測利用PsortlIProtein Sorting Prediction軟件對靜寧雞PPARa蛋白進行亞細胞定位預測，結果顯示，靜寧雞PPARa蛋白主要定位于細胞質（47.8%），其他定位于細胞核（26.1%）、線粒體（17.4%）、內質網（4.3%）和高爾基體（4.3%）。由此預測靜寧雞PPARa蛋白可能為胞質蛋白。

2.2.6 靜寧雞PPARa蛋白保守結構域預測利用在線軟件HMMER預測，結果發現靜寧雞PPARa蛋白在第100～168個氨基酸之間有1個核激素受體中的C4鋅指結構（ZnF_C4結構域），在第264～449個氨基酸之間有1個核激素受體配體結合結構域（HOLI結構域）。

2.2.7 靜寧雞PPARa蛋白二級結構預測使用在線軟件NPSA對靜寧雞PPARa蛋白二級結構進行分析，結果如圖6所示，176個氨基酸參與ax-螺旋，占比為37.61%;83個氨基酸參與β_折疊，占比為17.74%;209個氨基酸參與無規則卷曲，占比為44.66%。由此說明，在靜寧雞PPARa蛋白二級結構中，a-螺旋和無規則卷曲是主要的結構形式。

2.2.8 靜寧雞PPARa蛋白的三級結構預測利用在線軟件SWISS-MODEL的自動建模功能，對靜寧雞PPARx蛋白三級結構進行預測，三維模型如圖7所示。用于建立該模型的氨基酸殘基范圍為101～419位，該模型以3e00.1.B（Peroxisome proliferator-activated receptor gamma）蛋白為模板，序列同源性為64.39%。應用NCBI網站中Conserved domainsearch工具分析，證明PPARa蛋白是過氧化物酶體增殖劑激活受體。

2.2.9 靜寧雞PPARa多序列比對分析及系統進化樹對靜寧雞PPARa核苷酸序列和其他從NCBI下載的PPARa 核苷酸序列進行多序列比對。經BLAST在線分析，與靜寧雞PPARa核苷酸序列同源性高達99%的為原雞（Gallus gallus，NM001001464.1），同源性高達98%的分別為火雞（Meleagris gallo-pavo，XM_010716016.2）和吐綬雞（Numida melea-gris，XM_0213932-56.1），同源性高達97%的為鵪鶉（Coturnix japonica，XM_015870635.1），同源性高達96%的分別為綠頭鴨（Anasplatyrhynchos，XM_021273687.1）、阿德利企鵝（Pygoscelis adeliae，XM009323646.1）、卷羽鵜鶘（Pelecanus crispus，XM.009484690.1）、帝企鵝（Aptenodytes forsteri，XM_009283535.2）、倉鶚（Tyto alba，XM_009965270.1）、鴻雁（Anser cygnoides，KJ010765.1）、灰雁（Anseranser，AF481797.1）、白尾鷚（Phaethon lepturus，XM_010290610.1）、白鷺（Egretta gaettar，XM_009648790.1）、灰冠鶴（Balearica regulorum gibbericeps，XM_010300944.1）、朱鸚（Nipponia Nippon，XM_009465045.1）和加拿大金雕（Aquila chrysaetos Cana-densis，XM_011583046.1）。靜寧雞與原雞相比，PPARa基因編碼的氨基酸在98（亮氨酸/谷氨酸）和186（脯氨酸/亮氨酸）處發生突變，結果如圖8所示。利用MEGA6.10軟件構建NJ進化樹，如圖9所示。分析結果表明靜寧雞與原雞，火雞，吐綬雞和鵪鶉的PPARa基因親緣關系較近，其中與原雞的親緣關系最近。

3 討論

PPARx是脂肪酸氧化酶基因的主要轉錄調控子，在調節脂質代謝平衡中起主要作用”。在本研究中，基于NCBI公布的原雞的PPARa基因的CDS序列，通過Primer 5.0軟件設計特異性引物，從靜寧雞腎臟中成功克隆了該基因的全長CDS序列。該CDS序列全長1407 bp，所編碼的蛋白質包含468個氨基酸。其分子量為52300，理論等電點為5.85，分子式為C2301H3661N625O698S32，組成的原子總數為7317個。馬云等研究發現鴨PPARa基因的cDNA全長1430 bp[13]，最長開放閱讀框為1407 bp，共編碼468個氨基酸[12]孟和等研究發現鵝PPARa基因的cDNA序列長度為1407 bp[13]林森等從藏雞肝臟組織中成功克隆得到PPARa基因，該基因編碼區大小為1404 bp，編碼468個氨基酸[14]。因此本研究結果與查閱的文獻結果一致，說明已經成功克隆出靜寧雞PPARa基因。

磷酸化是蛋白質最常見和最重要的一種蛋白翻譯后的修飾[15]，并在許多生物反應中通過磷酸基團的變化，即添加或去除，發揮著“開關”作用[16]。糖基化具有調節蛋白質功能的作用，例如改善凝膠性和水合性等[17]。本研究中，靜寧雞PPARa蛋白為親水性蛋白質，且不存在跨膜區域，不包含信號肽，沒有N-糖基化位點，但存在16個0-糖基化位點，存在25個絲氨酸（Ser）磷酸化位點、12個蘇氨酸（Thr）磷酸化位點和6個酪氨酸（Tyr）磷酸化位點，靜寧雞PPARa蛋白主要位于細胞質中。據報道，PPARa通過鋅指結構域與視黃酸受體形成異源二聚體，然后與配體結合區（HOLI）連接。配體激活后的異源二聚體與靶基因啟動子上的PPARa反應元件結合以激活靶基因活化，調節轉錄表達，然后調節脂類的氧化[18-19]。本研究中，利用在線分析軟件預測結果顯示靜寧雞PPARx蛋白在第100～168個氨基酸之間有1個ZnF_C4結構域，在第264～449個氨基酸之間有1個HOLI結構域，與查閱的文獻結果[4]一致。

在靜寧雞PPARo蛋白的二級和三級結構中，a-螺旋和無規則卷曲是其主要的結構形式，β-折疊占比相對較少。多系列比對分析及系統進化樹結果表明，靜寧雞PPARa基因與原雞，火雞，吐綬雞和鵪鶉的PPARa基因親緣關系較近，其中與原雞的親緣關系最近。

家禽PPARo基因調節脂肪代謝，影響家禽脂肪的沉積，從而影響腹脂率和胴體品質[2]，因此，靜寧雞PPARo基因的成功克隆和功能預測為進一步研究PPARo基因在靜寧雞脂代謝中的作用和揭示PPARo基因功能提供了基礎，為低脂高質量雞育種提供基礎研究資料。

參考文獻：

[1]ISSEMANN L，GREEN S.Activation of a member of the steroid hormone receptor superfamily by perxisome proliferators[J].Nature，1990，347（6294）：645-650.

[2]黃春紅，肖調義，劉巧林，等.畜產動物過氧化物酶體增殖物激活受體基因的表達模式及功能研究進展[J].動物營養學報，2013，25（5）：923-931.

[3]BRAISSANT O，FOUFELLE F，SCOTTO C，et al.Differential expression of peroxisome proliferator activated receptors（PPARs）：Tissue distribution ofPPAR-a，β，and-y in the adult rat[J].Endocrinol，1996，137：345-349.

[4]ELBRECHT A，CHEN Y，CULLINAN C A，et al.Molecular cloning，expression and characterization of human peroxisome proliferator activated receptors γl and γ2[J].Biochemical and Biophysical Research Communicatons，1996，224：431-437.

[5]SUNDVOLD H，GRINDFLEK E，LIEN S.Tissue distribution of porcine peroxisome proliferator-activated receptor 0：detection of an alternatively spliced mRNA[J].Gene，2001，273：105-113.

[6]DIOT C，DUAIRE M.Characterization of a cDNA sequence encoding the peroxisome proliferator-activated receptor a in the chicken[J].Poultry Science，1999，78：119-1202.

[7]孟和，李輝，王宇祥.雞PPARs基因組織表達特性的研究[J].遺傳學報，2004，31（7）：682-687.

[8]田亞東，亢娟娟，孫桂榮，等.PPARa基因對安卡x固始雞資源群胴體品質的遺傳效應分析[J].華北農學報，2010，25（6）：14-18.

[9]孟和，王桂華，王啟貴，等.雞PPARa基因單核苷酸多態與脂肪性狀相關的研究[J].遺傳學報，2002，29（2）：119-123.

[10]李宏強.靜原雞的地方品種形成及其種質資源保護措施[J].中國家禽，2009，31（3）：52.

[11]張玥，姜寧，蘇麗，等.PPARa與運動改善脂質代謝的關系[J].中國康復醫學雜志，2008，23（6）：495-504.

[12]馬云，王云云，張曉婷，等.鴨PPARa基因結構及功能的生物信息學分析[J].浙江大學學報（農業與生命科學版），2011，37（4）：371-379.

[13]孟和，李輝，王宇祥.鵝PPARa基因全長eDNA的克隆和序列分析[J].遺傳，2004，26（4）：469-472.

[14]林森，徐亞歐，林亞秋，等.藏雞PPARax基因克隆與生物信息學分析[J].西南民族大學學報（自然科學版），2015，41（6）：661-666.

[15]KIMJH，LEEJ，OH B，et al.Prediction of phosphorylation sites using SVMs[J].Bioinformatics，2004，20（17）：3179-3184.

[16]EMANUELSSON O，NIELSEN H，BRUNAK S，et al.Predicting subcellular localization of proteins based on their N-terminal amino acid sequence[J].Journal of Molecular Biology，2000，300（4）：1005-l016.

[17]MARCHAL I，GOLFIER G，DUGAS 0，et al.Bioinformatics in glycobiology[J].Biochimie，2003，85（1/2）：75-81.

[18]NEVE B P，FRUCHART J C，STAELS B.Role of the peroxisome proliferator-activated receptors（PPAR）in atherosclerosis[J].Biochemical Pharmacology，2000，60（8）：1245-1250.

[19]LALLOYER F，STAELS B.Fibrates，glitazones and peroxisome proliferator-activated receptors[J].Arteriosclerosis Thrombosis andVascular Biology，2010，30（5）：894-899.