整合高通量表型组学与转录组学揭示水稻响应高浓度CO2与镉胁迫的转录调控机制

《iMetaOmics》：Integrated high-throughput phenomics and transcriptomics uncover the transcriptional mechanisms underlying rice responses to elevated CO2 concentration and cadmium stress

在高浓度CO₂（eCO₂）条件下降低作物镉（Cd）积累的品种选育，是当前全球农业面临的重大挑战。鉴定枢纽基因并阐明调控Cd胁迫与eCO₂响应的转录机制，需要整合性的、多组织组学研究策略。研究人员在本研究中针对水稻品种日本晴（Nipponbare），在Cd胁迫及3种eCO₂浓度条件下，跨越6个关键发育时期开展了高通量表型分析。基于图像的性状揭示了这两类胁迫下既存在共有、也存在差异性的表型响应。多组织转录组分析表明，差异表达基因（DEGs）与激素相关基因在Cd与eCO₂处理下表现出相似的表达模式。进一步地，研究构建了整合型共表达网络，并鉴定出5个胁迫响应子网络。其中，多重胁迫响应基因3（OsMSR3）被鉴定为枢纽基因，其单倍型2（Hap2）频率从低纬度地区向高纬度地区递减，并从中国南方向东北方向递减，这一格局与Cd污染土壤的地理分布趋势相呼应。总体而言，本研究揭示了水稻响应Cd与eCO₂胁迫的转录调控网络，为未来在持续升高CO₂水平背景下培育低Cd积累水稻品种提供了机制性认识与有价值的遗传资源。

该研究发表于《iMetaOmics》，聚焦全球变化与农田污染交织背景下水稻安全生产的关键科学问题，即持续升高的大气CO₂浓度与镉污染如何共同影响水稻生长、发育及分子调控。研究背景在于，气候变化导致大气CO₂浓度持续升高，而Cd作为高毒性重金属，可经土壤—植物—人体链条累积并危害健康。水稻又是典型的易积累Cd作物，因此在eCO₂条件下解析水稻Cd积累与适应反应的调控基础，对于保障粮食安全和开展低镉育种具有重要意义。既往研究虽已分别报道eCO₂对作物营养品质和Cd积累的影响，也鉴定出部分Cd转运蛋白及相关调控因子，但不同关键发育时期的动态表型响应、不同组织间的转录差异，以及Cd与eCO₂响应之间是否存在共享调控网络，仍缺乏系统阐明。正是在这一背景下，研究人员开展了高通量表型组学与多组织转录组学的整合研究。

研究人员以粳稻日本晴为材料，在自由空气CO₂富集系统（FACE）中设置3种eCO₂处理，同时利用高通量自动表型平台（HTP）开展Cd胁迫试验，对水稻从幼苗到成熟的6个关键生育时期进行连续成像和性状提取，并在拔节期与成熟期采集叶、茎及幼穗样本进行RNA-seq分析。在此基础上，研究进一步构建加权基因共表达网络（WGCNA），解析表型性状与基因模块的关联，并对候选枢纽基因OsMSR3进行群体单核苷酸多态性（SNP）与单倍型分析。研究最终表明，Cd与eCO₂虽在部分农艺后果上表现不同，但二者在转录层面存在高度重叠的响应程序；研究还鉴定出5个关键胁迫响应子网络，并指出OsMSR3是兼具分子调控意义与育种应用潜力的重要枢纽基因。该研究的重要意义在于，建立了水稻响应Cd与eCO₂的综合调控框架，为低Cd积累、适应未来高CO₂环境的分子设计育种提供了候选基因、网络资源和群体遗传依据。

研究所采用的主要技术方法包括：其一，依托位于中国河北廊坊的FACE系统和HTP平台，对日本晴在对照、Cd胁迫及3种eCO₂水平下进行全生育期高通量成像，提取54个基于图像的性状（i-traits）；其二，在拔节期和成熟期采集叶、茎及成熟期幼穗样本，开展多组织RNA测序（RNA-seq）与差异表达分析；其三，基于全部表达基因实施WGCNA，构建模块—性状关联网络并解析关键子网络；其四，结合3K Rice Genome数据库（RFGB）及RiceAtlas中国现代水稻品种队列，对OsMSR3进行SNP与单倍型地理分布分析。

在研究结果方面，论文首先在“Dynamic phenotypic changes in rice under Cd and eCO₂ conditions”部分系统揭示了Cd与eCO₂处理下水稻表型的时序变化。研究人员在31个时间点采集表型数据，覆盖幼苗生长、分蘖、拔节、孕穗、开花和成熟6个阶段。结果显示，i-A（基于图像的面积）和i-Cot（俯视紧致度）在不同处理间自中后期开始出现显著分化，并在成熟期表现出共同下降趋势，说明Cd与eCO₂均可影响植株空间构型与生长终点状态。但二者也存在明显差异：i-MARA（最小外接矩形面积）在Cd条件下升高、在eCO₂条件下降低，而i-Co（紧致度）则在eCO₂条件下升高、在Cd条件下降低。人工测定进一步表明，成熟期生物量和分蘖数在3种eCO₂处理下均显著增加，而在Cd胁迫下显著下降。由此可见，Cd与eCO₂诱导了既共享又分化的动态表型响应，高通量图像性状可有效量化胁迫强度及其发育阶段特异性效应。

在“Multi-tissue transcriptomes reveal shared transcriptional functions in response to Cd and eCO₂”部分，研究通过多组织转录组比较解析了两类胁迫下的分子共性与组织特异性。研究人员选择拔节期和成熟期作为代表性营养生长与生殖生长阶段，对叶、茎和幼穗进行RNA-seq。结果发现，eCO₂条件下不同组织和时期的差异表达基因数量差异显著，整体上拔节期多于成熟期，成熟期幼穗则表现出较高数量的响应基因。聚类分析表明，不同表达簇分别富集于光合作用、防御反应、氮代谢调控及碳水化合物代谢等功能，说明eCO₂影响具有明显的组织与发育阶段依赖性。更重要的是，Cd与eCO₂响应基因之间存在大范围重叠：Cd处理中85.64%的差异表达基因与eCO₂响应基因重合，且有1540个基因在两种胁迫及两个时期中共同差异表达。这些共有基因主要富集于“响应刺激”“光合作用—捕光”及“碳水化合物代谢过程”等生物学过程，说明水稻在应对Cd毒性与适应eCO₂时，激活了相当程度重叠的转录功能通路。

在“Dynamic changes in hormone-related genes under Cd and eCO₂ conditions”部分，研究进一步从植物激素角度解析胁迫应答的调控复杂性。基于MapMan注释，研究人员共识别出682个与8类主要激素相关的表达基因，包括生长素、乙烯、脱落酸（ABA）、油菜素内酯（BR）、细胞分裂素（CK）、茉莉酸（JA）、水杨酸（SA）和赤霉素（GA），并将其划分为合成与降解、信号转导、诱导响应3类功能模块。结果显示，大多数激素相关基因在Cd与eCO₂处理下具有很强的组织特异性，且不少基因在两种胁迫中的表达趋势相似。例如，SA和JA的部分合成降解基因、BR和CK的部分信号转导基因，以及生长素和ABA的部分诱导响应基因，在拔节期茎组织中特异高表达并呈现相近变化。OsOPR1等JA相关基因在Cd和eCO₂处理中均下调，提示两类胁迫可能共享部分激素调控路径。但乙烯相关基因呈现差异性模式，尤其成熟期茎组织中部分诱导响应基因在eCO₂下上调、在Cd下下调，表明乙烯信号可能参与区分两类环境因子的特异应答。

在“Construction of an integrative network in response to Cd and eCO₂ in rice”部分，研究人员利用WGCNA整合多组织转录组与表型数据，构建了水稻响应Cd与eCO₂的共表达网络。通过设置软阈值17，最终得到31个共表达模块，并将模块特征值与性状及处理条件进行相关分析。结果显示，蓝色、黄色和绿松石色模块分别与叶、幼穗和茎组织显著相关，体现出强烈的组织特异性。深红色模块则与Cd胁迫呈正相关，而与eCO₂呈负相关。功能富集分析表明，这些模块分别涉及光合作用、刺激响应、初级代谢调控、蛋白质代谢、氧化胁迫响应、微管过程、细胞死亡和防御反应等生物学功能。已知Cd响应基因如OsCLT1和OsMSR3位于黄色模块中，支持该网络具有较高的生物学可信度。该部分结果说明，多组织共表达网络能够有效连接胁迫表型与分子调控层次，为识别关键调控因子提供了框架。

在“Transcriptional regulatory subnetworks of rice respond to Cd and eCO₂”部分，研究在关键模块内进一步提取了以已知Cd响应基因为中心的调控子网络。总共识别出5个子网络，涵盖9个已报道基因。黄色模块中存在两个子网络，分别以OsMSR3和OsCLT1为枢纽；其中OsMSR3共表达100个基因，OsCLT1共表达4个基因，二者网络彼此独立，提示在生殖发育阶段可能承担不同的Cd解毒功能。蓝色模块中，OsMTI-3a与OsCADT1共同构成叶片相关子网络，关联60个共表达基因，指向叶组织中重金属耐受和活性氧（ROS）缓冲的协同调节。深红色模块中，OsZIP1作为枢纽基因连接100个靶基因，体现其与Cd特异响应密切相关。绿松石色模块则包含OsPCS1、OsIRT1、OsProT1和OsHIPP42四个金属相关基因，其中OsPCS1、OsIRT1与OsProT1共享多个靶基因，而OsHIPP42主要与OsPCS1共享靶标，说明这些金属稳态与解毒因子之间既存在共同调控，也有相对专一的连接方式。该结果构成了论文最核心的机制框架，即水稻对Cd与eCO₂的响应不是孤立基因事件，而是由组织特异性的枢纽网络共同驱动。

在“A natural OsMSR3 variation may contribute to subspecies differences in rice Cd and eCO₂ tolerance”部分，研究人员围绕OsMSR3开展了群体遗传学分析，以评估其育种应用潜力。首先，OsMSR3表达量与多个关键图像性状显著负相关，包括i-A、i-MARA和i-MECD，提示该基因与胁迫下植株形态变化密切相关。随后，研究在OsMSR3编码序列中筛选非同义SNP，确定一个位于第50位的T/C变异，该变异导致第17位氨基酸由L变为P，并据此将2887份水稻材料划分为Hap1和Hap2两种单倍型。Hap2主要分布于籼稻（indica）群体，表现出明显的亚种偏向性。在全球地理分布和中国现代品种分析中，Hap2频率均呈现从低纬度向高纬度递减、从中国南方向东北方向递减的趋势。这一变化格局与Cd污染土壤的空间分布趋势相呼应。基于原文表述，研究认为该单倍型变异可能是影响水稻Cd与eCO₂适应性的一个重要分子标记，并为未来低Cd积累水稻育种提供了候选位点。

讨论部分主要围绕四个层面展开。首先，论文强调整合高通量表型组学与多组织转录组学，能够构建水稻在Cd与eCO₂条件下的表型—转录组综合网络，从而突破过去仅关注单基因或单组织研究的局限。其次，高通量表型结果表明，图像性状比传统人工性状更能细致刻画不同胁迫在时序与空间构型层面的差异，为后续关联分析和功能基因挖掘提供了高分辨率表型基础。再次，差异表达基因与激素相关基因的大规模重叠，说明Cd与eCO₂在水稻中调动了相当程度共享的转录应答通路，但乙烯等信号途径又体现出特异性分化，提示胁迫应答兼具保守性与可塑性。最后，OsMSR3的网络中心地位及其单倍型的群体地理分布，凸显了该基因在理解Cd胁迫适应和开展分子育种中的潜在价值。论文同时指出，部分预测到的共表达关系仍有待后续实验验证，这反映出该研究提供的是高可信度的网络框架与资源基础。

研究结论部分可译述为：本研究整合高通量表型组学与多组织转录组学，阐明了水稻响应Cd与eCO₂胁迫的转录调控机制。研究记录了i-A和i-Co等关键图像性状在两类胁迫下的动态变化，并观察到eCO₂条件下分蘖数和生物量显著增加，而Cd条件下则明显下降。差异表达基因及激素相关基因呈现动态且具组织特异性的表达模式。进一步构建的共表达网络鉴定出与表型性状相关的枢纽模块和调控子网络，并发现了一个频率自低纬度向高纬度递减的优良OsMSR3等位变异。总体而言，这些结果揭示了一个调控水稻响应Cd与eCO₂的综合性转录调控框架，并为在持续升高的大气CO₂水平下培育低Cd积累水稻品种提供了机制依据和重要遗传资源。