Transcriptome Profiling Reveals Interplay of Multifaceted Stress Response in Escherichia coli on Exposure to Glutathione and Ciprofloxacin
转录组分析揭示大肠杆菌暴露于谷胱甘肽和环丙沙星时多层面应激反应的相互作用
来源：msystems.asm.org January/February 2018 Volume 3 Issue 1 e00001-18
1.论文摘要核心内容
本研究团队此前已报道外源性谷胱甘肽（GSH）可通过中和抗生素诱导的氧化应激并增强抗生素外排，促进大肠杆菌对环丙沙星的耐药性。本研究利用全基因组微阵列技术，系统分析了大肠杆菌在单独暴露于GSH、亚抑菌浓度环丙沙星以及两者联合处理下的转录组变化。微阵列数据显示，GSH补充显著影响了大肠杆菌的氧化还原功能、物质转运、酸休克和毒力相关基因的表达；环丙沙星则主要诱导DNA损伤修复和细胞分裂相关基因的表达。两者联合作用时，多个潜在的应激反应通路（包括上述基因及DNA损伤修复基因）发生相互作用，共同抵消环丙沙星的抗菌效应。研究通过逆转录定量PCR（RT-qPCR）对40个差异表达基因进行了大规模验证，结果与微阵列数据基本一致。进一步通过12个单基因缺失株的生长表型分析，证实了这些基因不仅在分子水平，更在功能水平参与了GSH介导的表型变化。基于转录组数据，研究还首次发现并验证了GSH补充能够显著提高大肠杆菌的酸休克存活率。综上，GSH除了在生理水平发挥作用外，还通过调控多个应激反应通路的基因表达，共同对抗环丙沙星的抗菌作用。
2.中文关键词（单行）
DNA损伤、酸抗性、抗生素抗性、细胞氧化还原状态、氧化还原功能、应激反应、转运蛋白、毒力
3.研究目的
在系统水平解析外源性GSH单独及与环丙沙星联合作用下大肠杆菌的全基因组表达谱，全面揭示GSH对细菌转录调控的影响。
鉴定GSH介导环丙沙星耐药性的关键基因和调控通路，补充此前仅发现的氧化应激中和与抗生素外排两种机制，构建完整的耐药性分子网络。
通过大规模基因缺失株的功能验证，建立差异基因表达变化与细菌生理表型之间的因果关系，区分功能性调控与非特异性表达波动。
探索GSH对大肠杆菌其他应激反应（如酸胁迫）的潜在影响，拓展对GSH在细菌应激适应中生物学功能的认知。
为理解细菌复合胁迫响应的整合机制提供新的视角，为开发克服抗生素耐药性的新策略提供理论基础。
4.研究思路
本研究采用"转录组学全局分析-分子水平验证-功能表型验证-新表型拓展"的系统性研究策略：
第一步，设置四个实验组：空白对照组（LB培养基）、10 mM GSH处理组、3 ng/ml亚抑菌浓度环丙沙星处理组、10 mM GSH+50 ng/ml环丙沙星联合处理组。培养大肠杆菌MG1655至对数中期，提取总RNA进行Affymetrix大肠杆菌全基因组微阵列分析，筛选差异表达基因并进行基因本体（GO）功能分类。
第二步，对所有表达变化≥5倍的41个基因（除ffs基因因片段过小无法设计引物外）进行RT-qPCR验证，确认微阵列数据的可靠性和准确性。
第三步，从Keio单基因缺失文库中选取36个可获得的缺失株，利用Bioscreen C全自动生长曲线分析仪，测定其在上述四种处理条件下的生长动力学，鉴定对GSH和/或环丙沙星表型有显著影响的功能基因。
第四步，基于转录组数据中酸休克基因显著上调的发现，设计酸抗性实验，测定大肠杆菌在不同处理条件下经pH 3.0酸休克1小时后的存活率，验证GSH对酸胁迫适应的影响。
第五步，整合所有实验数据，提出GSH通过多通路协同作用介导环丙沙星耐药性的综合机制模型。
5.研究亮点
首次报道了外源性GSH处理下大肠杆菌的全基因组转录组变化，系统揭示了GSH除经典抗氧化功能外的多重调控作用，包括酸抗性调控、物质转运重编程和毒力基因表达调节。
证明GSH介导的环丙沙星耐药性是一个多层面的系统过程，而非单一通路作用。除此前报道的氧化应激中和与TolC-AcrAB外排泵增强外，还涉及SOS DNA损伤反应抑制、酸应激通路激活和毒力网络重塑等多个机制。
通过大规模的基因缺失株表型筛选，明确了12个关键基因在GSH介导的生理表型中的功能作用，其中多个基因（如asr、yeiH、srlE）与GSH的关联此前从未被报道。
首次实验证实外源性GSH能够显著增强大肠杆菌在强酸性条件下的存活率，拓展了对GSH在细菌环境适应中作用的认知，为理解肠道致病菌的宿主适应机制提供了新线索。
发现超过一半的GSH调控蛋白定位于周质或细胞膜组分，提示细菌的膜系统和周质空间是GSH发挥作用的关键亚细胞位点。
6.可延伸的方向
深入解析GSH调控酸抗性的分子机制，明确其是通过直接调节胞内pH、激活谷氨酸依赖型酸抗性系统，还是通过代谢产生硫化氢等信号分子发挥作用。
研究GSH对其他类别抗生素（如β-内酰胺类、氨基糖苷类、四环素类）抗菌活性的影响，验证其是否具有介导广谱抗生素耐药性的能力。
利用ChIP-seq技术鉴定GSH响应的核心转录因子（如SoxR、OxyR、EvgA、RpoS）的全基因组结合位点，构建GSH介导的转录调控网络。
建立小鼠肠道感染模型，研究GSH在体内环境中对大肠杆菌致病性和抗生素治疗效果的影响，评估其临床相关性。
分析临床分离的多重耐药大肠杆菌菌株中GSH代谢相关基因的表达水平，探讨GSH调控通路在临床耐药性发展中的作用。
开发靶向GSH生物合成或其下游调控通路的小分子抑制剂，与环丙沙星等抗生素联合使用，探索逆转细菌耐药性的新策略。
研究GSH与其他宿主免疫因子（如活性氧、抗菌肽、补体系统）的相互作用，解析其在宿主-病原体互作中的复杂角色。
7.测量的数据、对应图表及研究意义
GSH介导环丙沙星耐药性的机制示意图，来自Figure 1。研究意义：直观总结了本研究团队此前发现的两种核心机制（中和抗生素诱导的氧化应激、增强TolC-AcrAB外排泵活性），为后续转录组研究提供了清晰的背景和研究基础。
全基因组微阵列表达数据，包括四种处理条件的基因表达热图、差异表达基因倍数分布饼图、不同处理组差异基因韦恩图，来自Figure 2。研究意义：全面展示了GSH、环丙沙星及联合处理对大肠杆菌转录组的全局影响，鉴定出609个表达变化≥2倍的基因，揭示了不同处理条件下基因表达模式的异同。
40个差异表达基因的RT-qPCR验证数据，来自Figure 3。研究意义：从分子水平验证了微阵列数据的可靠性，确认了大部分基因的表达变化趋势，为后续的功能研究提供了准确的分子依据。
野生型BW25113及5个关键基因缺失株（recA、yeiH、asr、srlE、manY）在四种处理条件下的生长曲线，来自Figure 4。研究意义：直观展示了基因缺失对细菌在GSH和/或环丙沙星胁迫下生长的影响，直接证明了这些基因在GSH介导的生理表型中的功能重要性。
大肠杆菌在不同处理条件下的酸休克存活率柱状图，来自Figure 5。研究意义：首次通过实验证明外源性GSH能够将大肠杆菌在pH 3.0条件下的存活率提高约2倍，验证了转录组数据中酸休克基因上调的生理意义。
41个表达变化≥5倍基因的微阵列表达数据及亚细胞定位信息，来自Table 1。研究意义：系统列出了最显著差异表达的基因及其功能和亚细胞定位，发现超过50%的GSH调控蛋白定位于周质或细胞膜，提示这些区域是GSH作用的关键位点。
用于RT-qPCR验证的40个基因的引物序列、功能及NCBI基因ID，来自Table 2。研究意义：提供了实验验证的详细方法学信息，确保了实验的可重复性和结果的可靠性。
本研究使用的所有大肠杆菌菌株列表及基因型，来自Table 3。研究意义：明确了实验菌株的来源和遗传背景，为实验结果的解读和比较提供了基础。
其余7个功能重要基因缺失株的生长曲线数据，来自Figure S1。研究意义：补充了Figure 4的结果，全面展示了12个关键基因缺失株在不同处理条件下的生长表型。
8.研究结论
外源性GSH处理显著重塑了大肠杆菌的转录组，影响了氧化还原平衡、物质转运、酸应激响应和毒力等多个生物学过程，其中膜系统和周质空间是GSH调控的主要亚细胞区域。
亚抑菌浓度的环丙沙星主要诱导DNA损伤修复和细胞分裂相关基因的表达，激活经典的SOS应激反应，这与氟喹诺酮类抗生素通过抑制DNA拓扑异构酶导致DNA断裂的作用机制完全一致。
GSH与环丙沙星联合处理时，基因表达模式主要与GSH单独处理相似，同时GSH能够部分逆转环丙沙星诱导的SOS反应，表明GSH通过多重通路协同抵消环丙沙星的抗菌效应。
12个单基因缺失株（asr、fecA、manX、manY、ogrK、recA、recN、srlB、srlD、srlE、yeiH、yjiY）在GSH和/或环丙沙星处理下表现出显著的生长缺陷，证明了这些基因是GSH介导的细菌应激适应所必需的。
外源性GSH能够显著增强大肠杆菌在强酸性条件下的存活率，这与转录组数据中酸休克核心基因（asr、frc、oxc、ydeO、yegR）的显著上调直接相关，揭示了GSH在细菌酸胁迫适应中的新功能。
GSH介导的环丙沙星耐药性是一个整合的系统响应，涉及氧化应激中和、抗生素外排增强、DNA损伤修复抑制、酸应激通路激活和毒力基因表达改变等多个机制的协同作用。
9.芬兰Bioscreen仪器测量的微生物生长曲线数据的研究意义详细解读
本研究使用芬兰Oy Growth Curves Ab Ltd.生产的Bioscreen C全自动微生物生长曲线分析仪，测定了野生型大肠杆菌BW25113及36个单基因缺失株在四种处理条件（LB对照、10 mM GSH、3 ng/ml环丙沙星、10 mM GSH+50 ng/ml环丙沙星）下的生长动力学。仪器设置为37℃恒温培养，连续振荡，每30分钟自动读取一次600nm光密度值（OD₆₀₀），连续监测15小时以上。所有生长曲线均接种至初始OD₆₀₀为0.01，培养体积为100μl，每个实验条件设置3个生物学重复。该仪器产生的生长曲线数据在本研究中具有以下六个核心研究意义：
第一，实现高通量功能基因筛选，高效建立分子变化与生理表型的关联。
转录组分析鉴定出数百个差异表达基因，传统的试管培养方法无法在短时间内完成如此多基因的功能验证。Bioscreen的100孔板设计使得可以同时测定多个菌株在多种条件下的生长曲线，本研究在短短一周内完成了超过144组（36株×4条件）独立的生长实验，筛选效率提高了数十倍。通过比较野生型和缺失株的生长曲线，研究人员快速从36个候选基因中鉴定出12个功能重要基因，直接将转录组的分子变化与细菌的生理表型联系起来，大大加速了研究进程。
第二，精确量化生长参数差异，客观评估基因缺失的表型效应强度。
传统的终点OD测量只能获得单一时间点的生物量数据，无法准确反映比生长速率、延迟期长度和平台期生物量等关键生长参数的细微差异。Bioscreen的连续监测能够生成完整的生长曲线，通过数学拟合精确计算每个菌株在不同条件下的生长参数。例如，数据显示recA缺失株在任何含有环丙沙星的培养基中都完全不能生长，而yeiH缺失株在GSH条件下的比生长速率降低了约40%，最大生物量降低了约50%。这种精确的定量差异，客观地评估了每个基因缺失对细菌适应不同胁迫条件的影响程度，为后续机制研究提供了量化依据。
第三，提供高时间分辨率的生长动力学数据，揭示基因功能的时间特异性。
不同基因可能在细菌生长的不同阶段发挥作用，Bioscreen每30分钟一次的高密度采样能够完整捕捉整个生长过程的动态变化。例如，manY缺失株在环丙沙星单独处理下的延迟期比野生型延长了约2小时，但对数生长期的生长速率与野生型基本一致；而在GSH+环丙沙星条件下，其整个生长周期的生长速率都显著降低。这种生长阶段特异性的表型差异，提示manY可能在细菌应对环丙沙星胁迫的早期适应阶段发挥重要作用，为深入理解基因的功能机制提供了重要线索。
第四，验证转录组数据的生理相关性，区分功能性调控与非特异性表达变化。
转录组分析会检测到大量差异表达基因，但其中很多变化可能是细胞的被动响应，而非具有生理功能的主动调控。Bioscreen的生长曲线数据能够验证这些基因表达变化是否具有实际的生理意义。例如，转录组数据显示srlB、srlD、srlE等山梨醇转运和代谢基因在GSH处理下显著下调，而对应的基因缺失株在GSH条件下的生长确实受到显著抑制，证明了这些基因的表达变化是GSH适应所必需的。相反，一些表达变化较小（2-3倍）的基因（如manY）也表现出明显的生长表型，提示即使是低水平的表达变化也可能具有重要的生理功能。
第五，标准化的生长数据确保实验的可重复性和可比性。
Bioscreen C是微生物学研究中广泛使用的标准生长曲线测定仪器，其标准化的操作流程和数据输出格式，使得本研究的生长数据能够与其他实验室的大肠杆菌生长研究结果进行直接比较。例如，本研究中野生型BW25113在LB培养基中的比生长速率约为0.8 h⁻¹，与已发表的文献数据高度一致，验证了实验体系的可靠性。同时，标准化的数据也为后续构建大肠杆菌的基因组规模代谢模型和应激响应模型提供了准确的生理参数。
第六，为后续机制研究提供明确的靶点和方向。
Bioscreen筛选出的12个功能重要基因，为后续深入研究GSH介导环丙沙星耐药性的分子机制提供了明确的靶点。例如，recA和recN是SOS DNA损伤反应的核心基因，它们的缺失株对环丙沙星极度敏感，且GSH无法恢复其生长，直接证明了GSH对SOS反应的抑制是其耐药机制的重要组成部分。asr是酸休克的关键诱导基因，其缺失株在GSH条件下生长缺陷，进一步支持了GSH通过激活酸应激通路增强细菌存活的结论。这些发现为后续的机制研究指明了方向，大大提高了研究效率。