副溶血性弧菌低细胞密度下ToxR调控子的转录组解析

《International Microbiology》：Transcriptomic dissection of the ToxR regulon in Vibrio parahaemolyticus at low cell density

副溶血性弧菌（Vibrio parahaemolyticus）是一种主要的海产品传播病原菌，其致病性与环境适应性依赖于调控因子ToxR，但该调控子的范围此前尚未明确。研究人员采用RNA-seq技术，在toxR表达峰值期比较了野生型RIMD2210633菌株与同基因ΔtoxR菌株的转录组，共鉴定出603个差异表达基因（突变株中348个下调、255个上调）。在所测试条件下，ToxR主要作为基因表达的激活因子发挥作用。该调控子包含32个调控基因、13个c-di-GMP代谢基因、III型分泌系统1（T3SS1）基因簇（在突变株中上调）、所有可检测到的T3SS2基因及tdh2（在突变株中下调）、差异调控的VI型分泌系统（T6SS1/T6SS2）基因、侧生鞭毛基因、scv基因、mfpABC操纵子及多种孔蛋白。上述发现确立了ToxR作为整合毒力、应激适应、生物被膜发育及第二信使信号转导相关基因表达的核心调控枢纽（hup）的地位，为理解副溶血性弧菌在宿主定殖与环境转换过程中如何协调其复杂生活方式提供了全面框架。

论文解读

一、研究背景与科学问题

副溶血性弧菌（Vibrio parahaemolyticus）是一种革兰氏阴性嗜盐菌，是全球范围内海产品源胃肠炎的首要致病菌。其致病性依赖于多种毒力因子，包括耐热直接溶血素（TDH）、TDH相关溶血素（TRH）、两套III型分泌系统（T3SS1和T3SS2）以及两套VI型分泌系统（T6SS1和T6SS2）。其中T3SS1主要介导对宿主细胞的细胞毒性，而T3SS2与TDH共同决定肠毒性及肠道定殖能力。这些毒力因子的表达受到温度、pH、盐度和胆汁盐等环境信号的严格调控。

ToxR是一种跨膜转录调控因子，最早在霍乱弧菌（V. cholerae）中被鉴定为霍乱毒素和毒素共调菌毛表达的主控因子。在副溶血性弧菌中，ToxR由toxRS操纵子编码，与周质伴侣蛋白ToxS共同作用。已有研究表明，ToxR可直接激活tdh2和vtrA的转录，间接激活T3SS2装置操纵子；同时通过多层次机制抑制T3SS1表达，包括直接激活LysR型调控因子CalR以及直接结合T3SS1启动子区域。ToxR还与群体感应（QS）调控因子AphA和OpaR协同抑制T6SS1，并通过CalR间接支持T6SS2表达。此外，ToxR还抑制c-di-GMP合成，参与生物被膜形成、鞭毛介导的运动性以及胆汁耐受和酸耐受。

尽管上述靶向研究已揭示了ToxR的多种调控功能，但对该调控因子在副溶血性弧菌中全基因组范围的调控网络仍缺乏系统性认识。因此，本研究旨在通过RNA-seq技术，在ToxR表达达到峰值的低细胞密度（LCD，OD₆₀₀=0.2）条件下，系统比较野生型菌株与ΔtoxR菌株的转录组差异，从而全面解析ToxR调控子的组成与功能。

二、研究设计与主要技术方法

研究人员以副溶血性弧菌RIMD2210633为野生型亲本菌株，使用此前构建的无标记非极性ΔtoxR菌株作为比较对象。两株菌均在HI培养基中于37°C振荡培养至OD₆₀₀=0.2（即toxR转录达到最大值的生长阶段），每株各取三个独立生物学重复进行RNA-seq分析。野生型RNA-seq数据来自已发表的NCBI BioProject PRJNA1235461数据集，ΔtoxR数据集为本研究新生成。测序在Illumina NovaSeq X Plus平台完成，以副溶血性弧菌RIMD2210633基因组（NC_004603.1和NC_004605.1，共4832个注释基因）为参考进行比对。差异表达基因（DEG）的判定标准为绝对倍数变化>2且q值（FDR, padj）≤0.05，分析使用DESeq2软件包完成。功能富集分析分别采用Gene Ontology（GO）、Clusters of Orthologous Groups（COG）和Kyoto Encyclopedia of Genes and Genomes（KEGG）数据库。此外，研究人员选取17个代表性基因进行qPCR验证，以recA为内参基因，采用2^?ΔΔCt方法计算相对表达量，并与RNA-seq数据进行Pearson和Spearman相关性分析。

三、研究结果

RNA-seq数据质量分析

六个Illumina文库各产生超过1420万条原始读段。经质控和过滤后，clean reads比对到参考基因组的总比对率介于98.9269%至99.3786%之间，唯一比对率介于95.2310%至95.4528%，多重比对率介于3.6740%至4.1456%，表明数据质量满足下游分析要求。

ToxR依赖性转录组概览

共评估了4729个注释基因，鉴定出603个DEG，其中348个在ΔtoxR菌株中下调、255个上调。COG注释显示，最大的功能类别为“一般功能预测”、“信号转导机制”和“功能未知”，其次为“无机离子转运与代谢”、“细胞壁/膜/包膜生物合成”和“氨基酸转运与代谢”。KEGG富集将大多数DEG归入代谢通路（189个基因），其次为细胞过程（22个）、环境信息处理（16个）和人类疾病（6个）。GO富集显示“致病性”和“离子转运”在生物学过程中包含的DEG数量最多。

调控子中的关键功能群

ToxR调控子包含32个调控基因（ΔtoxR菌株中18个下调、14个上调）、13个c-di-GMP代谢基因（9个下调、4个上调）、8个侧生鞭毛基因（全部上调）、scv胞外多糖（EPS）位点的5个基因（全部下调）、完整的T3SS1基因簇（34个基因，全部上调）、1个IV型菌毛基因mshF（下调）、6个T6SS1基因（主要上调，VP1418除外）以及5个T6SS2基因（下调）。此外还包括tdh2（下调）、8个T3SS2基因（全部下调）、两个抗氧化基因katG1和katG2（均下调）、mfpABC操纵子（全部上调）以及9个外膜蛋白基因（5个下调、4个上调）。在调控因子中，下调幅度最显著的包括calR（29.4倍）和scvE（31.5倍），而TetR/AcrR家族调控因子VP1229上调4.0倍。在c-di-GMP代谢方面，ToxR正调控多个磷酸二酯酶基因（如VP1423、VP1427、VP1768和VP1979），并抑制多个二鸟苷酸环化酶基因（如VP2888、VPA0184、VPA0360和VPA0556）。

qPCR验证

选取的17个基因的qPCR结果与RNA-seq数据在方向上基本一致。两个数据集经log₂转换后呈正相关（Pearson r=0.65，R²=0.42，P<0.01；Spearman ρ=0.72，P<0.01）。但相关性并非完全一致，部分RNA-seq鉴定出的DEG在qPCR中以2倍阈值衡量时未被重现，说明qPCR数据为RNA-seq结果提供了部分验证，既存在方向一致性也存在阈值水平的差异。

四、讨论与结论

本研究提供了副溶血性弧菌ToxR调控子在低细胞密度条件下的首个全基因组视图。由于突变株中大多数DEG为下调，ToxR在所测试条件下主要作为激活因子发挥作用。调控子中包含32个转录调控因子基因，表明ToxR通过直接启动子结合和下游调控级联双重机制控制基因表达。其中calR下调29.4倍，与已知的ToxR-CalR调控层级一致；opaR下调2.5倍，与OpaR抑制T3SS1和T6SS1表达、激活荚膜多糖表达的功能相符。

在c-di-GMP代谢方面，ToxR激活多个磷酸二酯酶基因并抑制多个二鸟苷酸环化酶基因，与ToxR抑制c-di-GMP合成的已知功能一致。在生物被膜相关基因中，scv EPS位点强烈下调表明ToxR激活EPS生物合成基因表达；mshF下调提示ToxR激活MSHA菌毛基因表达；而mfpABC操纵子上调则提示ToxR抑制其表达，这一发现与此前报道存在差异，可能反映了通过CalR或其他调控因子的间接调控，或菌株、培养基、生长阶段及标准化方法的差异。

T3SS1基因簇全部上调与ToxR抑制T3SS1表达的功能一致，这种抑制可能涉及ToxR依赖性激活calR和opaR以及ToxR直接结合T3SS1启动子。相反，T3SS2基因和tdh2下调与ToxR激活其表达的功能一致。T3SS1与T3SS2/TDH表达的相互调控可能反映了它们在环境生存与肠道致病中的不同角色。T6SS1基因主要上调而T6SS2基因下调，与ToxR直接抑制T6SS1表达并通过CalR支持T6SS2表达的已知功能一致，也与ToxR优先支持宿主定殖而非抗菌竞争的模型相符。

所有侧生鞭毛基因均上调，表明ToxR抑制其表达，这与ToxR激活生物被膜相关表达而抑制运动相关表达的模型一致，但与ToxR为运动性所需的表型观察存在表面矛盾，可能反映生长阶段特异性效应或游泳与群集运动的差异。ToxR还强烈激活chiP和VP2467表达，并抑制多个孔蛋白基因，同时激活katG1和katG2表达，提示ToxR参与氧化应激响应基因表达程序的激活。此外，硫、氨基酸和脂肪酸代谢基因的差异表达表明ToxR调控多种代谢通路的基因表达。

总体而言，本研究支持ToxR在低细胞密度下作为基因表达核心调控因子的模型：ToxR抑制T3SS1、T6SS1和侧生鞭毛基因表达；激活T3SS2和tdh2表达；激活scv EPS和MSHA菌毛基因表达而抑制mfpABC；调节孔蛋白表达；并协调代谢和氧化应激响应表达。通过32个ToxR依赖性调控因子基因，ToxR建立了整合群体感应、c-di-GMP信号和应激反应的多层级调控网络。该模型基于单一时间点的转录组数据，qPCR验证显示了方向一致性但存在阈值水平差异。未来研究应结合ChIP-seq或ChIP-exo、调控因子突变体的上位性分析、表型测定以及跨细胞密度和宿主相关条件的时间过程实验，以确定ToxR的直接靶标和功能后果。

结论：这项全基因组分析揭示了ToxR作为副溶血性弧菌中整合毒力、生物被膜相关、运动、代谢和应激响应程序的主控调控因子。调控子包括T3SS1、T3SS2、TDH、T6SS1、T6SS2、scv、mshF、侧生鞭毛、32个转录调控因子和c-di-GMP代谢基因。毒力基因表达的相互调控可能反映了对不同生态位的适应性。确定直接靶标和表型后果仍是理解该病原菌ToxR介导的生活方式转换的关键目标。