黏细菌甾体天然产物的基因组学与代谢组学整合发掘：从分散式生物合成基因架构到新型开环甾体发现

《Applied and Environmental Microbiology》：Integrated genomics and metabolomics survey of steroids from myxobacteria

　　

摘要专业翻译 与其它天然产物类别相比，来自黏细菌的甾体化合物在很大程度上仍未得到充分探索，这导致研究人员对其生物合成多样性、进化轨迹及药理潜力的理解存在显著空白。为弥补这一不足，研究人员对黏细菌基因组开展了系统的基因组学调查，以描绘甾醇生物合成基因的分布与系统发育背景，并发现这些基因以分散方式排列于染色体上，而非组织成典型的生物合成基因簇（BGCs）。尽管甾醇生物合成机制在黏细菌门（Myxococcota）中广泛分布，但其在不同谱系间的分布极不均衡，其中侏囊菌科（Nannocystaceae）展现出最高的生物合成潜力。对甾醇生物合成关键酶——氧化鲨烯环化酶（OSCs）——的序列相似性网络（SSN）挖掘揭示了广泛的序列分化和尚未探索的新型催化功能空间。受上述生物信息学预测的驱动，研究人员进一步证明，将色谱预分离与基于核磁共振（NMR）的代谢谱分析相结合，可显著简化甾体化合物的靶向分离流程。该方法从2株实验室自有黏细菌菌株中获得了12个（开环）甾体化合物，其中包括4个此前未描述的结构，这些结构展现出多样化的氧化修饰、骨架重排及环降解特征。生物活性测定揭示，共产生的代谢产物烷基糖苷可协同增强特定开环甾体的抗菌活性。综上，本研究建立了一条基因组学-代谢组学整合研究流程，不仅扩展了黏细菌甾体的结构库，也为未来隐秘甾醇代谢产物及其隐藏生物学功能的发现提供了蓝图。

论文解读文章

一、研究背景与科学问题

甾体化合物（steroids）是自然界中结构多样性最为丰富的天然产物类别之一，其核心骨架为环戊烷多氢菲四环体系，可通过氧化模式变化、双键位置调整、侧链结构修饰及骨架重排等方式实现广泛的结构多样化。甾体在真核生物中承担膜组分等基础生理功能，同时展现出广泛的药理活性，是药物发现中极具价值的优势骨架。从生态学视角看，甾体还可作为化学防御物质、共生信号分子以及环境胁迫下的膜适应调节因子。

黏细菌（Myxococcota门）是一类以复杂社会行为和卓越天然产物产生能力著称的革兰氏阴性菌。尽管黏细菌已贡献了埃博霉素（epothilones）等临床相关代谢产物，但其甾体生物合成长期以来仅见于零散报道，缺乏系统性认识。1983年Kohl等人首次从Nannocystis exedens中报道了甾体化合物；Müller研究组揭示了Stigmatella aurantiaca产生植物型甾醇环阿屯醇（cycloartenol），并克隆了首个细菌环阿屯醇合酶基因；Welander研究组通过系统发育基因组学、异源表达和脂质组学系统阐明了细菌甾醇生物合成。然而，与其它天然产物类别相比，黏细菌甾体的生物合成分布、结构多样性和生物学活性仍相对未被充分探索。

本研究旨在系统调查黏细菌门中甾醇生物合成基因的分布与系统发育背景，揭示其基因组组织特征，并在此基础上建立高效的甾体代谢产物靶向发现流程，以扩展黏细菌甾体的化学空间并理解其生物学功能。

二、主要技术方法

研究人员整合运用了以下关键技术方法开展研究：对369个高质量黏细菌完整基因组（其中280个来自NCBI公共数据库，89个为实验室自有测序基因组）进行生物信息学分析，基于16S rRNA序列构建贝叶斯系统发育树，并通过BLAST与profile hidden Markov model（profile HMM）双重策略筛选甾醇生物合成相关基因；利用EFI-EST网络资源构建氧化鲨烯环化酶（OSC）的序列相似性网络（SSN）；采用C₁₈反相色谱柱以20%梯度甲醇-水进行预分离，结合核磁共振（NMR）代谢谱分析实现甾体化合物的靶向检测；通过半制备高效液相色谱（RP-HPLC）分离纯化目标化合物，并综合运用高分辨质谱和二维NMR进行结构鉴定；采用CCK-8法评价抗增殖活性，采用肉汤微量稀释法评价抗菌活性，并通过棋盘稀释法（checkerboard assay）定量评估药物协同效应。

三、研究结果

黏细菌门中甾醇生物合成基因的系统发育分布

研究人员对369个黏细菌基因组进行甾醇生物合成潜力分析，构建的16S rRNA贝氏系统发育树实现了稳健的科级分类，产生五大主要分支：Myxococcaceae、Archangiaceae、Anaeromyxobacteraceae、Nannocystaceae和Polyangiaceae。基于甾醇生物合成关键基因存在的加权评分系统显示，甾醇生物合成机制在Myxococcota门内广泛分布，但各谱系间分布极不均衡。生物合成能力主要体现于Myxococcaceae（如Myxococcus和Corallococcus属）、Nannocystaceae（如Nannocystis和Enhygromyxa属）、Archangiaceae（如Archangium、Cystobacter和Stigmatella属）及Polyangiaceae（如Sorangium、Minicystis和Polyangium属）。其中Nannocystaceae科尤为突出，通常携带近乎完整的甾醇生物合成途径，提示甾醇生物合成可能在这些科中具有祖先性。值得注意的是，黏细菌中甾醇生物合成基因超过90%以分散方式排列于染色体上，而非集中于单一连续基因簇内，这与典型的细菌次级代谢产物BGCs（如PKS/NRPS基因簇）形成鲜明对比。

黏细菌OSC的序列相似性网络挖掘

通过BLAST与profile HMM联合搜索策略，在严格阈值（e值≤1×10^-10，一致性≥30%，覆盖率≥70%）下，从369个黏细菌基因组中鉴定出190个OSC阳性候选基因组。OSC分布统计显示Myxococcaceae科占比最高（52.0%），部分归因于该科测序菌株在数据集中的过度代表；Nannocystaceae科尽管代表性较低（11%），却展现出最高的甾醇生物合成潜力。基于EFI-EST构建的SSN网络（比对分数阈值360）显示，OSC聚类大体上重现了基于16S rRNA的黏细菌系统发育，但两者并非完全重叠，提示不同黏细菌谱系在OSC序列上经历了不同的进化轨迹。最大的OSC簇包含来自甾体丰富属的序列，所有参考OSC也位于该主要网络内。此外，从Archangium和Sandaracinus属中鉴定出多个与主要网络明显分离的未表征单例OSC，这些序列代表了通过异源表达发现新型甾体的首要候选。

NMR代谢组学结合色谱预分离实现甾体靶向分离

基于生物信息学预测，研究人员选取Myxococcus fulvus 09-6s3和Corallococcus interemptor SDUN13两株菌株开展靶向化学研究。将粗提物经C₁₈反相色谱柱以20%梯度甲醇-水预分离后，100%甲醇洗脱组分展现出甾体侧链双键的特征烯烃质子信号（δ_H 5.1–5.5）以及甾体骨架典型的高场区甲基共振信号（δ_H 0.6–1.2），这些特征在未分级的粗提物谱图中明显缺失或被严重掩盖。茴香醛-硫酸试剂显色的薄层色谱（TLC）分析显示特征性的紫色至蓝紫色斑点，指示甾体化合物存在。该策略最终从两株菌中分离获得12个（开环）甾体化合物1–12，其中化合物5–7和12为4个此前未描述的新结构。化合物1–7具有基本四环甾体骨架并含多种含氧官能团；化合物10–12为高度降解的甾体，具有七降麦角甾烷（heptanorergosterane）骨架和三环结构。

黏细菌甾体的生物合成途径

基于C. interemptor SDUN13和M. fulvus 09-6s3的基因组测序，研究人员提出了重排和未重排黏细菌甾体的生物合成模型。两株菌中均存在HMG-CoA合酶、甲羟戊酸激酶、甲羟戊酸-5-焦磷酸脱羧酶和IPP异构酶同源基因，提示DMAPP和IPP通过甲羟戊酸（MVA）途径生成。SQS同源基因的存在表明鲨烯产生在生化上可行。SQMO与OSC的协同作用建立基本四环三萜骨架（如羊毛甾醇），随后经CYP51和SdmA?C催化的C-14和C-4去甲基化产生中间体酵母甾醇。EBP、ERG3、DHCR24和SMT基因同源物的存在提示麦角甾醇可能作为化合物1–12的基本前体。麦角甾醇可经侧链和/或四环的连续氧化而不发生骨架重排，引入多种氧化修饰；另一方面，四环骨架可经氧化开环后发生重排产生多样化骨架。例如，B环C4-C5和C9-C10键的氧化断裂产生含共轭烯酮体系的开环甾体中间体，伴随A环的降解丢失，随后分子内内酯化形成三环incisterol型骨架（如化合物12）；C环C9-C11键的氧化断裂后经乙烯基醇醛加成产生重排的corasterol（化合物9）；D环C14-C15键的氧化断裂后与侧链烯烃发生乙烯基醇醛加成产生化合物8。该生物合成架构在SDUN13和09-6s3菌株中同样呈分散式分布，与植物次级代谢的组织特征高度相似。

生物活性测定

化合物1–7和9–11对五种人癌细胞系（A549、HepG2、MCF7、HCT116和THP-1）的抗增殖活性评价显示，除骨架重排的corasterol（9）对人单核细胞白血病THP-1细胞表现出中等活性（IC₅₀为22 μM）外，其余化合物均无活性（IC₅₀ > 50 μM）。抗菌活性测试中，仅开环甾体10和11对金黄色葡萄球菌ATCC 12600、白色念珠菌5312和尖孢镰刀菌P4表现出弱抑制活性（MIC为32–64 μg/mL），其余甾体在128 μg/mL浓度下均无活性。这一明确的结构-活性关系表明，甾体核心的氧化断裂和重排对于将普通惰性骨架转化为抗生素至关重要。进一步研究发现，M. fulvus 09-6s3产生的烷基糖苷（128 μg/mL）与化合物11联用时，可使化合物11对上述三种微生物的MIC降低达32倍。棋盘稀释实验定量验证了烷基糖苷与化合物11对尖孢镰刀菌P4的协同作用，FIC指数<0.5，表明产生菌内部可能存在复杂的化学生态关系，即特化的佐剂分子可能协同增强共产生甾体的活性。

四、讨论与结论

本研究最显著的发现是黏细菌甾醇生物合成基因 overwhelmingly 以分散方式排列于染色体上，超过90%的甾醇基因集不集中于单一基因座。研究人员提出这种分散架构并非简单的进化遗迹，而是主动促进结构多样化的功能创新。与典型BGC中酶物理共定位和协同调控不同，核心骨架形成酶（OSC、CYP51）与外围修饰氧化酶的脱偶联使各自能够独立进化：核心酶维持强纯化选择以保障基本骨架组装，而分散的修饰基因可通过复制、新功能化和水平转移自由多样化，而不干扰初级甾醇产生。这种"分裂途径"组织方式与植物次级代谢相似，可能解释了从简单羟基化到完全B环降解和乙烯基醇醛重排的显著结构多样性。

鉴于甾醇生物合成基因在黏细菌基因组中的分散组织，依赖连续BGC检测的传统基因组挖掘策略本质上难以预测结构新颖性。相比之下，NMR代谢组学恰好擅长于序列方法失效之处：它能够捕获分散修饰酶产生的全部氧化和骨架多样化程度，包括无法从碎片化基因组背景推断生物合成逻辑的开环甾体和重排甾体。甾体代谢产物具有特征性且易于识别的波谱特征，可作为甾体骨架的可靠生物标志物，使甾体代谢产物能够在无需任何基因组注释的情况下从高度复杂的粗提物中被标记识别。与色谱预分离整合后，NMR代谢组学的分辨能力显著增强。

黏细菌为甾体药物的化学酶法生物催化提供了有前景的生物技术解决方案。本研究揭示了黏细菌卓越的多步氧化修饰能力，包括骨架重排、环断裂和位点特异性羟基化。黏细菌细胞色素P450酶（如CYP260A1和CYP260B1）和非经典氧化鲨烯环化酶已展现出显著的底物混杂性和立体选择性，可作为药物化学酶法合成中后期C–H官能团化的强大生物催化剂。结合基因组调查揭示的大量未表征修饰酶库，黏细菌成为双重用途资源：既是发现具有前所未有生物活性的新型甾体化学型的沃土，也是工程化下一代甾体治疗药物的丰富酶工具库。

本研究存在以下局限性：369个基因组仅代表黏细菌多样性的一小部分，甾醇生物合成潜力的真实范围可能更大；基因组预测鉴定出具有甾醇生物合成能力的菌株并不保证在实验室培养条件下产生甾醇，因为途径表达可能受环境线索严格调控；所提出的重排和降解甾体生物合成模型仍属假设，需要通过候选加氧酶的异源表达和体外重构获得直接酶学证据加以验证。未来研究应优先关注三个方向：通过异源表达功能表征SSN挖掘中鉴定的未表征OSC；通过产生菌与非产生菌的比较基因组学鉴定负责骨架重排和环断裂的修饰加氧酶并进行生化验证；通过基因敲除研究和膜生物物理学研究黏细菌甾体的生态和生理功能。