牙鲆選育F3代家系的建立及遺傳效應分析

田永勝，陳紅林，2，齊文山，2，姜靜，2，張英平，2，劉萬軍，3，李寶玉，劉壽堂，孫德強，陳松林*

(1.中國水產科學研究院 黃海水產研究所 農業部海洋漁業可持續發展重點實驗室，山東 青島 266071；2.上海海洋大學 水產與生命學院，上海 201306；3.中國海洋大學 海洋生命學院，山東 青島 266003； 4.海陽市黃海水產有限公司，山東 海陽 265100)

牙鲆選育F3代家系的建立及遺傳效應分析

田永勝1，陳紅林1，2，齊文山1，2，姜靜1，2，張英平1，2，劉萬軍1，3，李寶玉4，劉壽堂4，孫德強4，陳松林1*

(1.中國水產科學研究院 黃海水產研究所 農業部海洋漁業可持續發展重點實驗室，山東 青島 266071；2.上海海洋大學 水產與生命學院，上海 201306；3.中國海洋大學 海洋生命學院，山東 青島 266003； 4.海陽市黃海水產有限公司，山東 海陽 265100)

牙鲆是太平洋西岸海水養殖、捕撈和自然增殖的主要魚類，但是人工養殖病害多、生長慢，針對這一問題，從2003年開始開展了牙鲆選擇育種。本文利用篩選和培育的F1、F2代牙鲆家系核心群體，構建了F3代家系52個，在生長到196 d時，對家系1 560尾魚苗全長、體寬、體質量、全長日增量和體質量日增量5個生長性狀進行測量和計算，獲得7 800個數據。采用“加性-顯性”線性模型、最小范數二階無偏估計法(MINQUE)對以上生長性狀的邊際平均值、加性方差分量、顯性方差分量、隨機方差分量、狹義遺傳力、廣義遺傳力、親本加性效應、雜交隨機效應等遺傳參數進行了評估。結果顯示：1202、1206、1226和1211家系以上生長性狀的邊際平均值同時顯著高于其他家系。以上5個生長性狀的加性方差分量比率為(0.099±0.01)～(0.24±0.02)，顯性方差分量比率為(0.30±0.02)～(0.31±0.02)，隨機方差分量比率為(0.48±0.02)～(0.61±0.01)，5個性狀狹義遺傳力為0.099～0.24，廣義遺傳力為0.38～0.54，以上遺傳參數均達極顯著性水平(p<0.01)。綜合比較18個親本5個性狀的加性效應預測值，發現F0917、F0908、F0905、F0915、F09125和KS(韓國牙鲆)親本在5個性狀上都具有極顯著的正向效應(p<0.01)，其他親本5個性狀都不同程度的表現為顯著(p<0.05)或極顯著(p<0.01)的負向加性效應。對42個雜交組合顯性隨機效應進行預測，D(2×6)、D(2×18)、D(3×4)、D(7×12)、D(7×14)和D(8×13) 6個雜交組合在5個性狀上都表現為極顯著(p<0.05)正效應，利用這些雜交組合建立家系將促進牙鲆以上5個性狀的生長。本文研究結果，為牙鲆優良家系的建立及新品種的培育初步篩選出了性狀優良的親本、雜交組合和家系，同時為牙鲆進一步選擇育種提供了理論和技術依據。

牙鲆；F3代家系；遺傳力；加性效應；顯性效應

1 引言

牙鲆(Paralichthysolivaceus)是中國、韓國、日本等太平洋西岸國家的主要經濟魚類，2009年世界牙鲆產量達到4.5萬噸，牙鲆的養殖產量26 477.5 t，占鲆鰈類產量的30%，在鲆鰈魚類養殖中處于第2位，僅次于大菱鲆(57 042 t)[1]，在漁業生產中具有重要的地位。我國牙鲆的大量養殖始于20世紀60年代，但直到21世紀初在國內從未有過關于牙鲆育種方面的研究報道。牙鲆繁殖和育苗親魚主要來源于養殖群體，經過多代的人工繁殖后，牙鲆的規模化養殖出現了種質退化現象，如生長速度慢、抗逆力下降、病害多(白化病、腸道白濁、鰻弧菌病、淋巴囊腫病、纖毛蟲感染等)[2]，導致牙鲆養殖效益降低，養殖者對牙鲆養殖信心不足，養殖規模相對縮小。同時也導致了牙鲆養殖中使用抗病藥物的種類和用量大增，使牙鲆商品魚品質無法保障；另一方面過多的使用藥物使養殖環境受到污染，造成惡性循環，因此培育生長快、抗逆力強的牙鲆養殖新品種已成為產業發展的迫切需求。

從21世紀初開始，國內在牙鲆養殖群體遺傳多樣性分析[3—4]，MHC抗病分子標記篩選及多態性分析[5—6]，遺傳連鎖圖譜構建[7]，數量遺傳參數評估[8—9]，雌核發育[10—11]，家系建立及選擇育種[12—13]等多方面開展了研究，利用家系選擇和雜交育種方法，初步培育出一個生長快牙鲆新品種“鲆優1號”[14]。日本通過抗淋巴囊腫相關微衛星標記輔助，篩選出抗淋巴囊腫病養殖新品種[15]。到目前國內已經建立了牙鲆育種核心種群及技術體系，為牙鲆養殖新品種的培育打下了基礎。

本文主要利用已經建立的牙鲆F1和F2代家系選育群體為基礎，通過回交、雜交和雌核發育等方式建立了牙鲆F3代家系，對其生長性狀進行測定，并對其數量遺傳參數進行分析，擬繼續篩選優良的牙鲆親本和家系，為培育出下一代生長快、抗病力強的牙鲆新品種提供育種群體和理論依據。

2 材料和方法

2.1 牙鲆優良繁殖群體構建

2007年利用實驗室建立的3個基礎群體：牙鲆抗鰻弧菌感染群體(RS)、日本群體(JS)和黃海野生群體(YS)，通過簡單多系交配方式建立牙鲆F1代63個家系，并通過遺傳參數估計等一系列方法，選育和篩選出生長快、抗逆性強的優良家系F0719、F0750、F0751、F0768和F0780家系[2，9]。2008年從韓國濟州島引進韓國牙鲆群體(KS)。2009年利用篩選出的F1代家系、韓國牙鲆、日本牙鲆群體，通過回交和雜交、雌核發育方式，構建F2代家系，經過3年的人工培育，生長和抗病力性狀測定、多重比較和遺傳參數估計，篩選到生長快、成活高的優良家系，如：F0905、F0908、F0915、F0917、F0921和F0927(F0750雌核發育系)、F0939、F0990、F09104、F09119、F09121和F09125等12個家系。2012年春天利用篩選到的F1、F2家系和韓國牙鲆群體親本，采用雜交、回交及雌核發育方式建立了F3代家系52個。

2.2 親魚培育和生殖調控

以上F1(253尾)、F2(443尾)代家系和韓國牙鲆群體(99尾)親魚共計795尾，分別飼養在6個30 m2的水泥池中，循環水培育，水流量4～5 m3/h，并利用純氧增氧系統向池水中加入純氧，溶氧量為7～8 mg/L。平時投喂人工加工的冷凍顆粒混合餌料(成分為4%小麥粉、1%魚粉，95%小雜魚，少量維生素C和E)，每天投喂2次，投喂量為親魚體質量的2%～3%。在進入繁殖季節后，以投喂野生雜魚為主(沙丁魚、玉筋魚等)，每天的投喂量為親魚體質量的5%～6%，并同時在飼料中添加500 mg/kg維生素E。

2012年3月初對牙鲆親魚進行生殖調控，工廠化培育親魚生殖調控必須利用人工控光、控溫、水流調控和餌料調控相結合的方式，才能達到成熟和產卵的目的。利用鍋爐加熱方式將親魚池水溫逐漸升溫到14℃，保持水溫至產卵，繁殖期最適宜水溫14～17℃。光照采用長日照模式，白天/黑夜時間比為16 h/8 h，光照強度300～500 lx。經過2個月的生殖調控牙鲆親魚達到性成熟。

2.3 家系的構建及標準化培育

在牙鲆繁殖季節，以不完全雙列雜交方法，利用以上18個牙鲆家系親本，構建牙鲆F3代全同胞和半同胞家系52個。親魚性成熟后，人工采卵，干法受精，受精后將卵置于50 cm×50 cm×100 cm的孵化網箱中，在14～17℃的水中進行孵化24 h左右，收取上浮卵，轉移到實驗車間，利用容積為3 m3的玻璃缸進行人工孵化，每缸初始放卵量為20 mL(20 000粒卵左右)。孵化后隨著魚苗成長，隨時做好苗種的標準化工作，在完成變態后每個家系魚苗數量保持在5 000尾左右，魚苗達5 cm時調整到2 000尾，10 cm時調整到1 000尾，盡量保證每個缸的家系魚數量保持一致。在魚苗長到15 cm左右時(日齡150 d)，利用熒光染料進行標記，每個家系標記150～250尾魚苗，集中放入容積為30 m3的水池中培育。

2.4 家系生長性狀數據測定

所有家系魚苗生長到196 d左右時，進行家系生長性狀進行測量。測量指標主要包括魚苗全長、體寬、體質量3個指標，每個家系測量30尾魚。共測定1 560尾魚苗的生長性狀，再分別利用體長/日齡、體質量/日齡計算出全長日增量和體質量日增量。共獲得數據7 800個。

2.5 數據分析方法

使用QGAStation軟件中的“加性-顯性”模型對F3代家系全長、體寬、體質量、全長日增量、體質量日增量5個性狀進行邊際均值及遺傳參數估計，其線性模型如下：

(1)

式中，y為所有觀察值構成的向量，常數向量1是系數為常數1的向量；μ是群體平均數；eA是加性效應向量；UA是加性效應向量關聯矩陣；eD是顯性效應向量；UD是顯性效應向量關聯矩陣；eε是隨機效應向量[16]。

采用最小范數二階無偏估計法[17]——MINQUE(1)法(Minimum Norm Quadratic Unbiased Estimation by setting 1 for all prior values)，估算各性狀的方差分量：即表型方差、加性方差、顯性方差和機誤方差。

3 結果

3.1 家系生長性狀的邊際均值

對牙鲆家系全長、體寬、體質量和日增量的估計邊際均值進行估算(表1)，根據全長邊際均值選出生長最快的前10個家系為 1202、1226、1211、1209、1206、1221、1219、1225、1218、1208。根據體寬估計邊際均值選出生長最好的前10個家系為1202、1226、1209、1211、1206、1223、1221、1225、1219、1207。根據體質量估計邊際均值選出生長最好的前10個家系為1202、1226、1211、1209、1206、1221、1219、1208、1203、1223，根據體質量日增量的估計邊際均值選出生長最好的前10個家系為1202、1226、1211、1209、1206、1221、1219、1208、1203、1223。綜合比較4個性狀所選擇的家系中，1202、1206、1226和1211等家系在4個性狀中均表現優秀。

表1 全長、體寬、體質量和日增體質量的估計邊際均值

續表1

注： 模型中出現的協變量在下列值處進行評估: Number of days=195.060 7。

3.2 生長性狀的方差分量和遺傳力

使用A-D模型對牙鲆F3代家系的全長、體寬、體質量、日增體長和日增體質量5個性狀進行遺傳分析，結果表明，全長、體寬、體質量、日增體長和日增體質量的加性、顯性方差分量、方差分量比率均達到極顯著水平(p<0.01)。以上5個性狀的加性方差分量比率為(0.099±0.01)～(0.24±0.02)，顯性方差分量比率為(0.30±0.02)～(0.31±0.02)，隨機方差分量比率為(0.48±0.02)～(0.61±0.01)。顯性方差分量均大于加性方差分量，機誤方差分量比率在0.48到0.61之間，達到了極顯著水平(p<0.01)。以上性狀狹義遺傳力在(0.099±0.01)到(0.24±0.02)之間，5個性狀都顯示低遺傳力，廣義遺傳力在(0.39±0.01)到(0.52±0.02)之間，均達到極顯著水平(p<0.01)，其中日增體質量的狹義、廣義遺傳力都是最低(表2)。

表2 利用“A-D”模型對牙鲆F3代家系5個性狀方差分量值、方差分量比率以及遺傳力估計結果(M±S.E)

注：** 表示差異極顯著(p<0.01)；*表示差異顯著(p<0.05)。

3.3 親本加性效應

表3列出了參與牙鲆F3代家系建立的18個牙鲆親本群體5個生長性狀的加性效應估計值，全長加性效應在(17.15±1.822)到(-12.86±1.352)之間，其中F0917、F0908、F0905、F0921、F09125、F0915、F0927和KS群體親本具有顯著或極顯著正向加效應，其他親本則為顯著的負向效應。體寬的加性效應在(5.37±0.687)到(-5.68±0.553)之間，同樣以上8個親本具有正向效應，其他家系為負向效應。體質量的加性效應在(19.62±2.427)到(-11.41±1.117)之間，F0917、F0908、KS、F0905、F0921和F0915六個親本具有極顯著的正向加性效應，其他親本具有不同程度的負向效應。全長日增量的加性效應在(0.04±0.006)到(-0.06±0.008)之間，其中F0917、KS、F0905、F0921、F0768、F09125、F09121 七個親本具有顯著或極顯著正向加性效應，其他親本具有不同程度的負向效應。體質量日增量的加性效應在(0.04±0.007)到(-0.06±0.007)之間，正向與負向效應與體長日增量相同。綜合以上5個性狀的加性效應預測值，可以看出F0917、F0908、F0905、F0915、F09125和KS群體的親本在5個性狀上都具有極顯著的正向效應(p<0.01)。其他親本5個性狀都不同程度的表現為顯著(p<0.05)或極顯著(p<0.01)的負向加性效應(見表3)。

表3 親本加性效應預測值(M±S.E)

注：** 表示差異極顯著(p< 0.01)；*表示差異顯著(p<0.05)。

3.4 親本不同雜交組合顯性隨機效應

對F3代家系建立中42個雜交組合的顯性隨機效應進行預測，結果顯示(表4)，全長的顯性效應在(45.55±2.585)到(-32.02±3.17)之間，有19個雜交給表現為不同程度的正效應，其中12個雜交組合表現出極顯著的正效應(p<0.01)。體寬的顯性效應在(17.01±1.031)到(-10.51±1.09)之間，有18個雜交組合表現為不同程度的正效應，其中有11個表現為極顯著正效應(p<0.01)。體質量的顯性效應在(41.77±3.154)到(-21.74±2.84)之間，有13個雜交組合表現為不同程度的正效應，其中有11個表現為極顯著正效應(p<0.01)。在全長、體寬和體質量3個性狀上都表現為極顯著正效應的雜交組合為D(6×6)、D(1×14)、D(2×6)、D(2×11)、D(3×4)、D(3×16)、D(3×17)、D(4×18)、D(7×12)、D(7×14)、D(8×13) 11個雜交組合。

全長日增量的顯性效應在(0.17±0.016)到(-0.10±0.017)之間，有22個雜交組合表現為不同程度的正向效應，其中有14個表現為極顯著正效應(p<0.01)。體質量日增量的顯性效應在(0.18±0.020)到(-0.12±0.018)之間，有22個雜交組合不同程度的表現為正效應，其中有13個表現為極顯著的正效應(p<0.01)。在這兩個性狀上都表現為極顯著正效應的雜交組組合為D(5×5)、D(1×3)、D(1×12)、D(2×6)、D(2×13)、D(2×18)、D(3×4)、D(4×11)、D(6×8)、D(7×12)、D(7×14)、D(8×11)、D(8×15)。對以上結果綜合比較，在5個性狀上都表現為正效應的雜交組合為D(2×6)、D(2×18)、D(3×4)、D(7×12)、D(7×14)、D (8×13)，利用這些雜交組合建立家系將促進牙鲆以上5個性狀的生長。

表4 各雜交組合顯性隨機效應值(M±S.E.)

續表4

注：D表示顯性效應；1～18為親本序號，與表3中親本序號及親本相同；** 表示差異極顯著(p<0.01)；*表示差異顯著(p<0.05)；+表示差異水平為p<0.10。

4 討論

牙鲆在我國的養殖面積相當廣泛，從渤海灣北部到福建南部沿海地區都有牙鲆的養殖[18]。養殖方式有工廠化、池塘和網箱養殖。近年來人們不斷總結養殖經驗，結合南北方海水溫度差異和互補優勢，逐漸發展了一種“池塘-網箱南北接力”養殖模式，當年可將牙鲆養殖到400～500 g，極大的提高了商品魚的養成速度和養殖效益，同時為牙鲆產業化探索了新的途徑。隨著牙鲆養殖產業化的進一步發展，對生長快、抗病力強、溫度適應性廣、抗逆力強牙鲆新品種的需求日益迫切。在牙鲆選擇育種方面，目前雖然已經培育出1個具有生長優勢雜交新品種，但需要對繁殖群體進行不斷的保持和延續，對優良性狀進行不斷的提純和復壯。另外在國內牙鲆養殖新品種的覆蓋率還相當低，為了解決養殖生產中出現的各種問題，還需要針對抗病、耐高溫、抗逆力等方面培育出更多的養殖新品種，來不斷滿足生產發展的需求。

選擇育種是目前牙鲆育種中最主要的育種手段，通過群體選育初步建立了牙鲆黃海野生群體、日本群體、韓國群體和抗鰻弧菌群體4個基礎育種群體[12—13]，從2007-2009年利用這些基礎群體，通過家系建立和選育方式建立了牙鲆F1、F2代家系，并從中篩選出生長快、成活率高的家系[13]。2012年利用F1、F2代優良家系和韓國牙鲆群體構建了F3代家系。本研究擬通過F3代家系的選育，為牙鲆新品種的培育準備具有豐富遺傳性狀的育種材料。

數量性狀選擇是選擇育種的方法之一，數量性狀通常受許多微效基因控制，同時受到環境因素的影響[17]。Fisher最早提出了數量性狀遺傳的“加性—顯性”遺傳模型，為數量遺傳模型的建立奠定了基礎[19]。1953年Henderson提出了3種適用于非均衡資料的Henderson方法，1971年Patterson和Thompson正式提出約束極大似然法(Restricted Maximum Likelihood，REML)用于遺傳方差組分估計，國內外的大量研究證明，REML方法是目前動物育種中遺傳參數估計比較理想的方法[20]。利用動物遺傳模型已經對銀大麻哈魚(Oncorhynchuskisutch)[21]、虹鱒(Oncorhynchusmykiss)[22]、尼羅羅非魚(Oreochromisnilogicus)[23]等魚類的體長和體質量遺傳參數進行了估計。

估計邊際均值是指剔除其他變量影響時算出的均值，能夠更加客觀地評估各家系的生長情況[24]，本文利用生長日齡為協變量對家系生長體長、體寬、體質量和日增質量幾個性狀的邊際平均值進行估計，結果直觀的反映了家系的生長情況。利用前3個性狀進行生長快家系的篩選，具有一定的趨同性，說明全長、體寬、體質量3個性狀的相關性比較好。而根據體質量日增量篩選出的優良家系與其他3個性狀的優良家系具有一定差別，主要是由于不同家系建立日期不同造成。綜合比較以上性狀的篩選結果，1202、1206、1226和1211等家系具有生長快特點。

牙鲆家系的表型性狀不僅受到遺傳因素的影響，同時受到環境因素的影響，其遺傳因素主要由基因的加性效應和顯性效應來決定，在“加性-顯性”線性模型中以加性方差和顯性方差來表現。本文中全長、體寬、體質量、全長日增量和體質量日增量5個生長性狀的顯性方差分量均大于加性分量，說明牙鲆家系在195日齡時以上5個性狀的遺傳效應主要取決于基因的顯性效應。但是利用230日齡F1代牙鲆家系生長數據估計全長和體質量的加性效應高于顯性效應[9]，是否可以推測在195日齡時有利于顯性效應的選擇，而在230日齡時有利于加性效應的選擇。另外這一時期環境影響造成的隨機效應值相對較高，說明環境因素對以上性狀有明顯影響。主要是由于牙鲆家系早期培育要受到家系建立時間、環境溫度、攝食餌料強度、不同培育水缸和操作人員的技術差異等因素的影響，隨著家系生長及標準化培育，這種影響會逐漸降低。

在魚類選擇良種中，育種值的估計是選種和配種方案制定的重要指標，育種值可以利用選擇指數、最佳線性無偏預測BLUP法和標記輔助BLUP等方法進行估計[25]，育種值是指親本的種用價值，在數量遺傳學中把決定數量性狀的基因加性效應值定義為育種值。利用育種值選擇和表型選擇方法分別對大菱鲆家系進行了選育評估，發現育種值的選擇效率高于表型值選擇的效率[26]。利用單性狀重復觀測值及動物模型對虹鱒體質量、體長、肥滿度的個體育種值和綜合育種值進行了估計[27]。本文利用“加性-顯性”模型估計了牙鲆全長等5個性狀的綜合育種值在0.099～0.24。而230日齡F1代牙鲆家系體長和體重育種值分別為0.20和0.26[9]，兩次評估所得到的體長和體質量的育種值相近。

遺傳力是衡量遺傳進展的主要標志，狹義遺傳力為加性效應方差占表型方差的比率，廣義遺傳力為加性方差和顯性方差之和在表型方差中所占比率[17]。一般認為具有0.2以上遺傳力的性狀就有選擇效果。鯉魚(Cyprinuscarpio)早期個體體長和體質量的遺傳力在0.11～0.16范圍內[28]；大黃魚(Larimichthyscrocea)40日齡時的體長和體質量遺傳力分別為0.31和0.29[29]；羅非魚(Oreochromisniloticus)在不同環境和不同世代的遺傳力估計結果具有較大差異(0.00～0.52)[30]，在同一養殖條件下同一世代魚類的收獲體重遺傳力為0.06～0.68[31]。田永勝等以3個牙鲆育種群體為研究對象，估計了體長和體質量的狹義遺傳率分別為0. 199 和0. 256，廣義遺傳率分別為0. 312 和0. 336[9]；Tian等通過完全雙列雜交方式估計了生長早期牙鲆生長性狀的遺傳力，體長遺傳力0.22～0.40、體質量遺傳力0.29～0.48[8]；Liu等估計牙鲆不同時期生長性狀的遺傳力結果為0.13～0.39[32]。本文中牙鲆5個生長性狀的狹義和廣義遺傳力分別在0.099～0.24和0.39～0.52，其中，只有全長的遺傳力達到了中等遺傳力水平(0.24±0.02)(h2>0.2)，其它性狀的遺傳力較低。環境效應造成的隨機效應方差達到了0.48～0.61，平均0.53，同時廣義遺傳力明顯大于狹義遺傳力，這一現象說明利用當前的牙鲆核心群體進行雜交育種，有望獲得較好的效果。另外反映出養殖環境對全長、體寬、體質量、全長日增量和體重日增量的影響較大。這主要由于家系的建立延時較長，繁殖前期和后期建立的家系在生長、環境適應、餌料應用等方面都有較大差異。此外1齡以內的魚苗在遺傳上受母體效應影響較大，從而體現為環境效應大于基因加性、顯性效應，給遺傳選擇方面造成一定的困難。因此對魚類家系性狀的選擇應當達到1～2齡較為適宜。

本次研究結果與已報道的關于牙鲆生長性狀遺傳力的估計結果相比，各性狀的遺傳力相對較小，主要原因可能是研究群體存在差異，目前家系已經為F3代家系，在F1、F2代家系的建立與篩選過程中，不同群體之間的遺傳基因不斷的雜合、家系所包含的遺傳多樣性越來越豐富，這方面的原因也會導致遺傳力變小。

生物性狀的加性效應通過遺傳選擇可以遺傳給后代。Pante等利用加性、加性-顯性、加性-顯性-上位動物模型對虹鱒3個核心群體體質量近交效應進行了分析[33]。之后又利用A、A+F、A+CE、A+CE+F、A+D+F、A+D+CE+F線型模型對3個虹鱒群體近交后代的體質量加性效應、顯性效應和環境效應等進行了分析，加性效應隨著近交效應、環境效應的增加有所降低。顯性方差比率在0%～22%，環境方差比率在0%～6%，有較大的變化[34]。王成輝等利用加性-顯性-上位性遺傳模型，對興國紅鯉(Cyprinuscarpiovar.singuonensis) 、荷包紅鯉(C.carpiovar.wuyuanensis) 和甌江彩鯉(C.carpiovar.color) 雙列雜交子二代的體質量10個形態性狀的加性和上位性遺傳效應進行了分析，結果顯示甌江彩鯉的加性和上位性效應能顯著增加后代的體質量、全長和體長，而興國紅鯉和荷包紅鯉的加性和上位性效應能顯著降低后代的體質量、全長和體長[35]。本文利用加性-顯性模型對19個親本群體的全長、體寬、體質量、體長日增量和體質量日增量5個性狀的加性效應進行了分析，結果顯示F0917、F0908、F0905、F0915、F09125和KS六個親本群體在5個性狀上都具有極顯著的正向效應，說明以上親本在育種中能夠促進牙鲆以上性狀的生長。

利用加性-顯性模型對41個雜交組合的顯性隨機效應進行預測，結果顯示，D(2×6)、D(2×18)、D(3×4)、D(7×12)、D(7×14)、D(8×13) 6個雜交組合的顯性效應在5個性狀中均表現為極顯著正效應，說明利用以上雜交組合進行家系的建立和優良品系的選育可以促進牙鲆體長、體質量的生長。雜交組合中2號親本為F0751、3為KS、4為F0750、6為F0915、7為F0917、8為F0768、12為F09121、13為F0905、14為F09104、18為F09125，其中5個親本具有顯著的正向加性效應。F0750和F0751為F1代家系中篩選出的具有生長快、成活率高的優良家系。F0768為具有抗鰻弧菌病能力的F1代家系。F0905、F0915、F0917、F09121、F09125五個親本為F2代家系中篩選出的具有生長快抗病力強的優良家系。KS為從韓國引進的優良牙鲆群體。可見以上篩選出的雜交組合具有優良的親本來源，因此，在以后的育種中可以充分利用這6個雜交組合來進行雜交育種和新品種選育。

綜上所述，利用牙鲆選育F3代家系生長數據和“加性-顯性”線型模型對牙鲆選擇育種結果進行分析和預測，可為牙鲆選擇育種操作提供可信的育種依據，通過以上遺傳效應分析，篩選出了F1202、F1206、F1226和F1211生長快的家系，F0917、F0908、F0905、F0915、F09125和KS優良親本，以及D(2×6)、D(2×18)、D(3×4)、D(7×12)、D(7×14)、D(8×13) 6個雜交組合，為牙鲆新品種的選育提供明確的技術路線。

[1] 楊正勇，王春曉. 全球視野下中國鲆鰈類養殖業的發展[J]. 中國漁業經濟，2009，6(27)：115-122．

Yang Zhengyong，Wang Chunxiao. On the development of flatfish aquaculture industry of China：a global perspective[J]. Chinese Fisheries Economics，2009，6(27)：115-122．

[2] 田永勝，陳松林，徐田軍，等. 牙鲆不同家系生長性能比較及優良親本選擇[J]. 水產學報，2009，33(6)：901-912．

Tian Yongsheng，Chen Songlin，Xu Tianjun，et al. The comparison of growth performances of Japanese flounder (Paralichthysolivaceus) families and selection of parents with good trait[J].Journal of Fisheries of China，2009，33(6)：901-912．

[3] Liu Y G，Chen S L，Li B F. Genetic differentiation among common and selected hatchery populations of flounder: Evidence from RAPD markers[J]. Biochemical Systematics and Ecology，2007，35: 689-695．

[4] 劉海金，朱曉琛，孫效文，等. 牙鲆5個養殖群體的遺傳多樣性分析[J]. 中國水產科學，2008，15(1)：30-37．

Liu Haijin，Zhu Xiaoshen，Sun Xiaowen，et al. Analysis of genetic structure among 5 cultured stocks of Japanese flounderParalichthysolivaceus[J]. Journal of Fishery Sciences of China，2008，15(1)：30-37．

[5] Zhang Y X，Chen S L，Liu Y G，et al. Major Histocompatibility Complex Class IIB Allele Polymorphism and Its Association with Resistance/Susceptibility to Vibrio anguillarum in Japanese Flounder (Paralichthysolivaceus) [J]. Marine Biotechnology，2006，8：600-610．

[6] Xu T J，Chen S L. Genomic structure of DAA gene and polymorphism within MHC-DAA alleles in Japanese flounder (Paralichthysolivaceus) [J]. Yi Chuan，2009，31: 1020-1028．

[7] Song W，Pang R，Niu Y，et al. Construction of high-hensity genetic linkage maps and mapping of growth-related quantitative trail loci in the Japanese flounder (Paralichthysolivaceus) [J]. PLoS One，2012，7(11): e50404．

[8] Tian Y S，Xu T J，Liang Y，et al. Estimates of genetic and phenotypic parameters for weight and length inParalichthysolivaceus(Temminck et Schlegel) [J]. Acta Oceanology Sinica，2011，30(6):58-64.

[9] 田永勝，徐田軍，陳松林，等.三個牙鲆育種群體親本效應及遺傳參數估計[J].海洋學報，2009，31(6):119-128．

Tian Yongshen，Xu Tianjun，Chen Songlin，et al.Parent effects and estimation of genetic parameters for three Japanese flounder breeding populations[J]. Haiyang Xuebao，2009，31(6):119-128．

[10] Zhu X P，You F，Zhang P J，et al. Effects of cold shock on microtubule organization and cell cycle in gynogenetically activated eggs of olive flounder (Paralichthysolivaceus) [J]. Mar Biotechnol (NY)，2006，8: 312-318．

[11] 劉海金，陸桂，王曉梅，等. 有絲分裂雌核發育牙鲆的微衛星鑒定[J].中國水產科學，2010，17(5)：889-894．

Liu Haijin，Lu Gui，Wang Xiaomei，et al. Identification of mitogynogenetic Japanese flounder (Paralichthysolivaceus) using microsatellite marker[J]. Journal of Fishery Sciences of China，2010，17(5)：889-894．

[12] 陳松林，田永勝，徐田軍，等.牙鲆抗病群體和家系的建立及其生長和抗病性能初步測定[J]. 水產學報，2008，32(5)：665-674．

Chen Songlin，Tian Yongsheng，Xu Tianjun，et al. Development and characterization for growth rate and disease resistance of disease-resistance population and family in Japanese flounder (Paralichthysolivaceus) [J].Journal of Fisheries of China，2008，32(5): 665-674．

[13] 田永勝，陳松林，徐田軍，等. 牙鲆不同家系生長性能比較及優良親本選擇[J]. 水產學報，2009，33(6)：901-912．

Tian Yongshen，Chen Songlin，Xu Tianjun，et al. The comparison of growth performances of Japanese flounder (Paralichthysolivacus) families and selecion of parents with good trait[J]. Journal of Fisheries of China，2009，33(6):901-912．

[14] 全國水產技術推廣總站. 2010水產新品種推廣指南[M]. 北京：中國農業出版社，2011：242-268．

National Fisheries Extension Center. 2010 Promotion of new varieties of aquatic guide[M]. Beijing: Agriculture Press of China，2011：242-268．

[15] Sakamoto T，Okamoto N. Marker-assisted breeding of a lymphocystis disease-resistant Japanese flounder (Paralichthysolivaceus) [J]. Aquaculture ，2007，272：291-295．

[16] 朱軍.遺傳模型分析方法[M].北京：中國農業出版社，1997：58-87．

Zhu Jun. Analysis method of genetic model[M].Beijing：Agriculture Press of China，1997：58-87．

[17] 盛志廉. 數量遺傳學[M]. 北京: 科學出版社，1999．

Sheng Zhilian. Quantitative genetics[M]. Beijing: Science Press，1999．

[18] 田永勝，陳松林，嚴安生，等.牙鲆的胚胎發育[J].水產學報，2004，28(6)：609-616．

Tian Yongsheng，Chen Songlin，Yan Ansheng，et al. Study on the embryonic development ofParalichthysolivaceus[J]. Journal of Fisheries of China，2004，28(6):609-616．

[19] 翟虎渠，王建康.應用數量遺傳[M].北京：中國農業科學技術出版社，2007：1-250．

Zhai Huqu，Wang Jiankang. Application of quantitative genetics[M]. Beijing: Agricultural Science and Technology Press of China，2007：1-250.

[20] 楊澤明，熊遠著，喻傳洲.影響豬遺傳評估效果的主要因素研究[J].華中農業大學學報，2001，20(6):598-602．

Yang Zeming，Xiong Yuanzhu，Yu Chuanzhou. The development of pig genetic evaluation[J]. Journal of Huazhong Agricultural University，2001，20(6):598-602．

[21] Neira R，Diaz N F，Gall G A E，et al. Genetic improvement in coho salmon (Oncorhynchuskisutch): Ⅰ. Selection response and inbreeding depression on harvest weight[J]. Aquaculture，2006，257: 9-17．

[22] Su G S，Liljedahl L E，Gall G A E. Genetic correlations between body weight at different ages and with reproductive traits in rainbow trout[J]. Aquaculture，2002，213: 85-94．

[23] Eknath A E，Bentsen H B，Ponzoni R W，et al. Genetic improvement of farmed tilapias: Composition and genetic parameters of a synthetic base population ofOreochromisniloticusfor selective breeding[J]. Aquaculture，2007，273: 1-14．

[24] 蔡一林，岳永生.水產生物統計[M]. 北京：中國農業出版社，2004．

Cai Yilin，Yue Yongshen. Aquatic biostatistics[M]. Beijing: Agriculture Press of China ，2004．

[25] 欒生，孔杰，王清印.水產動物育種值估計方法及其應用的研究進展[J].海洋水產研究. 2008，29(3)：101-108．

Luan Sheng，Kongjie，Wang Qingyin. Methods and application of aquatic animal breeding value estimation: A review[J]. Marine Fisheries Research，2008，29(3):101-108．

[26] 馬愛軍，王新安，雷霽霖.大菱鲆(Scophthalmusmaximus)不同生長階段體重的遺傳參數和育種值估計[J]. 海洋與湖泊，2009，40(2)：187-194．

Ma Aijun，Wang Xin’an，Lei Jilin.Genetic parameterization for turbotScophthalmusmaximus: implication to breeding strategy[J].Oceanologia et Limnologia Sinica，2009，40(2):187-194．

[27] 王炳謙，劉宗岳，高會江，等.應用重復力模型估計虹鱒生長性狀的遺傳力和育種值[J]. 水產學報，2009，33(2)：182-188．

Wang Bingqian，Liu Zongyue，Gao Huijiang，et al. Estimation of genetic parameters for growth traits in rainbow trout (Oncorhynchusmykiss) with differentmodels[J]. Journal of Fisheries China，2009，33(2):182-188．

[28] Nielsen H M，Odegard J，Olesen I，et al. Genetic analysis of common carp (Cyprinuscarpio) strains：Ⅰ. Genetic parameters and heterosis for growth traits and survival[J]. Aquaculture，2010，304：14-21．

[29] 王曉清，王志勇，何湘蓉. 大黃魚40 日齡體長和體重遺傳力估計[J]. 集美大學學報，2010，15(1)：7-10．

Wang Xiaoqing，Wang Zhiyong，He Xiangrong. Heritability estimate of 40-day-old body length and weight ofPseudosciaenacrocea[J].Journal of Jimei University(Natural Science)，2010，15(1):7-10．

[30] Thodesen J，Rye M，Wang Y X，et al. Genetic improvement of tilapias in China: Genetic parameters and selection responses in growth of Nile tilapia (Oreochromisniloticus) after six generations of multi-trait selection for growth and fillet yield[J]. Aquaculture，2011，(322/323)：51-64．

[31] Bentsen H B，Gjerde B，Nguyen N H，et al. Genetic improvement of farmed tilapias: Genetic parameters for body weight at harvest in Nile tilapia (Oreochromisniloticus) during five generations of testing in multiple environments[J]. Aquaculture，2012，(338/341): 56-65．

[32] Liu Y X，Wang G X，Wang Y F，et al. Estimation of genetic parameters for growth traits of Japanese flounderParalichthysolivaceususing an animal model[J]. Fisheries Science， 2011，77: 87-93．

[33] Pante M J R，Gjerde B，McMillan I. Effect of inbreeding on body weight at harvest in rainbow trout，Oncorhynchusmykiss[J]. Aquaculture，2001，192：201-211．

[34] Pante M J R，Gjerde B，McMillan I. Ignacy Misztal.Estimation of additive and dominance genetic variances for body weight at harvest in rainbow trout，Oncorhynchusmykiss[J]. Aquaculture，2002，204:383-392．

[35] 王成輝，李思發，劉志國，等.紅鯉生長性狀的上位性遺傳效應分析[J].中國水產科學，2006，13(4)：573-579．

Wang Chenghui，Li Sifa，Liu Zhiguo，et al.Genetic analysis of epistatic of growth traits in red common carp[J]. Journal of Fishery Sciences of China，2006，13(4):573-579．

Establishment and estimated genetic effects ofParalichthysolivaceusF3 families

Tian Yongsheng1，Chen Honglin1，2，Qi Wenshan1，2，Jiang Jing1，2，Zhang Yingping1，2，Liu Wanjun1，3，Li Baoyu4，Liu Shoutang4，Sun Deqiang4，Chen Songlin1

(1.KeyLaboratoryforSustainableUtilizationofMarineFisheryResources，MinistryofAgriculture，YellowSeaFisheriesResearchInstitute，ChineseAcademyofFisherySciences，Qingdao266071，China; 2.CollegeofFisheriesandLifeScience，ShanghaiOceanUniversity，Shanghai201306，China3.OceanUniversityofChina，Qingdao266003，China； 4.HuanghaiAquacultureLtd.，Haiyang265100，China)

Paralichthysolivaceusis the major species for marine aquaculture and fisheries on the western Pacific coast. Due to the high disease incidence and slow growth rate observed in artificial-breeding programs，the selective breeding of Japanese flounder has been pursued since the beginning of this century. Here，we established 52 F3 families based on the screened F1 and F2 core families that we bred previously. Approximately 196 days after hatching，we measured five growth traits of 1 560 fry from the 52 families，for a total of 7 800 data points. An “A-D” model and MINQUE were used to estimate the marginal mean，additive-variance component，dominance-variance component，random-variance component，narrow-sense heritability，broad-sense heritability，parental additive effect and random hybridization effect of each trait. The additive-variance component ratios of the five traits were between (0.099±0.01) and (0.24±0.02); the dominance-variance component ratios were between (0.30±0.02) and (0.31±0.02); the random-variance component ratios were between (0.48±0.02) and (0.61±0.01); the narrow-sense heritabilities were between 0.099 and 0.24; and the broad-sense heritabilities were between 0.38 and 0.54. Each of these genetic parameters showed highly significant differences (p<0.01). The random dominance effects of 42 crosses were predicted，and D(2×6)，D(2×18)，D(3×4)，D(7×12)，D(7×14) and D(8×13) all showed highly significant positive effects (p<0.01) on the five traits. The results of this study provide a preliminary set of optimal parental lines，crosses and families to establish improved Japanese flounder strains. These results also provide a theoretical and technical foundation for further selective breeding．

japanese flounder;Paralichthysolivaceus; F3 family; heritability; additive effect; dominance effect

10.3969/j.issn.0253-4193.2015.04.004

2014-07-15；

2014-12-17。

“863”高技術研究發展計劃(2012AA10A408)；國家自然科學基金(31372510)；上海海洋大學研究生科研基金(A1-0209-14-0900-37)；山東省泰山學者建設工程專項資助。

田永勝(1964—)，男，甘肅省會寧縣人，研究員，博士，主要從事魚類低溫生物學及遺傳育種研究。 E-mail:tianys@ysfri.ac.cn

*通信作者：陳松林，研究員。E-mail：chensl@ysfri.ac.cn

S917.4

A

0253-4193(2015)04-0040-12

田永勝，陳紅林，齊文山，等. 牙鲆選育F3代家系的建立及遺傳效應分析[J].海洋學報，2015，37(4)：40—51，

Tian Yongsheng，Chen Honglin，Qi Wenshan，et al. Establishment and estimated genetic effects ofParalichthysolivaceusF3 families[J]. Haiyang Xuebao，2015，37(4)：40—51，doi:10.3969/j.issn.0253-4193.2015.04.004