林 琳 高宇婷 程福山 辛本花 王貴春 夏富才 穆懷志*
(1.北華大學林學院,吉林 132013; 2.吉林省蛟河林業實驗區管理局,蛟河 132517)
賽黑樺(Betulaschmidtii)又稱鐵樺、遼東樺,主要分布于我國東北東部、朝鮮半島北部及俄羅斯錫霍特山脈南部地區,其木質厚重、堅硬,是東北林區最硬的木材之一[1]。由于賽黑樺生長緩慢,成熟林蓄積量低,先前人們并沒有將其作為主要的商品材樹種。隨著東北林區樹種資源的結構性調整,賽黑樺的采伐量逐年增加,目前其年采伐量可達10 000 m3,這使得賽黑樺成為一個單獨用材樹種[2]。連年不斷的采伐,致使賽黑樺資源日漸枯竭,吉林省已將其列為三級保護植物[3]。目前,關于賽黑樺的研究僅限于組織培養[4~5]、木材加工[6]、抗逆性測定[7~8]和生物信息學分析[9]等方面,對于其良種選育的研究尚未見報道。
種子萌發是苗木生長的基礎,種子活力是評價種子品質的重要指標,具有高活力的種子不僅發芽時間短、出苗整齊,而且生成的幼苗具有明顯的生長優勢和適應性[10~11]。先前人們開展了白樺[12]、光皮樺[13]、紅樺[14]等樺木屬植物的種子活力研究,篩選出多個樺木優良家系。本研究對3個采種地共30個賽黑樺半同胞家系的種子活力進行比較,對賽黑樺采種母樹進行初步評價,為賽黑樺優良家系的選擇提供參考。
2018年8月末,分別在吉林省蛟河林業實驗區管理局、琿春林業局三道溝林場、集安市五女峰國家森林公園有賽黑樺分布的林分,選擇樹干通直圓滿、胸徑粗大、無病蟲害的單株作為候選樹。以候選樹為中心,在其周圍10~15 m的范圍內,測量其他賽黑樺的樹高和胸徑,以100 m2內樹高不低于平均高、胸徑達到平均值120%為標準選擇優樹進行采種[15]。3個采種地共選擇30株優樹,分別命名為JH1~JH10、HC1~HC10、JA1~JA10,采集其半同胞家系種子經過室內晾干、篩去苞片后,置于-20℃冰箱中保存。
每個家系隨機選取300粒種子,每100粒種子作為1次重復,共重復3次,稱取質量后計算種子千粒重。將種子分別放入鋪有濾紙的培養皿中,加蒸餾水保持培養皿濕潤,然后置于30℃的恒溫培養箱中進行發芽試驗。當連續5 d沒有種子發芽時結束發芽試驗,參考顏啟傳等[16]的方法計算發芽率、發芽勢和發芽指數。
發芽率=發芽種子數/種子總數×100%
(1)
發芽勢=達到發芽高峰期累計的發芽數/種子總數×100%
(2)
發芽指數=∑第n天種子的發芽數/相應種子發芽所用天數
(3)
每個家系隨機選取30粒種子,每10粒種子作為1次重復,共重復3次,采用玻璃板直立法在30℃的恒溫培養箱中計算種子活力指數,與此同時,記錄幼苗發育狀況,將幼苗分為A、B、C、D四級:
活力指數=發芽指數×幼苗平均長度[17]
(4)

(5)

(6)

(7)

(8)

出苗率=(A級苗+B級苗+C級苗)/幼苗總數×100%
(9)
健壯率=(A級苗+B級苗)/幼苗總數×100%[12]
(10)
對獲得的百分比數據進行反正弦轉換后,利用SPSS16.0軟件進行方差分析、多重比較、相關分析和主成分分析,并計算家系遺傳力。
其中,各性狀的方差線性模型[18]為:
Xijk=μ+αi+βj(i)+Bk+αBik+βBj(i)k+eijk
(11)
式中:Xijk為第k重復第i采種地中第j家系觀測值;μ為試驗平均值;αi為第i采種地效應值;βj(i)為第i采種地中第j家系效應值;Bk為第k重復效應值;αBik為采種地與重復互作效應值;βBj(i)k為第i采種地中第j家系與區組互作效應值;eijk為機誤。
家系遺傳力采用公式[19]:
(12)

變異系數采用公式[20]:
(13)

相關系數采用公式[21]:
(14)

采用主成分分析法對家系進行綜合評定[22]:
Y=Σ(Ui·Zi)/n
(15)
式中:Y為主成分值;Ui為某一性狀特征向量;Zi為某一性狀標準分;n為主成分個數。
對參試的賽黑樺的3個采種地共30個半同胞家系種子的千粒重、發芽率、發芽勢、發芽指數、活力指數、出苗率和幼苗健壯率等7個性狀進行方差分析,集安和蛟河與琿春在幼苗健壯率方面差異顯著,在種子活力指數方面差異極顯著;7個性狀在不同家系間的差異均達到顯著或極顯著水平(表1)。通過計算家系遺傳力,發現除了出苗率和健壯率外,各性狀的家系遺傳力均高于0.68,說明種子千粒重和種子活力受環境影響較小、受遺傳控制程度較高。變異系數分析表明,千粒重的變異系數最小,僅為4.38%;發芽勢、發芽指數、活力指數和健壯率的變異系數均高于30%,其中活力指數的變異系數最大,高達45.22%,說明種子活力在家系間的變異十分豐富,具有很大的選擇潛力(表2)。
注:*表示在0.05水平差異顯著,**表示在0.01水平差異顯著。
Notes:*indicates a statistically significant at 0.05 level,**indicates a statistically significant at 0.01 level.

表2 不同家系種子千粒重及活力指標的遺傳力和主要統計參數


表4 賽黑樺種子各性狀的相關分析
注:*表示在0.05水平相關,**表示在0.01水平相關。
Notes:*indicates a statistically significant correlation at 0.05 level,**indicates a statistically significant correlation at 0.01 level.

表5 特征根及標準化特征向量
方差分析結果顯示,各性狀在不同家系間的差異均達到顯著或極顯著水平,因此對各家系進行多重比較。在種子千粒重方面,JH6表現最好,其千粒重達到426.7 mg,比各家系平均值高5.38%,比最差的JH8高13.48%;JH6、JA1和JA7的發芽率、發芽勢、發芽指數和活力指數表現較好,尤其在活力指數方面,JH6、JA1和JA7顯著高于其他家系,JH8的活力指數最低,僅為各家系平均值的41.12%;在出苗率和幼苗健壯率方面,JA1表現最好,其出苗率和健壯率分別高達100.00%和73.33%,分別比各家系平均值高24.32%和57.12%,HC9的出苗率和健壯率最差,分別僅為各家系平均值的78.73%和57.15%(表3)。
對30個半同胞家系的種子千粒重與活力指標進行相關分析,發現7個性狀間均呈顯著或極顯著的正相關(表4)。其中,發芽率與活力指數間的相關系數最大,高達0.991,二者之間為極顯著正相關;雖然千粒重與出苗率間的相關系數最小,僅為0.232,但二者之間的相關性也達到了顯著水平。上述結果說明種子質量與種子活力、幼苗的出苗率及健壯率均有相關性,種子千粒重越大,種子活力越強、幼苗越健壯。
上述分析就種子千粒重及活力指標等7個性狀進行了方差分析和多重比較。由于各性狀間存在顯著或極顯著的相關性,且優良家系的選擇需要綜合考慮多個性狀,因此采用主成分分析法求出特征根及各性狀的特征向量,通過計算主成分值(Y)對賽黑樺各家系進行綜合評價。特征根和特征向量的計算結果表明:只有1個主成分的特征根大于1,且各性狀的特征向量彼此相差不大,均介于0.304~0.408(表5),說明各性狀對家系綜合評價的影響較為接近。因此,可以只選用1個主成分進行家系綜合評價。主成分值的表達式為:Y=0.368X1+0.403X2+0.398X3+0.404X4+0.408X5+0.304X6+0.350X7。根據主成分值對各家系進行綜合評價,發現JH6的主成分值最大,達到4.555,JH8的主成分值最小,僅為-3.768,因此可以確定JH6為最優家系,JH8為最差家系。另外,JA1、JA7、JH10、JH3、HC1和HC10的各性狀標準分均大于0,且主成分值均在2以上,可以將其確定為優良家系(表6)。

表6 參試家系綜合評價
根據種子性狀評價采種母樹的優劣,對于優良家系的選擇十分重要。只有選擇種子活力高的家系,才能保證林木的產量和遺傳品質[23~24]。徐煥文等[12]通過比較白樺三倍體制種園中16個家系的種子活力,選出了3個三倍體優良家系。劉宇等[25]采用主成分分析法比較了白樺四倍體半同胞家系的種子活力和苗期生長性狀,選出了4個四倍體優良家系。本研究發現,雖然賽黑樺種子的千粒重、發芽率、發芽勢、發芽指數和出苗率在采種地間的差異不顯著,但在家系間的差異均達到了顯著或極顯著水平。JH6家系的種子千粒重、發芽率、發芽勢、發芽指數和活力指數分別比各家系平均值高5.38%、47.37%、45.80%、45.98%和86.78%,JA1家系的出苗率和健壯率分別比各家系平均值高24.32%和57.12%。這些結果說明根據種子千粒重和活力指標以及出苗率和幼苗健壯率進行采種母樹的初步選擇,可以為賽黑樺優良家系的選擇提供參考。
相關分析表明,賽黑樺種子千粒重與其他6個性狀均呈現顯著或極顯著正相關,這與白樺[26]、鹽松[27]、梓樹[28]的種子萌發試驗結果基本一致。對于小粒種子樹種而言,造成這種現象的原因可能是質量大的種子通常發育水平較高且含有較多的營養物質,進而具有較好的萌發能力。為提高育苗效率,可以對種子的千粒重、發芽率、發芽勢、活力指數、出苗率等性狀進行綜合利用,選擇高活力種子作為培育的對象,為生產優質苗木提供保障。由于種子千粒重的測定過程相對簡便快捷,因此如果選擇某一指標評價種子活力時,千粒重可以作為優先考慮的性狀。
由于賽黑樺各種子性狀間存在顯著或極顯著的正相關,在進行多性狀選擇時,很難確定各性狀所占的比重。主成分分析通過正交變換將一組存在相關性的變量轉換為一組線性不相關的變量,根據主成分值評價各家系綜合性狀的優劣,是一種比較理想的評價方法[22,25,29]。本研究通過主成分分析,將7個性狀轉化為1個主成分,根據主成分值選出JH6、JA1、JA7、JH10、JH3、HC1和HC10為優良家系,選出的7個優良家系分別來自蛟河、集安和琿春,說明每個采種地均可以選出優良家系。
上述研究僅對賽黑樺各半同胞家系的種子活力進行了比較,下一步我們將對這些家系進行育苗和造林,通過測定樹高、胸徑、保存率等性狀,進一步對母樹進行綜合評價,選擇種子活力高、子代生長迅速及適應性強采種母樹,為林業生產提供優良的賽黑樺種子。