rmf作为实验演化Escherichia coli跨化学胁迫趋同休眠的保守核糖体休眠标志物

《Microbiology Spectrum》：rmf, a conserved ribosome hibernation marker of convergent dormancy in experimentally evolved Escherichia coli

　　

摘要翻译：细菌能够适应极其多样的化学胁迫，但独立适应轨迹能在多大程度上趋同、以及趋同于何种终点，仍不清楚。研究人员重新分析了188个在47种不同化学胁迫下演化的大肠杆菌（Escherichia coli）谱系、4个M9对照和4个祖先重复的转录组。研究人员利用多个公共数据集检验趋异的谱系是否趋同于共享的生理终点。在化学性质多样的胁迫下演化的谱系趋同于一个共享转录轴，其表达变化在不同化学物之间呈正相关（r = 0.157；P = 0.001）。该轴对应于一个慢生长、由RpoS（稳定期σ因子）调控的稳定期与核糖体休眠程序，与RpoS/稳定期指数高度一致（r = 0.97），并与核糖体指数负相关（r = ?0.65）。核糖体调制因子（ribosome modulation factor, RMF）编码基因rmf可形成翻译沉默的100S核糖体二聚体，其表达沿PC1（第一主成分）轴升高（r = 0.65），并与100S休眠因子hpf、yfiA和sra特异性共表达。相比之下，其与持留菌毒素-抗毒素（toxin-antitoxin, TA）、多药外排和氧化应激基因的关联较弱或不一致。跨物种分析支持RMF受纯化选择（purifying selection），全基因组分析将rmf置于核心基因多样性最低的4%范围内，与通过调控而非编码改变实现趋同相一致。该轴和rmf表达均未追踪谱系的耐药水平，表明该程序是保守的适应状态而非耐药机制。这些发现将经典应激反应程序重新界定为适应性演化的趋同终点，并将rmf鉴定为该状态的保守读出标志。

论文解读：研究背景方面，细菌在抗生素及其他环境胁迫下可通过耐药、耐受和持留等策略存活；耐药通常表现为最低抑菌浓度（minimum inhibitory concentration, MIC）升高，耐受可在不升高MIC的情况下延长杀菌暴露中的存活，持留则是仅亚群存活的异质性耐受形式，休眠或生长停滞可贡献于耐受和持留但不等同于二者。尚待回答的演化问题是：当细菌适应化学和机制多样的胁迫时，是沿特异的胁迫路线演化，还是反复趋同于共享生理状态。适应性实验室演化（adaptive laboratory evolution, ALE）结合表达谱提示趋同真实存在，但其机制解释未定；既有研究多从耐药视角解读共享轴，尚不清楚该轴是否反映翻译下调或核糖体休眠、是否存在单一可解释基因标志、以及该标志是否受进化约束。核糖体调制因子（ribosome modulation factor, RMF）是55个氨基酸蛋白，在稳定期E. coli中促成100S核糖体形成，由σ⁷⁰依赖启动子转录并受(p)ppGpp调控，可与休眠促进因子（hibernation-promoting factor, HPF）共同连接两个70S核糖体形成翻译失活的100S二聚体。rmf突变株在丰富培养基中生长正常，但在饥饿、抗生素暴露和冷胁迫中存活力下降100至1,000倍。尽管如此，rmf跨化学无关胁迫的全基因组意义及其作为适应靶点、水平转移或与耐药基因共演化的情况仍未系统评估。研究人员因此重新分析多胁迫ALE转录组汇编，并结合独立公共数据集和细菌rmf演化记录，检验跨机制无关化学选择的转录组趋同、rmf是否标志主要趋同轴及其所属细胞程序、跨数据集可重复性和rmf的进化稳定性。总体结论是：化学多样选择在转录组水平趋同于以核糖体休眠和稳定期重塑为主的低维状态，rmf是该状态最可解释且高度受限的标志之一，并非化学特异性耐药决定因素；其意义在于把ALE趋同与持留菌趋同研究联系起来，将经典应激反应程序重新界定为适应性演化的趋同终点，并为耐受相关休眠状态提供保守读出候选，但其因果作用和治疗靶点价值仍需直接遗传与药理学验证。

关键技术与数据来源方面，研究使用GSE137348中196个转录组谱，包括188个在47种化学选择条件下独立演化的MDS42（无插入序列、基因组缩减的E. coli K-12衍生株）终点克隆、4个M9演化对照和4个祖先MDS42重复；另纳入18个手工整理的独立公共E. coli转录组数据集共149个样本，涵盖实验室、环境和临床或致病株；使用26个RMF同源物和200个完整E. coli基因组。主要方法包括主成分分析（principal component analysis, PCA）、全基因组表达变化向量的成对Pearson相关、1,000次基因标签置换检验、剔除rmf及核糖体蛋白和RpoS程序基因后的留出PCA、控制核糖体蛋白指数的偏相关；跨数据集z评分和残差化、GSE分层bootstrap、Fisher随机效应meta分析、方向性基因本体（Gene Ontology, GO）富集（Fisher精确检验、Benjamini-Hochberg校正、标准化富集分数）；MAFFT蛋白比对和回译密码子比对、PAML codeml自由ω与ω=1中性模型比较、似然比检验（likelihood-ratio test, LRT）和剖面似然置信区间（confidence interval, CI）；严格单拷贝核心基因过滤及核苷酸多样性（π）计算；NCBI Nucleotide记录人工审阅和去重。上述方法概括不超过250字，未涉及试剂和具体培养、质粒构建步骤。

研究结果部分保留原文小标题如下。Chemically diverse selection converges on a shared, low-dimensional transcriptomic state：通过PCA和成对相关分析，196个转录组未形成按药物靶点区分的离散簇，而是投射到单一共享轴；PC1解释21.7%总转录方差，PC2解释6.9%，不同无关靶点类别的PC1分布广泛重叠。同化学选择的重复谱系平均相关（mean r = 0.404；282对）高于异化学谱系（mean r = 0.212；17,296对），化学选择与M9演化谱系间也为正相关（mean r = 0.264；752对），M9对照间平均相关更高（mean r = 0.587），提示共享分批培养传代贡献显著。以4个M9对照平均谱为参照后，异化学选择谱系间相关仍为正（mean r = 0.157；P = 0.001），回归生长代理后异化学平均相关仍为正（r = 0.17），说明传代适应贡献趋同但不能完全解释，生长本身也不能单独解释。结论是化学多样选择诱导趋同的低维转录适应景观。

rmf marks the convergent hibernation state but not the resistance level of evolved lineages：通过将基因或模块表达与PC1关联，rmf表达沿PC1单调升高（Pearson r = 0.65）；核糖体蛋白指数与PC1强负相关（r = ?0.65），RpoS/稳定期指数与PC1几乎完全一致（r = 0.97），表明沿轴推进对应翻译下调和RpoS调控的稳定期程序。rmf与100S休眠因子hpf、yfiA和sra强而选择性共表达，与持留菌TA模块、多药外排基因和氧化应激基因基本无关联；按PC1排序后，rmf位于协同上调的休眠和稳定期模块中，核糖体蛋白基因构成相互下调模块。留出轴重算与原轴一致（agreement r = 1.00），rmf与该独立轴相关未减弱；控制生长指数后rmf仍追踪轴（partial r = 0.41，P = 3×10^?9）。rmf表达与谱系达到的相对IC₅₀（半数抑制浓度）基本无关（Pearson r = 0.11，P = 0.124），其他耐药指标均未超过弱关联（all |r| ≤ 0.17）。结论是rmf标志连接RpoS的趋同核糖体休眠程序，而非耐药水平。

The rmf hibernation program generalizes across independent E. coli data sets：通过18个独立公共转录组数据集149个样本的汇总相关和方向性GO富集，rmf与100S休眠因子hpf、yfiA及RpoS/稳定期程序基因dksA、rpoS、dps和sra稳健共表达；与多药外排泵、抗菌耐药调节因子和持留菌TA模块的关联弱、不稳定或为负。GO富集最显著术语的标准化富集分数为负，即富集于与rmf负相关的基因；翻译、核糖体结构与生物发生、rRNA和tRNA加工、核苷酸生物合成、氧化磷酸化和ATP合成以及质子转运均随rmf升高而协同下调，正相关基因中无类别富集。结论是该休眠程序在独立背景和临床株中可推广，对应代谢和翻译协同关闭，而非主动药物防御。

rmf is under purifying selection across species and lies within the low-diversity tail of the E. coli core genome：通过26个物种RMF比对和PAML codeml比较，24个可估计比较的ω（d_N/d_S）点估计均低于1，22个比较的95%剖面似然CI完全低于1，23/24个可估计比较的单侧LRT拒绝ω=1中性模型（P < 0.05）；Shigella flexneri 2a 301比较因无同义替换不可估计。通过200个完整E. coli基因组中1,519个严格单拷贝核心基因计算平均π，rmf的π为0.001193，按升序排名第61，位于严格核心分布多样性最低的4.0%内；预设小核糖体蛋白基因rpsU和rpmJ不变，rpmB多样性略低于rmf，休眠因子hpf和yfiA多样性略高。NCBI Nucleotide检索人工审阅去重后保留150条含明确rmf编码序列的记录，其中137条注释为染色体、13条注释为质粒，分别占91.3%和8.7%；这些计数描述NCBI来源注释，不等同于独特菌株、独特复制子或独立验证的质粒定位。结论是RMF受广泛纯化选择，rmf处于E. coli核心基因组低多样性尾，趋同通过调控而非编码改变实现。

讨论部分总结：讨论指出，化学和机制多样选择下的E. coli实验演化在转录组水平趋同于由核糖体休眠和稳定期重塑主导的共享低维状态，该状态由RpoS/稳定期模块读出，rmf是最可解释且高度受限的成员；同物种功能基因组学和跨物种群体基因组学两条证据支持rmf不是化学特异性耐药决定因素，而是100S核糖体休眠程序的保守组分，在持续毒性胁迫中激活。Maeda et al.的框架针对耐药而非核糖体生理、生长停滞或休眠，本研究无监督分析揭示其耐药监督投影未强调的互补主导轴；Fong et al.发现表型相似演化E. coli可通过不同表达状态到达共同终点，而本研究表明表达水平趋同，提示多样化学毒性下适应空间窄到足以将无关谱系导向共同调控终点。RMF与HPF促进翻译失活100S核糖体形成，与生长停滞和抗生素耐受相关而非耐药；rmf表达与定量耐药基本无关。跨物种codeml多支持ω<1，全基因组分析将rmf置于最低多样性4.0%核心基因内，但并非唯一不变，若干小核糖体蛋白同等或更保守，hpf和yfiA多样性略高，故rmf属于强约束核糖体相关基因群而非例外位点。Ekemezie and Melnikov提出阻断休眠因子结合可能促进休眠细胞保护性核糖体丢失，但休眠因子在不同细菌中不统一，RMF主要限于部分γ-变形菌，许多其他细菌使用长HPF，因此休眠过程保守不等于靶向RMF化合物具有广谱活性；已知抗生素结合位点重叠的因子是HPF、RaiA和Balon，而RMF结合位点距这些位点超过10 ?。研究支持rmf主要作为E. coli休眠状态的保守标志和功能测试候选，而非已验证直接药物靶点；若未来遗传和药理学研究显示破坏RMF依赖休眠增加抗生素敏感性，最合理应用是常规抗生素辅助而非单药治疗。作为标志，模块可互换，任何共调控成员（hpf、yfiA、sra和dps）均可报告该状态；突出rmf的原则在于其是机制上非冗余、特异且不可替代的100S形成启动者。局限包括趋同真实但部分（PC1解释21.7%方差，跨谱系平均相关中等），发现实验仅用无IS的MDS42，未捕获水平基因转移；标志非机制，观察性分析不能证明RMF介导休眠导致耐受状态；跨物种分析受短rmf序列和高d_S限制；核苷酸多样性分析限于E. coli，未扩展到Klebsiella和Enterobacter等其他ESKAPE-E γ-变形菌；共表达推断依赖富集胁迫和稳定期条件的公共转录组，可能高估模块一致性；rmf单独是否最佳标志或紧凑休眠模块签名更优仍待检验。研究结论翻译：这些发现将经典应激反应程序重新界定为适应性演化的趋同终点，并将rmf鉴定为该状态的保守读出；因此，rmf应被解释为E. coli核糖体休眠程序的保守标志物，而非可遗传耐药的决定因素或已验证的治疗靶点。