谷子耐盐性的双重进化策略：转录组学与重测序的整合分析

《Frontiers in Plant Science》：Transcriptomics and resequencing reveal dual evolutionary strategies for salt tolerance in foxtail millet (Setaria italica (L.))

　　

摘要 引言： 土壤盐渍化对谷物生产构成重大制约，然而谷子（Setaria italica (L.)）耐盐性的遗传基础在很大程度上仍不清楚。 方法： 为了剖析这些耐受机制，研究人员整合了表型分析、生理测定、转录组学和群体基因组学，从超过500份种质资源中筛选出18个极端耐盐群体（STP）和18个极端盐敏感群体（SSP）的材料，并在NaCl处理下进行评估。 结果： 在盐胁迫下，与SSP材料相比，STP材料表现出更好的生长状况和更高的叶绿素相对含量（RCC：第2天为5.9 vs. 3.0），同时膜损伤程度更低（相对电导率REC：14% vs. 24%）和活性氧（ROS）积累更少。转录组分析鉴定出一个与耐盐性相关的模块STP007，该模块富含AP2/ERF和IDD基序，并包含候选转录因子，如YABBY和WRKY。群体基因组学揭示，尽管遗传背景相对同质，但第3和第6号染色体上的局部选择性清除与耐盐性相关。通过整合转录组互作网络与群体进化特征，研究人员鉴定了五个参与脱落酸（ABA）生物合成和氧化还原稳态的核心候选基因，突显了耐盐性可能的遗传基础。 讨论： 总之，这些发现揭示了谷子耐盐性可能存在一种双层机制，可能由核心基因组节点的局部适应性选择和转录网络内的动态可塑性协同调控。该模型为缓解生长与胁迫之间的权衡提供了见解，并为培育抗逆性强的谷类作物提供了潜在的候选靶点。

谷子耐盐性的双重进化策略：转录组学与重测序的整合分析解读

1. 研究背景、问题与研究意义

土壤盐渍化是全球范围内限制作物生产力的主要非生物胁迫之一。过量的盐分会对植物造成离子毒害和渗透胁迫，导致生长抑制、代谢紊乱和产量下降。因此，阐明植物耐盐性的分子机制对于提高作物抗逆性和保障粮食安全至关重要。

谷子作为一种新兴的C₄作物模型，具有基因组小（约510 Mb）、遗传资源丰富、生命周期短等优点，且中国拥有丰富的谷子种质资源。尽管在功能基因组学方面取得了显著进展，例如通过时间序列RNA-seq和加权基因共表达网络分析（WGCNA）揭示了耐盐性可能源于协调的调控网络，但在群体水平上驱动耐盐性分化的进化驱动力仍未得到充分探索。一个关键的科学问题是：耐盐性分化主要是由编码序列的适应性进化、胁迫诱导的调控网络重编程，还是两者的协同作用驱动的？为了解决这一问题，本研究旨在通过整合表型组、生理学、转录组学和群体基因组学数据，系统性地阐明谷子耐盐性的遗传架构。

2. 研究方法概述

研究人员从超过500份谷子种质资源中，筛选并构建了由18个极端耐盐群体（STP）和18个极端盐敏感群体（SSP）组成的关联群体。采用的主要关键技术方法包括：

1. 1.
   生理生化分析： 通过DAB和NBT组织化学染色检测活性氧（ROS）积累，并测量相对电导率（REC）和相对叶绿素含量（RCC）评估细胞膜损伤和光合能力。
2. 2.
   转录组测序与分析： 对36份材料的叶片组织进行RNA-seq，利用DESeq2进行差异表达基因（DEGs）分析，并通过多WGCNA（multiWGCNA）构建基因共表达网络。
3. 3.
   全基因组重测序与群体遗传学分析： 对36份材料进行全基因组重测序，利用ADMIXTURE和主成分分析（PCA）进行群体结构分析，并通过固定指数（F_ST）和跨群体复合似然比（XP-CLR）方法进行选择性清除扫描。

3. 研究结果

3.1 STP与SSP的生理分化

通过田间实验和实验室筛选，研究人员从500多份材料中选出了18个高度耐盐（STP）和18个极度盐敏感（SSP）品系。在300 mM NaCl处理下，STP幼苗能保持绿叶和直立形态，而SSP则表现出严重的叶片黄化和萎蔫。生理指标测定表明，STP能维持更高的RCC和更低的REC，显示出更强的光合能力和细胞膜稳定性。

3.2 STP与SSP幼苗叶片组织化学染色

通过DAB和NBT染色评估抗氧化能力。正常条件下，两组间H₂O₂和超氧阴离子（O₂^-）含量无显著差异。盐胁迫后，STP叶片染色较浅，表明其ROS积累较少，具有更强的活性氧清除能力；而SSP则染色深，表明ROS大量积累，耐盐性较差。

3.3 盐胁迫下STP与SSP的生理响应

RCC测量结果显示，处理前两组基线值相似。处理后第2天和第4天，STP的RCC均显著高于SSP。REC测量结果同样显示，处理前两组无差异，但处理后SSP的REC显著高于STP，表明STP能更好地维持细胞膜完整性。

3.4 STP与SSP的转录差异及调控网络分析

RNA-seq分析鉴定了2365个具有显著基因型-处理交互作用的基因。主成分分析（PCA）显示盐处理是转录变异的主要驱动力。GO富集分析表明，耐盐群体的DEGs富集于广谱非生物胁迫响应通路，而盐敏感群体的DEGs则偏向生理重塑。交互作用DEGs特异性地富集于茉莉酸介导的防御反应、苯丙烷类次生代谢等通路。

3.5 耐盐候选基因的表达模式聚类与功能富集

将交互作用基因的表达模式分为五类，其中类别I（594个基因）在耐盐群体中特异性响应，被确定为关键候选基因簇。该类基因主要富集于两个核心功能模块：核RNA监控和脂质磷酸化驱动的气孔运动，表明转录后调控和气孔防御是STP的主要耐盐机制。

3.6 WGCNA揭示群体特异性共表达模块和调控网络

WGCNA分析鉴定出STP特有的STP007模块（含81个基因），该模块在SSP网络中不保守。该模块包含三个核心转录因子（TF）：一个YABBY和两个WRKY。转录因子结合位点（TFBS）分析显示，AP2/ERF和IDD家族基序在该模块启动子区显著富集，暗示乙烯和ABA信号是该模块的主要调控因子。

3.7 STP与SSP的遗传分化和群体结构

全基因组重测序获得约250万个高质量SNP。群体结构分析和PCA均未显示STP和SSP之间存在显著的遗传分化，表明两者遗传背景高度同质。尽管如此，F_ST和XP-CLR分析一致地鉴定出第3和第6号染色体上存在强烈的选择性清除信号。

3.8 候选基因的交叉验证与功能整合

通过将转录组交互网络与选择性清除信号交叉比对，研究人员鉴定出5个受强烈多样性选择的关键基因，涉及ABA生物合成（Seita.6G138300，CNX2直系同源基因）、氧化还原稳态（Seita.2G037700，硫氧还蛋白H型）和胁迫感知（Seita.7G058300，CRK同源基因）。相反，STP007模块中的81个共表达基因在相同阈值下完全不存在选择性清除信号。等位基因频率和单倍型分析进一步证实了这些基因在STP中存在显著的突变等位基因富集。

4. 总结与讨论

讨论总结： 研究表明，谷子的盐适应可能由复杂的多基因调控网络驱动，而非组成型表达差异。耐盐性与转录响应重编程密切相关，其中转录后质量控制与物理屏障响应的协同激活是STP的主要防御策略。STP007模块作为一个耐盐特异性共表达单元，可能通过WRKY和YABBY转录因子形成调控级联，并由乙烯和ABA信号主导。尽管遗传背景高度同质，但第3和第6号染色体上的强选择性清除表明，盐驱动的定向选择精准靶向了少数局部基因组区域。这种由核心节点适应性变异与下游网络动态可塑性组成的双层架构，可能有助于缓解生长稳定性与胁迫灵活性之间的权衡。

研究结论翻译：

总之，本研究阐明了谷子耐盐性的一种双层进化和调控机制。该机制由两种不同的策略协同驱动：首先，核心基因组节点上适应性等位基因的方向性选择和富集，维持了基础的生理稳态；其次，下游转录调控网络的高度动态可塑性，介导了对环境胁迫快速且可逆的响应。这种双机制模型可能有助于缓解植物生长与抗逆性之间的进化权衡。同时，这些发现为工程化培育高产、稳定且气候适应性强的谷类作物提供了有前景的候选靶点。