王憲軍，向 華*，張煥容，任玉鵬，朱江江

(1.西南民族大學(xué)畜牧獸醫(yī)學(xué)院，成都 610041；2.青藏高原動物遺傳資源保護與利用四川省重點實驗室，成都 610041)

透明質(zhì)酸結(jié)合蛋白(hyaluronan binding protein 4,HABP4)又被稱為Ki-1/57，廣泛分布于哺乳動物組織中。HABP4蛋白含有多個RGG/RXR盒族，這些被看做是精氨酸甲基化的靶標(biāo)[1]。透射電鏡分析顯示，該蛋白不僅存在于細(xì)胞質(zhì)，還存在于細(xì)胞核和核孔，推測該蛋白亞細(xì)胞定位的多變性可能與該蛋白中存在精氨酸和賴氨酸殘基相關(guān)，這些殘基能與一些帶負(fù)電荷的生物大分子相結(jié)合，包括葡糖氨基葡聚糖、RNA、透明質(zhì)酸，因此，它也被稱為細(xì)胞內(nèi)透明質(zhì)酸結(jié)合蛋白，這也暗示了HABP4蛋白可能在細(xì)胞內(nèi)發(fā)揮著多種功能[2]。

HABP4蛋白最初在細(xì)胞內(nèi)被發(fā)現(xiàn)可作為單克隆抗體Ki-1的細(xì)胞內(nèi)交叉反應(yīng)物，且是第一個用于特異性檢測霍奇金淋巴瘤惡性細(xì)胞的抗體[3]。盡管HABP4和霍奇金病之間關(guān)系尚不清楚，但研究發(fā)現(xiàn)，HABP4蛋白參與了腫瘤的發(fā)生，如Gray-Mcguire等[4]在分析家族性結(jié)腸癌的風(fēng)險相關(guān)SNP中發(fā)現(xiàn)，其和HABP4的SNP存在強烈的不平衡性。此外，HABP4蛋白作為轉(zhuǎn)錄因子復(fù)合體，可通過與多種蛋白相互作用，在細(xì)胞轉(zhuǎn)錄控制與基因表達的調(diào)控中起著至關(guān)重要的作用。Lemos等[5-7]發(fā)現(xiàn)，HABP4蛋白與染色質(zhì)重塑蛋白(chromodomain helicase DNA binding protein 3,CHD3)、轉(zhuǎn)錄因子(myocyte enhancer factor 2C,MEF2C和p53)相互作用后在生物體轉(zhuǎn)錄調(diào)控方面直接或者間接地發(fā)揮作用。Bressan等[8-11]報道，HABP4蛋白與hnRNPQ(heterogeneous nuclear ribonucleoprotein)和SFRS9(arginine/serine-rich 9)蛋白作用后調(diào)控mRNA前體的拼接，以及與 RACK1(receptor for activated C kinase 1)、CIRP(cold inducible RNA binding protein)、FMRP(fragile X mental retardation protein)和RPL38(ribosomal protein L38)蛋白作用后調(diào)控蛋白的翻譯。Bressan等[8]進一步研究發(fā)現(xiàn)，HABP4蛋白通過與富U RNA探針的結(jié)合參與RNA的新陳代謝過程，且具有深入調(diào)節(jié)mRNA前體拼接和翻譯方面的功能。2013年，Braunschweig等[12]發(fā)現(xiàn)，HABP4蛋白在mRNA的拼接、核染色質(zhì)重建、轉(zhuǎn)錄和翻譯中發(fā)揮作用，暗示此蛋白在調(diào)控整合基因的表達過程中具有重要作用。同時，HABP4蛋白參與了細(xì)胞內(nèi)RNA的代謝作用[8]。

近年研究發(fā)現(xiàn)，通過微陣列數(shù)據(jù)分析顯示，HABP4蛋白的表達調(diào)控了人源細(xì)胞周期、凋亡、增殖相關(guān)基因表達水平[13]。HABP4蛋白在真核細(xì)胞內(nèi)本質(zhì)上屬于非結(jié)構(gòu)蛋白家族的成員之一，其結(jié)構(gòu)的靈活性也暗示了該蛋白能與許多不同的蛋白質(zhì)發(fā)生相互作用[14]，推測HABP4蛋白可能通過對生物體轉(zhuǎn)錄翻譯過程的調(diào)節(jié)從而參與細(xì)胞周期、細(xì)胞凋亡、mRNA代謝等[15-17]。由此看出，HABP4蛋白在生物機體內(nèi)發(fā)揮著重要的作用?；诖?，研究人員對該蛋白產(chǎn)生了濃厚的興趣并進行重點關(guān)注。然而，目前尚未見HABP4蛋白在調(diào)控山羊細(xì)胞凋亡中作用的相關(guān)報道。

簡州大耳羊是我國自主培育的第二個肉用山羊品種，具有生長發(fā)育快、適應(yīng)性好、抗病力強等特點，充分挖掘簡州大耳羊的品種優(yōu)勢對推動肉用山羊產(chǎn)業(yè)的發(fā)展具有重要意義。本試驗以簡州大耳羊為研究對象，利用TA法克隆得到山羊HABP4基因，并利用生物信息學(xué)方法預(yù)測其結(jié)構(gòu)及功能；利用實時熒光定量PCR(RT-qPCR)技術(shù)分析其在不同組織中的表達；同時，使用RT-qPCR技術(shù)檢測與HABP4蛋白可能存在相互作用的10個蛋白的mRNA表達特性并分析其相關(guān)性，進行HABP4基因及相關(guān)基因的GO功能富集，以期為進一步驗證HABP4在調(diào)控山羊細(xì)胞凋亡中的作用奠定基礎(chǔ)。

1 材料與方法

1.1 樣本采集及總RNA提取

試驗于2020年在西南民族大學(xué)畜牧獸醫(yī)學(xué)院完成。于四川省簡陽大哥大牧業(yè)有限公司隨機選取10月齡左右，體重為55 kg左右的健康簡州大耳羊公羊16只，作為16個生物學(xué)重復(fù)樣本。清晨空腹屠宰后，迅速采集各只山羊心、肝、脾、肺、腎、小腸、大腸、瘤胃和胰腺9個組織樣品，DEPC水清洗后將組織樣品用錫箔紙包好，立刻置于液氮中保存，待用。使用TRIzol法提取組織總RNA，并用1%的瓊脂糖凝膠電泳檢測RNA完整性。采用NanoDropND-1000(Agilent)分光光度計測定RNA樣品的濃度以及OD260 nm/OD280 nm值，并控制在1.9～2.1之間。利用gDNAse去除基因組DNA后，使用QuantiTect Reverse Transcription Kit進行cDNA鏈的合成。1%的瓊脂糖凝膠電泳檢測cDNA的質(zhì)量。

1.2 主要試劑及儀器

TRIzol、Revert Aid First Strand cDNA Synthesis Kit 試劑盒、QuantStudio3購自美國Thermo Fisher Scientific公司；Premix Taq、DL2000 DNA Marker、pMD-19T Simple載體、DH5α感受態(tài)細(xì)胞購自寶生物工程(大連)有限公司(TaKaRa)；SYBR?Premix Ex TaqTM(2×)購自TaKaRa公司；膠回收試劑盒、質(zhì)粒提取試劑盒購自美國Omega公司；氯仿、異丙醇、DEPC等均為國產(chǎn)試劑。

1.3 山羊HABP4基因的克隆及鑒定

參照GenBank 中已公布的山羊HABP4基因預(yù)測序列(登錄號：XM_018052483.1)，利用Primer Premier 5軟件設(shè)計1對特異性引物，上游引物:5′-GGCTTCGCTAAGACTCGC-3′，下游引物：5′-G-TGTCTGGAGAAATCCCTGTT-3′，引物由上海生工生物工程有限公司合成。以山羊脾cDNA為模板，使用上、下游引物進行PCR擴增。PCR反應(yīng)體系：Premix Taq 12.5 μL，上、下游引物各1 μL，模板1 μL，ddH2O 9.5 μL。反應(yīng)條件：95 ℃預(yù)變性4 min；95 ℃變性30 s，60 ℃退火30 s，72 ℃延伸1 min，共35個循環(huán)；72 ℃延伸7 min。PCR產(chǎn)物經(jīng)1%瓊脂糖凝膠電泳進行檢測正確后，切膠回收產(chǎn)物連接至pMD-19T Simple載體，并轉(zhuǎn)化DH5α感受態(tài)細(xì)胞，提取質(zhì)粒后陽性重組質(zhì)粒送上海生工生物工程有限公司測序正確后保存?zhèn)溆?，將鑒定正確的質(zhì)粒命名為pMD19-HABP4。

1.4 山羊HABP4基因的生物信息學(xué)分析

使用clustalw2 在線軟件(http://www.ebi.ac.uk/Tools/clustalw2/index.html)對本試驗克隆獲得的山羊HABP4基因CDS序列與綿羊預(yù)測序列(XM_004005238.4)、牛(NM_001081523.1)、豬(NM_001244681.1)、人(NM_014282.4)、小鼠(NM_019986.3)和大鼠(NM_001079940.1)的HABP4基因進行同源性分析。使用BioXM2.6軟件對山羊與綿羊(預(yù)測序列)、牛、人和豬進行HABP4基因序列比對。使用Protparam在線軟件(https://web.expasy.org/protparam/)對HABP4的氨基酸序列一級結(jié)構(gòu)和理化性質(zhì)進行分析。使用ExPASy ProtScale在線軟件(https://web.expasy.org/protscale/)計算HABP4的疏水性。使用 Tmpred 在線軟件(http://www.ch.embnet.org/software/TMPRED_form.html)分析蛋白質(zhì)的跨膜區(qū)。利用EMBL-EBI在線軟件(https://www.ebi.ac.uk/Tools/pfa/phobius)進行信號肽預(yù)測。使用phyre2在線軟件 (http://www.sbg.bio.ic.ac.uk/phyre2/phyre2)預(yù)測山羊HABP4蛋白質(zhì)二級結(jié)構(gòu)。使用SWISS MODEL在線軟件(https://swissmodel.expasy.org/interactive)分析蛋白的三級結(jié)構(gòu)。使用Net Phos 3.1 (http://www.cbs.dtu.dk/services/NetPhos/)預(yù)測HABP4蛋白的磷酸化位點。使用在線軟件PsortⅡ(https://www.genscript.com/tools/psort)進行蛋白亞細(xì)胞定位預(yù)測。使用 MEGA7.0 軟件對HABP4 基因序列構(gòu)建系統(tǒng)進化樹。利用STRING在線軟件(https://string-db.org/)預(yù)測可能與HABP4蛋白發(fā)生相互作用的蛋白。

1.5 山羊HABP4基因在不同組織的表達分析

選取10月齡左右健康山羊(n=16)的9個組織(心、肝、脾、肺、腎、小腸、大腸、瘤胃和胰腺)，分別對每個組織提取總RNA并反轉(zhuǎn)錄為cDNA。根據(jù)“1.3”中克隆的HABP4基因CDS序列，運用Primer Premier5設(shè)計1對特異性引物，RT-qPCR檢測HABP4基因在各組織中的表達差異。以實驗室前期保存的Ubiquitously-Expressed Transcript(UXT)基因作為內(nèi)參基因，內(nèi)參基因的引物序列見表1。RT-qPCR的反應(yīng)體系：SYBR?Premix Ex TaqTM10 μL，上、下游引物各1 μL，模板1 μL，ddH2O 6.8 μL，ROX 0.2 μL。反應(yīng)條件：95 ℃ 預(yù)變性3 min；95 ℃變性10 s，60 ℃退火30 s，同時采集熒光信號，72 ℃延伸30 s，共35個循環(huán)，熔解曲線：95 ℃變性15 s，60 ℃冷卻1 min，使DNA雙鏈充分結(jié)合，最后從60 ℃逐步加熱至95 ℃，每一步增加0.15 ℃并保持1 s，同時采集熒光信號。

1.6 RT-qPCR檢測HABP4互作蛋白的mRNA表達水平

根據(jù)STRING數(shù)據(jù)庫預(yù)測可能與HABP4蛋白發(fā)生相互作用的蛋白為基礎(chǔ)，利用RT-qPCR檢測與HABP4蛋白發(fā)生相互作用的蛋白的mRNA表達。各基因引物見表1。RT-qPCR反應(yīng)體系：SYBR?Premix Ex TaqTM10 μL，上、下游引物各1 μL，模板1 μL，ddH2O 6.8 μL，ROX 0.2 μL。反應(yīng)條件：95 ℃預(yù)變性3 min；95 ℃變性10 s，60 ℃退火30 s，同時采集熒光信號，72 ℃延伸30 s，共35個循 環(huán)，熔解曲線：95 ℃變性15 s，60 ℃冷卻1 min，使DNA雙鏈充分結(jié)合，最后從60 ℃逐步加熱至95 ℃，每一步增加0.15 ℃并保持1 s，同時采集熒光信號。每個待測樣本設(shè)置3個重復(fù)。

表1 引物信息Table 1 Primers information

1.7 數(shù)據(jù)統(tǒng)計與分析

使用Metascape在線分析軟件(http://metascape.org/)對HABP4進行GO功能富集。分析HABP4基因在山羊不同組織中表達水平時，采用2-△△Ct方法分析HABP4基因相對于內(nèi)參基因UXT的表達差異，其中△Ct=CtHABP4-CtUXT，△△Ct=△Ct組織-△Ct心。計算結(jié)果在SPSS 22.0中進行單因素方差分析和顯著性檢驗。采用2-△Ct方法分析HABP4基因互作蛋白在胰腺中的表達，其中△Ct=Ct預(yù)測互作蛋白目的基因-CtUXT。并將計算結(jié)果在SPSS 22.0中進行相關(guān)性分析，P<0.05表示差異顯著，P<0.01表示差異極顯著。

2 結(jié) 果

2.1 山羊HABP4基因的克隆及鑒定

本試驗以山羊組織cDNA為模板，通過PCR方法成功擴增HABP4基因，經(jīng)1%的瓊脂糖凝膠電泳檢測可見與預(yù)期目的片段大小(1 496 bp)相符合(圖1A)，包含5′UTR 61 bp，CDS區(qū)1 254 bp，3′UTR 181 bp，編碼417個氨基酸，起始密碼子為ATG，終止密碼子為TGA。上傳HABP4基因CDS區(qū)序列到NCBI獲得GenBank登錄號：MW748177。膠回收產(chǎn)物克隆轉(zhuǎn)化后，進行菌落PCR，產(chǎn)物經(jīng)1%瓊脂糖凝膠電泳檢測，與預(yù)測HABP4基因的大小相符合(圖1B)。

A.HABP4基因的PCR擴增檢測；B.重組質(zhì)粒pMD19-HABP4的PCR擴增。M.DNA相對分子質(zhì)量標(biāo)準(zhǔn)；N.陰性對照；1～2.HABP4基因的擴增條帶；3～6.重組質(zhì)粒pMD19-HABP4擴增條帶A.HABP4 gene amplified by PCR;B.Recombinant plasmid pMD19-HABP4 amplified by PCR.M.2 000 bp DNA ladder;N.Negative control;1-2.HABP4 gene amplified by PCR;3-6.Recombinant plasmid pMD19-HABP4 amplified by PCR圖1 HABP4基因克隆及重組質(zhì)粒鑒定Fig.1 HABP4 gene cloning and recombinant plasmid identification

2.2 HABP4基因的生物信息學(xué)分析

本研究獲得的HABP4基因CDS區(qū)的同源性分析表明，山羊HABP4基因的CDS區(qū)與NCBI中牛(NM_001081523.1)的同源性高達98.41%，與綿羊預(yù)測序列(XM_004005238.4)的同源性為100%，與豬(NM_001244681.1)、人(NM_014282.4)、小鼠(NM_019986.3)和大鼠(NM_001079940.1)的同源性分別為94.03%、89.49%、84.64%、85.05%。使用BioXM2.6軟件對山羊與綿羊(預(yù)測序列)、牛、豬和人進行HABP4基因序列比對，結(jié)果見圖2。

利用Protparam在線軟件對克隆獲得的山羊HABP4基因CDS區(qū)的氨基酸序列進行分析，預(yù)測其蛋白分子式為C1 957H3 108N638O644S8，分子質(zhì)量為46 134.57 u；理論等電點(theoretical PI)為6.54，不穩(wěn)定指數(shù)為58.05，該蛋白為不穩(wěn)定蛋白。ExPASy ProtScale在線軟件計算HABP4的疏水性(圖3A)，推測該蛋白為親水性蛋白。使用Tmpred Server v.2.0 對HABP4蛋白的跨膜區(qū)進行分析，推測該蛋白不屬于跨膜蛋白(圖3B)。使用 SIGNALP 4.0分析HABP4蛋白，預(yù)測該蛋白不存在信號肽結(jié)構(gòu)，推測該蛋白為非分泌性蛋白 (圖3C)。Net Phos 3.1在線軟件預(yù)測HABP4蛋白可能的磷酸化位點，結(jié)果顯示該蛋白有19個絲氨酸磷酸化位點、11個蘇氨酸磷酸化位點和3個絡(luò)氨酸磷酸化位點，無O-糖基化位點和N-糖基化位點。亞細(xì)胞定位分析顯示，HABP4蛋白主要定位于細(xì)胞核(73.9%)和細(xì)胞質(zhì)(8.7%)中。Phyre2在線軟件對山羊HABP4蛋白進行二級結(jié)構(gòu)分析，預(yù)測結(jié)果顯示HABP4蛋白中α螺旋占19%，無規(guī)則卷曲占83%，β折疊占4%。利用 Phyre2 Server預(yù)測了山羊HABP4蛋白的三級結(jié)構(gòu)(圖3D)。同時分析了與山羊HABP4蛋白氨基酸同源性最高的綿羊(預(yù)測序列)HABP4蛋白三級結(jié)構(gòu)(圖3E),分析比對發(fā)現(xiàn)兩個物種HABP4蛋白三級結(jié)構(gòu)無明顯差異，與豬HABP4蛋白三級結(jié)構(gòu)(圖3F)對比也無明顯差異，預(yù)測其功能具有一定相似性；與牛HABP4蛋白三級結(jié)構(gòu)(圖3G)對比分析發(fā)現(xiàn)有差異性，預(yù)測其功能相似度較低。利用 MEGA7.0軟件構(gòu)建各物種HABP4基因序列系統(tǒng)進化樹。如圖4A顯示，山羊與牛親緣關(guān)系最近，與豬的親緣關(guān)系較遠，與小鼠親緣關(guān)系最遠。利用Sting數(shù)據(jù)庫檢索到MRPS7、CHD4、UBA52(ubiquitin A-52 residue ribosomal protein fusion product 1)、RPL32(ribosomal protein L32)、GNB2L1(guanine nucleotide-binding protein subunit beta-2-like 1,)、RPS9(ribosomal protein S9)、RPS3(ribosomal protein S3,)、RPS16(ribosomal protein S16)、RPL31(ribosomal protein L31,)、RPL38共10個蛋白可能和HABP4蛋白存在相互作用(圖4B)。

A.山羊HABP4蛋白的疏水性分析；B.山羊HABP4蛋白的跨膜區(qū)分析；C.山羊HABP4蛋白的信號肽預(yù)測；D.山羊HABP4蛋白的三級結(jié)構(gòu)預(yù)測；E.綿羊HABP4蛋白的三級結(jié)構(gòu)預(yù)測；F.豬HABP4蛋白的三級結(jié)構(gòu)預(yù)測；G.牛HABP4蛋白的三級結(jié)構(gòu)預(yù)測A.Hydrophobicity analysis of goat HABP4;B.Transmembrane region analysis of goat HABP4 protein;C.Predicted signal peptide of goat HABP4 protein;D.Predicted tertiary structure of goat HABP4 protein;E.Predicted tertiary structure of sheep HABP4 protein;F.Predicted tertiary structure of pig HABP4 protein;G.Predicted tertiary structure of cattle HABP4 protein圖3 山羊HABP4蛋白的生物信息學(xué)分析Fig.3 Bioinformatics analysis of goat HABP4 protein

A.HABP4基因進化樹分析；B.與HABP4可能存在互作的蛋白預(yù)測A.Evolutionary tree analysis of HABP4 gene;B.Prediction of proteins interacting with HABP4 protein圖4 山羊HABP4蛋白的生物信息學(xué)分析Fig4 Bioinformatics analysis of goat HABP4 protein

2.3 山羊HABP4基因不同組織的mRNA表達水平分析

以UXT為內(nèi)參基因，采用RT-qPCR技術(shù)檢測山羊9個不同組織中HABP4基因mRNA的表達水平。結(jié)果顯示，在山羊的心、肝、脾、肺、腎、胰腺、小腸、大腸和瘤胃等9個組織中均能檢測到HABP4基因的mRNA表達，其中HABP4基因mRNA表達水平在胰腺中最高，極顯著高于其他組織(P<0.01)，腎組織次之，肝和瘤胃最低(圖5)。

圖5 HABP4基因mRNA在不同組織中相對表達量Fig.5 Relative expression of HABP4 mRNA in different tissues

2.4 HABP4互作蛋白預(yù)測及各基因表達量間的相關(guān)性分析

根據(jù)STRING在線軟件預(yù)測到HABP4蛋白與MRPS7、CHD4、UBA52、RPL32、GNB2L1、RPS9、RPS3、RPS16、RPL31、RPL38共10個蛋白可能存在相互作用，因此推測，這些基因可能與HABP4的表達水平具有一定的相關(guān)性。以UXT為內(nèi)參基因，利用RT-qPCR檢測這10個基因在16只羊胰腺中的mRNA相對表達水平。結(jié)果見表2，HABP4與UAB52(r=0.642,P<0.01)具有極顯著正相關(guān)，與RPL32(r=0.564,P<0.05)、RPS9(r=0.589,P<0.05)具有顯著正相關(guān)性；MRPS7和RPS3(r=0.910,P<0.01)、RPL32(r=0.862,P<0.01)、RPS9(r=0.850,P<0.01)具有極顯著正相關(guān)，與RPS16(r=0.570,P<0.05)具有顯著正相關(guān)性；CDH4和RPL31(r=0.792,P<0.01)、RPS16(r=0.808,P<0.01)、GNB2L1(r=0.959,P<0.01)具有極顯著正相關(guān)，與RPS3(r=0.501,P<0.05)具有顯著正相關(guān)性；UBA52與RPL32(r=0.596,P<0.05)具有顯著正相關(guān)性；RPS3與RPS16(r=0.796,P<0.01)、RPL32(r=0.710,P<0.01)、GNB2L1(r=0.653,P<0.01)、RPS9(r=0.672,P<0.01)具有極顯著正相關(guān)；RPL38和RPL31(r=0.710,P<0.01)具有極顯著正相關(guān),與RPS16(r=0.619,P<0.05)具有顯著正相關(guān)；RPL31與RPL32(r=-0.395,P<0.05)、RPS9(r=-0.399,P<0.05)具有顯著負(fù)相關(guān)；RPS16與RPL32(r=0.298,P<0.01)、RPS9(r=0.277,P<0.01)具有極顯著正相關(guān)；RPL32與GNB2L1(r=0.021,P<0.01)具有極顯著正相關(guān)；GNB2L1與RPS9(r=0.040,P<0.01)具有極顯著正相關(guān)，相關(guān)性網(wǎng)絡(luò)分析發(fā)現(xiàn)，各個基因表達量之間均具有較強的相關(guān)性，如RPS9與7個基因具有較強的相關(guān)性，其次RPL32、RPS16、RPS3各與6個 基因具有較強相關(guān)性，GNB2L1、RPL31、CDH4、MRPS7各與4個基因具有較強相關(guān)性，UBA52、HABP4各與3個基因具有較強相關(guān)性，最少的為RPL38與2個基因具有較強相關(guān)性。

表2 與 HABP4互作蛋白的基因表達量之間相關(guān)性分析Table 2 Correlation analysis of genes expression among proteins interacting with HABP4

2.5 HABP4蛋白及相關(guān)性基因的GO功能富集

從2.4結(jié)果中得知HABP4基因與UAB52、RPL32、RPS9共3個基因具有顯著相關(guān)性，對其進行GO分析，以富集結(jié)果P<0.01 為參考標(biāo)準(zhǔn)。GO 分析包括生物過程(圖6A)和分子功能(圖6 B)兩個方面。在生物過程方面主要涉及細(xì)胞定位、新陳代謝、多生物過程以及生物過程負(fù)調(diào)節(jié)；分子功能方面與翻譯起始、核糖體轉(zhuǎn)錄、細(xì)胞質(zhì)中的代謝有關(guān)。

A.生物學(xué)過程分析；B.分子功能分析A.Biological process analysis;B.Molecular function analysis圖6 山羊HABP4及相關(guān)基因功能富集分析Fig.6 Functional enrichment analysis of goat HABP4 and related genes

3 討 論

本研究獲得山羊HABP4基因的CDS區(qū)域以及預(yù)測其生物信息學(xué)相關(guān)信息對揭示該基因的功能將起到重要的作用。HABP4基因CDS區(qū)為1 254 bp，其中GC含量達60%以上，由于HABP4基因的GC含量較高，在擴增過程中導(dǎo)致多條非特異性條帶的產(chǎn)生。本試驗通過多次優(yōu)化反應(yīng)條件，同時使用PrimeSTAR Max Premix高保真酶(TaKaRa)克隆獲得了條帶單一的PCR產(chǎn)物。測序結(jié)果表明，本試驗克隆得到的HABP4基因的CDS區(qū)與NCBI上公布的牛、豬、人、小鼠、大鼠的同源性分別為98.41%、94.03%、89.49%、86.64%和85.05%，可以看出，HABP4在哺乳動物間存在較高的保守性并且與牛的同源性最高。不同的生物信息學(xué)分析軟件預(yù)測基因功能的結(jié)果可能有所不同，因此，篩選好的生物信息學(xué)分析軟件對精確預(yù)測基因功能是至關(guān)重要的。如MHMM-2.0和TMPred兩款在線比對軟件常用于蛋白質(zhì)跨膜區(qū)域的預(yù)測[18]。本研究采用ProtScale在線軟件預(yù)測蛋白的親疏水性，發(fā)現(xiàn)采用試驗克隆獲得的HABP4基因序列和其氨基酸序列得到不同的結(jié)果，其中基因序列結(jié)果分析顯示為疏水性蛋白，氨基酸序列分析結(jié)果為親水性蛋白。李武峰等[19-21]均采用目的基因的氨基酸序列來分析該蛋白對應(yīng)的疏水性。可以看出，利用氨基酸序列分析蛋白的疏水性的可靠度更高。因此，本試驗使用HABP4基因的氨基酸序列通過ProtScale在線軟件，預(yù)測結(jié)果顯示HABP4蛋白為親水性蛋白。HABP4蛋白不屬于跨膜蛋白，無明顯的跨膜結(jié)構(gòu)，這也與其在細(xì)胞質(zhì)中發(fā)揮作用的研究結(jié)果是一致的。蛋白的磷酸化位點和糖基化位點可為研究HABP4蛋白的活性區(qū)域，并進行定點突變具有重要的指導(dǎo)意義，且這些位點常常對生物進化的研究具有重要的參考價值[22]。蛋白的亞細(xì)胞定位往往能一定程度上反映出該蛋白的生物學(xué)功能[23]。本研究亞細(xì)胞定位分析顯示，HABP4蛋白主要定位于細(xì)胞核(73.9%)和細(xì)胞質(zhì)(8.7%)中，這與Rohde等[24]的研究報道結(jié)果一致，推測山羊HABP4蛋白可能與人源HABP4蛋白的功能相似。相互作用蛋白分析顯示，HABP4蛋白可能與MRPS7、CHD4、UBA52等蛋白存在相互作用。

RT-qPCR分析mRNA表達水平時，合適內(nèi)參基因的選擇能得到相對準(zhǔn)確的RT-qPCR數(shù)據(jù)[25]。內(nèi)參基因是一種標(biāo)準(zhǔn)化核酸樣品數(shù)量的基因，從理論上其表達量應(yīng)保持穩(wěn)定。但有研究得出，常用的內(nèi)參基因存在波動性表達，穩(wěn)定表達只存在于特定的試驗條件、細(xì)胞或組織中[26]，因此，篩選表達水平相對穩(wěn)定的內(nèi)參基因顯得尤為重要[27-28]。常用于研究的內(nèi)參基因主要包括GAPDH、ACTIN、18SrRNA、RPLP0、TBP、UXT等，其中GAPDH長期作為公認(rèn)的內(nèi)參基因被大家廣泛使用[29]，如2017年，張蓉等[27]研究發(fā)現(xiàn)，GAPDH在蛋雞組織表達中為最優(yōu)選擇的內(nèi)參基因。池永東等[30]報道，相同的內(nèi)參基因并不適合于所有的物種，且發(fā)現(xiàn)GAPDH基因在山羊組織中并非穩(wěn)定性最強。2013年，Bonnet等[31]研究發(fā)現(xiàn)，UXT、EIF3K和RPLP0基因是?；蛏窖蚋巍⒓∪?、乳腺和脂肪組織中最為穩(wěn)定的參考基因，可使mRNA豐度正常化。同時，2020年，謝光杰等[32]以UXT為內(nèi)參基因，利用RT-qPCR技術(shù)檢測了山羊APOL6基因在簡州大耳羊 背最長肌、皮下脂肪和腹間脂肪的表達量。許晴等[33-35]在山羊基因功能研究中以UXT為內(nèi)參基因矯正目的基因的表達。因此，本研究結(jié)合前人研究結(jié)果，選擇UXT作為內(nèi)參基因。

分析基因在不同組織中的表達規(guī)律可為研究該蛋白的功能提供試驗依據(jù)。從山羊HABP4的組織分布來看，該基因在各組織中均有表達，其中在胰腺中的表達量最高，腎組織次之，肝和瘤胃最低。HABP4在各組織間的表達差異可能和其在不同組織中的功能差異相關(guān)。研究顯示，胰腺一直以一種特殊的實質(zhì)器官被研究，胰腺的內(nèi)分泌及外分泌功能向來被認(rèn)為是相互獨立存在[36],在細(xì)胞轉(zhuǎn)錄調(diào)控、mRNA代謝、RNA拼接以及mRNA穩(wěn)定性等過程中發(fā)揮著作用，而HABP4蛋白同樣在基因轉(zhuǎn)錄調(diào)控、mRNA代謝等水平發(fā)揮作用，從而推測HABP4蛋白可能在山羊的胰腺中發(fā)揮著類似重要的生物學(xué)功能。

目前，關(guān)于HABP4蛋白具體與哪種或者哪些蛋白之間相關(guān)性更強的研究較少。參照文獻[37]的方法，本研究利用RT-qPCR技術(shù)檢測可能與HABP4蛋白存在相互作用的10個蛋白在16只羊胰腺中的mRNA相對表達水平并進行GO 功能富集。結(jié)果顯示，HABP4基因和各個基因之間均具有相關(guān)性，其中主要與UAB52、RPL32、RPS9 3個基因具有顯著相關(guān)性。UBA52是由N端泛素和C端核糖體蛋白L40組成的融合蛋白[38]。在細(xì)胞內(nèi)，被泛素化修飾的UBA52將泛素轉(zhuǎn)移到靶蛋白，使靶蛋白被26S蛋白酶體識別而降解，從而有效調(diào)節(jié)細(xì)胞內(nèi)蛋白質(zhì)的周轉(zhuǎn)，參與細(xì)胞內(nèi)熱休克反應(yīng)、DNA損傷修復(fù)、細(xì)胞周期調(diào)控等一系列細(xì)胞生物學(xué)反應(yīng)[39-40]。RPL以及RPS家族作為核糖體的主要成分，在調(diào)控蛋白質(zhì)合成過程中起著重要作用[41]。同時，RPL以及RPS家族還可以參與調(diào)節(jié)DNA損傷修復(fù)、凋亡、細(xì)胞增殖和分化、耐藥性、細(xì)胞遷移和侵襲[42]。研究報道顯示，細(xì)胞內(nèi)RPL32被敲低后會導(dǎo)致核糖體應(yīng)激，而核糖體是細(xì)胞內(nèi)翻譯的重要細(xì)胞器[43]，直接影響rRNA的成熟[44]。同時，RPS9蛋白與腫瘤的發(fā)生存在一定的關(guān)聯(lián)[45-46]。因而，根據(jù)HABP4蛋白和UAB52、RPL32、RPS9共3個蛋白相對表達的正相關(guān)，推測HABP4蛋白可能在細(xì)胞內(nèi)蛋白質(zhì)的翻譯、合成以及細(xì)胞周期、凋亡等方面發(fā)揮一定的調(diào)控作用，具體作用有待進一步研究。GO功能富集分析常被用于檢測差異基因在GO term的富集程度，能一定程度上對基因的功能進行限定和描述[47]。2007年，Uenishi等[48]利用GO功能富集分析了豬源HABP4基因可能存在的分子功能并繪制了樹形結(jié)構(gòu)圖。本研究通過GO功能富集發(fā)現(xiàn)，HABP4基因的生物學(xué)過程分析主要富集于細(xì)胞定位、新陳代謝、多生物過程以及生物過程負(fù)調(diào)節(jié)；分子功能方面與翻譯起始、核糖體轉(zhuǎn)錄、細(xì)胞質(zhì)中的代謝有關(guān)。

4 結(jié) 論

山羊HABP4基因cDNA全長1 496 bp，包含5′UTR 61 bp，CDS區(qū)1 254 bp，3′UTR 181 bp，編碼417個氨基酸，在胰腺中表達水平最高。HABP4與UAB52表達量呈極顯著正相關(guān)，與RPL32、RPS9呈顯著正相關(guān)，推測其在細(xì)胞定位、新陳代謝、多生物過程及生物過程負(fù)調(diào)節(jié)以及翻譯起始、核糖體轉(zhuǎn)錄、細(xì)胞質(zhì)中的代謝等方面發(fā)揮著重要作用，為進一步探討HABP4在調(diào)控細(xì)胞凋亡過程中的作用提供了數(shù)據(jù)參考。