《Evolutionary Applications》：Modeling Phenotypic Trait Variation and Plasticity in Elymus elymoides to Guide Climate-Informed Seed Transfer
编辑推荐：
这篇综述以重要生态恢复草种Elymus elymoides（瓶刷松鼠尾草）为研究对象，通过三个同质园实验，系统评估了其生长、繁殖、形态、物候和存活等表型性状的定量遗传变异，并揭示了其与种源气候的复杂关联。研究创新性地结合多元回归树（multivariate regression trees）和随机森林（random forests）等非标准建模方法，构建了具有明确层级结构的固定边界种子区（fixed-boundary seed zones）和动态的焦点种子转移模型（focal-point seed transfer models），为全球变化背景下制定灵活、精确的种子调运（seed-sourcing）决策提供了新框架。
本文的研究对象是一种名为瓶刷松鼠尾草（Elymus elymoides (Raf.) Swezey）的禾本科植物。这种植物是美国西部山间地区（Intermountain Region）生态恢复中的重要物种，常被称为“主力”草种，广泛应用于野火后植被恢复等项目。它耐干扰、能一定程度耐受火烧，并对入侵性草种有潜在的竞争力。瓶刷松鼠尾草分布范围广泛，生境多样，至少被鉴定出五个地理上重叠的亚种，并且此前的研究已证实了其在亚种和种群水平上存在性状变异，且这种变异与种源气候相关联，对恢复效果有潜在影响。
本研究核心目标是刻画该物种在核心分布区内的定量遗传变异，并利用性状与种源气候之间的关联，开发一套基于实证的种子转移指导模型，以应对全球变化的挑战。研究团队在2007至2010年间，从美国西部五个州的98个野生种群收集了种子。种子先在生长室中萌发，并在俄勒冈州立大学的温室中培养，随后在2011年秋季移植到三个位于华盛顿州、俄勒冈州和内华达州的同质园（Central Ferry, Powell Butte, Reno）。每个种群在两个母系家庭中各取一个个体，随机分配到每个花园的六个区块中，每个花园共种植了1176株测试植物，总计3528株。为了最小化边缘效应，每个花园的边缘种植了一圈由多余植株构成的缓冲区。
在2013年的生长季，研究人员系统测量了与植物生长、叶片形态、物候、繁殖和存活相关的多个性状，包括地上生物量、株高、花序数量、叶片长度和宽度、叶片长宽比、抽穗期、种子成熟期以及植株存活率。他们利用从ClimateWNA获取的种源地30年（1981-2010）气候常态数据，来分析与性状的关联。通过方差分量模型分析，研究发现，尽管大部分性状变异来自花园效应（即表型可塑性），但种群水平的变异在绝大多数性状中依然显著，是后续建模的基础。然而，叶片长宽比在Reno花园的种群水平变异可忽略不计，因此在后续分析中被排除。
研究还量化了种群水平的性状可塑性，并将其与种源气候进行了相关性分析。结果显示，生物量、株高、叶片宽度和存活率的可塑性通常与一个或多个降水变量呈正相关，与表明气温更高、无霜期更长或更干旱的变量呈负相关。而种子成熟期可塑性的模式则完全相反。这表明来自不同气候的种群，其性状对环境变化的响应能力（即可塑性）存在系统性差异。
研究的核心创新在于其种子转移建模方法。传统的标准方法通常基于主成分分析（Principal Component Analysis, PCA）降维后，利用多元线性或多项式回归建立转移函数，再进行空间投影和分区。本文则采用了非标准方法，将多元回归树用于定义固定边界种子区，将随机森林用于构建焦点种子转移模型。
首先，研究人员对26个性状-花园组合（9个性状×3个花园，减去被排除的叶片长宽比）的数据矩阵进行主成分分析，得到了6个显著的主成分轴，共同解释了77%的总方差。这6个轴代表了性状在不同花园中表达的复杂模式对比，例如轴1上多个花园的晚花、晚熟、高生物量和高株高与Reno花园的高存活率形成对比。
固定边界种子区：研究人员以6个主成分轴为响应变量，以种源气候变量为预测变量，建立了一个多元回归树模型。该模型最终产生13个终端节点（即13个种子区），每个节点包含5-17个种群。这个决策树模型呈现了一个清晰的层级结构：第一层由季节性温差（TD）和最暖月均温（MWMT）划分出三个主区（A, B, C）；每个主区进一步被降水、温度、霜冻等多种气候变量划分为2-3个二级区；部分二级区还可再细分为三级区。这三个层级的划分分别解释了模型22%、14%和12%的方差。将这些分区空间投影到研究区域，就形成了一套具有层级结构的种子区划图。例如，A区（低季节性温差）主要位于研究区西部的内华达山脉、喀斯喀特山脉等地；B区和C区则大致沿内陆山脉从高海拔、较冷湿气候向低海拔、较温暖干旱气候呈梯度分布。来自较温和气候区（如低季节性温差、较凉爽或较湿润）的种群，其植株通常更大、花序更多、叶片更宽、物候更晚；而来自较严酷气候区的种群则相反。
焦点种子转移模型：研究人员为上述6个主成分轴分别构建了随机森林模型，以建立每个性状轴与气候变量之间的关联函数。模型揭示了不同性状轴所关联的最重要气候变量各不相同，例如轴1与季节性温差、最暖月均温和冷却度日（DD18）关系密切，轴2则与冷度日（DD_0）、无霜期结束日期和年均降水量更相关。通过部分依赖图可以观察这些关键气候变量的边际效应。将这些随机森林模型空间投影，就得到了每个主成分轴的连续分布图。将所有这些轴的预测值结合起来，可以计算地图上任意两点在性状空间中的多维欧几里得距离，从而动态评估任意一个“焦点”（如一个计划进行生态恢复的地点）与整个区域内所有潜在种源之间的性状相似性。距离越近，表示性状越相似，种子转移的适宜性理论上越高。这构成了一个灵活的、无预设阈值的焦点种子转移模型。
此外，研究还探讨了亚种分化和性状可塑性对种子转移指导的潜在影响。在分析的98个种群中，鉴定出四个亚种：ssp. elymoides(57个)、ssp. brevifolius(16个)、ssp. californicus(9个) 和 ssp. hordeoides(8个)。当在多元回归树模型中加入“亚种”作为一个分类预测变量时，模型的第一级划分就基于亚种身份，将ssp. brevifolius和 ssp. californicus与其他亚种区分开。这导致为这两个亚种组单独建立了与之前不同的种子区模型，表明不同亚种组的性状-气候关联模式存在明显差异。因此，在制定种子转移指南时，考虑种内谱系分化是必要的。
综上所述，本研究通过创新的建模框架，成功刻画了瓶刷松鼠尾草复杂的表型变异及其与气候的关联，并开发了兼具实用性（层级化、可调整的固定边界种子区）和精确性（连续的焦点转移模型）的种子转移指南。研究结果强调了性状-气候关联的多面性，并指出在全球变化时代，将亚种分化、性状可塑性等进化因素纳入考量的种子转移建模，对于做出明智的种子调运决策至关重要。这套方法不仅适用于瓶刷松鼠尾草，也为其他生态恢复物种的适应性管理提供了可借鉴的范例。
