《Metabolites》：DNA Barcoding and Untargeted Metabolomics Reveal Phylogenetic and Chemotaxonomic Relationships Among 11 Ferula Species from Xinjiang, China
编辑推荐：
背景/目的：阿魏属（Ferula）物种在形态上的相似性，尤其是在营养生长期，给叶片材料的可靠鉴定带来了困难。野生与栽培阿魏资源之间系统发育关系与叶片化学分类学模式之间的对应关系尚未得到充分表征。本研究针对中国新疆11种阿魏属物种探讨了三个问题：组合DNA条形码
背景/目的：阿魏属（Ferula）物种在形态上的相似性，尤其是在营养生长期，给叶片材料的可靠鉴定带来了困难。野生与栽培阿魏资源之间系统发育关系与叶片化学分类学模式之间的对应关系尚未得到充分表征。本研究针对中国新疆11种阿魏属物种探讨了三个问题：组合DNA条形码是否支持物种鉴定与宽泛的系统发育分组；叶片代谢谱是否揭示类群相关的分化；系统发育分组与代谢分组在多大程度上相互对应。方法：将ITS2、trnH–psbA、rps11和atpA–trnR四个DNA条形码位点组合用于分子鉴定和系统发育分析。采用响应面法优化超声辅助提取条件，以建立标准化的样品制备流程，随后对代表11个物种的33份叶片样品进行基于液相色谱-质谱联用（LC–MS）的非靶向代谢组学分析。采用多元统计分析表征类群相关的代谢模式及化学分类学关系。结果：组合DNA条形码数据集总体上支持对大多数所检样品的鉴定，并揭示了四个宽泛的系统发育组，尽管若干近缘种质的分辨率有限。非靶向代谢组学共推定注释了1972种代谢物，主要为萜类、黄酮类、酚酸类和香豆素类。主成分分析（PCA）揭示了类群相关的代谢模式，而层次聚类分析（HCA）将样品分为三个主要代谢簇。值得注意的是，新疆阿魏（F. sinkiangensis）、阜康阿魏（F. fukanensis）、羊食阿魏（F. lehmannii）、多伞阿魏（F. teterrima）和圆锥茎阿魏（F. conocaula）出现在同一宽泛系统发育组和同一代谢组内，表明分子关系与代谢表型之间存在部分对应。结论：DNA条形码与非靶向代谢组学为阿魏属叶片的鉴定和比较表征提供了互补的遗传与化学证据。所观察到的系统发育关系与代谢谱之间的部分对应，为理解所研究阿魏属物种间的化学分化提供了探索性基础。
论文解读
一、研究背景与科学问题
阿魏属（Ferula）隶属于伞形科（Apiaceae），全球约有180–185种，广泛分布于欧洲和亚洲，是重要的药用、芳香及食用植物资源。在中国约记录有30种，其中新疆是该属的重要多样性分布中心。传统上，阿魏属的利用集中于根部及根茎分泌的油胶树脂，这些树脂在中医、伊朗传统医学、阿育吠陀及阿拉伯医学中广泛用于治疗消化系统疾病、呼吸系统疾病、炎症性疾病、风湿病及皮肤病等。然而，既往研究主要聚焦于药用树脂，对可食用叶片资源的遗传鉴定和化学多样性研究相对不足。在新疆，多种阿魏属植物的鲜叶长期被作为蔬菜食用，如新疆阿魏（F. sinkiangensis）、阜康阿魏（F. fukanensis）、羊食阿魏（F. lehmannii）和多伞阿魏（F. teterrima）等，传统上被认为有益于消化健康。
然而，阿魏属物种在营养生长期叶片形态高度相似，仅凭形态特征难以可靠区分不同物种。新鲜叶片通常经过焯水冷冻保存，进一步削弱了形态特征的诊断价值。此外，不同物种的代谢物组成和积累模式可能存在差异，物种混杂会掩盖类群相关的化学特征，影响阿魏叶资源的化学评价一致性。因此，亟需能够准确鉴定形态相似的阿魏叶片并表征种间化学差异的可靠方法。DNA条形码通过靶向标准化遗传位点，所需生物材料少、重现性高、不依赖完整形态特征，特别适用于形态高度相似的植物类群鉴定。但传统DNA条形码对近缘或近期分化类群的区分能力可能受限。代谢组学通过全面表征代谢物组成和相对丰度，可提供互补的化学信息，用于评估种间化学变异和化学分类学关系。因此，本研究将多基因座DNA条形码与LC–MS非靶向代谢组学整合，对新疆11种阿魏属植物的33份叶片样品进行系统研究，旨在回答三个核心问题：组合DNA条形码是否支持物种鉴定和宽泛系统发育分组；叶片代谢谱是否揭示类群相关的分化；系统发育分组与代谢分组在多大程度上对应。
二、主要技术方法
研究人员于2023至2025年间从新疆多个地点采集了11种阿魏属植物的33份叶片样品，涵盖野生种群和栽培资源，所有材料经新疆医科大学徐海燕和张雪佳根据生殖期形态特征鉴定，凭证标本存放于新疆医科大学标本馆。研究采用四个DNA条形码位点（ITS2、trnH–psbA、rps11和atpA–trnR）进行分子鉴定和系统发育分析，分别构建邻接法（NJ）和最大似然法（ML）系统发育树。代谢组学方面，以阿魏酸为代表性酚类标记物，通过单因素实验和响应面法优化超声辅助提取条件，建立标准化样品制备流程，随后对全部叶片样品进行超高效液相色谱-四极杆-静电场轨道阱质谱联用（UHPLC–Q Exactive HF-X Orbitrap MS）非靶向代谢组学分析，采用Progenesis QI v3.0处理数据，MetaboAnalyst 6.0进行多元统计分析。
三、研究结果
形态特征观察：所研究的阿魏属物种呈现该属的一般形态特征，包括直立茎、高度分裂的复叶、复伞形花序、黄色花和背腹压扁的果实。物种间在株高、分枝方式、叶片分裂程度、花序结构和果实形态上存在较大变异，但多个近缘种的营养器官特别是叶片形态高度相似。可靠的形态鉴定依赖于生殖性状，而这些性状在营养期和加工后的叶片材料中不可获得，因此叶片形态单独用于物种鉴定的区分能力有限。
系统发育分析：基于四个条形码位点串联序列构建的NJ树将所研究的阿魏属种质与三个外类群（Angelica sinensis、Bupleurum commelynoideum、Chamaesium jiulongense）明确分开，并可将阿魏属种质初步分为四个宽泛系统发育组。第I组包含多伞阿魏、新疆阿魏、圆锥茎阿魏、阜康阿魏和羊食阿魏，被视为宽泛的蒜味阿魏类群；第II组包含F. akitschkensis和F. syreitschikowii；第III组包括F. ovina、F. dissecta、F. kirialovii、F. olivacea、未鉴定阿魏种质及F. feruloides 1号种质；第IV组主要包括F. songarica及F. feruloides 2号和3号种质。部分节点自展支持率较低，个别同种质未与同种其余样品聚在一起。ML分析恢复了NJ树中若干末端关系，但深层次分支模式存在差异。野生与栽培种质在系统发育树上未出现按栽培状态分离的模式。
提取条件优化与响应面分析：单因素实验表明乙醇浓度、液料比和提取时间均影响阿魏酸得率。响应面分析建立的二次多项式回归模型高度显著（F = 27.39，p = 0.0001），失拟项不显著，R2 = 0.9724，校正R2 = 0.9369。三因素对响应的影响顺序为提取时间 > 乙醇浓度 > 液料比。数值优化预测最佳条件为49.507%乙醇、液料比22.309 mL/g、提取时间34.735 min，预测阿魏酸得率7.962 mg/g。调整为50%乙醇、22 mL/g、35 min后，三次验证实验平均得率7.91 ± 0.06 mg/g，与预测值相对偏差仅0.52%，表明模型预测与实验验证高度一致。
代谢物注释概况：基于精确前体质量、MS/MS碎裂模式和数据库匹配，在11种阿魏叶片中共推定注释1972种代谢物。萜类为最大化学类别，共797种，占40.42%，其中倍半萜类225种（11.41%）；黄酮类485种（24.59%）；酚酸类及其衍生物253种（12.83%）；香豆素类及其衍生物177种（8.98%）。多种特征性倍半萜香豆素如farnesiferol A、colladonin、conferone和feropolidin被归入香豆素类。
主成分分析与偏最小二乘判别分析：PCA模型中PC1和PC2分别解释24.6%和19.7%的总方差，累计44.3%。同种样品总体趋向聚集，部分物种占据相对独立的得分图区域，另一些则部分重叠。基于PCA得分空间的PERMANOVA支持物种间显著整体分离（F = 126.75，R2 = 0.9829，p = 0.005）。PLS-DA的样品分布与PCA总体一致，交叉验证显示两成分模型R2 = 0.2658，Q2 = 0.0965，预测性能有限，但置换检验表明观察到的组间分离显著大于随机置换结果（p < 0.001）。VIP得分最高的代谢物包括dihydroxyacidissiminol、kiwiionoside、4-carvomenthenol、feropolidin、farnesiferol A、umbelliferone、colladonin、hesperetin、nerolidol、conferone和(E)-1-propene-1-sulfenic acid等，涵盖香豆素类、萜类、黄酮类和含硫化合物。
层次聚类分析：基于标准化LC–MS数据集的HCA将11种阿魏分为三个主要代谢组。A组包含F. songarica、F. akitschkensis、F. syreitschikowii和F. olivacea，其中F. akitschkensis与F. syreitschikowii代谢谱相对相似。B组由F. feruloides和F. ovina组成，两物种各自形成紧密亚簇。C组包含圆锥茎阿魏、新疆阿魏、羊食阿魏、阜康阿魏和多伞阿魏，这些物种形成独立的主要代谢组，具有食用记录的四种阿魏均在此组内。
野生与栽培种质比较：在NJ树、PCA和HCA中，同种野生与栽培种质均未出现一致的按栽培状态分离的模式。这可能与新疆阿魏栽培历史较短、栽培种质多源自野生资源有关，栽培环境对代谢物的影响尚不足以产生显著分化。
韦恩图分析：在C组五种阿魏中，1190种代谢物为五物种共有，多伞阿魏、圆锥茎阿魏、新疆阿魏、阜康阿魏和羊食阿魏分别有73、41、21、11和6种独有代谢物。B组中F. feruloides和F. ovina共有1596种代谢物，分别有58和204种独有。A组四物种共有1339种代谢物，F. songarica、F. akitschkensis、F. syreitschikowii和F. olivacea分别有51、9、40和31种独有代谢物。
热图分析：基于前50个VIP排名代谢特征的热图显示，圆锥茎阿魏、新疆阿魏、阜康阿魏、多伞阿魏和羊食阿魏的样品在大部分VIP代谢物上呈现大致相似的积累模式，多种香豆素和倍半萜香豆素相关代谢物如colladonin、farnesiferol A、feropolidin、umbelliferone和conferone在这些样品中相对丰度较高。相反，F. songarica、F. olivacea、F. syreitschikowii和F. akitschkensis在多种黄酮类代谢物如quercetin、isoquercitrin、myricetin、gossypin和quercetin-3-gentiobioside上相对丰度较高。F. feruloides和F. ovina也呈现特征性积累模式。
四、讨论与结论
讨论部分指出，仅凭叶片形态难以准确鉴定阿魏属物种，而组合DNA条形码总体上支持了形态鉴定并揭示了宽泛的系统发育结构。第I组中五种阿魏（新疆阿魏、阜康阿魏、羊食阿魏、多伞阿魏和圆锥茎阿魏）在系统发育和代谢分组中均聚在一起，表明较近的进化关系可能伴随叶片整体代谢表型的相似性，这从生物学上是合理的，因为次生代谢物生物合成依赖遗传编码的酶和调控通路。然而，这种对应并非普遍存在，环境条件、发育阶段和栽培历史也可能造成差异，因此支持系统发育关系与代谢分组之间存在部分对应。含硫代谢物与阿魏资源特征性葱蒜气味的关联值得关注，但现有数据尚不能建立单一代谢物与食用接受度、风味或生物活性之间的直接关系。研究局限包括每物种仅三个样品、地理来源与采集时间及栽培状态存在混杂、PLS-DA预测性能有限、代谢物为推定注释等。未来需扩大采样规模、系统评估环境和栽培因素、使用标准品确认代谢物身份。
结论部分翻译：本研究整合多基因座DNA条形码与LC–MS非靶向代谢组学，对中国新疆11种阿魏属植物的33份叶片样品进行了物种鉴定、系统发育关系和化学多样性评价。组合ITS2、trnH–psbA、rps11和atpA–trnR条形码数据集为大多数所检样品的鉴定提供了分子支持，并揭示了阿魏属类群间的宽泛系统发育关系，尽管若干近缘种质分辨率有限。非靶向代谢组学共推定注释1972种代谢物，其中萜类、黄酮类、酚酸类和香豆素类为主要化学类别。多元分析进一步揭示了类群相关的代谢模式和种间化学分类学分化。值得注意的是，代谢分组与DNA系统发育关系存在部分对应，尤其是系统发育第I组与代谢C组之间的对应，表明较近的系统发育关系可能伴随整体代谢谱的相似性。同时，两种分析框架之间的不完全对应表明系统发育亲缘关系与代谢表型提供相关但不完全相同的信息。所观察到的模式代表当前采样框架内的探索性类群相关代谢趋势。扩大采样、系统评估环境和栽培因素以及使用标准品确认代谢物注释为进一步验证提供了方向。总体而言，DNA条形码与非靶向代谢组学的整合为形态相似的阿魏叶片鉴定、比较表征和化学分类学评价提供了互补的分子与化学信息，为未来资源鉴定和质量评价奠定了基础。
