矮鸢尾不同生境种群叶片性状对实验增温响应的种内变异研究

《Plants》：Intraspecific Variation in Leaf Trait Responses to Experimental Warming in Iris pumila Populations from Contrasting Habitats

　　

摘要 来自不同生境的植物种群对气候变暖的表型响应可能存在差异。为评估矮鸢尾（Iris pumila）叶片性状对增温响应的种内变异，研究人员在受控条件下，将来自沙丘（Dune）和树林（Woods）种群的40个基因型依次暴露于环境温度、+1 °C和+4 °C的温度处理。功能叶片性状（比叶面积，SLA；叶片干物质含量，LDMC）和机制叶片性状（比叶含水量，SLWC；叶厚，LT，由SLA和LDMC估算；气孔密度，SD；以及RGB衍生指数I1，此前已在矮鸢尾中校准为与叶绿素浓度呈负相关）随温度处理序列发生变化。SLA、SLWC、LT和I1升高，而LDMC和SD降低。沙丘与树林基因型在SLA、LDMC和I1上存在差异，其中沙丘基因型的SLA和I1较高，树林基因型的LDMC较高。表型可塑性指数（PIV）在每个温度对比中因叶片性状而异，总体上SD和I1的值较高，而SLA的值较低。种群间PIV的差异在环境温度至+1 °C的LDMC和SLWC（沙丘 > 树林）以及+1 °C至+4 °C的LDMC（树林 > 沙丘）中表现明显。性状相关性结构随温度序列而变化，尽管整体相关矩阵相似性的证据对性状组成的改变并不稳健。在不同温度处理或种群之间，未检测到表型整合程度的统计学显著差异。总体而言，矮鸢尾的叶片响应揭示了跨温度处理序列的性状特异性变化和种群依赖性模式，反映了种内对增温响应的复杂性。

论文解读：矮鸢尾不同生境种群叶片性状对实验增温的响应及其种内变异

一、研究背景与科学问题

全球气候变暖正持续改变自然植物种群所经历的热环境条件。2011—2020年间全球地表温度较工业化前水平升高约1.1 °C，且未来数十年在温室气体持续排放情景下预计将进一步升温。温度变化会影响植物代谢、生长与发育，进而改变物候、生理表现和生产力。由于植物固着生长，其原位表型调整能力尤为重要。表型可塑性——即同一基因型在不同环境条件下表达不同表型的能力——可在种群通过遗传进化作出响应之前提供快速的代内响应。然而，可塑性并不天然具有适应性，其存在或大小本身不能作为适应能力的证据。

叶片性状为研究植物对环境变化的响应提供了有价值的框架。比叶面积（SLA）和叶片干物质含量（LDMC）作为功能性状，反映了资源获取与组织投资策略之间的权衡。然而，这些功能性状的整合性可能限制其解析植物对气候变化更直接响应组分的能力。因此，研究者将SLA和LDMC视为功能性状，而将比叶含水量（SLWC）、估算叶厚（LT）、气孔密度（SD）和RGB衍生叶绿素相关指数（I₁）视为机制性状，分别代表叶片表型的水合—结构、气孔—解剖和叶绿素相关组分。

除单个性状的响应外，功能性状与机制性状可能因共享功能、结构或发育关系而协同变化。表型整合——即有机体内功能相关性状间的协变——为检验这些关联提供了框架。升温不仅可能影响单个性状值和可塑性，还可能影响性状间的关联模式。此外，物种内不同种群间的性状表达和可塑性可能存在差异，来自对比环境的种群可能在平均性状表达和表型可塑性上有所不同。此前研究仅包含单一来源种群和两个温度条件，尚不清楚中等增温的响应是否持续或随温度进一步升高而变化，以及其幅度和协调性是否在不同生境种群间存在差异。本研究旨在通过联合考虑表型可塑性、种群来源和性状整合，更全面地理解叶片响应。

二、研究设计与主要技术方法

本研究选取来自塞尔维亚Deliblato沙地两个对比生境的矮鸢尾自然种群：开阔沙丘种群（Dune）和遮荫林地种群（Woods），每个种群选取20个克隆基因型，共40个基因型。实验材料在受控生长室中维持三年以最小化母体环境效应。三个温度处理（环境温度、+1 °C和+4 °C）在同一生长室中依次施加，光周期、光源、基质、浇水计划保持一致。叶片性状测量包括：SLA（叶面积/干质量）、LDMC（干质量/鲜质量）、SLWC（鲜质量与干质量之差/叶面积）、LT（由SLA和LDMC估算）、SD（表皮印迹法计数）以及RGB指数I₁（R + G ? 2B）。统计方法包括重复测量方差分析、剖面分析、表型可塑性指数（PI_V）计算、Friedman ANOVA、Mann–Whitney U检验、非参数纵向分析、Pearson相关分析、Mantel检验、主成分分析和置换检验。

三、主要研究结果

3.1 功能与机制叶片性状随温度序列的响应

跨温度序列，SLA、SLWC、LT和I₁总体升高，LDMC和SD总体下降。SLA在+4 °C时两个种群均较环境温度增加约9%，LDMC分别下降15%和16%。重复测量方差分析显示所有性状均随温度序列显著变化。种群效应对SLA（ηp2 = 0.385）、LDMC（ηp2 = 0.255）和I₁（ηp2 = 0.313）显著。温度×种群交互作用对SLWC和LT显著，表明这两个机制性状的响应剖面在种群间非平行。

3.2 表型可塑性跨温度序列的变化

反应规范图显示SLA、SLWC、LT和I₁总体呈上升趋势，LDMC呈下降趋势，SD表现出先升后降的复杂模式。PI_V值在0.027至0.087之间。沙丘基因型中SD在环境温度vs. +1 °C和+1 °C vs. +4 °C对比中最高，LDMC在环境温度vs. +4 °C对比中最高；树林基因型中I₁在所有三个温度对比中均最高。种群比较显示，LDMC在环境温度vs. +1 °C（沙丘PI_V = 0.058 vs. 树林0.038）和+1 °C vs. +4 °C（树林0.046 vs. 沙丘0.027）对比中以及SLWC在环境温度vs. +1 °C对比（沙丘0.068 vs. 树林0.033）中存在显著差异。

3.3 性状相关与整合模式

相关分析显示，沙丘基因型在+1 °C时LDMC与SD（r = 0.488）和I₁（r = 0.512）正相关，+4 °C时SLA与LDMC负相关（r = ?0.622），SD与I₁正相关（r = 0.506）。树林基因型在+1 °C时LT与SD正相关（r = 0.527），SLA与LDMC（r = ?0.533）和SD（r = ?0.609）负相关，+4 °C时SLA与LDMC（r = ?0.500）和SD（r = ?0.521）负相关。Mantel检验显示六性状集下相关矩阵显著相似，但四性状敏感性分析中不再显著。机制复合性状组在跨温度比较中保持正相关，而功能—机制复合性状组间的负相关在排除SLWC和LT后不再显著。r_VE分析未检测到温度处理间或种群间表型整合程度的显著差异。

四、讨论与结论

研究表明，矮鸢尾叶片响应涉及功能性状向更资源获取型策略的转变，同时伴随水合—结构、气孔—解剖和叶绿素相关组分的改变。不同性状的响应轨迹各异：SLA和LDMC跨两个温度增量渐进变化，I₁主要在初始增温阶段响应，SD的主要下降发生在后期。两个种群在总体性状水平和部分响应动态上存在差异，沙丘基因型维持较高的SLA和I₁，树林基因型维持较高的LDMC，最明显的种群特异性响应涉及SLWC和LT。性状相关结构在种群和温度处理间存在描述性变异，但整体表型整合程度未检测到显著差异。由于温度处理为顺序施加，时间因素与温度序列共变，+4 °C下的响应应在先前暴露于+1 °C的背景下解释。研究结论强调，预测矮鸢尾对未来增温的响应需联合考虑性状表达、表型可塑性和性状协变，同时注意多变量模式对性状组成选择的依赖性。