荷斯坦奶牛FGFR2基因組織表達及生物信息學分析

鄧惜緣，劉麗元，周子航，溫 萬，高旭紅，顧亞玲

（1. 寧夏大學農學院，銀川 750021；2.寧夏畜牧站，銀川 750004）

荷斯坦奶牛原產于荷蘭，為一種風土馴化能力極強的品種，適應性強、體型大、產奶量高，是我國存欄量最多品種，全年平均產乳量超7 000 kg。我國荷斯坦奶牛又分為乳用型荷斯坦奶牛和兼用型荷斯坦奶牛。

成纖維細胞生長因子受體家族（Fibroblast growth factor receptor，FGFR）包括5 個基因，其中FGFR1、FGFR2、FGFR3 和FGFR4，為跨膜酪氨酸激酶受體，其胞外區具有2 或3 個免疫球蛋白（Ig）樣結構域，而FGFR5 缺乏細胞內蛋白酪氨酸激酶結構域，該結構域被一個短的富含組氨酸基序的胞內尾部取代[1]。FGFR2 全稱為成纖維細胞生長因子受體2，位于奶牛第26 號染色體，全長107 412 bp，具有不同亞型，在上皮細胞中編碼FGFR2b 亞型，在間質細胞中編碼FGFR2c 亞型。FGFR2 基因分別結合不同FGF 激活信號通路，介導FGFs 信號傳遞入細胞質，為b FGF 最主要受體。b FGF與其結合后可激活細胞內絲裂原活化蛋白激酶、細胞外信號調節激酶（MAPK/ERK）和磷脂酰肌醇-3 激酶和蛋白激酶B（PI3K/AKT）等信號通路調節血管內皮細胞、上皮細胞、成纖維細胞、成骨細胞和神經膠質細胞等多種起源于中胚層、神經外胚層的細胞分化增殖，還參與機體多種生理病理學過程（如胚胎發育、組織愈合、腫瘤發生與轉移）[2]。Cunningham 等最新研究擴展FGF 信號已知領域，并為功能研究確定許多新靶點[3]。FGF信號通路結構靈活，在典型的RAF/MAPK/ERK/RSK和PI3K/AKT/PDK/mTOR/S6K 通 路 中，FGFR2 激 酶活性抑制反應被識別。同時觀察到磷酸化依賴性負反饋途徑被抑制，定義FGFR2 抑制的內在抗性機制。該發現對FGFR 抑制劑治療應用有指導意義，因FGFR抑制劑識別具有相同遺傳損傷和耐藥途徑細胞的反應。

本課題組以寧夏地區1 280頭荷斯坦母牛為研究對象，利用Illumina 牛150K Bovine 全基因組SNP 芯片作基因分型，對5 個泌乳性狀（產奶量、乳脂率、乳脂量、乳蛋白率、乳蛋白量）與體細胞評分，運用Farm CPU 統計方法作全基因組關聯分析[4]。在奶牛第26號染色體上存在P值4.68×10-8且名為ARS-BFGL-NGS-106897 的SNP 位點與奶牛體細胞數有關，與其相距30 kb即41 956 121 bp存在FGFR2基因且最小等位基因頻率為0.55。通過GO分析發現，FGFR2 參與FGF 信號通路（P00021），也參與多種生物學過程，如MAPK 級聯、細胞分化、細胞增殖和凋亡過程負調控。因此本研究分析前期GWAS 篩選出的FGFR2 基因不同組織表達量，再分析荷斯坦奶牛FGFR2 基因生物信息學，以期為進一步研究荷斯坦奶牛FGFR2 基因生物學功能提供參考，為后期FGFR2 基因功能驗證提供依據。

1 材料與方法

1.1 材料

1.1.1 試驗樣品

荷斯坦奶牛來自寧夏回族自治銀川市賀蘭縣茂盛奶牛場，選擇3頭具有相同父本且飼養管理條件相同的荷斯坦奶牛，取屠宰后奶牛組織（肝臟、卵巢、乳腺、腎臟、心臟、子宮），DEPC 水洗凈，液氮速凍后送回實驗室，-80 ℃保存，備用。

1.1.2 主要試劑與儀器

主要試劑：LA Taq（DRR02AG）(TaKaRa)，反轉錄試劑盒（TUREscript 1st Stand cDNA SYNTHESIS Kit）（艾德來）；2×SYBR?Green 預混液（中國DF）；pcDNA3.1-EGFP-C 載體（優寶生物）；DH5α感受態細胞（天根）。

主要試驗儀器：analytikjena-qTOWER 2.2型熒光定量PCR 儀（德國）、SCILOGEX D3024R 離心機（美國）、普通PCR 儀analytikjena-Easycycler（德國）、超微量核酸蛋白測定儀scandrop100（德國）、移液器（美國bio-rad）、U-3010 紫外-可見分光光度計（Hitachi）、凝膠成像系統（Gene Genius）、H6-1微型電泳槽（上海精益有機玻璃制品儀器廠）、SW-CJ-1D 潔凈工作臺（上海琪特分析儀器有限公司）、DYY-8 型穩壓穩流電泳（江蘇蘇潔凈化設備廠）、PCR 反應儀擴增儀（加拿大BBI）、YXJ-2 離心機（湘儀離心機儀器有限公司）。

1.2 方法

1.2.1 總RNA提取及反轉錄

使用Trizol 提取荷斯坦奶牛不同組織總RNA，通過凝膠電泳試驗檢測提取效果，超微量核酸蛋白測定儀（scandrop100）檢測RNA其OD260/OD280值。

采用Aidlab 公司反轉錄試劑盒（TUREscript 1st Stand cDNA SYNTHESIS Kit）反轉錄，采用20 μL 反應體系：總RNA 500 ng、5×RT Reaction Mix 4 μL、Rondam primer/oligodT 0.8 μL、 TUREscript HRTase/RI Mix 0.8 μL、RNase Free dH2O 加至20 μL。反轉錄反應條件如下：42 ℃40 min，65 ℃10 min反應結束后，得到cDNA，-80 ℃保存。

1.2.2 奶牛FGFR2基因克隆

根據GenBank 中已知牛FGFR2 基因序列（XM_024985484.1）及GADPH 序列（NM_0010 34034.2），利用Primer Premier 5.0 設計特異性熒光定量引物，由陜西致研生物科技有限公司合成。

由于CDS 區序列較長，分六段擴增，由上海生工生物工程股份有限公司合成。具體引物序列如表1所示。

PCR 擴增體系（25 μL）：2X GC Buffer I 12.5 μL，上下游引物各0.5 μL，dNTP（10 mmol·L-1）0.2 μL，ddH2O 10.1 μL，Template 1 μL（cDNA）、Taq 酶(5 U·μL-1) 0.2 μL。PCR 反應條件：95 ℃預變性3 min，94 ℃變性30 s，58 ℃退火30 s，72 ℃延伸2 min，72 ℃修復延伸7 min，循環33次。PCR產物用1%瓊脂糖凝膠電泳檢測，特異性好擴增產物保存于-20 ℃冰箱備用。純化后產物通過TA 克隆將該序列在16 ℃連接至pcDNA3.1-EGFP-C 載體，轉化到大腸桿菌DH5α 感受態細胞，涂板于LB 固體培養基，37 ℃培養12 h，藍白斑篩選，選擇白色菌落菌液送至上海生工生物工程有限公司測序。

1.2.3 組織表達差異分析

根據序列設計實時熒光定量PCR 引物, 檢測FGFR2 基因在心臟、肝臟、乳腺、子宮、卵巢和腎臟組織表達量。反應體系（10 μL）：2×SYBR?Green Supermix 5 μL，上下游引物各0.5 μL，cDNA 1 μL，ddH2O 3 μL。反應程序：95 ℃PCR 預變性3 min，95 ℃變性10 s，58 ℃組合延伸30 s，循環39 次。共設置3 個生物學重復以及3 個技術學重復，以GAPDH為內參，利用2-△△Ct計算各樣品基因相對表達量，利用SPASS 22.0 分析荷斯坦奶牛FGFR2基因不同組織表達量單因素差異。

1.2.4 生物信息學分析

基 本BLAST 檢 索（http：//www.ncbi.nlm.nih.gov/blast）；利用Prot-Param 在線軟件分析理化性質（http：//web. expasy.org/protparam /）；利用DNAMAN比對序列；運用SignalP 4.1 軟件分析信號肽（http：//www.ebs.dtu.dk/services/SignalP/）；運用ProtScale 軟件分析親疏水性（http：//web.expasy.org /protscale/）；通過在線軟件TMHMM（http：//www.cbs.dtu.dk/services/TMHMM/）預測蛋白質跨膜結構；保守結構域通過NCBI 上CDD 預測（http：//www.ncbi.nlm.nih.gov/cdd）；運用Psort 軟件分析亞細胞定位（https：//www.genscript.com/psort.html）；通過在線軟件Prabi 服務器SOPMA 程序（https：//npsa-prabi.ibcp.fr/）預測蛋白質二級結構，通過在線軟件SWISS-MODEL（https：//npsa-prabi.ibcp.fr/cgi-bin/npsa-automat.pl）預測蛋白質三級結構，研究確定FGFR2 蛋白質各項生物信息，了解該基因部分功能。

2 結果與分析

2.1 FGFR2基因RNA提取與擴增及克隆測序

不同組織RNA 提取后，用1.5%瓊脂糖凝膠電泳，電泳效果見圖1，產物用超微量核酸蛋白測定儀（scandrop100）檢測RNA OD 值（見表2），可看到明顯的28S與18S條帶，全部符合RNA要求。以荷斯坦奶牛乳腺組織為模板，分六段擴增，引物特異性較好，六段目的基因片段長度符合要求，結果見圖2。

2.2 荷斯坦奶牛FGFR2基因組織表達分析

目的基因FGFR2擴增曲線圖3，圖4為溶解曲線，產物TM 值為80 ℃。內參基因GAPDH 為擴增曲線圖5，圖6 為溶解曲線，產物TM 值為83 ℃。綜合溶解曲線與擴增曲線圖，試驗結果可靠。

表2 不同組織RNA提取濃度與純度檢驗Table 2 Test of RNA extraction concentration and purity in different tissues

以GAPDH為內參，利用2-△△Ct計算各樣品基因相對表達量，利用SPSS 25.0 一般線性模型（GLM）作平均值下單因素方差分析，Duncan's多重檢驗結果顯示，乳腺組織FGFR2 基因mRNA 表達量最高，并極顯著高于肝臟、卵巢、子宮和心臟，與腎臟組織mRNA表達量差異不顯著（見圖7）。

2.3 荷斯坦奶牛FGFR2蛋白理化性質分析

ProtParam 在線服務器分析結果如圖8 和表3所示，奶牛FGFR2 基因編碼790 個氨基酸，20 種氨基酸中亮氨酸（Leu）數目最多，共64 個，占全部氨基酸8.1%。負電荷殘基（Asp+Glu）總數為100，正電荷殘基（Arg + Lys）總數為89，分子式為C3923H6159N1067O1190S44，相對分子質量為12 383。半衰期為30 h，脂肪指數為78.71，親水性平均值為-0.409?；蚓幋a產物不穩定指數為46.09，顯示編碼產物不穩定。

表3 荷斯坦奶牛FGFR2蛋白氨基酸組成Table 3 Amino acid composition of FGFR2 protein in Holstein cows

2.4 荷斯坦奶牛FGFR2基因序列相似性及系統進化樹分析

利用Meg Align 軟件將克隆得到的荷斯坦奶牛FGFR2 基因序列與家貓、黃牛、家雞、獼猴、黑猩猩、野豬、家犬、綿羊、熱帶爪蟾、蘇門答臘猩猩基因序列作相似性比對，其中荷斯坦奶牛與黃牛序列相似性高達100%，與熱帶爪蟾相似性最低為67.2%（見表4）?；赥amura-Nei 模型最大似然法推導進化樹模型，選擇具有最高對數可行性的樹（-6 150.52），將Neighbor-Join 和BioNJ算法應用于最大復合似然（MCL）方法估計成對距離矩陣，選擇具有優越對數似然值拓撲，最后獲得用于啟發式搜索的初始樹。分析涉及10 個核苷酸序列，所有包含空白和缺失數據的位置均刪除，所有進化分析在MEGA7.0 中完成，結果顯示荷斯坦奶牛與黃牛親緣性最近，其次為野豬；而與熱帶爪蟾較低等脊椎動物相距最遠（見圖9）。

2.5 荷斯坦奶牛FGFR2蛋白疏水性分析

運用ExPASy 服務器中ProtScale 程序分析奶牛FGFR2 蛋白親疏水性發現，第358 位丙氨酸（Ala）疏水性最強為3.211，第167 位賴氨酸（Lys）疏水性最弱為-2.622，在負值區域有較多氨基酸（見圖10）。

表4 荷斯坦奶牛與其他10 個物種FGFR2基因序列相似性比對Table 4 Sequence similarity comparison of FGFR2 gene between Holstein cows and other 10 species

2.6 荷斯坦奶牛FGFR2蛋白跨膜結構預測及信號肽分析

運用TMHMM 2.0 分析奶牛FGFR2 基因蛋白跨膜結構，預測結果如圖11 所示。結果表明，FGFR2 蛋白存在1 個跨膜結構，位于第344～366 氨基酸殘基之間。使用SignaIP 4.1在線軟件預測蛋白信號肽發現，編碼產物C 值為0.112，Y 值為0.107，S 值為0.116，荷斯坦奶牛FGFR2 蛋白不具有信號肽（見圖12）。結合軟件SecretomeP2.0（http：//www.cbs.dtu.dk/services/SecretomeP/）分 析， 得 到SecP score分數高于0.6，則荷斯坦奶牛FGFR2基因編碼蛋白屬于非經典分泌蛋白。

2.7 荷斯坦奶牛FGFR2蛋白二級結構和三級結構預測

由圖13 可知，α-螺旋占23.80%，β-轉角占4.18%，無規則卷曲占51.27%， 延伸鏈占20.76%。由此推測，無規卷曲為荷斯坦奶牛FGFR2 最大量二級結構元件，α-螺旋、β-轉角和延伸鏈散布其中。由圖14可知，FGFR2蛋白以無規卷曲和α-螺旋為主，與二級結構分析一致。

2.8 荷斯坦奶牛FGFR2蛋白亞細胞定位分析

運用PSORTⅡ預測FGFR2 編碼產物亞細胞分布情況：細胞質21.7%，細胞核43.5%，線粒體17.4%，質膜4.3%，液泡4.3%，過氧化物酶體4.3%，分泌系統小泡4.3%。

2.9 荷斯坦奶牛FGFR2蛋白保守結構域預測

運用NCBI 的CDD 數據庫預測荷斯坦奶牛FGFR2 蛋白保守結構域。如圖15 所示，荷斯坦奶牛FGFR2蛋白有3個特定匹配和2個匹配的超家族保守結構域。

3 討 論

FGFR2 基因已確定與胃癌[5-6]、乳腺癌[7-8]、和膽管癌[9]等多種癌癥密切相關，且Pemazyre 治療FGFR2 基因融合型膽管癌已獲得美國FDA 首批認證[10-12]。在骨骼發育過程中, FGFR2 基因可調節位于軟骨膜、骨膜和顱縫處前骨生成細胞和骨生成細胞增生和分化[13]。Brajadenta 等研究發現FGFR2基因中至少有7 個突變引起Apert（尖頭并指綜合征）綜合征[14]。Kim 等研究發現FGFR2 信號為導致面中部前后向空間異常關鍵因素，研究結果為預防FGFR2基因突變引起的顱面畸形提供新療法[15]。為證明FGFR2 主要影響軟骨細胞增生和分化，尹良軍等應用基因敲入技術建立FGFR2 功能獲得型基因突變小鼠模型，研究發現FGFR2 不僅具有調控成骨細胞系發育功能，還具有調控軟骨細胞系發育重要功能，從而同時影響膜內成骨和軟骨內成骨，在骨骼發育、骨折修復等方面發揮重要作用[16]。可知FGFR2在調節骨相關生長發育中發揮作用。

分析克隆的荷斯坦奶牛FGFR2 基因序列發現，該基因CDS 區序列長2 372 bp，編碼790 個氨基酸，其中亮氨酸含量最多，亮氨酸增加蛋白質合成，抑制蛋白質分解。亮氨酸代謝產物α- 酮戊己酸（α-KIC）和β- 羥基-β- 甲基丁酸（HMB）具有調節蛋白質代謝作用。帶負電荷殘基總數（Asp+Glu）大于帶正電荷殘基總數（Arg+Lys），故FGFR2基因編碼蛋白質帶負電荷， 其不穩定指數為46.09，為不穩定蛋白，預測該蛋白為親水性蛋白。由于具有跨膜結構但無信號肽，則推測為非經典分泌蛋白，與FGF2分泌途徑：由質膜直接易位介導，依賴于質膜內磷脂酰肌醇（4,5）二磷酸和質膜外層的硫酸乙酰肝素蛋白多糖結論一致[17-18]。

畢艷楠等研究發現奶牛乳腺炎金黃色葡萄球菌誘導小鼠乳腺組織中FGFR2 表達量顯著升高[19]，表明FGFR2 基因與奶牛乳腺纖維化分子作用機制關系密切，與前期GWAS 篩選結果一致，且與不同組織表達差異分析得出乳腺組織表達量最高結果也一致。王澤鑫等通過瞬時轉染pcDNA3.1-TOPO-FGFR2 真核表達質粒進入人膽管癌Hucct1 和HEK293T 細胞，在細胞水平驗證FGFR2基因過表達具有促進細胞生長作用，與GWAS 結論相符[20]。同時，畢艷楠等發現，FGFR2與bFGF共同作用促進奶牛乳腺成纖維細胞中PI3K 和Akt mRNA和蛋白表達[21]。

基于前期GWAS 及不同組織差異分析得出乳腺組織相對表達量高的結果，再結合其他學者研究可推斷荷斯坦奶牛FGFR2 基因具有促進乳腺細胞生長作用，后期將開展細胞水平功能驗證。

4 結 論

本研究成功克隆荷斯坦奶牛FGFR2基因, 得到2 373 bp 的CDS 區序列，編碼790 個氨基酸殘基；生物信息學分析結果顯示，FGFR2 蛋白為非經典分泌蛋白，編碼產物主要分布在細胞核；實時熒光定量PCR 結果顯示，荷斯坦奶牛FGFR2 基因在乳腺組織中相對表達量最高。研究結果為荷斯坦奶牛FGFR2 基因在奶牛產奶性狀相關調控機制研究提供理論依據。