本研究由深圳大学生命科学与海洋学院、深海科学与工程研究所（中国科学院）等机构的陈静上、周乔、Jackson Nkoh Nkoh、刘静、王俊杰、胡章立、Quaid Hussain等作者共同完成。论文发表于Elsevier旗下的期刊 ecotoxicology and environmental safety 第281卷（2024），文章编号116609，在线发表日期为2024年6月20日。
学术背景
本研究属于环境生态毒理学与植物生理学的交叉领域，具体聚焦于红树林植物对重金属胁迫的响应机制。随着工业、农业和城市活动的增加，重金属污染日益严重，对包括红树林在内的生态系统构成威胁。铜（Cu）是植物必需的微量营养元素，参与多种生理生化过程，但过量铜会产生毒害作用，影响植物生长、光合作用，并诱导氧化应激。
木榄（*Bruguiera gymnorhiza*）是一种典型的红树植物，在海岸带生态系统中扮演着重要角色，具有净化污染物的潜力。然而，关于木榄如何响应长期铜胁迫，尤其是在气孔反应、叶片蜡质合成及其与分子、生理生化水平的协同作用方面，研究尚不充分。叶片表皮蜡质和气孔行为是植物应对外界胁迫（如减少水分流失）的重要策略。因此，本研究旨在通过整合生理、生化和转录组学分析，探究长期铜胁迫下木榄叶片的适应机制，特别是对气孔大小、蜡质晶体和成分的影响，以期为理解红树林对重金属胁迫的响应提供新见解，并为开发新型管理实践提供参考。
详细研究流程
本研究是一项为期两年的户外控制实验，流程系统且全面，主要包括以下几个关键环节：
植物材料与胁迫处理：实验于2017年11月至2019年12月在广西北海的红树林保护区进行。实验使用一年生木榄幼苗，栽培于装有特定比例沙土混合物的盆栽中。设置三种铜处理浓度：0 mg L⁻¹（对照组，Cu0）、200 mg L⁻¹（Cu200）和 400 mg L⁻¹（Cu400），以CuCl₂溶液形式通过预埋管道每半年施加一次。每个处理设3个重复小区，每个小区10-12株植物。实验为期两年，以确保长期胁迫效应。实验结束后，采集植株上部第二或第三片叶片用于各项分析。
生理与形态指标测定：
- 生长与元素含量：扫描叶片并计算叶面积、长度、宽度和长宽比。将叶片干燥研磨后，采用半微量凯氏定氮法和比色法测定总氮（N）和总磷（P）含量；使用原子吸收光谱仪测定叶片铜含量。
- 光合色素与相关指数：使用80%丙酮提取叶片光合色素，通过分光光度法测定叶绿素a（Chl a）、叶绿素b（Chl b）、总叶绿素和类胡萝卜素（Car）含量。同时，使用植物氮平衡指数仪测定叶片相对叶绿素、花青素、类黄酮含量及氮平衡指数（NBI）。
- 生化指标分析：使用商用试剂盒，通过分光光度法测定抗氧化酶活性，包括超氧化物歧化酶（SOD）、过氧化物酶（POD）、过氧化氢酶（CAT）和抗坏血酸过氧化物酶（APX）。同时测定活性氧（过氧化氢H₂O₂、丙二醛MDA）含量以及渗透调节物质（脯氨酸、可溶性糖）含量。
叶片表面蜡质结构分析：
- 扫描电子显微镜（SEM）观察：将叶片样品干燥后镀金，使用Hitachi S-4800场发射扫描电镜在5 kV加速电压下观察叶片近轴面（上表皮）和远轴面（下表皮）的蜡质晶体形态和气孔结构。
- 蜡质成分分析：结合气相色谱-质谱联用技术（GC-MS），分析叶片蜡质负载量及其化学成分，包括烷烃、醇类、脂肪酸、伯醇、正构烷烃等的含量变化。
转录组学分析：
- RNA提取与测序：从三个处理组（Cu0, Cu200, Cu400）的叶片中各取3个生物学重复，共9个样本，进行总RNA提取。使用Illumina平台进行RNA测序（RNA-seq）。对测序数据进行质量评估，确保Q20和Q30比例合格。
- 基因注释与差异表达分析：使用DESeq2软件包（参数：Fold Change ≥ 1.5，p-value ≤ 0.05）进行差异表达基因（DEGs）的筛选。将DEGs比对到NR、NT、Pfam、KOG、Swiss-Prot、GO和KEGG等数据库进行功能注释。
- 功能富集分析：利用GO和KEGG数据库对DEGs进行功能富集分析，识别显著富集的生物学过程和代谢通路。
- 实时定量PCR验证：从DEGs中选取13个关键基因（涉及光合作用、角质/木栓质/蜡质生物合成、过氧化物酶体、植物激素信号转导等通路），进行实时定量PCR（qRT-PCR）验证，以确保RNA-seq数据的可靠性。
数据分析：所有生理生化数据使用Statistix 8.1软件进行单因素方差分析（One-way ANOVA），并采用最小显著性差异法（LSD）进行多重比较（p < 0.05）。使用GraphPad Prism 9软件进行图表绘制。
主要研究结果
铜胁迫对植物生长、营养吸收和叶片形态的影响：与对照组（Cu0）相比，Cu200和Cu400处理显著抑制了木榄的生长。叶片铜含量随处理浓度升高而增加，而叶片氮和磷含量则显著降低，表明铜胁迫干扰了植物对必需营养元素的吸收。叶片形态发生显著变化：Cu200处理下叶面积减小，而Cu400下叶面积反而有所增加但叶片颜色变深（呈紫/深紫色），叶长宽比也发生改变，表明不同浓度铜胁迫对叶片发育产生了差异化影响。
铜胁迫对光合色素和生理生化指标的影响：
- 光合色素：Cu200处理下，叶绿素a、b、总叶绿素和类胡萝卜素含量均显著高于对照；但在Cu400处理下，这些色素含量均显著下降，表明高浓度铜严重损害了光合系统的功能。花青素和类黄酮含量在Cu200下积累最多，而氮平衡指数在铜胁迫下均下降。
- 抗氧化系统和渗透调节：铜胁迫显著提高了POD和APX的活性，但降低了CAT的活性；SOD活性未发生显著变化。活性氧方面，H₂O₂含量在胁迫下显著升高（Cu400 > Cu200），而MDA含量却显著降低，这可能与抗氧化酶（特别是POD和APX）活性的上调有效缓解了膜脂过氧化损伤有关。此外，脯氨酸和可溶性糖含量在铜胁迫下（尤其是Cu200）显著积累，发挥了重要的渗透调节和抗氧化保护作用。
铜胁迫对叶片蜡质结构和成分的影响：
- SEM观察：随着铜浓度增加，叶片近轴面和远轴面的蜡质晶体沉积增多。在远轴面（气孔主要分布面），高浓度铜（Cu400）导致蜡质晶体显著覆盖气孔，迫使气孔关闭，这直观地解释了铜胁迫下可能的气孔限制机制。
- GC-MS分析：与对照相比，Cu200和Cu400处理下，叶片蜡质负载总量、烷烃、醇类、伯醇（C26:0, C28:0, C30:0, C32:0）以及正构烷烃（C29, C30）的含量显著更高，而脂肪酸（C12, C16, C18）含量降低。这些成分变化表明，木榄通过重塑叶片蜡质组成和增加蜡质沉积来应对铜胁迫。
转录组学分析结果：
- 差异表达基因：比较分析共鉴定出大量DEGs。其中，Cu0 vs Cu200有1240个（771上调，469下调），Cu0 vs Cu400有1000个（723上调，277下调），Cu200 vs Cu400有1476个（808上调，668下调）。Cu200与Cu400之间的DEGs数量最多，说明不同胁迫强度引发了差异显著的分子响应。
- 通路富集分析：KEGG通路富集分析显示，铜胁迫主要影响了八条代谢通路，包括：光合作用、苯丙烷生物合成、植物激素信号转导、角质/木栓质/蜡质生物合成、淀粉和蔗糖代谢、脂肪酸延伸、类黄酮生物合成以及内质网中的蛋白质加工。其中，角质/木栓质/蜡质生物合成通路在Cu200 vs Cu400组中显著富集，与上述蜡质观察结果直接对应。
- 关键基因表达：qRT-PCR验证了RNA-seq数据的可靠性。分析发现：
- 与光合作用相关的基因（如 psbc, atpi, psbs, atpa 等）在Cu0 vs Cu200和Cu0 vs Cu400中多为上调，这与Cu200下光合色素增加的现象部分吻合，但高浓度铜（Cu400）下光合功能的最终受损可能是多因素综合结果。
- 与抗氧化相关的多个*POD*（过氧化物酶）基因（如 per3, per10, *per12*等）显著上调，与POD酶活性升高结果一致。
- 在角质/木栓质/蜡质生物合成通路中，HTH (hothead)、*FAR2*、FAR3 (脂肪酸酰基-CoA还原酶) 和 CER1 (eceriferum) 等基因在Cu200 vs Cu400中上调，而 HHT1 基因下调。这些基因的表达变化为蜡质成分和沉积量的改变提供了分子解释。
- 植物激素信号转导通路中，涉及生长素（Auxin）、细胞分裂素（Cytokinin）、赤霉素（Gibberellin）、脱落酸（ABA）、乙烯（Ethylene）等激素的多个受体、响应因子及转录因子基因（如 ARF, AUX/IAA, SAUR, TIFY 等）发生差异表达，表明植物激素网络深度参与了铜胁迫的调控。
结论
本研究通过整合生理、生化和转录组学分析，系统揭示了木榄（*Bruguiera gymnorhiza*）响应长期铜胁迫的复杂机制。主要结论如下： 1. 长期铜胁迫抑制木榄生长，减少氮磷吸收，损害光合作用（尤其在400 mg L⁻¹浓度下）。 2. 木榄通过激活特定的抗氧化酶系统（POD, APX）、积累渗透调节物质（脯氨酸、可溶性糖）以及降低MDA含量，有效缓解了铜诱导的氧化损伤。 3. 铜胁迫诱导木榄叶片蜡质晶体沉积显著增加，蜡质化学成分发生重塑（烷烃、伯醇增加，脂肪酸减少），并且高浓度铜导致气孔被蜡质覆盖而关闭。这是木榄减少水分流失和可能限制重金属进一步进入的物理适应策略。 4. 转录组学分析从分子层面证实了上述生理生化表现。铜胁迫差异性地调控了涉及光合作用、抗氧化（苯丙烷生物合成）、蜡质合成、激素信号和初级代谢（淀粉蔗糖代谢）等一系列关键通路的基因表达。特别是角质/木栓质/蜡质生物合成通路相关基因的上调，直接关联于叶片蜡质沉积的增加和气孔行为的改变。
研究价值与意义
本研究的科学价值在于首次从气孔形态、蜡质物理结构、化学成分与分子调控网络联动的角度，全面阐释了红树植物木榄应对长期铜胁迫的适应性策略。它不仅丰富了植物重金属胁迫生理学的内容，也为红树林生态系统中重金属污染的生态毒理学评估提供了新的生物标志物（如特定蜡质成分、相关基因表达谱）和理论依据。应用价值方面，研究结果有助于评估红树林在重金属污染海岸带的生存潜力与修复功能，为利用红树林进行生态修复时的物种选择和管理（如关注其蜡质屏障特性）提供科学指导。
研究亮点
- 研究视角新颖：将气孔大小、蜡质晶体结构（SEM观察）、蜡质化学成分（GC-MS分析）与高通量转录组学（RNA-seq）深度整合，从微观形态到分子机制全方位解析了木榄对铜胁迫的响应，突破了以往研究多局限于生理或分子单一层面的局限。
- 长期胁迫实验：进行了为期两年的户外长期铜胁迫处理，更贴近自然条件下红树林面临的慢性污染情景，所获数据更能反映植物的真实适应过程。
- 机制阐释深入：不仅观察到了蜡质沉积增加和气孔关闭的现象，还通过转录组学找到了驱动这一现象的潜在关键基因（如 FAR2, FAR3, *CER1*），并将激素信号转导等调控网络与植物的整体生长响应相联系，构建了相对完整的“胁迫信号感知-分子调控-生理生化适应-形态结构改变”响应链条。
- 发现矛盾与深入探讨：研究注意到了高浓度铜下MDA含量反而降低这一看似矛盾的现象，并将其与抗氧化酶（POD, APX）活性的显著上调以及相关基因的表达关联起来，进行了合理解释，体现了分析的深度。