无性系选育作为现代良种选育的主要方式,在提升人工林遗传质量和经营效率方面具有显著优势。早期的无性系选育主要侧重树高、胸径和单株材积等生长性状,而对主干材高、通直度等与木材品质和生态适应性密切相关的性状关注不足
[1]。随着木材加工业的转型升级以及生态林体系建设的推进,单一性状已难以全面评估和反映无性系的综合育种价值。林木育种逐步由以单一经济性状为目标的选育模式,向兼顾高生长性和优良干形结构的多性状协同改良方向发展
[2]。无性系的系统性评价和综合选择,对于构建优质高产的良种体系,提升人工林育种效率与质量,具有重要的理论意义与实践价值
[1]。目前已对杉木(
Cunninghamia lanceolata)、马尾松(
Pinus massoniana)、白皮松(
Pinus)和西南桦(
Betula alnoides)
[3-6]等开展了无性系多性状综合评价与优选研究,筛选出一批优良无性系,显著提升了人工林的产量和质量。
柚木(
Tectona grandis L.f.)属落叶半落叶高大乔木,材质优异,被广泛应用于造船、建筑、雕刻及高档家具制造领域,是世界著名的珍贵用材树种
[7-8]。我国自2000年以来已开展了3批次柚木无性系区域试验。黄桂华等
[9]对云南景谷柚木23个无性系生长数据进行分析评价,从中筛选出7029、7531、Z408和8301 4个优良无性系。潘栋康等
[10]对广西凭祥柚木19个无性系开展综合评价,筛选出3个兼具速生性与良好分枝特性的无性系。然而始终缺乏柚木多无性系、多指标、长林龄综合选择与评价。本研究基于定植于海南省定安县的14年生25个柚木无性系,开展生长性状和形质性状综合评价研究,旨在筛选出多性状优良的柚木无性系,为柚木良种的繁育与推广提供参考。
1 材料与方法
1.1 试验地概况
试验地位于海南省定安县海南省林业科学研究院龙州基地(110°15' 17″ E,19°38' 29″ N),属热带海洋性季风气候,年平均气温24 ℃,年平均降水量2 000 mm。试验地地势平坦,土壤为褐色砖红壤,pH值5.3。
1.2 试验设计
参试无性系25个,采用6株行式小区,随机完全区组排列,6次重复。穴规格60 cm×60 cm×40 cm,穴施基肥1.0 kg钙镁磷肥+1.0 kg石灰。株行距2.5 m ×4 m,2010年7月造林。
1.3 测定指标与方法
于2024年10月对柚木无性系林各指标进行测定。
1.3.1 生长性状
每木检尺,测量树高(H)、胸径(DBH)和主干材高(clear bole height, CBH)。采用测高杆测量树高和主干材高,精度为0.1 m;采用胸径尺测量胸径,精度0.1 cm。保存率(survival rate)计算公式为:,其中为每小区调查存活总数。
单株材积(V)计算公式为:V=0.478 7×DBH2×H。
1.3.2 形质性状
树干通直度、干材位置、分枝模式与主干分枝采用分级评价法
[11-12],具体指标及得分评估标准如下。
(1) 通直度(straightness,STR)。根据树干主轴在树高方向的直立程度,通直度分为5个等级:树干从基部有大分叉或有2段以上明显弯曲,得1分;树干有两段以上稍弯曲,得2分;树干有1~2段稍弯曲,得3分;树干无明显弯曲,得4分;树干通直,无弯曲,得5分。
(2) 干材位置(bole position,BP)。从下到上将树的总高度分为4个相等的部分,然后据此判断干材位置,可分为4个等级:干材位置在树的第1个1/4处,得1分;干材位置在树的第2个1/4处,得2分;干材位置在树的第3个1/4处,得3分;干材位置在树的第4个1/4处,得4分。
(3) 分枝模式(mode of branching,MB)。根据树木分枝的特征,将分枝模式分为5个等级:树干从下半部分开始分枝,得1分;树干从上半部分开始分枝,得2分;树干无大分枝且侧枝扩散不规则,得3分;树干无大分枝且侧枝扩散较规则,得4分;树干无大分枝且侧枝扩散规则,得5分。
(4) 主干分枝(main stem branching,MSB)。根据主干上小枝的数量和分布情况,将主干分枝分为4个等级:主干小枝极密,均匀分布于整个树干,得1分;主干小枝较多,约50%的部位有分布,得2分;主干小枝较少,约25%的部位有分布,得3分;主干几乎无小枝,得4分。
为了减少偏见,以上指标的定性分数均由同一位评估者进行评分并记录。
1.4 数据统计与分析
利用Excel软件整理数据,用SPSS 22软件进行方差分析和邓肯多重比较。采用R语言sommer包拟合多性状混合线性模型
[13-14],计算遗传相关系数。方差分析线性模型、表型变异系数、重复力、遗传相关系数、表型相关系数和遗传增益等参数的计算参照文献[
15-
16]进行。
2 结果与分析
2.1 柚木无性系各性状遗传变异分析
保存率大小反映无性系对环境的适应性水平。
表1表明,不同柚木无性系间保存率差异明显。参试的25个无性系的平均保存率为80.77%,其中15个无性系保存率超过平均值。无性系8301和7549表现优良,保存率分别高达97.22%和94.44%;而无性系J730的保存率最低,仅为47.22%。
表1表明,不同柚木无性系间的树高、胸径和单株材积差异较大。25个柚木无性系的树高为8.54~12.46 m,平均树高为10.85 m,其中无性系7514和7509树高较大,分别达12.46和12.43 m,较平均值分别增加了14.84%和14.56%;而无性系7112的树高最低(8.54 m),较平均值降低了21.29%。25个柚木无性系的胸径为10.67~17.47 cm,平均胸径为14.02 cm,其中无性系7122和7559胸径较高,分别达17.47和16.24 cm,较平均值提高了24.60%和15.83%;无性系7531的胸径最低(10.67 cm),较平均值降低了26.03%。25个柚木无性系的单株材积为0.05~0.17 m³,平均单株材积为0.11 m³,其中无性系7122、7514和7029的单株材积较高,分别为0.17,0.16和0.16 m³,较平均值提高了45.5%~54.5%;无性系7531和7112的单株材积均较低(均为0.05 m³),均较平均值降低了54.5%。25个柚木无性系的主干材高为3.94~7.85 m,平均主干材高为5.70 m,其中无性系7514和7219主干材高较大,分别达7.85和7.30 m,较平均值分别增加了37.72%和28.07%;而无性系7112的主干材高最低(3.94 m),较平均值降低了30.88%。
表2表明,各无性系干材位置、通直度、分枝模式和主干分枝等指标均存在明显差异。25个柚木无性系干材位置得分为2.38~3.07,平均干材位置得分为2.66,其中无性系7514、7146和7219干材位置得分较高,分别为3.07,3.00和3.00,较平均值提高了12.78%~15.41%;无性系7112干材位置得分最低(2.38),较平均值降低了10.53%。25个柚木无性系通直度得分为2.16~3.06,平均通直度得分为2.57,其中无性系7029和3078-5通直度得分较高,分别达3.06和2.97,较平均值分别增加了19.07%和15.56%;而无性系Z408通直度得分最低(2.16),较平均值降低了15.95%。25个柚木无性系分枝模式得分为2.00~2.89,平均分枝模式得分为2.49,其中无性系7219和71-5分枝模式得分较高,分别达2.89和2.83,较平均值分别增加了16.06%和13.65%;而无性系J731分枝模式得分最低(2.00),较平均值降低了19.68%。25个柚木无性系主干分枝得分为2.34~2.93,平均主干分枝得分为2.74,其中无性系7122和3078-5主干分枝得分较高,分别达2.93和2.93,均较平均值增加了6.93%;而无性系7137主干分枝得分最低(2.34),较平均值降低了14.60%。
2.2 柚木无性系遗传参数分析
表3表明,不同柚木无性系在生长性状和形质性状方面表现出显著的遗传多样性。
2.2.1 生长性状变异
由
表3可以看出,树高、胸径和单株材积生长性状表型变异系数分别为19.53%,22.79%和52.89%;遗传变异系数依次为9.19%,14.10%和32.11%。两个变异系数均呈现单株材积>胸径>树高的趋势,表明参试无性系具有丰富的生长性状变异。在遗传稳定性上,树高、胸径和单株材积的重复力均在0.80以上,具体为胸径(0.94)=单株材积(0.94)>树高(0.88),表明生长性状主要受遗传因素控制,且各无性系内个体生长表现高度一致。方差分析进一步表明,参试的25个无性系之间、区组之间以及无性系×区组互作效应对生长性状的影响总体均达到极显著水平(
P<0.001),其中无性系间的变异贡献率最高,凸显了遗传因素在柚木生长性状调控中的主导作用。
2.2.2 形质性状变异
由
表3可以看出,参试无性系形质性状表现出丰富的表型变异,各性状的表型变异系数由大到小依次为主干材高(43.46%)>干材位置(32.46%)>分枝模式(29.58%)>通直度(29.00%)>主干分枝(18.00%),可知主干材高变异幅度最大,而主干分枝相对较低。遗传稳定性分析表明,主干材高的重复力最高;除干材位置外,其余性状的重复力均在0.50以上。各形质性状重复力由大到小依次为主干材高(0.78)>通直度(0.71)>分枝模式(0.69)>主干分枝(0.61)>干材位置(0.46)。方差分析结果表明,参试柚木无性系在形质性状上存在显著或极显著差异。对于通直度、分枝模式和主干分枝3个指标而言,各效应的影响程度依次为区组效应>无性系效应>无性系×区组互作效应,其中区组效应对这些性状的影响最为突出,表明环境因素在决定这些性状表现中起主导作用。主干材高和干材位置同时受到无性系效应和区组效应的显著或极显著影响。
2.3 柚木无性系各性状间相关性分析
图1表明,多数性状间在表型和遗传方面均存在显著或极显著相关性,这些性状可能存在共同的遗传基础。在生长性状方面,各指标间表现出较强的相关性,其中树高与胸径在表型水平呈中度正相关(
rp=0.57),而在遗传水平上相关性较高(
rg=0.78),表明通过遗传改良提升树高可有效促进胸径增长。单株材积作为综合指标,与树高(
rp=0.71,
rg=0.85)和胸径(
rp=0.90,
rg=0.96)均呈极显著正相关,显示出优异的遗传潜力。表型水平上,树高、胸径和单株材积与除干材位置外的所有形质性状均呈极显著相关。主干材高与其余各形质指标存在显著表型相关性,表明形质性状之间存在协同发育机制,通过选育高主干材高的无性系可同步改良其他形质性状。
2.4 柚木无性系综合评价与筛选
对25个无性系的8个性状(树高、胸径、单株材积、主干材高、干材位置、通直度、分枝模式和主干分枝)进行主成分分析,并计算各主成分得分,结果显示,第一主成分的方差贡献率达65.24%,其主要反映柚木无性系整体生长量,得分较高的无性系如7514和7219在多个生长及形质性状上表现优异。第二主成分方差贡献率为14.01%,主要反映柚木无性系生长与形质性状的平衡特征。前2个主成分的累计方差贡献率为79.18%,以前2个主成分的方差贡献率计算综合得分,结果(
表4)显示,在兼顾生长和形质性状的基础上,采用15%作为优良无性系入选率,最终筛选出了4个无性系,分别为7514、7219、3078-5和7122。
入选柚木无性系不同性状的遗传增益如
表5所示。
表5表明,入选的优良无性系在生长性状和形质性状方面整体表现优异,且表现出一定的遗传改良潜力,其中7514在多个性状上具有较高的遗传增益,特别是在单株材积(49.70%)和主干材高(37.93%)指标上表现突出,在生长与形质性状方面具备较强的育种潜力。
3 讨 论
本研究发现,柚木不同无性系在保存率、生长性状和形质性状上均存在显著差异和丰富的遗传变异,且多数性状间存在表型和遗传相关性。保存率是衡量林木无性系适应性与稳定性的重要指标
[17]。本研究中,海南定安14年生25个柚木无性系测定林的保存率变幅在47.22%~97.22%,明显高于梁坤南等
[18]对滇西南地区29年生柚木种源/家系的保存率。这可能源于两地生态条件、种植密度以及种植材料遗传背景等方面的差异,也反映出基因型与环境条件之间存在的密切互作关系。此外,本研究发现柚木25个无性系之间保存率差异明显,说明无性系对环境的适应能力存在较强的遗传分化。因此,在开展无性系试验或推广造林时,应重视保存率的评估,以确保材料在特定环境中的成活率及稳定性。
研究表明,不同树种的无性系群体普遍存在极显著的遗传变异,但变异模式常表现出明显的树种特异性
[19]。研究发现,柚木无性系在树高、胸径等性状上具有高遗传力
[20]。本研究中,柚木无性系在生长性状(树高、胸径和单株材积)和形质性状(主干通直度和分枝模式)方面均表现出显著的遗传分化,揭示了柚木种质资源丰富的基因型多样性及其巨大的遗传改良潜力,这与前人的研究结果
[20]一致。本研究中,无性系7029、3078-5、7514和7122在多个性状指标上表现优异,其中7029的优良特性与黄桂华等
[9]对11年生柚木无性系的研究结果一致,表明该无性系具有稳定的遗传优势。
遗传变异是选择的基础,只有充分了解林木的遗传变异规律,才能提高林木改良的准确性
[21]。本研究发现,柚木单株材积的表型变异系数(52.89%)和重复力(0.94)均最高,这与李运兴等
[22]对柚木家系评价研究结论一致,验证了单株材积为林木遗传改良的重点性状。本研究的表型与遗传相关系数分析结果显示,单株材积与树高、胸径等生长性状呈极显著正相关,可综合反映两者的生长特性。单株材积与主干材高和通直度等形质性状的遗传相关性也较高,表明优选形质性状优良的无性系有望同步提升其生长潜力
[23-25]。
主成分分析法将不同无性系的性状表现置于统一标准下进行对比分析。本研究筛选出的4个优良无性系(7514、7219、3078-5和7122)不仅在单个性状上表现突出,且在多性状间发展协调,为柚木无性系在海南的大面积推广利用提供高增益遗传材料支持。优良无性系各性状的遗传增益中,单株材积的遗传增益(21.84%~60.07%)最高,凸显了单株材积在多性状综合选择中的主导地位,反映出其在衡量无性系生长潜力与选育效益方面的重要作用
[4,26-27]。本研究仅对海南单一地区柚木进行分析和评价,筛选结果具有一定局限性,尚需要结合多点试验进行全面分析
[28],以完善研究结论。
4 结 论
海南省定安县柚木无性系的生长和形质性状具有丰富的变异,综合评价后从中筛选出7514、7219、3078-5和7122 4个无性系,其保存率高,在多个目标性状上表现出较大遗传与现实增益,综合表现优异,为海南省柚木良种繁育与推广种植提供了重要的种质材料。
“十四五”国家重点研发计划项目(2024YFD2200102)