大肠杆菌在亚致死浓度利福平胁迫下早期与长期适应的转录因子网络动态差异

《mSystems》：Distinct transcription factor network dynamics underlie early and long-term adaptation to rifampicin in Escherichia coli

　　

当大肠杆菌（Escherichia coli）暴露于亚致死浓度利福平（rifampicin）胁迫时，敏感细胞可快速启动适应过程以提高存活率。研究人员探究了转录因子网络（transcription factor network, TFN）是否参与该适应过程。首先，709个响应基因相对对照组出现偏离，提示强烈扰动；尽管这很大程度上是利福平的直接效应，但研究人员提供证据表明RNA聚合酶（RNA polymerase, RNAP）、σ70、σ38、Fnr及(p)ppGpp具有影响力，提示早期适应机制。进一步证据显示这些响应可依赖启动子序列，并与单细胞变异性的变化部分解耦；互作基因的响应呈相关性，提示协调调控。后期转录组与对照组部分重对齐，但272个基因表现与长期适应一致，其中若干基因的相关动态可由互作解释，形成分别由gcvB、cdaR和lldR控制的三个独立星状适应性基因模块。最后研究人员发现，进化距离遥远的病原体结核分枝杆菌（Mycobacterium tuberculosis）的同源基因对利福平的响应强度具相关性，提示上述现象可能跨菌种保守。总体而言，亚致死利福平对大肠杆菌全基因组转录水平的初始效应随时间被削弱，继而进入受基因网络影响的特异状态，可能促进适应；该发现有助于开发干扰细菌利福平适应的新策略。

研究背景与立项依据

在自然环境中，细菌常遭遇低于最小抑菌浓度（minimum inhibitory concentration, MIC）的亚致死抗生素（antibiotic, AB）胁迫，此类暴露可促进耐药机制涌现。利福平通过结合RNAP全酶、阻断启动子逃逸而抑制转录，其初始效应为全基因组mRNA迅速下降，但后续基因调控网络（gene regulatory network, GRN）如何驱动从早期扰动向长期适应态转变尚不清楚。既往研究多聚焦耐药突变或致死浓度下mRNA降解率测定，缺乏对转录因子网络动态贡献的时序解析。为此，研究人员以大肠杆菌K-12 MG1655为模型，在低于MIC（2.5 μg/mL，MIC=7.5 μg/mL）浓度下开展时序转录组与单细胞蛋白水平联合分析，旨在界定早期（80 min）与长期（160 min）适应相关的GRN驱动机制，并检验跨菌种保守性。该文发表于《mSystems》。

主要关键技术方法

研究以大肠杆菌K-12 MG1655野生型及系列报告株（YFP融合库、RpoC-GFP内源标记RL1314、Fnr-GFP单拷贝质粒库、染色体整合pFAB120合成启动子LacI位点突变库）为样本队列；在OD₆₀₀~0.3加入2.5 μg/mL利福平，于80 min与160 min采集样本，以无药对照校正生长时相差异。采用双时间点RNA-seq（Illumina NovaSeq 6000，DESeq2算log₂折叠变化（LFC））界定早期响应（709基因）与适应性基因（272基因，需在80与160 min间行为切换）；用流式细胞术定量单细胞YFP/GFP/mVenusNB水平；以半数诱导（250 μM IPTG）pFAB120文库检验启动子序列效应；利用RegulonDB v13.0互作关系与随机置乱零模型判定TF互作对LFC相关性的贡献；跨种比较借用NCBI GEO中结核分枝杆菌六株亚致死利福平RNA-seq（TPM），通过NCBI Gene正交基因比对评估响应强度相关性。

研究结果

预期利福平触发的GRN事件时序

基于利福平作用模式（mode of action, MoA），研究人员划分三阶段：直接RNA下降（~15 min内）、早期GRN驱动相（80 min）、晚期适应相（160 min）。选用80与160 min捕捉后两阶段。

亚致死利福平胁迫影响RNAP但不影响ATP、生长与形态

2.5 μg/mL不抑生长速率、ATP（QUEEN-2m探针）与细胞大小（脉冲宽度），但RpoB-YFP与RpoC-GFP单水平升高，RNA-seq亦佐证RNAP亚基转录上调，提示游离RNAP增加可部分抵消利福平阻塞。

直接/间接响应的时序特征

15 min LFC分布最广；互作基因输入输出LFC相关性在15 min最弱、80与160 min增强，且强于置乱零模型，表明GRN间接效应随时间放大。

操纵位点核苷酸改变扭转利福平效应

pFAB120文库中LacI位点单碱基突变使mVenusNB均值LFC与CV²变化解耦，证明启动子序列可微调利福平对逃逸步的阻碍程度。

亚致死利福平诱发转录组适应进程

80 min检出709早期响应基因（激活或抑制），160 min仅272基因表现行为切换（适应性基因），提示初期广扰后被压缩为特异重编程。YFP库验证astC、bhsA、hisJ、phoH、lldD等蛋白动态与RNA-seq趋势一致。

控制早期响应基因的调控机制

全局调控因子中仅fnr及其靶标一致响应；σ⁷⁰与σ³⁸靶标80 min普降、160 min恢复；spoT-YFP上调致(p)ppGpp水解增强，早期被(p)ppGpp激活基因多下调、被抑制基因不下调；早期互作基因输入输出LFC显著正相关（R²=0.33–0.38），置乱后消失，但早期→适应性基因间无相关。

控制适应性基因的调控机制

适应性基因互作仅160 min呈显著正相关，置乱零模型不成立；网络含34节点，分三独立星状模块，中枢为lldR（有氧乳酸利用）、cdaR（二酸糖摄取代谢TF）、gcvB（sRNA），伴8个自调控（6负2非正）与3孤立对，无三角环，提示有限独立程序与负自稳偏好。

大肠杆菌与结核分枝杆菌同源基因响应强度相关

567对正交基因在80（大肠杆菌）与60 min（结核分枝杆菌）LFC/TPM对数显著相关，虽弱于结核内株间相关，但支持利福平响应策略跨门级保守。

讨论与结论翻译

许多耐药株可能源于亚致死抗生素下的种群；亚致死浓度即可引发大规模转录重编程，且结核中TF可调节利福平耐受。本文揭示两阶段态：~700基因早期扰动，272基因长期切换（仅15%曾早期响应），构成阻尼初始效应并转入稳定适应态的基因组程序。早期阻尼由RNAP上调、σ⁷⁰/Fnr/(p)ppGpp局部重配及启动子序列贡献；早期互作基因相关但不延伸至适应性基因。晚期态由适应性基因内部互作维持，三星状模块（lldR/cdaR/gcvB）加负自稳，构成可靶向脆弱点。跨种同源基因响应相关，提示工程化抗适应策略或具广谱性。遗留问题包括TF是否依调控模改输出、枢纽删除株表型等。综上，亚致死利福平对大肠杆菌转录的初始冲击被时序削弱，继受GRN塑造为特异适应态，靶向上述TFN节点有望干扰细菌适应。