李偉業(yè),羅海忠,殷小龍,張川,柳敏海,章霞,邱豪軍,毛志增,油九菊,徐志進*

(1.浙江省舟山市水產研究所，浙江 舟山 316000；2.舟山市海洋與漁業(yè)局，浙江 舟山 316000；3.浙江萬里學院，浙江 寧波 315100)

大黃魚Larimichthyscrocea為暖水性魚類，主要分布于黃海南部、東海，以及瓊州海峽以東的南海北部沿岸，是中國重要的經濟魚類之一。大黃魚的水溫適應范圍為10～32 ℃，最適溫度為18～25 ℃，低于14 ℃和高于30 ℃攝食下降，低于7 ℃時大黃魚死亡[1]。近年來，由于生態(tài)環(huán)境的破壞，極端天氣頻繁發(fā)生，中國南方沿海地區(qū)也常受寒流侵襲，致使養(yǎng)殖海區(qū)水溫長期偏低，對于不耐低溫的大黃魚造成了嚴重的損害。據報道，2016年南方漁業(yè)因寒潮造成經濟損失超過20億元，僅大黃魚損失就超過3億元[2]。因此，開展低溫對大黃魚影響的研究，對開發(fā)抗寒大黃魚新品種、促進大黃魚產業(yè)可持續(xù)發(fā)展具有重要意義。

目前，有關低溫對大黃魚影響的研究較多，但多集中于生長特性、生理生化指標等方面，而對大黃魚耐寒分子機制研究較少。金屬硫蛋白(Metallothionein, MT)是由微生物和動植物產生的金屬結合蛋白，是一類可被多種因素誘導的內源性蛋白，在水生動物體內各生長發(fā)育階段發(fā)揮著重要作用，具有較強的清除自由基能力[3]。水通道蛋白(Aquaporin，AQP)，又名水孔蛋白，是一種位于細胞膜上的蛋白質(內在膜蛋白)，不僅在水分子的快速跨膜轉運中起主要作用，還具有陽離子通道的作用，同時有助于CO2和NH3的交換運輸[4]。目前，已發(fā)現的水通道蛋白有200多種亞型, 其中對AQP1的研究最為廣泛和深入[5]。AQP1是由Agre等從紅細胞膜分離純化Rh血型多肽時偶然發(fā)現的，相對分子質量為 28 000。它不僅可以調節(jié)滲透壓, 還與小分子氣體交換和運輸等有關, 也與某些疾病有著重要的相關性[6]。熱休克蛋白(Heat shock protein，HSPs)，是一組在生物進化過程中高度保守的蛋白分子，HSP60 家族是HSP家族中重要的成員之一，可以有效保護機體細胞免受逆境脅迫[7]。TCP1 (t-com-plex polypeptide-1)屬于HSP60家族Ⅱ類亞族，是一種廣泛存在于細胞漿中的異型寡聚蛋白，對提高細胞耐受應激能力、維持細胞內環(huán)境的穩(wěn)定具有重要作用[8]。

研究發(fā)現，MT、AQP1、TCP1 3個基因與魚類的抗寒性能相關[3-8]。為此，本研究中選取了上述基因，并比較分析了耐低溫品系與非耐低溫品系大黃魚中3個基因在低溫條件下的表達特點，以驗證3個候選基因與大黃魚耐寒能力的相關性，旨在從分子層面闡釋大黃魚在低溫脅迫下基因調控的相關機理，以期為今后研究大黃魚的耐寒分子機制和抗寒育種提供科學參考。

1 材料與方法

1.1 材料

試驗用魚取自浙江省舟山市水產研究所繁育的F5代耐低溫品系大黃魚和非耐低溫品系大黃魚幼魚，耐低溫品系大黃魚體質量為(55.35±3.52)g，體長為(14.72±0.25)cm，非耐低溫大黃魚體質量為(56.45±4.41)g，體長為(14.91±0.32)cm。將試驗魚置于圓形玻璃鋼桶(直徑為200 cm、高為100 cm)內，納米氣管增氧，于每天9:00、16:00各投喂一次配合飼料(天邦大黃魚配合飼料3號料)，并于投喂后1 h吸污，日換水率100%。試驗用水為過濾后的自然海水。

1.2 方法

1.2.1 低溫脅迫對比試驗 試驗設2個組，分別為耐低溫品系大黃魚組(N組)和非耐低溫品系大黃魚組(F組)，每個組設2個平行。每個試驗桶各放入50尾大黃魚幼魚。試驗桶水溫為16 ℃(自然溫度)。試驗前，將試驗魚移至試驗桶中，暫養(yǎng)7 d，待試驗魚的狀態(tài)穩(wěn)定后開始正式試驗。試驗起始溫度為16 ℃，按照1 ℃/d的速率降溫，分別在16、13、10、7 ℃時取樣，并觀察試驗魚的行為狀態(tài)。在各溫度下隨機抽取3尾魚，取其肝臟、肌肉，并于液氮中速凍后置于-80 ℃超低溫冰箱中保存?zhèn)溆谩?/p>

1.2.2 熒光定量PCR 采用天根生化科技(北京)有限公司試劑盒提取總RNA，經瓊脂糖電泳檢測RNA 的完整性，用RNA作為模板，以Oligo dT18作為引物合成第一鏈cDNA。根據GenBank數據庫中MT、AQP1、TCP1基因序列，使用Primer 5.0軟件設計3對特異性引物(表1)，用于熒光定量分析(qPCR)，選取β-actin作為內參基因[9]。擴增反應在ABI 7300plus熒光定量擴增儀上進行，熒光定量檢測采用SYBR定量檢測試劑盒(TaKaRa)。反應體系(共20 μL)包括：SYBR Premix Ex TaqⅡ10 μL，上、下游引物各0.8 μL，Rox Reference Dye 0.4 μL，cDNA模板2 μL, ddH2O 6 μL。反應程序為：95 ℃下預變性30 s; 95 ℃下變性 5 s，60 ℃ 下退火31 s, 共進行40個循環(huán)；最后在72 ℃下延伸10 min。溶解反應：95 ℃ 下反應15 s，60 ℃ 下反應60 s (數據采集)，95 ℃ 下反應15 s。采用2-ΔΔCT法計算基因的相對表達量。

表1 大黃魚基因定量分析所用引物

1.3 數據處理

試驗結果均以平均值±標準差(mean±S.D.)表示，采用SPSS 19.0軟件對數據進行單因素方差分析，采用Duncan法進行組間多重比較，顯著性水平設為0.05，極顯著性水平設為0.01。

2 結果與分析

2.1 不同溫度下耐低溫和非耐低溫品系大黃魚的行為狀態(tài)

從表2可見：N組大黃魚和F組大黃魚在16、13、10 ℃下，行為狀態(tài)較為相似；在7 ℃下，N組大黃魚和F組大黃魚的行為狀態(tài)出現差異，N組大黃魚盡管活力一般，也不攝食，但未出現翻肚側游的現象，在水體各層分布較均勻，而F組大黃魚活力不佳，部分出現翻肚側游的現象，且多分布于水體的中上層。

表2 不同溫度下耐低溫和非耐低溫品系大黃魚的行為狀態(tài)

2.2 不同溫度下耐低溫和非耐低溫品系大黃魚肝臟和肌肉中MT基因表達量的變化

從圖1可見：隨著溫度的降低，N、F組大黃魚肝臟和肌肉中MT基因表達量均呈上升趨勢，且各試驗組均有顯著性差異(16 ℃時除外)(P<0.05)；在7 ℃時，大黃魚肝臟和肌肉中MT基因表達量均達到最高，N組、F組大黃魚肝臟中的表達量分別為16 ℃時的10.8倍和22.3倍，N組、F組大黃魚肌肉中的表達量分別為16 ℃時的9.8倍和18.3倍，且F組大黃魚表達量也顯著高于N組(P<0.05)。

標有不同字母者表示組間有顯著性差異(P<0.05)，標有相同字母者表示組間無顯著性差異(P>0.05)，下同。The means with different letters are significantly different in the groups at the 0.05 probability level, and the means with the same letter are not significant differences, et sequentia.圖1 不同溫度下耐低溫和非耐低溫品系大黃魚肝臟和肌肉中MT基因相對表達量的變化Fig.1 Relative expression levels of MT gene in liver and muscle of large yellow croaker Larimichthys crocea exposed to different temperature

2.3 不同溫度下耐低溫和非耐低溫品系大黃魚肝臟和肌肉中AQP1基因表達量的變化

從圖2(a)可見：N、F組大黃魚肝臟中AQP1表達量均隨溫度的降低呈先升高后降低的趨勢，且各試驗組間均有顯著性差異(P<0.05)，F組表達量均顯著高于N組(P<0.05)；在13 ℃時N組、F組大黃魚的表達量均達到最高，上調幅度分別為16 ℃時的1.73倍和1.82倍；當溫度降到10 ℃時，各試驗組大黃魚肝臟中的AQP1表達量顯著下降；在7 ℃時N組、F組大黃魚表達量達到最低，分別為16 ℃時的0.23%和0.39%。

圖2 不同溫度下耐低溫和非耐低溫品系大黃魚肝臟和肌肉中AQP1基因相對表達量的變化Fig.2 Relative expression levels of AQP1 gene in liver and muscle of large yellow croaker Larimichthys crocea exposed to different temperature

從圖2(b)可見：N、F組大黃魚肌肉中AQP1表達量也均隨著溫度的降低呈先升高后降低的趨勢；在13 ℃時N組、F組大黃魚的表達量均達到最高，上調幅度分別為16 ℃時的1.49倍和1.64倍，且F組表達量顯著高于N組(P<0.05)；當溫度降到10 ℃時，各試驗組大黃魚肌肉中的AQP1表達量顯著下降；在7 ℃時N組、F組大黃魚表達量達到最低，分別為16 ℃時的0.14%和0.26%，但N組和F組表達量間無顯著性差異(P>0.05)。

2.4 不同溫度下耐低溫和非耐低溫品系大黃魚肝臟和肌肉中TCP1基因表達量的變化

從圖3可見：隨著溫度的降低，N、F組大黃魚肝臟和肌肉中TCP1基因表達量均呈上升趨勢，且各試驗組間有顯著性差異(16 ℃時除外)(P<0.05)；在7 ℃時，肝臟和肌肉中表達量均達到最高，N組、F組大黃魚肝臟中TCP1基因表達量分別為16 ℃時的28.3倍和19.3倍，N組、F組大黃魚肌肉中TCP1基因表達量分別為16 ℃時的9.3倍和6.3倍，且N組表達量顯著高于F組(P<0.05)。

圖3 不同溫度下耐低溫和非耐低溫品系大黃魚肝臟和肌肉中TCP1基因相對表達量的變化Fig.3 Relative expression levels of TCP1 gene in liver and muscle of large yellow croaker Larimichthys crocea exposed to different temperature

3 討論

3.1 低溫脅迫對耐低溫和非耐低溫品系大黃魚MT基因表達的影響

MT為一類可被多種因素誘導的內源性蛋白，在水生動物體內各生長發(fā)育階段均發(fā)揮著重要作用，具有較強的清除自由基能力，其清除羥基自由基的能力是SOD的1萬倍，清除氧自由基的能力是GSH的25倍。眾多研究表明，MT的合成與溫度有明顯的溫度效應關系。Oh等[10]研究發(fā)現，冷、熱應激和耐力鍛煉等均可引起大鼠肝臟MT含量升高；Baer等[11]也發(fā)現，MT含量與小鼠對寒冷的耐受性呈正相關(r=0.78)。另有研究認為，MT可能在冷應激中起到穩(wěn)定細胞膜、保護細胞器的重要作用[12-14]。

本研究中發(fā)現，隨著溫度降低，大黃魚肝臟和肌肉中MT基因表達量均呈顯著上升，且非耐低溫品系大黃魚的表達量顯著高于耐低溫品系，這表明低溫脅迫會促進MT基因的表達，這與Sardella等[15]的研究結果相似。此外，這或許從側面反映了耐低溫品系大黃魚比非耐低溫品系大黃魚具有更強的低溫適應能力。

3.2 低溫脅迫對耐低溫和非耐低溫品系大黃魚AQP1基因表達的影響

AQP基因不僅在水分子的快速跨膜轉運中起到主要作用，還具有陽離子通道的作用，同時有助于CO2和NH3的交換運輸，它屬于水選擇性跨膜通道[16]。研究表明，AQP1基因與動物機體的抗寒能力關系密切，在樹蛙Hylachrysoscelis和面包酵母Saccharomycescerevisiae研究中均發(fā)現，AQP1 基因與這些物種的抗寒性能相關[17-18]。此外，朱華平等[19]研究表明，AQP1基因表達量與溫度顯著相關，且隨著溫度的降低，AQP1基因的表達量逐漸下調。本研究中發(fā)現，AQP1基因表達量與溫度也具有顯著的相關性，溫度降至10 ℃后，耐低溫品系和非低溫品系大黃魚肝臟和肌肉的AQP1基因表達均顯著下調，在7 ℃時表達量最低，且耐低溫品系大黃魚肝臟和肌肉的AQP1表達量均低于非耐低溫品系。結合耐低溫品系和非耐低溫品系大黃魚在7 ℃下的行為狀態(tài)(耐低溫品系大黃魚緩慢游動，非耐寒品系大黃魚部分側游)，筆者認為，耐低溫品系大黃魚具有更強的耐低溫能力，這也可能與其AQP1表達量下調幅度更大有關，這也與朱華平等[16]的研究結果相似。其機制可能是大黃魚機體為適應低溫條件，AQP1 基因表達量大幅下調，降低細胞膜上AQP1 的含量或活性，減少水分子的跨膜運輸，使細胞內環(huán)境保持相對穩(wěn)定，降低外界低溫對細胞的損害。

3.3 低溫脅迫對耐低溫和非耐低溫品系大黃魚TCP1基因表達的影響

TCP1伴侶蛋白是一大類在生物大分子折疊、組裝及降解過程中起重要協(xié)同作用且與細胞應激損傷關系密切的蛋白。對提高細胞耐受應激的能力、維持細胞內環(huán)境的穩(wěn)定具有重要作用[8,20]。眾多研究表明，TCP1與生物體抗寒能力相關，如酵母SaccharomycescerevisiaeTCP-1-alpha 和TCP-1-beta 的低溫誘導表達模式顯示，TCP-1是一種冷激表達蛋白[21]，肉蠅Sarcophagacrassipalpis滯育期的低溫存活能力與熱休克蛋白如TCP-1、HSP70家族等的上調表達關系密切[22]；蔥蠅Deliaantique幼蟲耐寒能力增強也與其TCP-1的上調表達密切相關[23]。在水生生物中也有研究表明，TCP基因與抗寒密切相關，如殷勤等[24]在低溫脅迫試驗中證實，TCP基因與凡納濱對蝦耐寒性狀存在相關性，不同基因型的個體在耐寒能力上具有差異；謝建麗等[25]對尼羅羅非魚的研究中發(fā)現，隨著低溫脅迫程度的加強，TCP基因表達上調幅度增大，這表明TCP-1-beta和TCP-1-eta是潛在的耐寒相關基因。本研究中發(fā)現，TCP1基因表達量與溫度也具有顯著的相關性，隨著溫度的降低，各試驗組大黃魚肝臟和肌肉中的TCP1基因表達量均顯著上調，在溫度降至7 ℃后，TCP1基因表達量上調幅度最大。此外，耐低溫品系大黃魚肝臟和肌肉中的TCP1表達量均顯著高于非耐低溫品系。本研究結果與謝建麗等[25]的研究結果相似，這表明TCP1與大黃魚的抗寒能力密切相關，同時，耐低溫品系大黃魚的TCP1表達量顯著高于非耐低溫品系，這也從側面反映了耐低溫品系大黃魚比非耐低溫品系具有更強的抗寒能力。

4 結論

1)2個品系的大黃魚MT、AQP1、TCP1 3個基因在低溫下表現出顯著差異性，表明這3個基因與大黃魚抗寒能力具有一定的相關性。

2)耐低溫品系大黃魚在低溫下存活率和活力要優(yōu)于非耐低溫品系，結合MT、AQP1、TCP1 3個與抗寒相關的基因的表達差異，證明兩者在低溫適應能力和抗寒能力上存在一定區(qū)別。

3)在大黃魚耐低溫品系的篩選過程中，尋找出抗寒相關因子及其基因，獲得抗寒相關基因的分子標記，是低溫選育耐低溫品系的重要手段。而找到與性狀顯著相關的SNP位點是進行分子標記輔助育種的前提和基礎，因此，今后應重點進行抗寒基因編碼序列SNP位點的篩查。