一种杉木优良个体的选择方法

摘要

本发明公开了一种杉木优良个体的选择方法，利用基于空间分析方法估计空间残差方差与空间自相关、限制性极大似然法估计遗传方差和最佳线性无偏预测杉木个体育种值的混合线性模型方法，分析杉木主要经济性状包括树高、胸径、木材基本密度和微纤丝角等的遗传变异、遗传控制和性状间的遗传相关，预测杉木每个个体每个性状的育种值，并根据杉木的多性状指数选择公式选出杉木生长与材性在遗传上兼优的优良个体。该方法，在遗传参数估计时剔除了环境影响的同时，基于育种值的多性状指数选择能有效选出综合性状在遗传上表现优良的个体，准确性高，所选出的优良个体可以作为杉木无性快繁的材料，也可以作为杉木下一代育种群体的育种材料。

说明书

技术领域

本发明涉及林木育种技术领域，特别涉及一种杉木优良个体的选择方法。

背景技术

目前，杉木优良单株常按表现型进行选择的，而表现型是基因型和环境相互作用的 结果，良好的表现型取决于由环境引起的位置效应和遗传因素，它能不能传递到后代以 及有多大的强度是未知的，因此单凭表现性的优劣来决定后代的好坏和建立高世代种子 园是十分困难的。在杉木高世代育种中，目前主要以表型选择为主，不能真正选出遗传 上优良的个体，最佳线性无偏预测的最重要优点是所有来源于多世代遗传测定的数据都 能被整合到一个分析中，在所有世代中最佳线性无偏预测的育种值被调整后来解释在随 后的育种循环中逐渐增加的遗传平均值，并且是在同一个水平上来进行预测，这样就可 以不考虑候选个体的世代，在考虑亲缘关系的同时可只根据最佳线性无偏预测育种值来 进行选择，即在第1个育种循环中通过表型选择表现优良的个体也可能包括在第2、3 和4个育种循环的育种群体中，子代中表现优良的个体也可以被选择到育种群体中。在 选择指数中，2个或2个以上具有很高正相关性状只要取其中的一个进行计算即可。因 为对于一个性状的选择将对另外具有很高遗传相关的性状产生很强的间接增益，在这些 性状上分配较高的权重系数在数学上没有任何意义，并且把2个或2个以上具有很高正 相关性状计算得到的选择指数通常会产生很小的遗传价值的增益。根据育种值和选择指 数得到的优良个体往往会集中在一些表现好的家系上，这使得遗传多样性大大减小。尤 其对于遗传力介于0.1至0.2的这些性状，在决定遗传价值上来源于子代测定的家系表 现是最有用的，因为遗传力低的性状由于个体树之间的变异受环境影响较大，个体树表 现通常给予一个较小的权重来预测它的育种值。

发明内容

发明目的：针对现有技术中存在的不足，本发明的目的就是提供一种杉木优良个体 的选择方法，为筛选生长量大、木材品质表现稳定和抗逆性强的杉木高世代种子园中的 优良家系和子代个体，挖掘良种潜力提供可行的方法。

技术方案：为了实现上述发明目的，本发明采用的技术方案为：

一种杉木优良个体的选择方法，包括以下步骤：

(1)测定杉木经济性状；

(2)利用基于空间分析方法估计空间残差方差与空间自相关、限制性极大似然法 估计遗传方差和最佳线性无偏预测杉木个体育种值的混合线性模型方法，分析杉木经济 性状的遗传变异、遗传控制和性状间的遗传相关，预测杉木每个个体每个性状的育种值；

(3)根据杉木的多性状指数选择公式，选出杉木生长与材性在遗传上兼优的优良 个体；多性状指数选择公式如下：

$I_i=0.7{\hat{a}}_{i1}+0.3{\hat{a}}_{i2}$

式中，I_i是基因型i的选择指数值；0.7和0.3分别是杉木生长性状和材性性状的权 重；和是基因型i每个性状的最佳线性无偏预测育种值。

步骤(1)中，所述的主要经济性状包括树高、胸径、木材基本密度、木材微纤丝 角度、木材红心材率和偏心率。

步骤(2)中，具体过程如下：

(1)数据整理与标准化，使数据符合正态分布；

(2)空间分析；空间分析模型如下：

Y＝Xb+Zf+ξ+η

其中，Y是数据向量，固定效应向量为均值向量，b和f是随机效应向量，X和Z 是观测值与随机效应相对应的设计矩阵，ξ是空间随机误差自变量向量，η是随机残差 向量；

(3)遗传参数估计，利用限制性极大似然法估计方差分量，并根据各方差组分计 算遗传力、遗传相关；

(4)育种值计算，建立亲缘关系矩阵，根据计算得到的家系方差和环境方差计算 育种值，公式为：

$\left(\begin{array}{c}\hat{b}\\ \hat{f}\end{array}\right)={\left(\begin{array}{cc}{X}^{\prime }X& X^{\prime }Z\\ Z^{\prime }X& Z^{\prime }Z+A^{-1}\frac{{\sigma }_e^2}{{\sigma }_f^2}\end{array}\right)}^{-1}\left(\begin{array}{c}{X}^{\prime }y\\ Z^{\prime }y\end{array}\right)$

其中y为性状表型观测值构成的向量；b为固定效应构成的向量，且第一个元素为 群体均值；f为随机效应构成的向量；e为剩余误差向量；X为固定效应的关联矩阵；Z 为随机效应的关联矩阵；假设随机误差和随机效应向量独立且服从相同的分布，即

E(f)＝0，E(e)＝O，E(y)＝Xb，Var(y)＝ZGZ′+R，Cov(f，e)＝0

，其中R和G分别为误差方差协方差矩阵和随机效应方差协方差矩阵，A为分子亲缘 关系矩阵，为误差方差，随机效应方差。

步骤(3)中，选择I_i指数值从高到低排序并以1%的选择强度选出优良个体用于无 性繁殖。

有益效果：与现有技术相比，本发明的杉木优良个体的选择方法，在遗传参数估计 时剔除了环境影响的同时，基于育种值的多性状指数选择能有效选出综合性状在遗传上 表现优良的个体，准确性高，为筛选生长量大、木材品质表现稳定和抗逆性强的杉木高 世代种子园中的优良家系和子代个体，挖掘良种潜力提供可行的方法。所选出的优良个 体可以作为杉木无性快繁的材料，也可以作为杉木下一代育种群体的育种材料。具有很 好的实用性。

附图说明

图1是木材基本密度的空间残差方差空间趋势图；

图2是微纤丝角的空间残差方差空间趋势图；

图3是胸径的空间残差方差空间趋势图；

图4是树高的空间残差方差空间趋势图。

具体实施方式

下面结合具体实施方式对本发明作进一步的说明。

实施例1

材料来源：福建省沙县官庄国有林场杉木第2代种子园位于池村工区石景山40林 班54大班2小班，面积195亩。1982年2月定砧，1984年3月至4月1嫁接。种子园 前茬为天然林，经砍伐、炼山清杂，整个种子园共分17个组34个小区，其中7、14、 26、27小区分两个组，区号到30小区，256个无性系，采用顺序错位设计配置方式。 其中的256个家系分别来自于种源试验林、福建当地选择的优树和南平地区早期选择优 树的杂交子代。

本发明用的苗木全部由福建省沙县官庄国有林场提供，参试家系80个，编号从1 到80，其中参试家系亲本包括50个半同胞家系和30个全同胞家系，以官庄国有林场2 代种子园混合种子育的苗木为共同对照，编号为81，共有81个处理。种子于2001年 10月采自福建省沙县官庄国有林场杉木第2代种子园，其中1到80号为杉木第2代种 子园家系代码，81号为官庄国有林场杉木第2代种子园混合代码，2002年春在福建省 沙县官庄国有林场场部附近苗圃分家系育苗。

试验点选择在立地类型和气候条件不同的福建省沙县官庄国有林场、福建省邵武卫 闽国有林场、福建省仙游溪口国有林场和福建省国有来舟林业试验场等4个试验点，每 个试验点约50亩。试验林地选择要求具有当地代表性的山场，地形变化小，开阔平坦， 坡面较大，交通方便但不易遭人畜破坏的山地。

造林试验采用平衡格子设计(balanced lattice)，即使每一对家系都能在同一区组中 相遇一次，这种设计的优点是各对家系间比较的精确性相对一致，分析方法相对简单。 参试家系81个，处理数81个，重复10次，区组数90个，区组内设置9个小区，采用 4株单列小区，从上坡到下坡顺坡排列。

区组沿山坡水平方向横坡排列，区组内各小区顺坡单行排列，区组内的立地条件要 求一致，重复内的立地条件基本一致，9个区组连在一起构成一个重复，区组之间可以 相连，也可以不相连，视山场变化情况而定，区组内小区一定相连。区组在山上的排列 顺序面向山场遵循“由上到下、由左往右”的原则，同一重复内的9个区组连在一起。 试验区外围设置普通杉木为保护行。

树高、树皮厚度、冠幅和胸径测量方法：苗期用游标卡尺测量地径，精确至0.1cm； 每木测定树高和胸径，精确至0.1m和0.1cm；冠幅测量树木的两个垂直方向上的宽度取 平均值，精确至0.1m；离地1.3m的胸径处于每木上取通过髓心的完整木芯，树皮厚度 为胸径减去树皮木芯的长度，精确至0.1cm。

木材基本密度测定：杉木木材基本密度的测定材料，去皮的完整的一半木芯，余下 一半木芯用作其他指标测定。测定方法采用最大饱和含水量法，即将木芯浸入盛净水的 容器中，每天定时换水，一个星期后开始每隔1天测定其重量，至其达饱和含水量后， 然后置于105℃烘箱中烘至恒重，取出后用电子天平快速称重，再按如下公式计算木材 基本密度：

$G_f=\frac{1}{\frac{m_m-m_0}{m_0}+\frac{1}{G_{S_0}}}\times 1000$

其中G_f表示木芯的基本密度，m_m表示最大饱和水时木芯的重量，单位为克(g)，m₀表示完全烘干后木芯的重量，单位为克(g)，表示木材的比重，其值介于1.50和1.56 之间，取其平均值1.53，由于水的密度为1g/cm³，故木芯的基本密度单位为kg/m³，结 果精确至0.1kg/m³。

木材红心材率和偏心率的测定：木材的红心材是指红心材指自髓心向外木材颜色呈 红褐色部分的心材。木材红心材率(red-heart ratio，RHR)采用木芯中红心部分的长度 与去树皮后木芯的全长比例表示，即：红心材率＝红心全长÷去皮木芯全长×100%。

偏心率是指椭圆两焦点间的距离和长轴长度的比值。即某一椭圆轨道与理想圆环的 偏离，长椭圆轨道偏心率高，而近于圆形的轨道偏心率低。木材的偏心率(eccentricity  ratio，ER)计算公式为：

木材的偏心率＝|髓芯到树皮之间的木芯长度÷去皮木芯全长-0.5|×100%。

木材微纤丝角度的测定：切取每个木芯最外一个生长轮晚材即靠近树皮一侧的薄 片，厚度0.1～0.3cm左右的，切片时尽量保持切片与弦切面平行。用X射线衍射法进 行测定，实验使用具有超能检测器的X射线衍射仪X’PERTPRO，设定其电压为40kV， 电流为40mA，X射线束截面积为4m×2m，2θ角为22.6°，试样旋转角度为360°，扫描 一个试样的时间为3min。将试样竖直粘贴在试样台上进行测试，扫描得到试样的衍射 强度曲线，在强度分布曲线上确定两拐点，并分别作切线。它们与横轴的两交点之间的 1/2距离定义为角距离(T)。平均微纤丝角为

$\overline{\mathrm{MFA}}=0.6T$

其中MFA表示微纤丝角即Micro Fiber Angle，系数0.6表示0.6法即平均微纤丝角 与峰值高度60%处的宽度相对应，通常假定峰值60%处弧的角宽度的一半等于平均微纤 丝角^[3]。

遗传分析的线性混合模型：y＝Xb+Zf+e

式中y为性状表型观测值构成的向量；b为固定效应构成的向量，且第一个元素为群 体均值；f为随机效应构成的向量；e为剩余误差向量；X为固定效应的关联矩阵；Z为 随机效应的关联矩阵。假设随机误差和随机效应向量独立且服从相同的分布，即

E(f)＝0，E(e)＝0，E(y)＝Xb，Var(y)＝ZGZ′+R，Cov(f，e)＝0

，其中R和G分别为误差方差协方差矩阵和随机效应方差协方差矩阵，A为分子亲缘 关系矩阵，为误差方差，随机效应方差。

Herderson在1953年提出了求解上式的混合模型方程组法(MME，Mixed model  equations)：

$\left(\begin{array}{cc}{X}^{\prime }{R}^{-1}X& X^{\prime }{R}^{-1}Z\\ Z^{\prime }{R}^{-1}X& Z^{\prime }{R}^{-1}Z+G^{-1}\end{array}\right)\left(\begin{array}{c}\hat{b}\\ \hat{f}\end{array}\right)=\left(\begin{array}{c}{X}^{\prime }{R}^{-1}y\\ Z^{\prime }{R}^{-1}y\end{array}\right)$

由上式可以计算得到固定效应的估计值和随机效应的预测值分别为：

$\left(\begin{array}{c}\hat{b}\\ \hat{f}\end{array}\right)={\left(\begin{array}{cc}{X}^{\prime }X& X^{\prime }Z\\ Z^{\prime }X& Z^{\prime }Z+A^{-1}\frac{{\sigma }_e^2}{{\sigma }_f^2}\end{array}\right)}^{-1}\left(\begin{array}{c}{X}^{\prime }y\\ Z^{\prime }y\end{array}\right)$

似然比检验(likelihood ratio tests，LRT)被用于检验模型的有效性和估计参数的显 著性。似然函数取得最大值表示相应的参数能够使得统计模型最为合理，如果参数限制 是正确的，那么加入这样一个参数应当不会造成似然函数最大值的大幅变动，一般使用 两者的比例来进行比较，利用对数似然标准进行模型的对比，其显著性检验公式为：

LRT＝2×(logL.offull model-logL.of reduced model)

其中L.ofull model表示复杂模型的似然值，L.of reduced model表示简单模型的似然值。 LRT检验中，自由度等于在复杂模型中增加的模型参数的数目。由于LR近似的服从卡 方分布，这样根据卡方分布临界值表，就可以判断模型差异是否显著，即如果q为自由度，则表示估计的方差分量是显著的。在ASReml中标准误通过泰勒展开得到。

在本发明中，EXCEL用于数据的整理，R语言用于作图，ASReml2.0用于遗传分 析。

测定的部分杉木经济性状数据如表1所示。

表1部分杉木经济性状数据

实施例2

以实施例1测定的数据为原始数据，利用基于空间分析方法估计空间残差方差与空 间自相关、限制性极大似然法估计遗传方差和最佳线性无偏预测杉木个体育种值的混合 线性模型方法，分析杉木经济性状的遗传变异、遗传控制和性状间的遗传相关，预测杉 木每个个体每个性状的育种值；具体过程如下：

(1)数据整理与标准化，使数据符合正态分布。

(2)空间分析。空间分析模型如下：

Y＝Xb+Zf+ξ+η

其中，Y是数据向量，固定效应向量为均值向量，b和f是随机效应向量，X和Z 是观测值与随机效应相对应的设计矩阵，ξ是空间随机误差自变量向量，η是随机残差 向量。

空间分析结果为杉木木材基本密度、微纤丝角、胸径和树高的空间残差方差空间趋 势图。如图1、图2、图3和图4所示，分别为木材基本密度的空间残差方差空间趋势 图、微纤丝角的空间残差方差空间趋势图、胸径的空间残差方差空间趋势图和树高的空 间残差方差空间趋势图。

(3)遗传参数估计，利用限制性极大似然法(REML)估计方差分量，并根据各方差 组分计算遗传力、遗传相关。

(4)育种值计算，建立亲缘关系矩阵，根据计算得到的家系方差和环境方差计算 育种值，公式为：

$\left(\begin{array}{c}\hat{b}\\ \hat{f}\end{array}\right)={\left(\begin{array}{cc}{X}^{\prime }X& X^{\prime }Z\\ Z^{\prime }X& Z^{\prime }Z+A^{-1}\frac{{\sigma }_e^2}{{\sigma }_f^2}\end{array}\right)}^{-1}\left(\begin{array}{c}{X}^{\prime }y\\ Z^{\prime }y\end{array}\right)$

其中y为性状表型观测值构成的向量；b为固定效应构成的向量，且第一个元素为 群体均值；f为随机效应构成的向量；e为剩余误差向量；X为固定效应的关联矩阵；Z 为随机效应的关联矩阵。假设随机误差和随机效应向量独立且服从相同的分布，即 E(f)＝0，E(e)＝0，E(y)＝Xb，Var(y)＝ZGZ′+R，Cov(f，e)＝0

，其中R和G分别为误差方差协方差矩阵和随机效应方差协方差矩阵，A为分子亲缘 关系矩阵，为误差方差，随机效应方差。计算得到的胸径的育种值如表2所示。

表2胸径的育种值(部分数据)

个体树号 胸径育种值 木材密度育种值 3101 0.11 -0.91 3102 0.03 -0.53 3103 0.10 -0.53 3104 0.01 -0.53 3205 0.37 -0.93 3206 0.33 -0.76 3207 0.31 -0.76 3208 0.32 -0.76 3309 -0.11 -0.24 33010 -0.02 -0.40 33011 -0.18 -0.24

33012 -0.14 -0.24 34013 0.16 -0.44 34014 0.07 -0.44 34015 0.04 -0.44 34016 0.26 -0.55 35017 0.55 0.16 35018 0.57 0.16 35019 0.77 0.31

实施例3

根据杉木的多性状指数选择公式，选出杉木生长与材性在遗传上兼优的优良个体； 多性状指数选择公式如下：

$I_i=0.7{\hat{a}}_{i1}+0.3{\hat{a}}_{i2}$

式中，I_i是基因型i的选择指数值；0.7和0.3分别是杉木生长性状和材性性状的权 重；和是基因型i每个性状的最佳线性无偏预测育种值(即为实施例2计算出 的胸径育种值和木材密度育种值)。

(1)福建省沙县官庄国有林场

共调查2243株，各性状的平均值、遗传方差估计值、加性遗传变异系数和单株狭 义遗传力见表3。

表3福建省沙县官庄国有林场数据表

注：**P<0.01；ns不显著；NE表示没有估计值或为0；单尾似然比检验显著性水平； “±”后的数值为估计的标准误。

胸径与树高性状间的加性遗传相关估计值为0.92±0.05(P<0.01)，预测的育种值和 根据选择指数并以1%的选择强度中选的优良个体见表4。

表4福建省沙县官庄国有林场选中的优良个体

个体编号 个体所在位置 育种值(cm) 树高(m) 胸径(cm) G.32102678 R9B76P670F21T2 1.06 11.2 17.6 G.32102892 R10B82P723F21T4 1.04 11.9 19.0 G.32202002 R7B57P501F22T2 0.85 12.0 19.9 G.32202004 R7B57P501F22T4 0.91 11.7 20.3 G.34301551 R5B44P388F43T3 1.20 9.2 15.3 G.34302043 R7B58P511F43T3 1.07 11.4 18.9 G.3501104 R4B32P276F5T4 1.03 10.5 16.8 G.3502387 R8B68P597F5T3 1.07 10.1 17.1 G.36002621 R9B74P656F60T1 0.98 10.7 17.4 G.3600525 R2B15P132F60T1 1.24 10.2 16.8 G.36601131 R4B32P283F66T3 0.77 10.7 15.6 G.36601808 R6B51P452F66T4 0.81 9.6 16.2 G.36701273 R4B36P319F67T1 0.89 10.4 18.0 G.36701843 R6B52P461F67T3 1.09 10.2 18.6 G.37101417 R5B40P355F71T1 1.29 10.6 16.9 G.37102802 R9B79P701F71T2 1.31 10.2 17.7 G.37502130 R7B60P533F75T2 1.10 11.2 19.6 G.3750924 R3B26P231F75T4 1.50 9.2 18.9 G.37801206 R4B34P302F78T2 0.84 10.5 16.2 G.37801636 R6B46P409F78T4 0.88 10.5 17.7 G.3801211 R4B35P303F8T3 1.05 11.8 20.7 G.3802492 R8B71P623F8T4 0.87 13.4 18.1 G.3901888 R6B54P472F9T4 1.36 12.9 21.7 G.3902528 R8B72P632F9T4 1.01 11.6 18.5

注：R表示重复，B表示区组，P表示小区，F表示家系，T表示株号。

(2)福建省邵武卫闽国有林场

共调查2398株，树高的空间行相关为0.88，空间列相关为0.97，空间残差方差为 0.86m²；胸径的空间行相关为0.86，空间列相关为0.96，空间残差方差4.33cm²。各性 状的平均值、遗传方差估计值、加性遗传变异系数和单株狭义遗传力见表5。

表5福建省邵武卫闽国有林场数据表

注：*P<0.05；**P<0.01；ns不显著；NE表示没有估计值或为0；单尾似然比检验显著 性水平；“±”后的数值为估计的标准误。

胸径与树高性状间的加性遗传相关估计值为0.83±0.01(P<0.01)，树高与冠幅的加 性遗传相关估计值为0.67±0.18(P<0.01)，胸径与冠幅的加性遗传相关估计值为0.78±0.11 (P<0.01)，树高与红心材比率的加性遗传相关估计值为0.19±0.33(不显著)，胸径与红 心材比率0.42±0.25的加性遗传相关估计值为(不显著)。预测的育种值和根据选择指数 并以1%的选择强度中选的优良个体见表6。

表6福建省邵武卫闽国有林场选中的优良个体

个体编号 个体所在位置 育种值(cm) 树高(m) 胸径(cm) W.31905415 R17B153P1354F19T3 1.20 11.5 21.0 W.31905700 R18B161P1425F19T4 1.56 12.6 25.7 W.32203640 R12B103P910F22T4 1.12 9.9 18.4 W.32204774 R15B135P1194F22T2 0.93 10.3 17.9

W.32403921 R13B111P981F24T1 1.13 9.2 18.2 W.32403923 R13B111P981F24T3 1.11 10.4 18.0 W.34303756 R12B106P939F43T4 1.20 11.3 20.2 W.34305744 R18B162P1436F43T4 0.95 10.3 19.7 W.34905069 R16B143P1268F49T1 1.01 9.0 17.1 W.34905428 R17B153P1357F49T4 1.22 11.2 21.4 W.3504304 R14B122P1076F5T4 1.23 10.5 22.0 W.3505908 R19B167P1477F5T4 1.06 10.6 20.9 W.35804188 R14B118P1047F58T4 1.16 11.9 20.4 W.35806284 R20B177P1571F58T4 1.19 10.0 19.7 W.36404117 R13B116P1030F64T1 1.04 10.2 18.9 W.36405863 R19B165P1466F64T3 1.02 9.8 19.8 W.36504155 R13B117P1039F65T3 1.25 11.4 20.3 W.36506359 R20B179P1590F65T3 1.26 9.5 22.1 W.36604332 R14B122P1083F66T4 1.37 10.5 21.8 W.36606396 R20B180P1599F66T4 1.39 8.7 22.2 W.37503628 R12B102P907F75T4 1.47 9.7 23.2 W.37505329 R17B150P1333F75T1 0.91 10.6 18.2 W.3804409 R14B125P1103F8T1 1.34 9.4 20.8 W.3804412 R14B125P1103F8T4 1.25 9.1 20.0

注：R表示重复，B表示区组，P表示小区，F表示家系，T表示株号。

(3)福建省国有来舟林业试验场

共调查2536株，各性状的平均值、遗传方差估计值、加性遗传变异系数和单株狭 义遗传力见表7。

表7福建省国有来舟林业试验场数据表

性状 平均值 家系间方差 表型方差 加性遗传变异系数 单株狭义遗传力 树高 8.0m 0.0668±0.0207** 1.0868 6.4% 0.23±0.07** 胸径 10.5cm 0.20±0.06** 3.96 8.7% 0.20±0.06**

性状 平均值 家系间方差 表型方差 加性遗传变异系数 单株狭义遗传力 树皮厚度 1.5cm 0.0119±0.0054** 0.1243 14.7% 0.35±0.37** 基本密度 351.5kg/m³NE 32.0 NE NE 微纤丝角 14.9° 0.67±0.30** 5.52 10.6% 0.48±0.64** 红心比率 41.2% 0.0002±0.0001^ns0.0015 6.0% 0.40±0.56^ns偏心率 2.4% 0.12±0.11^ns16.9 21.0% 0.09±0.11^ns

注：*P<0.05；**P<0.01；ns不显著；NE表示没有估计值或为0；单尾似然比检验显著 性水平；“±”后的数值为估计的标准误。

胸径与树高性状间的加性遗传相关估计值为0.82±0.07(P<0.01)，树高与树皮厚度 的加性遗传相关估计值为0.60±0.22(P<0.01)，胸径与树皮厚度的加性遗传相关估计值 为0.40±0.23(不显著)，胸径与微纤丝角的加性遗传相关估计值为0.10±0.24(不显著)， 树高与微纤丝角的加性遗传相关估计值为-0.09±0.25，树皮厚度与微纤丝角的加性遗传 相关估计值为0.25±0.33(不显著)。预测的育种值和根据选择指数并以1%的选择强度 中选的优良个体见表8。

表8福建省国有来舟林业试验场选中的优良个体

个体编号 个体所在位置 育种值(cm) 树高(m) 胸径(cm) L.31606941 R22B196P1736F16T1 1.33 10.3 15.7 L.31607616 R24B215P1904F16T4 1.18 10.1 15.8 L.3208999 R29B254P2250F2T3 1.32 9.7 16.0 L.3209317 R30B263P2330F2T1 1.18 8.4 13.6 L.32208402 R27B237P2101F22T2 1.42 10.0 16.4 L.32209217 R29B260P2305F22T1 1.09 11.8 16.3 L.32407656 R24B216P1914F24T4 1.73 12.5 17.8 L.32407942 R25B224P1986F24T2 1.85 12.9 19.5 L.33907625 R24B215P1907F39T1 1.15 11.0 15.5 L.33908158 R26B230P2040F39T2 0.99 11.9 17.4 L.34306572 R21B185P1643F43T4 1.25 9.5 14.0 L.34308942 R28B252P2236F43T2 1.29 11.3 16.5

L.3507826 R25B221P1957F5T2 1.40 12.0 17.1 L.3508462 R27B239P2116F5T2 1.48 10.3 15.9 L.35808950 R28B252P2238F58T2 1.23 11.8 15.8 L.35809482 R30B267P2371F58T2 1.81 9.7 17.7 L.36508170 R26B230P2043F65T2 1.31 11.6 19.0 L.36508850 R28B249P2213F65T2 1.17 11.4 15.5 L.36606664 R21B188P1666F66T4 1.39 10.3 15.9 L.36607963 R25B224P1991F66T3 1.44 10.8 17.6 L.36806669 R21B188P1668F68T1 1.07 8.5 15.0 L.36806670 R21B188P1668F68T2 1.09 9.4 15.1 L.3706939 R22B196P1735F7T3 1.29 12.3 15.1 L.3709176 R29B259P2294F7T4 1.20 11.6 16.1

注：R表示重复，B表示区组，P表示小区，F表示家系，T表示株号。

对中选的优良个体经现场踏查，其生长量和材性均表现优良，明显优于未选中的个 体。选出的优良个体可以作为杉木无性快繁的材料，也可以作为杉木下一代育种群体的 育种材料。