Clerodendrum trichotomum（海州常山）染色体水平基因组组装揭示其复杂全基因组复制历史

《Plant Diversity》：Chromosome-level genome assembly of Clerodendrum trichotomum sheds light on its complex whole genome duplication history

　　

Clerodendrum trichotomum Thunb.（海州常山）是唇形科（Lamiaceae）大青属（Clerodendrum）的一种灌木或小乔木，广泛分布于中国、韩国、日本及菲律宾北部。该物种兼具观赏价值，其叶片和茎被揉碎后释放的特殊不愉快气味是众多大青属植物的典型特征。除园艺用途外，海州常山长期被用作治疗高血压和炎症的传统中药，其嫩叶在中国西南地区作为野生蔬菜食用。植物化学研究表明，该植物富含萜类、黄酮类、苯丙素类、苯乙醇苷类、酚苷类及生物碱等次生代谢产物，其中萜类最为丰富。然而，该物种的挥发性萜类化合物尚未得到系统分析。为此，研究人员组装了高质量染色体水平基因组，并对海州常山挥发性萜类进行了系统分析，旨在揭示该物种基因组进化与特化代谢之间的联系。研究结果首次提供了大青属的染色体水平基因组（约464.40 Mb，锚定到26条假染色体，BUSCO完整性98.39%，QV 33.70）。系统发育和K_s分析表明，海州常山经历了核心真双子叶植物共同六倍化（ECH/γ）、与其他唇形科共享的β-WGD事件（约56.7 Mya），以及谱系特异性的Ctri-WGD事件（约24.8 Mya）。基于亚基因组间基因保留的平衡性、无显著表达优势及较低的P-index（0.271和0.258），支持Ctri-WGD为同源四倍体化事件。通过HS-SPME-GC-MS共鉴定688种挥发性化合物，其中146种为萜类；结合差异表达分析，发现6个萜类合成酶（TPS）基因与萜类积累显著正相关，大多数与WGD相关拷贝有关。这些结果表明，WGD保留及随后的局部/分散重复，特别是TPS-a亚家族的谱系特异性扩增，促进了挥发性萜类的多样化。该研究为唇形科进化与功能研究提供了重要基因组资源。

论文解读：海州常山染色体水平基因组揭示复杂全基因组复制历史与挥发性萜类多样性

一、研究背景与问题
Clerodendrum trichotomum Thunb.（海州常山）是唇形科大青属的灌木或小乔木，广泛分布于中国、韩国、日本及菲律宾北部。该种具有观赏价值，其破碎叶片和茎释放的特殊不愉快气味也是大青属许多物种的典型特征。在传统医药中，海州常山长期用于治疗高血压和炎症，嫩叶在中国西南地区被用作野生蔬菜。植物化学研究已表明该物种富含萜类、黄酮类、苯丙素类、苯乙醇苷类、酚苷类及生物碱等多种次生代谢产物，其中萜类含量最为丰富。挥发性萜类是唇形科植物中占主导地位的特化代谢物，在植物防御、传粉和环境适应中发挥关键化学调控作用。然而，尽管该科植物已知化学多样性丰富且挥发性物质具有生态重要性，海州常山的挥发性化合物却尚未得到系统分析。因此，研究人员组装了该物种的高质量染色体水平基因组，并对其挥发性萜类进行系统分析，以揭示基因组进化与特化代谢之间的联系。

二、研究总体内容与结论
研究人员首次完成了大青属染色体水平的基因组组装，并基于比较基因组学解析了海州常山复杂的全基因组复制（whole genome duplication, WGD）历史。结果表明，海州常山在经历核心真双子叶植物共同六倍化（ECH/γ）和唇形科共享的β-WGD之后，又发生了一次谱系特异性的近期同源四倍体化事件（Ctri-WGD）。亚基因组进化的平衡特征支持该事件为同源四倍体化。进一步结合代谢组和转录组分析，研究人员揭示了WGD保留及后续局部/分散重复对挥发性萜类多样化的贡献。该研究为唇形科进化与功能基因组学研究提供了重要资源。

三、关键技术方法
研究人员利用Oxford Nanopore Technologies（ONT）长读长（约66×）、Illumina短读长和Hi-C染色质构象捕获（约322×）数据进行基因组组装，并采用purge_haplotigs去除单倍型冗余，最终获得染色体水平基因组。组装质量通过BUSCO完整性和QV值评估。转座子注释采用基于结构的方法，基因预测结合从头预测和同源证据。进化分析使用多物种基因家族聚类、共线性分析、同义替换率（K_s）分布、系统发育定年及P-index计算。挥发性物质通过顶空固相微萃取结合气相色谱-质谱联用（HS-SPME-GC-MS）分析，TPS基因表达通过RNA-seq差异表达和Pearson相关分析评估。

四、主要研究结果
（1）基因组组装与注释特征
基于26.12 Gb Illumina数据的k-mer分析估计基因组大小约446.96 Mb，流式细胞术测定为468.38 Mb。最终组装基因组为464.40 Mb，contig N50为8.03 Mb，其中461.46 Mb（99.37%）被锚定到26条假染色体。BUSCO评估显示完整性为98.39%，QV值为33.70，表明组装高质量。转座子（TE）占基因组的45.83%，其中长末端重复逆转录转座子（LTR）占15.65%，以Gypsy（7.17%）和Copia（6.58%）超家族为主；LINE占15.29%，主要来自L1类型；Helitron DNA转座子占7.02%。LTR插入时间分析显示近期爆发，主峰约0.46百万年前（Mya），中位插入时间0.88 Mya。基因注释共预测34,778个蛋白编码基因，96.4%获得功能注释。

（2）系统发育与基因家族进化
研究人员将海州常山与9个代表物种（包括唇形科近缘种、其他唇形目物种及核心真双子叶外群）进行基因家族聚类，共获得24,826个直系同源群。海州常山中29,907个基因（85.99%）聚类到14,604个基因家族，其中3,341个基因属于718个物种特异家族。进化分析鉴定到2,660个扩张和1,498个收缩基因家族（p<0.05）。扩张家族显著富集于硫转运、半胱氨酸和蛋氨酸代谢通路，推测与其先前报道的二氧化硫耐受和高吸收能力相关；收缩家族显著富集于跨膜运输通路。

（3）全基因组复制事件鉴定
通过精细的重复基因分类框架，研究人员将19,854个基因归类为典型WGD/片段重复，6,994个基因归类为源于古老WGD区域的次级重复（包括分散重复、近端重复和串联重复），4,263个基因为谱系分化后的新建重复。基因组自比对鉴定出24,335对旁系同源基因对应1816个共线性区块。K_s分布显示海州常山存在三个峰值（0.438、0.948、2.155），分别对应近期谱系特异WGD（Ctri-WGD）、共享β-WGD和核心真双子叶植物六倍化（ECH/γ）。与葡萄（Vitis vinifera）和黄芩（Scutellaria baicalensis）的共线性比较支持海州常山经历了额外的基因组复制。系统发育定年估计β-WGD约为56.7 Mya，Ctri-WGD约为24.8 Mya。

（4）亚基因组进化分析
以S. baicalensis和S. barbata为参考，海州常山的26条染色体被划分为亚基因组A和亚基因组B。两种参考下亚基因组A/B保留的基因数分别为9,191/8,217和9,215/8,592，基因保留率差异仅为3.5%和1.7%。在叶、花、果实组织中没有观察到明显的亚基因组表达优势。P-index值分别为0.271和0.258，均低于0.3。这些平衡的基因保留、无显著表达优势以及低P-index共同表明Ctri-WGD最可能源于同源四倍体化（autotetraploidization）。

（5）挥发性萜类与TPS基因关联分析
利用HS-SPME-GC-MS，研究人员在海州常山的叶、花和果实中共鉴定688种挥发性化合物，其中146种为萜类，是主要挥发性组分。单萜在叶和花中占优势，倍半萜在果实中占优势。在25个萜类合成酶（terpene synthase, TPS）基因中，15个在组织间差异表达。Pearson相关分析（经FDR校正）发现6个TPS基因与萜类含量显著正相关（r=0.890–0.981），其中TPS-a亚家族的CtriTPS18和CtriTPS1关联最广泛。这些正相关基因大多可归因于WGD/片段重复或源于古老WGD区域的次级重复。基因得失分析显示TPS-a亚家族在海州常山中表现出最强的谱系特异性拷贝数增加（+7/?3）。结果提示WGD保留是TPS基因进化的主要基因组库，后续局部或分散重复进一步促进了拷贝特异性表达与萜类积累的关联。

五、总结讨论与结论
研究人员认为，海州常山经历了古老的核心真双子叶植物六倍化、唇形科共享的β-WGD以及约24.8 Mya发生的近期谱系特异WGD事件。亚基因组间的基因保留平衡、无明显表达优势以及低P-index共同支持该近期WGD为同源四倍体化事件。WGD保留以及后续的局部或分散重复，尤其是TPS-a亚家族的扩张，很可能促进了挥发性萜类的多样化。研究结论翻译如下：本研究组装了海州常山的染色体水平基因组，首次为大青属的同源四倍体化和进化提供了基因组学见解。进化分析表明海州常山经历了ECH/γ、与其他唇形科共享的β-WGD，以及约24.8 Mya的谱系特异性Ctri-WGD事件。亚基因组进化、无偏基因保留和低P-index的综合证据支持Ctri-WGD最可能源于同源四倍体化。WGD保留和次级局部/分散重复，特别是TPS-a扩张，可能促进了挥发性萜类多样化。总之，这些发现阐明了海州常山的多倍化历史、适应性进化和挥发性萜类多样性，为唇形科进化与功能研究提供了关键基因组资源。