多组学整合分析揭示田间高粱光合作用与抗旱性之间的调控关联

《Plant, Cell & Environment》：Integrated Multi-Omic Analyses Uncover a Regulatory Link Between Photosynthesis and Drought Tolerance in Field-Grown Sorghum

　　

全球变暖与日益加剧的水资源短缺对农业构成重大挑战，凸显了将逆境生物学洞见转化为抗旱品种培育的迫切需求。本研究评估了6份遗传背景多样的高粱（Sorghum bicolour）材料在美国亚利桑那州中部田间干旱条件下的表现，在7周时间内整合了生理学、转录组学及氧化胁迫测定数据。研究人员利用基于共表达相关性、转录因子（TF）结合基序特征及蛋白质-蛋白质相互作用（PPI）的网络分析，鉴定出一个与光合能力强烈相关的干旱关联模块。该模块包含一个胁迫响应转录因子SbDof8（本文中称为SbCDF2/3-like，或SbCDF2/3L）作为高连接度枢纽基因，且CDF2/3相关结合基序在共表达模块成员的启动子区域显著富集。这些共表达成员富集于胁迫响应、代谢及光合作用相关过程，且在耐旱材料中较敏感材料在干旱下持续维持更高表达水平。对包含两个独特材料的独立高粱干旱时序数据集的分析显示，SbCDF2/3L与光合作用相关模块基因在耐旱与敏感基因型间呈现一致的表达模式，增强了该调控模块稳健性的证据。综上所述，研究结果突出了一类光系统I（PSI）相关基因——包括捕光蛋白、PSI亚基以及一个重要的干旱响应转录调控因子——在抗性材料中协同上调以应对干旱响应。这些数据为改良高粱田间条件下抗旱性与生物量产量的育种靶点提供了重要参考。

全球气候变暖与水资源短缺对农业生产构成严峻挑战，将逆境生物学研究成果转化为抗旱品种培育策略已成为当务之急。高粱作为世界第五大谷类作物，虽具有适应干旱环境的固有能力，但不同材料间抗旱性存在显著遗传变异，其调控机制尚不明确。现有研究多基于单一材料或受控环境条件，缺乏田间多头材料整合多组学数据的系统性分析，尤其难以解析上游调控因子与下游防御基因及产量性状之间的关联。此外，Dof转录因子家族虽被报道参与光合 machinery 调控与环境信号整合，但具体成员功能及其在作物田间抗旱中的调控网络仍有待阐明。

本研究发表于《Plant, Cell & Environment》。研究人员以6份遗传背景多样、抗旱性各异的高粱材料为对象，在美国亚利桑那州田间干旱条件下开展系统研究，整合转录组、氧化胁迫状态、光合速率及全株表型等多维度时序数据，构建共表达网络并结合转录因子结合基序分析与蛋白质互作数据，鉴定出以SbCDF2/3L为核心的干旱-光合调控模块，揭示了其在耐旱材料中维持光合基因表达、保障产能的分子机制，为高粱抗旱育种提供了关键靶标基因。

研究采用的关键技术方法包括：田间试验设计与样本队列（6份高粱材料，设置充分灌溉与限量灌溉处理，7周时序采样）；LI-6800便携式光合测定系统获取光合参数；氧化胁迫相关生化指标测定；RNA-seq转录组测序与数据处理；加权基因共表达网络分析（WGCNA）构建共表达网络；基于iTAK的转录因子与蛋白激酶预测；利用Arabidopsis DAP-seq数据库进行启动子基序富集分析（AME）；STRING数据库蛋白质-蛋白质相互作用整合；GO与KEGG功能富集分析；以及独立公共数据集（Varoquaux等，2019）的验证分析。

研究结果部分如下：

**3.1 高粱材料对干旱的多样表型响应**

研究首先确认了田间干旱处理的有效性。土壤体积含水量（SVWC）监测显示限量灌溉组维持在约16%，显著低于充分灌溉组的24%。表型分析表明，材料'A'（M1）在干旱下呈现最明显的叶片褐变，且其根际内生真菌群落组成在控水与干旱间差异较小。生物量测定显示，'A'和'D'（N290-B）在两种条件下总生物量均最低；'A'的根冠比（RSR）在干旱下降幅最大，而'C'（80M）和'F'（SC391）则呈上升趋势。光合效率参数Fv'/Fm'与净CO₂同化速率（Anet）测定表明，'A'在第五周和第七周表现出明显的干旱与高温胁迫特征，而'B'（DL/59/1530）和'C'更为稳健。综合形态与生理数据，将'B'和'C'归为耐旱型，'A'和'F'归为敏感型。

**3.2 高粱材料干旱胁迫下的分子响应变异**

转录组主成分分析（PCA）显示基因型是样本聚类的最大驱动因素，其次是处理条件和时间点。差异表达基因（DEG）分析揭示不同材料响应模式高度异质：'A'和'E'（SC1271）DEG数量最多，而'D'和'F'最少；'B'、'C'、'D'几乎不富集应激响应GO条目，与'A'、'E'形成鲜明对比。氧化胁迫代谢物与酶活相关性分析显示，抗氧化代谢物与部分抗氧化酶呈正相关，而氧化胁迫标记物与抗氧化酶呈部分负相关。第五至第七周，'A'和'F'表现出更高的羟基丙酮酸还原酶（HPR）和乙醇酸氧化酶（GOX）积累，以及更低的抗氧化酶/代谢物水平，表明其遭受更严重氧化胁迫，与持续上调的应激响应DEG模式一致。

**3.3 干旱相关共表达转录本群的鉴定**

基于WGCNA构建的共表达网络包含22,314个转录本，分为44个模块。13个模块与性状显著相关，其中'Darkmagenta'模块与光合参数及氧化胁迫标记物相关性最强，且基因显著性（GS）与模块成员度（MM）呈正相关，表明性状相关基因在模块中处于核心位置。该模块富集于非生物胁迫响应与代谢过程，AME基序分析鉴定出10个转录因子结合基序显著富集，其中4个注释为非生物胁迫调控因子（bZIP、ERF及CDF家族）。SbCDF2/3L展现出最高的MM和GS，提示其为连接转录调控与光合性状变异的关键枢纽。

**3.4 SbCDF2/3L结合基序关联基因的干旱响应模式**

系统发育分析将SbCDF2/3L（原注释为SbDof8）归入Cycling Dof（CDF）亚家族，与拟南芥CDF2（AT5G39660）和CDF3（AT3G47500）为姊妹群。'Darkmagenta'模块1,093个基因中，681个（62%）启动子区域含有SbCDF2/3L预测结合位点，其核心基序为5'-T/AAAAG-3'，与玉米及拟南芥已报道的CDF结合基序一致。表达模式分析显示，'A'材料中SbCDF2/3L在第三至第七周急剧下降。SbCDF2/3L与其模块内预测靶基因的相关系数显著高于模块内非靶基因及模块外随机基因，且其靶基因在干旱下的log₂倍变变化显著高于非靶基因。整合TF结合、PPI数据构建的186基因子网络中，129个为SbCDF2/3L预测靶标，涉及净CO₂同化、类黄酮含量、叶温及总抗氧化能力等性状，多个基因与已报道的CO₂同化数量性状位点（QTL）重叠。

**3.5 公共 field 胁迫数据中SbCDF2/3L调控响应的验证**

利用Varoquaux等的公共数据集（后开花期干旱，8周时程），研究人员比较了耐旱材料BT×642与敏感材料RT×430中'Darkmagenta'核心基因的表达模式。结果显示，两类材料中耐旱型均维持更高的光系统I、胁迫响应及代谢过程相关基因表达，而敏感型显著下调。SbCDF2/3L在耐旱型中表达水平显著高于敏感型，且含CDF2/3结合基序的基因在耐旱型中呈现更高的干旱诱导表达。这些跨实验验证结果支持SbCDF2/3L作为高粱光合与胁迫响应基因正向调控因子的作用。

**讨论**

研究人员综合讨论指出，干旱胁迫通过复杂表型与分子响应降低作物生物量与产量，而高粱虽适应干旱环境，其材料间遗传变异背后的调控机制尚不清晰。本研究通过多头材料田间多组学整合，鉴定出以SbCDF2/3L为核心的关键调控模块，其特色在于：该转录因子将环境信号与光合 machinery 转录调控耦合，在耐旱材料中维持PSI相关基因（LHCA蛋白、PsaO/K/L亚基等）的高表达，从而保障电子传递与光合复合体稳定性。讨论中提到，这类光合作用相关基因的多层协调调控（涉及捕光、电子传递、环式电子流及叶绿体调控）并非单一通路受限，而是整个光合系统的协同响应。SbCDF2/3L的独特性在于它是该模块中唯一的CDF/Dof家族成员，且其靶基因启动子富集、表达相关性及跨实验稳健性均支持其作为耐旱正调控因子的功能。研究人员也指出，尽管拟南芥CDF的昼夜节律特性已在少数禾本科Dof转录因子（如RDD1）中得到证实，SbCDF2/3L是否具有类似节律调控特征及其与胁迫下生长整合的关系仍需探究。最终，研究提出通过调控SbCDF2/3L在干旱条件下的表达以同时提升抗旱性与光合能力，为高粱抗旱高产育种提供了可操作的分子框架。

**研究结论**

研究人员提出概念模型：SbCDF2/3L在干旱条件下结合并调控模块内胁迫响应基因，包括一系列核编码的光合作用相关基因，这些基因与LI-6800测定性状相关，且已知在碳固定（CO₂/HCO₃^?代谢）、光合复合体组装、电子传递及叶绿体水平光合性能调控中发挥作用。增强SbCDF2/3L介导的转录调控有望在未来成为维持限量灌溉条件下光合性能、提升抗旱性、生物量积累及产量稳定性的重要途径。