桉树无性系生长性状基因型与环境互作效应及优良无性系选择.pdf

桉树无性系生长性状基因型与环境互作效应及优良无性系选择王楚彪1，熊涛2，张磊2，何文亮1，卢万鸿1，林彦1，罗建中1（1.中国林业科学研究院速生树木研究所，广东湛江524022；2.广西国有东门林场，广西扶绥532199）摘要：【目的目的】对 20 个桉树Eucalyptus 无性系，包括 1 个赤桉 E.camaldulensis 无性系和 19 个杂交桉无性系，分别在广西南宁、广西钦州、广东雷州、广西扶绥和海南定安的立地条件下进行测试，探讨基因型、环境及其互作对无性系选育的影响，旨在为桉树的无性系筛选适生环境并推广提供借鉴。【方法方法】以各地点所测试的 4.5 年生桉树无性系为研究对象，基于各无性系的树高、胸径、单株材积、保存率的调查数据，比较各地点各无性系的生长表现。对单地点、多地点建立线性混合模型，计算各性状的方差分量和无性系重复力，同时计算单地点的最佳线性无偏预测(BLUP)值和多地点各性状的方差效应量。最后通过材积平均值和 BLUP 值进行基因型和基因型与环境互作双标图(GGE)分析。【结果结果】单地点分析中，各地点的最优无性系生长排序不同，无性系各性状受到环境的影响也不同。多地点分析中，桉树无性系基因型、环境以及基因型与环境互作存在显著差异，各性状的无性系重复力为 0.78160.8685，其中胸径的无性系重复力最高且表现最稳定，树高易受环境影响。基于单株材积性状，筛选出 EC186、EC188、EC184 为优良无性系。基于材积平均值和 BLUP 值进行 GGE 模型分析得出 EC183 和 EC184 为高材积、高稳定性的无性系。【结论结论】20 个无性系中，EC183 和 EC184 为高材积、高稳定性的无性系。5 个地点相比，广西南宁地区适合这批无性系的快速生长，广西钦州为这批无性系理想的选材环境。图 1 表 5 参 38关键词：桉树无性系；BLUP-GGE 联合分析；重复力；基因型与环境互作中图分类号：S722文献标志码：A文章编号：2095-0756(2023)05-0951-10GenotypeenvironmentinteractionongrowthtraitsofEucalyptusclonesandselectionofeliteclonesWANGChubiao1，XIONGTao2，ZHANGLei2，HEWenliang1，LUWanhong1，LINYan1，LUOJianzhong1（1.ResearchInstituteofFast-growingTrees,ChineseAcademyofForestry,Zhanjiang524022,Guangdong,China;2.GuangxiDongmenStateForestFarm,Fusui532199,Guangxi,China）Abstract:Objective20Eucalyptusclones,including1Eucalyptus camaldulensisand19Eucalyptushybridclones,are tested at Nanning,Qinzhou,Fusui in Guangxi,Leizhou in Guangdong and Dingan in Hainan,respectively.And,theinfluenceofgenotypeandenvironmentinteractiononEucalyptusclonesweremainlydiscussed.TheresultscanprovideabasementforscreeningsuitableenvironmenttothetestedclonesandreferenceforthepromotionofexcellentclonesofEucalyptus.MethodThe4.5-year-oldEucalyptusclonestestedatdifferentsiteswereselectedastheresearchobjects.Basedonthesurveydataoftreeheight,diameteratbreastheight(DBH),individualvolumeandsurvivalrateofeachclone,thegrowthperformanceofeachclone收稿日期：2022-11-28；修回日期：2023-05-17基金项目：中国林业科学研究院基本科研业务费专项资助(CAFYBB2021MB002)；广东省林业科技创新项目(2017KJCX031,2019KJCX014)作者简介：王 楚 彪(ORCID:0000-0001-5417-4629)，副 研 究 员，博 士，从 事 林 木 遗 传 育 种 研 究。E-mail:。通信作者：罗建中(ORCID:0000-0001-8656-9254)，研究员，博士，从事林木遗传育种研究。E-mail:浙江农林大学学报，2023，40（5）：951960Journal of Zhejiang A&F Universitydoi:10.11833/j.issn.2095-0756.20220734atdifferentsiteswascompared.Alinearmixingmodelwasestablishedforsinglesiteandmulti-sitetocalculatethevariancecomponentsandclonalrepeatabilityofeachtrait.Thebestlinearunbiasedpredictive(BLUP)valueofsinglesiteandthevarianceeffectsizeofeachtraitatmulti-sitewerealsocalculated.Finally,theaveragevolumeandBLUPvalueswereusedforgenotypemaineffects(G)andgenotypeenvironmentinteraction(GE)analysismodel(GGE).ResultInsinglesiteanalysis,theorderofoptimalclonesineachsiteweredifferent,andthetraitsindifferentsiteswerealsoaffectedbydifferentenvironments.Inthemulti-siteanalysis,thereweresignificantdifferencesingenotype,environmentandGEinteractionamongEucalyptusclones.Theclonerepeatabilityofeachtraitrangedfrom0.7816to0.8685,amongwhichtheDBHtraithadthehighestandmoststablerepeatability,thetreeheighttraitwasmosteasilyaffectedbyenvironment.Basedontheindividualvolumetrait,EC186,EC188andEC184wereselectedassuperiorclones.BasedontheaveragevolumeandBLUPvalue,theGGEmodelanalysisshowedthatEC183andEC184werecloneswithhighvolumeandhighstability.ConclusionAmongthe20clones,EC183andEC184areclonewithhighindividualvolumeandhighstability.Comparedwiththefiveregions,NanninginGuangxiissuitablefortherapidgrowthoftheseclones,QinzhouinGuangxiistheidealenvironmentforthematerialselectionoftheseclones.Ch,1fig.5tab.38ref.Key words:Eucalyptusclone;BLUP-GGEanalysis;repeatability;genotypeenvironmentinteraction桉树Eucalyptus 生长迅速，作为一种重要工业用材树种，在满足日益增长的全球木材供应需求方面发挥了重要作用1。在以往桉树育种策略中，培育速生桉树无性系能短时间内满足木材需求，前提需在无性系区域试验基础上，选出适宜本区域的优良无性系以推广应用2。材积是桉树重要的经济性状，对木材产量有显著影响。然而，推广的无性系在不同生境生长难免出现经济性状不良的问题，因此桉树育种策略需兼顾培育适应范围广的优良无性系。区域试验为选育适应性范围广的品种提供依据，且能反映植物生长状况及性状、品种稳定性的问题。区域试验也存在着基因型与环境互作(genotypeenvironmentinteraction,GE)影响品种在不同环境下的排序。GE 指基因型的相对表现在不同环境下缺乏稳定性3，即基因型与环境互作效应越强，稳定性越低，是影响品种稳定性的主要原因。基因型与环境互作研究的重要性包括 3 点：其一，可认识不同性状在不同生长环境中的遗传规律，如赤桉E.camaldulensis 的保存率受土壤盐浓度的影响4，桉树无性系的胸径受株行距的正向影响5，翅茎西番莲Passiflora alata 的生长性状受海拔的影响6等；其二，可为特定生境选择经济性状表现最优的基因型，如在特定的恶劣环境中，选育出抗寒7、抗盐抗旱8、抗风9的桉树品种，使跨区域(或气候)种植成为可能；其三，可为多种生境选择稳定性好的优良基因型，如在不同土壤特性、气候类型的环境中选择较优的桉树品种1012。通过基因型与环境互作的研究，可选择适应范围广、品种产量稳定的优良无性系。几乎所有重要的林木品种的基因型和环境互作都有报道，包括杨属Populus13、松属 Pinus14、柳属Salix15等。对桉树无性系与环境互作的研究，目前大多数使用的分析方法有稳定性分析、主效可加互作可乘(additivemainmultiplicativeinteraction,AMMI)分析、主成分分析以及近期流行的基因型和基因型与环境互作双标图 genotypemaineffects(G)andgenotypeenvironmentinteraction(GE),GGE 分析方法等10。其中 GGE 的优势是在 AMMI 分析基因型与环境互作效应的基础上，又更多考虑基因型的作用16。此外，对遗传和表型参数的估计，包括各性状方差分量、无性系重复力、遗传增益等，对育种策略的成功至关重要17。如仅从一个环境中估算遗传和表型参数，未考虑较大的 GE 影响，会导致结果不准确，因此，从多个环境中获得的参数具有参考价值18。目前，仍缺乏对桉树基因型与环境互作效应的研究，已有的研究也仅在一个地点开展1921，因此评估栽培桉树不同基因与环境互作的效应具有重要意义。本研究通过研究不同地点中桉树无性系的生长性状，筛选各地点中的适生品种；对单地点、多地点建立线性混合模型以计算遗传参数；最后以材积平均值和最佳线性无偏预测(bestlinearunbiased952浙江农林大学学报2023 年 10 月 20 日prediction,BLUP)值进行 GGE 分析，选出优良无性系及其理想种植地点，为桉树优良品种的推广提供借鉴。1材料和方法1.1材料参试无性系共 21 个，4.5 年生，详见表 1，以种植面积较大的无性系 DH32-29 作为对照。1.2试验设计及环境概况采用完全随机区组设计，每个小区 4 行，每行5 株，株行距为 2m3m，样地设置 4 个重复，1 个重复内的各个无性系重复 20 次。于 2016 年 6 月种植，2020 年 12 月调查树高、胸径、保存率。试验地分布在广西、广东及海南，详见表 2。1.3数据统计与分析1.3.1生长表现树高(H,m)和胸径(D,cm)分别使用瑞典 Vertex超声波林木测高仪和卷尺测量，计算单株材积(V,m3株1)22，公式为：V=D2H/30000。各无性系保存率为活立木株数占总株数的百分比。表型材积变异系数(Cv)公式为：Cv=/100%，其中 表示无性系材积标准差，表示无性系材积平均值。1.3.2方差分析和遗传参数单地点方差分析使用线性混合模型(linearmixedmodel)，模型23为：Yijk=+Bi+Cj+BCij+ijk。其中 Yijk为单点试验第 i 重复第 j 无性系第 k 单株的观测值，为群体的总体平均值，Bi为第 i 重复效应，Cj为第 j 无性系效应，BCij为第 j 无性系和第 i 重复的互作效应，ijk为剩余项。其中以无性系、无性系与重复的互作效应为随机效应，获得单地点各无性系的 BLUP 值用于GGE 分析。Yhij=+Lh+B(L)i(h)+Cj+ILChj+hij多地点方差分析同样建立线性混合模型，计算各性状的方差分量和方差效应量，模型23为：。其中，Yhij为单点试验第2(A)=SSASSA+SSEh 地点内第 i 重复第 j 无性系的小区平均值，为群体的总体平均值，Lh为第 h 地点效应，B(L)i(h)为第h 地点内第 i 重复效应，Cj为第 j 无性系效应，ILChj为第 j 无性系和第 h 地点的互作效应，ijk为剩余项。此外各性状的均方(MS)是由第类平方和与自由度(DF)的比值计算。偏效应量(ES，2)基于第类平方和计算，公式为：。其中，A 表示某个变异因子(或因素)。SSA表示某个变异因子的平方和，SSE表示误差平方和。基于线性混合模型，单地点以无性系、无性系与重复互作效应作为随机因子；多地点以无性系、无性系与地点互作效应作为随机因子，其他变异因子作为固定效应。获得的各性状方差分量用于计算遗传参数。无性系重复力(Rc)的计算公式21为：Rc=VgVg+VgeL+VeNL。表1参试的无性系Table1Clonestested树种和杂交组合()无性系赤桉SX95粗皮桉巨桉EC199粗皮桉尾叶桉EC195巨桉粗皮桉EC183巨桉巨桉EC184巨桉尾叶桉EC182尾叶桉赤桉EC189尾叶桉粗皮桉EC186、EC187、EC190、EC191、EC192、EC193、EC194、EC196、EC200尾叶桉巨桉EC181、DH32-29(ck)细叶桉粗皮桉EC198细叶桉巨桉EC188细叶桉尾叶桉EC180说明：赤桉E.camaldulensis，粗皮桉E.pellita，巨桉E.grandis，尾叶桉E.urophylla。表2各试验地基本情况Table2Basicinformationofeachtestarea试验地编码纬度(N)经度(E)土壤类型年均气温/年均降水量/mm广西南宁E1228210837赤红壤23.51350广西钦州E2219810866砖红壤24.01600广东雷州E3209211010沙壤23.21864广西扶绥E4226410791赤红壤24.01250海南定安E5196811036砖红壤26.01639说明：气候信息来源参考中国气象数据网(http:/)。第 40 卷第 5 期王楚彪等：桉树无性系生长性状基因型与环境互作效应及优良无性系选择953其中，Vg为无性系基因型方差，Vge为无性系与重复(或地点)交互作用方差，Ve为残差，L 为重复(或地点)个数，NL为重复(或地点)内各无性系个数，其中每个重复(或地点)内各无性系的个数不同，NL的值需使用调和平均数。1.3.3BLUP-GGE 分析模型为了鉴别各试点中表现最好的无性系、分析各试点对无性系评价的相似性、筛选理想选材环境和无性系，分别基于材积平均值和 BLUP 值进行 GGE 模型分析，该模型公式2425为：YijYij=1i1j1+2i2j2+ij。Yij其中，Yij为无性系 i 在环境 j 下的预测产量；为所有无性系 i 在环境 j 下的总体平均值；1和 2为主成分单值分解，i1和 i2是无性系 i 在主成分上的特征向量，j1和 j2为环境 j 在主成分上的特征向量；ij是随机误差。使用 R4.1.1、Rstudiov.2022.12.0+353、SPSS26、Excel2019 进行数据分析。2结果和分析2.1不同无性系生长表现各无性系的材积、材积变异系数及保存率的数据统计见表 3，单地点方差分析表明各无性系的材积差异极显著(P0.01,未列表)。5 个试验地相比，南宁的平均材积(0.1010m3株1)较其他地点大，定安的平均材积(0.0646m3株1)较小，说明定安对于参试无性系的适应性相对差。在各地中选择前 3 个优良无性系(除对照)，其中南宁优良无性系为 EC186、EC184 和 EC189；钦州为 EC184、EC183 和EC191；雷州为 EC192、EC186 和 EC188；扶绥为 EC188、EC184 和 EC183；定安为 EC180、EC187 和EC190。可见 EC184 在南宁、钦州和扶绥 3 地表现优良。各无性系(多地点分析)中，EC186、EC188、EC184 的 平 均 材 积 大 于 DH32-29。EC182、EC183、EC184 和 EC189 的 保 存 率 均 大 于 80.0%，但EC182 与 EC189 的材积表现略低。此外，EC183、E184、EC189 材积的变异系数(38.51%39.82%)相对较低，表明这些无性系稳定性相对较好；EC186 的材积变异系数(45.61%)则较高。2.2不同地区无性系重复力2.2.1不同地区联合方差分析由表 4 可见：对于树高，地点无性系的交互作用占总表型贡献率(即方差效应量)的 14.3%，地点的效应量为 22.0%，无性系的效应量为 17.0%；对于胸径，地点无性系的效应量为 12.2%，地点的效应量为 6.8%，无性系的效应量为 21.2%；对于材积，地点无性系的效应量为13.9%，地点的效应量为 8.7%，无性系的效应量为 20.7%。可见树高受到地点影响较其他性状大，胸径受地点影响较小。2.2.2无性系各性状在不同地点下的重复力基于线性混合模型，获得各性状(除保存率)的方差结果和遗传参数值，结果如表 5。树高、胸径、材积的无性系重复力分别为 0.71140.9186、0.67810.9346、0.64030.9207，说明各性状无性系之间的生长变异主要受遗传控制。广西南宁试验地材积无性系重复力较高，为 0.8615，胸径较低，为 0.7907，该地残差较其他地大。广西钦州试验地各性状的无性系重复力几乎相等(趋于 0.89)，基因型方差分量占较大比例，说明各无性系内单株之间生长变异受到遗传控制较高，不同重复、重复无性系之间的胸径性状差异不显著，表明不同重复的种植环境对胸径生长无影响。广东雷州试验地树高性状无性系重复力(0.7511)较低，可能其对环境更加敏感；而该地区无性系内胸径、材积性状重复力都较高。广西扶绥试验地各性状无性系重复力大于 0.90，可能受环境因素的影响较小，而且方差结果表明该试验地各性状的无性系效应量较其他变异来源高。海南定安试验地各性状重复无性系方差所占的分量比基因型大，导致无性系重复力相对较低，同时该地区的 GE 大小(通过重复无性系的方差分量占无性系与重复无性系两者之和的百分比来衡量)较其他地点高。5 个地点联合分析中，无性系与地点的互作方差比单一地点重复互作方差要高，说明地点很大程度影响无性系的生长且该混合模型更有合理性。此时性状之间无性系重复力由大到小为胸径、材积、树高。2.3不同无性系的基因型与环境互作对材积原始数据进行 GGE 双标图分析见图 1，PC1 与 PC2 共同解释总 G+GE 效应的 78.96%，954浙江农林大学学报2023 年 10 月 20 日图 1A 中多边形由离原点最远的无性系连接而成，通过原点到多边形各边的垂线将双标图分成若干扇区，每个扇形区域中顶角的品种是该扇区中表现最好的无性系，靠近原点的品种对环境不敏感。根据落点的扇区可将广西南宁、海南定安划为一组，其中 EC186 在该环境材积最高；广西钦州、扶绥划为一组，其中 EC184 在该环境材积最高；广东雷州划为一组，其中 EC183 在该环境材积最高。图 1A 中两地线段间夹角的余弦值是它们的相关系数，夹角小于 90表示正相关，说明两地对品种排序相似26。夹角较小说明试验点是重复设置的，去掉 1 个试验点不影响对品种的评价。各地点之间存在正相关，其中广西南宁和海南定安夹角最小，说明这 2 个环境对无性系的选择具有一致性。广西南宁和钦州线段较长，说明广西南宁和钦州有较强的无性系区分能力。然而，海南定安线段较短，说明海南定安的不同无性系之间的生长差异小。表3不同参试无性系生长性状Table3Growthtraitsofdifferentclonestested无性系材积/(m3株1)平均材积/(m3株1)变异系数/%保存率/%广西南宁广西钦州广东雷州广西扶绥海南定安EC1800.0712hi0.0467fg0.1019bcde0.0846def0.0813a0.085947.2068.5EC1810.0750ghi0.0539defg0.0779g0.0622hi0.0477g0.065043.9574.3EC1820.1184abcde0.0289gh0.0585h0.0505ij0.0722abc0.064245.5583.1EC1830.1163bcdef0.1168a0.0967def0.1119b0.0657cde0.098938.5283.8EC1840.1378ab0.1212a0.0977cdef0.1132b0.0697bc0.102138.8580.7EC1860.1519a0.0849bc0.1149ab0.0967cd0.0771ab0.104845.6166.7EC1870.1129bcdef0.0797bcde0.0777g0.0870cdef0.0797ab0.085133.4253.7EC1880.1040bcdefgh0.0665cdef0.1126bc0.1262a0.0644cdef0.103245.2257.7EC1890.1348abc0.0722bcdef0.0843fg0.0815efg0.0685bcd0.082338.4984.0EC1900.0946defghi0.0517efg0.0935defg0.0603hi0.0794ab0.075646.7878.5EC1910.1285abcd0.0991ab0.0866efg0.0837ef0.0565efg0.085444.5176.6EC1920.1132bcdef0.0846bc0.1276a0.0909cde0.0615cdef0.100049.0063.0EC1930.1105bcdefg0.0816bcd0.0850fg0.0831ef0.0561efg0.080550.9756.7EC1940.0797fghi0.0467fg0.0974cdef0.0751fg0.0654cde0.078237.4173.2EC1950.0796fghi0.0727bcdef0.0993cdef0.0485jk0.0705abc0.073946.2178.0EC1960.1000cdefgh0.0775bcde0.0964def0.0971c0.0714abc0.090142.9769.9EC1980.0837efghi0.0574cdef0.0929defg0.0693gh0.0614cdef0.075850.2166.7EC1990.0585ij0.0748bcdef0.0529hi0.0778fg0.0573defg0.063351.2573.3EC2000.0902efghi0.0790bcde0.1089bcd0.0816efg0.0531fg0.084240.7156.4SX950.0370j0.0142h0.0415i0.0380k0.0327h0.036651.6176.2DH32-29(ck)0.13470.09170.100332.9894.0平均0.10100.07050.07050.08100.0646说明：同列不同字母表示在0.05水平下差异显著。表示数据未检测。表4地点与无性系的方差汇总Table4Summaryofvariancesofsitesandclones变异因子自由度树高胸径材积SSMSES/%SSMSES/%SSMSES/%修正模型11010992.999.9*7491.968.1*2.1320.019*截距1927809.3927809.3*567606.6567606.6*25.96525.965*地点44411.01102.8*22.01116.6281.9*6.80.3780.095*8.7无性系193201.2168.5*17.04104.9215.5*21.21.0370.055*20.7地点无性系762617.334.4*14.32108.627.8*12.20.6440.008*13.9重复(地点)11763.569.4*4.7161.814.7*1.10.0720.007*1.8误差375915632.84.215244.14.13.9750.001总计3870954435.0590342.632.071修正后总计386926625.722736.06.107说明：SS、MS、ES分别表示第1类平方和、均方、方差效应量。*与*分别表示在0.01、0.05水平下F检验差异显著。第 40 卷第 5 期王楚彪等：桉树无性系生长性状基因型与环境互作效应及优良无性系选择955图 1B 中带箭头的直线为平均环境轴，直线上圆圈表示平均环境值，试验地线段和平均环境轴的角度越小，代表性越强27。其中广西钦州与平均环境轴角度小，且线段较长，表明广西钦州具有较强的代表性和区分力，有区分力却没代表性的地区可以淘汰稳定性差的品种，既有区分力又有代表性的地区才能选择高产稳产的品种28。即在广西钦州选择的高产稳产品种较为可靠；广西南宁适于无性系的快速生长，同时也淘汰不稳定品种。图 1B 中环境平均轴所指方向是无性系平均值的走向，无箭头的绿直线表示该性状的总平均值27。材积由大到小依次列出 EC184、EC186、EC183、EC192、EC191、EC188、EC189、EC196、EC193、EC200(只列出大于材积总平均值的无性系)。无性系与平均环境轴之间的垂线越长，表示无性系越不稳定，其中 EC184 具有较强的稳定性，即该无性系在各环境下生长表现差异不大；其次 EC182 的稳定性最低，说明该无性系在各环境下生长表现差异大，由表 3 数据：EC182 在广西南宁生长较好，广西钦州生长较差。综合选择 EC184 为高材积量、高稳定性的无性系。对材积的 BLUP 值进行 GGE 双标图分析可知(因结果基本相似图未列出)：主成分共同解释G+GE 效应的 80.68%。选出无性系 EC184 在广西南宁、广西钦州、广西扶绥、海南定安环境下材积的表现最优，EC188 在广东雷州环境下表现最好。除定安地区外，其他地区有较强的无性系区分能力。广西钦州(E2)与前文分析相同，为理想的材积试验环境。无性系材积由大到小依次为 EC184、EC183、EC186、EC192、EC188、EC191、EC196、EC189、EC187、EC193、EC200(只列出大于 BLUP 值的总平均值的无性系)。尽管 EC184 材积量最大，但 BLUP-GGE 模型得出 EC183 具有较强的材积稳定性，与前面分析有差异，因此综合选择 EC183 为高材积量、高稳定性的无性系。表5不同地点无性系各性状的方差结果和遗传参数Table5Varianceresultsandgeneticparametersofeachcharacterinclonesfromdifferentlocations性状方差结果/遗传参数地点广西南宁广西钦州广东雷州广西扶绥海南定安5地点结合树高重复13.409*143.413*101.418*8.696*无性系34.142*80.532*57.043*101.538*16.978*重复无性系9.511*17.439*8.793*5.235*误差5.3923.8944.7804.0381.310Ve5.38933.91454.78473.85871.32014.4072Vg1.76253.47430.88631.66310.38360.7533Vge0.56110.83520.33180.43210.9360Rc0.83050.88110.75110.91860.71140.7816胸径重复7.15030.401*15.180*6.115*无性系37.191*56.974*117.873*88.210*13.528*重复无性系6.4808.302*5.599*4.505*误差7.3484.7875.2483.1271.150Ve7.34514.86365.45033.17391.18424.1470Vg1.85282.94221.82721.42030.28061.0333Vge0.10710.13150.16470.36220.6725Rc0.79070.88830.93460.93780.67810.8685材积重复0.0010.014*0.006*0.002*无性系0.013*0.014*0.028*0.028*0.003*重复无性系0.0010.002*0.002*0.001*误差0.0020.0010.0010.0010.000Ve0.0017780.0010320.0013650.0008760.0002960.001058Vg0.0007380.0006400.0004410.0004696.26E050.000247Vge3.10E055.51E055.83E059.80E050.000219Rc0.86150.88660.92070.94130.64030.8333说明：广西南宁无性系、误差自由度分别为 19、304，其他地的重复、无性系、重复无性系、误差的自由度分别为广西钦州(2、19、25、325)；广东雷州(3、19、54、1149)；广西扶绥(3、19、54、1174)；海南定安(3、19、43、632)。Ve表示残差，Vg表示无性系基因型方差，Vge表示重复(或地点)与基因型互作方差，Rc表示无性系重复力。*与*分别表示在0.01、0.05 水平下 F 检验差异显著。表示缺少重复数据。956浙江农林大学学报2023 年 10 月 20 日0.050.0500.050.100.020.040.0600.020.040.060.08第 2 主成分(13.27%)0.020.040.0600.020.040.060.08第 2 主成分(13.27%)第 1 主成分(65.69%)第 1 主成分(65.69%)EC180EC181EC182EC183EC184EC186EC187EC188EC189EC190EC191EC192EC193EC194EC195EC196EC198EC199EC200SX95E1E2E3E4E5A00.050.10EC180EC181EC182EC183EC184EC186EC187EC188EC189EC190EC191EC192EC193EC194EC195EC196EC198EC199EC200SX95E1E2E3E4E5BA 表示试验点分组及无性系选择；B 表示试验点代表值和无性系稳定性。图中 E1E5 分别表示广西南宁、广西钦州、广东雷州、广西扶绥、海南定安地区。图1基于材积平均值的 GGE 分析图Figure1GGEanalysisdiagrambasedonaveragewoodvolume3讨论3.1不同因子对无性系生长的影响无性系生长不仅受到遗传控制，而且受许多环境因子的影响，如降水量、土壤类型、气温、风力等。本研究表明：南宁地区的无性系生长表现最好。ROCHA 等29研究表明：潮湿环境更有利于桉树无性系材积增长，然而广西南宁的年均降水量相比广西钦州、广东雷州、海南定安较低，意味着这 3 个点地的无性系本应该有较好的材积生长量，但是无性系材积都低于广西南宁，推断有其他因素影响这 3 个地区的无性系生长。本研究无性系品种以尾叶桉和巨桉为母本居多，不排除因品种因素导致广西南宁在较低年均降水量下有良好生长表现的可能29。另外，由于气候条件相差不大，还可能是土壤类型所致。广西南宁的土壤类型为赤红壤，广西钦州和海南定安的土壤为浅海沉积物发育的砖红壤，一般认为赤红壤比浅海沉积物砖红壤更有高导水性、持水和供水的能力30。其次，海南定安材积量较其他地小，可能与其年均气温较高、地理位置靠近赤道有关，这与 deARAUJO 等24认为桉树无性系接近赤道的生产力在下降的结果有一定关联。此外，广西钦州、广东雷州、海南定安沿海，常年受到台风的影响，也可能导致无性系生长受抑制。由此可见，无性系的生长是遗传与环境相互作用的结果。3.2各性状方差效应量及重复力无性系生长过程中，各性状在不同环境下表现程度有所不同，有的性状受环境影响较大，不同性状的重复力也有所差异。本研究表明：树高易受到环境影响，重复力较低。方差效应量是衡量各变异因子大小的指标，树高性状中地点的方差效应量(22.0%)较无性系、无性系与地点的交互作用大，说明树高主要受地点影响，而胸径、材积主要受遗传的影响，与 FREITAS 等31、dosSANTOS 等32的研究结果一致。本研究无性系各性状的重复力(0.64030.9346)处于高水平，说明桉树无性系生长性状受较高程度的遗传控制。多点结合的各性状重复力中，树高重复力较低，与朱显亮等33得出尾细桉无性系树高重复力高于胸径的结果不一致，可能是其研究仅涉及单个试验地，同时本研究单地点试验中广西南宁与海南定安的树高重复力较胸径高，说明单地点重复之间影响树高的环境因子较小，导致树高重复力高于胸径。可见重复力大小与环境因子有一定关联，影响树高的环境因子有待进一步研究。3.3基因型与环境互作本研究表明：GGE 模型得出 EC184 为高产稳产的无性系，BLUP-GGE 模型则为 EC183，2 个模型一致得出广西钦州为理想选材环境。倘若没有应用基因型与环境互作相关分析模型，以传统手段从多个地点综合选出优于 DH32-29 的优良无性系，那么结果可能与预期不一致。本研究优于 DH32-29 的无性系包括 EC184、EC186、EC188，但 EC186 的变异系数大，EC188 则保存率低。周家维等7在贵州地区中筛选出 6 个优良尾细桉优于 DH32-29，李光友等34在桂北地区筛选出优于 DH32-29 的无性系均为尾巨桉杂交所得，并认为尾巨桉杂交值得关注，但都局限于单一地点，单一地区选择的优良无性系并不能第 40 卷第 5 期王楚彪等：桉树无性系生长性状基因型与环境互作效应及优良无性系选择957说明该无性系在其他环境具有代表性。以往研究表明：桉树生长性状受显著的 GE 影响35，本研究亦如此，不同地点的无性系重复之间是由环境造成，每个重复甚至单株之间所生长的空间环境独立，然而本研究广西钦州胸径的重复无性系相互作用不显著，其他地区却显著，表明同一地点内往往会表现微环境、块状或梯度型的空间差异，而广西钦州的胸径则没有表现出环境差异。为了提高无性系遗传效应估算的准确性，有学者利用各个树体的行列数据进行空间分析36，前提是种植布局的行列需整齐。对于 GE 大小，本研究通过方差分量的比