《Current Plant Biology》：Integrative Multi-Omics Reveals Genome Evolution and CYP719-Mediated BIA Biosynthesis in
Tinospora sagittata
编辑推荐：
本研究针对药用植物青牛胆（Tinospora sagittata）的基因组进化与苄基异喹啉类生物碱（BIA）生物合成机制展开。通过整合比较基因组学、转录组学、代谢组学等技术，研究人员解析了该植物基因组进化历史（发现约1.5百万年前的串联重复和约86.9百万年前的全基因组复制事件），并鉴定出CYP719基因候选TsA02G014550可催化(S)-canadine的形成，从而解析了jatrorrhizine的生物合成通路。该研究为深入理解植物BIA生物合成、基因组进化及后续可持续生产药用次生代谢产物提供了关键基础。
在传统中医药典籍中，青牛胆（Tinospora sagittata(Oliv.) Gagnep.）的干燥块根——“金果榄”（Radix Tinosporae）——因其清热、解毒、抗炎等功效而被广泛应用，并已收录于《中华人民共和国药典》2020年版。其药效主要源于两类重要的次生代谢产物：二萜类化合物和苄基异喹啉类生物碱（Benzylisoquinoline Alkaloids, BIAs）。后者包含多种高活性的化合物，如小檗碱、巴马汀、药根碱（jatrorrhizine）等。药根碱作为青牛胆中的主要活性BIA成分，已被研究证实具有抗糖尿病、抗菌、抗癌、抗肥胖、降血脂及有益于中枢神经系统等多种药理活性，显示出巨大的医疗应用潜力。
然而，尽管其药用价值显著，我们对青牛胆的认知仍存在巨大的空白。首先，其基因组资源的匮乏严重阻碍了对该物种在进化树上的精准定位及其与近缘物种关系的理解。其次，关于BIA，尤其是药根碱的生物合成途径，仍有许多关键酶和调控步骤尚未阐明。这种“知其然而不知其所以然”的困境，限制了科学家利用现代生物技术手段，对青牛胆进行定向遗传改良或通过合成生物学方法大规模、可持续地生产这些珍稀药用成分。为此，一个系统性的、整合多组学技术的研究势在必行。
为了填补这些知识空白，一个来自华中农业大学的研究团队开展了一项雄心勃勃的综合性研究。他们利用先进的比较基因组学、转录组学、代谢组学（靶向与非靶向）技术，对青牛胆的基因组进化、基因家族扩张以及BIA生物合成通路进行了深入解析。该研究成果以“Integrative Multi-Omics Reveals Genome Evolution and CYP719-Mediated BIA Biosynthesis in Tinospora sagittata”为题，发表在《Current Plant Biology》期刊上，为我们揭开了这种神奇药用植物生命蓝图的一角。
为了开展这项研究，研究人员采用了多项关键技术方法。首先，他们从NCBI数据库下载了青牛胆及其它16种植物的基因组数据进行比较基因组学分析。研究材料来自中国湖北省利川市的青牛胆种子，并在华中农业大学研究站的受控环境中种植。通过转录组学分析，他们对青牛胆的嫩叶、茎、根和块茎等不同组织进行了RNA测序，以解析基因表达模式。在代谢组学方面，团队结合了靶向和非靶向代谢物分析，对青牛胆中的BIA等次生代谢物进行了定性和定量分析，特别是利用高效液相色谱（HPLC）对药根碱含量进行了测定。此外，通过同源搜索、系统发育分析和基因簇挖掘，鉴定出与BIA生物合成相关的候选基因。最后，利用酵母（Saccharomyces cerevisiae WAT11菌株）异源表达系统和体外酶活测定，对关键候选基因TsA02G014550进行了功能验证。
3.1. 青牛胆与16个其他物种的基因组比较
研究人员将青牛胆的基因组与其他16种植物的基因组进行了系统性比较。分析发现，在所有研究的物种中，青牛胆与黄连（Coptis chinensis）的亲缘关系最近，两者共享了181个基因家族。共线性分析进一步证实了青牛胆与黄连、罂粟（Papaver somniferum）等物种之间存在广泛的同源基因关系，为后续的进化分析奠定了基础。
3.2. 独特基因家族的功能分类
对青牛胆独有的766个基因家族进行功能富集分析（GO和KEGG）发现，这些家族与RNA-DNA杂交核糖核酸酶活性、转录调控区DNA结合等分子功能相关，并在苯并噁嗪类生物合成、ABC转运蛋白、单萜类生物合成等通路中富集。这表明青牛胆在次生代谢调控方面可能具有独特的遗传基础。
3.3. 基因家族扩张、系统发育特征与全基因组复制
基因家族进化分析显示，在青牛胆基因组中，有186个基因家族发生了扩张，6个发生了收缩。扩张的家族功能富集于DNA解旋酶活性、UDP-糖基转移酶活性等，并且在苯并噁嗪类生物合成、2-氧代羧酸代谢以及异喹啉生物碱生物合成（isoquinoline alkaloid biosynthesis）通路中显著富集，其中就包含了与BIA生物合成密切相关的细胞色素P450（CYP）和BIA桥酶（BBE）基因家族。
系统发育树构建和分化时间估算结果表明，青牛胆与黄连大约在89-109百万年前（白垩纪时期）从一个共同祖先分化出来。更重要的是，通过分析同义替换率（Ks）和四倍简并位点置换（4DTv）的分布，研究人员发现了青牛胆基因组经历两次重大复制事件的证据：一次是约1.5百万年前发生的串联重复（tandem duplication），另一次是约86.9百万年前发生的古全基因组复制（ancient Whole-Genome Duplication, WGD）。与葡萄（Vitis vinifera）基因组的共线性深度比例（2:3）分析表明，青牛胆经历了一次与核心真双子叶植物六倍化事件相关的古老WGD。
进一步分析显示，BIA生物合成通路中的多个基因（如PAT, TAT, CAS, BBE等）在青牛胆基因组中扩张为四个拷贝，且这些复制事件与近期的WGD事件相关。微共线性分析也证实，BIA生物合成基因在青牛胆、黄连和罂粟的基因组中是保守的。
3.4. 青牛胆中次级代谢物基因簇及BIA生物合成通路基因的确定
利用PLANTiSMASH软件，研究人员在青牛胆基因组中预测了58个次级代谢物基因簇，其中包含7个与生物碱合成相关的簇。通过差异基因表达和共表达网络分析，筛选出7个与BIA生物合成通路密切相关的基因，其中包括位于A02染色体上的TsA02G014550，该基因被注释为(S)-stylopine/(S)-canadine/(S)-nandinine合成酶，属于CYP719家族。
3.5. 药根碱生物合成关键候选基因的鉴定
通过整合转录组和代谢组数据，研究人员发现54个差异表达基因和8个差异积累代谢物在BIA生物合成通路中富集。在青牛胆基因组中鉴定出543个CYP450基因，其中TsA02G014550与来自罂粟、日本黄连（C. japonica）和黄连的已知(S)-canadine合成酶（CAS）基因聚为一支，强烈暗示其参与药根碱的生物合成。
非靶向代谢组学在青牛胆中注释了124种生物碱，包括32种异喹啉生物碱。靶向HPLC-MS分析则明确显示，药根碱在叶片中的含量最高。
3.6. TsA02G014550的体外催化分析
基因表达分析（qRT-PCR和RNA-seq）均表明TsA02G014550在叶片中表达量最高，与药根碱的积累模式一致。为了验证其功能，研究人员将该基因在酵母中进行异源表达，并提取微粒体蛋白进行体外酶活测定。结果证实，以(S)-四氢非洲防己碱为底物时，TsA02G014550编码的重组蛋白能够催化生成(S)-canadine，从而明确了其作为(S)-canadine合成酶（CAS）的功能。
本研究通过多组学整合分析，系统解析了青牛胆的基因组进化历史与BIA生物合成机制，得出了一系列重要结论。
首先，在基因组进化层面，研究证实青牛胆经历了一次古老的WGD事件（约86.9百万年前）和一次近期的串联重复事件（约1.5百万年前）。这次古老的WGD发生在毛茛目（Ranunculales）内部、防己科（Menispermaceae）与罂粟科（Papaveraceae）、毛茛科（Ranunculaceae）分化之后，与黄连、罂粟中发现的WGD事件相似，共同促进了BIA生物合成相关基因家族的扩张与功能分化。共线性分析支持了青牛胆A、B亚基因组之间的同源关系，并揭示了B亚基因组中存在大量的结构变异（如易位、倒位、拷贝数变异等），这些变异可能影响了淀粉与蔗糖代谢、蛋白质加工、异喹啉生物碱生物合成等关键通路。
其次，在BIA生物合成层面，研究成功挖掘出青牛胆基因组中与BIA合成相关的多个基因家族，并发现其中CYP450和BBE等基因家族显著扩张。特别重要的是，通过序列挖掘、系统发育分析和多组学关联，鉴定出CYP719家族成员TsA02G014550为关键候选基因。随后的体外功能实验确凿无疑地证明，TsA02G014550编码的蛋白具有(S)-canadine合成酶活性，能够催化(S)-四氢非洲防己碱转化为(S)-canadine，这是药根碱生物合成通路中的关键步骤。
这项研究的意义深远。它不仅首次在染色体水平上解析了重要药用植物青牛胆的基因组，厘清了其在毛茛目中的进化地位和分化时间，为植物比较基因组学和进化生物学提供了宝贵数据；更重要的是，它通过整合多组学策略，成功地从复杂的基因组中定位并验证了控制关键药用成分合成的“靶点”基因。这为未来通过基因工程或合成生物学手段，在微生物或模式植物中重建并优化BIA（尤其是药根碱）的生产途径奠定了坚实的遗传基础，有望实现这些高价值药用化合物的可持续、规模化生产，从而降低对野生植物资源的依赖，具有重要的科学价值和应用前景。
