棘腹蛙生長激素/類胰島素樣生長因子軸關鍵基因的克隆與序列特征分析

羅 潔 姜玉松 楊 帆 樊汶樵 李曉英

(重慶文理學院林學與生命科學學院，重慶珍稀瀕危水產(chǎn)資源保護與開發(fā)研究中心，重慶402168)

棘腹蛙生長激素/類胰島素樣生長因子軸關鍵基因的克隆與序列特征分析

羅 潔 姜玉松 楊 帆 樊汶樵*李曉英

(重慶文理學院林學與生命科學學院，重慶珍稀瀕危水產(chǎn)資源保護與開發(fā)研究中心，重慶402168)

機體的生長激素(GH)/類胰島素樣生長因子(IGFs)軸由GH系統(tǒng)和IGFs系統(tǒng)構成，可促進細胞增殖、調(diào)節(jié)生長發(fā)育、調(diào)控生理代謝，在機體生長發(fā)育調(diào)控方面有著重要作用。為明確棘腹蛙GH/IGFs軸的功能結構和進化特征，為棘腹蛙生長發(fā)育調(diào)控方面的研究提供理論依據(jù)，本試驗對棘腹蛙GH、類胰島素生長因子-Ⅰ(IGF-Ⅰ)和類胰島素生長因子-Ⅱ(IGF-Ⅱ)進行克隆并對其序列特征進行分析。結果顯示：1)與兩棲類模式動物的多重序列比對發(fā)現(xiàn)，棘腹蛙GH、IGF-Ⅰ和IGF-Ⅱ的功能結構域嚴格保守，具有一定的遺傳多態(tài)性；IGF-Ⅱ的N端呈簡縮進化趨勢。2)遺傳進化聚類分析發(fā)現(xiàn)，棘腹蛙IGFs與兩棲動物聚為一支，并與硬骨魚相對近緣，說明IGF-Ⅰ和IGF-Ⅱ進化地位相對原始；棘腹蛙GH則與蛙類、魚類等水生動物相對近緣，暗示該基因具有趨同進化趨勢。3)為進一步明確上述基因的特異性功能位點，利用Swiss-model server軟件解析其蛋白質結構，最終確定IGF-Ⅰ的THR52、LEU53、PHE72、PHE73、SER74為潛在的功能分化位點，IGF-Ⅱ的TYR81、LYS82、LYS83為潛在的功能分化位點，GH的PHE208為潛在的功能分化位點。由此可知，棘腹蛙GH/IGFs軸的主要基因相對保守，但與已知模式物種相比，存在潛在的功能分化位點，可作為后期棘腹蛙GH/IGF軸功能研究和遺傳進化特征分析的分子靶標。

棘腹蛙；類胰島素生長因子；生長激素；序列分析；結構特征

生長激素(GH)是腺垂體細胞分泌的肽類激素，是影響動物生長和發(fā)育的主要內(nèi)源因素。類胰島素生長因子(IGFs)分為類胰島素生長因子-Ⅰ(IGF-Ⅰ)和類胰島素生長因子-Ⅱ(IGF-Ⅱ)2種類型，是存在于血漿內(nèi)的一類既有促生長作用又有胰島素樣作用的多肽。目前已闡明機體的生長是由GH/IGFs軸所控制，GH從垂體分泌后，通過GH受體的介導進而刺激肝臟和其他組織合成并分泌IGFs，后者通過IGFs受體的介導發(fā)揮生物學功能[1]。

棘腹蛙(Paaboulengeri)又名石坑、石蛙等，隸屬于無尾目，蛙科，棘蛙屬，是我國西部中高海拔區(qū)域特有的珍稀兩棲動物[2-3]。由于環(huán)境污染、生態(tài)破壞等原因，野生棘腹蛙的規(guī)模日趨萎縮，現(xiàn)已被《中國瀕危物種紅皮書》[4]、《中國物種紅色目錄》[5]等收錄。棘腹蛙目前市場需求旺盛，人工飼養(yǎng)的子二代出欄價達到45～70元/kg，其人工養(yǎng)殖已經(jīng)成為山區(qū)農(nóng)民脫貧致富的重要途徑。為更好地保護野生資源，同時滿足經(jīng)濟需求，棘腹蛙的人工養(yǎng)殖熱度逐年上升。

本課題組前期對棘腹蛙的人工飼養(yǎng)條件進行探索時發(fā)現(xiàn)，溫度過高或過低都會影響其生長發(fā)育，而溫度對GH/IGFs軸的影響尤為明顯。對棘腹蛙的遺傳特性研究發(fā)現(xiàn)棘腹蛙種群屬于一個單系分支[6-7]，在進化過程中相對獨立[8-9]，并無其他模式物種的相關數(shù)據(jù)可以參考。本研究擬從轉錄組測序獲得的信息入手，克隆GH/IGFs軸的關鍵基因(IGF-Ⅰ、IGF-Ⅱ和GH)并對其進行生物信息學分析，深入了解其基因特性，發(fā)掘物種特異性的潛在功能作用位點，為探索以棘腹蛙為代表的兩棲動物的生理反應和生長發(fā)育調(diào)節(jié)機制奠定理論基礎。

1 材料與方法

1.1 試劑與耗材

2歲齡健康棘腹蛙，飼養(yǎng)于重慶珍稀瀕危水產(chǎn)資源保護與開發(fā)研究中心兩棲動物流水養(yǎng)殖系統(tǒng)。Trizol、焦碳酸二乙酯(DEPC)購自上海生工生物工程有限公司，cDNA合成試劑盒、Taq酶、PCR純化試劑盒購自Promega公司，凝膠回收試劑盒購自Omega公司，DL2000 DNA Marker購自寶生物工程(大連)有限公司。

1.2 引物設計與合成

基因序列主要基于前期本實驗室建立的棘腹蛙蝌蚪的Illumina solexa第二高通量轉錄組數(shù)據(jù)庫，利用實時熒光定量PCR(RT-qPCR)分別獲取了棘腹蛙GH、IGF-Ⅰ、IGF-Ⅱ的全長基因序列，并經(jīng)過比對驗證其序列與轉錄組數(shù)據(jù)結果一致。最后，利用Primer 5.0軟件設計基因的特異性引物，所用引物如表1所示，由蘇州金唯智生物科技有限公司合成。

表1 引物設計Table 1 Primer design

1.3 cDNA制備

取雙毀髓法處死的棘腹蛙肝臟，立即放入液氮中。隨后總RNA按照試劑盒說明書進行提取，并利用Promega公司的cDNA合成試劑盒說明書完成單鏈cDNA(ss cDNA)的制備。

1.4 序列克隆與測定

利用高保真DNA Taq酶進行PCR擴增，RT-qPCR反應體系為：cDNA模板1 μL、10×Ex Taq PCR buffer 2.5 μL、Ex Taq 0.65 U、dNTP(100 μmol/L)1 μL、MgCl21 μL、上游引物1 μL、下游引物1 μL，加ddH2O至總體積25 μL。PCR反應程序：1)94 ℃預變性4 min；2)94 ℃變性90 s；3)58 ℃退火90 s；4)72 ℃延伸90 s；5)將步驟2～4重復循環(huán)30次；6)72 ℃終延伸10 min；7)12 ℃保存。最后，PCR產(chǎn)物經(jīng)1%瓊脂糖凝膠電泳進行檢測后切取與基因長度一致的片段，經(jīng)TaKaRa凝膠回收試劑盒純化后連接pMD 19-T載體，并轉化至大腸桿菌(E.coli)JM109感受態(tài)細胞中，利用藍白斑篩選陽性質粒，經(jīng)雙酶切驗證后送蘇州金唯智生物科技有限公司測序。

1.5 生物信息學分析

1.5.1 遺傳多態(tài)性分析方法

在NCBI數(shù)據(jù)庫中檢索GH、IGF-Ⅰ、IGF-Ⅱ的同源序列，利用ClustalX 1.83軟件進行多重序列比對和Boxshade Server序列對齊，最后結果由Photoshop CS6進行處理。

1.5.2 遺傳進化分析

在GenBank數(shù)據(jù)庫中檢索近緣物種序列，利用ClustalW對齊后，利用MEGA 5.0軟件對GH、IGF-Ⅰ、IGF-Ⅱ同源基因以鄰近(neighbor-joining)法構建系統(tǒng)進化樹，取1 000次重復檢驗以估算各分支的置信值。

1.5.3 蛋白質結構分析

首先利用Phyre2軟件(http://www.sbg.bio.ic.ac.uk/)檢索蛋白質模型，然后使用Swiss-model server軟件(https://swissmodel.expasy.org/)預測目的蛋白質的三維結構，最后采用Swiss-Pdb-viewer軟件進行圖像處理和蛋白質結構還原。

2 結果與分析

2.1IGF-Ⅰ的序列特征分析

以2歲齡棘腹蛙肝臟組織cDNA為模板進行

IGF-Ⅰ的基因擴增(圖1-A)及測序(圖1-B)，并將棘腹蛙IGF-Ⅰ與已發(fā)布的2種蟾蜍的近緣基因進行比較，發(fā)現(xiàn)前者保守的功能結構域存在較多的氨基酸位點突變(圖1-C)。

A:棘腹蛙IGF-Ⅰ的基因擴增；B:棘腹蛙IGF-Ⅰ的序列分析；C:已知兩棲動物的IGF-Ⅰ遺傳多態(tài)性分析。NP_001156865.1注釋為非洲爪蟾IGF-Ⅰ；XP_002936875.1注釋為熱帶爪蟾IGF-Ⅰ；IGF__1-like-Pb注釋為棘腹蛙IGF-Ⅰ。黃色表示棘腹蛙IGF-Ⅰ潛在功能分化位點。

A: the gene amplification ofIGF-Ⅰ fromPaaboulengeri; B: the sequence analysis ofIGF-Ⅰ fromPaaboulengeri; C: the genetic polymorphism analysis ofIGF-Ⅰ from known amphibians. NP_001156865.1 represents forIGF-Ⅰ fromXenopuslaevis; XP_002936875.1 represents forIGF-Ⅰ fromXenopustropicalis; IGF__1-like-Pb representsIGF-ⅠfromPaaboulengeri. The potential differentiation amino acid sites ofIGF-Ⅰ fromPaaboulengeriare shadowed in yellow.

圖1 棘腹蛙IGF-Ⅰ序列特征分析

Fig.1 The sequence characteristic analysis ofIGF-Ⅰ formPaaboulengeri

取不同代表物種IGF-Ⅰ進行聚類分析并構建系統(tǒng)進化樹(圖2)，發(fā)現(xiàn)IGF-Ⅰ主要分為兩大支，一支為高等哺乳動物，一支為爬行動物、鳥類和兩棲動物；其中，棘腹蛙IGF-Ⅰ與無尾目兩棲動物、部分水棲性爬行動物聚為一支，與其他陸生脊椎動物近緣，說明棘腹蛙IGF-Ⅰ與兩棲類模式動物爪蟾IGF-Ⅰ存在較大程度的遺傳分化。值得注意的是，揚子鱷IGF-Ⅰ有3個拷貝，其中2個拷貝與鳥綱近緣，另1個拷貝與中華鱉單獨聚為一支，與緬甸蟒蛇和龜殼花蛇近緣，推測脊椎動物的IGF-Ⅰ可能是由原始的爬行兩棲動物起源而來。

為深入分析棘腹蛙IGF-Ⅰ的蛋白質結構，我們利用Swiss-model server軟件進行結構預測，如圖3所示。相較于人類IGFs蛋白質結構3lri.1A[10]而言，兩者序列相似度達84.6%。其氨基酸序列中，氨基酸31～48處存在1個跨膜域，氨基酸51～109處為IGF-Ⅰ-like superfamily功能結構域。IGF-Ⅰ的蛋白質結構成分主要由2個β折疊和4個α螺旋構成，且兩端存在未知功能區(qū)域，其中THR52、LEU53、PHE72、PHE73、SER74為物種特異的氨基酸突變位點，暗示該蛋白質存在一定功能分化。

圖2 棘腹蛙IGF-Ⅰ的聚類分析Fig.2 The clustering analysis of IGF-Ⅰ form Paa boulengeri

A：棘腹蛙IGF-Ⅰ的蛋白質三維結構分析；B：棘腹蛙IGF-Ⅰ的特異功能結合位點預測；C：棘腹蛙IGF-Ⅰ與人類IGFs模板3lri.1.A[10]的蛋白質構象差異分析。

A: the 3D protein structure analysis of IGF-Ⅰ fromPaaboulengeri; B: the specific functional target prediction of IGF-Ⅰ fromPaaboulenger; C: the protein conformation difference analysis between IGF-Ⅰ fromPaaboulengeriand the IGFs template 3lri.1.A from human beings[10].

圖3 棘腹蛙IGF-Ⅰ的蛋白質結構分析

Fig.3 The protein structure analysis of IGF-Ⅰ from Paa boulengeri

2.2IGF-Ⅱ的序列特征分析

擴增棘腹蛙的IGF-Ⅱ基因并測序，利用ClustalX 1.83軟件對其遺傳多態(tài)性進行分析，結果如圖4所示。其氨基酸序列中，氨基酸54～114處為IGF-Ⅱ的功能結構域，氨基酸147～204處為IGF-Ⅱ的C端功能結構域，暗示該基因可能主要參與棘腹蛙的免疫防御。與非洲爪蟾和熱帶爪蟾2種已知的兩棲模式動物相比，棘腹蛙IGF-Ⅱ的拷貝數(shù)明顯減少，并且該基因的N端呈簡縮進化趨勢，且其功能結構域呈現(xiàn)插入突變，說明棘腹蛙IGF-Ⅱ的功能分化較為明顯。

A：棘腹蛙IGF-Ⅱ的基因擴增；B：棘腹蛙IGF-Ⅱ的序列分析；C：已知兩棲動物的IGF-Ⅱ遺傳多態(tài)性分析。XP_012816009、XP_012816008、NP_001107144、OCA36474.1注釋為熱帶爪蟾IGF-Ⅱ；NP_001085129、AAH72153.1、NP_001082128、AAH70545.1注釋為非洲爪蟾IGF-Ⅱ；IGF-2-Pb.aa注釋為 棘腹蛙IGF-Ⅱ。黃色表示棘腹蛙IGF-Ⅱ潛在功能分化位點，紅色下劃線表示N端功能分化區(qū)域。

A: the gene amplification ofIGF-Ⅱ fromPaaboulengeri; B: the sequence analysis ofIGF-Ⅱ fromPaaboulengeri; C: the genetic polymorphism analysis ofIGF-Ⅱ from known amphibians. NP_001085129, AAH72153.1, NP_001082128 and AAH70545.1 are stand forIGF-Ⅱ fromXenopuslaevis; XP_012816009, XP_012816008, NP_001107144 and OCA36474.1 are stand forIGF-Ⅱ fromXenopustropicalis; GF-2-Pb.aa representsIGF-ⅡfromPaaboulengeri. The potential differentiation amino acid sites ofIGF-Ⅱ fromPaaboulengeriare shadowed in yellow, and the N terminal differentiation region is underline in red.

圖4 棘腹蛙IGF-Ⅱ序列特征分析

Fig.4 The sequence characteristic analysis ofIGF-Ⅱ formPaaboulengeri

棘腹蛙IGF-Ⅱ同源序列的系統(tǒng)進化分析結果(圖5)顯示，整個進化樹分為三大支，一支為高等哺乳動物類，一支為兩棲動物類，另一支則包括硬

骨魚、陸生動物及部分鳥類，說明IGF-Ⅱ基因的進化相對保守。

利用Swiss-model server軟件模擬IGF-Ⅱ蛋白質的三維結構，結果顯示，其與棘腹蛙IGF-Ⅰ蛋白質的結構類似，主要由6個α螺旋構成，其中TYR81、LYS82、LYS83存在特異的氨基酸位點突變，暗示該蛋白質在物種進化過程中存在一定的功能分化。

對比IGF-Ⅰ和IGF-Ⅱ作用位點和結合方式差異的預測，我們利用PDB-Viewer還原了這2種蛋白質的蛋白質骨架，如圖3-A和圖6-A所示，黃色部分為IGF-Ⅰ和IGF-Ⅱ作用位點相同的區(qū)域，結果二者僅有3處區(qū)域有明顯差異，其蛋白質修飾方式的對比顯示所有差異位點均位于功能結合區(qū)域(圖3-B和圖6-B)，隨后對二者功能差異位點的比對結果顯示二者差異部分有2處為IGFs保守的功能結構域，其中一處為棘腹蛙IGF-Ⅰ特異的功能區(qū)域。由此可知，深入挖掘IGFs的功能對于研究棘腹蛙生理代謝的物種特異性具有重大意義。

2.3GH的序列特征分析

通過克隆測序(圖7-A)，獲得了棘腹蛙GH基因序列(圖7-B)。遺傳多態(tài)性分析結果發(fā)現(xiàn)棘腹蛙GH與已知的兩棲動物同源基因存在較少的氨基酸位點突變，并且其中大多數(shù)為同義氨基酸突變，對其功能結構影響不大(圖7-C)。

與GH同源基因進行系統(tǒng)進化樹聚類分析，結果如圖8所示。棘腹蛙的GH與牛蛙近緣，并與草魚、斑馬魚和淡水鰉等水生動物單獨聚為一支，暗示水生動物GH具有共同的祖先，鳥類與陸生動物聚為一支，高等哺乳動物單獨聚為一支，證實了棘腹蛙GH基因的進化相對保守。

對棘腹蛙GH蛋白質結構進行分析，結果如圖9所示。棘腹蛙GH具有典型的Hormone_1結構域(圖9-C)，屬于典型的4-螺旋細胞因子(4-helix cytokines)(圖9-A)。其中，PHE208與模式脊椎動物相比為突變位點(圖9-B和圖9-C)，而與兩棲模式動物相比則相對保守，說明棘腹蛙GH功能結構域相對保守。

3 討 論

在GH-IGFs軸中，GH處于上游位置，而IGFs則處于下游位置[12]。GH-IGFs軸的作用并不單一，是由多種因子和多種調(diào)控模式組成了GH合成與分泌的調(diào)控網(wǎng)絡[13]。GH/IGF-Ⅰ軸作為調(diào)控機體生長發(fā)育的主線[14]，其信號轉遞通路主要是通過啟動2條信號傳遞鏈，即磷酯酰肌醇-3激酶(PI3-K)激活途徑和絲裂原活化蛋白激酶(MAPK)激酶激活途徑，把有絲分裂和代謝信號傳遞到細胞核內(nèi)，從而啟動IGFs分泌，促進細胞增殖、分化以及抑制細胞的凋亡。IGF-Ⅰ與其受體結合后，首先導致胰島素底物-1(IRS-1)磷酸化，IRS-1被磷酸化后，PI3-K和生長因子結合蛋白2(Grb2)才能夠與其結合，由此啟動2條信號傳遞鏈。一條途徑是PI3-K被激活并形成磷酸化磷酸肌醇(PIP3)，PIP3就是細胞生長的信號，且PIP3途徑是抑制細胞凋亡的最經(jīng)典途徑。另一條途徑是激活細胞外信號調(diào)節(jié)激酶(ERK)，ERK把信號傳遞到細胞核內(nèi)啟動有絲分裂[15]。也就是說，IGF-Ⅰ一方面通過增加細胞的有絲分裂來促進機體的生長，另一方面則通過抑制細胞的凋亡來促進機體的生長。

A：棘腹蛙IGF-Ⅱ的蛋白質三維結構分析；B：棘腹蛙IGF-Ⅱ的特異功能結合位點預測；C：棘腹蛙IGF-Ⅱ與人類IGFs模板3lri.1.A[10]的蛋白質構象差異分析。

A： the 3D protein structure analysis of IGF-Ⅰ fromPaaboulengeri; B: the specific functional target prediction of IGF-Ⅰ fromPaaboulenger; C: the protein conformation difference analysis between IGF-Ⅰ fromPaaboulengeriand the IGFs template 3lri.1.A from human beings[10].

圖6 棘腹蛙IGF-Ⅱ的蛋白質結構分析

Fig.6 The protein structure analysis of IGF-Ⅱ fromPaaboulengeri

目前針對兩棲動物GH功能的報道較少，在僅有的部分早期研究中發(fā)現(xiàn)，兩棲動物變態(tài)前外源GH可刺激其幼體——蝌蚪的生長且在幼體的生長發(fā)育(尤其是肢的生長)中發(fā)揮特殊的作用[16]。外源性GH能促進蟾蜍幼體的生長，也能促進牛蛙幼體下肢的生長[17]。在蝌蚪發(fā)育的后期，血液中內(nèi)源性GH水平會不斷升高[16]。但內(nèi)源性GH在兩棲動物幼體生長、發(fā)育和變態(tài)中的確切作用仍然沒有相應的報道。本研究針對棘腹蛙內(nèi)源GH、IGF-Ⅰ和IGF-Ⅱ進行分析，結果發(fā)現(xiàn)三者與已知兩棲模式生物相比均存在較多氨基酸突變位點，說明其存在一定程度的功能分化。遺傳進化樹聚類分析顯示，與已知兩棲模式動物相比，棘腹蛙IGF-Ⅰ和IGF-Ⅱ均存在特異的物種分化趨勢，而GH的進化則相對保守。蛋白質結構分析顯示，與人類IGFs相比，棘腹蛙IGF-Ⅰ存在5處功能作用位點的分化區(qū)域，IGF-Ⅱ存在3處功能作用位點的分化區(qū)域；而棘腹蛙GH與模式脊椎動物GH相比，則僅有PHE208處存在氨基酸突變，說明該氨基酸位點對于棘腹蛙的功能分化具有重要研究價值。值得注意的是，該位點在兩棲動物中則相對保守，暗示PHE208是兩棲動物與高等哺乳動物功能分化的氨基酸位點。

綜上所述，盡管功能結構域和遺傳進化特征相對保守，但棘腹蛙IGF-Ⅰ的THR52、LEU53、PHE72、PHE73、SER74，IGF-Ⅱ的TYR81、LYS82、LYS83以及GH的PHE208為潛在的功能作用位點和結合位點。鑒于本課題組前期研究發(fā)現(xiàn)GH/IGFs軸易受到溫度影響，導致棘腹蛙蝌蚪無法變態(tài)甚至死亡，成蛙生長異常。因此，本研究為揭示溫度如何介導GH/IGFs軸代謝機制進而調(diào)控棘腹蛙生長發(fā)育篩選了分子靶標，并為今后棘腹蛙優(yōu)質品種選育提供理論依據(jù)。

A：棘腹蛙GH的擴增；B：棘腹蛙GH的序列分析；C：已知兩棲動物的GH遺傳多態(tài)性分析。XP_002933088.2注釋為熱帶爪蟾GH；AAP4356.1和CAA31038.1注釋為牛蛙GH；GH pre__Pb注釋為棘腹蛙GH。黃色表示棘腹蛙GH潛在功能分化區(qū)域。

A： the gene amplification ofGHfromPaaboulengeri; B: the sequence analysis ofGHfromP.boulengeri; C: the genetic polymorphism analysis ofGHfrom known amphibians. AAP4356.1 and CAA31038.1 both stand forGHfromRanacatesbeiana; XP_012816009 and XP_002933088.2 both stand forGHfromXenopustropicalis; GH pre__Pb stand forGHfromPaaboulengeri. The potential differentiation amino acid sites ofGHfromPaaboulengeriare shadowed in yellow.

圖7 棘腹蛙GH的序列特征分析

Fig.7 The sequence characteristic analysis ofGHformPaaboulengeri

圖8 棘腹蛙GH的聚類分析Fig.8 The clustering analysis of GH from Paa boulengeri

A：棘腹蛙GH的蛋白質三維結構分析；B：棘腹蛙GH的特異功能結合位點預測；C：棘腹蛙GH與人類GH模板1bp3.1.A[11]的蛋白質構象差異分析。

A： the 3D protein structure analysis of GH fromPaaboulengeri; B: the specific functional target prediction of IGF-Ⅰ fromPaaboulenger; C: the protein conformation difference analysis between GH fromPaaboulengeriand the GH template 1bp3.1.A from human beings[11].

圖9 棘腹蛙GH的蛋白質結構分析

Fig.9 The protein structure of GH fromPaaboulengeri

4 結 論

棘腹蛙GH/IGFs軸的主要基因相對保守，但與已知模式物種相比，存在潛在的功能分化位點，可作為后期棘腹蛙GH/IGFs軸功能研究和遺傳進化特征分析的分子靶標。

[1] POLLAK M.The insulin and insulin-like growth factor receptor family in neoplasia:an update[J].Nature Reviews Cancer,2012,12(3):159-169.

[2] 費梁,葉昌媛,黃永昭,等.中國兩棲動物檢索及圖解[M].成都:四川科學技術出版社,2005:137-139.

[3] YAN F,ZHOU W W,ZHAO H T,et al.Geological events play a larger role than Pleistocene climatic fluctuations in driving the genetic structure ofQuasipaaboulengeri(Anura: Dicroglossidae)[J].Molecular Ecology,2013,22(4):1120-1133.

[4] 汪松,趙爾宓.中國瀕危動物紅皮書(兩棲類和爬行類)[M].北京:科學出版社,1998:80-83.

[5] 汪松,解焱.中國物種紅色名錄:第1卷:紅色目錄[M].北京:高等教育出版社,2004:1-224.

[6] SHAN X,XIA Y,ZHENG Y C,et al.The complete mitochondrial genome ofQuasipaaboulengeri(Anura:Dicroglossidae)[J].Mitochondrial DNA, 2014,25(2):83-84.

[7] XIA Y,HU L J,SHAN X,et al.Isolation and characterization of eleven polymorphic tetranucleotide microsatellite loci forQuasipaaboulengeri(Anura:Dicroglossidae)[J].Conservation Genetics Resources,2013,5(1):5-7.

[8] YUAN S Q,XIA Y,ZHENG Y C,et al.Development of microsatellite markers for the spiny-bellied frogQuasipaaboulengeri(Anura:Dicroglossidae) through transcriptome sequencing[J].Conservation Genetics Resources,2014,7(1):229-231.

[9] QING L Y,XIA Y,ZHENG Y C,et al.A de novo case of floating chromosomal polymorphisms by translocation inQuasipaaboulengeri(Anura,Dicroglossidae)[J].PLoS One,2012,7(10):e46163.

[10] LAAJOKI L G,FRANCIS G L,WALLACE J C,et al.Solution structure and backbone dynamics of long-[Arg(3)]insulin-like growth factor-Ⅰ[J].The Journal of Biological Chemistry,2000,275(14):10009-10015.

[11] SOMERS W,ULTSCH M,DE VOS A M,et al.The X-ray structure of a growth hormone-prolactin receptor complex[J].Nature,1994,372(6505):478-481.

[12] PERRINI S,LAVIOLA L,CARREIRAM M C,et al.The GH/IGF1 axis and signaling pathways in the muscle and bone:mechanisms underlying age-related skeletal muscle wasting and osteoporosis[J].Journal of Endocrinology,2010,205(3):201-210.

[13] CANOSA L F,CHANG J P,PETER R E.Neuroendocrine control of growth hormone in fish[J].General and Comparative Endocrinology,2007,151(1):1-26.

[14] BARTKE A.Impact of reduced insulin-like growth factor-1/insulin signaling on aging in mammals:novel findings[J].Aging Cell,2008,7(3):285-290.

[15] JONES J I,CLEMMONS D R.Insulin-like growth factors and their binding proteins:biological actions[J].Endocrine Reviews,2002,16(1):3-34.

[16] MCLEAN E,DONALDSON E M.The role of growth hormone in the growth of poikilotherms[M]//SCHREIBMAN M.The endocrinology of growth,development,and metabolism in vertebrates.Amsterdam:Elsevier,1993:43-71.

[17] KIKUYAMA S,KAWAMURA K,TANAKA S,et al.Aspects of amphibian metamorphosis: hormonal control[J].International Review of Cytology,1993,145:105-148.

*Corresponding author, associate professor, E-mail: wonderbreeze@126.com

(責任編輯 菅景穎)

Cloning and Sequence Characteristic Analysis of Growth Hormone/Insulin-Like Growth Factors Axis Key Genes fromPaaboulengeri

LUO Jie JIANG Yusong YANG Fan FAN Wenqiao*LI Xiaoying

(ChongqingResearchCentersofConservationandDevelopmentonRare&EndangeredAquaticResources,CollegeofLifeScience&Forestry,ChongqingUniversityofArt&Science,Chongqing402168,China)

Growth hormone (GH)/insulin-like growth factors (IGFs) axis was composed of GH system and IGFs system, which can promote cell proliferation, regulate growth and development, regulate physiological metabolism, as well as played an important role in the regulation of growth and development. In order to clear the functional structure and evolution characteristics of GH/IGFs axis fromPaaboulengeriand to provide the theoretical basis for regulation study of growth and development forPaaboulengeri, theGH, insulin-like growth factor-Ⅰ (IGF-Ⅰ) and insulin-like growth factor-Ⅱ (IGF-Ⅱ) fromPaaboulengeriwere cloned and the sequence characteristics of them were analyzed. The results showed as follows: 1) multiple sequence alignment with amphibian model animals revealed thatGH,IGF-Ⅰ andIGF-Ⅱ displayed respective genetic polymorphisms within conserved functional domains, and the N-terminal domain ofIGF-Ⅱ showed a abbreviated evolution trend. 2) The phylogenic clustering analysis found thatIGFsfromPaaboulengeriand amphibians were classified to one cluster, and they had a near relationship with bony fish, which indicatedIGF-Ⅰ andIGF-Ⅱ were comparatively ancient genes.IGFsfromPaaboulengerihad a near relationship with aquatic animals such as frogs and fishes, which impliedGHhad a trend of convergent evolution. 3) To further define the specific functional sites of above genes, the protein structure of them were analyzed by Swiss-model server software, and some potential functional differentiation sites were found including THR52, LEU53, PHE72, PHE73 and SER74 from IGF-Ⅰ, TYR81, LYS82 and LYS83 from IGF-Ⅱ, and PHE208 from GH. The results indicate that the major genes of GH/IGFs axis fromPaaboulengeriare relatively conserved, but their exist some potential functional differentiation sites compared with known model species, which can be as molecular targets for the functional study and genetic characteristics analysis of GH/IGFs axis fromPaaboulengeriin future.[ChineseJournalofAnimalNutrition, 2017, 29(7)：2492-2501]

Paaboulengeri; IGFs; GH; sequence analysis; structure characteristics

10.3969/j.issn.1006-267x.2017.07.034

2016-12-07

重慶市基礎科學與前沿技術研究(重點)專項(cstc2015jcyjBX0013)；重慶市前沿與應用基礎研究(一般)項目(cstc2014jcyjA80042)；永川區(qū)自然科學基金項目(Ycstc,2015nc1001)；重慶文理學院人才引進項目(R2013LS14，R2014LX07)

羅 潔(1983—)，女，河南周口人，博士，從事微生物與免疫研究。E-mail: luojiecq@sina.cn

*通信作者：樊汶樵，副教授，E-mail: wonderbreeze@126.com

S917.4

A

1006-267X(2017)07-2492-10