《Microorganisms》：Computational Identification of Novel Transcriptional Regulators and Functional Gene Clusters in Lactococcus lactis Using Integrated Bioinformatics Approaches
编辑推荐：
乳酸乳球菌（Lactococcus lactis）是一种重要的工业“细胞工厂”，广泛应用于食品发酵和生物技术领域。然而，关于该微生物适应过程的调控机制，一个关键生物学问题仍未解决：在环境胁迫和感受态期间，细菌如何通过转录协调初级代谢与细胞表面重塑之间的权衡？迄
乳酸乳球菌（Lactococcus lactis）是一种重要的工业“细胞工厂”，广泛应用于食品发酵和生物技术领域。然而，关于该微生物适应过程的调控机制，一个关键生物学问题仍未解决：在环境胁迫和感受态期间，细菌如何通过转录协调初级代谢与细胞表面重塑之间的权衡？迄今为止，尚缺乏一个统一的、全局性的基因调控网络（GRN）模型来解释这一转变。为解决这一碎片化问题并消除选择偏差，研究人员整合了截至2025年夏季NCBI数据库中公开可用的L. lactis全部转录组数据集。该全面数据集涵盖多种条件，包括热胁迫、酸胁迫、噬菌体诱导胁迫以及自然感受态，为稳健的网络推断提供了必要的转录变异。研究人员采用集成生物信息学流程，使用GENIE3算法推断所有测试条件下共有的核心调控网络，并辅以DeepTFactor、Entraf和p2TF工具的集成组合来预测L. lactis的菌株特异性潜在转录因子（TFs）。共表达网络被划分为50个功能簇，特别突出了可能参与细胞表面修饰的独特WxL操纵子的假定调控因子。此外，研究人员提出了主感受态调控因子ComX的候选调控靶标，并通过计算预测CpsY为支链氨基酸代谢的潜在LysR家族调控因子。这些发现提供了一个基于转录组学的L. lactis调控计算模型。通过明确区分已建立的调控通路与纯粹的计算预测，研究人员提出若干未表征蛋白作为细菌响应环境挑战的假定关键节点。虽然需要直接实验验证以建立物理相互作用，但该计算方法生成了高置信度假设，并为靶向代谢工程提供了精选的候选资源。
**研究背景与问题**
乳酸乳球菌（*Lactococcus lactis* subsp. *lactis*）是一种重要的革兰氏阳性工业菌株，广泛用作乳制品发酵的起始培养物，并因其GRAS（一般认为安全）状态和缺乏内毒素而被用于重组蛋白、维生素和治疗分子的生产。然而，环境胁迫（如快速酸化、温度波动和渗透压）会限制细菌生长和代谢活性，导致工业效率下降。细菌在胁迫下生存需要快速协调的遗传响应。尽管比较基因组学已勾勒出乳酸菌的核心保守调节子，且一些全局调控因子（如CodY和CcpA）在氮代谢和碳分解代谢物阻遏中的作用已明确，但*L. lactis*调控图谱中仍存在大量“暗节点”或功能注释不足的基因。例如，利用ESM2蛋白质语言模型的大规模预测研究仅识别出四个主全局调控因子（CcpA、Spx、Nnr和CmbR），而RefSeq数据库已注释约181个位点特异性转录调控因子。这一巨大差距表明，需要基于动态表达的网络推断来揭示未表征调控节点的功能。传统实验方法鉴定转录因子及其靶标资源密集，导致转录组数据与调控逻辑理解之间的鸿沟。系统微生物学和机器学习的进步使得通过分析跨条件共表达模式进行基因调控网络（GRN）的计算机重建成为可能。本研究旨在解决的核心调控问题是：识别控制细胞壁通透性和胁迫适应的核心调控枢纽，这些枢纽无法从单条件实验中分离。因此，研究人员整合了截至2025年底NCBI数据库中38个*L. lactis* subsp. *lactis*转录组数据集，涵盖热胁迫、酸胁迫、噬菌体诱导胁迫和自然感受态，以捕获最大可能的表型和转录变异，从而迫使GENIE3算法识别通用调控模块而非条件特异性伪影。
**主要技术方法**
研究人员整合了38个公开转录组数据集（源自NCBI BioProject，涵盖酸胁迫、热胁迫、噬菌体感染和遗传扰动条件），使用标准化的计算流程重新处理原始数据。通过集成深度学习工具DeepTFactor、基序导向算法Entraf和结构域搜索方法p2TF，预测了98个潜在转录因子（TFs）。应用GENIE3随机森林算法（1000棵树）推断核心调控网络，保留前4000条边，并使用Ward.D2层次聚类将网络划分为50个功能簇。利用Operon-mapper预测操纵子结构，通过XSTREME服务器进行顺式调控元件和基序发现，并整合COG、KEGG、GO和PFAM数据库进行功能注释。
**研究结果**
- **全局拓扑与网络重建**：GENIE3算法构建的共表达网络包含1791个基因（占基因组约70%）和98个潜在转录因子，覆盖了38个样本和4种条件（pH4、热休克、基因操作和噬菌体感染）。与RefSeq注释的181个调控因子相比，仅98个在跨条件中具有足够表达变异的TFs被保留，代表活跃的调控核心。网络被划分为50个功能模块，其中Cluster 18完美富集核糖体蛋白，作为内部验证标准。
- **WxL操纵子（Cluster 11）的鉴定**：Cluster 11包含6个基因，其中4个编码含WxL结构域的细胞表面细胞壁结合蛋白，可能形成一个操纵子。该簇与一个XRE家族转录因子（LL14B4_RS03315）和另一个XRE蛋白（LL14B4_RS08280）相连，并连接至肽酶、黏蛋白结合蛋白、*dapB*和*ilvA*等基因。在*comX*过表达条件下，WxL基因、*dapB*和*ilvA*下调，而其他基因不变。基序发现显示所有WxL基因共享一个共同结合基序（ATTATAAAAG），与*B. subtilis* PerR类似，提示与氧化应激响应和感受态调控的关联。
- **ComX调节子扩展与自然感受态**：Cluster 5包含几乎所有*com*基因（包括*comX*、*comGD*、*comGC*、*comGB*、*comGA*、*comEC*、*comFC*、*comFA*、*comEA*和*comC*）以及ATP合酶操纵子（*atp*）、23个未知功能基因和*rec*操纵子基因。在*comX*过表达下，该簇90%以上基因上调（log2FoldChange > 1）；在*covRS*缺失下，多数基因下调；在*codY*缺失下，表达变化不大。这表明CovRS是此模块的主要调控因子。基序分析发现一个与CodY结合位点（ATTCATAAA）匹配的基序。同时，多个生物合成通路（如组氨酸、含硫氨基酸、支链氨基酸合成）在*comX*激活后显著下调，其中Cluster 48的*ilv*基因（BCAA合成）平均log2FoldChange < -5，且与CodY和CpsY转录因子拓扑相连。
- **CpsY在代谢重编程中的调控作用**：CpsY被鉴定为Cluster 48中的潜在转录因子，与*ilv*基因高度互连，并与未表征的Cluster 29（编码*opp*基因的ABC转运蛋白）相连。在*comX*过表达下，*cpsY*被抑制（log2FoldChange = -4.64），其相关基因也下调。Cluster 29、*cpsY*和*ilv*基因的多数（>80%）含有CodY和CcpA结合位点，提示存在一个CpsY–CodY–Ilv调控轴，可能在感受态期间通过降低BCAA合成来解除CodY对*comX*的抑制。
- **操纵子感知的再评估**：通过将预测的调控边映射到Operon-mapper预测的操纵子结构，发现关键模块（如*opp*和*com*操纵子）的调控因子精确靶向先导基因，而*ilv*操纵子等则内部基因也存在强边，反映真实的亚操纵子调控动态。
**讨论与结论**
讨论部分指出，重建网络的准确性和可靠性通过内部基准（如Cluster 18的核糖体蛋白）和与已发表实验数据的比较得到验证，成功恢复了CodY对*ilv*操纵子的抑制、CovRS对ComX感受态调节子的控制以及CcpA结合基序的识别。多层次的ComX调控涉及CodY、CcpA和CovR的抑制以及MecA-ClpCP复合物的翻译后降解。CpsY–CodY轴在感受态期间通过降低BCAA合成来解除CodY对*comX*的抑制，从而维持感受态状态。WxL操纵子可能通过下调WxL蛋白暂时改变细胞壁蛋白组成，促进DNA转运，其调控可能受PerR样氧化应激响应因子控制。研究还发现ATP合酶操纵子与*com*基因共上调，满足DNA摄取的能量需求。比较基因组学与转录组学方法的互补性得到了强调：序列模型描绘物理接线图，而共表达网络揭示动态生理串扰。研究揭示了多个新节点，如XRE家族转录因子LL14B4_RS08280可能参与细胞壁结构组织，CpsY在*L. lactis*中参与BCAA合成（不同于其同源物CmbR在链球菌中调控半胱氨酸/甲硫氨酸代谢）。局限性包括样本数量有限、缺乏物理相互作用验证以及无法捕捉转录后调控和多组学整合的需求。
**研究结论**：本研究代表了一个基于转录组推断的*Lactococcus lactis* subsp. *lactis*基因调控网络的系统级重建，预测了胁迫适应和自然感受态的假定调控模块。基于计算分析，研究人员提出了一个新型细胞壁重塑WxL操纵子（Cluster 11），揭示了CpsY–CodY–Ilv代谢轴（Cluster 48），并展示了受CovRS控制的DNA转化能量耦合（Cluster 5）。这些调控枢纽虽需直接实验验证，但为工程化具有特定生物技术特性的稳健工业菌株提供了有前景的遗传靶点。例如，调控新发现的CpsY–CodY–Ilv轴可优化支链氨基酸代谢（直接负责乳制品发酵中挥发性风味化合物生成），修饰WxL操纵子可改变表面特性以增强噬菌体抗性，而调整CovRS介导的胁迫网络可提高培养物对抗工业胁迫（如快速酸化或冷冻干燥）的存活能力。计算预测提供了一个数据驱动框架和高质量假设基础，用于后续体内验证和靶向代谢工程。
