韓龍海，潘鳳剛，劉洪志，宿效忠，于連軍，張含國

(1.葦河林業(yè)局青山林木種子園，黑龍江 尚志 150600；2.東北林業(yè)大學林學院，黑龍江 哈爾濱 150040)

紅松(Pinuskoraiensis)主要分布于長白山、老爺嶺、張廣才嶺、完達山和小興安嶺，是東北地區(qū)重要的用材與堅果樹種，是天然林主要的森林組成樹種[1-3].由于營養(yǎng)豐富、含油量高，關(guān)于紅松種子在健康食品[4-5]、食用油加工[6-7]等方面的研究日益增多.已經(jīng)證實，紅松種子產(chǎn)量、營養(yǎng)品質(zhì)在產(chǎn)地、群體間均存在較大變異[8-11].目前，紅松種子園無性系間種實形態(tài)和營養(yǎng)成分變異的研究報道較少[12-13]，而紅松無性系偏雌、偏雄現(xiàn)象普遍存在[14]，因此，開展種實性狀變異分析及優(yōu)良無性系選擇十分必要.本次研究通過調(diào)查葦河1984年營建的紅松種子園近81個無性系多年結(jié)實量和種實性狀，研究無性系種實變異規(guī)律；通過穩(wěn)定性和主成分分析選擇出高產(chǎn)、穩(wěn)產(chǎn)的優(yōu)良無性系，為當?shù)丶t松堅果林營建及遺傳改良提供參考.

1 材料與方法

1.1 試驗地點與材料來源

試驗地點位于黑龍江省尚志市葦河林業(yè)局(44°40′N，128°03′E)，平均海拔300 m，年降水量666.1 mm，年日照時間2 552.3 h，7月均溫21.6 ℃，年均溫2.3 ℃，大于5 ℃積溫2 753.2 ℃，年蒸發(fā)量1 084.4 mm.1984年，按順序5錯位法建立紅松無性系種子園36.6 hm2，共238個無性系，每個無性系29～37株.優(yōu)樹選自鶴北林業(yè)局天然母樹林.

1.2 數(shù)據(jù)調(diào)查

2019年，調(diào)查8、9、15區(qū)81個無性系的單株球果數(shù)、單株球果重、單果重、球果長、球果寬、種子長、種子寬、種仁重、種皮重、種子重、千粒重、種子粒數(shù)、出種率、出仁率、空殼率.2000年、2011年、2014年調(diào)查球果數(shù)量、球果產(chǎn)量、種子產(chǎn)量、千粒重、出種率、出仁率.其中，球果數(shù)為無性系單株球果數(shù)量；球果產(chǎn)量為無性系單株球果總重；種子產(chǎn)量為無性系單株球果脫粒后種子總重；出種率為種子(含外種皮)質(zhì)量占單個球果質(zhì)量的比；出仁率為單個球果種仁質(zhì)量與帶殼種子質(zhì)量比；空殼率為種子空殼數(shù)量與所有種子數(shù)量的比.調(diào)查40個無性系球果長、寬，種子長、寬，用游標卡尺測定，種子的縱軸長為種長，垂直種臍、種面的最大橫向?qū)挾葹榉N寬；采用四分法調(diào)查千粒重，取平均值.每個無性系測30粒種子，重復(fù)3次.

1.3 數(shù)據(jù)處理與分析

利用SPSS 22.0和DPS V14.10進行數(shù)量性狀統(tǒng)計和分析，百分率數(shù)據(jù)經(jīng)過反正弦平方根轉(zhuǎn)換.

2 結(jié)果分析

2.1 無性系球果產(chǎn)量

球果調(diào)查結(jié)果顯示：結(jié)實大于2株的無性系共有51個.單株平均球果數(shù)為6.3個，排名前10的無性系分別為44、148、42、180、11、200、112、50、62、41號無性系，平均為10.7個，大于無性系均值70.7%；單株球果重排名前10的無性系為44、42、11、148、180、200、62、50、112、115號，平均為1 818.8 g，大于無性系均值69.5%.單株球果數(shù)與單株球果重較高的前10個無性系基本一致，球果產(chǎn)量無性系間差異不顯著(P球果數(shù)為0.391，P球果重為0.322).

2.2 種實性狀

81個無性系種實性狀方差分析結(jié)果見表1.由表1可知：15個種實性狀在81個無性系間差異極顯著.從整體變異系數(shù)看，除空殼率外，其余性狀變異相對較小，變幅為2.36%～13.39%.球果重變異較大的無性系有131、68、137、63、40號，平均為18.8%；變異系數(shù)較小的無性系有74、64、118、2、65號，平均為1.3%.種子重變異系數(shù)平均為7.4%，變異系數(shù)較大的無性系有60、68、146、56、113號，平均為17.9%；變異系數(shù)較小的無性系有53、103、66、42、125號，平均為1.3%.出種率變異系數(shù)平均為8.4%，變異系數(shù)較大的無性系有120、39、61、146、129號，平均為21.0%；變異系數(shù)較小的無性系有62、45、88001、119、46號，平均為1.6%.球果長變異系數(shù)平均為12.2%，變異系數(shù)較大的無性系有38、48、74、58、179號，平均為18.3%；變異系數(shù)較小的無性系有119、191、63、121、126號，平均為6.8%.

從優(yōu)到劣對81個無性系的9個種實質(zhì)量性狀和4個形態(tài)性狀排序，前15位排序結(jié)果見表2.由表2可知：140號無性系的9個質(zhì)量性狀中有7個表現(xiàn)優(yōu)良，42、58、64號有6個性狀表現(xiàn)優(yōu)良，45、48、49、50、200、88010號有5個表現(xiàn)優(yōu)良，38、72、116號有4個表現(xiàn)優(yōu)良.種實較重、球果種子較大的無性系有140、42、58、48、49、72號.

表1 81個無性系種實性狀方差分析

表2 各性狀排名前15位的無性系與均值比較

2.3 性狀相關(guān)分析

對81個無性系球果、種實11個性狀進行皮爾遜相關(guān)分析，結(jié)果見表3.由表3可知：球果重的無性系種仁重、種皮重、種子重較大，種子粒數(shù)較多，這些性狀與千粒重、出種率之間的相關(guān)性均達到顯著水平.由此表明，球果較重的無性系，其千粒重和出種率較高，直接影響到單株種子產(chǎn)量.種子長與單果重、種子重、千粒重、出種率，以及種子寬與千粒重、空殼率正相關(guān)顯著或極顯著.

表3 紅松種實性狀相關(guān)性

表3(續(xù))

2.4 多性狀綜合選擇

對2019年調(diào)查的無性系球果重、種仁重、種皮重、種子重、千粒重、出種率、出仁率、空殼率、種子粒數(shù)、種子質(zhì)量10個性狀進行主成分分析，結(jié)果見表4.由表4可知：前3個主成分累計貢獻率達83.0%，能夠反映種實性狀的絕大部分信息.F1中球果重、種仁重、種皮重、種子重、出種率、種子粒數(shù)權(quán)重較大，代表了紅松結(jié)實的“質(zhì)量性狀”；F2與出仁率和空殼率的相關(guān)性較大，代表了“籽粒飽滿性狀”；F3代表千粒重.3個主成分劃分與性狀相關(guān)性分析結(jié)果基本一致.

表4 紅松種實性狀主成分分析

根據(jù)成分得分系數(shù)矩陣分別計算得到3個主成分計算公式：

y1=0.154x1+0.18x2+0.177x3+0.186x4+0.069x5+0.147x6+0.047x7-0.068x8+0.161x9+0.08x10，

y2=0.103x1-0.131x2+0.151x3+0.053x4-0.093x5-0.018x6-0.473x7+0.444x8+0.115x9+0.175x10，

y3=-0.033x1-0.02x2+0.061x3+0.037x4+0.837x5+0.039x6-0.18x7+0.04x8-0.42x9+0.077x10，

其中：x1為球果重；x2為種仁重；x3為種皮重；x4為種子重；x5為千粒重；x6為出種率；x7為出仁率；x8為空殼率；x9為種子粒數(shù)；x10為種子質(zhì)量.以每個主成分對應(yīng)的特征值占所提取主成分總特征值之和的比例為權(quán)重，各系數(shù)乘以權(quán)重之和為綜合模型系數(shù)，得出主成分綜合模型：

y=0.118x1+0.083x2+0.156x3+0.136x4+0.132x5+0.095x6-0.101x7+0.063x8+0.075x9+0.101x10.

對81個無性系各性狀測定值進行標準化轉(zhuǎn)化，得到性狀值標準化向量矩陣.將標準化向量分別帶入主成分計算公式，得到無性系各性狀主成分得分和排名.以10%的入選率選出綜合得分前10%的無性系，共8個，分別為58、50、88010、42、140、64、45、48號無性系.入選無性系與總體無性系均值見表5.由表5可知：入選的8個無性系各性狀表現(xiàn)均優(yōu)于無性系平均值.其中，58號無性系單果重、種仁重、種子粒數(shù)、出種率表現(xiàn)優(yōu)良；88010和42號無性系單株產(chǎn)量高；8個無性系種子重、種仁重、種子粒數(shù)、出種率及種子產(chǎn)量分別大于全部無性系平均值的44.2%、51.8%、25.1%、20.8%及85.0%.

表5 入選無性系種實性狀分析

2.5 結(jié)實量多年穩(wěn)定性

利用DPS軟件對2000、2011、2014、2019年共有的23個無性系種子產(chǎn)量、球果產(chǎn)量和球果數(shù)量進行年度穩(wěn)定性分析，結(jié)果見表6.由表6可知：種子產(chǎn)量豐產(chǎn)、穩(wěn)產(chǎn)較好的無性系是58、60、74、11、47、200、56號.其中，58、60、74號產(chǎn)量較高，56、200、47號穩(wěn)定性較好.球果數(shù)量豐產(chǎn)、穩(wěn)產(chǎn)的無性系為58、60、74、47號.

表6 種子產(chǎn)量豐產(chǎn)性及穩(wěn)定性

3 結(jié)論與討論

研究結(jié)果表明，葦河林業(yè)局青山紅松種子園無性系間種實性狀差異顯著，單株球果數(shù)量、球果產(chǎn)量差異不顯著，可能是由于紅松結(jié)實存在大、小年的規(guī)律[15]，在某些年份容易出現(xiàn)無性系間結(jié)實量差異不顯著現(xiàn)象.

多性狀綜合分析篩選出的8個無性系分別為58、50、88010、42、140、64、45、48號，入選無性系的種子重、種仁重、種子粒數(shù)、出種率及種子產(chǎn)量分別大于全部無性系平均值的44.2%、51.8%、25.1%、20.8%及85.0%，空殼率小于全部無性系均值.產(chǎn)量性狀多年綜合表現(xiàn)優(yōu)良且穩(wěn)定的有58、60、74號無性系，3個無性系種子產(chǎn)量分別大于無性系均值32.2%、46.1%、38.0%.

性狀相關(guān)性分析和主成分分析均表明：紅松單株產(chǎn)量與球果重、單果種子重、出種率顯著正相關(guān)，與出仁率、空殼率和千粒重等性狀相關(guān)不顯著.由此說明，紅松種子大小和飽滿程度與種子產(chǎn)量無關(guān)，這與那冬晨[16]的研究結(jié)果一致.紅松的產(chǎn)量性狀和種仁性狀可以分別進行選擇，有利于選出產(chǎn)量高、品質(zhì)好的無性系.種子長與單果重、種子重、千粒重、出種率正相關(guān)顯著，這為紅松果用無性系間接選擇提供了簡便、快捷的選擇方法.