红锥Castanopsis hystrix又名红黎、赤黎、黎木、刺栲、红栲、红掾栲，为壳斗科栲属常绿乔木，是华南地区重要的乡土阔叶珍贵用材和高效多用途树种^[1-3]。自20世纪70年代末以来，广西、广东、福建等省区已开展红锥天然林资源的调查、材性试验、良种选育、栽培试验、混交林技术、生物学特性、木材利用技术以及病虫害防治等方面的研究^[4-6]。我国阔叶树种子园各项技术研究虽起步较晚，但可借鉴杉木Cunninghamia spp.和松树Pinus spp.等的种子园营建及经营技术上的研究成果^[7]，包括园址选择、亲本开花结实特性、控制授粉、树体管理、土壤管理等方面^[8]。本文以广东省龙眼洞林场2代红锥种子园为研究对象，对种子园内的56个红锥无性系的生长性状树高、胸径和冠幅，以及形质性状尖削度、通直度、分枝角度、分枝数、分枝大小、叶片疏密和树形进行调查，比较红锥生长性状和形质性状的差异性及其相关性，采用主成分分析法和聚类分析法对不同无性系的生长性状和形质性状进行综合评价，旨在为红锥常规杂交育种研究和优异基因的挖掘提供参考。

1 材料与方法 1.1 调查材料

试验地设在广东省龙眼洞林场墓园工业区内(东经113°21′~113°27′，北纬23°10′~23°18′)。该地年降雨量1 694 mm，年均温21.8 ℃，海拔120 m。种子园林地为深厚的壤质赤红壤，土层厚40~60 cm，pH4.8~5.2。该园为2008年嫁接的2代种子园，56个无性系均为广西、广东和福建种源的种子经过子代测定选择的优良个体。种子园试验采用完全随机区组设计，7次重复，重复内的家系采用不相邻小区种植，各重复内家系植株数量为3~20株，株行距5.0 m×5.0 m。2015年12月对红锥56个无性系1 030个单株进行了树高、胸径、冠幅、通直度、分枝角度和分枝数的测量，具体测定方法见文献[9]。此外，还有4个指标：尖削度=树干直径(1/2树高位置)/树干直径(1.3 m 树高位置)；分枝大小，总体侧枝粗＜1.5 cm为1级，总体侧枝粗[1.5, 2.0) cm为2级，总体侧枝粗≥2.0 cm为3级；叶片疏密，疏松计为1，紧密计为2；树形，分为圆满型、三角形和伞形3种类型，依次计为1、2 和 3。

1.2 数据统计分析

利用Excel对原始数据进行整理，采用R软件“psych”包进行主成分分析和样品聚类分析，其中聚类分析采用系统聚类法。应用SAS软件进行不同红锥无性系间各生长性状和形质指标的方差分析、相关性分析。其中，调用GLM过程进行方差分析并选择TYPEⅢ的平方和类型^[10]。红锥种子园生长和形质性状观测数据按下列模型进行方差分析：

${Y_{i\!jk}} = u + {B_i} + {C_j} + {e_{i\!jk}},$

式中，Y_ijk为无性系的性状观测值，u为总体平均值，B_i为区组效应，C_j为无性系效应，e_ijk为剩余误差，并计算参数：无性系重复力(R²)，R²=1–1/F，其中，F为方差分析TYPEⅢ的F值；变异系数(CV)，CV=SD/M，其中，SD和M分别表示表型标准差和均值；遗传相关系数(r)， $r={\rm{cov}}XY/\sqrt {\sigma _X^2 \sigma _Y^2} $ ，其中cov为协方差，σ²为方差。

2 结果与分析 2.1 2代红锥种子园生长和形质性状变异分析

对56个红锥无性系生长和形质的10个性状指标的变异分析结果见表1。红锥平均树高 7.8 m、胸径 9.1 cm、尖削度 0.80、冠幅 3.90 m、通直度 3.9、分枝角度 56°、分枝数 5.8个、分枝大小为 2.20级、叶片疏密为 1.4、树形为 1.60。除尖削度的变异系数较小外，其他均高于15%，说明这些性状个体间的差异较大，性状表现不稳定，因此在建立红锥高世代种子园时有必要对建园材料进行再选择。红锥的生长和形质性状变异系数排序为：分枝数>树形>叶片疏密>分枝大小>通直度>分枝角度>胸径>冠幅>树高>尖削度。

2.2 2代红锥种子园生长和形质性状方差分析

按单株小区试验设计进行方差分析，分别对56个红锥无性系生长和形质性状的差异显著性进行检测。由表2可知，除尖削度的区组效应未能达到显著水平外，其他性状的区组效应均达显著或极显著水平。不同性状的无性系效应检测F值介于1.27~2.14，除分枝角度外，其他性状的F值显著性水平均达到显著或极显著。不同性状的重复力为0.21~0.53，表明这些性状在单株水平上受到中等水平的遗传控制，进一步开展无性系选择潜力较大。

2.3 2代红锥种子园生长和形质性状的相关性分析

对不同红锥无性系生长和形质性状的相关性分析结果(表3)表明，不同红锥无性系树高与胸径、冠幅、分枝数之间存在极显著正相关性(P<0.001)；胸径与冠幅、分枝数、分枝大小和叶片疏密之间存在极显著正相关性(P<0.001)；叶片疏密与分枝角度、分枝数和分枝大小间存在明显的正相关性(P<0.05或P<0.01)。冠幅与分枝角度，尖削度与分枝大小间存在显著的负相关(P<0.001或P<0.05)，但相关值较低, 介于–0.13~ –0.28。可见红锥无性系生长量对分枝数、分枝大小和叶片疏密具有较大的影响，但对尖削度、通直度、分枝角度和树形影响不大。

2.4 2代红锥种子园生长和形质性状的主成分分析

主成分分析结果(表4)表明：第1主成分的特征值为3.18，方差贡献率为44%，代表了全部信息的44%。第2主成分的特征值为1.77，方差贡献率为28%，代表了全部信息的28%。第3主成分的特征值为1.29，方差贡献率为13%，代表了全部信息的13%。前3个主成分的特征值均大于1，且累积方差贡献率为85%，包含了全部指标的大部分信息，因此可以选取前3个主成分作为红锥生长和形质性状选择的综合指标。

在第1主成分中，胸径、树高、冠幅占有较高载荷，说明第1主成分反映的是与生长量有关的指标信息，可称为生长量因子。在第2主成分中，尖削度、通直度占有较高的载荷，说明第2主成分反映了与树形有关的信息，可称为树形因子。在第3主成分中，分枝角度占有较高的载荷，反映的是与侧枝有关的指标信息，可称为侧枝因子。由此可知，红锥无性系之间的差异主要表现在生长量、树形和侧枝。

2.5 2代红锥种子园生长和形质性状的聚类分析

以所测的10个生长和形质性状的平均值为参数，采用系统聚类法对56个红锥无性系进行聚类分析，结果见图1。10个生长和形质性状采用系统聚类能够将56个无性系完全区分开，说明所研究的生长和形质性状能够反映各个样品间的差异。取阈值为15时，56份红锥无性系可分为4组：第1组包含8份材料，约占总样本数的14.3%；第2组包含4份材料，约占总样本数的7.1%；第3组包含8份材料，约占总样本数的14.3%；第4组包含36份材料，分为3个亚组，约占总样本数的64.3%。由此可见，56个红锥无性系具有各自的生长和形质表型特征，这为开展红锥种内杂交育种的亲本选择奠定基础。

3 讨论与结论

通过对2代红锥种子园56份无性系单株生长和形质性状调查表明，10个生长和形质性状的变异系数介于11%~47%，方差分析也表明各无性系不同性状间差异达显著或极显著水平，不同性状的重复力为0.21~0.53。大量研究表明林木种子园不同无性系间的生长性状和形质性状间存在显著差异^[11-12]。廖焕琴等^[9]发现红锥1.5代种子园无性系间的不同性状表现出显著或极显著差异。陈苏英等^[13]的研究发现1.5代杉木种子园不同无性系的生长性状(树高、胸径、单株材积)间存在显著差异。孙文生^[14]的研究发现红松Pinus koraiensis种子园各无性系间的结实量存在极显著差异，立地条件的不同、管理方式差异、树龄及不同无性系对种子园结实量具有显著影响。因此，有必要对红锥种子园内各无性系植株进行调查与清理，淘汰不适合本地生长的无性系植株。

各生长和形质性状间存在不同程度的相关性。树高与胸径、冠幅与树高、胸径间极显著正相关，相关系数分别为0.51、0.60和0.63。胸径与分枝数、分枝大小和叶片疏密间也存在极显著正相关，相关系数介于0.51~0.60。这与廖焕琴等^[9]的研究结果一致。冠幅与分枝角度、尖削度与分枝大小存在显著的负相关，相关值较低。可见，对相关性系数较高的性状进行选择时，对一个性状的改良，可能同时影响另一个表型性状。因此，在对生长和形质等综合性状进行改良时，应关注性状间的正负相关性并充分利用这些相关性提高林木育种效率。

对56份红锥无性系生长和形质性状的主成分分析结果表明，前3个主成分累积贡献率达85%，因此所选的前3个性状具有较强的代表性。大量研究表明主成分分析结果与参试资源和性状指标的选择有关系^[15]。从对56份红锥无性系的主成分分析可以看出，不同性状在3个主成分中具有明显不同的载荷值，载荷值的绝对值越大，说明该因子对当前变量的影响程度越大^[16]，因此可以将9个生长和形质性状归类为：第1主成分以胸径性状为主，第2主成分以尖削度为主，第3主成分以分枝角度性状为主。由此可知，红锥无性系之间的差异主要表现在生长量、尖削度和分枝角度。