思茅松無性系生長性狀早晚相關研究及早期選擇

李思廣，李宏斌，姜遠標，楊 斌

（1.云南省林業和草原科學院，云南昆明 650201；2.國家林業和草原局思茅松工程技術研究中心，云南思茅650201；3.普洱市職業教育中心，云南普洱 665000）

思茅松（Pinuskesiyavar.langbianensis）是云南省主要的材脂兼用樹種，其生態適應性強，生長迅速，產脂量較高，主要分布在云南省普洱市、紅河哈尼族彝族自治州、西雙版納傣族自治州景洪市和臨滄市的部分縣，分布海拔為600～1 600 m[1]。云南省思茅松林地總面積約為56.7 萬公頃，立木總蓄積量為6 072.9萬立方米[1]。中國是世界上最大的松香和松節油生產國。2016年，云南省松香產量為7.8萬噸，松節油產量約為3.5 萬噸，思茅松松節油產量占云南省松節油產量的90%左右。

現代林木育種的主要研究內容之一為縮短林木育種周期，加速育種進程。已有研究者開展一些樹種主要生長性狀遺傳規律及其相關性的研究。上世紀80年代以來，國內科研工作者對長白落葉松（Larixolgensis）、杉木（Cunninghamialanceolata）、蒙古櫟（Quercusmongolica）、油松（Pinustabuliformis）、火炬松（Pinustaeda）、馬尾松（Pinusmassoniana）和尾葉桉（Eucalypturophylla）等樹種開展早晚相關及早期選擇研究[2-10]，通過早期性狀對成熟期的經濟性狀作出判斷，提出最佳早期選擇林齡，達到早期選擇目的，使造林和營林不盲目進行。

思茅松遺傳改良研究主要集中在對松脂產量的改良上，初步揭示松脂產量的遺傳變異特征，并選出一批脂用思茅松優良家系和無性系[11-12]，但關于思茅松無性系早晚相關性的研究尚未開展。本研究以思茅松無性系為研究材料，對其樹高、胸徑和單株材積進行分析，研究早晚期相關關系，建立生長早期預測模型，研究早期選擇的效率和適宜的選擇林齡，為思茅松無性系早期選擇工作提供參考。

1 材料與方法

1.1 試驗地概況

試驗地位于云南省普洱市景谷縣的思茅松松脂基地（100°03′E，22°51′N），海拔1 520 m，年均降水量1 235 mm；年均氣溫21.1 ℃，最冷月（1月）平均氣溫13.0 ℃，最熱月（7月）平均氣溫24.6℃，≥10 ℃活動積溫7 360.9 ℃；5—10月為雨季；土壤類型主要為紅壤，西北坡向，平均坡度15°。

1.2 試驗設計

選擇40 株思茅松高產脂優樹。2001年3月整地，7月定植砧木；2002年5月，采集優樹穗條嫁接，營建無性系嫁接試驗林。采用隨機區組設計，6 次重復，6株為1小區，順坡排列，株行距為3 m×3 m[1]。2002—2016年（2004、2011、2013 和2015年除外）12月采用塔尺測量樹高，2005—2016年（2011、2013 和2015年除外）12月采用胸圍尺測量胸徑。

1.3 指標測定

單株材積（V，m3）計算公式[13]為：

式中，D為胸徑（cm）；H為樹高（m）。

方差分析線性模型[8]為：

式中，μ為總體平均值；αi為重復效應值；βi為無性系效應；εij為誤差。

遺傳變異系數（GCV，%）計算公式[8]為：

式中，δ2為遺傳方差；為性狀均值。

重復力（h2）計算公式[14]為：

式中，F為方差檢驗值。

為檢驗樹高、胸徑和單株材積的穩定性，根據其數據的排列位次，計算不同林齡思茅松無性系樹高、胸徑和單株材積與15年生思茅松無性系的秩次相關系數（rs），計算公式[15]為：

式中，di為第i對無性系位次的差值；n為無性系組成對數。

遺傳相關系數（rgjm）、表型相關系數（rpjm）和環境相關系數（rejm）計算公式[16]為：

式中，j、m為前期各林齡和第15年；σgm、σpm和σem分別為第15年遺傳方差平方根、表型方差平方根和環境方差平方根；σgj、σpj和σej分別為前期各林齡遺傳方差平方根、表型方差平方根和環境方差平方根；covgjm、covpjm和covejm分別為樹齡間遺傳協方差、表型協方差和環境協方差。

采用線性回歸模型，以早晚林齡比值（t/T）的自然對數為自變量，以相應的表型相關系數（rpjm）為因變量，建立生長性狀早晚相關系數估算模型[3]。

式中，A為常數；B為系數。

早期選擇效率（E）計算公式[17]為：

式中，Tm和Tj分別為輪伐期林齡和早期選擇林齡。

1.4 數據處理

采用SPSS、DPS和Excel軟件進行數據分析。

2 結果與分析

2.1 思茅松無性系生長性狀表型變異

樹高、胸徑和單株材積均值分別為0.37～11.48 m、5.12～14.08 cm 和0.004 1～0.101 3 m3；遺傳變異系數分別為9.7%～24.4%、14.0%～16.8%和33.2%～49.5%（表1）。樹高的遺傳變異系數隨樹齡增加而減少；胸徑和單株材積的遺傳變異系數均隨樹齡增加有波動。

2.2 思茅松無性系生長性狀方差分析及重復力

不同林齡思茅松無性系樹高、胸徑和單株材積均差異極顯著（P<0.01）（表2）。樹高、胸徑和單株材積的重復力均較高，分別為0.75～0.94、0.82～0.88和0.85～0.89。

表2 不同林齡思茅松無性系生長性狀方差分析及重復力Tab.2 Analysis of variance and repeatability forces of growth traits of P.kesiya var.langbianensis clones with different stand ages

2.3 思茅松無性系生長性狀的秩次相關分析

不同林齡思茅松無性系與15年生思茅松無性系樹高、胸徑和單株材積的秩次相關系數均呈極顯著水平（P<0.01）（表3）。隨樹齡增加，樹高、胸徑和單株材積的秩次相關系數均呈增加趨勢。6～13年生思茅松無性系與15年生思茅松無性系樹高、胸徑和單株材積的秩次相關系數分別為0.88～0.92、0.88～0.98和0.90～0.98，均較穩定。

表3 不同林齡思茅松無性系與15年生思茅松無性系生長性狀的秩次相關系數Tab.3 Rank correlation coefficients of growth traits among P.kesiya var.langbianensis clones with different stand ages and 15-year-old P.kesiya var.langbianensis clones

2.4 思茅松無性系生長性狀表型、遺傳和環境相關系數

不同林齡思茅松無性系與15年生思茅松無性系樹高、胸徑和單株材積的3 個相關系數均表現為遺傳相關系數>表型相關系數>環境相關系數（表4）。隨林齡增加，樹高、胸徑和單株材積的遺傳相關系數均呈增加趨勢。

表4 不同林齡思茅松無性系與15年生思茅松無性系生長性狀的表型、遺傳和環境相關系數Tab.4 Phenotypic,genetic and environmental correlation coefficients of growth traits among P.kesiya var.langbianensis clones with different stand ages and 15-year-old P.kesiya var.langbianensis clones

2.5 思茅松無性系樹高和胸徑早期預測模型及選擇效率

樹高和胸徑生長早晚相關系數估算模型方程為樹高RH早-H晚=1.009 9 + 0.220 8ln(x)，相關系數為0.984；胸徑RD早-D晚=1.021 4 + 0.161 8ln(x)，相關系數為0.928。兩個模型均存在極顯著線性相關，可用于估算早晚相關系數和早期選擇效率。

思茅松的主伐林齡為31年，更新期為5年，Tm值為36年。將樹高和胸徑早晚相關系數估算模型方程代入早期選擇效率公式，分別計算樹高和胸徑的早期選擇效率（表5）。

胸徑的早期選擇效率隨林齡增加而減少，第5年時達最高；第5年胸徑均值僅為5.12 cm，剛達起測徑階，不能進行早期選擇。第6年的選擇效率排名為第2，從第6年起胸徑的遺傳相關系數均達0.93以上，秩次相關系數均大于0.88；綜合考慮生長狀況、選擇效率、遺傳和秩次相關系數，胸徑適宜的早期選擇林齡為第6年。

樹高的早期選擇效率隨林齡增加呈先增后減的趨勢，第4年時達到最高。嫁接前5年因受嫁接和補接的影響，幼樹尚不能發揮其真正的生長特性。綜合考慮選擇效率、遺傳和秩次相關系數，樹高適宜的早期選擇林齡為第6年。

2.6 樹高、胸徑早期選擇正確率分析

以15年生思茅松無性系的生長性狀為目標性狀，對樹高和胸徑進行早期選擇驗證。入選率為5.0%時，樹高性狀的入選正確率為100%；入選率為12.5%時，入選正確率較低，僅為40%；入選率為25.0%時，入選正確率為70%。入選率為5.0%時，胸徑性狀的入選正確率為100%；入選率為12.5%時，入選正確率降為60%；入選率為25.0%時，入選正確率為80%（表6）。樹高和胸徑在6年生時，入選率為5.0%和25.0%的入選正確率較高，其早期選擇效果較為顯著。以胸徑性狀開展早期選擇的入選正確率高于樹高性狀。

表6 生長性狀早期選擇正確率Tab.6 Accuracy rates of early selection of growth traits

3 討論與結論

本研究中，樹高的遺傳變異系數總體隨林齡增加而減少，與賴猛[18]對落葉松和馬常耕等[4]對杉木的研究結果相似，原因可能是試驗林中的弱勢單株死亡，其他單株適應環境的能力不斷增強。

不同林齡思茅松無性系樹高、胸徑和單株材積均差異極顯著，說明其生長性狀變異豐富；重復力均較高，表明生長性狀均受遺傳控制，利用重復力進行良種選育可靠性較高。

生長性狀的年-年相關性可顯著提高間接選擇效率，是進行早期選擇的重要依據。不同林齡思茅松無性系各生長性狀的3個相關系數均表現為遺傳相關系數>表型相關系數>環境相關系數，說明生長性狀主要受遺傳控制。隨林齡增加，樹高、胸徑和單株材積的遺傳相關系數均呈增加趨勢，表明樹體越大，各生長性狀表現越接近成熟期。這與賴猛[18]對落葉松及馬常耕等[4]和胡德活等[14]對杉木的研究結果相似。隨林齡增加，樹高、胸徑和單株材積的秩次相關系數均呈增加趨勢；不同林齡思茅松無性系與15年生思茅松無性系的秩次相關系數均達極顯著水平；說明可以進行無性系早期選擇。

依據樹高和胸徑的早期選擇效率，結合重復力、遺傳和秩次相關系數，確定樹高和胸徑較適宜的早期選擇林齡為第6年。大部分無性系生長排序較穩定，但一部分無性系生長較為特殊，被誤選或漏選。59 號無性系的胸徑、樹高生長不一致，按胸徑排名（第2），可入選；按樹高排名（中游），會漏選。47 號無性系前期和后期生長速度不一致，按前期樹高、胸徑排名（均排名20 位后），會漏選；按后期樹高、胸徑生長情況，可入選。89 號無性系前期生長快，后期生長變慢；6年生時其樹高排名第3，胸徑排名第6；15年生時其樹高排名第10，胸徑排名第6，會造成誤選。進行思茅松無性系早期選擇時，應綜合考慮樹高和胸徑，才能提高選擇準確率。

利益沖突：所有作者聲明無利益沖突。

作者貢獻聲明：李思廣負責研究計劃制定、試驗調查、數據收集與分析、論文撰寫與修改和文獻檢索；李宏斌負責研究計劃制定、論文修改；姜遠標負責研究計劃制定、數據收集和文獻檢索；楊斌負責試驗調查、數據收集。