《Plants》：Transcriptome-Based Dissection of the Molecular Mechanisms Underlying Flooding Stress Responses of Eastern Cottonwood in the Floodplains of the Middle and Lower Reaches of the Yangtze River
Guowei Huang,
Xueli Zhang,
Xinye Zhang,
Ning Liu,
Changjun Ding,
Jinhua Li,
Fenfen Liu,
Kailian Long,
Chengcheng Gao and
Qinjun Huang
+ 2 authors
编辑推荐：
洪水胁迫严重影响长江中下游滩地杨树生长与生存。研究人员通过整合多年生长调查与转录组测序，解析了两个不同耐涝性的东部棉杨（Populus deltoides）品种对季节性淹水的表型与分子响应，揭示耐涝品种HBI通过激活乙烯（ETH）信号、协调光合恢复、类黄酮（Flavonoid）生物合成与抗氧化防御等途径维持生长优势的适应机制，为耐涝林木品种选育及林业可持续发展提供理论依据。
论文解读
在全球气候变化加剧的背景下，极端降雨和洪水事件频发，已成为制约植物生长和陆地生态系统健康的重大威胁。其中，洪水引发的缺氧（Hypoxia）胁迫是限制作物和树木生长的最关键非生物胁迫之一，据统计可导致全球作物产量年损失高达15%至80%。长江中下游地区拥有广阔的季节性淹没滩地，夏季年淹没期可达0-60天，这严重限制了适宜树种的种植选择。杨树是中国重要的速生用材树种，广泛种植于该地区，但其在长期季节性淹水下的生存能力和生长表现存在显著差异。为了探究树木应对淹水胁迫的内在分子机制，并筛选出优良的耐涝品种，研究人员对两个具有不同表现的关键东部棉杨（Populus deltoides）品种进行了深入研究。
研究人员在长江中游的李家洲林场建立了长期定位观测试验，比较了耐涝的‘江汉1号’（HBI）与对淹水敏感的‘楚林2号’（CL）这两个东部棉杨品种在长达七年（2016-2022年）季节性淹水胁迫下的生长表现与存活率。通过整合多年的生长调查数据与叶片转录组测序（RNA-seq）分析，他们系统地解析了两个品种在淹水胁迫下的表型差异和分子响应网络。相关成果发表在期刊《Plants》上，为林木耐涝性研究及品种选育提供了宝贵的分子见解。
主要技术方法
研究团队在长江中下游滩地典型季节性淹没区建立了包含HBI和CL在内的杨树品种试验林，对株高和胸径进行了长期监测。在2024年洪水季（淹水前5月29日和水退后8月20日），研究人员采集了叶片样品，随后进行了RNA提取、文库构建和高通量测序。通过对测序数据进行质量控制和比对，利用差异表达基因（DEGs）分析、基因本体（GO）富集和京都基因与基因组百科全书（KEGG）通路分析，系统比较了两个品种的转录组响应。最后，通过定量实时聚合酶链式反应（qRT-PCR）对关键差异表达基因的表达模式进行了验证。
研究结果
2.1. 淹水对东部棉杨生长和存活的影响
经过7年季节性淹水，耐涝品种HBI的存活率为73.91%，其树高（23.1米）和胸径（26.3厘米）均显著优于CL品种，证实HBI是一个耐涝性强的品种。
2.2. RNA测序数据质量评估
对样本的转录组测序数据质量进行评估，结果表明，所有样本的清洁读段比例均超过98.20%，整体比对率高达96.45%至98.05%，Q20和Q30值均很高，证实测序数据质量优异，可用于后续分析。
2.3. 淹水前后HBI和CL的差异表达基因
淹水胁迫后，在HBI叶片中鉴定到3195个差异表达基因（DEGs），其中1030个上调，2165个下调；在CL叶片中鉴定到2021个DEGs，其中516个上调，1505个下调。两品种共同存在1098个DEGs，其中971个下调，127个上调。qRT-PCR验证结果与RNA-seq数据高度一致，证实了转录组分析的可靠性。
2.4. HBI和CL共享差异表达基因的功能注释
对两品种共享的DEGs进行GO和KEGG富集分析发现，这些基因主要富集于信号转导、胁迫响应、谷胱甘肽代谢过程等生物学过程，以及脂肪酸降解、昼夜节律、次级代谢物生物合成等通路。这表明两品种可能通过共同的应激感知、防御反应、代谢重编程和生物钟协调机制来应对淹水胁迫。
2.5. HBI和CL特异性差异表达基因的功能注释
品种特异性DEGs的富集分析揭示了关键差异。在耐涝品种HBI中，其特异性DEGs显著富集于植物激素信号转导、MAPK（促分裂原活化蛋白激酶）信号通路、淀粉和蔗糖代谢、光合作用、类黄酮生物合成等通路。而在不耐涝的CL中，其特异性DEGs则主要富集于过氧化氢分解过程、细胞对DNA损伤刺激的反应、对活性氧（ROS）的反应等氧化应激和损伤修复通路。
2.6. 参与东部棉杨淹水响应的类黄酮生物合成通路
在耐涝品种HBI中，鉴定出12个在类黄酮生物合成通路中显著富集的DEGs。这些基因编码花青素合成酶（ANS）、无色花青素还原酶（LAR）、黄烷酮3-羟化酶（F3H）、查尔酮异构酶（CHI）、咖啡酰辅酶A O-甲基转移酶（CCoAOMT）、羟基肉桂酰辅酶A奎宁酸/莽草酸羟基肉桂酰转移酶（HCT）和反式肉桂酸4-单加氧酶（C4H/CYP73A）等关键酶。这些酶的差异表达表明，HBI在淹水条件下对类黄酮积累和花青素生物合成进行了特定的调控调整。
2.7. 淹水胁迫下HBI和CL的不同响应模式
基于差异表达分析结果，研究人员提出了一个耐涝调控模型。耐涝品种HBI在淹水胁迫下，通过特异性激活乙烯（ETH）及其他激素相关基因，构建了一个协同的调控网络。这个网络能够同时协调缺氧适应、抵御再充氧诱导的氧化应激、恢复光合功能以及类黄酮介导的抗氧化保护，从而维持生长活力，提高在季节性淹没环境中的存活率。相比之下，淹水敏感品种CL则缺乏这种ETH-类黄酮调控网络，主要通过激活针对ROS积累的响应通路来应对胁迫。
结论与意义
本研究通过长期生长调查与转录组学分析相结合，揭示了两个东部棉杨品种应对季节性淹水胁迫的不同分子机制。耐涝品种HBI表现出优越的生长和存活能力，这与其转录组层面激活的特定调控网络密切相关。该网络的核心特征包括：显著激活了乙烯（ETH）等植物激素信号转导通路，从而协调了整体的胁迫响应；通过调节光合作用相关基因表达，应对淹水造成的光合抑制；并特异性激活了类黄酮生物合成通路，利用类黄酮强大的抗氧化能力清除活性氧（ROS），减轻氧化损伤。而敏感品种CL的响应则更多地聚焦于氧化应激的直接应对和损伤修复，这可能是其在淹水后遭受严重氧化胁迫的表现。
此项研究不仅在转录组水平明确了杨树品种间耐涝性差异的分子基础，更重要的是，它为林木耐涝品种的分子标记辅助育种提供了关键靶点，例如那些参与乙烯信号和类黄酮合成的关键基因。研究结果不仅拓宽了木本植物淹水胁迫生物学的研究前沿，也为长江中下游等易受洪水影响地区的林业可持续发展提供了重要的理论基础和技术支撑。未来，整合代谢组学、蛋白质组学分析，以及对根部等关键胁迫感知器官的研究，将有助于更全面地阐明杨树耐涝性的分子调控网络，并推动其遗传改良实践。
