Genotype × environment interaction on growth traits of Eucalyptus clones and selection of elite clones

参试无性系共21个，4.5年生，详见表1，以种植面积较大的无性系DH32-29作为对照。

采用完全随机区组设计，每个小区4行，每行5株，株行距为2 m×3 m，样地设置4个重复，1个重复内的各个无性系重复20次。于2016年6月种植，2020年12月调查树高、胸径、保存率。试验地分布在广西、广东及海南，详见表2。

树高 (H, m) 和胸径 (D, cm) 分别使用瑞典Vertex Ⅳ超声波林木测高仪和卷尺测量，计算单株材积 (V, m³·株⁻¹) ^[22]，公式为：V=D²H/30 000。各无性系保存率为活立木株数占总株数的百分比。表型材积变异系数(C_v)公式为：C_v=σ/μ×100%，其中σ表示无性系材积标准差，μ表示无性系材积平均值。

单地点方差分析使用线性混合模型 (linear mixed model)，模型^[23]为：${{Y}}_{{ijk}}{=}{\mu}{+}{{B}}_{{i}}{+}{{C}}_{{j}}{+}{{B}}_{{C}{ij}}{+}{{ \varepsilon }}_{{ijk}}$。其中Y_ijk为单点试验第i重复第j无性系第k单株的观测值，μ为群体的总体平均值，B_i为第i重复效应，C_j为第j无性系效应，B_Cij为第j无性系和第i重复的互作效应，ε_ijk为剩余项。其中以无性系、无性系与重复的互作效应为随机效应，获得单地点各无性系的BLUP值用于GGE分析。

多地点方差分析同样建立线性混合模型，计算各性状的方差分量和方差效应量，模型^[23]为：${{Y}}_{{hij}}{=}{\mu}{+} {{L}}_{{h}}{+}{{B}{(}{L}{)}}_{{i}{(}{h}{)}}{+}{{C}}_{{j}}{+}{{I}}_{{{\rm{LC}}hj}}{+}{{ \varepsilon }}_{{hij}}$。其中，Y_hij为单点试验第h地点内第i重复第j无性系的小区平均值，μ为群体的总体平均值，L_h为第h地点效应，B(L)_i(h)为第h地点内第i重复效应，C_j为第j无性系效应，I_LChj为第j无性系和第h地点的互作效应，ε_ijk为剩余项。此外各性状的均方 (MS) 是由第Ⅰ类平方和与自由度 (DF) 的比值计算。偏效应量 (ES，η²) 基于第Ⅰ类平方和计算，公式为：${{\eta}}_{\left({A}\right)}^{{2}}{=}\dfrac{{{S}}_{{{\rm{S}}}{A}}}{{{S}}_{{{\rm{S}}}{A}}{+}{{S}}_{{{\rm{S}}}{{\rm{E}}}}}$。其中，A表示某个变异因子 (或因素)。S_SA表示某个变异因子的平方和，S_SE表示误差平方和。

基于线性混合模型，单地点以无性系、无性系与重复互作效应作为随机因子；多地点以无性系、无性系与地点互作效应作为随机因子，其他变异因子作为固定效应。获得的各性状方差分量用于计算遗传参数。无性系重复力(R_c)的计算公式^[21]为：

其中，V_g为无性系基因型方差，V_ge为无性系与重复 (或地点)交互作用方差，V_e为残差，L为重复 (或地点) 个数，N_L为重复 (或地点) 内各无性系个数，其中每个重复 (或地点) 内各无性系的个数不同，N_L的值需使用调和平均数。

为了鉴别各试点中表现最好的无性系、分析各试点对无性系评价的相似性、筛选理想选材环境和无性系，分别基于材积平均值和BLUP值进行GGE模型分析，该模型公式^[24−25]为：

其中，Y_ij为无性系i在环境j下的预测产量；$\overline Y_{ij} $为所有无性系i在环境j下的总体平均值；λ₁和λ₂为主成分单值分解，ξ_i1和ξ_i2是无性系i在主成分上的特征向量，η_j1和η_j2为环境j在主成分上的特征向量；ε_ij是随机误差。

使用R 4.1.1、Rstudio v.2022.12.0+353、SPSS 26、Excel 2019进行数据分析。

各无性系的材积、材积变异系数及保存率的数据统计见表3，单地点方差分析表明各无性系的材积差异极显著 (P＜0.01, 未列表)。5个试验地相比，南宁的平均材积 (0.101 0 m³·株⁻¹) 较其他地点大，定安的平均材积 (0.064 6 m³·株⁻¹) 较小，说明定安对于参试无性系的适应性相对差。在各地中选择前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%) 则较高。

由表4可见：对于树高，地点×无性系的交互作用占总表型贡献率 (即方差效应量) 的14.3%，地点的效应量为22.0%，无性系的效应量为17.0%；对于胸径，地点×无性系的效应量为12.2%，地点的效应量为6.8%，无性系的效应量为21.2%；对于材积，地点×无性系的效应量为13.9%，地点的效应量为8.7%，无性系的效应量为20.7%。可见树高受到地点影响较其他性状大，胸径受地点影响较小。

基于线性混合模型，获得各性状 (除保存率) 的方差结果和遗传参数值，结果如表5。树高、胸径、材积的无性系重复力分别为0.711 4~0.918 6、0.678 1~0.934 6、0.640 3~0.920 7，说明各性状无性系之间的生长变异主要受遗传控制。

广西南宁试验地材积无性系重复力较高，为0.861 5，胸径较低，为0.7907 ，该地残差较其他地大。广西钦州试验地各性状的无性系重复力几乎相等 (趋于0.89)，基因型方差分量占较大比例，说明各无性系内单株之间生长变异受到遗传控制较高，不同重复、重复×无性系之间的胸径性状差异不显著，表明不同重复的种植环境对胸径生长无影响。广东雷州试验地树高性状无性系重复力 (0.751 1) 较低，可能其对环境更加敏感；而该地区无性系内胸径、材积性状重复力都较高。广西扶绥试验地各性状无性系重复力大于0.90，可能受环境因素的影响较小，而且方差结果表明该试验地各性状的无性系效应量较其他变异来源高。海南定安试验地各性状重复×无性系方差所占的分量比基因型大，导致无性系重复力相对较低，同时该地区的G×E大小 (通过重复×无性系的方差分量占无性系与重复×无性系两者之和的百分比来衡量) 较其他地点高。

5个地点联合分析中，无性系与地点的互作方差比单一地点重复互作方差要高，说明地点很大程度影响无性系的生长且该混合模型更有合理性。此时性状之间无性系重复力由大到小为胸径、材积、树高。

对材积原始数据进行GGE双标图分析见图1，PC1与PC2共同解释总G+G×E效应的78.96%，图1A中多边形由离原点最远的无性系连接而成，通过原点到多边形各边的垂线将双标图分成若干扇区，每个扇形区域中顶角的品种是该扇区中表现最好的无性系，靠近原点的品种对环境不敏感。根据落点的扇区可将广西南宁、海南定安划为一组，其中EC186在该环境材积最高；广西钦州、扶绥划为一组，其中EC184在该环境材积最高；广东雷州划为一组，其中EC183在该环境材积最高。

Figure 1.  GGE analysis diagram based on average wood volume

图1A中两地线段间夹角的余弦值是它们的相关系数，夹角小于90°表示正相关，说明两地对品种排序相似^[26]。夹角较小说明试验点是重复设置的，去掉1个试验点不影响对品种的评价。各地点之间存在正相关，其中广西南宁和海南定安夹角最小，说明这2个环境对无性系的选择具有一致性。广西南宁和钦州线段较长，说明广西南宁和钦州有较强的无性系区分能力。然而，海南定安线段较短，说明海南定安的不同无性系之间的生长差异小。

图1B中带箭头的直线为平均环境轴，直线上圆圈表示平均环境值，试验地线段和平均环境轴的角度越小，代表性越强^[27]。其中广西钦州与平均环境轴角度小，且线段较长，表明广西钦州具有较强的代表性和区分力，有区分力却没代表性的地区可以淘汰稳定性差的品种，既有区分力又有代表性的地区才能选择高产稳产的品种^[28]。即在广西钦州选择的高产稳产品种较为可靠；广西南宁适于无性系的快速生长，同时也淘汰不稳定品种。

图1B中环境平均轴所指方向是无性系平均值的走向，无箭头的绿直线表示该性状的总平均值^[27]。材积由大到小依次列出EC184、EC186、EC183、EC192、EC191、EC188、EC189、EC196、EC193、EC200 (只列出大于材积总平均值的无性系)。无性系与平均环境轴之间的垂线越长，表示无性系越不稳定，其中EC184具有较强的稳定性，即该无性系在各环境下生长表现差异不大；其次EC182的稳定性最低，说明该无性系在各环境下生长表现差异大，由表3数据：EC182在广西南宁生长较好，广西钦州生长较差。综合选择EC184为高材积量、高稳定性的无性系。

对材积的BLUP值进行GGE双标图分析可知 (因结果基本相似图未列出)：主成分共同解释G+G×E效应的80.68%。选出无性系EC184在广西南宁、广西钦州、广西扶绥、海南定安环境下材积的表现最优，EC188在广东雷州环境下表现最好。除定安地区外，其他地区有较强的无性系区分能力。广西钦州 (E2) 与前文分析相同，为理想的材积试验环境。无性系材积由大到小依次为EC184、EC183、EC186、EC192、EC188、EC191、EC196、EC189、EC187、EC193、EC200 (只列出大于BLUP值的总平均值的无性系)。尽管EC184材积量最大，但BLUP-GGE模型得出EC183具有较强的材积稳定性，与前面分析有差异，因此综合选择EC183为高材积量、高稳定性的无性系。

无性系生长不仅受到遗传控制，而且受许多环境因子的影响，如降水量、土壤类型、气温、风力等。本研究表明：南宁地区的无性系生长表现最好。ROCHA等^[29]研究表明：潮湿环境更有利于桉树无性系材积增长，然而广西南宁的年均降水量相比广西钦州、广东雷州、海南定安较低，意味着这3个点地的无性系本应该有较好的材积生长量，但是无性系材积都低于广西南宁，推断有其他因素影响这3个地区的无性系生长。本研究无性系品种以尾叶桉和巨桉为母本居多，不排除因品种因素导致广西南宁在较低年均降水量下有良好生长表现的可能^[29]。另外，由于气候条件相差不大，还可能是土壤类型所致。广西南宁的土壤类型为赤红壤，广西钦州和海南定安的土壤为浅海沉积物发育的砖红壤，一般认为赤红壤比浅海沉积物砖红壤更有高导水性、持水和供水的能力^[30]。其次，海南定安材积量较其他地小，可能与其年均气温较高、地理位置靠近赤道有关，这与de ARAUJO等^[24]认为桉树无性系接近赤道的生产力在下降的结果有一定关联。此外，广西钦州、广东雷州、海南定安沿海，常年受到台风的影响，也可能导致无性系生长受抑制。由此可见，无性系的生长是遗传与环境相互作用的结果。

无性系生长过程中，各性状在不同环境下表现程度有所不同，有的性状受环境影响较大，不同性状的重复力也有所差异。本研究表明：树高易受到环境影响，重复力较低。方差效应量是衡量各变异因子大小的指标，树高性状中地点的方差效应量 (22.0%) 较无性系、无性系与地点的交互作用大，说明树高主要受地点影响，而胸径、材积主要受遗传的影响，与FREITAS等^[31]、dos SANTOS等^[32]的研究结果一致。本研究无性系各性状的重复力 (0.640 3~0.934 6) 处于高水平，说明桉树无性系生长性状受较高程度的遗传控制。多点结合的各性状重复力中，树高重复力较低，与朱显亮等^[33]得出尾细桉无性系树高重复力高于胸径的结果不一致，可能是其研究仅涉及单个试验地，同时本研究单地点试验中广西南宁与海南定安的树高重复力较胸径高，说明单地点重复之间影响树高的环境因子较小，导致树高重复力高于胸径。可见重复力大小与环境因子有一定关联，影响树高的环境因子有待进一步研究。

本研究表明：GGE模型得出EC184为高产稳产的无性系，BLUP-GGE模型则为EC183，2个模型一致得出广西钦州为理想选材环境。倘若没有应用基因型与环境互作相关分析模型，以传统手段从多个地点综合选出优于DH32-29的优良无性系，那么结果可能与预期不一致。本研究优于DH32-29的无性系包括EC184、EC186、EC188，但EC186的变异系数大，EC188则保存率低。周家维等^[7]在贵州地区中筛选出6个优良尾细桉优于DH32-29，李光友等^[34]在桂北地区筛选出优于DH32-29的无性系均为尾巨桉杂交所得，并认为尾巨桉杂交值得关注，但都局限于单一地点，单一地区选择的优良无性系并不能说明该无性系在其他环境具有代表性。

以往研究表明：桉树生长性状受显著的G×E影响^[35]，本研究亦如此，不同地点的无性系重复之间是由环境造成，每个重复甚至单株之间所生长的空间环境独立，然而本研究广西钦州胸径的重复×无性系相互作用不显著，其他地区却显著，表明同一地点内往往会表现微环境、块状或梯度型的空间差异，而广西钦州的胸径则没有表现出环境差异。为了提高无性系遗传效应估算的准确性，有学者利用各个树体的行列数据进行空间分析^[36]，前提是种植布局的行列需整齐。对于G×E大小，本研究通过方差分量的比值来权衡，得出海南定安各性状的G×E最大，各性状中树高G×E较大。同样地，OSORIO等^[37]通过两两地点的B型遗传相关 (实际上是方差分量的比值) 表明巨桉的树高G×E较材积、木材密度大。

GGE分析中，基于BLUP数据的GGE模型比原始数据具有更高的变异解释能力，因此BLUP-GGE模型更为可靠，这与程玲等^[38]的研究一致。为了解决试验地环境异质性问题并使分析结果更可靠，程玲等^[38]基于平均值、BLUP值与GGE双标图联合分析火炬树Rhus typhina的基因型与环境互作效应，发现原始平均值与BLUP值对应的理想选材环境不同，但是本研究平均值和BLUP值对应的理想选材环境都是广西钦州，可能是前者多考虑了空间分析造成选材环境有所差异。其次，基于平均值、BLUP值的GGE模型分别筛选出EC184、EC183为高材积量、高稳定性品种，2个无性系都是以巨桉为母本，说明巨桉在生长方面具有较强的基因稳定性以及适应能力，李光友等^[34]也认为巨桉为母本具有良好的表现。

本研究综合5个地点的数据，桉树各无性系的生长表现存在显著差异。5个地点联合方差分析及遗传参数中，胸径的地点效应量 (6.8%) 较其他性状低，胸径的无性系重复力较其他性状高，说明胸径性状较其他更稳定，受到遗传控制较高。

基于材积平均值的GGE分析，广西钦州为理想选材环境，EC184为材积较高且稳定的无性系；基于材积BLUP值的GGE分析，广西钦州同样为材积理想选材环境，EC183为材积较高且稳定的无性系。所以GGE模型筛选出EC183和EC184为高材积和高稳定性的无性系。