六种工业大肠杆菌菌株转录调控网络(Transcriptional Regulatory Network, TRN)差异分析

《Metabolic Engineering Communications》：Transcriptional Regulatory Network Differences between Six Industrial Escherichia coli Strains

　　

大肠杆菌(Escherichia coli, E. coli)菌株广泛应用于众多工业与生物技术过程中，但其性能存在显著差异，这种差异无法从基因组注释中预测。由于转录调控网络(Transcriptional Regulatory Network, TRN)调控细胞运动性、胁迫响应、代谢灵活性及生产效率等功能，TRN 组织与使用的差异可能是观察到的表型差异的重要基础。为探究菌株间 TRN 差异，研究人员生成了 6 种常用工业 E. coli 菌株（BL21、C、Crooks、MG1655、W 和 W3110）共 433 套匹配 RNA-Seq 表达谱汇编，并应用 iModulon 分析比较相似生长条件下各菌株 TRN 状态。分析显示，功能相似的核心调控程序在不同菌株间存在不同的调控连接(wiring)，且相同环境胁迫下各菌株以不同方式启动这些程序。上述发现揭示了工业 E. coli 菌株间表型表达的转录调控多样性。通过提供常用宿主菌 TRN 差异的整体视图，本研究为解读菌株特异性行为提供了基础，并为更合理的菌株选择与优化提供支持。

六种工业大肠杆菌菌株转录调控网络差异的iModulon分析——论文解读

本文发表于《Metabolic Engineering Communications》。工业大肠杆菌（Escherichia coli, E. coli）K?12衍生物及B族菌株因重组蛋白生产、代谢工程和高附加值化学品合成而被广泛使用，但即便基因含量高度相似，不同菌株在生产效率、胁迫耐受性和代谢特征上仍存在明显差异。以往的比较多集中于基因组、代谢物和生理表型层面，而转录组水平的系统比较十分有限。基因组仅定义了细胞的遗传潜能，无法反映动态转录调控及资源分配逻辑，因此难以解释近缘菌株的表型分化。本研究以六种代表性工业E. coli菌株（BL21、C、Crooks、MG1655、W、W3110）为对象，通过构建匹配条件的RNA?Seq转录组汇编，利用Multi?View独立成分分析(Independent Component Analysis, ICA / Multi?View ICA) 提取iModulon（独立调控模块），从转录调控网络(Transcriptional Regulatory Network, TRN)的组织结构（基因成员组成）和条件依赖性激活（活性）两个维度，系统比较菌株间TRN差异，为菌株选型与工程改造提供机制层面的依据。

主要关键技术方法

研究人员选取6种常用工业E. coli菌株（BL21、C、Crooks、MG1655、W、W3110），整合新生成及已发表RNA?Seq数据，构建含433个转录组的匹配条件汇编（涵盖胁迫、碳源、培养基及基因敲除/进化株）；经质控与log?TPM标准化后，采用Multi?View ICA框架联合分解各菌株表达矩阵，识别跨菌株共有（核心）iModulon及菌株特异iModulon；iModulon基因成员由绝对权重经分位数过滤+Otsu阈值确定，功能注释借助eggNOG及RegulonDB（仅MG1655）；同时基于全基因组SNP系统发育树、平均核苷酸一致性(Average Nucleotide Identity, ANI)及Biolog表型微阵列量化基因组相似性与底物利用差异，并与iModulon组成/活性做关联分析。

研究结果

Exploring genomic similarities and substrate?utilization variation among strains（探索菌株间基因组相似性与底物利用差异）

通过对六菌株全基因组SNP系统发育、ANI分析及泛基因组(pan?genome)构建，发现所有菌株对ANI > 98%，MG1655与W3110 ANI > 99.9%，核心基因占59%，证实基因组高度相似。但Biolog底物利用谱显示菌株间Pearson相关系数差异显著，MG1655与W3110底物利用最接近，BL21与C较接近，而部分基因组近缘菌株（如C与MG1655）底物利用却分歧明显。基因组距离(1?ANI)与底物利用差异(1?r)无线性关联，说明基因组信息不足以解释生理表型差异，需转录组水平分析。

Mapping transcriptional regulatory architecture across industrial strains using transcriptomic analysis（利用转录组分析解析工业菌株转录调控架构）

研究人员获得每菌株69–76个高质量RNA?Seq样本（覆盖26个共有条件）。Multi?View ICA分解得到15个所有菌株共有核心iModulon及少数菌株特异iModulon。核心iModulon代表在各菌株中可重复提取的共调控基因集，反映保守生物学程序；菌株特异iModulon捕捉独特调控信号。此数据驱动方法不依赖先验调控网络，可比较有无Curated TRN图谱的菌株。

Strain?specific TRN structure revealed by differential core iModulon gene memberships（核心iModulon基因成员差异揭示菌株特异性TRN结构）

计算各菌株–iModulon保守得分（该iModulon基因在全部菌株iModulon交集中的比例），发现即使功能相同的核心iModulon，基因组成存在菌株差异。以Flagella（鞭毛）iModulon为例：BL21因基因组缺失fli/mot/che基因致iModulon缩小，符合B族菌株丧失运动性已知特征；菌株C虽保有这些基因却未纳入同一iModulon，提示调控连接(wiring)差异；W3110的Flagella iModulon异常包含LPS核心及O抗原合成基因(waa/wbb/rfb)，暗示该菌株鞭毛组装与包膜生物合成存在独特调控偶联。Copper（铜稳态）iModulon在Crooks菌株额外包含染色体整合的pco/sil抗铜/银基因岛，而菌株C虽携带同源基因却不表达也未纳入iModulon，说明相同基因含量未必导致相同转录响应。Anaerobic（厌氧）iModulon中，C、Crooks、W包含formate hydrogen lyase(FHL)复合体hyc基因，与这些菌株较高厌氧生长速率及甲酸分泌相符；BL21独有发酵型乳酸脱氢酶基因ldhA，与其独有厌氧产乳酸能力一致。表明核心iModulon可通过不同基因编组反映菌株特异调控策略。

Core iModulon activity captures differential engagement of regulatory programs across strains（核心iModulon活性反映菌株对调控程序的不同调动）

计算共有条件下各iModulon活性的株间方差，OXPHOS、Flagella、Fnr iModulon株间变异最大。酸性pH 5条件下：Iron iModulon除MG1655外均激活（MG1655需在更低pH激活）；Crooks中Fnr iModulon强烈抑制而OXPHOS iModulon强烈激活，W3110中OXPHOS被抑制，显示呼吸/厌氧回路菌株差异调控。H₂O₂氧化胁迫下：MG1655 Iron iModulon受抑制（符合Fur?OxyR/SoxRS经典途径），其余5株Iron iModulon激活；Crooks与W强烈抑制Flagella iModulon（节能避险），MG1655反而激活Flagella iModulon（可能轻度胁迫下趋避反应）。证明相同环境刺激引发菌株特异iModulon激活模式，经典E. coli调控表型不能笼统推广至所有工业菌株。

Unique iModulons reflect strain?specific regulatory signals（独特iModulon反映菌株特异调控信号）

BL21具3个特异iModulon（其一富集噬菌体相关基因），W具1个（质粒编码基因富集），W3110具1个（菌毛/部分鞭毛基因，腺嘌呤处理下受抑，与该株腺嘌呤敏感性呼应），MG1655与Crooks未发现菌株特异iModulon，可能与测试条件未覆盖菌株特有诱导因子（如W株蔗糖利用）有关。

讨论与结论（翻译浓缩）

本研究建立的比较框架强调了调控多样性是工业E. coli菌株间有意义且此前未被充分重视的差异来源。通过匹配转录组建档与Multi?View ICA，研究人员提供了六种常用生产菌株调控组织与激活差异的统一系统视角。尽管基因组高度相似，其TRN在调控程序的基因编组方式及条件依赖性启动上存在差异。核心iModulon可由不同基因组合构成，说明近缘E. coli菌株以不同方式协调共有基因，调控组织本身即是工业宿主多样性的来源。各菌株在相同环境挑战下以不同iModulon组合响应，表明调控调动(engagement)具菌株特异性。本研究的iModulon结构与活性谱可辅助后续基因调控网络建模及菌株工程设计，为理解工业宿主基因型?表型间隙提供转录组层面的中间层解释，并可指导特定工艺条件下的菌株优选与改造。