整合生理与转录组分析揭示辣椒（Capsicum annuum L.）耐涝性的多层调控网络及WRKY/HB/C2H2候选调控因子

《Scientia Horticulturae》：Integrative physiological and transcriptomic analysis reveals multi-layered regulatory mechanisms underlying waterlogging tolerance in pepper ( Capsicum annuum L.)

整合生理与转录组分析揭示辣椒（Capsicum annuum L.）耐涝性的多层调控机制；整合生理与转录组分析揭示辣椒在水涝胁迫下的多层调控网络；耐涝辣椒基因型协调抗氧化系统与渗透调节以减轻氧化损伤；耐涝基因型中不饱和脂肪酸、玉米素（zeatin）、支链氨基酸（branched-chain amino acid, BCAA）降解和ABC转运蛋白（ATP-binding cassette transporters, ABC transporters）代谢途径被激活；WRKY、HB和C2H2转录因子被鉴定为连接基因表达与生理适应的候选调控因子；提出整合抗氧化防御、渗透调节和代谢重编程的辣椒耐涝工作模型。

研究背景：辣椒（Capsicum annuum L.）原产墨西哥，是全球第二大茄科蔬菜，中国种植面积和经济产出较大。辣椒根系浅、吸收根少，耐涝性较弱；贵州降雨集中且喀斯特地貌排水差，水涝常导致代谢紊乱、生长受抑甚至死亡。水涝使根际缺氧，有氧呼吸受阻、ATP合成下降，厌氧呼吸增强，养分吸收和光合产物运输受损，气孔关闭、光合下降，叶片萎蔫衰老。耐涝植物可形成通气组织、不定根或皮层皮孔，调整能量代谢，激活活性氧（reactive oxygen species, ROS）清除系统，并由ERF-VII、MYB、WRKY等转录因子（transcription factor, TF）及乙烯、ROS、Ca²⁺信号参与。已有辣椒研究提示细胞壁生物合成和钙信号，但不同遗传背景的保守性、上游转录调控和特定代谢通路尚不清楚。因此，研究人员以耐涝纯系GC41和敏感纯系GC39为材料，整合生理与转录组，解析多层次调控，为耐涝育种提供依据。

方法上，材料来自贵州农科院辣椒研究所品种资源实验室的纯系GC41和GC39；6叶心期持续水涝9天，水位保持土表以上1–2 cm，设0、3、6、9 d和恢复3 d五个时间点，取叶片进行表型、叶绿素相对含量（soil and plant analyzer development, SPAD）、ROS、丙二醛（malondialdehyde, MDA）、渗透调节物质、抗氧化酶和AsA-GSH循环测定；RNA测序（RNA sequencing, RNA-seq）按2基因型×2处理×5时间点×3重复共60样本，Illumina NovaSeq 6000测序，HISAT2比对Zunla-1，FeatureCounts计数，DESeq2筛选差异表达基因（differentially expressed genes, DEGs），并进行基因本体（Gene Ontology, GO）/京都基因与基因组百科全书（Kyoto Encyclopedia of Genes and Genomes, KEGG）、定量实时聚合酶链反应（quantitative real-time polymerase chain reaction, qRT-PCR）和加权基因共表达网络分析（weighted gene co-expression network analysis, WGCNA）；统计用SPSS 26.0。

结果：3.1 表型分析：持续水涝下GC39中下部叶片黄化脱落、茎基吸水膨大、根基部水浸状腐烂，GC41叶色保持、叶萎蔫脱落少、根基部产生侧生不定根；恢复后GC41株高显著高于GC39，存活率95.56%显著高于GC39，表明GC41耐涝。3.2 叶片黄化与SPAD：9天水涝后GC39 SPAD由49.62降至45.48（降8.31%），GC41由51.89升至53.75（升3.59%），GC41全期显著较高，叶绿素保持能力强。3.3 H₂O₂、O₂^?和MDA动态：GC39的H₂O₂第9天峰值且显著高于GC41，O₂^?全期较高，MDA随胁迫逐步升高且显著高于GC41，说明GC41 ROS积累慢、膜脂过氧化轻。3.4 ROS清除系统活性：GC41早期超氧化物歧化酶（superoxide dismutase, SOD）快速升高并保持，过氧化氢酶（catalase, CAT）和过氧化物酶（peroxidase, POD）中后期协同增强，AsA-GSH循环关键酶活性多高于GC39，形成协调ROS清除系统。3.5 生理响应比较：GC41快速积累抗坏血酸（ascorbic acid, AsA）和谷胱甘肽（glutathione, GSH）并维持，脯氨酸（proline, PRO）、可溶性糖和可溶性蛋白显著高于GC39，说明其抗氧化和渗透调节能力强。3.6 高通量测序与组装：60样本获2 498 909 476条clean reads，Q20>95.70%、Q30>89.77%，平均比对率93.77%，数据可靠。3.7 主成分分析（principal component analysis, PCA）：同基因型和处理重复聚类紧密，PC1区分处理、PC2区分基因型，实验设计有效。3.8 差异基因表达：共鉴定4894个DEGs（3239上调、1655下调），WT vs TCK 995个，WS vs SCK 1139个，WT vs WS 2760个，耐涝与敏感基因型响应程度不同。3.9 水涝相关DEG分析：WT和WS分别有404和533个特异DEGs，42个共表达，提示共有与基因型特异响应。3.10 GO功能注释：WT vs TCK富集几丁质代谢/降解、细胞壁大分子降解、氨基糖代谢、微管相关复合体、线粒体膜间隙、序列特异性DNA结合和几丁质酶活性；WS vs SCK富集几丁质代谢/降解、药物代谢、金属离子转运、糖苷水解酶、金属离子跨膜转运和DNA结合TF活性；WT vs WS富集DNA整合、辅酶A代谢、生物刺激响应、细胞壁修饰、端粒维持、解剖结构稳态、细胞壁、胞间连丝、tRNA（m¹A）甲基转移酶复合体、ADP结合、转移酶、肽酶调节、内肽酶抑制、丝氨酸肽酶/水解酶和分子功能调节。3.11 KEGG通路富集：WT vs TCK富集11条通路，包括MAPK信号、植物激素信号转导、植物-病原互作、不饱和脂肪酸生物合成、玉米素生物合成、缬氨酸/亮氨酸/异亮氨酸降解、氨基糖和核苷酸糖代谢、酪氨酸代谢；WS vs SCK富集9条，包括MAPK、氨基糖和核苷酸糖、生物碱、甘油磷脂、植物-病原、二萜、甘油酯、维生素B6和半乳糖代谢；WT vs WS富集15条，包括玉米素、半乳糖、花青素、倍半萜和三萜、ABC转运蛋白、苯丙氨酸、植物-病原、植物激素、异黄酮、谷胱甘肽、硫代谢、酪氨酸、半胱氨酸和甲硫氨酸、色氨酸代谢及MAPK信号。结论是耐涝基因型优先能量重组、维持膜结构和协调信号-代谢-防御，敏感基因型膜脂降解和糖代谢紊乱。3.12 差异表达转录因子：WT vs TCK、WS vs SCK、WT vs WS分别鉴定61、69和68个TF，分属22、20和27个家族；WT vs WS差异家族最多，ARF、B3、GATA、TCP、WOX等见于WT vs TCK，G2-like、GeBP、HD-zip、MYB-related、Nin-like等见于WT vs WS，并包含MYB、WRKY、ERF、NAC等已知耐涝相关家族。3.13 qRT-PCR验证：随机选7个DEGs，qRT-PCR与RNA-seq表达趋势基本一致，验证转录组数据可靠。3.14 候选TF与生理指标相关性：CaWRKY11、CaWRKY32与H₂O₂负相关，CaWRKY11与单脱氢抗坏血酸还原酶（monodehydroascorbate reductase, MDHAR）极显著正相关；CaC2H2–72与SOD、POD、GSH、PRO正相关，CaC2H2–87与AsA正相关、与MDA负相关，CaHB71与PRO、AsA正相关；错误发现率（false discovery rate, FDR）校正后CaC2H2–72与SOD、GSH、PRO、SPAD仍显著，提示其为稳健候选。3.15 WGCNA：β=12构建网络，黑色模块与WT3（GC41水涝3天）显著相关（r=0.66，P<0.01）；GC41第3天快速瞬时上调、第6天回基线，GC39延迟减弱并在6–9天小峰，与抗氧化酶和渗透调节物质变化时间吻合。

讨论认为，GC41形态适应、抗氧化防御和渗透调节协同；转录组显示不饱和脂肪酸生物合成维持膜流动性、玉米素维持分生组织活性、BCAA降解供能、ABC转运蛋白促进解毒和离子稳态；敏感基因型富集膜脂降解和糖代谢紊乱。WRKY、HB、C2H2为候选上游调控因子。ERF-VII/N-end rule未显著富集，可能与合并时间点稀释瞬态信号、辣椒保守性未知及长期淹水有关；候选因子为补充而非替代。局限为仅叶片、苗期和转录证据，需功能验证和代谢组确认。模型整合抗氧化防御、渗透调节和代谢重编程。

结论部分翻译：本研究对两种耐涝性不同的辣椒基因型进行生理、生化和转录组整合分析，揭示胁迫适应的多层调控机制。耐涝基因型GC41具叶片黄化脱落少、不定根形成和顶端分生组织活性维持等形态优势，最终存活和恢复更好。生理上，GC41早期快速激活SOD、POD、CAT，长期维持AsA-GSH循环和氧化还原稳态，降低MDA；PRO、可溶性糖和可溶性蛋白积累参与渗透调节和基础代谢。转录组显示GC41重塑更广泛，富集MAPK和植物激素信号、维持膜和细胞分裂相关的不饱和脂肪酸和玉米素生物合成、能量相关BCAA降解、结构适应和解毒相关细胞壁修饰和ABC转运蛋白。Pearson相关和qRT-PCR表明WRKY、HB、C2H2转录因子与抗氧化酶和PRO积累相关，CaC2H2–72经FDR校正仍与SOD、GSH、PRO和SPAD显著相关，为核心候选。基于此提出多层工作模型，为耐涝辣椒分子育种提供候选基因和通路资源。