银胁迫下大肠杆菌全基因组DNA甲基化组特征解析

《Microorganisms》：Characterization of the Global DNA Methylome in Escherichia coli Under Silver Stress

　　

细菌持续暴露于各种胁迫因子，其感知、适应和应对环境变化的能力对生存适应性至关重要。双组分响应系统（two-component response systems，TCRS）在介导这些响应中发挥主要作用。然而，表观遗传修饰在适应和驯化中的作用尚不清楚。研究人员利用单分子实时测序（Single-Molecule Real-Time sequencing，SMRT）技术，表征了野生型大肠杆菌（Escherichia coli）K-12 MG1655、四个cisS突变株及六个银适应突变株在有无50 ng/mL AgNO3暴露1小时条件下的N6-甲基腺嘌呤（N6-methyladenine，6mA）和N4-甲基胞嘧啶（N4-methylcytosine，4mC）DNA甲基化模式。选择cisS突变株旨在研究突变对DNA甲基化的影响。所有菌株在有无银暴露条件下，全局平均6mA比例均较高，且主要由各菌株中占主导地位的GATC基序以及野生型和银适应突变株中的AACN6GTGC/GCACN6GTT基序驱动；而全局平均推定4mC比例则较低。在不同条件下，中等至低甲基化水平的扩展简并基序存在变异。研究结果表明，基序特异性甲基化模式随银暴露和遗传背景的不同而变化，而大肠杆菌的全局甲基化水平则表现出显著的韧性。本研究加深了对微生物胁迫适应与胁迫响应的理解，提示DNA甲基化可能参与TCRS介导的胁迫耐受。

论文解读

一、研究背景与科学问题

细菌在自然环境中持续面临多种胁迫因子，其感知、适应和应对环境变化的能力直接关系到生存适应性。双组分响应系统（two-component response systems，TCRS）是细菌介导胁迫响应的核心机制，由感知外界刺激的组氨酸激酶和接受磷酸基团并启动转录响应的响应调节蛋白组成。在大肠杆菌（Escherichia coli）中，银适应与抗性最稳定地与CusS/R双组分系统相关联：CusS感知Ag⁺后发生自磷酸化，将磷酸基团转移至CusR，进而激活cisRS和cisCFBA基因的转录，介导Ag⁺外排。银作为一种广谱抗菌剂已被使用数百年，但细菌对银的耐药性已有报道，且实验进化研究显示银抗性可快速产生，使其成为研究细菌胁迫、适应与驯化响应的理想抗菌环境。

尽管TCRS介导的转录响应已被广泛研究，但DNA甲基化在细菌胁迫响应中的作用，以及其对应变特异性差异或适应相关差异的贡献，尤其是与TCRS的关联，仍不清楚。DNA甲基化是一种表观遗传修饰，通过向特定核苷酸碱基共价添加甲基基团，在不改变DNA序列的前提下影响基因组功能。细菌中三种主要DNA甲基化形式为N6-甲基腺嘌呤（6mA）、N4-甲基胞嘧啶（4mC）和5-甲基胞嘧啶（5mC）。已有研究表明，细菌在胁迫下可发生局部或基序特异性的甲基化可塑性变化，而非全基因组甲基化组的破坏。例如，大肠杆菌在营养限制下、鼠伤寒沙门氏菌（Salmonella enterica Typhimurium）在氧化胁迫下以及幽门螺杆菌（Helicobacter pylori）在酸胁迫下均观察到局部甲基化变化，而全基因组甲基化水平基本保持稳定。

基于上述背景，本研究旨在利用单分子实时测序（SMRT）技术，系统表征野生型大肠杆菌K-12 MG1655、四个cisS突变株和六个银适应突变株在低水平银胁迫条件下的全基因组DNA甲基化图谱，重点识别主导序列基序、基序分布变化及跨菌株和基因型的比较甲基化模式。低水平银暴露（50 ng/mL AgNO₃，1小时）被选为亚致死胁迫条件，在维持细菌生长的同时评估甲基化组在非致死暴露条件下的响应。

二、研究设计与技术方法

研究人员使用野生型大肠杆菌K-12 MG1655（ATCC #10798D-5）、四个cisS突变株（R15L、T14P、T17P和N279H）和六个银适应突变株（SAM1至SAM6）。cisS突变株通过pORTMAGE重组工程在大肠杆菌K-12 MG1655中引入cisS单核苷酸替换构建；银适应突变株通过硝酸银实验进化产生。所有菌株在LB培养基中复苏后，经Davis Minimal Broth（DMB）三步驯化培养，在OD₆₀₀=0.2（对数中期）时进行银暴露处理，共收集66个沉淀样本（3个生物学重复×11个菌株×2个处理条件）。

主要技术方法包括：PacBio Revio平台的SMRT测序（HiFi化学，平均覆盖度>400×）；基因组组装采用Flye 2.9.2并利用Circlator 1.5.5环化；甲基化检测使用ipdSummary 3.0算法，基于聚合酶动力学信号（interpulse duration ratios）识别6mA和4mC修饰位点；统计比较采用卡方检验、Mann–Whitney U检验和Kruskal–Wallis检验（Dunn事后多重比较校正），显著性阈值p<0.05。

三、研究结果

3.1 测序输出、甲基化组数据质量与甲基化图谱概览

基因组大小约4.64 Mb，符合大肠杆菌预期大小；测序深度>400×平均覆盖度，% bp ≥ Q30至少97%。检测到的基序仅在平均IPD比值>1时才用于分析和数据解读。除6mA和推定4mC外，还检测到额外的聚合酶动力学信号，但这些信号的化学身份无法通过标准SMRT甲基化流程确定，因此未作进一步推断。SAM6（银处理，重复1）组装为两个重叠群；SAM4（无银，重复2）组装长度正常，而其余五个SAM4组装因缺乏e14区域而短约15.2 kb。

3.2 跨菌株或条件的比较甲基化模式

全局基因组甲基化分析显示，6mA是数据集中检测到的主要修饰类别。卡方检验显示，所有菌株的合并基序计数与对照组（无银野生型）存在显著差异（p<0.05），但效应量（Cramer's V）在cisS突变株中较弱（0.13–0.14），在SAMs中可忽略（<0.1），表明基序分布差异可检测但幅度微弱。cisS突变株效应量较大，反映其与野生型对照的基序组成差异更大，特别是AACN₆GTGC/GCACN₆GTT基序的完全缺失，而该基序在SAMs中仍可检测。4mC基序的卡方检验显示选择性显著变化，在银适应突变株中比cisS突变株更明显，但效应量均为可忽略水平。

Mann–Whitney U检验显示，所有菌株在有无银条件下全局甲基化水平均稳定，与对照相比无显著差异（p>0.05）。全局加权平均6mA甲基化比例范围为约0.69至接近1.00；全局加权平均推定4mC甲基化比例范围为约0.08至接近0.09。Kruskal–Wallis检验显示全局6mA甲基化比例存在显著菌株效应（p<0.05），Dunn事后检验鉴定SAM1在银暴露下与未处理野生型对照显著不同（校正p=0.0344），但对应的Mann–Whitney U比较不显著（p=0.1）。全局推定4mC也显示组间显著总体差异（p<0.05），但无菌株在Dunn事后检验中与对照显著不同。

在无银条件下，SAM3的一个重复偏离主导6mA全局甲基化模式，经敏感性分析确定该偏差主要由低比例基序GACYNVNBND（甲基化比例约0.06）驱动，而非全基因组甲基化丢失。

6mA基序中，GATC和AACN₆GTGC/GCACN₆GTT基序在多个菌株中显示完全或接近完全甲基化。GATC基序在所有菌株和重复中均检测到，不受处理影响；AACN₆GTGC/GCACN₆GTT基序在野生型和银适应群体的所有重复中均存在，但在所有cisS突变株中均缺失，归因于所有cisS重组克隆中共享的hsdS单碱基对缺失。银暴露后，野生型和大多数银适应群体中大多数简并基序信号未被检测到，仅SAM3和SAM6中检测到ACYKGATCW和VAGGKGATC基序。

推定4mC基序包括CNK、CCAGGBAD、CCAGGYAD和CGAGGAYCD。无银条件下，CNK基序存在于所有菌株和重复中，但甲基化比例低（约0.08–0.09）；CCAGGBAD在所有菌株中检测到，但未在所有重复中普遍检测到；CCAGGYAD仅在CCAGGBAD未检测到的重复中出现。银暴露后，CNK基序保持稳定，CCAGGYAD未被检测到，而CCAGGBAD在所有菌株和重复中一致存在。CGAGGAYCD仅在SAM3的一个银处理重复中检测到，敏感性分析显示其不影响整体结果。

总体而言，研究结果支持基序水平甲基化再分布与全局甲基化比例相对稳定并存。卡方检验敏感于甲基化基序计数的分布变化，Mann–Whitney和Kruskal–Wallis检验评估重复水平甲基化比例是否在组间存在差异。热图显示处理前甲基化基序谱更异质，处理后更同质，主导甲基化系统未改变。

四、讨论总结与结论翻译

讨论总结：本研究提供了大肠杆菌K-12 MG1655的全基因组DNA甲基化表征。6mA是最主要的甲基化形式，GATC基序几乎完全甲基化（甲基化比例>0.99但<1），在所有菌株和条件下保守，反映Dam甲基转移酶的高活性，与文献报道的GATC基序在多种胁迫条件下的高度稳定性一致。AACN₆GTGC/GCACN₆GTT基序在野生型和所有SAMs中保守且接近饱和甲基化（0.99–1），是EcoKI相关I型限制修饰系统（restriction-modification system，R-M系统）的经典特征。cisS突变株中该基序的完全缺失可归因于pORTMAGE工程中hsdS基因的脱靶单碱基对缺失，hsdS编码DNA识别亚基，其完整性对EcoKI系统的功能至关重要。SAM4中五个组装缺失e14相关前噬菌体区域（约15 kb），而一个重复保留该区域，提示原始甘油菌株库中存在异质性，银适应可能促成或贡献了e14切除。推定4mC基序（CNK、CCAGGBAD、CCAGGYAD）在文献和REBASE数据库中缺乏报道，其中CCAGGBAD和CCAGGYAD含CCAGG核心，对应Dcm识别序列CCWGG的W=A变体，其4mC分类应视为推定而非确定性证据。银处理后低丰度简并基序检测减少可能反映胁迫降低甲基化组复杂性，但也可能受银处理样本较低测序深度的影响。

结论翻译：本研究通过表征十一种大肠杆菌K-12 MG1655菌株背景下的6mA和推定4mC甲基化，扩展了当前对细菌DNA甲基化组的理解。据研究人员所知，这是首个在有无银暴露条件下，比较野生型大肠杆菌K-12 MG1655、cisS突变株和银适应群体全基因组DNA甲基化的研究。主导且可重复检测的甲基化基序的鉴定，为比较遗传变异体、适应群体和胁迫暴露条件建立了基线甲基化组资源。在本研究条件下，还鉴定了非保守和菌株相关基序，以及银处理条件下基序检测的额外变异。尽管菌株间存在遗传差异且暴露于低水平胁迫，主要甲基化基序在生物学重复间基本保守，表明大肠杆菌核心甲基化系统高度稳定且具有韧性。低丰度或推定基序中的变异凸显了可能反映生物学和/或技术来源的更细微基序特异性差异的潜力。结合文献报道，研究人员推测在低水平银胁迫下，大肠杆菌K-12的主导DNA甲基化特征表现出显著韧性，与核心甲基化组的结构性而非条件响应性作用一致。这些发现为未来研究细菌DNA甲基化在环境胁迫存在与缺失条件下的生物学意义及其对细菌适应和驯化的影响奠定了基础。