川藏黑豬NR5A2基因克隆及生物信息學分析

（四川省畜牧科學研究院養豬研究所，四川省動物遺傳育種重點實驗室，四川成都 610000）

NR5A2屬于核受體超家族 NR5A亞家族，是一種重要的孤兒受體，與體細胞編程、動物胚胎發育及分化、類固醇激素生成和癌癥等有密切關系，在成年哺乳動物中NR5A2 主要在肝、腸和卵巢等組織中表達，特別是在卵巢組織中高表達。在卵巢組織中，NR5A2 在顆粒細胞和黃體細胞中特異表達，而在卵泡膜細胞和間質細胞中不表達，說明 NR5A2 在卵巢顆粒細胞和黃體細胞中發揮重要作用。牛卵巢優勢卵泡和排卵前卵泡顆粒細胞中 NR5A2 基因的表達水平極顯著高于小卵泡、黃體細胞，說明卵巢顆粒細胞中的 NR5A2 在哺乳動物卵泡發育和優勢化過程中發揮重要的作用。NR5A2 在調控哺乳動物卵泡發育、成熟和排卵等方面具有重要作用，是雌性動物繁殖力的關鍵調節因子。本研究以我國自主培育豬種川藏黑豬為研究對象，采用克隆測序技術獲得川藏黑豬NR5A2 基因編碼區序列，利用生物信息學方法分析其序列特征。

1 材料與方法

1.1 試驗動物

川藏黑豬母豬來自四川省畜牧科學研究院養豬科研基地，杜長大商品豬來自龍泉宏業屠宰場，各3頭，所有試驗豬經屠宰后立即采取卵巢組織樣，置于液氮中帶回實驗室，用于提取組織總RNA。

1.2 RNA提取和反轉錄

采用 Trizol 試劑（Invitrogen，USA）一步法提取卵巢組織總RNA，用反轉錄試劑盒（Promega，USA）反轉錄成 cDNA，具體步驟按說明書操作。反轉錄產物-30 ℃ 保存備用。

1.3 NR5A2基因ORF擴增

利用ORF Finder工具，尋找NR5A2基因可能的ORF，發現了全長為1506bp的完整ORF。在NCBI上尋找包含NR5A2基因ORF上下游mRNA的序列為模板設計引物引物如下：Forward primer：GCTGCACTTTTTCTTACGCCA

Reverse primer：TAGTGTTGCTTAGAGCCGTTCA

1.4 NR5A2 基因編碼區克隆和測序

以川藏黑豬卵巢組織 cDNA 為模板，擴增 NR5A2 基因編碼區序列。PCR 反應體系：1 μL 模板，1 μL dNTPs，1 μL 上、下游引物（10 μmol·L -1），2 μL 10 × Buffer，0.2 μL Taq 聚合酶（TaKaRa），加雙蒸水至終體積 20 μL。PCR 反應條件：94 ℃ 5 min；94 ℃ 1 min，退火 30 s，72 ℃ 1 min，35 個循環；72 ℃ 10 min。PCR產物經 0.8% 瓊脂糖凝膠電泳檢測，用回收試劑盒（Axygen，USA）純化，純化的目的片段與 pMD19-T 載體（TaKaRa）連接，連接產物轉化 DH5α 感受態細胞；挑取陽性克隆，用質粒提取試劑盒（Axygen，USA）提取質粒，由Invitrogen 公司進行測序。

1.5 NR5A2 基因的序列分析

川藏黑豬 NR5A2 基因編碼區序列的整理、翻譯和比對使用DNAStar 和 DNAMAN 軟件。用 ExPASy Pro-teomics Server在線軟件預測氨基酸的理化性質。利用BioEdit對川藏黑豬NR5A2氨基酸疏水性進行分析。采用SWISS-MODEL 在線分析軟件來進行預測NR5A2 蛋白的三級結構。

1.6 哺乳動物 NR5A2 基因進化分析

利用 BLAST 在線軟件 從 GenBank 數據庫中獲取其他哺乳動物NR5A2 基因編碼區序列，包括虎鯨（NW_004438658.1）、綿羊（JN182926）、大熊貓（NW_003217343.1）、黑猩猩（X M_001142695）、人（NM_205860）、牛（NM_001206816）、小鼠（NM_030676）、大鼠（NM_021742）等。利用 DNAStar 軟件進行序列的整理、比對。采用 DNASP 5.0軟件統計多態位點，分析核苷酸多樣度（π），進行中性檢驗，計算序列間的非同義替代數（dN）、同義替代數（dS）和 dN /dS 值。利用 MEGA 5.0 軟件計算序列間的轉換數（Ts）、顛換數（Tv）、轉換 / 顛換比值（Ts / Tv），并以雞（NM_205078）為外類群，采用最大似然法構建哺乳動物系統發育樹。

2 結果與分析

2.1 RNA電泳圖

川藏黑豬和商品豬組織總RNA電泳圖見圖1，根據設計的引物進行NR5A2基因PCR擴增，擴增產物電泳圖見圖2。

圖1 川藏黑豬和商品豬組織總RNA電泳圖

圖2 川藏黑豬和商品豬NR5A2基因PCR擴增產物電泳圖

2.2 NR5A2核酸堿基同源度分析

序列經過克隆、比對，剪切，獲得了全長為1506bp的完整ORF，通過測序得知其序列組成特征見表1。Pos#1、Pos#2、Pos#3分別代表密碼子第一、二、三位。從表1可看出，A+T含量低于G+C。同源性分析發現NR5A2基因的核苷酸序列與虎鯨、山羊、大熊貓、黑猩猩、人、牛、小鼠、大鼠等哺乳動物的一致性分別為93.0%、91.4%、83.9%、83.7%、90.4%、91.6%、75.4%、77.9%。川藏黑豬的NR5A2基因與引物源序列（XM_021064001）一致性99.87%，共有2個非保守位點T1209C、G1451C。

表1 川藏黑豬基因NR5A2編碼區核苷酸組成特征

2.3 NR5A2氨基酸序列同源度分析和頻率分析

NR5A2基因編碼區序列編碼包含了501個氨基酸殘基的蛋白質，如圖3所示。氨基酸使用頻率見表2。氨基酸序列同源比對發現川藏黑豬與虎鯨、牛、山羊、人、黑猩猩、大熊貓、大鼠、小鼠的一致性分別為86.5%、85.6%、85.6%、86.9%、86.9%、86.5%、79.8%和79.4%。

2.4 NR5A2蛋白理化性質和疏水性分析

ExPASy顯示NR5A2基因的分子量為57.3KD，等電點理論值為8.09，Met在N端，體外半衰期30h，蛋白的不穩定系數為51.60，預測顯示該蛋白結構不穩定。脂溶指數81.18，疏水性均值Grand average of hydropathicity（GRAVY）-0.385，具有較強的疏水性，見圖4。

圖3 川藏黑豬 NR5A2 基因的核苷酸序列和預測的氨基酸序列

圖4 NR5A2基因蛋白的疏水性

表2 NR5A2基因氨基酸頻率

圖5 NR5A2蛋白三級結果預測

2.5 NR5A2蛋白質三級結果預測

三級結構預測顯示川藏黑豬 NR5A2 蛋白主要由 11 個 α 螺旋和 2個 β 折疊組成，形成 ααααααββααααα 結構（圖5）。

表3 NR5A2基因堿基對組成頻數

2.6 NR5A2基因進化與系統發育分析

利用得到的川藏黑豬NR5A2基因完整編碼區和NCBI上下載的牛、山羊、大熊貓、人、黑猩猩、虎鯨、大鼠以及小鼠的NR5A2基因的編碼區序列進行比對和剪齊，其中保守位點（C）1125個，占總位點的74.7%；變異位點（V）381個，占總位點的25.3%；簡約信息位點（Pi）277個，占變異位點72.7%；單態位點（S）104個，占總變異位點的27.3%。堿基替換變異類型中轉換數106次，顛換數57次，轉換/顛換比值R為1.85。

圖6 基于NR5A2基因構建的哺乳動物的系統發育樹

根據哺乳動物 NR5A2 基因核甘酸序列，以鳥類中的雞為外類群，構建了哺乳動物 8 個物種的系統發育樹，結果見圖 6。哺乳動物 9個物種明顯地聚為 4 大類：偶蹄目中的牛、羊、川藏黑豬和鯨魚首先聚為一類；奇蹄目中的馬單獨聚為一類；靈長目中的人和黑猩猩聚為一類；嚙齒目中的大鼠和小鼠聚為一類，熊貓單獨聚為一類。在偶蹄目中，牛科中牛和綿羊首先聚為一類，然后與鯨魚聚在一起，再與豬聚在一起，可見系統發育樹與經典的分類結果是一致的。

3 討論

NR5A2目前發現的核受體超家族基因是核受體（nuclear receptors，NRs）超家族中的重要一員，成員大致可以分為 7 個亞家族（即 NR0 - NR6），它們在結構上均比較保守，一般由以下 4 個結構域構成：①A /B 結構域，含有一個 AF － 2 基序；②C 結構域，即 DNA 結合域（DNA-binding domain，DBD），這個結構域高度保守，由 2 個鋅指組成；③D 結構域，這個結構域充當 DBD 和 LBD 結構域間的鉸鏈；④E 結構域，即配體結合域（ligand-binding domain，LBD）。NR5A2 隸屬于 NR5 亞家族 A組，在 NR5A 的所有成員中，與 DNA 結合的特異性是由 DBD 結構域 C 末端 26 個氨基酸殘基組成的 Ftz － F1 盒（Ftz － F1 box）所決定研究發現 NR5A2 通過 CYP19A1 基因調節雌激素（E2）的合成，從而促進大鼠顆粒細胞的增殖，本文通過克隆測序獲得了川藏黑豬 NR5A2 基因編碼區序列，發現其與哺乳動物其他物種核苷酸的序列的同源性為 75%～93%，氨基酸序列的同源性為79%～87%，可見 NR5A2 基因在不同哺乳動物中具有一定的保守性。進一步分析發現，川藏黑豬 NR5A2 蛋白與哺乳動物其他物種一樣，都具有 DNA 結合域、配體結合域和 Ftz －F1 盒等保守結構域，推測川藏黑豬 NR5A2 蛋白與哺乳動物其他物種一樣在卵泡發育、成熟和排卵等方面具有重要作用。