-
泡桐(Paulownia)是我国栽培历史最悠久的优质速生用材树种,在我国25个省(市、自治区)均有分布,具有优质、速生、分布和用途广泛以及适用于农林复合经营等众多优良特性。20世纪70年代以来,泡桐产业得到迅猛发展,学者们分别从泡桐的种类、形态特征、生态特性、分布状态以及性状稳定性等开展了多方面的综合研究[1-2],将泡桐分为11个种、2个变种、6个变型,并针对泡桐的传统栽培区(黄淮海平原)先后选育出优良无性系30余个。
近年来随着林权制度的改革,泡桐因其较高的经济和生态收益日益受到重视,加之南方低山丘陵区丰富的土地资源、优越的水热和气候等条件,泡桐种植业得到迅速发展。泡桐在优良无性系选育以及高效栽培技术等方面虽取得一定的成效,但前期研究主要集中在北方,针对适宜南方低山丘陵区域栽植的泡桐优良无性系缺乏系统研究,导致南方栽植的泡桐用材林存在诸多问题,比如盲目引种、品系混杂、主干低矮不通直等,严重阻碍了泡桐产业的健康发展;而且前期的泡桐无性系选择多采用单性状直接选择法(速生、干形、抗丛枝病等)或者单性状排列选择法(逐步实现多个性状的遗传改良),对其进行多性状综合改良研究较少[2-7],大大降低了聚合多性状育种的选择效果[8-9];白花泡桐(Paulownia fortunei (Seem.) Hemsl.)是我国南方泡桐种类的代表树种,分布范围较广,遍布长江流域以南各地,除了一部分人工栽培外,多为天然分布;其树形多为长卵形或塔形、生长速度较快、树干通直高大、自然接干能力较强;木材的亮度和白度较高、材质均匀、色泽淡雅、材性较优[4, 10-11]。本研究针对18个白花泡桐无性系的6个主要性状指标进行观测和统计(选择目前栽培范围较广、市场前景广阔的主栽品种9501为对照),并对各性状的表型变异、遗传变异、表型相关和遗传相关等进行综合比较分析,采用多性状选择指数法对各无性系进行综合评价,以期选择出适宜南方低山丘陵区域栽植的泡桐优良无性系,为当地泡桐种植业的健康发展提供科学指导。
HTML
-
2012年2月对试验地进行整地挖穴(规格为60 cm×60 cm×60 cm),18个参试无性系均按照完全随机区组设计,6株小区,4次重复,株行距为4 m × 5 m,1年生苗木造林(平均苗干高度为2.92 m,平均胸径为6.55 cm),底肥为3 kg有机肥。其中9501为对照(CK),四周均设有保护行,试验林采用一致的抚育管理措施。试验材料的来源详见表 1。
无性系Clone 种源Provenance 种类Specie 来源Origin 01-1 江西抚州 白花泡桐P. fortunei 实生选种seeding selection 01-4 湖南郴州 白花泡桐P. fortunei 实生选种seeding selection 01-11 重庆涪陵 白花泡桐P. fortunei 实生选种seeding selection 01-18 湖北宜昌 白花泡桐P. fortunei 实生选种seeding selection 01-20 浙江衢州 白花泡桐P. fortunei 实生选种seeding selection 01-22 湖南郴州 白花泡桐P. fortunei 实生选种seeding selection 01-23 湖南郴州 白花泡桐P. fortunei 实生选种seeding selection 1-3 湖南张家界 白花泡桐P. fortunei 优树plus tree 1-38 贵州凯里 白花泡桐P. fortunei 优树plus tree 1-39 广西柳州 白花泡桐P. fortunei 优树plus tree 1-43 广西河池 白花泡桐P. fortunei 优树plus tree 1-45 广西河池 白花泡桐P. fortunei 优树plus tree 1-48 湖南衡阳 白花泡桐P. fortunei 优树plus tree 1-58 广西河池 白花泡桐P. fortunei 优树plus tree -1 江西宜春 白花泡桐P. fortunei 实生选种seeding selection -4 湖南邵阳 白花泡桐P. fortunei 实生选种seeding selection -9 广西桂林 白花泡桐P. fortunei 实生选种seeding selection 9501 湖南长沙 白花泡桐P. fortunei 实生选种seeding selection Table 1. The origin of Paulownia clones used in this study
-
2012年3月试验林营建后,调查各无性系初始的苗干高度和胸径。于每年落叶后(11月份左右),分别对18个参试无性系的生长情况进行观察和测定,其中包括胸径(D1.3)、2.6 m直径(D2.6)、苗干高度(H1)、接干高度(H2)、接干0.5 m处直径(De0.5)和接干1.5 m处直径(De1.5),根据公式(1)~(8)分别计算其主干高(H)、苗干削度(β1)、接干削度(β2)、主干削度(β)、苗干材积(V1)、接干材积(V2)、总材积(V)和形数(f)[12]。以已达到早期选择年龄[13]的5年生泡桐的调查数据为基础进行统计分析。为了减少初始苗干高度和胸径对试验结果的影响,分别采用协方差分析方法对高、径生长指标进行调整[14]。
(1) H=H1+H2
(2) β1=(D1.3-D2.6)/(2.6-1.3)
(3) β2=(De0.5-De1.5)/(1.5-0.5)
(4) β=(β1+β2)/2
(5) V1={2.6×D1.32+[D2.6-β1×(H1-2.6)/2]2×(H1-2.6)}×π/40000
(6) V2={[De0.5-β2×(H2/2-0.5)]2}×π×H2/40000
(7) V=V1+V2
(8) f=V1/[(π/40000)×D1.32×H1]
-
方差分析和重复力计算具体参照续九如[15]的林木数量遗传学,用方差分析线性模型xijk=μ+αi+βj+αβijk+eijk对参试的18个泡桐无性系进行方差分析,均以单株测定数据为统计单元(式中μ:总体均值;αi:无性系效应;βj:区组效应;αβijk:无性系和区组的交互效应;eijk:随机误差)。无性系重复力计算公式为H2=1-1/F(式中F:无性系方差分析的F值)。
表型变异系数$ CV = \frac{S}{{\bar x}}$×100%(式中S:某性状的表型标准差;$ {\bar x}$:某性状的群体平均值);遗传变异系数$ GCV = \frac{{\sqrt {{\sigma _a}^2} }}{{\bar x}}$×100% (式中σa2:无性系方差分量;x:某性状的群体平均值)[16];表型相关系数$ {r_{p12}} = \frac{{{{{\mathop{Cov}\nolimits} }_{p12}}}}{{\sqrt {\sigma _{p1}^2\sigma _{p2}^2} }}$(式中Covp12:2个性状的表型协方差;σp12、σp22:2个性状的表型方差);遗传相关系数$ {r_{g12}} = \frac{{{{{\mathop{Cov}\nolimits} }_{g12}}}}{{\sqrt {\sigma _{g1}^2\sigma _{g2}^2} }}$(式中Covg12:2个性状的遗传协方差;σg12、σg22:2个性状的遗传方差)。
利用Smith-Hazel指数选择法进行多性状的综合选择[17],采用等权重法估算无性系各性状的经济权重w,具体计算公式为$ {w_i} = \frac{1}{{{\sigma _i}}}$[6](式中σi:各性状的表型标准差),指数选择的计算公式为I= b′X(式中b′:各性状的指数系数;X:对应性状的平均值),指数系数计算公式为b′=P-1Gw(式中P:选择性状的表型方差协方差矩阵;G:选择性状的遗传方差协方差矩阵)[6, 18-19]。
参照赵承开等[20]方法进行遗传增益(ΔG)和实际增益(G)的计算,但又做了适当的调整,具体计算公式如下:ΔG=(H2M/x)×100%、G=(M/x)×100%,式中M为入选无性系群体某一性状的平均值与CK平均值之差。
2.1. 试验材料
2.2. 试验方法
2.3. 统计分析方法
-
对参试无性系的各性状分别进行方差分析和遗传参数估算,结果表明:胸径、主干高、总材积、接干高、主干削度和形数等6个性状指标在18个无性系间的差异均达到极显著水平(表 2)。其胸径、主干高、总材积、接干高、主干削度和形数的平均值分别为18.282 1、7.452 1、0.139 2、4.611 9、1.258 3和0.964 6,其最高值分别是最低值的2.41、2.88、5.50、9.63、103.47和1.76倍,在立地环境一致的条件下不同无性系间表型的差异客观反映了它们基因型存在差异,表明开展无性系选择是必要的。
性状Traits 胸径D1.3 /cm 主干高H/m 总材积V/m3 接干高H2/m 主干削度β 形数f 平均值Means 18.282 1 7.452 1 0.139 2 4.611 9 1.258 3 0.964 6 无性系间F值F values of clones 12.399 2** 3.444 5** 11.672 7** 8.664 4** 2.824 4** 7.196 9** 无性系间方差分量Variance 37.042 8 3.302 4 0.012 7 8.091 3 0.888 9 0.019 5 互作方差分量
Interaction variance component1.313 8 0.420 3 0.001 3 0.603 0 0.429 9 0.004 2 误差Error 2.987 6 0.958 8 0.001 1 0.933 9 0.314 7 0.002 7 无性系重复力Repeatability of clones 0.964 6 0.872 7 0.895 0 0.925 5 0.516 3 0.785 8 表型变异系数PCV/% 12.840 3 13.869 8 32.780 1 25.113 6 48.093 1 6.495 9 遗传变异系数GCV/% 33.290 9 24.385 8 80.938 9 61.678 4 74.926 2 14.486 8 Table 2. Analysis of variance and genetic parameter estimation
18个参试无性系各性状的Duncan’s多重比较结果如表 3所示,结果表明胸径、主干高、总材积均显著高于CK的有01-22、01-23和1-58,接干高显著高于CK的有01-22、01-23、1-38、1-45和1-58,主干削度显著低于CK的有01-22、01-23和1-58,形数显著高于CK的是01-23和-9。
无性系Clones 胸径D1.3/cm 主干高H/m 总材积V/m3 接干高H2/m 主干削度β 形数f 9501(CK) 18.386 7 cdef 7.403 5 abc 0.143 0 cde 3.642 1 fg 1.292 3 abcd 0.891 9 i 01-1 17.333 5 fgh 6.534 2 d 0.124 4 defg 3.653 6 fg 1.694 6 a 0.948 9 defgh 01-4 16.887 5 gh 7.317 4 bcd 0.124 4 defg 4.220 9 efg 1.488 0 abc 0.940 9 fgh 01-11 18.969 2 cd 7.578 9 ab 0.157 8 bc 4.782 0 bcde 1.044 3 cd 0.998 7 bc 01-18 17.907 5 cdefg 7.576 0 ab 0.139 1 cdef 4.538 2 cde 1.261 8 abcd 0.956 7 cdefg 01-20 17.501 3 efg 6.605 8 cd 0.117 7 efg 3.541 1 g 1.489 3 abc 0.909 0 hi 01-22 20.420 4 b 8.206 7 a 0.199 6 a 5.406 7 b 0.970 6 d 0.964 5 bcdefg 01-23 21.916 3 a 7.855 4 ab 0.187 2 a 6.505 4 a 0.971 6 d 1.044 2 a 1-3 18.100 0 cdefg 7.163 6 bcd 0.128 9 defg 4.296 9 defg 1.450 8 abcd 0.989 8 bcd 1-38 16.852 5 gh 7.783 6 ab 0.124 4 defg 5.056 2 bcd 1.655 4 ab 0.950 1 defgh 1-39 17.664 4 defg 7.675 0 ab 0.129 7 defg 4.804 3 bcde 1.037 2 cd 0.985 0 bcde 1-43 17.484 2 efg 7.210 7 bcd 0.118 0 ef 4.452 7 cde 1.258 0 abcd 0.984 2 bcdef 1-45 18.896 7 cde 7.681 1 abc 0.139 0 cdef 5.138 9 bc 1.116 5 cd 0.944 3 efgh 1-48 16.863 5 gh 7.715 9 ab 0.114 7 fg 4.738 9 bcde 1.160 5 cd 0.951 0 defg 1-58 20.403 8 b 7.837 2 ab 0.176 0 ab 5.006 4 bcde 0.985 1 d 0.994 4 bc -1 16.078 3 h 7.406 5 abc 0.109 1 g 4.331 7 cdef 1.406 0 abcd 0.934 4 gh -4 18.227 5 cdefg 7.283 3 bcd 0.118 7 efg 4.366 8 cdef 1.139 9 cd 0.949 4 defgh -9 19.270 8 bc 7.298 2 bcd 0.146 7 cd 4.498 9 cde 1.189 6 bcd 1.004 0 b Table 3. Duncan's test of the Paulownia clones performance
-
表现型是遗传效应和环境效应共同决定的。重复力是指基因型方差和一般环境方差之和在表型方差中所占的比例。本研究18个泡桐无性系各性状的重复力如表 2所示,其中胸径的重复力最大(0.964 6),主干高、总材积和接干高的重复力也均达到了0.872 7以上,表明泡桐无性系各性状受遗传因素的影响较大,受外界环境的制约较小,对其进行表型选择获得优良无性系的概率较大[21]。形数的重复力次之(0.785 8),重复力最小的是主干削度(0.516 3)。
较大的变异可以提供更大的选择空间[22],在18个参试无性系的6个选择性状中,总材积和主干削度的遗传变异系数最大,接干高和胸径的遗传变异系数次之,主干高和形数的遗传变异系数较小,遗传变异系数的变化范围为14.486 8%~ 80.938 9%;表型变异系数最大的为主干削度和总材积,其次是接干高、主干高和胸径,最小的是形数,表型变异系数的变化范围为6.495 9%~ 48.093 1%;该结果表明参试的18个无性系之间存在着比较丰富的遗传变异和表型变异,变幅最大的性状是主干削度,其次是接干高。
-
泡桐无性系各性状间相关系数如表 4所示,胸径、主干高、总材积和接干高4个生长性状之间均呈现出极显著的表型和遗传正相关关系;除与主干高、接干高分别呈不显著和显著的表型负相关关系外,主干削度与其它性状间的表型和遗传负相关关系均达到极显著水平;形数与主干高表型正相关关系不显著外,与其它性状间的表型和遗传相关关系均呈极显著水平;遗传相关对于育种来说更为重要,它是指表型性状中可以稳定遗传部分的相关,表 4结果显示:18个参试无性系的4个生长性状指标中遗传相关系数最高的可达到0.949 3(胸径与总材积间的相关系数),其两两之间的遗传相关系数均大于0.583 3,变化范围为0.583 3~ 0.949 3。主干削度与4个生长性状间的遗传相关系数的绝对值均接近于1,形数与其它5个性状间的遗传相关系数的变化范围为-0.561 6~ 0.784 2。
性状Traits 胸径D1.3/cm 主干高H/m 总材积V/m3 接干高H2/m 主干削度β 形数f 胸径D1.3 0.259 8** 0.816 5** 0.411 3** -0.165 5** 0.248 3** 主干高H 0.583 3** 0.521 2** 0.869 8** -0.098 9 0.023 2 总材积V 0.949 3** 0.717 6** 0.527 7** -0.243 7** 0.152 3** 接干高H2 0.760 0** 0.825 1** 0.756 6** -0.146 0* 0.286 9** 主干削度β -0.972 5** -1.150 6** -0.827 0** -1.044 1** -0.504 3** 形数f 0.649 4** 0.371 3** 0.486 7** 0.784 2** -0.561 6** 注:左下角为遗传相关系数,右上角为表型相关系数。*表示相关关系达到显著水平(P<0.05),**表示相关关系达到极显著水平(P<0.01)。
Notes: Phenotypic correlation coefficients in the upper triangle, and genetic correlation coefficients in the below triangle. * and ** represent correlation is significant at 5% level and 1% level, respectively.Table 4. Correlation coefficients between traits of Paulownia clones
-
为了从多性状角度来评价泡桐无性系,本研究利用18个参试无性系的胸径、主干高、总材积、接干高、主干削度和形数等6个性状指标对其进行综合评价。首先按照等权重法确定6个性状指标的经济权重(表 5),后分别用等权重法、强调生长改良和强调干形改良三种方法建立多性状综合指数方程进行整理分析。根据各性状间的相关关系,选择不同的性状组合,分别构建等权重的指数方程(I1-I3)、强调生长的指数方程(I4-I6)和强调干形的指数方程(I7-I9),并分别对9个指数方程的综合选择进展(表 6)和各性状的遗传进展进行统计分析(表 7)。表 6和表 7的结果表明:在等权重的各指数方程中,以胸径、主干高、总材积、接干高、主干削度和形数这6个性状组合的指数方程I1的综合选择进展最大(13.272 4),此时胸径和接干高的遗传进展也最大;强调生长改良的3个指数方程中,虽然I5的主干高和总材积的遗传进展最大,但I4的综合选择进展、胸径和接干高的遗传进展均最大;强调干形改良的3个指数方程中,I7的综合选择进展、胸径和接干高的遗传进展均最大。
方法Method 指数方程
Selection indices经济权重Economic weight 胸径D1.3 主干高H 总材积V 接干高H2 主干削度β 形数f 等权重
Equal weightI1 0.426 0 0.967 5 21.915 2 0.863 4 -1.652 5 15.960 0 I2 0.426 0 0.967 5 21.915 2 0.863 4 -1.652 5 I3 0.426 0 0.967 5 21.915 2 0.863 4 15.960 0 强调生长
Strengthening growthI4 1.278 0 2.902 5 65.745 5 2.590 2 -1.652 5 15.960 0 I5 1.278 0 2.902 5 65.745 5 2.590 2 -1.652 5 I6 1.278 0 2.902 5 65.745 5 2.590 2 15.960 0 强调干形
Strengthening stem formI7 0.426 0 0.967 5 21.915 2 0.863 4 -4.957 4 47.880 1 I8 0.426 0 0.967 5 21.915 2 0.863 4 -4.957 4 I9 0.426 0 0.967 5 21.915 2 0.863 4 47.880 1 Table 5. Estimation method and results of 6 traits economic weight
序号No. 性状配合及指数方程式Selection indices equation with various traits 综合选择进展ΔI I1 I=+1.646 9×1+3.551 3×2-41.965 0×3+0.405 5×4-5.985 6×5+13.389 8×6 13.272 4 I2 I=+1.538 1×1+3.813 7×2-39.618 4×3-0.046 2×4-6.337 3×5 11.662 9 I3 I=+1.096 1×1+2.797 9×2-8.128 6×3+0.221 4×4+32.872 9×5 10.980 8 I4 I=+3.959 2×1+9.079 1×2-86.517 6×3+0.622 5×4-16.360 0×5+22.267 7×6 33.341 7 I5 I=+3.862 8×1+8.774 3×2-85.086 9×3+0.590 2×4-17.312 0×5 31.745 1 I6 I=+3.057 4×1+9.678 2×2-18.822 7×3-0.798 4×4+83.502 0×5 29.386 8 I7 I=+2.628 4×1+5.126 2×2-81.342 5×3+0.999 4×4-7.582 5×5+31.291 6×6 19.952 3 I8 I=+2.289 4×1+6.480 4×2-73.386 6×3-0.775 1×4-8.037 3×5 15.003 4 I9 I=+1.326 8×1+1.513 3×2-13.691 8×3+1.684 0×4+47.989 7×5 14.825 7 Table 6. Combined selection progress and selection indices equation with growth and stem form
序号No. 胸径D1.3 主干高H 总材积V 接干高H2 主干削度β 形数f I1 6.745 6 1.887 6 0.114 0 3.240 2 -0.967 5 0.105 1 I2 6.683 2 1.985 7 0.115 2 3.165 1 -0.990 8 I3 6.316 6 1.540 0 0.108 7 3.000 8 0.114 4 I4 6.770 5 1.923 7 0.115 4 3.214 3 -0.965 8 0.100 0 I5 6.741 1 1.961 6 0.115 7 3.175 3 -0.971 1 I6 6.307 6 1.599 5 0.110 7 2.966 2 0.108 0 I7 6.634 9 1.808 1 0.110 5 3.250 2 -0.960 5 0.112 5 I8 6.517 6 2.023 7 0.113 4 3.123 1 -1.026 2 I9 6.211 4 1.392 1 0.102 8 3.010 9 0.124 9 Table 7. Genetic progress on the traits of all selection indices
综上所述,等权重的指数方程I1、强调生长改良的指数方程I4和强调干形改良的指数方程I7比较理想。按照20%的选择强度依据指数方程I1从参试无性系中选择出排在前4位的无性系依次是01-23、01-22、1-58和1-45,指数方程I4排在前4位的无性系依次是01-23、01-22、1-58和1-45,指数方程I7排在前4位的无性系依次是01-23、1-58、01-22和01-11。3个理想指数方程分别选出的4个无性系都包含01-23、01-22和1-58,可作为入选的生长、干形兼优的优良无性系,入选率为17%。表 8是入选的3个无性系和CK各性状的平均值,其中,01-23、01-22和1-58这3个入选无性系的胸径、主干高、总材积、接干高和形数与CK相比均不同程度的得到了提高,幅度分别为10.97%~ 19.20%、5.86%~ 10.85%、23.08%~ 39.58%、37.46%~ 78.62%和8.14%~ 17.08%;其主干削度与CK相比降低的幅度为23.77%~ 24.89%;3个入选无性系的6项指标均得到不同程度的改良。
无性系Clone 胸径D1.3 主干高H 总材积V 接干高H2 主干削度β 形数f 01-23 21.916 3 7.855 4 0.187 2 6.505 4 0.971 6 1.044 2 01-22 20.420 4 8.206 7 0.199 6 5.406 7 0.970 6 0.964 5 1-58 20.403 8 7.837 2 0.176 0 5.006 4 0.985 1 0.994 4 9501 (CK) 18.386 7 7.403 5 0.143 0 3.642 1 1.292 3 0.891 9 Table 8. Means of 3 traits of Paulownia clones and its value of index formula
3个入选群体各性状的遗传增益和实际增益如图 1所示。在所有性状中,接干高的实际增益最高,达到43.31%,其次是总材积(32.04%),二者的遗传增益也较高,分别达到40.08%和28.68%,表明入选群体接干高和总材积在获得较大实际增益的基础上,也显著提高了其遗传增益;其次是主干削度和胸径,其实际增益分别达25.16%和13.82%,但由于主干削度的重复力较小(0.516 3),导致其遗传增益与实际增益相差稍大[23];入选群体形数和主干高的遗传增益和实际增益稍小。