大黃魚選育家系F1代生長性狀研究

章 霞，油九菊，柳敏海，傅榮兵，李偉業，殷小龍，徐志進，張 川

( 浙江省舟山市水產研究所，浙江 舟山 316000 )

大黃魚選育家系F1代生長性狀研究

章 霞，油九菊，柳敏海，傅榮兵，李偉業，殷小龍，徐志進，張 川

( 浙江省舟山市水產研究所，浙江 舟山 316000 )

通過選取舟山藍科海洋生物研究所(岱山)大黃魚和舟山市水產研究所大黃魚親本，構建12個家系：岱山♀×舟山所♂(A1、A2、A3)，舟山所♀×岱山♂(B1、B2、B3)，岱山♀×岱山♂(C1、C2、C3)，舟山所♀×舟山所♂(D1、D2、D3)，并對其3月齡、6月齡、9月齡、12月齡的體長(L)、體質量(m)進行測量分析。試驗結果證明家系A1、A2、A3、B1、B2、B3、C1、C2、C3的體長體質量在各個月齡均大于家系D1、D2、D3，且差異顯著(P<0.05)，并擬合其中兩個生長特性差異較大的家系生長曲線如：A1，m=0.0212L2.9248(r2=0.9987)，D3，m=0.0362L2.6816(r2=0.9971)，擬合度較好，表明該生長階段網箱養殖的大黃魚為等速生長，但不同家系間體長與體質量的關系差異顯著。綜合結果表明，通過岱山大黃魚和舟山所大黃魚產生F1代大黃魚的生長性能優于同代舟山所大黃魚子代，可應用于下一步的育種工作。

大黃魚；體長；體質量

大黃魚(Pseudosciaenacrocea)，屬石首魚科、黃魚屬，曾是我國海洋四大主捕對象之首，為中下層暖溫性海水魚類[1]。到了20世紀80年代，大黃魚因過度捕撈而資源匱乏[2]，由此大黃魚養殖業的快速發展應運而生。但由于世代的近交繁殖導致大黃魚生長速度慢，抗病力差，體色發暗等種質問題瀕發[3-6]，因此大黃魚良種選育以及種質改良工作的開展具有重要實際意義。

在水產育種中，家系育種在新品種培育、品種優化等生產實踐中有著重大的實踐意義。目前國內研究者已對部分水產經濟動物進行較大規模的家系構建及生長性能研究，并取得較好的選育效果，如半滑舌鰨(Cynoglossussemiliaevis)[7]、石斑魚(Epinephelus)[8]、大菱鲆(Scophthalmusmaximus)[9]、虹鱒(Oncorhynchusmykiss)[10]等均已有諸多報道，大黃魚[11-12]也有部分報道。

本研究擬構建較大數量大黃魚全同胞及半同胞家系，利用舟山藍科海洋生物研究所(岱山)大黃魚和舟山市水產研究所大黃魚構建12個家系，分別對這12組家系進行3月齡、6月齡、9月齡、12月齡的生長性能測定并比較分析，旨在獲得大黃魚家系不同生長發育階段生長性狀表現情況，為進一步的育種工作提供參考資料。

1 材料與方法

1.1 親魚來源

2013年12月，取舟山市水產研究所養殖的成熟個體(個體大、性成熟、活力強)300尾作為親本備用，其中雌魚220尾，體長約26.4～32.6 cm，體質量354～545 g，雄魚80尾，體長約25.5～32.4 cm，體質量252～504 g。另從岱山的舟山藍科海洋生物研究所購進大黃魚親魚300尾，其中雌魚120尾，體長約32.1～35.5 cm，體質量501～650 g，雄魚180尾，體長約28.5～32.4 cm，體質量350～513 g，作為親本備用。

1.2 子代繁育

親魚經暫養、營養強化后，于2014年3月采用人工授精方法構建家系，共構建12個家系，岱山♀×舟山所♂(A1，A2，A3)，舟山所♀×岱山♂(B1，B2，B3)，岱山♀×岱山♂(C1，C2，C3)作為試驗組，舟山所♀×舟山所♂(D1，D2，D3)作為對照組。

將12個家系相同質量的受精卵分別放在3 m3的玻璃鋼桶內進行魚苗孵化和前期養殖，培育溫度(22±0.5) ℃。40 d時，倒桶，隨機留3000尾大黃魚進行養殖，到3月齡時進行倒桶，并每個家系挑選出生長較快的魚苗900尾各自放置于3 m3的玻璃鋼桶中進行養殖，各個玻璃鋼桶聯通。6月齡時挑選出生長較快的魚苗450尾，9月齡時挑選出生長較快的魚苗350尾繼續放置在玻璃鋼桶中進行養殖，每日7:00與16:00投喂相同質量的天邦大黃魚配合飼料，飼料投喂量一般控制缸內魚體總質量的2%～3%，兼以其不再攝食為標準。若出現病害狀況，減少投喂或者停食，待病情好轉，繼續投喂。6—11月份進行流水養殖，11月到次年3月進行換水養殖，每日刷桶2次，換水量為200%，每半個月倒桶1次。

1.3 大黃魚生物學測定

每3個月每個群體隨機抽取30個樣本，用直尺(精確至0.1 cm)測量大黃魚的體長，用電子天平(精確至0.1 g)測量體質量，求取平均值，結果以平均值±標準差表示。

1.4 數量參數分析

擬合各群體體長(L)與體質量(m)冪函數曲線關系，其關系式為：m=aLb[13]

式中a為生長的條件因子，b為冪指數函數。

其他參數計算公式如下：

增積量=(m2-m1)(L2-L1)/(t2-t1)2

肥滿度/%=m/L3×100%

絕對質量增加率/g·d-1=(m2-m1)/(t2-t1)

絕對增長率/cm·d-1=(L2-L1)/(t2-t1)

變異系數=SD/MN[11]

式中，m為體質量(g)，L為體長(cm)，m1、m2和L1、L2分別為時間t1、t2時的體質量(g)和體長(cm)；SD為標準差，MN為平均值。

數據統計分析用獨立樣本T檢驗不同樣本間體長、體質量相關系數差異，采用統計軟件Excel 2007 和SPSS 17.0進行單因素方差分析并采用Duncan氏法進行顯著性差異檢驗。

2 結果與分析

2.1 體長參數對比

12個大黃魚家系各時期平均體長和變異系數見表1～2，家系D1、D2、D3的體長在各個月齡均小于家系A1、A2、A3、B1、B2、B3、C1、C2、C3，差異顯著(P<0.05)。3月齡時，家系A2的平均體長約為D1的1.45倍，家系A1、A2、A3、C1、C2、C3顯著大于家系B1、B2、B3(P<0.05)；6月齡時，家系 C2體長生長顯著優于家系A3、B1、D1、D2、D3；9月齡時，家系A1、A2、A3體長生長優于其他家系，顯著優于家系C1、C2、C3和家系D1、D2、D3；12月齡時家系A1、A2、A3和家系B1、B2、B3體長大于家系C1、C2、C3和家系D1、D2、D3，且差異顯著。

表1 12個家系不同時期的平均體長 cm

注：同列上標小寫字母不同代表顯著差異(P<0.05)，下同.

表2 12個家系不同時期平均體長的變異系數

2.2 體質量參數對比

12個大黃魚家系各時期平均體質量見表3，變異系數見表4，家系D1、D2、D3的平均體質量在各個月齡均小于家系A1、A2、A3、B1、B2、B3、C1、C2、C3，差異顯著(P<0.05)，家系A2的平均體質量約為D1的2.79倍。在3月齡時，家系A1、A2、A3、C1、C2、C3顯著大于家系B1、B2、B3，差異顯著(P<0.05)。這可能由魚苗的孵化率、仔魚期的卵黃囊體積及吸收率、仔魚期的成活率、食餌相遇密度等相關因素的差異造成[14-15]。6月齡時，家系A1、A2、A3、B1、B2、B3、C1、C2、C3顯著優于家系D1、D2、D3(P<0.05)，9月齡、12月齡時，家系A1、A2、A3、B1、B2、B3體質量優于家系C1、C2、C3、D1、D2、D3，且差異顯著(P<0.05)。12個家系的體長、體質量變異系數分別為0.023～0.127，0.089～0.391，體質量變異量比體長大。

表3 12個家系不同時期的平均體質量 g

表4 12個家系不同時期平均體質量的變異系數

2.3 肥滿度及增積量

12個大黃魚家系各時期肥滿度見表5，各個時期各個大黃魚家系的肥滿度差異較大，且按照養殖時間延長而有所下降，家系A3、B2、B3、D2、D3肥滿度在3月齡優于其他家系，差異顯著(P<0.05)，6月齡時，肥滿度家系A3>家系A1、A2、D1、D3>D2>家系B1、B2、B3，家系C1>家系C2、C3，且差異顯著(P<0.05)。9月齡時，各家系肥滿度均優于C2，且差異顯著(P<0.05)，12月齡時，各家系之間肥滿度差異不顯著。

12個大黃魚家系各時期增積量見圖1，在3～6月齡時，家系A1、A2、A3、B1、B2、B3、C1、C2、C3，顯著大于家系D1、D2、D3(P<0.05)，其中家系B2、B3顯著大于其他家系。在6～9月齡時，家系A2顯著高于其他家系，家系A1、A3顯著高于家系B1、B2、B3(P<0.05)，家系B1、B2、B3顯著高于家系C1、C2、C3(P<0.05)，家系D1、D2、D3、C2顯著高于家系C1、C3(P<0.05)。在9～12月齡時，日增積量家系B1>家系A2、B2、B3>家系A1、D1、D3>家系A3、C1、C2、C3、D2，差異顯著(P<0.05)。在3～12月齡養殖期間，日增積量家系B1、B2、B3>家系A1、A2、A3>家系C1、C2、C3、D1、D2、D3，差異不顯著。

表5 12個家系不同時期的肥滿度

2.4 絕對質量增加率及絕對增長率

12個大黃魚家系各時間段的絕對質量增加率見圖2，隨養殖時間增加而減少，養殖前期3～6月齡，家系A1、A2、A3、B1、B2、B3、C1、C2、C3日質量增加率約達0.4 g/d，約是家系D1、D2、D3的2倍，家系A1、A2、A3、B1、B2、B3、C1、C2、C3絕對質量增加率顯著大于家系D1、D2、D3(P<0.05)；6～9月齡養殖期間，絕對質量增加率家系A1、A3>家系C1、C2、C3、D1、D2、D3，且差異顯著(P<0.05)；在9～12月齡時期，絕對質量增加率家系由大到小依次為家系A2、A1、B1、D2、C2、C3、D3、C1、B2、B3、A3、D1，差異不顯著。在3～12月齡養殖期間，絕對質量增加率家系A1、A2、A3、B1、B2、B3>家系C1、C2、C3、D1、D2、D3，其中家系A1、A2、A3、B1、B2、B3絕對質量增加率顯著大于家系D1(P<0.05)。

各個大黃魚家系在不同時期的絕對增長率見圖3。隨養殖時間增加而減少，養殖前期3～6月齡絕對增長率家系B1、B2、B3>家系C1、C2、C3>家系A1、A2、A3>家系D1、D2、D3，為0.066～0.088 cm/d；6～9月齡養殖期間，絕對增長率家系A1>家系A2>家系A1、B1、D1、D2、D3>家系B2、B3家系C1、C2、C3，差異顯著(P<0.05)；9～12月齡養殖期間，絕對增長率家系A1、D3>家系A2、A3、B1、B2、B3、C2、D2>家系C1、C3、D1，差異顯著(P<0.05)；3～12月齡，各家系之間絕對增長率家系B1、B2、B3>家系A1、A2、A3、D2、D3>家系C1、C2、C3、D1，差異不顯著。

2.5 生長型比較

以Keys氏公式m=aLb擬合3～12月齡2個生長性能差異較大的大黃魚選育家系體質量(m，g)與體長(L，cm)的關系，回歸方程分別：

A1：m=0.0212L2.9248(r2= 0.9987)

D3：m=0.0362L2.6816(r2=0.9971)

這兩個家系均表現為b≈3，即該生長階段網箱養殖的大黃魚為等速生長，但不同家系間體長和體質量的關系有差異(圖4)。

圖1 12個家系各個時期的增積量

圖2 12個家系各個時期的絕對質量增加率

圖3 12個群體各個時期的絕對增長率

圖4 2個家系的體長與體質量的關系

3 討 論

3.1 大黃魚選育家系F1代生長性狀結果

本研究表明，12個家系在各個時期的體質量變異系數超過體長，因此在大黃魚生長性能的篩選上，選擇體質量作為直接指標，體長作為間接指標，這與許益銨[16]的研究一致。經擬合兩個生長性能差異最大的家系的生長曲線關系，符合硬骨魚類體長體質量之間的冪函數關系：m=aLb，b一般為2.5～4.0[17-18]，且各回歸方程中b值均接近于3，r2值較高，擬合度較好，這表明這兩個家系在養殖過程中可排除環境條件和營養狀況的影響而導致的生長差異，呈勻速生長狀態，符合李明云等[19-23]的研究。本試驗中日增積量、絕對增長率和絕對質量增加率等數據顯示家系A1、A2、A3、B1、B2、B3的生長性狀優于家系C1、C2、C3、D1、D2、D3，肥滿度處于中等，此試驗結果符合李明云等[19,24-25]的研究結果，造成此現象的原因可能是因為岱山大黃魚親本為2010年捕撈野生大黃魚所產生F1代，遺傳多樣性豐富，而舟山所大黃魚親本是舟山所在2000年捕撈的野生大黃魚經多代近親交配所得的子代，可能存在種質退化的因素，從而造成岱山大黃魚和舟山所大黃魚所產的子代在生長性能上優于舟山所大黃魚所產的子代，但具體的遺傳機制和相關因素仍需要通過不同的方法比如分子標記、蛋白鑒定、酶聯免疫等手段來進行探索和驗證。

3.2 基于大黃魚快速生長性狀進行家系選育的可行性

家系選育在水產動物優良性狀選育上應用廣泛，韋信鍵等[26]利用不平衡槽設計構建大黃魚生長快速家系選育，獲得快速生長家系3個，生長較快家系11個；陳松林等[27]通過查明半滑舌鰨不同家系雌、雄魚遺傳和生理性別比例及其與生長速率的關系，篩選出生理雌魚比例高且生長快的家系；劉旭東等[28]通過褐牙鲆(Paralichthysolivaceus)家系選育出第二代群體的平均體質量和平均全長比未經選育的對照組提高2.91 g、1.33 cm，差異顯著，且具有高遺傳力；唐章生等[29]對吉富羅非魚(GIFTOreochromisniloticus)家系進行生長快速選育，獲得生長速度比群體的均值高18.00%以上的家系等；本研究結果表明家系A1、A2、A3、B1、B2、B3的生長速度均快于家系C1、C2、C3、D1、D2、D3，且差異較大，可繼續應用于后期的良種選育。由此可見，基于大黃魚快速生長性狀進行家系選育可行，并且鑒于家系選育的遺傳背景較為統一，后期可再采用基于電子標記或分子標記的手段輔助家系選育，提高對大黃魚生長性狀遺傳選擇的準確性，并在育種進程中將具有性狀優勢的等位基因或基因型進行富集，應用于QTL定位、遺傳連鎖圖譜構建和全基因組關聯分析[30-32]等進一步提高育種家系的生長性能。

[1] 蘇永全, 張彩蘭, 王軍, 等. 大黃魚養殖[M]. 北京:海洋出版社, 2004:1-10.

[2] 劉家富, 劉招坤. 福建閩東大黃魚Larimichthyscrocea(Richardson)產業展望[J]. 現代漁業信息, 2008, 23(12):3-5.

[3] 黃振遠, 蘇永全, 張建設, 等. 閩粵群和岱衢群養殖大黃魚(Pseudosciaenacrocea)及其雜交子代遺傳差異的SSR分析[J]. 海洋與湖沼, 2011, 42(4):592-595.

[4] 王軍, 全成干, 蘇永全, 等. 官井洋大黃魚遺傳多樣性的RAPD分析[J]. 海洋學報, 2001, 23(3):87-91.

[5] 王娟, 封永輝, 蔡立勝, 等. 來自大黃魚(Pseudosciaenacrocea)腸道的弧菌拮抗菌的篩選與鑒定[J]. 海洋與湖沼, 2010, 41(5):707-713.

[6] 劉賢德, 韋信鍵, 蔡明夷, 等. 大黃魚22 個微衛星標記在F1家系中的分離方式及與生長性狀的相關分析[J]. 水產學報, 2012, 36(9):1322-1330.

[7] 陳松林, 杜民, 楊景峰. 半滑舌鰨家系建立及其生長和抗病性能測定[J]. 水產學報，2010, 34(12):1669-1704.

[8] 劉付永忠, 趙會宏, 劉曉春, 等. 赤點石斑魚與斜帶石斑魚雜交初步研究[J]．中山大學學報:自然科學版, 2007, 46(3):72-75.

[9] 王新安, 馬愛軍, 雷霽霖, 等. 大菱鲆不同家系生長性能的比較[J]. 海洋科學, 2011, 35(4):1-8.

[10] 盧國, 谷偉, 白慶利. 電子標記輔助虹鱒家系建立及快速生長家系篩選[J]. 中國水產科學, 2012, 19(1):77-83.

[11] 劉賢德, 蔡明夷, 王志勇, 等. 閩—粵東族大黃魚生長性狀的相關與通徑分析[J]. 中國海洋大學學報：自然科學版, 2008, 38(6):916-920.

[12] 劉招坤. 大黃魚群體選育世代F1和F2生長比較分析[J]. 水產科學, 2015, 34(1):38-42.

[13] 葉金清. 官井洋大黃魚的資源和生物學特征[D]. 上海:上海海洋大學, 2012.

[14] 殷名稱. 魚類仔魚期的攝食和生長[J]. 水產學報, 1995,19(4):335-342.

[15] 王宏田, 張培軍. 環境因子對海產魚類受精卵及早期仔魚發育的影響[J]. 海洋科學, 1998,22(4):50-52.

[16] 許益銨. 舟山附近海域大黃魚遺傳多樣性及家系的生長性狀研究[D]. 舟山:浙江海洋學院,2014.

[17] 華元渝, 胡傳林. 魚種重量與體長相關公式(W=aLb)的生物學意義及其應用[G]∥中國魚類學會.魚類學論文集，北京:科學出版社, 1981:125-131.

[18] Brown M E. Experimental studies on growth. in the Physiology of Fishes [M]. London:Academic Press, 1957:361-400.

[19] 李明云, 胡玉珍, 苗亮, 等. 岱衢洋和官井洋大黃魚自交與雜交子代生長性能及雜交優勢分析[J]. 2010, 34(6):679-684.

[20] 黃偉卿, 張藝, 柯巧珍, 等. 大黃魚選育子二代生長性狀研究[J]. 南方水產科學, 2013, 9(3):14-19.

[21] 陳慧, 陳武, 林國文, 等. 官井洋種群網箱養殖大黃魚的形態特征與生長式型[J]. 海洋漁業, 2007, 29(4):331-336.

[22] 陳成進. 人工養殖大黃魚主要生長特征觀察[J]. 現代漁業信息, 2011, 26(3):24-29.

[23] 徐恭昭, 羅秉征, 黃頌芳. 大黃魚生殖季節體長體質量關系的種內變異[J]. 海洋科學集刊, 1984, 22(1):1-8.

[24] 鄭煒強, 黃偉卿, 韓坤煌, 等. 大黃魚選育群體與野生群體雜交F1生長性狀研究[J]. 水產科學, 2014, 33(11):667-672.

[25] 黃偉卿, 韓坤煌, 陳仕璽, 等. 海捕野生大黃魚選育子代生長性能及現實遺傳力分析[J].水產科學, 2016, 35(3):204-209.

[26] 韋信健, 劉賢德, 王志勇. 32個大黃魚家系早期階段生長形狀比較及遺傳參數估計[J]. 集美大學學報:自然科學版, 2013, 18(5):321-328.

[27] 陳松林, 李仰真, 張靜, 等. 半滑舌鰨快速生長及高雌性家系的篩選[J]. 水產學報, 2013, 347(4):481-488.

[28] 劉旭東, 劉志鵬, 王亞楠, 等. 對牙鲆進行一代選擇之后的育種效果分析[J]. 海洋科學進展, 2012, 30(4):548-555.

[29] 唐章生, 林勇, 黎筠, 等. 吉富羅非魚不同家系的生長性狀差異[J]. 廣西師范大學學報:自然科學版, 2011, 29(3):74-79.

[30] 陳宏. 基因工程實驗技術[M]. 北京:中國農業出版社, 2005:19-36.

[31] 張玲. 微衛星DNA 標記研究進展及應用[J]. 安徽農業科學, 2007, 35(4):972-975.

[32] 竺俊全, 錢偉平. 海水養殖種類種質退化原因及對策[J]. 寧波大學學報:理工版, 2000, 13(2):87-91.

GrowthCharacteristicsofF1CrossinLargeYellowCroakerPseudosciaenacrocea

ZHANG Xia, YOU Jiujü, LIU Minhai, FU Rongbing, LI Weiye, YIN Xiaolong, XU Zhijin, ZHANG Chuan

( Zhoushan Fisheries Research Institute of Zhejiang Province，Zhoushan 316000，China )

The body length and body weight of 3-12 month old were compared among large yellow croakerPseudosciaenacroceaF1 crossed between Zhoushan Blue Institute of (Daishan) and Zhoushan fisheries research institute self-reproducing large yellow craker including Daishan♀×Zhoushan Research Institute♂ (Family A1, A2, A3); Zhoushan Research Institute♀×Daishan♂ (Family B1, B2, B3); Daishan♀×Daishan♂ (Family C1, C2, C3); Zhoushan Research Institute♀×Zhoushan Research Institute♂ (Family D1, D2, D3), total 12 families. The results showed that the fish in families A1, A2, A3, B1, B2, B3, C1, C2, and C3 had larger body length and body weight than the fish in families D1, D2, and D3 at every period (P<0.05). The fitting relationship between body length (L) and body weight (m) was expressed as the following:Family A1,m=0.0212L2.9248(r2= 0.9974), Family D3,m=0.0362L2.6816(r2=0.9952), indicating that the growth was the same in a net cage culture. However, the relationship between body length and body weight was different among different families. The findings suggest that F1 crossed between Daishan and Zhoushan Research Institute large yellow croaker be applied to the breeding of large yellow croaker in the future.

Pseudosciaenacrocea; body length; body weight

10.16378/j.cnki.1003-1111.2017.06.007

S965.323

A

1003-1111(2017)06-0735-06

2016-10-27；

2017-01-22.

浙江省公益技術研究農業項目(2015C32111)；浙江省重大科技專項(2016C02055-7-3)；2016浙江省海洋與漁業綜合管理與發展專項資金資助項目(舟財農[2016]33號)；浙江省海洋與漁業綜合管理和產業發展專項資金資助項目(舟財農[2017]40號);舟山市海洋經濟創新發展示范工作專項資金資助項目(國海科字[2016]496號).

章霞(1989-)，女，工程師；研究方向：水生動物增養殖. E-mail：yufan414515@163.com.通訊作者：柳敏海(1979-)，男，高級工程師；研究方向：水產動物增養殖.E-mail：Okso1125@aliyun.com.