《Current Research in Microbial Sciences》：Modelling
Fusobacterium lifestyles transitions by integrating transcriptomics and growth data
编辑推荐：
本研究通过构建并实验验证具核梭杆菌的基因组尺度代谢模型（GEM），并结合不同生理状态（自由生长、粘附、侵袭）的转录组数据，系统揭示了该菌在宿主微环境中代谢重编程的动态特征。研究发现，粘附阶段伴随支链氨基酸分解代谢抑制和特定营养摄取增加，而侵袭阶段则部分恢复了中心碳氮代谢通路活性，同时短链脂肪酸生成和氧化还原平衡的改变可能促进细菌持久定植并调控肿瘤微环境。这项研究为理解具核梭杆菌在结直肠癌等疾病中的作用提供了系统的代谢视角，并为靶向其代谢脆弱性提供了新思路。
在我们身体的微观世界里，居住着数以万亿计的细菌，它们中的大多数与我们和平共处，甚至互利互惠。然而，其中一些“坏邻居”却与疾病的发生发展密切相关。近年来，科学家们发现，生活在我们口腔中的一种名为具核梭杆菌（Fusobacterium nucleatum）的厌氧菌，竟然在结直肠癌组织中异常富集，并与癌症的进展和不良预后紧密关联。这种细菌就像一个“机会主义者”，它不仅与牙周炎、阑尾炎等多种感染性疾病有关，更被发现能粘附并侵入人类肠道上皮细胞，在肿瘤内部“安家落户”，直接或间接地促进肿瘤生长、帮助癌细胞逃避免疫攻击，甚至导致化疗耐药。
尽管人们已经认识到具核梭杆菌在癌症中的重要作用，但对于它如何在复杂的肿瘤微环境中生存、适应并施加影响，其内在的代谢机制仍是一团迷雾。细菌的代谢是其在多变环境中生存和发挥功能的核心，理解具核梭杆菌在从自由生活到粘附宿主细胞，再到侵入细胞内部这一系列“生活方式”转变过程中的代谢变化，对于揭示其致病机理和寻找干预靶点至关重要。然而，目前尚缺乏一个系统性的工具来模拟和解析这一动态过程。
为了解决这一难题，一项发表在《Current Research in Microbial Sciences》上的研究应运而生。研究人员致力于通过整合多组学数据，构建一个能够模拟具核梭杆菌在不同生理状态下代谢网络的“数字孪生”模型。这项研究的核心目标是：构建一个经过实验验证的具核梭杆菌基因组尺度代谢模型（GEM），并利用公开的转录组数据，模拟该菌在自由生长、粘附及侵袭人类细胞三种不同状态下的代谢通量变化，从而揭示其适应宿主环境的代谢重编程策略。
为开展这项研究，作者团队运用了几个关键的技术方法：首先，他们对具核梭杆菌ATCC 25586菌株进行了生理学表征，包括在不同稀释度下绘制生长曲线以确定指数生长期，并在特定时间点（T1, T2, T3）取样进行RNA测序，分析差异表达基因。其次，他们提取了生物质组分（蛋白质、RNA、DNA、脂质、糖原），为构建菌株特异的生物质方程提供实验数据。研究的核心是构建并验证基因组尺度代谢模型（Fn-GEM），他们从KBase平台的自动草稿模型出发，通过手动比对KEGG、BioCyc等数据库以及参考E. coli的同源关系，进行精细的模型重构与验证，最终得到一个包含1041个反应、938种代谢物和452个基因的模型，并验证了其在不同限定培养基中预测生长速率的能力。最后，为了模拟生活方式转变，他们进行了转录组数据的整合与通量分析，利用来自F. animalis菌株（因其与F. nucleatum核心代谢基因高度保守）在控制、粘附和侵袭状态下的公开RNA-seq数据，通过GIMME算法将表达数据映射到Fn-GEM上，生成情境特异性模型，并使用OptGPSampler进行通量取样分析，通过统计检验识别不同状态间发生显著变化的代谢反应。
研究结果
3.1. 具核梭杆菌在复杂培养基中的生理特性
研究人员首先在富含营养的YCFA培养基中培养了具核梭杆菌，确定了其生长动力学，并选择了1/10的稀释度作为后续实验条件。他们在指数生长期的三个时间点（T1, T2, T3）取样，进行了RNA测序和生物质组分分析。结果显示，蛋白质是细胞干重的主要成分（约40%），其次是RNA（20%）、脂质（9%）和糖原（2.1%）。转录组分析揭示了随着生长阶段变化，基因表达谱发生显著改变，其中氨基酸代谢和转运相关基因的表达变化尤为突出，强调了氨基酸代谢对该菌的重要性。
3.2. 具核梭杆菌基因组尺度代谢模型（Fn-GEM）的重构
通过“自下而上”的手动重构和广泛验证，研究团队成功构建了具核梭杆菌ATCC 25586的GEM模型。该模型能够较好地预测菌株在多种化学限定培养基（CDM, CDMG, CDMF, CDMFII）中的生长速率，尽管存在一定程度的低估或高估，但总体趋势一致，证明了模型在模拟不同营养条件下生长表型的可靠性。
3.3. 梭杆菌属保守核心网络的情境特异性建模
为了探究具核梭杆菌在生活方式转变中的代谢变化，研究人员将F. animalis菌株（与F. nucleatum核心代谢基因高度保守）在控制（自由生长）、粘附和侵袭状态下的转录组数据整合到Fn-GEM中。通过GIMME算法构建情境特异性模型并进行通量取样分析，他们发现从自由生长到侵袭的过程中，模型网络规模逐渐缩小（反应数量从1041个减少至609个），但显著变化的反应主要涉及特定通路的精细调控，而非整个代谢过程的完全开启或关闭。
3.4. 从自由生活到宿主相关状态的整体网络重编程
综合分析表明，粘附状态的特征是支链氨基酸（BCAA）分解代谢的抑制，以及丝氨酸和蛋氨酸摄取增加，同时碳通量向混合酸发酵途径转移。相比之下，向侵袭阶段的过渡则涉及BCAA生物合成的部分重新激活、一碳代谢的重塑以及氮通量的更广泛再分配。
3.5. 从自由生活到宿主粘附代谢状态的转变
粘附诱导了深刻的代谢重编程。氨基酸利用方面，多种必需氨基酸（如亮氨酸、谷氨酰胺、缬氨酸）的摄取减少或逆转，而蛋氨酸和丝氨酸的摄取显著增加。丝氨酸分解为丙酮酸和铵（SERD_L）的活性也增强。中心碳代谢方面，丙酮酸甲酸裂解酶（PFL）通量增加，导致甲酸盐分泌增多，这被认为是一种关键的致癌代谢物。此外，柠檬酸摄取和TCA循环部分酶（如ACONT, ICDHyr）的活性增强，同时乙酰辅酶A ACP转酰基酶（ACOATA）通量增加，暗示脂质生物合成可能增强，这可能与细菌外膜囊泡（OMVs）的产生有关，OMVs在粘附过程中起关键作用。果糖摄取和利用也增加。氧化应激防御相关的硫氧还蛋白过氧化物酶（THIORDXi） 活性升高。
3.6. 在人结肠上皮细胞内定植期间的代谢网络适应
侵袭阶段的代谢变化反映了细菌对宿主细胞内环境的适应。丝氨酸利用减少。色氨酸代谢（TRPS2, TRPAS2）活性降低，恢复到接近自由生长状态的水平。异亮氨酸转氨酶（ILETA）通量反向，表明转向BCAA的生物合成以满足需求。谷氨酰胺摄取（GLNt4）和NAD+依赖性谷氨酸脱氢酶（GLUDxi）活性增强，可能用于应对肿瘤微环境中的营养需求。发酵途径方面，D-乳酸脱氢酶（LDH_D）通量增加，表明乳酸被重新利用。丙酮酸甲酸裂解酶（PFL）活性仍保持高位但略有下降。丁酸盐代谢相关反应（BUTKr, PBUTT, BUTCT）通量小幅下降，显示出动态平衡。TCA循环早期步骤（如柠檬酸摄取、ACONT, ICDHyr）活性降低。叶酸代谢中，亚甲基四氢叶酸脱氢酶（MTHFD）活性大幅增强，而MTHFC反应几乎失活。
研究结论与意义
本研究成功构建并验证了具核梭杆菌的基因组尺度代谢模型（Fn-GEM），并利用跨物种转录组数据模拟了其在自由生长、粘附和侵袭三种状态下的代谢网络重编程。研究揭示，粘附阶段伴随着支链氨基酸分解代谢的抑制、特定营养物（如蛋氨酸、丝氨酸）摄取的增加，以及向混合酸发酵的代谢转变，这可能与能量生产、氧化还原平衡以及支持外膜囊泡生成的脂质生物合成有关。侵袭阶段则表现出中心碳氮代谢的部分重新激活、谷氨酰胺利用增强以及发酵产物的再利用，这些变化可能有助于细菌在细胞内环境中的生存和增殖。
该研究的重要意义在于：
- 1.
提供了系统的代谢视角：首次通过整合实验验证的GEM与多状态转录组数据，系统揭示了具核梭杆菌在宿主相互作用过程中的动态代谢适应策略。
- 2.
揭示了潜在致病机制：研究预测的代谢变化，如甲酸盐和短链脂肪酸产量的改变，与具核梭杆菌促进结直肠癌进展的已知机制（如激活HIF1α信号、调节免疫微环境）相吻合，为理解其致病性提供了新的代谢层面解释。
- 3.
奠定了靶向治疗的基础：识别出的关键代谢脆弱点（如特定的氨基酸转运、发酵途径、一碳代谢酶），为未来开发针对具核梭杆菌的代谢干预策略（例如小分子抑制剂或营养干预）以预防或治疗相关疾病（特别是结直肠癌）提供了潜在靶点。
- 4.
构建了可扩展的研究平台：所构建的Fn-GEM是一个宝贵的资源，未来可进一步与宿主细胞代谢模型整合，用于研究宿主-微生物互作；或嵌入微生物群落模型，以探究其在复杂肠道菌群环境中的代谢互作。
总之，这项研究为深入理解具核梭杆菌这一重要病原体在疾病背景下的代谢灵活性提供了关键工具和见解，是连接细菌代谢与宿主病理生理学的重要桥梁。
