喜马拉雅百里香（Thymus linearis）花粉数量–质量权衡与耐胁迫形态的海拔梯度位移

《Discover Plants》：Elevational shift in pollen quantity–quality trade-off and stress-tolerant morphology in Himalayan Thymus linearis

　　

摘要：高海拔环境施加多种非生物胁迫因子，塑造了植物的繁殖和形态响应；然而，同一物种内花粉数量–质量关系沿海拔梯度的检验仍较罕见。研究人员调查了克什米尔喜马拉雅地区五个样点（1800–3499 m）的Thymus linearis的16个营养与繁殖性状。每个样点取15个个体测定营养与花部性状，每样点3个个体测定花粉性状（分别为n = 75和n = 15）。所有测定性状在样点间均差异显著（p ≤ 0.003），除节间长度外其余性状在p < 0.001水平显著。叶长和叶宽表现出显著的样点水平差异，若干叶片维度随海拔升高而增加。花粉–胚珠比（P/O）随海拔显著升高（R² = 0.741, p = 0.0354），而花粉可染色性显著降低（R² = 0.780, p = 0.0227）。二者的直接关系呈强负相关但不显著（r = ?0.841, R² = 0.707, p = 0.0743），与提出的花粉数量–质量关系一致。高海拔种群表现出更大的茎围（+71%）、根围（+76%）、更宽的花（+16%）和更少的花序数（?64%）。主成分分析解释了总变异的88.17%，其中PC1（65.62%）代表营养结构分化，PC2（22.55%）代表植株架构和繁殖形态的额外变异。总体而言，各群落在繁殖性状、叶片形态和植株生长形态上表现出显著的种内分化，高海拔结构性状与耐胁迫预期大体一致。这些发现为理解喜马拉雅海拔梯度上协同性状变异提供了功能基础。

论文解读：喜马拉雅百里香（Thymus linearis）沿海拔梯度的花粉数量–质量权衡与耐胁迫形态协同变异

一、研究背景与科学问题

高海拔环境对植物生长和繁殖施加多种非生物胁迫，包括低温、强紫外辐射、干燥风、短生长季以及不可预测的传粉者 availability。植物通常通过局部适应或表型可塑性在营养和繁殖性状上响应这些胁迫。在群落水平，花粉–胚珠比（P/O）作为繁殖系统保守指标，一般随海拔升高而增加，提示在传粉者稀缺条件下植物增加花粉投资或转向异交；叶片功能性状也随海拔显著变化，高海拔环境常选择更保守的叶片表型；生物量在叶、茎、根之间的分配对环境条件敏感。然而，群落水平模式并不能完全预测单个物种沿同一环境梯度的响应，因为群落性状均值既包含物种更替也包含种内变异。尽管已有全球尺度分析显示某些功能叶性状随海拔一致变化，但叶面积并不存在普遍的种内响应；地中海灌木群落也表现出显著的叶性状种内和种间变异。因此，亟需开展单物种研究以区分物种特异性响应与群落平均模式。

Thymus linearis Benth.（唇形科）是西喜马拉雅特有的多年生半灌木，分布海拔范围广，主要虫媒传粉，兼具有性生殖和营养繁殖。其宽阔的海拔分布和易于测定的繁殖、营养性状使其成为研究喜马拉雅海拔梯度种内功能响应的理想材料。本研究针对三个未解决的科学问题：第一，花粉数量与花粉质量的关系尚不明确。群落研究显示P/O比随海拔升高，但花粉数量单独并不能完全代表繁殖表现，花粉质量（如可染色性）可能也随环境梯度变化。研究人员推测高海拔种群可能表现出更高的P/O比但更低的花粉可染色性，即存在花粉数量–质量权衡。第二，叶片形态的海拔响应存在争议。尽管高海拔常与资源保守型叶片相关，但种内叶片响应并不普遍。同属物种Thymus mongolicus的地理分布表现出显著的叶片维度环境相关变异，提示T. linearis可能并不符合叶片随海拔减小的一般预期。研究人员假设较大叶片可增加光合面积，在受限的高海拔生长季中促进碳获取（光合补偿假说）。第三，植物生长形态沿海拔的种内转变尚不清楚。环境胁迫可改变植株架构和结构投资模式，Grime的CSR框架预测环境约束增强时耐胁迫特征表达增加。研究人员假设T. linearis种群从中间海拔相对较高、分枝较多的生长型，向高海拔更大的结构投资（更粗茎和更发达根系）转变。

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

本研究在克什米尔喜马拉雅山谷选取5个样点（海拔1800–3499 m），涵盖不同气候、植被和土壤条件。每样点随机选取15个健康成熟个体（n = 75）测定营养和花部性状；因花粉计数工作量大，每样点取3个个体（n = 15）测定花粉性状。采样于2024年7–8月盛花期进行。主要技术方法包括：（1）形态性状测定：对株高、一级/二级分枝数、节间长、茎围、根围、上/下叶长宽等10个营养性状，以及花序数、每花序花数、花长、花宽、花粉可染色性、P/O比等6个繁殖性状进行标准化测量；（2）花粉分析：采用1%醋酸洋红染色法评估花粉可染色性，每样本计数至少300粒花粉；P/O比通过每花总花粉数除以胚珠数（T. linearis每花4个胚珠）计算；（3）统计分析：使用单向方差分析（ANOVA）和Tukey HSD事后检验比较样点间差异；用线性混合效应模型检验花粉性状的海拔趋势；基于种群均值的普通最小二乘（OLS）回归分析形态性状与海拔关系；对14个形态性状（排除两个花粉变量）进行主成分分析（PCA），另以全部16个性状进行补充PCA验证结构稳健性；并通过相关矩阵分析性状间协变模式。

三、主要研究结果

（一）营养性状沿海拔的变异

所有10个营养性状在5个种群间差异显著。最高海拔样点（Pir ki Gali，3499 m）茎围最粗（1.85 cm）、根围最大（2.17 cm），而中海拔样点（Drung 2146 m，Dara 2276 m）植株更高（达16.73 cm）、二级分枝更多（达33.2）。与叶片随海拔减小的一般预期相反，所有4个叶片维度在种群间均差异显著（ANOVA F值40.4–800.0，p < 0.001），其中下叶宽表现出最强的种群间分化（F_4,70 = 800.0），从1800 m的0.18 cm增至3499 m的0.46 cm，差异达156%。线性回归显示茎围每100 m增加0.024 cm（R² = 0.950，p = 0.004，总增幅71%），根围每100 m增加0.029 cm（R² = 0.836，p = 0.031，总增幅76%），下叶宽每100 m增加0.009 cm（R² = 0.950，p = 0.004）。二级分枝数随海拔显著减少（每100 m减少0.238个，R² = 0.803，p = 0.040）。

（二）繁殖性状沿海拔的变异

繁殖性状表现出分歧的海拔模式。P/O比随海拔显著升高（R² = 0.741，p = 0.035），从1800 m的3268.9增至3499 m的4740.4，增幅45%。相反，花粉可染色性随海拔显著降低（R² = 0.780，p = 0.023），从91.3%降至70.8%，降幅22%。P/O比与花粉可染色性呈强负相关（r = ?0.841，R² = 0.707，p = 0.074），虽未达到α = 0.05显著水平，但方向与花粉数量–质量权衡一致。花宽每100 m增加0.002 cm（R² = 0.870，p = 0.021），最高海拔种群花宽增加16%。花序数随海拔呈近显著下降趋势（R² = 0.714，p = 0.071），最高海拔种群花序数比Kangan少64%。

（三）性状整合与协变模式

相关分析显示茎投资与根投资紧密耦合（r = 0.93），且均与叶片性状正相关：根围与上叶长（r = 0.97）、下叶宽（r = 0.99）高度正相关，表明高海拔存在粗茎、大根、大叶的协同耐胁迫综合征。P/O比与根围显著正相关（r = 0.94，p < 0.05），与茎围正相关但不显著（r = 0.88）；花粉可染色性与茎围显著负相关（r = ?0.95，p < 0.05），与根围和叶性状的负相关不显著。花序数与根围显著负相关（r = ?0.97，p < 0.01）。这些模式表明繁殖性状与结构性状沿海拔梯度存在协调的种群水平变异。

（四）多变量性状空间中的种群分化

基于14个形态性状的PCA显示前两个主成分解释总变异的88.17%（PC1 = 65.62%，PC2 = 22.55%）。PC1主要代表营养结构投资梯度，正向载荷最大的是根围（0.329）、下叶宽（0.329）、上叶长（0.323）、一级分枝（0.304）和茎围（0.295），负向载荷包括花序数（?0.324）、花长（?0.302）、节间长（?0.277）和二级分枝（?0.268）。PC2代表植株架构和繁殖形态的次级梯度。PCA双标图清晰分离5个种群：Pir ki Gali位于PC1正极端，与粗茎、大根围、宽叶相关，表现出最强的耐胁迫表型；Gulmarg位于PC1正区域但沿PC2与Pir ki Gali区分；Kangan和Dara位于PC1负侧；Drung以强烈负PC2得分为特征。包含全部16个性状的补充PCA结果总体结构一致（前两轴解释87.87%），且P/O比与花粉可染色性在PC1上载荷方向相反（P/O = 0.293，PV = ?0.248），进一步支持数量–质量权衡。

四、讨论与结论

本研究发现三个值得注意的模式：第一，P/O比与花粉可染色性在种群水平呈强负相关（虽然不显著），且二者沿海拔呈现显著相反趋势，这与花粉数量–质量权衡假说一致。高海拔种群P/O比显著升高，提示在传粉者访问可能稀少且不可预测的条件下，植物增加雄配子生产以提高异交成功率；而花粉可染色性降低可能与高海拔UV-B辐射损伤、低温干扰花粉发育或贫瘠土壤资源限制有关。第二，叶片维度随海拔增加而增大，与常见叶面积随海拔减小的模式相反，支持短生长季下光合补偿假说——较大叶片可在有限时间内最大化碳获取，超过蒸腾和热负荷成本。第三，植物生长形态从中间海拔的相对高大、多分枝型，向高海拔的粗茎、大根、少花序型转变，与CSR框架下耐胁迫预期大体一致。然而，所测性状不允许正式分类CSR策略，仅作为解释框架。与群落水平研究比较，本研究的P/O结果与西南中国84物种的海拔增加模式定性一致，但补充了花粉可染色性数据，提示群落水平P/O升高可能部分反映数量–质量权衡而非单纯异交增加。研究局限性包括花粉性状采样强度较低（每种群3个体）、醋酸洋红染色反映可染色性而非直接萌发能力、海拔观测设计无法区分环境驱动因素及表型可塑性 vs. 遗传分化。未来需通过同园实验、移栽实验、增加花粉采样、直接测定传粉者访问和繁殖成功率以及花粉萌发实验来验证机制。在气候变化背景下，高海拔种群已表现降低的花粉可染色性（70.8%），若温度持续上升，花粉质量可能进一步下降，威胁高海拔种群繁殖成功；而低海拔种群可能面临上行迁移物种的竞争。保护应注重维持海拔连通性和多个种群，以保留物种的功能多样性。

研究结论：本研究证明了T. linearis沿喜马拉雅海拔梯度的营养和繁殖性状存在协调的种群水平变异。高海拔种群具有更大的茎围和根围、若干更大的叶片维度，同时花序数和花粉可染色性较低。这些形态模式与CSR框架下的耐胁迫预期大体一致，但所测性状不允许正式CSR分类。P/O比随海拔升高而花粉可染色性降低，二者直接负相关虽强但不显著，因此与花粉数量–质量权衡保持一致而非确证。结果显示了显著的种内分化，但未区分表型可塑性与遗传基础的局部适应。论文发表于《Discover Plants》。