Blumea axillaris（菊科：Blumea属）的完整叶绿体基因组及其系统发育分析

首次完成B. axillaris叶绿体基因组测序，揭示其151,043 bp基因组具有典型四分区结构，含134个基因（90蛋白编码、36 tRNA、8 rRNA），并发现LSC区含大倒位（13 kb）嵌套小倒位（4.9 kb）。系统发育分析证实其与B. oxyodonta为姐妹群，支持属内B. lacera分支的系统地位。

　　
Blumea axillaris的叶绿体基因组研究揭示了该物种在玄参科中的系统发育地位及基因组特征。作为药用植物，Blumea axillaris在传统医学中应用广泛，其全基因组测序填补了该属基因组数据的空白。研究团队通过采集广西合浦县的野生植株叶片，采用Illumina测序平台完成基因组组装，最终获得151,043bp的完整叶绿体基因组序列，该长度与同科其他物种保持一致，属于典型四区体结构（LSC-IRa-SSC-IRb-LSC）。基因组包含134个功能基因，其中包含90个蛋白质编码基因、36个tRNA和8个rRNA，这一基因组成在现存研究数据中属于中等偏高水平。

在基因组结构分析中发现，该物种存在两个倒位结构。其中一个位于LSC区的大规模倒位（约13kb）内，嵌套着另一个4.9kb的小型倒位。这种双倒位结构在玄参科其他属种中并不常见，但与近缘物种Barnadesia caryophylla的基因组存在关联性特征。倒位区域可能影响基因表达调控或维持基因组稳定性，具体机制仍需进一步研究。

系统发育分析方面，研究团队选取了13个近缘物种作为对比，包括Aster属和Carpesium属等代表类群。通过ML和BI两种方法构建的系统发育树均显示，Blumea axillaris与B. oxyodonta形成姐妹群关系，且该结论在两次独立分析中得到高度支持（ML bootstrap≥95%，BI posterior probability≥0.95）。值得注意的是，B. axillaris与B. lacera亚族的进化关系存在分化，其系统发育位置更接近B. densiflora亚族，这与基于ITS和trnL-trnF片段的研究结果存在差异。这种矛盾可能源于传统分子标记与全基因组数据的进化信息差异，提示需要更全面的数据支持属内分类。

基因组成分析显示，该物种存在独特的基因拷贝现象。rps12基因存在两份拷贝，每份包含三个外显子和一个内含子。这种多拷贝现象在叶绿体基因组中较为罕见，可能与特定代谢途径的进化压力相关。此外，15个基因包含内含子结构，其中atpF、ndhA等关键基因的内含子分布符合植物叶绿体基因组的典型特征。

研究特别关注了基因组结构变异对系统发育的影响。通过对比B. axillaris与B. oxyodonta的叶绿体基因组，发现两者在倒位结构、重复序列分布等方面存在显著差异。其中B. axillaris特有的LSC区倒位结构，可能成为区分该物种与其他近缘种的重要分子标记。这种结构变异不仅影响基因组的物理结构，还可能改变调控元件的空间排列，从而影响基因表达模式。

在应用价值方面，研究团队指出完整叶绿体基因组的建立为该属物种的DNA条形码开发奠定了基础。通过比对不同物种的基因组特征，可以建立更精准的物种鉴定体系。例如，B. axillaris与B. oxyodonta的亲缘关系在基因水平上得到确认，这为中药制剂中区分不同Blumea属药材提供了分子证据。此外，基因组中发现的重复序列和倒位结构，可能为基因编辑技术提供靶点区域。

未来研究方向建议结合多组学数据深入探讨基因组变异的机制。例如，通过比较转录组数据与基因组结构，可以揭示倒位区域对特定基因表达的影响。同时，扩大物种采样范围至整个Blumea属，结合核基因组数据（如ITS区域）的分析，将有助于完善该属的系统发育框架。在资源开发方面，研究建议利用基因组的重复序列设计特异性分子标记，用于田间种植材料的快速鉴定。

该研究在方法学上实现了重要突破，首次在Blumea属中完成全基因组测序。通过优化测序深度（平均4118.86倍覆盖）和组装算法（GetOrganelle v1.7.7），成功克服了植物叶绿体基因组中常见的高GC区域和重复序列带来的组装难题。基因组注释工具Plastid Genome Annotator的应用，显著提高了基因功能的鉴定准确度，特别是对tRNA基因的识别率达到100%。

在生态适应方面，研究发现该物种的GC含量（37.49%）显著低于同科其他物种平均水平，这可能与高原生境下的紫外线辐射适应机制相关。通过比较不同海拔样本的基因组变异，未来研究可深入探讨植物在山地生态系统中的适应性进化策略。

综上所述，该研究不仅为Blumea axillaris提供了完整的分子遗传学档案，更在玄参科系统发育学领域取得重要进展。其揭示的基因组结构特征和系统发育关系，为药用植物的资源保护与开发提供了关键分子依据，同时为理解叶绿体基因组在被子植物进化中的功能提供了新案例。后续研究可结合地理种群采样和表型数据分析，建立更全面的物种鉴定体系，促进该药用植物的产业化应用。