|
猪鼻支原体SK76甲基化组与限制-修饰系统（R-M系统）的多组学整合表征摘要 限制-修饰（R-M）系统在支原体（Mycoplasma）物种中广泛存在，参与抵御外源DNA和表观遗传调控。然而，柔膜体纲（Mollicutes）有限的遗传可操作性阻碍了R-M系统的功能表征。本研究通过整合比较基因组学、结构建模、定量蛋白质组学和全基因组甲基化组分析，对猪鼻支原体（Mycoplasma hyorhinis）SK76菌株的R-M系统谱系进行了表征。研究人员发现SK76编码I型、II型和III型R-M系统，且部分组分存在重复。对于I型系统，预测HsdS2识别一个二部分的TA-N?-ATC基序，该基序被检测为次要的6mA相关甲基化特征，而其余HsdS变体则表现出与功能减弱一致的特征。在II型R-M系统中，Mhy3及其同源核酸内切酶End3与5mC基序GATC相关联，而其余预测的甲基转移酶则缺乏蛋白质表达和相应的甲基化特征。关于III型系统，Mod1被鉴定为负责不对称GGAT中心甲基化模式的主要活性腺嘌呤甲基转移酶，而其他Mod同源物可能代表相变或条件调控元件，活性有限。总体而言，猪鼻支原体（M. hyorhinis）的功能性甲基化组仅代表其基因组R-M系统谱系的一个子集，这凸显了整合基因组学、蛋白质组学和甲基化组数据以定义活性表观遗传系统的价值。 论文解读文章一、研究背景与科学问题猪鼻支原体（Mycoplasma hyorhinis）是一种无细胞壁的柔膜体纲（Mollicutes）细菌，主要被认识为猪病原体，同时也是哺乳动物细胞培养的常见污染物。除兽医相关性外，该菌因在人类肿瘤中频繁检出而日益受到生物医学关注，包括胃癌（56%）、结直肠癌（55%）和肺癌（53%）。持续感染和不同细菌蛋白已被牵涉到炎症信号传导、增强细胞迁移和肿瘤进展中。重要的是，实验证据表明猪鼻支原体的限制-修饰（R-M）系统和解旋酶有助于细菌毒力和宿主相互作用。此外，来自猪鼻支原体的II型DNA甲基转移酶已被证明在入侵后修饰宿主基因组中的胞嘧啶残基，改变与人类增殖和癌变相关基因的甲基化模式。这些发现表明DNA甲基化不仅在细菌基因组防御中发挥重要作用，还可能参与细菌-宿主相互作用和致病性。 尽管有这些相关证据，猪鼻支原体编码的DNA甲基转移酶的全谱仍未完全明确。在细菌中，DNA甲基化主要由R-M系统介导，这是区分自身与非自身DNA的最广泛机制之一。这些系统通常包含切割未甲基化外源DNA的限制性核酸内切酶（R）和修饰宿主基因组内特定序列基序的相应DNA甲基转移酶（M），从而保护自身DNA免于切割。除防御外源DNA的经典作用外，R-M系统还牵涉基因表达调控、基因组稳定性、水平基因转移和细菌适应。R-M相关甲基转移酶修饰特定胞嘧啶或腺嘌呤残基，产生5-甲基胞嘧啶（5mC）、4-甲基胞嘧啶（4mC）和N6-甲基腺嘌呤（6mA）等表观遗传标记。具有甲基化活性的R-M系统根据亚基组成、DNA识别和催化要求的差异分为三大类（I-III型）。 尽管有证据表明猪鼻支原体编码几种推定的DNA甲基转移酶，但其DNA甲基化系统的多样性和序列特异性仍未完全表征。以往研究通过在大肠杆菌（Escherichia coli）中异源表达三种推定的猪鼻支原体SK76甲基转移酶（Mhy）酶，随后进行基于亚硫酸盐的分析，鉴定出CATG和CG DNA基序上的5mC甲基化。然而，亚硫酸盐测序能有效检测5mC，但不能可靠区分4mC（可能被误认为5mC），且不检测6mA。因此，仅通过亚硫酸盐测序表征的细菌甲基化组可能低估DNA甲基化标记的多样性。第三代测序技术（包括单分子实时测序和牛津纳米孔测序）的进展，现在能够在单分子分辨率下直接检测多种DNA碱基修饰（即6mA、5mC和4mC），无需化学转化或扩增偏差。这凸显了利用能够检测多类DNA修饰的技术对猪鼻支原体甲基化组进行全面、基序分辨表征的必要性。 二、研究内容与主要结论本研究提供了猪鼻支原体SK76 DNA甲基化景观的全面表征。通过整合基因组注释、结构建模、蛋白质组学和牛津纳米孔测序，研究人员定义了其限制-修饰系统的谱系和组织结构，鉴定了在分析条件下表达的组分，并将这些系统与全基因组检测到的甲基化模式联系起来。这一整合分析还使现有注释得以完善，并提供了该菌株中三类限制-修饰系统的统一视图。 研究发现SK76基因组编码完整的I型R-M系统，包含典型的HsdR-HsdM-HsdS结构，包括一个HsdM甲基转移酶、一个HsdR限制性亚基和三个推定的HsdS特异性亚基。其中，HsdS2因其较高的蛋白质丰度和结构完整性，被确定为主要特异性亚基，最可能识别二部分TA-N?-ATC基序。II型系统中，三个推定的胞嘧啶甲基转移酶（Mhy1-Mhy3）中仅Mhy3及其同源限制性蛋白End3被蛋白质组学可靠检测，且与GATC中心的5mC基序相关。III型系统中，Mod1被鉴定为主要活性腺嘌呤甲基转移酶，负责不对称GGAT中心的6mA甲基化模式，而其他Mod同源物（Mod2A、Mod2B、Mod3A、Mod3B）可能代表相变或条件调控元件，活性有限。 三、关键技术与方法研究人员主要采用了以下关键技术方法：（1）基于Pfam数据库的隐马尔可夫模型（HMM）结构域搜索，结合Clustal Omega多序列比对，鉴定和验证R-M系统同源物；（2）AlphaFold 3和AlphaFold-Multimer结构预测，用于建模蛋白质结构及蛋白质-DNA复合物，评估亚基间结构兼容性；（3）定量蛋白质组学，采用数据非依赖采集（DIA）结合DIA-NN分析，检测指数期和稳定期培养条件下R-M组分的蛋白质表达；（4）牛津纳米孔测序（GridION Mk1平台，R10.4.1流动池），结合MicrobeMod和STREME进行甲基化碱基检测和基序发现；（5）比较基因组学分析，涵盖23株猪鼻支原体菌株的R-M基因保守性评估。样本来源为猪鼻支原体SK76菌株（由Kim S. Wise提供），基因组DNA提取自30 mL稳定期培养物，三个生物学重复的DNA合并后建库测序。 四、研究结果猪鼻支原体SK76中I型R-M组分的计算机表征与蛋白质组学证据 通过结构域搜索在SK76基因组中鉴定出完整的I型R-M系统，包括一个HsdM、一个HsdR和三个推定的HsdS特异性亚基。hsdS1和hsdS2基因以紧凑的反向配置排列，有25 bp重叠，而hsdS3位于远端染色体位点。通过序列分析发现，SK76 NC_019552.1和CP003914.1注释中hsdS1的起始密码子定义错误，实际起始密码子位于两个变体上游。AlphaFold 3结构建模显示三个变体具有相似的整体折叠，包含两个球状DNA靶标识别结构域（TRD1和TRD2），由α-螺旋连接区相连。序列同一性分析显示三个蛋白质间存在显著差异：HsdS2和HsdS3在中心保守区和TRD2共享高同一性（96.55%和92.37%），但TRD1仅48.15%；HsdS1在所有区域更为分化。HsdS2包含完整的卷曲螺旋中心螺旋区，而HsdS1含有较短的螺旋，HsdS3缺少该区域的第二个α-螺旋，提示其可能是假基因。HsdM序列分析鉴定出保守的NPPY催化基序和Rossmann样AdoMet结合核心，支持I型甲基转移酶的催化能力。HsdR含有PF04313 I型限制酶R结构域以及DEAD/DEAH-box解旋酶基序，与ATP依赖性DNA转位和限制活性一致。蛋白质组学检测到HsdS1和HsdS2在指数期和稳定期均表达，而HsdS3无独特肽段可确认。HsdS1的平均蛋白质强度比HsdS2低约3个log?单位。AlphaFold-Multimer M?S组装模型显示，所有三种HsdS变体与HsdM的组合均显示相似的pTM/ipTM分数，高于阳性对照（EcoKI复合物），远高于阴性对照（非同源2Mhy1-HsdS1）。综合HsdS1低丰度和HsdS3截短的特征，支持HsdS2为I型系统最可能的主要特异性亚基。 猪鼻支原体SK76中II型R-M组分的计算机表征与蛋白质组学证据 在独立染色体位点鉴定出三个推定的II型胞嘧啶甲基转移酶（Mhy1、Mhy2和Mhy3），每个甲基转移酶基因均位于一个限制性相关蛋白（End1、End2和End3）旁，其中mhy3和end3基因部分重叠。系统发育分析显示，每个限制性相关蛋白与其相邻甲基转移酶的聚类一致：End1和Mhy1与来自无乳支原体（M. agalactiae）的MAG4270相关蛋白聚类，End2和Mhy2与来自微小脲原体（Ureaplasma parvum）的UPV_230和UPV_229聚类，End3和Mhy3与来自金黄色葡萄球菌（Staphylococcus aureus）的Sau3AIR和Sau3AIM相关蛋白聚类。蛋白质组学分析显示显著表达差异：Mhy3及其相关限制性蛋白End3被高丰度可靠检测；Mhy1和Mhy2均未检测到；End1未检测到，End2仅在置信度阈值处检测到且支持肽段有限。结构分析聚焦于Mhy3-End3对，AlphaFold 3模型显示与金黄色葡萄球菌Sau3AI甲基转移酶和限制酶高度结构相似（RMSD分别为0.4 ?和1.7 ?）。Mhy3序列分析鉴定出C5胞嘧啶甲基转移酶特征的保守PCQ催化基序和参与S-腺苷-L-甲硫氨酸结合的FxGxG基序。End3保留PD-(D/E)XK核酸酶基序，与催化活性的II型限制性核酸内切酶一致。 猪鼻支原体SK76中III型R-M组分的计算机表征与蛋白质组学证据 III型R-M系统组织为紧凑的基因组位点，包含多个串联聚集的Mod相关基因（mod1、mod2B、mod2A、mod3B和mod3A），紧邻预测的Res限制性亚基（res1）。在Mod序列中，仅Mod1保留甘氨酸丰富的FxGxG AdoMet结合基序和活性III型甲基转移酶特征的催化DPPY基序。mod1基因含有一个基因内11重复AG丰富短序列重复（SSR），SSR可在DNA复制过程中发生滑链错配，产生可改变阅读框的重复长度变体，这一过程称为相变。比较RefSeq和GenBank注释发现Mod2A和Mod3A预测起始密码子的差异：根据GenBank注释，Mod2A和Mod3A是截短蛋白，保留FxGxG AdoMet结合基序但缺乏催化DPPY基序，不太可能作为独立甲基转移酶发挥作用。在mod2A-mod2B（Mod2）和mod3A-mod3B（Mod3）之间鉴定出SSR，序列重建显示分别缺失两个碱基和一个碱基可恢复连续的mod2和mod3开放阅读框，编码包含两个保守基序的蛋白质。系统发育分析将Mod1与来自其他支原体物种的全长Mod甲基转移酶归为一组，而Mod2和Mod3形成密切相关但较短的分化分支。Res1与来自肺炎支原体（M. pulmonis）的特征明确的III型限制酶显示结构相似性，保留典型结构域架构，并具有保守的Walker A和Walker B ATPase基序。蛋白质组学分析鉴定出分配给Mod1、Mod2B和Mod3B的独特肽段，而Mod2A、Mod3A或推定的SSR恢复全长产物因高序列相似性无法可靠检测。Mod1和Mod2B丰度相似，Mod3B丰度较低。mod2-mod3区域的截短组织在多个猪鼻支原体菌株中保守。Res1丰度与Mod1相似。 I、II和III型R-M系统基因在猪鼻支原体菌株中广泛保守且大多必需 以SK76为参考，分析了22株猪鼻支原体菌株中预测的I、II和III型R-M基因的保守性。I型HsdM、HsdR和HsdS亚基的同源物在大多数菌株中鉴定到，特异性亚基通常显示>70%序列同一性。三个II型胞嘧啶甲基转移酶（Mhy1-Mhy3）在菌株间高度保守，相应限制性核酸内切酶大多保留；End1在JF5820及其衍生菌株中缺失，IVB14_OD_0817含End2基因但缺Mhy2同源物。III型R-M系统大体保守，仅22株中4株缺失相关基因。来自JF5820的必需性数据显示，三种R-M类型的所有甲基转移酶均被分类为必需，End1和End3为非必需，End2为必需。III型系统中Mod1为必需，Res1为非必需。Mod1在菌株间的高保守性及其必需性，凸显其在维持猪鼻支原体SK76 III型甲基化系统中的核心作用。 猪鼻支原体SK76具有多层DNA甲基化景观 通过抗6mA抗体的斑点印迹分析确认了SK76基因组DNA中存在6mA，信号强于大肠杆菌DH5α野生型DNA，而大肠杆菌dam?/dcm?菌株DNA无信号。全基因组6mA、5mC和4mC修饰通过牛津纳米孔测序和MicrobeMod分析。在STREME阈值0.5下，鉴定出两个主要的6mA相关基序和一个低覆盖基序：第一个为非回文序列5′-WGGATNWNNW-3′，检测到2,333个基因组位点，77.2%甲基化覆盖；第二个为二部分基序5′-TAWNNWNNATCW-3′，检测到890个位点，47.2%覆盖；第三个为低覆盖基序WGATCW，A主要在正链甲基化。对于5mC，主要基序为回文序列5′-WGATCW-3′（2,320个位点，90.1%覆盖），另有低覆盖的5′-WNNCGNWNNWNNW-3′基序。4mC基序CATCC仅在极低频率检测到。在STREME阈值0.9下，主要6mA信号保留GGAT中心核心（1,956个位点，77.4%覆盖），主要5mC信号解析为延伸的回文GATC中心基序（824个位点，93.2%覆盖）。二部分6mA相关基序（TA-N6-ATC配置）在较高阈值下未恢复。甲基化信号广泛分布于基因组，6mA比5mC更广泛，仅检测到7个4mC位点，无稳健的基序相关4mC信号证据。 活性甲基转移酶组分和推定同源限制酶的整合表征 由于猪鼻支原体遗传可操作性有限，研究人员整合基因组组织、保守催化基序、蛋白质组学证据、结构建模和纳米孔甲基化组谱，将候选甲基化基序分配给SK76中鉴定的R-M组分。I型系统中，HsdS2是最可能识别二部分5′-TAWNNWNNATCW-3′基序的HsdS变体，基于其结构完整性和蛋白质组检测。II型系统中，Mhy3和End3是唯一被蛋白质组学可靠检测的甲基转移酶-限制酶对，Mhy3靶向GATC包含基序，End3属于同样识别GATC的Sau3AI家族，综合MicrobeMod结果，次要甲基化基序更佳表示为5′-CGATCG-3′。III型系统中，Mod1被确定为主要候选负责5′-GGAT-3′中心6mA信号，基于保守催化基序的存在和蛋白质组检测。HsdR1、End3和Res1被蛋白质组学检测并保留与限制活性一致的保守序列或结构特征，但DNA切割活性未直接测试，因此被视为推定同源限制组分。 五、讨论与结论总结讨论部分指出，本研究通过整合比较基因组学、结构建模、蛋白质组学和全基因组甲基化谱，提供了猪鼻支原体SK76 R-M景观的综合视图，展示了如何在遗传难处理的生物（如某些柔膜体纲）中重建功能性表观遗传系统。SK76揭示了由I型、II型和III型R-M模块组成的多层表观遗传架构，各模块显示不同程度的保守性、表达和对全局甲基化组的贡献。这些模块可能构成分层表观遗传架构，其功能输出受基因保留、蛋白质丰度和潜在相变机制调节。 I型系统中，HsdS2因较高蛋白质丰度和结构完整性成为最可能的特异性亚基，HsdS1和HsdS3可能处于假基因化过程。与通过SSR驱动HsdS相变的其他柔膜体纲不同，SK76缺乏此类重复，表明更稳定的架构，支持其在基因组防御和表观遗传调控中的可能作用。II型系统中，仅Mhy3得到蛋白质组检测和相应甲基化特征支持，而Mhy1和Mhy2虽此前证明在大肠杆菌中活性，但在测试条件下未检测到蛋白质表达和相应甲基化信号，表明SK76中胞嘧啶甲基化由单一活性II型系统驱动。End2的极低丰度检测及其在JF5820中的必需性预测之间的差异，提示该位点的功能相关性可能因菌株而异。III型系统中，Mod1是最可能的GGAT中心6mA甲基化负责酶，而Mod2和Mod3变体可能代表相变库，具有在亚群中产生甲基转移酶多样性的潜力，但缺乏融合产物的蛋白质组学证据及截短架构的跨菌株保守性，不支持在测试条件下有效恢复替代全长酶。关于4mC信号（CATCC），仅两个甲基化位点支持且覆盖极低（0.016），最可能代表背景噪声。推定的6mA基序AAGATC虽覆盖较低（0.088），但由56个甲基化位点支持，不能完全排除其存在，可能反映测试条件下不活跃或弱活性的额外腺嘌呤甲基转移酶的残留或条件依赖性活性。 研究结论表明，猪鼻支原体SK76的功能性甲基化组仅代表其基因组R-M系统谱系的一个子集。整合基因组学、蛋白质组学和甲基化组数据对于定义活性表观遗传系统具有重要价值，为理解该物种的表观遗传调控、宿主相互作用和致病性提供了框架和基础。 |
