《Plant Biotechnology Journal》：Cell-Specific Expression and Cellular Compartmental Regulation in Camptothecin Biosynthesis
编辑推荐：
植物次生代谢产物，如单萜吲哚生物碱（MIAs），其生物合成和积累在特定细胞类型的特定细胞器中受到严格调控。然而，抗癌MIA喜树碱在茜草科植物喜树（Ophiorrhiza pumila）中的生物合成基因的细胞特异性表达模式，以及通过细胞区室化调控喜树碱生物合成的
植物次生代谢产物，如单萜吲哚生物碱（MIAs），其生物合成和积累在特定细胞类型的特定细胞器中受到严格调控。然而，抗癌MIA喜树碱在茜草科植物喜树（Ophiorrhiza pumila）中的生物合成基因的细胞特异性表达模式，以及通过细胞区室化调控喜树碱生物合成的机制仍知之甚少。本研究中，通过对喜树叶片进行单细胞转录组测序（scRNA-seq），获得了9181个细胞，利用标记基因注释和原位杂交将其分为四类。基因表达谱分析结合拟时间（pseudotime）分析表明，喜树碱生物合成基因优先在表皮细胞和幼嫩细胞中表达。值得注意的是，编码喜树碱关键酶——文朵灵合酶（strictosidine synthase, STR）和色氨酸脱羧酶（tryptophan decarboxylase, TDC）的基因在喜树的STR-TDC基因簇中发生了扩增，这支持了其细胞特异性功能分化。同时，该基因簇进化出了一个新功能基因OpAVT1，其被注释为液泡转运蛋白，负责介导色氨酸从液泡向细胞质的运输，从而调控喜树碱及其他MIAs的积累。总之，这些发现为喜树碱的生物合成与调控建立了细胞层面的框架。
**基于单细胞测序揭示喜树碱生物合成的细胞基础与区室化调控机制**
**研究背景与问题**
喜树碱是一种来源于喜树（Ophiorrhiza pumila）的单萜吲哚生物碱（monoterpenoid indole alkaloids, MIAs），其衍生物是拓扑异构酶I（topoisomerase I）的特异性抑制剂，在临床上具有重要的抗癌疗效。然而，喜树碱在源植物中的完整生物合成途径尚未完全阐明，特别是其中涉及的多个细胞色素P450（cytochrome P450, CYP450）酶的催化功能仍不清楚。此外，植物次生代谢产物通常表现出严格的细胞类型特异性积累模式，其生物合成和积累发生在特定的细胞和细胞器中。尽管喜树碱在多种组织中均有积累，但其生物合成在细胞层面如何被组织和调控，尤其是在哪种特定细胞中进行、该过程如何随细胞发育而动态变化、以及关键前体物质如何被调控以支持生物合成，这些问题尚不明确。因此，从细胞层面解析喜树碱的生物合成与调控机制，对于理解其生产规律和进行代谢工程改造至关重要。
研究人员为了在单细胞分辨率上系统解析喜树碱的生物合成与调控网络，利用单细胞转录组测序（scRNA-seq）技术对喜树叶片进行了深入研究。该研究揭示了喜树碱生物合成基因的细胞特异性表达模式和发育动态，鉴定了一个新的、位于核心代谢基因簇中的氨基酸液泡转运蛋白OpAVT1，并阐明了其通过调控色氨酸的跨液泡膜运输来协调喜树碱及其他MIAs生物合成的功能。
这项研究发表于《Plant Biotechnology Journal》，其核心意义在于首次为喜树碱的生物合成提供了一个精细的细胞图谱，并揭示了一个由基因簇编码的、集底物运输与催化功能于一体的细胞特异性代谢生产模块，为通过细胞工程提高喜树碱产量提供了新的靶点和理论依据。
**关键技术方法**
本研究采用的主要技术方法包括：1）单细胞转录组测序（scRNA-seq），利用10× Genomics平台对喜树叶片原生质体进行高通量测序，以获得细胞分辨的转录组数据；2）细胞类型鉴定与验证，通过计算聚类、同源标记基因注释以及原位杂交（in situ hybridisation, ISH）技术相结合，准确划分并验证了叶片的四种主要细胞类型；3）拟时间（pseudotime）分析，使用Monocle2软件重建表皮细胞的发育轨迹，以研究基因表达的动态变化；4）系统发育与共线性分析，对不同物种间的STR-TDC基因簇进行比较基因组学研究，以追溯其进化历程；5）功能验证，包括酵母互补实验验证OpAVT1的转运功能，以及利用CRISPR-Cas9技术构建基因敲除（knockout, KO）转基因毛状根系以进行体内功能分析；6）代谢组学与转录组学整合分析，通过液相色谱-串联质谱（LC-MS/MS）和qRT-PCR等技术，研究OpAVT1敲除对代谢物和基因表达的影响。研究样本主要来源于两个月龄的无菌喜树苗叶片。
**研究结果**
**2.1 单细胞RNA测序揭示喜树叶片细胞图谱**
研究人员从两个月龄的喜树叶片中分离原生质体，利用10× Genomics平台进行scRNA-seq，最终获得一个由9181个细胞、21448个基因构成的高质量数据集。通过主成分分析（principal component analysis, PCA）和无监督聚类，这些细胞被分为19个簇。结合拟南芥和长春花的已知叶片细胞标记基因进行同源比对，并利用原位杂交（ISH）技术进行验证，最终将这19个簇明确归类为四种主要的叶片细胞类型：表皮细胞（epidermal cells, ECs）、栅栏组织薄壁细胞（palisade mesophyll cells, PMCs）、海绵组织薄壁细胞（spongy mesophyll cells, SMCs）和维管细胞（vascular cells, VCs）。这一分析为后续研究提供了可靠的细胞分类基础。
**2.2 喜树碱生物合成在喜树叶片中的细胞类型特异性分布**
分析scRNA-seq数据发现，大部分已知的喜树碱生物合成基因在表皮细胞中特异性高表达。通过ISH实验对五个功能已验证的基因（OpG10H1、OpLAMT1、OpSLS1、OpTDC1和OpSTR1）进行验证，证实了它们确实在表皮细胞中定位。进一步利用多组学数据（包括不同细胞类型、组织和培养系的转录组数据）进行共表达网络分析，鉴定出31个可能调控喜树碱合成的候选转录因子，它们主要属于MYB、HD-ZIP和WRKY家族，且在表皮细胞、根和毛状根中高表达，与喜树碱的积累模式相吻合。
**2.3 拟时间分析筛选参与喜树碱生物合成的候选CYP450基因**
为了刻画喜树碱生物合成基因在表皮细胞发育过程中的动态表达，研究人员进行了拟时间分析。该分析将表皮细胞按照发育轨迹排序，并划分为7个状态。结果显示，包括OpG10H1、OpLAMT1、OpTDC1和OpSTR1在内的多个喜树碱生物合成基因在幼嫩的表皮细胞中高表达，其表达水平随着细胞成熟而下降。基于这一发育动态模式，结合全转录组数据分析了所有192个CYP450基因的表达谱，最终筛选出26个候选基因，它们的表达模式与已知的喜树碱合成基因相似。进一步的拟时间分析显示，其中10个CYP450基因在幼嫩表皮细胞中上调，表明它们可能参与喜树碱生物合成的早期步骤。
**2.4 STR-TDC基因簇的进化关系和功能验证**
基因组分析在喜树的10号染色体上发现了一个包含3个TDC和2个STR编码基因的代谢基因簇。通过跨物种比较基因组学和系统发育分析发现，虽然STR-TDC基因簇在多种MIA产生植物中保守存在，但基因的串联重复扩增事件仅发生在喜树中，产生了OpTDC3、OpTDC4、OpTDC5以及OpSTR1和OpSTR2。ISH验证和异源表达功能实验证实，这些重复的基因获得了细胞类型特异性的表达模式。例如，OpTDC3和OpTDC5特异地在表皮细胞中表达，而OpTDC4在表皮细胞、栅栏组织细胞和气孔中均有表达。所有三个OpTDC蛋白均能将L-色氨酸转化为色胺，OpSTR2也能催化严格苷（strictosidine）的形成，但其催化活性低于OpSTR1。这表明基因重复在保留原始催化功能的同时，驱动了功能的特化和分化。
**2.5 STR-TDC基因簇中的OpAVT1介导色氨酸从液泡向细胞质运输**
在进化过程中，该基因簇持续招募新基因以获得新功能。研究人员在喜树的STR-TDC基因簇中鉴定出一个保守的氨基酸液泡转运蛋白基因OpAVT1。系统发育分析显示它属于AVT1家族。ISH实验证实OpAVT1在表皮细胞中特异性高表达，与喜树碱合成基因的表达位置一致。亚细胞定位显示OpAVT1定位于液泡膜。在酵母互补实验中，单独表达OpAVT1不能支持酵母在以色氨酸为唯一氮源的培养基上生长，但与已知的色氨酸转运蛋白AtAAP6共表达时，会抑制酵母生长，表明OpAVT1在酵母中介导色氨酸的外排。进一步测试发现，OpAVT1还能外排苯丙氨酸、酪氨酸等多种氨基酸。通过CRISPR-Cas9技术构建的OpAVT1敲除毛状根系中，OpAVT1转录水平显著降低，根的生长表型发生改变，喜树碱含量显著下降，而游离色氨酸大量积累。这表明OpAVT1对于将液泡中储存的色氨酸转运至细胞质以供喜树碱合成至关重要。
**2.6 OpAVT1调控喜树碱及其他MIAs的生物合成**
代谢组学分析显示，在OpAVT1敲除系中，除了色氨酸积累和喜树碱减少外，其下游中间产物（如色胺、严格苷、喜树苷、喜果苷）以及其他多种经典MIAs（如毛钩藤碱、16-羟基它波宁、帽柱木碱、阿吗碱、阿枯米定和长春花碱）的水平均显著降低。同时，一些色氨酸衍生物和其它氨基酸相关的代谢物水平也发生了变化。转录组分析显示，敲除OpAVT1导致表皮细胞特异性表达的喜树碱合成基因（如OpTDC1、OpTDC3、OpTDC4和OpSTR1）表达下调。这些结果表明，OpAVT1通过两种机制调控喜树碱生物合成：一是直接调控色氨酸等前体物质的液泡-细胞质跨膜运输；二是影响关键生物合成基因的表达。
**讨论与结论**
本研究利用scRNA-seq技术，首次构建了喜树叶片的高分辨率细胞图谱，并从中解析了喜树碱生物合成的细胞特异性基础。研究发现，已知的喜树碱生物合成基因高度富集于表皮细胞，特别是幼嫩的表皮细胞中，这与长春花中MIA合成基因的多细胞分布模式不同。拟时间分析进一步揭示了这些基因在表皮细胞发育过程中的动态表达规律，并据此筛选出10个候选的、可能参与早期喜树碱合成的CYP450基因。这些发现建立了一个细胞层面的研究框架，有助于加速对喜树碱完整生物合成途径的解析。
本研究还深入表征了喜树中保守的STR-TDC代谢基因簇。该基因簇在喜树中发生了物种特异性的基因重复，产生了具有细胞类型特异性表达的TDC和STR基因家族成员，实现了功能的多样化。更重要的是，该基因簇进化出了一个关键的新成员——液泡转运蛋白OpAVT1。功能研究表明，OpAVT1定位于液泡膜，负责将储存在液泡中的色氨酸（以及其它氨基酸）转运至细胞质。在细胞质中，色氨酸被细胞质中的OpTDC酶转化为色胺，进而被OpSTR催化形成严格苷，后者是喜树碱合成的核心中间体。因此，OpAVT1、OpTDC和OpSTR在表皮细胞中形成了一个物理上紧密连锁、功能上高度协同的代谢模块，共同确保了喜树碱的高效合成。OpAVT1基因的敲除导致前体供应中断，不仅严重削弱了喜树碱的合成能力，也影响了其他次级代谢产物的生产，并短暂改变了植物的生长表型，这体现了初生代谢与次生代谢之间的平衡与权衡。
**研究结论**
综上所述，本研究通过单细胞转录组学揭示了喜树碱生物合成具有严格的表皮细胞特异性，并在发育过程中呈现动态变化。研究证实了喜树的STR-TDC基因簇是一个经历了谱系特异性扩张的功能性代谢基因簇，其编码的OpTDC和OpSTR酶具有细胞类型特异性的表达模式。研究鉴定并验证了该基因簇中一个新型液泡转运蛋白OpAVT1的功能，阐明了其通过调控液泡中色氨酸的输出来协调喜树碱生物合成的分子机制。这项工作建立了一个从细胞类型到基因簇功能、再到代谢物调控的完整模型，为理解植物次生代谢的细胞区室化调控提供了范例，并为通过细胞工程和基因编辑策略定向改良喜树碱生产奠定了理论基础。
