## 中文标题 揭示欧洲黑松（*Pinus Nigra* J.F. Arnold）幼苗干旱响应的种内变异

《Plant, Cell & Environment》：Unravelling the Intraspecific Variation in Drought Responses in Seedlings of European Black Pine (Pinus Nigra J.F. Arnold)

　　

理解干旱耐受性的种内变异对于预测气候变化下森林树种的适应能力至关重要。然而，在生态和经济上具有重要意义的针叶树种中，这种变异的分子基础仍鲜为人知。研究人员将高通量表型分析与代谢组学和转录组学相结合，在标准化土壤干旱条件下研究了针叶树种*Pinus nigra*九个气候迥异种源的干旱响应。研究人员检验了干旱耐受性（以光系统II最大量子产额（Fv/Fm）的下降来衡量）是否在种源间存在差异、是否遵循气候渐变群（clinal cline）、以及是否涉及与生长的权衡。为揭示其潜在分子基础，研究人员对代表对比性干旱耐受性的四个种源进行了代谢组学和转录组学分析。干旱耐受性在种源间存在显著差异，且与生长解偶联，但未沿气候渐变群表现出分化。耐旱种源与敏感种源在类黄酮化合物松属素（pinostrobin）和二萜类化合物羟基脱氢枞酸（hydroxydehydroabietic acid）的组成型丰度及干旱诱导积累方面均存在差异。转录组图谱进一步突显了与类黄酮相关基因表达中种源特异性的差异。总体而言，这些发现将表型干旱响应与种源间的分子变异联系起来，为理解这一重要针叶树种干旱响应的分子基础提供了见解。

## 论文解读

### 一、研究背景与科学问题

干旱是限制植物生存和生产力的重要环境胁迫因子。随着气候变化加剧，干旱事件日益频繁、持久和剧烈，对全球生态系统稳定性构成重大威胁。在针叶树种中，极端干旱事件与森林死亡率上升和生产力下降密切相关。理解物种内（种内）干旱耐受性的变异对于预测森林树种的适应潜力、制定辅助迁移、种源选择和保育规划等策略至关重要。欧洲黑松（*Pinus nigra* J.F. Arnold）是一种在欧洲广泛分布但遗传碎片化的针叶树种，具有重要的生态和经济价值，其广泛的地理和气候分布使其成为研究针叶树种干旱耐受性种内变异的理想模型。然而，此前基于生长、存活和生理性状的研究结果并不一致——有的研究报告种源间存在显著差异，有的则未发现差异。更重要的是，关于种源水平干旱耐受性差异背后的分子响应和机制知之甚少，尤其是在选择压力可能最强的幼苗阶段。因此，亟需在严格控制的环境条件下，结合多层次的分子分析，系统评估欧洲黑松种源间的干旱耐受性变异及其分子基础。

### 二、主要技术方法

研究人员利用配备机器人称重灌溉系统的高通量植物表型平台（HTPP），在标准化干旱条件下对来自九个遗传多样种源的欧洲黑松幼苗进行表型分析，通过暗适应叶绿素荧光（Fv/Fm）评估干旱耐受性。随后，选取四个代表对比性干旱耐受性的种源（两个耐旱、两个敏感），进行靶向代谢组学（定量36种已知干旱相关代谢物）和非靶向代谢组学（UHPLC-MS）分析，并利用mRNA-seq（QuantSeq 3' mRNA-Seq文库制备，Illumina NextSeq 2000平台测序，映射至*P. tabuliformis*参考基因组）进行转录组分析。样品来自奥地利维也纳生物中心植物科学设施（VBCF Vienna BioCenter Core Facilities GmbH）培育的幼苗。

### 三、研究结果

**3.1 干旱对光系统II效率的影响**：干旱处理下，各品种幼苗Fv/Fm在第3天前保持接近最大值（~0.83），随后下降并达到品种特异性的最低值。包含品种作为固定效应的GAMM模型显著优于不含品种的模型（p < 0.0001），表明品种间Fv/Fm下降存在显著差异。将冠层面积作为协变量并未改善模型拟合，排除植物大小对干旱响应的干扰。基于PLCF（叶绿素荧光损失百分比）估算的干旱耐受性在品种间存在显著差异（p = 0.001），PN1（奥地利）和PN3（西班牙）表现出最低PLCF值（耐旱），而PN6（克罗地亚）PLCF最高（干旱敏感）。

**3.2 干旱对生长的影响**：干旱处理显著降低了所有品种的生长，平均冠层面积增量（CAI）在干旱条件下约为4%，而在充分浇水条件下约为60%（P < 2.2e?16）。品种间CAI差异显著（P < 4.4e?07），PN1在两种处理下均表现出最高CAI。处理与品种的交互作用不显著（p = 0.075），但品种对相对冠层面积增量（rCAI）影响显著。

**3.3 生长与干旱耐受性之间缺乏权衡**：PLCF与CAI在充分浇水（R = ?0.73，p = 0.031）和干旱（R = ?0.81，p = 0.011）条件下均呈显著负相关，表明较高生长与较高干旱耐受性相关，不存在权衡关系。rCAI与充分浇水条件下的CAI关系不显著，表明生长较快的品种不一定更易受干旱影响。

**3.4 黑松干旱耐受性沿气候梯度的渐变群变异**：PLCF与19个生物气候变量之间均无显著相关性，但PLCF与最湿月降水量（BIO13）和年均温（BIO1）呈弱正趋势（Spearman ρ ≈ 0.3），表明来自更湿润、更温暖地区的品种倾向于表现出略高的PLCF。生长相关变量与气候的相关性同样较弱且不显著，未发现强渐变群变异。

**3.5 靶向代谢组学揭示物种广泛和品种特异的干旱响应**：在36种靶向代谢物中，13种在至少一个品种中干旱处理后显著增加（FDR < 0.05），表明存在物种广泛的干旱响应。品种间这些化合物水平差异显著，但无一种在耐旱品种中持续高于敏感品种。干旱诱导方面，脱落酸（ABA）、脯氨酸和玉米黄质在敏感品种PN6中表现出显著更强的诱导倍数变化，而类黄酮没食子儿茶素和原花青素B1在耐旱品种PN1中表现出显著更强的诱导。萜烯类（α-蒎烯、β-蒎烯和总萜烯）虽无显著处理效应，但耐旱品种在干旱条件下总体水平趋于更高。

**3.6 非靶向代谢组学鉴定与增强干旱耐受性相关的候选代谢物**：非靶向分析初始检测到3205个代谢特征。PCA显示各品种在PC1（24.3%方差）和PC2（20.7%方差）上分别聚类，表明品种间存在独特的基线代谢差异。聚焦PC3（13.7%方差）后，鉴定出84个精选特征。层次聚类显示，独立于处理，耐旱组（PN1、PN3）与敏感组（PN4、PN6）沿PC2（24.9%）明显分离，主要由Cluster 3代谢物驱动，包括奎宁酸（酚类前体）、松属素（类黄酮）、羟基脱氢枞酸（二萜类）、阿朴类胡萝卜素和侧金盏花黄质（四萜类）等。这些代谢物在耐旱品种中组成型丰度更高且在干旱下诱导更强。羟基脱氢枞酸表现出显著的品种×处理交互作用（p < 0.0001），耐旱品种基态水平更高（log2FC_DT/DS > 1），且干旱下积累显著更多（log2FC_D/WW > 1），而敏感品种几乎无变化（log2FC_D/WW = ?0.13）。松属素在耐旱品种PN1中诱导显著更强。

**3.7 转录组分析揭示物种广泛和品种特异的干旱响应**：PCA沿PC1（74.5%方差）分离处理，沿PC2（7.5%）在干旱下分离耐旱与敏感组。耐旱品种的差异表达基因（DEG）数量约为敏感品种的一半，且DEG数量与Fv/Fm呈显著负相关（R2 ≈ 1，p < 0.05）。所有品种共享906个干旱响应基因，富集于渗透胁迫、光保护、原花青素生物合成、光合作用和水分运输等GO条目，包括脯氨酸生物合成（P5CS）、糖转运（SWEET）、ABA信号（SNRK2、PYL、PP2C）和LEA蛋白等基因。耐旱品种共享的DEG中包含组蛋白修饰（H3）、钙信号（CIPK3）、热激响应（HSP20）和类黄酮生物合成（查尔酮合酶CHS、花青素还原酶ANR）相关基因。敏感品种共享的DEG中，上调基因富集于糖转运、苯丙烷代谢和氧化应激响应，下调基因与光合机制相关（psbS、psaK、petF、ATP合酶），而这些基因在耐旱品种中相对稳定。

**3.8 DEG层次聚类识别共调控基因模块**：对所有5137个DEG进行层次聚类，鉴定出三个主要簇。Cluster 1（1671个基因）和Cluster 2（3025个基因）分别包含干旱下上调和下调的基因，两者在敏感品种中变化更强。Cluster 3（441个基因）包含干旱诱导且在耐旱品种PN3中表达显著更强的基因，富集于类黄酮生物合成、细胞分裂素信号和应激响应。Cluster 3包含编码CHS、F3H、DFR、ANS、ANR和COMT-like的类黄酮通路基因拷贝，这些基因在耐旱品种PN3中诱导更强，其中CHS（Pt3G29450）和ANR（PtJG50190）为耐旱品种共享DEG。

**3.9 干旱下黑松萜烯生物合成的转录调控**：未发现富集于萜烯生物合成的表达簇。虽然脱落二烯合酶（abietadiene synthases）非DEG，但上游MEP通路基因在敏感品种中表现出相对更强的抑制（平均log2FC: ?1.2），而耐旱品种抑制较弱（平均log2FC: ?0.41）。

### 四、讨论与结论

研究发现欧洲黑松在幼苗阶段即表现出明显的干旱响应种内变异，体现在光合效率、代谢和转录多个生物学层面。耐旱品种（PN1、PN3）在干旱下保持较低PSII损伤（低PLCF），且转录组相对稳定，DEG数量约为敏感品种的一半，表明其具有更强的维持生理功能的能力。敏感品种（PN4、PN6）则表现出更强的光合机制基因下调和广泛的转录重编程。干旱耐受性与生长之间未发现权衡，且与起源地气候无显著关联，这提示气候来源不能作为预测干旱耐受性的唯一指标，也表明耐受性和生长可同时进行遗传改良。非靶向代谢组学和转录组学综合分析将松属素和羟基脱氢枞酸鉴定为潜在的干旱耐受性代谢标记物：松属素作为一种甲基化黄烷酮，其较高的干旱诱导积累可能通过清除脂溶性活性氧（ROS）保护光系统II稳定性；羟基脱氢枞酸作为一种枞酸型二萜树脂酸，其组成型和干旱诱导的高积累可能通过ROS清除作用限制氧化损伤。这些候选代谢物及其共表达的类黄酮生物合成基因（CHS、COMT等）为后续功能验证和分子标记开发提供了基础。论文发表在《Plant, Cell》。

研究结论：该研究提供了欧洲黑松干旱响应种内变异的明确证据，这种变异在幼苗阶段即可检测，并表现在功能、代谢和转录水平。耐旱品种表现出较少的生长和光化学效率下调及较小的代谢和转录重编程。综合转录组和代谢组分析揭示了共享的核心响应和种源特异的分子特征，将松属素和羟基脱氢枞酸鉴定为潜在的代谢标记物，为未来研究其在欧洲黑松及其他针叶树种干旱耐受性中的作用奠定了基础。