張躍環(huán) 秦艷平 張 揚(yáng) 李 軍 肖 述 向志明馬海濤 喻子牛①
(1. 中國科學(xué)院南海海洋研究所熱帶海洋生物資源與生態(tài)重點(diǎn)實(shí)驗(yàn)室 廣東省應(yīng)用海洋生物學(xué)重點(diǎn)實(shí)驗(yàn)室 廣州 510301;2. 南海生物資源開發(fā)與利用協(xié)同創(chuàng)新中心 廣州 510275)
對(duì)于貝類選擇育種而言, 連續(xù)幾代的選擇均可以獲得積極的選擇反應(yīng)(張國范等, 2009)。Haskin等(1979)對(duì)美洲牡蠣(Crassostrea virginica)進(jìn)行了連續(xù)四代的抗單孢子蟲病選擇, 使其抗病力比野生群體提高了 8—9倍; Newkirk(1982)認(rèn)為歐洲牡蠣(Ostrea edulis)生長選擇每代提高 10%—20%是合理的;Paynter等(1990)報(bào)道了養(yǎng)殖人員經(jīng)過18代選擇的美洲牡蠣, 其生長速度顯著快于野生群體; Nell等(1999)對(duì)悉尼巖牡蠣(Saccostrea commercialis)體重進(jìn)行了連續(xù) 2代選擇, 獲得了 17.8%的遺傳改進(jìn)量; Nell等(2006)進(jìn)行悉尼巖牡蠣三代的抗單孢子蟲病(QX疾病,Queensland Unknown Disease)選育, 獲得了較強(qiáng)的抗病系; Li等(2011)和Wang等(2012)對(duì)長牡蠣(C. gigas)進(jìn)行了連續(xù)6代的混合選擇, 獲得了積極的改良效果,培育出“海大1號(hào)”長牡蠣新品種; Dégremont等(2015)對(duì)長牡蠣進(jìn)行了連續(xù)4代的抗皰疹病毒選育, 獲得了抗性較強(qiáng)的牡蠣新品系, 存活力提高了61.8%。
至于其它常見經(jīng)濟(jì)貝類, Zheng等(2004, 2006)對(duì)海灣扇貝(Argopecten irradians)進(jìn)行了連續(xù)4代的混合選擇, 發(fā)現(xiàn)基礎(chǔ)群體對(duì)其遺傳效果具有顯著影響,培育出“中科紅”海灣扇貝; Zheng等(2013)還對(duì)華貴櫛孔扇貝(Chlamys nobilis)進(jìn)行了橘紅色閉殼肌品系連續(xù)4代的選育, 培育出“南澳1號(hào)”華貴櫛孔扇貝新品種; 閆喜武等(2010)和劉振等(2014)對(duì)菲律賓蛤仔(Ruditapes philippinarum)進(jìn)行了連續(xù) 8代混合選擇,培育出“斑馬蛤”菲律賓蛤仔新品種; You等(2009)對(duì)九孔鮑(Haliotis diversicolor)進(jìn)行了連續(xù) 4代混合選擇, 培育出“東優(yōu)1號(hào)”雜色鮑新品種。
先前, 我們對(duì)熊本牡蠣(C. sikamea)的多嵴、無嵴品系進(jìn)行了一次上選, 經(jīng)過一周年的養(yǎng)成, 獲得了比較好的選擇反應(yīng), 并且獲得了8.41%、7.71%遺傳改進(jìn)量(張躍環(huán)等, 2016)。本研究是以兩個(gè)品系子一代選擇組為材料, 繼續(xù)觀測(cè)兩個(gè)品系的選擇反應(yīng)、現(xiàn)實(shí)遺傳力, 計(jì)算二次選擇的遺傳改進(jìn)量, 為熊本牡蠣遺傳育種提供理論基礎(chǔ)和實(shí)踐經(jīng)驗(yàn)。
如圖1所示, 本研究以2012年5月采集的湛江野生群體熊本牡蠣為材料, 構(gòu)建了核心基礎(chǔ)群。于2013年5月, 從核心基礎(chǔ)群中篩選出多嵴和無嵴品系,按照 10%留種率, 進(jìn)行了兩個(gè)品系熊本牡蠣子一代選擇反應(yīng)及其現(xiàn)實(shí)遺傳力研究。在此基礎(chǔ)上, 于2014年 5月, 以多嵴、無嵴品系熊本牡蠣上選的子一代作為材料, 隨機(jī)篩選每個(gè)品系若干個(gè)個(gè)體(剔除掉外部形態(tài)不整齊個(gè)體)作為育種材料進(jìn)行子二代連續(xù)上選。
利用截頭法分別對(duì)多嵴和無嵴品系子一代進(jìn)行混合選擇, 留種率均為10%, 選擇強(qiáng)度(i)均為1.755。從兩種類型中各隨機(jī)抽取 1200個(gè)個(gè)體, 用游標(biāo)卡尺進(jìn)行測(cè)量(精確至 0.02mm), 繪制個(gè)體殼高大小的頻率分布圖(圖 1)。從每個(gè)群體中選擇 10%的最大個(gè)體(120個(gè)), 作為選擇組2(selected group 2, SG2), 在剩余的個(gè)體中隨機(jī)抽取120個(gè)個(gè)體作為對(duì)照組2(control group2, CG2)。其中, 多嵴品系親本的截點(diǎn)為62.37mm, 雌雄比例為 87:33; 無嵴品系親本的截點(diǎn)為57.46mm, 雌雄比例為81:39(圖1)。

圖1 兩個(gè)品系熊本牡蠣子二代選擇育種技術(shù)路線圖Fig.1 Pedigree of selected strains indicating founder strain and selection intensity(i) for a successive two-generation selection
采用解剖方式獲取配子, 為了保證所有雄性個(gè)體的繁殖機(jī)會(huì), 每個(gè)實(shí)驗(yàn)組均采用單個(gè)雄性個(gè)體與所有卵子配對(duì)的原則進(jìn)行(1♂:n♀)。實(shí)驗(yàn)中為了防止發(fā)生組間污染, 利用300目篩絹網(wǎng)進(jìn)行幼蟲選育, 各實(shí)驗(yàn)組均保留D型幼蟲約1.2億個(gè), 多余幼蟲全部放掉; 每個(gè)實(shí)驗(yàn)組放置在3個(gè)10m3水泥池中培養(yǎng), 幼蟲培養(yǎng)密度控制在4ind./mL。
幼蟲培育期間, 鹽度控制在 20—25; 稚貝培養(yǎng)的鹽度控制在15—20。室內(nèi)培育期間, 6日齡以前, 投喂云微藻(Chlorella sp.), 之后混合投喂云微藻和牟氏角毛藻(Chaetoceros muelleri)(體積比=1:1), 投餌量視攝食情況而定, 日投喂量在(5—10)×104cell/mL。當(dāng)幼蟲出現(xiàn)鰓原基、眼點(diǎn)、足時(shí), 利用波紋板采苗。當(dāng)稚貝殼高在10mm左右時(shí)裝入5mm網(wǎng)孔塑料袋在竹林鹽廠生態(tài)池中進(jìn)行中間育成, 每袋密度控制在240—300個(gè)。隨著稚貝生長, 定期調(diào)整密度, 最后移至孔徑 10mm的扇貝籠中養(yǎng)成, 密度控制在每層50—60個(gè)。
殼高小于 400μm 的幼蟲在顯微鏡下用目微尺(10×)測(cè)量, 殼長大于 3.0mm 的子代和親本用游標(biāo)卡尺測(cè)量, 精確值0.02mm。每個(gè)實(shí)驗(yàn)組設(shè)3個(gè)重復(fù), 每個(gè)重復(fù)隨機(jī)測(cè)量30個(gè)個(gè)體。
為了保證方差齊性, 所有的殼高進(jìn)行以 10為底的對(duì)數(shù)轉(zhuǎn)化。用SPSS19.0統(tǒng)計(jì)軟件對(duì)數(shù)據(jù)進(jìn)行分析處理, 不同實(shí)驗(yàn)組間數(shù)據(jù)的比較采用雙因素方差分析方法(Two-way ANOVA), 用Turkey法進(jìn)行組間多重比較, 差異顯著性設(shè)置為P<0.05。
參照 Zheng等(2006)使用的方法, 將兩個(gè)品系熊本牡蠣子二代的選擇反應(yīng)(Selected reaction, SR)、現(xiàn)實(shí)遺傳力(hR2)及遺傳改進(jìn)量(Genetic gain, GG)用以下公式計(jì)算

公式(1)、(2)和(3)中: XSG2和 XCG2分別是選擇組(SG2)和對(duì)照組(CG2)個(gè)體的平均殼高, SCG2是對(duì)照組(CG2)的標(biāo)準(zhǔn)差, i是選擇強(qiáng)度。
為了確定品系來源及其選擇效應(yīng)對(duì)生長性狀的影響, 參照 Zhang等(2007)方法, 利用雙因子分析模型檢測(cè)品系來源(strain origin, SO)及選擇效應(yīng)(selective effect, SE)對(duì)各實(shí)驗(yàn)組子代殼高的影響

式中: Yijk為k個(gè)重復(fù)i個(gè)品系j種選擇下的殼高(存活率); μ為常數(shù); SOi為殼高的品系來源(i=1, 2); SEj為殼高的選擇效應(yīng)(j=1, 2); (SO×SE)ij為品系與選擇效應(yīng)的交互作用; eijk為隨機(jī)誤差(k=1, 2, 3)。
幼蟲期, 3日齡時(shí), 多嵴品系選擇組殼高顯著大于其它實(shí)驗(yàn)組(P<0.05); 9—15日齡期間, 多嵴品系選擇組殼高顯著大于無嵴品系對(duì)照組, 而多嵴品系對(duì)照組與無嵴品系無顯著差異(P>0.05)。經(jīng)過雙因子分析模型分析, 品系來源及其選擇效應(yīng)均對(duì)子二代幼蟲殼高造成了顯著影響, 15日齡時(shí)二者間的交互作用也造成了顯著影響(P<0.05)。稚貝期(30—60日齡),多嵴品系選擇組殼高顯著大于無嵴品系對(duì)照組, 而其它兩個(gè)實(shí)驗(yàn)組間無顯著差異(P>0.05)。在此期間,選擇效應(yīng)對(duì)子二代殼高具有顯著影響。養(yǎng)成期,90—180日齡期間, 多嵴品系選擇組殼高一直顯著大于無嵴品系對(duì)照組(P<0.05); 360日齡時(shí), 多嵴品系選擇組顯著大于其它三個(gè)實(shí)驗(yàn)組, 無嵴品系選擇組顯著大于相應(yīng)對(duì)照組(P<0.05), 表現(xiàn)出明顯的選擇效果。180—360日齡期間, 品系來源及其選擇效應(yīng)均對(duì)子二代殼高具有顯著影響(P<0.05), 二者間交互作用尚未發(fā)生顯著影響(表1, 2)。
從遺傳參數(shù)上看, 兩個(gè)品系選擇反應(yīng)、現(xiàn)實(shí)遺傳力均表現(xiàn)為幼蟲期<稚貝期<養(yǎng)成期。多嵴、無嵴品系幼蟲期的選擇反應(yīng)分別為0.47±0.12、0.34±0.06, 稚貝期為 0.65±0.05、0.40±0.05, 養(yǎng)成期為 0.82±0.07、0.56±0.07; 幼蟲期的現(xiàn)實(shí)遺傳力為 0.27±0.07、0.19±0.04, 稚貝期為 0.37±0.03、0.23±0.03, 養(yǎng)成期為 0.47±0.04、0.32±0.04。從其遺傳改進(jìn)量上看, 多嵴、無嵴品系品系幼蟲期為 4.83%±1.01%、3.61%±0.71%, 稚貝期為7.55%±0.40%、5.76%±0.22%, 養(yǎng)成期為9.95%±0.18%、6.47%±0.13%。經(jīng)過一周年的養(yǎng)成, 多嵴、無嵴品系的選擇反應(yīng)、現(xiàn)實(shí)遺傳力和遺傳改進(jìn)量分別為0.88、0.62, 0.50、0.35, 10.13%、6.61% (表 3, 圖 2)。

表1 兩個(gè)品系子二代選擇組和對(duì)照組殼高Tab.1 Shell height of the control and selective groups in different ages

表2 品系來源(SO)及選擇效應(yīng)(SE)對(duì)其子二代生長性狀影響的方差分析Tab.2 Analysis of variance showing the strain origin (SO) and selective effects (SE) on growth for each experimental group

表3 兩個(gè)品系子二代的選擇反應(yīng)(SR)、現(xiàn)實(shí)遺傳力(hR2)及遺傳改進(jìn)量(GG)Tab.3 Response to selection (SR), realized heritability (hR2)and genetic gain (GG) of the two oyster strains
本研究中, 通過設(shè)置親本子代對(duì)照組來計(jì)算遺傳參數(shù), 非用親本本身, 這主要是因?yàn)樨愵惖纳L經(jīng)常受到環(huán)境效應(yīng)影響, 如水質(zhì)、餌料、水流等影響著貝類的生長性能(Zhenget al, 2004, 2006; 李冬春等,2013)。研究中發(fā)現(xiàn), 兩個(gè)品系子一代的殼高就存在著顯著差異, 即多嵴品系子一代顯著大于無嵴品系, 且其選擇反應(yīng)均較高(張躍環(huán)等, 2016), 所以子二代還存在著一定程度的遺傳改良效應(yīng)。為此, 我們進(jìn)行了兩個(gè)品系的子二代連續(xù)上選。
經(jīng)過兩代連續(xù)上選, 兩個(gè)品系遺傳改進(jìn)量存在著分歧, 多嵴品系遺傳改進(jìn)量一直大于無嵴品系, 尤其是到了360日齡時(shí), 多嵴品系和無嵴品系遺傳改進(jìn)量分別為10.13%、6.61%。這種分歧可能與品系內(nèi)遺傳多樣性變化水平有關(guān), 也就是說多嵴品系內(nèi)遺傳多樣性水平較高, 而無嵴品系內(nèi)遺傳多樣性水平較低(Liet al, 2011; Wanget al, 2012)。在貝類連續(xù)選擇中, Newkirk等(1982)認(rèn)為歐洲牡蠣生長選擇每代提高 10%—20%是合理的; Nell等(1999)對(duì)悉尼巖牡蠣體重進(jìn)行了連續(xù)兩代選擇, 獲得了17.8%的遺傳改進(jìn)量; Zheng等(2006)對(duì)海灣扇貝進(jìn)行了兩代連續(xù)上選,發(fā)現(xiàn)新引進(jìn)的海灣扇貝群體可以獲得 10%左右的遺傳改進(jìn)量; 張國范等(2010)對(duì)菲律賓蛤仔進(jìn)行中日韓三個(gè)群體的連續(xù)兩代的混合選擇, 獲得了 10%左右遺傳改進(jìn)量; 王慶志(2011)對(duì)長牡蠣進(jìn)行了中日韓三個(gè)地理群體的三代混合選擇, 每代均獲得了 10%的遺傳改進(jìn)量。由此可見, 本研究中子二代獲得遺傳改進(jìn)量處于較低水平, 這可能與熊本牡蠣自身個(gè)體較小, 生長緩慢有著直接關(guān)系。

圖2 兩個(gè)品系不同階段平均的選擇反應(yīng)(SR)、現(xiàn)實(shí)遺傳力(hR2)及遺傳改進(jìn)量(GG)Fig.2 Response to selection (SS), realized heritability (hR2) and genetic gain (GG) of the two oyster strains in different stages
從遺傳參數(shù)看, 在整個(gè)生活史中, 多嵴品系的選擇反應(yīng)及現(xiàn)實(shí)遺傳力均高于無嵴品系, 尤其是經(jīng)過周年養(yǎng)成的成體, 二者間差異明顯。我們分析, 這可能是兩個(gè)品系間的遺傳多樣性水平差異不同, 多嵴品系中仍存在著較大的遺傳方差, 大的遺傳方差暗示蘊(yùn)藏著大的遺傳力, 遺傳力大, 遺傳改進(jìn)量相對(duì)較多, 這也是多嵴品系遺傳改進(jìn)量較多的原因(Zheng et al, 2004, 2006; 馬愛軍等, 2012)。此外, 在熊本牡蠣野生群體中, 多嵴類型占據(jù)大多數(shù), 無嵴類型相對(duì)較少, 可能導(dǎo)致無嵴品系遺傳多樣性水平較低。
與其子一代選擇相比, 多嵴品系的現(xiàn)實(shí)遺傳力較高, 說明還具有較大的選擇空間; 而無嵴品系的現(xiàn)實(shí)遺傳力相對(duì)居中, 較一代選擇時(shí)結(jié)果略有下降, 這可能是遺傳多樣性水平變化小, 或者不存在較大的效應(yīng)基因?qū)е碌?張國范等, 2009; 張躍環(huán)等, 2014)。經(jīng)過雙因子模型分析發(fā)現(xiàn), 兩個(gè)品系經(jīng)過一周年的養(yǎng)成, 仍存在著顯著的生長分歧, 這也說明熊本牡蠣多嵴品系子二代具有較好的遺傳改良效應(yīng), 連續(xù)上選的效果是顯著的。
本文發(fā)現(xiàn)經(jīng)過連續(xù)兩次的上選, 兩個(gè)品系生長性狀仍具有較高的選擇反應(yīng), 且多嵴品系的遺傳參數(shù)及遺傳改進(jìn)量均明顯大于無嵴品系, 出現(xiàn)了遺傳分歧, 這也說明多嵴品系仍具有較好的遺傳改良前景。本文為熊本牡蠣遺傳改良和新品種培育奠定了理論與實(shí)踐基礎(chǔ)。
馬愛軍, 王新安, 黃智慧等, 2012. 大菱鲆(Scophthalmus maximus)家系選育 F2早期選擇反應(yīng)和現(xiàn)實(shí)遺傳力估計(jì).海洋與湖沼, 43(1): 57—61
王慶志, 李 琪, 劉士凱等, 2011. 長牡蠣不同地理群體選育系數(shù)量性狀的比較. 中國海洋大學(xué)學(xué)報(bào), 41(7—8): 36—41
劉 振, 閆喜武, 張躍環(huán)等, 2014. 菲律賓蛤仔不同品系斑馬蛤的選擇反應(yīng)及現(xiàn)實(shí)遺傳力. 海洋與湖沼, 45(2): 259—264
閆喜武, 張躍環(huán), 霍忠明等, 2010. 不同地理群體菲律賓蛤仔的選擇反應(yīng)及現(xiàn)實(shí)遺傳力. 水產(chǎn)學(xué)報(bào), 34(5): 704—710
李冬春, 閆喜武, 霍忠明等, 2013. 蛤仔(Ruditapes philippinarum)養(yǎng)成期殼長遺傳力與育種值估計(jì). 海洋與湖沼, 44(6): 1550—1556
張國范, 閆喜武, 2010. 蛤仔養(yǎng)殖學(xué). 北京: 科學(xué)出版社, 287—291
張國范, 鄭懷平, 2009. 海灣扇貝養(yǎng)殖遺傳學(xué). 北京: 科學(xué)出版社, 52—73
張躍環(huán), 王昭萍, 喻子牛等, 2014. 香港牡蠣(Crassostrea hongkongensis)與長牡蠣(Crassostrea gigas)種間雜種遺傳力評(píng)估. 海洋與湖沼, 45(6): 1367—1373
張躍環(huán), 張 揚(yáng), 李軍等, 2016. 熊本牡蠣多嵴和無嵴品系 F1的選擇反應(yīng)及現(xiàn)實(shí)遺傳力研究. 中國水產(chǎn)科學(xué), 已錄用
Dégremont L, Nourry M, Maurouard E, 2015. Mass selection for survival and resistance to OsHV-1 infection in Crassostrea gigas spat in field conditions: response to selection after four generations. Aquaculture, 446: 111—121
Haskin H H, Ford S E, 1979. Development of resistance to Minchinia nelsoni (MSX) mortality in laboratory-reared and native oyster stocks in Delaware Bay. Marine Fisheries Review, 41(1): 54—63
Li Q, Wang Q Z, Liu S K et al, 2011. Selection response and realized heritability for growth in three stocks of the Pacific oyster Crassostrea gigas. Fisheries Science, 77(4): 643—648
Nell J A, Smith I R, Sheridan A K, 1999. Third generation evaluation of Sydney rock oyster Saccostrea commercialis(Iredale and Roughley) breeding lines. Aquaculture,170(3—4): 195—203
Nell J A, Perkins B, 2006. Evaluation of the progeny of third-generation Sydney rock oyster Saccostrea glomerata(Gould, 1850) breeding lines for resistance to QX disease Marteilia sydneyi and winter mortality Bonamia roughleyi.Aquaculture Research, 37(7): 693—700
Newkirk G F, Haley L E, 1982. Phenotypic analysis of the European oyster Ostrea edulis L.: relationship between length of larval period and postsetting growth rate. Journal of Experimental Marine Biology and Ecology, 59(2—3): 177—184
Paynter K T, Dimichele L, 1990. Growth of tray-cultured oysters(Crassostrea virginica Gmelin) in Chesapeake Bay.Aquaculture, 87(3—4): 289—297
Wang Q Z, Li Q, Kong L F et al, 2012. Response to selection for fast growth in the second generation of Pacific Oyster(Crassostrea gigas). Journal of Ocean University of China,11(3): 413—418
You W W, Ke C H, Luo X et al, 2010. Divergent selection for shell length in two stocks of small abalone Haliotis diversicolor. Aquaculture Research, 41(6): 921—929
Zhang H B, Liu X, Zhang G F et al, 2007. Growth and survival of reciprocal crosses between two bay scallops, Argopecten irradians concentricus Say and A. irradians irradians Lamarck. Aquaculture 272(S1): 88—93
Zheng H P, Zhang G F, Liu X et al, 2004. Different responses to selection in two stocks of the bay scallop, Argopecten irradians irradians Lamarck (1819). Journal of Experimental Marine Biology and Ecology, 313(2): 213—223
Zheng H P, Zhang G F, Liu X et al, 2006. Sustained response to selection in an introduced population of the hermaphroditic bay scallop Argopecten irradians irradians Lamarck (1819).Aquaculture, 255(1—4): 579—585
Zheng H P, Zhang T, Sun Z W et al, 2013. Inheritance of shell colours in the noble scallop Chlamys nobilis (Bivalve:Pectinidae). Aquaculture Research, 44(8): 1229—1235