撒哈拉来源 Sinorhizobium meliloti IRAMC:0087 的基因组可塑性、胁迫耐受性与独特 pSymA 共生适应

《Frontiers in Sustainable Food Systems》：Genomic plasticity, stress tolerance, and a distinctive pSymA in the symbiotic bacterium Sinorhizobium meliloti IRAMC:0087

　　

干旱生态系统促进根瘤菌菌株（rhizobial strains）的出现，这些菌株适应严重非生物胁迫（abiotic stress）同时保留共生能力（symbiotic competence）。在此，研究人员表征了 IRAMC:0087，该根瘤菌菌株分离自突尼斯南部撒哈拉灌木 Calobota saharae 的根瘤。表型试验显示其对热、盐度和渗透胁迫具有高耐受性，并具有植物促生性状（plant growth-promoting traits），包括吲哚-3-乙酸（indole-3-acetic acid, IAA）产生、磷酸盐溶解和铁载体（siderophore）分泌。体内试验证明其对 Vachellia tortilis 的有效固氮结瘤，而 Medicago 和 Trifolium 物种虽通过非共生机制被促进生长，却未被结瘤。共聚焦显微镜确认了根部定殖和根瘤占据，表明宿主依赖性侵染行为。为更好理解该菌株的生理行为，研究人员生成了其完整基因组序列，揭示了一个 7.37 Mb 的环化基因组，由五个复制子（replicon）组成：一个 3.65 Mb 染色体、两个对应于 pSymA（0.59 Mb）和 pSymB（1.67 Mb）的大质粒，以及两个附属质粒（accessory plasmids），总体 GC 含量为 62%。系统发育基因组学（phylogenomics）、平均核苷酸一致性（Average Nucleotide Identity, ANI）和数字化 DNA–DNA 杂交（digital DNA–DNA hybridization, dDDH）分析将 IRAMC:0087 归属于 Sinorhizobium meliloti。然而，比较基因组学揭示 pSymA_IRAMC:0087 与 S. meliloti 1,021 存在显著结构和功能分化，而与 S. garamanticum 和 S. numidicum 相似性更高。NodC 系统发育将该菌株置于 symbiovar tunisiaense 内，与其独特宿主范围一致。基因组注释进一步鉴定出与胁迫耐受、根际竞争力、生物活性代谢物产生、信号传导和共生相关的决定因子，包括位于附属质粒上的基因。这些发现突出 IRAMC:0087 作为一种沙漠适应根瘤菌，具有独特共生特性和强生态适应性，为理解干旱环境中根瘤菌进化提供了见解。

干旱与半干旱生态系统中，本土豆科植物及其根瘤菌（rhizobia）是未曾充分挖掘的遗传与功能多样性资源。根瘤菌可与豆科植物建立共生关系，将大气氮（N₂）转化为氨，减少合成氮肥依赖并改善土壤肥力。然而，附属复制子多样性对胁迫适应和共生性状分布与进化的贡献仍理解不足。突尼斯南部撒哈拉灌木 Calobota saharae 根瘤中分离的 Sinorhizobium meliloti IRAMC:0087 所处生态位具有极端温度、低水可利用性和高土壤盐分。研究人员结合表型、共生和基因组分析，旨在表征其适应特征，鉴定与生存、竞争力和胁迫下共生表现相关的基因组决定因子，并比较其与近缘 S. meliloti 菌株的基因组结构，以评估染色体外变异对环境适应和宿主特异性的贡献。该研究为理解干旱环境根瘤菌基因组进化及筛选耐逆生物接种剂提供了依据。

关键技术方法概括：菌株 IRAMC:0087 于 2018 年 3 月分离自突尼斯南部撒哈拉灌木 Calobota saharae 根瘤。研究人员开展热、盐、渗透、抗生素和重金属耐受表型测定，植物促生性状筛选，六种豆科植物结瘤试验，mRFP（单体红色荧光蛋白）标记与共聚焦显微镜观察；采用 Illumina NextSeq500 和 Oxford Nanopore MinION 测序，Flye、Medaka、Polypolish、POLCA 组装，PGAP 和 RAST 注释；比较基因组分析使用 TYGS dDDH、EzBioCloud ANI、Roary、PGPg_finder、Mauve、eggNOG-mapper、BlastKOALA、ISfinder、nucmer、TYGS 和 MEGA11；统计采用 SPSS 的 ANOVA、Tukey 和 t 检验。以上方法概括约 250 字以内。

3.1 S. meliloti IRAMC:0087 表型特征。通过生长测定，IRAMC:0087 在 40 °C 连续热胁迫下保留 77.16% 生长，47 °C 和 50 °C 热激后分别保留 90.63% 和 72.06%；在 340 mM NaCl 下保留 79.17%，510 mM NaCl 下保留 62%；20% PEG 6000 使生长较对照降低 46%。该菌株耐庆大霉素 2.5 μg/mL、卡那霉素 50–100 μg/mL、氨苄青霉素 100 μg/mL、新霉素 50–150 μg/mL、氯霉素 20 和 30 μg/mL，但四环素 10 μg/mL 及链霉素 20 和 30 μg/mL 完全抑制生长。重金属耐受方面，其耐受 MnCl₂ 至 2 mM、CoCl₂ 0.05 mM、CuCl₂ 0.05 mM、NiCl₂ 0.025 mM 和 Pb(CH₃COO)₂ 0.025 mM。促生性状中，IAA 产量为 21.05 ± 0.16 μg/mL，铁载体产生指数为 1.46 ± 0.03，溶磷阳性，自由生活条件下未检测到 ACC 脱氨酶活性。

3.2 共生互作与根定殖能力。宿主范围试验显示，IRAMC:0087 可结瘤 Vachellia tortilis，形成粉红色根瘤，并显著提高地上部干重，支持其有效固氮；Calicotome villosa、Medicago arborea、M. sativa、M. truncatula 和 Trifolium subterraneum 均未结瘤。在 M. sativa 中，接种使根干重较未接种增加约 48%；在 M. truncatula 中，地上部干重和根干重分别增加 68% 和 153%，提示通过非共生机制促进生长。共聚焦显微镜显示，该菌株可定殖不同豆科植物根系：在 V. tortilis 中先定殖主根表面并富集于次生根发生处；在 M. sativa 中粘附于主根表面和根毛基部但未见侵入；在 T. subterraneum 中定殖根毛和胞间空间，提示启动内生互作而非结瘤共生。根瘤表面灭菌后回收菌落经 16S rRNA 基因测序确认 IRAMC:0087 占据根瘤。

3.3.1 一般基因组特征与分类归属。长读长 Nanopore 与短读长 Illumina 联合测序获得 IRAMC:0087 完整环化基因组，总大小 7,367,083 bp，G+C 含量约 62%，由 1 个染色体和 4 个质粒组成，包括 pSymA、pSymB、p_A1 和 p_A2。注释得到 6,950 个编码 DNA 序列，其中 6,881 个蛋白编码基因，326 个假基因，6 个 rRNA 基因，56 个 tRNA 基因和 4 个 ncRNA 基因。16S rRNA 基因与 S. meliloti LMG 6133^T 100% 一致；ANI 为 98.96%，与 S. meliloti USDA1002^T 的 dDDH 为 91.60%，超过物种界定阈值。

3.3.2 泛基因组与植物促生活性比较。对包含 71 株 S. meliloti 及两个外群的分析，共鉴定 37,010 个直系同源基因，其中核心基因 497 个（1.3%），软核心 3,453 个（9.3%），壳基因 2,878 个（7.8%），云基因 30,182 个（81.6%），显示高度开放泛基因组。IRAMC:0087 聚类于 S. meliloti 主分支，并具有自身附属基因子集。植物促生性状基因家族比较显示，IRAMC:0087 携带 4,520 个 PGPT 基因家族，比 S. meliloti 内群均值 4,399 ± 217 多约 121 个，并高于外群。PCA 支持其整体 PGPT 谱与 S. meliloti 一致，但在内群中形成独立短支。通路水平上，IRAMC:0087 富集群体感应/生物膜形成（131 vs. 102.1 ± 10.7）、细菌分泌系统（79 vs. 49.8 ± 13.1）、细菌适应性基因（534 vs. 437.6 ± 64.8）和胞外多糖产生（48 vs. 40.7 ± 2.1），并富集定殖表面附着、根定殖、通用胁迫响应及若干植物信号/植物激素类别；氮获取减少（161 vs. 191.6 ± 13.8），重金属解毒、异源物降解等也降低。结果表明其功能重组偏向竞争和定殖，部分经典共生基因减少。

3.3.3 与模式菌株 S. meliloti 1,021 的基因组结构与功能分化。比较显示，S. meliloti 1,021 基因组为 6.7 Mb、3 个复制子，而 IRAMC:0087 为 7.3 Mb、5 个复制子；染色体和 pSymB 类染色体大小及 G+C 含量相似，但 pSymA_IRAMC:0087 较 pSymA_1,021 小约 2.5 倍。IRAMC:0087 还携带 p_A1（1,247,198 bp）和 p_A2（197,378 bp）。Mauve 比对显示染色体和 pSymB 保守，pSymA 核苷酸组织差异明显，p_A1 与 pSymA_1,021 共享多个局部共线块。重排断点沿 p_A1 呈热点分布，插入序列家族拷贝数差异提示 IS 介导重组可能参与 pSymA 来源序列的重排和再分布。COG 分析中，IRAMC:0087 在碳水化合物运输与代谢、转录、复制重组修复、细胞壁/膜/包膜生物发生、信号转导、胞内运输分泌和防御机制等类别基因数较高。KEGG 分析显示其富集环境信息处理、信号与细胞过程、氨基酸代谢、聚糖生物合成与代谢；S. meliloti 1,021 则富集能量代谢、辅因子和维生素代谢、异源物降解，并具有完整反硝化途径而 IRAMC:0087 缺失。总体表明 IRAMC:0087 在环境感知、细胞响应、结构防御和根际适应方面投入更多。

3.3.4 胁迫耐受与植物相关性状的基因组决定因子。基因组注释发现热休克蛋白及调节因子，如 HSPQ、HSP20、HSP33、HSP70、HrcA、GrpE，以及 DnaK/DnaJ、ClpA/ClpB、GroEL-GroES 伴侣系统；rpoH 编码的 σ³² 位于染色体。冷适应相关基因包括 ltrA 和冷休克蛋白基因。盐和低水势响应中，存在胆碱、甘氨酸甜菜碱和海藻糖合成与摄取基因，如 betABCI、choVWX、betB、betI 及 otsA、treY、treZ、treS；pSymB 携带 exoFYXUVWTILAMOP 胞外多糖合成操纵子，mucR 有 4 个拷贝。氧化胁迫相关基因包括过氧化氢酶、过氧化物酶、超氧化物歧化酶、MsrA-MsrB 及 oxyR、soxR。抗生素抗性相关决定因子包括 MacAB-TolC 外排泵、β-内酰胺酶家族基因、cml、氯霉素磷酸转移酶和 sbmA；重金属耐受相关基因包括 MerR 样调节因子、ABC 型铁锰转运体、镍钴转运体、CorA、HlyC/CorC、P 型 ATP 酶和 ArsR/SmtB 家族调节因子。植物促生相关基因包括色氨酸依赖 IAA 合成基因 trpAB、trpC、trpD、trpEG；磷酸动员基因 phoD、ppx、ppk2、phnFGIJKLM、phoB、phoR、phoU、pstABC 及 pqq 操纵子；铁载体合成与摄取基因主要位于附属质粒，p_A1 携带 rhbABCDEF、rhtA、rhrA、RhtX，p_A2 携带 lucA/lucC 家族基因和 TonB 依赖性铁载体转运体；pSymA 携带 acdS，但自由生活条件下未检测到 ACC 脱氨酶活性。根际竞争与通讯相关基因包括趋化 cheA、cheB、cheD、cheR、cheW，运动 fli/flh 基因簇，群体感应 TraI/TraR 系统，LysR 家族调节因子及 Sec、T1SS、T2SS、T3SS、T4SS 分泌系统。

3.3.5 pSymA 进化与共生基因组织。系统发育基因组分析显示，pSymA_IRAMC:0087 与 S. garamanticum LMG 24692^T 和 S. numidicum CIP 109850^T 的 pSymA 最近，二者均分离自突尼斯 Argyrolobium uniflorum 根瘤。Mauve 比对显示高度共线，repABC 操纵子系统发育将 IRAMC:0087 与上述两株菌聚为强支持分支，bootstrap = 100，并与 S. meliloti 1,021 明显分开。共生基因方面，pSymA 携带 nodD 和 syrM，nodD1 靠近 nodABCIJ，nodD2 靠近 nifHDKEX；宿主特异性结瘤基因 nodH、nodE、nodF、nodG、nodL、nodM、nodN、nolF、nolG 存在，nodQ 和 nfeD 位于 pSymB。固氮相关基因包括 nifH、nifD、nifK、nifE、nifN、nifQ、nifX、nifT、nifZ、nifB、4Fe–4S 铁氧还蛋白、nifS、nifU、nifA、fixAB 和 fixJLH；fixAB 有 3 个拷贝。nod/nol 和 nif 簇与 S. garamanticum、S. numidicum 共享保守多块基因顺序，而 S. meliloti 1,021 为单一连续 nif/fix 排列；IRAMC:0087 缺少 nodD3、noeA/noeB 和 ccoNOP。nodA 和 nodC 与 S. garamanticum LMG 24692^T 相似性分别为 98.82% 和 97.11%；nodC 系统发育将该菌株置于 symbiovar tunisiaense 分支，bootstrap = 100，与其独特宿主范围一致。

讨论部分总结：Sinorhizobium 属成员具有复杂多分复制子基因组，反映自由生活与共生生活方式适应。IRAMC:0087 的五个复制子与 S. meliloti 基因组高可塑性一致，附属质粒可能贡献环境适应、根际竞争力，有时参与共生功能。其撒哈拉来源与耐热、耐盐表型相符，基因组中热休克蛋白、分子伴侣、渗透保护物合成与摄取、胞外多糖和抗氧化基因支持这些表型。泛基因组高度开放，核心基因极少，云基因占绝大多数，符合 S. meliloti 由水平基因转移和复制子特异性进化塑造的塑性基因组特征。IRAMC:0087 的 PGPT 功能富集于群体感应、生物膜、分泌系统和表面附着，提示对自由生活持久和竞争性根定殖的投入增加；氮获取基因减少则与宿主范围常由可水平转移的离散共生基因决定而非物种分类一致。pSymA 结构分化、较小且更接近 S. garamanticum 和 S. numidicum，支持水平基因转移或共享质粒祖先。IRAMC:0087 虽归为 S. meliloti，但不能结瘤测试的 Medicago 物种，却有效结瘤 V. tortilis，且粉红根瘤和地上部干重增加支持有效固氮；其共生特性更符合 symbiovar tunisiaense。研究局限在于未能用 Calobota saharae 和 Argyrolobium uniflorum 活种子测试结瘤，且部分适应性性状由基因组推断，尚需实验验证。

研究结论翻译：理解根瘤菌对极端环境的适应，对于解析塑造共生效率、宿主特异性和生态适应性的进化机制至关重要。干旱生态系统提供了环境约束与宿主互作共同驱动微生物多样化和基因组可塑性的自然场景。本研究对撒哈拉 S. meliloti 菌株 IRAMC:0087 进行了综合表型、共生和基因组表征。其多分复制子基因组结构，以及共生质粒 pSymA 的显著结构分化，凸显该菌株对极端环境条件的进化灵活性。尽管系统发育基因组学将其归为 S. meliloti，IRAMC:0087 的宿主范围不同于经典 Medicago 相关生态型，而与 symbiovar tunisiaense 内与野生豆科植物相关的菌株一致。这些发现强调，决定宿主特异性和共生功能的是 symbiovar 身份，而非仅物种归属。除独特共生特征外，IRAMC:0087 结合了多种胁迫耐受机制和植物促生性状，包括耐热、耐盐和耐渗透胁迫，以及与养分动员和根际竞争力相关的性状。这种多功能组合可能为干旱土壤中的选择优势提供基础，使其既能作为自由生活土壤细菌持续存在，又能作为植物共生体发挥作用。总体而言，这些发现推进了对干旱生态系统中根瘤菌进化的理解，并说明基因组多样化，尤其是共生质粒和附属质粒内的多样化，如何促成 Sinorhizobium 属内的生态特化。IRAMC:0087 因此是研究沙漠适应根瘤菌的有价值模型，也是开发耐胁迫生物接种剂、改善生物固氮、减少合成氮输入并支持胁迫易发土壤可持续农业的有前景候选菌株。