基因改良火炬松生长表现与基因型特异性气候敏感性——基于西湾海岸平原20年遗传试验的树轮证据

《Journal of Forestry》：Climate Sensitivity and Growth Response of Genetically Improved Loblolly Pine: A Case Study from a 20-Year Genetic Trial in the West Gulf Coastal Plain

　　

识别基因型特异性气候敏感性可为改良火炬松（Pinus taeda L.）的选择和部署策略提供指导。研究人员利用2004年建立的20年生火炬松遗传试验中的增量芯，重建年基面积增量（basal area increment，BAI），并评估美国南部干旱易发区西湾海岸平原（West Gulf Coastal Plain，WGCP）不同遗传材料生长与时间气候参数之间的关系。基因型包括3个开放授粉（open-pollinated，OP）材料（OP 7–56、OP 8–103和OP LA）和2个品种（V9和V93），均种植于同一地点。OP 7–56在20年生时表现出最高的平均BAI（17.3 cm²），生长峰值出现在6–10年；V9和V93则表现出更中等但稳定的生长轨迹。响应函数分析表明，V9和V93对当年气候最敏感，尤其是春季降水和夏末至秋季干旱指数；OP 7–56和OP 8–103对前一年初夏温度响应更强；OP LA的气候信号最弱。采用错误发现率（false discovery rate，FDR）程序控制多重比较后，关键基因型的生长—气候关系仍可检测，提供了差异气候敏感性的初步证据，值得进一步研究。小波分析揭示持续的年际（2–4年）和多年（3–6年）生长周期，以及基因型对较长时间尺度气候变率的特异性敏感性。总体而言，改良火炬松基因型可能在生长表现和时间气候敏感性上存在差异，在提高南部松树人工林对气候变化韧性时需加以考虑，尤其是在西湾海岸平原。

研究背景与目的：火炬松（Pinus taeda L.）是美国南部重要商品林树种，长期树木改良已形成开放授粉（open-pollinated，OP）家系、控制授粉家系和克隆品种等不同遗传材料，其生长潜力与环境胁迫抗性存在差异。现有问题在于，南部地区干旱频率和高温风险增加，而木材生产、碳汇和可再生生物质需求同时上升；WGCP是火炬松栽培中最干旱、干旱易发区域，土壤水分有限且气候变率大，基因型×环境互作更明显。已有树轮气候研究增进了对生长—气候关系的认识，但将其用于评价不同改良遗传材料气候敏感性仍有限；小波分析虽适于非平稳、多尺度气候信号检测，但在树轮气候研究中应用仍不足。为此，研究人员利用2004—2024年树轮数据，评估WGCP改良火炬松基因型、OP家系和克隆品种对月及季节降水、温度和干旱指数的BAI响应差异，以检验气候—生长敏感性是否因基因型而异，并展示无损树轮分析作为基因型筛选和部署决策工具的价值。研究得出，改良基因型在生长表现和时间气候敏感性上存在差异，需在提升南部松树人工林气候韧性时加以考虑。

关键技术方法：研究地点位于美国路易斯安那州Claiborne Parish的Hill Farm Research Station（LSU AgCenter），属WGCP；2004年建立20年生遗传试验，8.1 ha、3个随机完全区组、54个0.1 ha样地，含OP 7–56、OP 8–103、OP LA、V9和V93，OP LA仅Block 1（3个样地）作区域对照。2024年测DBH和树高，按径阶每样地选3株（OP LA 5株/样地），共取168根生长芯；经固定、打磨、1200 dpi扫描，WinDENDRO测年轮宽，dplR交叉定年并构建BAI年表；采用PRISM 800 m气候数据与SPEI 3/6/9月指标，treeclim响应相关和FDR校正，计算抗旱性、恢复力和韧性指数，并进行连续小波变换。

研究结果：
1）基面积增量随年龄变化：基于2006—2024年BAI和累积BA分析，所有基因型呈S形生长，5—9年生长最快，12—20年下降；OP 7–56生长最高，20年累积BA近400 cm²，但个体变异大；V9和V93较中等且稳定，V93年BAI随时间最低。Welch方差分析显示基因型对年BAI有显著影响，OP 7–56显著高于V9和V93，但与OP 8–103、OP LA无显著差异。OP家系变异系数多高于50%，V9最稳定（28.1%—38.6%），OP 7–56和OP LA成熟期仍高于85%。结论是快速生长与稳定性存在基因型差异。
2）不同基因型BAI年表：基于去趋势BAI年表统计，克隆品种年表一致性较高（V93和V9的Rbar分别为0.68和0.66），OP 7–56对前一年生长依赖最弱；各基因型EPS均≥0.95，克隆品种SNR较高，OP LA最低；OP LA平均敏感性最高（0.46），Gini系数最高（0.35），MSI以V93（0.52）和OP 8–103（0.50）较高，OP 7–56最低（0.37）。结论是克隆品种年表更一致稳定，OP材料变异更大。
3）研究期气候趋势：基于Mann-Kendall趋势检验，年降水无显著趋势（1453.62±345.10 mm，τ=0.11，p=0.49），但年际变率大，2010年极端干旱（872.23 mm）；年均温呈边缘显著正趋势（τ=0.27，p=0.097，17.88±0.66 °C）。SPEI显示2011年干旱最严重，尤其2010—2012年；2004、2016和2019年水分盈余明显，2017年短期干旱加剧。结论是长期增温可能较降水更突出，干旱变率显著。
4）基因型特异性生长—气候关系：基于treeclim响应函数相关和FDR校正，V9和V93对当年气候更敏感，尤其3月降水、V93的8月温度和夏末至秋季SPEI；OP家系对前一年6月平均和最高温度响应较强；OP LA除前一年6月均温外无显著相关，且样本较小、统计功效有限。FDR校正后名义显著关系未保持显著，应作探索性解释。基于2011和2017年干旱事件的Rt、Rc、Rs指数，基因型间无显著差异；2011年Rt、Rc、Rs多大于1，2017年Rt和Rs多小于1。结论是克隆品种对当年水分和热胁迫更敏感，OP材料对前一年温度更敏感，事件尺度差异不显著。
5）不同基因型径向生长的多尺度周期性：基于标准化年轮指数RWI连续小波变换，4个基因型均呈持续年际（约2—4年）和多年（约3—6年）周期，早期以短周期为主，2008—2018年3—6年信号最强；OP家系6—8年较长周期更明显。RWI约2015年达峰后下降，2018年后小波功率减弱。结论是径向生长受短期干扰与多年水分/胁迫恢复过程共同影响，OP材料对低频气候变率更敏感。

讨论总结：讨论认为，OP 7–56早期生长最快但气候敏感性较高、树间变异较大，在干旱易发地风险较高；V9和V93峰值BAI较OP 7–56低11%—12%，但年际稳定性和树间一致性更高，体现快速生长与长期韧性之间的权衡。苗木类型可能部分解释差异：OP为裸根苗，克隆品种为容器苗，容器苗早期根系生长潜力较大，但该效应具有地点依赖性，不能仅归因于遗传。单点试验应视为案例研究，区域外推需谨慎。气候敏感性方面，克隆品种对当年水分和干旱指数较敏感，V93对晚季高温表现为负响应；OP材料对前一年6月温度响应较强，可能涉及储存碳水化合物的滞后效应。FDR校正后信号减弱，且PRISM 800 m数据难以解析样地微气候，故需谨慎解释。小波结果说明生长受短周期扰动和多年气候过程共同控制，2015年后RWI下降及2018年后功率减弱可能与累积生理胁迫、碳平衡限制、竞争及大小/年龄限制有关；2011年抗旱指数较高可能反映幼龄快速生长，而非单纯耐旱。管理上，OP家系适合短期生长最大化且干旱风险较低或有适应管理的地点，克隆品种适合土壤水分较可靠或采用降低水分胁迫措施的地点；维持OP与品种材料组合有助于分散风险、保持遗传多样性和长期生产力。

研究结论（翻译）：本研究考察了干旱易发的WGCP地区基因改良火炬松的生长—气候关系，该区域生产力和气候韧性差异显著，并关系到可持续商品林经营。在给定地点，没有单一基因型能同时最大化生产力、生长稳定性和对气候极端事件的抗性，这些权衡很可能随地点变化，提示应采用气候智慧型组合部署。OP家系（如OP 7–56）可用于具备适应性管理条件的地点；克隆品种（V9和V93）生长较保守，适合优先追求稳定生长、较低经营风险和较少韧性干预的地点。由于OP LA样本量较小，其推论为初步结果，尚不宜指导生产部署。小波分析识别出持续的年际和多年生长周期，表明径向生长不仅受逐年气候条件影响，也受较长期水分和胁迫动态影响。这些基因型—气候关系为降低可持续经营的生物和经营风险提供依据：将耐旱或气候缓冲型基因型匹配到水分受限地点，可降低多年干旱下生长崩溃风险、稳定碳汇并减少以高成本经营投入“抢救”林分的压力；将种植密度和间伐计划与基因型特异性生长轨迹和气候敏感性匹配，有助于将BAI维持在限制水分胁迫和死亡率的阈值内。总体而言，研究支持将遗传选择纳入更广泛的可持续规划，部署OP家系和克隆品种材料的混合可分散遗传基础、缓冲未来气候不确定性，并保持木材和碳目标的灵活选择。未来应延长年表长度以增强统计功效，并改变地点条件验证本研究的探索性信号；将小波衍生的气候敏感性纳入生长与产量模型和部署指南，将有助于土地所有者和经营者在气候变率加剧下建立和经营高产、具韧性的火炬松林分。