王玲軍, 陶妍, 黨靜, 王孫勇

(上海海洋大學食品學院,上海201306)

在魚類肌肉蛋白質的組成中,肌球蛋白約占肌原纖維蛋白的50%以上,它主要決定了魚肌肉蛋白質的加工適性。肌球蛋白分子含有兩個相對分子質量約200 000的重鏈亞基和4個相對分子質量約20 000的輕鏈亞基,每個重鏈的氨基端形成一個球狀的頭部(S1)和頸狀結構(S2),余下的羧基端形成具α-螺旋結構的被稱之為輕酶解肌球蛋白(light meromyosin,LMM)的尾部;兩對輕鏈分別非共價地結合在重鏈的頸部[1]。肌球蛋白分子的頭部具有ATPase的活性;而LMM具有使肌球蛋白在生理條件下形成絲狀的能力,它的結構性質在很大程度上決定了魚肌肉蛋白質的凝膠形成特性[2]。

自從1990年 Gerlach等[3]首次報道了鯉魚(Cy prinus carpio)肌球蛋白重鏈基因的多樣性后,Kikuchi等[4]進一步證明了鯉魚肌球蛋白重鏈屬于“多基因家族”。3種與馴化溫度相關的肌球蛋白重鏈同工型基因從溫度馴化的鯉魚骨骼肌中被分離到,并且這些同工型基因伴隨馴化溫度的表達模式也被確定[5];近年來,3種草魚(Ctenopharyngodon idellus)肌球蛋白重鏈同工型基因也被分離,它們之間在LMM的一級結構上展現(xiàn)了明顯的差異[6-7]。

綜觀世界淡水漁業(yè)的發(fā)展,低成本、高產量的低值淡水魚占淡水魚總產量的50%以上[8]。鰱魚(Hy pophthalmichthys molitrix)屬我國淡水魚資源中具有代表性的魚種之一,具有低值、高產的特點,因此也是加工利用的主要對象。為了探明鰱魚肌肉蛋白質的加工適性,對來自不同季節(jié)的鰱魚骨骼肌肌球蛋白Ca2+-A TPase的熱穩(wěn)定性和凝膠形成特性進行了測定和比較,發(fā)現(xiàn)春、夏季鰱魚肌球蛋白的熱穩(wěn)定性明顯優(yōu)于秋、冬季鰱魚的[9,10],并且它們在凝膠形成模式上也不同[11]。然而,迄今為止,國內外關于鰱魚肌球蛋白分子水平上的研究報道很少。據此,作者以夏季和冬季鰱魚為研究對象,對編碼其骨骼肌LMM的部分基因進行cDNA克隆,以期分離不同的肌球蛋白重鏈同工型基因,并比較不同同工型之間在LMM的一級結構氨基酸組成上的差異,初步探明鰱魚肌球蛋白的功能性質發(fā)生季節(jié)性變化的分子機理。

1 材料與方法

1.1 實驗材料

鰱魚(體質量730 g,體長40 cm),取自上海市青浦區(qū)養(yǎng)殖場,于2005年8月和2006年1月購自上海市圖們路農貿市場,分別作為夏季和冬季的實驗材料。鮮活鰱魚運至實驗室后,用鈍器擊斃,立即取其背部白色骨骼肌切成1 cm×1 cm×1 cm見方小塊,用液氮快速凍結后,保存于-80℃冰箱待用。

1.2 總RNA提取和第一股cDNA合成

取-80℃保藏的鰱魚肌肉約0.4 g于50 mL離心管中,立即加入 Trizol試劑(100 mg∶1 mL)(Invitrogen Corp.,Carlsbad,CA,USA),用高速組織搗碎機處理2～3 min;取1 mL搗碎液于1.5 mL離心管中,28℃下靜置5 min,加入0.2 mL氯仿,來回振蕩15 s,28℃下放置2～3 min;采用高速冷凍離心機離心15 min(4℃,11 900g);上清液轉移至新的離心管中,加入0.5 mL異丙醇,28℃下靜置10 min,離心10 min(條件同上);沉淀加75%酒精1 mL,離心15 min(4 ℃,≤7 500g);棄去上清液,沉淀在鋁箔上自然干燥 30 min后,溶解于10μL DEPC·H2O中,-80 ℃保存。

第一股cDNA的合成使用Super RT Kit(Invitrogen Corp)。2μL總 RNA、5μL AP通用引物(Adapter Primer)、2μL ddH2O混勻,70 ℃保溫10 min,冰浴 1～2 min;依次加入 4μL 5×First strand buffer、4μL dNTP mixture(2.5 mmol/L)、2μL DTT(0.1 mol/L),42 ℃保溫 2 min,再加入1μL Reverse Transcriptase(RNase H-)(200 U/L)保溫70 min;70℃保溫15 min、冰上冷卻后得到第一股cDNA。

1.3 輕酶解肌球蛋白基因的PCR擴增

分別以冬季和夏季鰱魚的第一股cDNA為模板,通過3′-RACE法對編碼鰱魚肌球蛋白重鏈羧基端的部分輕酶解肌球蛋白基因進行擴增,見圖1。正向引物為MYH-12F (5′-CTGCACTCTCAGAACACAAGCCT-3’),該引物根據鯉魚肌球蛋白重鏈同工型基因的保守序列設計;反向引物為AUAP通用引物(Abridged Universal Amplification Primer)(5′-GGCCACGCGTCGACTA GTAC-3′)。PCR反應體系如下:5μL第一股cDNA、各5 μL 10μmol/L的正向和反向引物、2.5μLEx TaqDNA polymerase(5 U/μL)(Takara,Otsu,Japan)、10μL 10 ×ExTagBuffer、16μL dN TP Mixture,用無菌水將反應液調至100μL。PCR反應條件如下:94℃預變性4 min、94℃變性30 s、59℃退火1 min、72℃延伸1 min、最終72℃終延伸5 min,共進行30個循環(huán)。

圖1 鰱魚部分輕酶解肌球蛋白基因的PCR擴增策略圖Fig.1 Schematic representation of PCR amplification for partial light meromyosin from silver carp

1.4 序列分析

通過PCR擴增得到的目的片段經DNA純化試劑盒(天根生物科技公司,北京)純化后,與p GMT質粒載體(天根生物科技公司,北京)按1∶1在16℃下用 T4DNA連接酶連接過夜,然后轉化感受態(tài)細胞 Top10;通過藍白篩選后,采用 PCR法進行插入檢驗,每個平板挑選8～10個陽性克隆,由上海生工生物工程有限公司進行DNA測序。采用Chromas Version 1.62和BioEdit Version 7.0.5對測序結果進行分析。

1.5 氨基酸的親水性分析和分子系統(tǒng)樹的構建

通過BioEdit Version 7.0.5軟件,根據 Hopp等[12]的方法對所得鰱魚兩種肌球蛋白重鏈同工型基因的氨基酸序列進行親水性分析。分子系統(tǒng)樹根據 Saitou等[13]提出的鄰接法(Neighbor-joining method)構建。作者將鰱魚兩種同工型的部分LMM區(qū)域的氨基酸序列與其他鯉科魚類和溫血動物的相應序列相比較,用CLUSTAL W進行多重比對后,通過 MEGA version 3.1[14]軟件對獲得的結果進行聚類分析,系統(tǒng)樹中各結點的置信度由自引導值(Bootstrap value)估計,重復次數為1 000。用于構建分子系統(tǒng)樹的序列除表1中涉及的3種鯉科魚類的各種肌球蛋白重鏈同工型之外,其它的序列來自于:Cad-雞的成熟骨骼肌(U87231)[15]、Rperi-圍產期大鼠的肌肉 (X04267)[16]、Rα-大鼠α心臟肌 (X15938)[17]、Rβ-大 鼠β心臟肌(X15939)[18]、Argo-扇貝A rgopecten irradians橫紋肌(X55714)[19]、Placo-扇貝Placopecten magellanicus橫紋肌 (U59294)[20]、Cgiz-雞的胃組織(X06546)[21]、Cnon-雞的非肌肉組織(M26510)[22]。

2 結果與討論

2.1 鰱魚肌球蛋白重鏈同工型基因的分離及部分輕酶解肌球蛋白基因的核苷酸序列比較

以基于冬季和夏季鰱魚的第一股cDNA為模板,均擴增得到了大約750 bp的DNA片段;將這兩個片段純化后插入p GM-T質粒載體,分別隨機挑選8～10個陽性克隆用于測序。結果表明,基于冬季鰱魚的所有克隆均顯示了相同的DNA序列,而基于夏季鰱魚的所有克隆亦顯示了相同的DNA序列。但是,這兩種克隆之間在DNA序列上顯示了明顯的差異,它們之間的核苷酸同源性僅79.5%(表1、圖2)。前人的研究工作已經證明了哺乳動物肌球蛋白重鏈同工型基因之間在3′-和5′-區(qū)域會顯示大的序列差異[23];Imai等[5]和 Tao等[6]分別從鯉魚和草魚骨骼肌中分離到的3種肌球蛋白重鏈同工型基因之間在終止密碼子的上游和下游區(qū)域顯示了較低的核苷酸同源性。作者從鰱魚骨骼肌中分離到的兩種cDNA克隆之間在終止密碼子的上游和下游區(qū)域亦顯示了較大的核苷酸序列差異,據此,可以認為這兩種cDNA克隆代表了兩個不同的肌球蛋白重鏈同工型基因,它們分別被命名為低溫型(sc-w)和高溫型(sc-s)。Yuan等[24]曾在蛋白質水平上報道了與季節(jié)溫度變化相關的兩種鰱魚肌球蛋白同工型,他們發(fā)現(xiàn)兩種同工型之間在蛋白質的熱穩(wěn)定性、凝膠形成特性等方面顯示了明顯的差異。將鰱魚兩種cDNA克隆與草魚和鯉魚的相應核苷酸序列進行比對,發(fā)現(xiàn)鰱魚低溫型(scw)與草魚 10℃型(gc10)之間的同源性高達96.1%、鰱魚高溫型(sc-s)與草魚中間(gcI)和30℃型(gc30)之間的同源性則較高;另一方面,鰱魚高溫型(sc-s)較其低溫型(sc-w)與鯉魚3種同工型之間顯示了更高的同源性。

表1 鰱魚、草魚和鯉魚輕酶解肌球蛋白cDNA克隆的核苷酸序列及推斷的氨基酸序列的同源性Tab.1 DNA nucleotide and deduced amino acid sequences identities between light meromyosin cDNA clones isolated from silver carp and those of grass carp and common carp

圖2 鰱魚、草魚和鯉魚部分輕酶解肌球蛋白基因的核苷酸序列比較Fig.2 The comparison in the nucleotide sequences of partial light meromyosin genes from silver carp and those of grass carp and common carp

2.2 鰱魚部分輕酶解肌球蛋白的氨基酸組成分析

根據編碼鰱魚兩種肌球蛋白重鏈同工型的部分輕酶解肌球蛋白的cDNA序列推斷其氨基酸序列,并進行比對,同時與草魚和鯉魚的相應氨基酸酸序列進行比較。結果顯示,sc-w與sc-s的氨基酸同源性僅為88.1%,見表1;sc-w與gc10之間顯示了97.2%的最高同源性。而sc-s與gcI和gc30之間顯示較高的同源性,分別為96.8%和94.9%,與cc30顯示94.9%的同源性。根據Tao等[6]的報道,gcI也是從30℃馴化的草魚骨骼肌中分離到的。另外,sc-w、gc10和ccI的氨基酸序列較其它同工型的氨基酸序列均縮短了一個氨基酸殘基,見圖3。根據Imai等[5]的報道,ccI不同于gcI,它主要是從10℃馴化的鯉魚骨骼肌中分離到的。結果似乎表明了來自于相似溫度環(huán)境的淡水魚在肌球蛋白重鏈的結構上更具有相似性。然而,sc-w與鯉魚的cc10之間反而顯示了較之與ccI和cc30之間略低的氨基酸同源性,此結果亦出現(xiàn)在草魚3種與鯉魚3種全長LMM同工型的氨基酸比對中[7],可能與cc10的基因歧化有關,此推測被以下的分子系統(tǒng)樹進一步證明。

由圖3的陰影部分可以發(fā)現(xiàn),sc-w和gc10在His-21、Ser-24、Thr-29、Leu-49、Ser-62、Asn-87、Asp-140、Ala-147、Ser-167、Gly-176、Tyr-177、Lys-180、Leu-181、Glu-192、Glu-217 的氨基酸殘基位點上顯示了與其它同工型不同的特有變異,并且其中His-21、Ser-24、Leu-49、Ser-62、Asn-87、Asp-140、Ser-167、Gly-176、Tyr-177、Lys-180、Leu-181 顯示了高度的保守。根據 Hopp等的方法[12],對兩種鰱魚同工型的部分輕酶解肌球蛋白的氨基酸序列進行親水性分析,結果見圖4。兩種同工型的氨基酸基本都顯示了親水性模式。然而,這兩種同工型的氨基酸的親水周期還是有些明顯的不同,比如sc-w相對于 sc-s在第 31-43、55-58、68-70、79-87、141-145和189-195的氨基酸區(qū)域展示了更高的親水性。

圖3 鰱魚、草魚和鯉魚部分輕酶解肌球蛋白的氨基酸序列比較Fig.3 The comparison in the deduced amino acid sequences of partial light meromyosins from silver carp and those of grass carp and common carp

圖4 鰱魚兩種肌球蛋白重鏈同工型基因部分輕酶解肌球蛋白氨基酸序列的親水性分析Fig.4 Hydrophilicity analysis of partial deducedamino acid sequences for light meromyosin from the two types of silver carp myosin heavy chain isoform

我們已經在蛋白質水平上通過差示掃描量熱儀(DSC)和流變儀對冬季和夏季鰱魚的肌球蛋白熱穩(wěn)定性和凝膠形成特性進行了測定,發(fā)現(xiàn)前者的熱穩(wěn)定性明顯低于后者,且兩者之間在凝膠形成模式上顯示了明顯的差異[11];Tao等[25-26]報道了10℃馴化和冬季草魚的肌球蛋白熱穩(wěn)定性較30℃馴化和夏季草魚的明顯要低,彼此之間亦顯示了不同的凝膠形成模式。據此,可以認為不同肌球蛋白重鰱同工型在輕酶解肌球蛋白特定區(qū)域氨基酸殘基上的變異可能導致它們之間在結構上的差異,進而引起功能性質上的變化。這樣,似乎基于低溫馴化的淡水魚表達的肌球蛋白重鏈同工型必須在它們的初級結構上發(fā)生某些氨基酸殘基的變異才能有助于其適應低溫的棲息環(huán)境。

2.3 分子系統(tǒng)樹分析

我們根據3種鯉科魚類和兩種溫血動物的骨骼肌、心臟肌或非肌肉組織的部分LMM氨基酸序列構建了分子系統(tǒng)樹。由圖5所示,系統(tǒng)樹分為3個主要的分支:一個最大的分支由草魚、鰱魚和鯉魚的骨骼肌、溫血動物雞的成年骨骼肌以及大鼠的心臟肌的肌球蛋白重鏈同工型組成,其自引導值為99%;另兩個分支分別由非脊椎動物的扇貝肌和雞的非肌肉組織的肌球蛋白重鏈同工型組成,它們的自引導值均為100%。在大分支中,3種鯉科魚類的8個肌球蛋白重鏈同工型形成了一個大組,其中,gcI、gc30、sc-s和cc30均從高溫棲息的魚肌肉中分離,因此,它們之間顯示了更近的基因距離;而ccI、sc-w、gc10均從低溫棲息的魚肌肉中分離,它們之間則顯示了近的基因距離。雖然cc10也是來自于低溫馴化的鯉魚,但在這個大組中,它形成了獨立于其它類型的分支,其自引導值為72%。來自于雞的非肌肉組織的肌球蛋白重鏈同工型與所有其它來自于骨骼肌的肌球蛋白重鏈同工型之間顯示了低的同源性。分子系統(tǒng)樹的分析結果進一步證明了廣溫性淡水魚肌球蛋白重鏈的基因歧化和功能進化在很大程度上受棲息環(huán)境溫度的影響。這些結果完全符合了前人報道的對鯉魚和草魚所作的相關分析結果[7,27]。

圖5 根據鰱魚和其它生物的輕酶解肌球蛋白氨基酸序列構建的分子系統(tǒng)樹Fig.5 Phylogenetic analysis based on the amino acid sequences of partial light meromyosins from silver carp and those of other species

3 結 語

作者通過RT-PCR法和3′-RACE法從冬季和夏季鰱魚的骨骼肌中分別分離到一種低溫型(sc-w)和一種高溫型(sc-s)肌球蛋白重鏈同工型基因,兩種同工型基因的cDNA序列包含了部分3′-翻譯區(qū)和全部3′-非翻譯區(qū),編碼了鰱魚輕酶解肌球蛋白(LMM)羧基端的一部分氨基酸序列。兩種同工型之間的核苷酸或氨基酸序列之間顯示了較低的同源性;與草魚和鯉魚的各種同工型之間的比對結果,基本上顯示了來自低溫和高溫棲息環(huán)境的魚其各自肌球蛋白重鏈同工型之間具有更高的同源性,表明了廣溫性的淡水魚必須在它們的肌球蛋白重鏈的初級結構上發(fā)生某些氨基酸殘基的變異才能有助于其適應低溫的棲息環(huán)境。分子系統(tǒng)樹的分析結果進一步證明了淡水魚肌球蛋白重鏈同工型的基因距離與它們所棲息的環(huán)境溫度有關。

[1]Harrington W F,Rodgers M E.Myosin[J].Annu Rev Biochem,1984,53:35-73.

[2]Sano T,Noguchi S F,Matsumoto J J,et al.Thermal gelation characteristics of myosin subfragments[J].J Food Sci,1990,55(1):55-58.

[3]Gerlach G,Turay L,Malik T A,et al.Mechanism of temperature acclimation in the carp:a molecular approach[J].Am J Physiol,1990,259:237-244.

[4]Kikuchi K,Muramatsu M,Hirayama Y,et al.Characterization of the carp myosin heavy chain multigene family[J].G ene,1999,228,189-196.

[5]Imai J I,Hirayama Y,Kikuchi K,et al.cDNA cloning of myosin heavy chain isoforms from and their gene expression associated with temperature acclimation[J].J Exp Biol,1997,200:27-34.

[6]Tao Yan,Kobayashi M,Liang C S,et al.Temperature-dependent expression patterns of grass carp fast skeletal myosin heavy chain genes[J].Comp Biochem Physiol,2004,B139:649-656.

[7]Wang Sun Yong,Tao Yan,Liang C S,et al.cDNA cloning and characterization of temperature-acclimation-associated light meromyosins from grass carp fast skeletal muscle[J].Comp Biochem Physiol,2008,B149:378-387.

[8]張慜,張駿.國內外低值淡水魚加工與下腳料利用的研究進展[J].食品與生物技術學報,2006,25(5):115-120.ZHANG Min,ZHANG Jun.A research review of low value freshwater fishes processing and their discards utilization[J].Journal of Food Science and Biotechnology,2006,25(5):115-120.(in Chinese)

[9]韓黎,陶妍,汪之和.淡水魚肌球蛋白Ca2+-ATPase熱穩(wěn)定性的季節(jié)變化[J].上海水產大學學報,2007,16(4):362-367.HAN Li,TAO Yan,WANG Zhi-he.Seasonal changes in thermal stability of Ca2+-ATPase for myosins from freshwater fish[J].Journal of Shanghai Fisheries University,2007,16(4):362-367.(in Chinese)

[10]韓黎,陶妍,汪之和.草魚和鰱魚肌球蛋白Ca2+-ATPase熱力學活化參數的季節(jié)變化[J].水生生物學報,2008,32(2):22-27.HAN Li,TAO Yan,WANG Zhi-he.Seasonal changes in thermodynamic activation parameters of myosin Ca2+-ATPase from grass carp and silver carp[J].Acta Hydrobiologica Sinica,2008,32(2):22-27.(in Chinese)

[11]董永,陶妍.淡水魚骨骼肌肌球蛋白熱穩(wěn)定性及凝膠特性的季節(jié)變化[J].食品工業(yè)科技,2008,29(7):124-127.DONG Yong,TAO Yan.Seasonal changes in thermal stability and gel-forming property of myosin from grass carp fast skeletal muscle[J].Science and Technology of Food Industry,2008,29(7):124-127.(in Chinese)

[12]Hopp T P,Woods K R.Prediction of protein antigenic determinants from amino acid sequences[J].Proc Natl Acad Sci USA,1981,78:3824-3828.

[13]Saitou N,Nei M.The neighbor-joining method:a new method for reconstructing phylogenetic trees[J].Mol Biol Evol,1987,4:406-425.

[14]Kumar S,Tamura K,Nei M.MEGA3:Integrated software for molecular evolutionary genetics analysis and sequence alignment[J].Briefings in Bioinformatics,2004,5:150-163.

[15]Chao T H,Bandman E.Cloning,nucleotide sequence and characterization of a full-length cDNA encoding the myosin heavy chain from adult chicken pectoralis major muscle[J].G ene,1997,199(1-2):265-270.

[16]Karsch-Mizrachi I,Feghali R,Shows T B,et al.Generation of a full-length human perinatal myosin heavy-chain-encoding cDNA[J].G ene,1990,89:289-294.

[17]Mcnally E M,Gianola K M,Leinwand L A.Complete nucleotide sequence of full length cDNA for rat alpha cardiac myosin heavy chain[J].Nucleic Acids Res,1989,17:7527-7528.

[18]Kraft R,Bravo-Zehnder M,Taylor D A,et al.Complete nucleotide sequence of full length cDNA for rat beta cardiac myosin heavy chain[J].Nucleic Acids Res,1989,17:7529-7530.

[19]Nyitray L,Goodwin EB,Szent-Gyorgyi A G.Nucleotide sequence of full length cDNA for a scallop striated muscle myosin heavy chain[J].Nucleic Acids Res,1990,18:7158.

[20]Perreault-Micale C,Kalabokis V,Nyitray L,et al.Sequence variations in the surface loop near the nucleotide binding site modulate the ATP turnover rates of molluscan myosins[J].J Muscle Res Cell Motil,1996,17:543-553.

[21]Yanakisawa M,Hamada Y,Katsuragawa Y,et al.Complete primary structure of vertebrate smooth muscle myosin heavy chain deduced from its complementary DNA sequence.Implications on topography and function of myosin[J].J Mol Biol,1987,198:143-157.

[22]Shohet R V,Conti M A,Kawamoto S,et al.Cloning of the cDNA encoding the myosin heavy chain of a vertebrate cellular myosin[J].Proc Natl Acad Sci USA,1989,86:7726-7730.

[23]Weiss A.The mammalian myosin heavy chain gene family[J].Annu Rev Cell Dev Biol,1996,12:417-439.

[24]Yuan C,Kaneniwa M,Wang X,et al.Seasonal expression of two types of myosin with different thermostability in silver carp muscle(Hypophthalmichthys molitrix)[J].J Food Sci,2005,70:C326-C331.

[25]Tao Yan,Kobayashi M,Fukushima H,et al.Changes in enzymatic and structural properties of grass carp fast skeletal myosin induced by the laboratory-conditioned thermal acclimation and seasonal acclimatization[J].Fish Sci,2005,71:195-204.

[26]Tao Yan,Kobayashi M,Fukushima H,et al.Changes in rheological properties of grass carp fast skeletal myosin induced by seasonal acclimatization[J].Fish Sci,2007,73:189-196.

[27]Nihei Y,Kobiyama A,Ikeda D,et al.Molecular cloning and mRNA expression analysis of carp embryonic,slow and cardiac myosin heavy chain isoforms[J].J Exp Biol,2006,209:188-198.