藏山羊FGF21基因克隆及生物信息學分析

沈閱　林亞秋　梅寒　李想　張小玉

摘要：為探討藏山羊FGF21基因的分子結構特征，利用RT-PCR技術對該基因進行克隆、測序及相關生物信息學分析。結果表明：獲得藏山羊FGF21基因序列650 bp，ORF長度為634 bp，編碼209個氨基酸。氨基酸序列分析顯示，藏山羊FGF21編碼蛋白屬于不穩定酸性蛋白，分子式為C₁₀₂₂H₁₅₉₆N₂₇₆O₂₉₆S₅，分子質量為22.65 ku。預測該蛋白含有13個磷酸化位點，并與多種蛋白存在相互作用；無跨膜結構域，有信號肽結構。在47～164位氨基酸殘基存在FGF結構域。亞細胞定位顯示，FGF21蛋白大部分位于細胞質內。同源性分析指出山羊FGF21氨基酸序列保守性較高，與GenBank中藏羚羊、牛、印度水牛及豬等的同源性在90%-99%。本研究結果為闡明FGF21基因在藏山羊中的作用提供重要的基礎數據。

關鍵詞：藏山羊；FGF21基因；克隆；生物信息學分析

中圖分類號：S827.2 文獻標志碼：A 文章編號：1002—1302（2016）01—0044—04

成纖維細胞生長因子（fibroblast growth factor，FGF）是一類由FGF基因家族編碼的結構相關的蛋白質。人類的FGF家族包括22個成員，按種系和序列又分為7個亞科，大部分FGF因子對細胞的生長和分化具有調節作用。成纖維細胞生長因子FGF21（fibroblast growth factor 21）是FGF家族的一個新成員，屬于FGF家族中FGF19亞家族成員之一。FGF21主要在肝臟、脂肪和其他組織表達，它對機體的能量平衡調節和糖代謝具有重要作用，已經成為全球糖脂代謝研究的熱點基因。

研究指出生酮飼喂小鼠可使小鼠體質量減輕，并且用生酮飼喂FGF21基因敲除大鼠，則老鼠出現體質量明顯上升現象。Coskun等報道FGF21過表達的轉基因鼠對多食誘導的肥胖有明顯的改善作用，并可逆轉肝臟脂肪變性及胰島素抵抗，使體質量下降約20%。學者在給患糖尿病的恒河猴注射FGF21后發現其不僅對血糖有控制作用，還具有輕度的減重作用。研究發現肥胖患者空腹血清FGF21表達水平顯著高于正常對照組，FGF21表達水平與腰圍以及WHR、FPG、HbA1c、FINS、TG、hsCRP的表達水平呈顯著正相關，與HDL-C呈顯著負相關。上述研究結果都反映出，FGF21在脂肪代謝和控制體質量方面都起了非常重要的作用，但上述研究均集中在人和小鼠等動物及細胞系上進行，在動物脂肪沉積方面的研究尚不深入。在山羊上尚未見相關報道，僅見山羊FGF21的預測序列，因此需要進行更深入的研究。

藏山羊是我國青藏高原特有畜種，又稱為“克什米爾”山羊，其肌肉細嫩、風味好，是藏族同胞重要的生產和生活資料。且藏山羊一般不需要專門肥育，是研究山羊優良肉質性狀形成的理想材料。因此本研究以藏山羊為研究對象，應用RT—PCR技術獲得FGF21基因序列，同時進行生物信息學分析，為進一步研究FGF21基因在藏山羊脂肪代謝及能量平衡中的作用機制，對提高藏山羊的育肥效率、提高肉品質提供新的理論依據。

1材料與方法

1.1材料

1.1.1試驗動物與組織的采集 以四川省若爾蓋種羊場飼養的成年健康公藏山羊為試驗動物，現場屠宰，迅速采集背最長肌組織樣本，立即置于液氮中保存帶回實驗室備用。

1.1.2主要試劑 Trizol試劑盒、熒光定量試劑盒SYBR@Premix Ex TaqTM（2×）、pMD-19T Vector均購自大連TaKa-Ra公司，反轉錄試劑盒Revert Aid First Strand cDNA SynthesisKit均購于Thermo公司，Taq聚合酶、膠回收試劑盒與E.ColiDH5a感受態細胞均購于天根生化科技有限公司。

1.2方法

1.2.1 RNA的提取及cDNA合成按照Trizol試劑盒說明書提取藏山羊背最長肌的總RNA，測定D_260nm與D_280nm比值在1.8～2.0之間，檢測RNA濃度和質量，然后取2μg按照反轉錄試劑盒說明書來合成cDNA第一鏈，-80℃保存備用。

1.2.2引物設計與合成根據GenBank上山羊FGF21基因預測序列（XM_005692688），利用Primer Premier 5.0軟件設計PCR克隆引物（正：5′-ACTGATGGGCTGGGACGA-3′；反：5′-TGCACAGCGGACGTCTTC-3′），產物長度為650 bp。引物送交四川成都擎科梓熙生物技術有限公司進行合成。

1.2.3藏山羊FGF21基因克隆 利用RT-PCR技術克隆臧山羊FGF21基因序列，PCR反應總體系為：2×Taq PCRMaster Mix 12.5μL，模板cDNA 1μL，上、下游引物各1μL（10μmol/L），最后加水至25μL。擴增條件：94℃4 min；94℃30 s，60℃30 s，72℃90 s，36個循環；72℃10 min。然后1%瓊脂糖凝膠電泳檢測擴增產物，膠回收試劑盒純化目的片段，然后再將產物與pMD-19T載體連接，并轉化至感受態細胞中，經菌液PCR鑒定后送至四川成都擎科梓熙生物技術有限公司測序。

1.2.4藏山羊FGF21基因和蛋白的生物信息學分析 Prot-Param（http：∥web.expasy.org/protparam/）分析蛋白理化性質；ProtParam（http：∥Web.expasy.org/protparam/）分析蛋白理化性質；NetPhos 2.0分析預測磷酸化位點；Signal P 4.1Server信號肽預測；TMHMM預測跨膜結構域；NPS預測二級結構，SWISS-MODEL三級結構預測；NCBI在線程序CDD預測保守功能域；通過PSORTⅡ進行亞細胞定位，利用STRING交互式數據庫進行搜索蛋白相互作用分析；NCBI中Blast在線比對同源性分析，Clustalx 1.83和MEGA 5.0構建進化樹；

2結果與分析

2.1山羊FGF21基因的克隆

用藏山羊背最長肌組織cDNA為模板經PCR擴增后獲得FGF21目的片段，經瓊脂糖凝膠電泳檢測與預期相符（圖1），測序后經Blast比對后確定其為藏山羊FGF21基因序列650 bp，分析表明FGF21基因序列起始密碼子ATG位于5 bp處，終止密碼子位于634 bp處，其中核苷酸序列的堿基組成分別為：A 17.2%，C 34.3%，G 29.8%，T 18.6%。

2.2藏山羊FGF21生物信息學分析

2.2.1蛋白理化性質分析 蛋白理化性質分析表明，藏山羊FGF21蛋白由209個氨基酸殘基組成（圖2），分子式為C₁₀₂₂H₁₅₉₆N₂₇₆O₂₉₆S₅，分子質量為22.65 ku。帶負電荷的氨基酸殘基總數（Asp+Glu）為22個，帶正電荷的氨基酸殘基總數（Arg+Lys）為20個，即該蛋白質可能帶負電荷。結合不穩定系數（52.g7）、理論等電點（6.08）和平均疏水指數（0.293）可知，FGF21蛋白為不穩定、酸性不溶類蛋白。且分析顯示藏山羊FGF21的47～164位氨基酸序列殘基存在FGF結構域（圖2、圖3）。

2.2.2蛋白磷酸化位點分析 經磷酸化位點預測，在76、79、122、152、195、200、204和206位共有8個絲氨酸（Ser）磷酸化位點，在68和98位共有2個蘇氨酸（Thr）磷酸化位點，在132、201和207位共有3個絡氨酸（Tyr）磷酸化位點。

2.2.3蛋白跨膜結構分析信號肽 分析結果指出山羊FGF21蛋白無跨膜螺旋結構。在1～28位氨基酸處存在信號肽。

2.2.4亞細胞定位及蛋白相互作用分析 通過PSORTⅡ預測，主要在細胞質（55.6%）和細胞核（22.2%）發揮生物學作用，少量作用于液泡（11.1%）和線粒體（11.1%）。利用STRING交互式數據庫進行搜索FGF21蛋白相互作用，條件限制在10個蛋白以內。搜索結果表明，FGF21蛋白可能和FGFR1、IGF1R、FGFR2、INSR、FGFR、PPARA、CSF1R、CYP7A1、FGFR3和FGFR4蛋白存在相互作用。構建蛋白相互作用網絡（圖4）。

2.2.5蛋白二級結構與三級結構預測 通過SOPMA在線軟件預測藏山羊FGF21蛋白二級結構（圖5），結果顯示，在藏山羊FGF21的209個氨基酸中，33個（15.7%）氨基酸可能形成α-螺旋（h），48個（22.97%）氨基酸可能形成β-折疊占（e），19個（9.09%）氨基酸可能形成β-轉角（t），109個（52.15%）氨基酸可能形成無規則卷曲（c）。通過I-TSSER軟件預測山羊FGF21蛋白三級結構（圖6），該三級結構與二級結構預測結果基本一致。

2.3藏山羊FGF21預測的氨基酸序列與其他物種同源性比較

通過NCBI中BlastN程序和BlastP程序分別進行同源性在線分析，結果表明，藏山羊FGF21推測的氨基酸序列與藏羚羊（XP_005955412）、牛（XP_005895443）、印度水牛（XP_006050969）、豬（NP_001156882）、馬（XP_001489202）、人（NP_061986）、小鼠（NP_064397）、大鼠（NP_570108）、家犬（XP_541510）和野駱駝（XP_006173890）的同源性分別為99%、99%、99%、91%、86%、85%、75%、74%、90%和90%，表明藏山羊與藏羚羊的親緣關系最近（圖7）。

3討論與結論

FGF21是近年來發現的一種與糖脂代謝關系密切的內源性脂肪細胞因子，有助于脂肪的代謝，是饑餓時抑制生長的重要因子。并且體外研究表明它可以降低血糖、抑制胰高血糖素的分泌、改善脂代謝。Varady等指出FGF21在豬的脂肪控制中具有調控作用，并且生長激素可以直接刺激牛肝臟中FGF21的表達。基于此，畜牧工作者對FGF21在動物脂代謝中的作用也非常關注。因此本試驗通過克隆得到藏山羊FGF21基因序列，為后續在山羊中的研究奠定基礎。

本試驗通過RT-PCR技術克隆得到藏山羊FGF21基因序列650 bp，編碼209個氨基酸，具有FGF家族典型的結構。推測的氨基酸序列與藏羚羊的同源性最近，這符合物種的進化規律。同時利用生物信息學分析了藏山羊的磷酸化位點、信號肽、跨膜結構、亞細胞定位、二級和三級結構以及蛋白相互作用等。結果發現，藏山羊FGF21沒有跨膜結構，具有信號肽結構，含有13個磷酸化位點，并主要發生在絲氨酸殘基上（8個），推測其翻譯后修飾主要是通過此來進行的。二級結構預測發現在藏山羊FGF21氨基酸序列主要有33個α-螺旋、48個β-折疊、19個β-轉角和109個無規則卷曲結構，并且三級結構預測也證明了二級結構的預測，有利于進一步闡明藏山羊FGF21的功能。

本研究發現經過蛋白質相互作用分析預測與藏山羊FGF21相互作用的蛋白可能是FGFR1、IGF1R、FGFR2、INSR、FGFR、PPARα、CSF1R、CYP7A1、FGFR3和FGFR4等。

有研究指出FGF21的啟動子上有PPARα結合位點，具有上調FGF21表達的功能，禁食和飼喂可以通過PPARα促進FGF21在肝臟中的大量表達，本試驗的預測結果與之一致。但關于FGF21在藏山羊中的具體的作用機制需要試驗來進一步證明。