植物相关鞘氨醇单胞菌属（Sphingobium）谱系对喜树碱（CPT）的转化与质粒附属区域GAP2的关联性研究

《Current Research in Microbial Sciences》：Camptothecin transformation in plant-associated Sphingobium lineages is linked to a plasmid accessory region

　　

植物次生代谢物在植物相关微生物栖息地中扮演着“化学通货”的核心角色。这些次生代谢物可作为防御分子、化感物质或信号分子，其在土壤和根际环境中的持久性不仅受非生物过程影响，也受微生物转化和降解作用影响。微生物对植物来源芳香族化合物的代谢已被广泛研究，尤其是相对简单的酚类化合物和类外源物质底物。相比之下，关于细菌如何代谢结构复杂的植物生物碱，目前所知甚少。这一知识空白在微生物生态学中具有重要意义，因为生物碱可显著影响微生物群落结构；同时该信息缺失在生物技术领域同样重要，因为微生物酶可在温和条件下选择性转化复杂天然产物。 喜树碱（Camptothecin，CPT）是一种五环喹啉类生物碱，是临床使用的拓扑异构酶I抑制剂（如伊立替康和拓扑替康）的先导骨架。然而，CPT也是一种特化代谢物，存在于产CPT植物（包括Ophiorrhiza pumila Champ. ex Benth.）的根相关组织中。此类植物释放的CPT或相关化合物可使根际微生物暴露于具有生物活性的含氮芳香骨架。这一生态背景提出了一个独特的微生物学问题：植物相关细菌能否转化或利用CPT？如果能，哪些生物体、代谢物和基因介导了该生物碱进入细菌代谢途径？ 关于CPT的细菌代谢细节在很大程度上仍未解决。本研究中，研究人员从植物相关环境中分离出转化CPT的Sphingobium菌株，并整合了表型测定、代谢物化学、基于基因组的分类学、比较基因组学、跨谱系保守性以及RNA测序，以界定与CPT转化相关的候选基因组框架。从O. pumila根际分离的菌株OP1将CPT转化为一种瞬时产物M1，该产物保留了E环内酯并在孵育过程中随后减少。当CPT作为唯一添加碳源时，OP1的活菌数也增加。基于基因组的分类学分析将OP1、NN19和KJ09归入一个推定的新物种级谱系，并将IROP2归入一个相关但不同的谱系。将NN19与不转化CPT的NN19衍生突变体NTG18进行比较，鉴定出GAP2——一个富含加氧、开环、氧化还原、水解酶、转运和调控候选基因的质粒相关区域。整合缺失作图、跨谱系保守性、功能注释以及在CPT生长期间的表达数据，优先确定了GAP2内的一个紧凑加氧/开环模块。研究人员提出一个可检验的模型：质粒相关的附属代谢使植物相关Sphingobium谱系能够转化结构复杂的植物生物碱，并可能有助于其在化学选择性植物相关栖息地中的适应。

论文解读文章

一、研究背景与科学问题

植物次生代谢物在植物与微生物互作中扮演“化学通货”的重要角色，可作为防御分子、化感物质或信号分子影响根际微生物群落结构。微生物对植物来源芳香族化合物的代谢已有大量研究，但关于细菌如何代谢结构复杂的植物生物碱，学界认知极为有限。喜树碱（Camptothecin，CPT）是一种五环喹啉类生物碱，是临床抗肿瘤药物伊立替康和拓扑替康的先导骨架，以往研究主要集中于其药理学、生物合成及生物技术应用。然而，CPT也是产CPT植物根相关组织中的特化代谢物，其释放可使根际微生物暴露于具有生物活性的含氮芳香骨架。此前，CPT在植物相关环境中的微生物转化机制、参与转化的菌株及遗传决定因子均未被阐明。环境研究虽考察了CPT的降解，但海水中的生物降解据报道可忽略不计。因此，细菌如何代谢CPT这一科学问题长期悬而未决。Sphingobium属细菌以降解芳香族或难降解有机化合物著称，其降解功能常由质粒相关的附属复制子编码，这使其成为研究植物生物碱微生物代谢的理想模式类群。

二、研究内容与总体结论

研究人员从产CPT植物Ophiorrhiza pumila根际及日本西表岛相关土壤中分离获得多株可转化CPT的Sphingobium菌株（OP1、NN19、KJ09、IROP2），并通过表型测定、代谢物化学分析、全基因组分类学、比较基因组学、跨谱系保守性分析及转录组测序，系统解析了CPT细菌转化的遗传基础。研究发现：菌株OP1可将CPT转化为保留E环内酯的瞬时产物M1，且以CPT为唯一添加碳源时活菌数显著增加；通过比较CPT转化菌株NN19与其不转化突变体NTG18，鉴定出位于质粒pNN19-1上的97.8-kb附属区域GAP2；GAP2内包含一个由加氧酶和开环酶基因组成的核心候选模块，该模块在另一独立转化谱系IROP2中高度保守，且在CPT生长条件下高表达。研究提出一个可检验的模型：质粒相关的加氧/开环模块使植物相关Sphingobium谱系能够进入结构复杂植物生物碱的代谢途径。该研究首次将CPT的环境归趋与特定细菌谱系及其候选遗传决定因子联系起来，具有重要的化学生态学和生物催化意义。

三、主要关键技术方法

研究人员采用了以下核心技术手段：从植物相关土壤进行CPT富集培养并分离转化菌株；利用365-nm紫外荧光筛选法检测CPT转化活性；采用高效液相色谱（HPLC）、液相色谱-电喷雾电离-飞行时间质谱（LC-ESI-TOF-MS）及核磁共振（NMR）波谱进行代谢产物M1的结构鉴定；利用N-甲基-N'-硝基-N-亚硝基胍（NTG）诱变获得不转化突变体NTG18；采用PacBio HiFi长读长测序和DNBSEQ短读长测序进行全基因组组装，并通过DFAST进行基因组注释；利用Type (Strain) Genome Server（TYGS）进行基于全基因组的分类学分析；采用MUMmer4 dnadiff进行成对基因组比对以界定GAP2区域；利用BLASTp进行跨谱系同源蛋白搜索；采用RNA-seq（HISAT2比对和featureCounts计数）评估候选基因在CPT生长条件下的表达水平。

四、研究结果

4.1 菌株OP1的分离与CPT转化特性鉴定

研究人员从O. pumila根际土壤经富集培养分离获得菌株OP1，该菌在365-nm紫外照射下表现出CPT荧光局部减弱，表明其具有CPT消耗或转化能力。16S rRNA基因序列初步将其归入Sphingobium属。在以CPT为唯一添加碳源的无机盐培养基中，OP1活菌数从0 h的5.3 ± 0.7 log₁₀ CFU/mL增至48 h的8.6 ± 0.1 log₁₀ CFU/mL，而对照S. japonicum、S. chlorophenolicum和S. yanoikuyae几乎不消耗CPT。培养体系pH基本不变（7.31–7.33），排除了pH依赖性水解的干扰。

4.2 基于基因组的分类学分析

TYGS分析将OP1、NN19和KJ09归入同一物种级聚类（dDDH值为100.0%），IROP2则构成一个相关但不同的姐妹谱系（与前三者的dDDH值为57.0%，低于70%物种界限）。OP1、KJ09和NN19的基因组大小分别为3,054,081 bp、3,054,235 bp和3,060,055 bp，G+C含量均为62.0%左右；IROP2基因组较小（2,858,869 bp，62.2% G+C）。这些菌株代表两个相关的推定新物种级Sphingobium谱系。

4.3 OP1将CPT转化为瞬时产物M1

静息细胞实验证实OP1可主动转化CPT，且该活性不依赖于预先的CPT暴露。HPLC分析显示CPT峰下降伴随新峰M1出现，细胞-free对照中无M1生成。高分辨质谱确定M1分子式为C₂₀H₁₆N₂O₆（[M+H]⁺实测m/z 381.1107），NMR分析表明M1保留了CPT的E环内酯。时间过程HPLC显示M1在反应早期瞬时积累后随孵育时间延长而减少，提示其被进一步转化而非稳定终产物。

4.4 CPT类似物的底物选择性转化

OP1可转化10-羟基喜树碱（HCPT），但对7-乙基-10-羟基喜树碱（EHCPT）的转化延迟或受限。LC-ESI-TOF-MS检测到EHCPT的单加氧产物（C₂₂H₂₁N₂O₆ [M+H]⁺，m/z 409.14136）和双加氧产物（C₂₂H₂₁N₂O₇ [M+H]⁺，m/z 425.13287），表明OP1介导的转化对CPT骨架的特定结构特征具有识别性，加氧是早期转化步骤。

4.5 GAP2候选区域的界定

通过NN19与其不转化突变体NTG18的全基因组比对，鉴定出NTG18缺失了pNN19-1上的两个区域：R2（10,083 bp，位置401,277–411,359）和R3（86,997 bp，位置411,463–498,459）。该缺失区间（pNN19-1位置401,171–498,940，共97.8-kb）被命名为GAP2。HiFi读长验证确认R2和R3缺失，同时NTG18存在23个高置信度单核苷酸变异（SNVs），其中19个为NTG18特异，但未涉及明显的CPT加氧或开环候选蛋白。

4.6 GAP2基因编码加氧和开环候选酶

GAP2内鉴定出多个与芳香族氧化、开环、水解、氧化还原和转运相关的基因。其中四个基因被指定为核心候选模块：LOCUS_040（NN19_28260，编码外二醇双加氧酶）、LOCUS_070（NN19_28290，编码NtaA/DmoA家族FMN依赖性单加氧酶）、LOCUS_080（NN19_28300，编码芳香环羟化双加氧酶α亚基）和LOCUS_090（NN19_28310，编码外二醇型假定蛋白）。该组合与代谢物数据推断的化学过程一致：先对CPT芳香/多环骨架进行初始加氧，随后进行开环型加工。

4.7 GAP2位于更广泛的质粒附属区域背景中

GAP2区间含77个编码序列，功能分类显示29个基因（37.7%）与氧化还原功能相关，7个（9.1%）与水解酶相关，10个（13.0%）与转运或受体相关，6个（7.8%）与转录调控相关。R2含6个编码序列，R3含71个编码序列。该组成支持将GAP2解释为一个类分解代谢岛附属区域，而非单一酶基因座。

4.8 IROP2携带相关的GAP2样候选模块

BLASTp分析显示，IROP2蛋白组中存在12个功能筛选的GAP2候选蛋白中11个的同源物，包括核心候选LOCUS_040（97.5%一致性）、LOCUS_070（94.7%）、LOCUS_080（98.8%）和LOCUS_090（97.0%）。这些同源基因位于pIROP2的两个簇中，表明IROP2携带一个相关但重排的质粒模块。IROP2也能从HCPT和EHCPT产生加氧产物，为其CPT类似物加氧能力提供了表型支持。

4.9 RNA-seq表达分析优先确定候选基因

RNA-seq数据显示，在CPT生长条件下，LOCUS_080表达量最高（平均TPM 9316.9，全基因组表达排名第13位），LOCUS_070平均TPM 1428.5（排名113），LOCUS_040平均TPM 1394.8（排名116），LOCUS_090平均TPM 982.6（排名166）。相比之下，编码酰胺水解酶家族蛋白的LOCUS_060转录丰度较低。表达模式突出显示了GAP2内的紧凑早期氧化/开环模块。

4.10 细菌CPT转化的工作假说

综合代谢物、NTG18缺失、IROP2保守性、功能注释和RNA-seq数据，研究人员提出：CPT代谢的起始步骤是CPT骨架的加氧，随后进行开环。LOCUS_080是早期加氧步骤的主要候选基因；LOCUS_070可能在平行或顺序途径中介导加氧；LOCUS_040和LOCUS_090可能作用于羟化或儿茶酚类CPT中间体的开环反应；LOCUS_050和LOCUS_060等基因可能参与下游加工。

五、讨论与结论

讨论部分指出，本研究将CPT代谢置于植物-根际微生物功能的背景下。OP1从产CPT植物根际的分离与其转化植物特化代谢物的能力一致，M1的瞬时积累和消失表明OP1进一步转化M1而非将其作为稳定终产物积累。CPT转化菌株代表两个此前未描述的Sphingobium物种级谱系，扩展了CPT微生物学的研究范畴。GAP2作为质粒相关附属区域，其基因组成与芳香族化合物降解途径的典型架构一致。四个独立证据线汇聚支持GAP2模块为CPT转化的主要候选基因组区域：不转化突变体NTG18缺失GAP2；独立转化谱系IROP2携带相关pIROP2模块；核心加氧/开环基因在CPT生长期间高表达；预测酶功能与观察到的产物形成和转化一致。

研究结论指出：本研究鉴定出转化CPT的植物相关Sphingobium谱系，并将该表型与质粒相关GAP2区域联系起来。菌株OP1形成保留E环内酯的瞬时产物M1，且以CPT为唯一添加碳源时活菌数增加。GAP2在不转化突变体NTG18中的缺失、pIROP2上重排同源模块的保守性以及优先基因在CPT生长期间的转录，共同指向一个紧凑的加氧/开环候选模块。这些发现支持一个可检验的模型：质粒相关的附属代谢使细菌能够进入结构复杂植物生物碱的代谢途径，并可能有助于其在化学选择性植物相关栖息地中的适应。优先确定的基因为后续遗传和生化验证提供了聚焦靶点。