尼佩塔属植物（Nepeta glomerulosa）的遗传多样性、种群结构及环境影响：对伊朗-图兰地区生物保护工作的启示

《Ecology and Evolution》：Genetic Diversity, Population Structure, and Environmental Influences in Nepeta glomerulosa: Implications for Conservation in the Irano-Turanian Region

　　**摘要**

本研究首次全面评估了Nepeta glomerulosa的遗传多样性、种群结构和进化历史。这种半特有药用植物分布于伊朗-图兰地区，具有多种治疗用途。我们使用序列相关扩增多态性（SRAP）标记，分析了来自伊朗主要山脉（Alborz、Zagros和Khorassan-Kopet Dagh）的18个种群的188个样本。共扩增出251条带，其中99.20%（249条）为多态性条带。该物种表现出中等程度的遗传多样性（Nei基因多样性指数He = 0.2605），其中Alborz山脉的多样性最高，这表明它可能是该物种的起源中心。种群间的遗传分化显著（Gst = 0.4199），分子方差分析（AMOVA）显示30.17%的总变异可归因于种群间的差异，而69.83%的变异来自种群内部。通过UPGMA聚类分析和主坐标分析（PCoA），样本被分为两个主要群体；而贝叶斯STRUCTURE分析则识别出三个具有显著混合的遗传群体，尤其是在Zagros山脉的种群中。这一模式表明Zagros山脉可能是北部和东南部种群之间重要的遗传交流区域。空间和环境因素，包括地理距离、海拔和土壤条件，对观察到的遗传结构有显著影响。野外观察进一步表明，N. glomerulosa种群受到过度放牧、栖息地退化和不可持续采伐的威胁。鉴于其分布范围有限且种群规模较小，迫切需要采取有针对性的保护措施。

**1 引言**

特有性是指一个物种仅限于一个明确的地理区域，这通常是由历史、生态和进化过程共同作用的结果（Anderson 1994；Noroozi等人2019）。伊朗位于西南亚，拥有丰富的地貌多样性，包括广阔的干旱至半干旱山区、山间盆地和高原。这种明显的地形和气候异质性在促进物种多样化和特有类群的出现中发挥了核心作用。伊朗拥有极其丰富的生物多样性，包括超过8000种维管植物（Khalvati等人2025；Noroozi等人2016, 2018）。此外，该国还栖息着1000多种大陆脊椎动物（Firouz 2005），进一步凸显了其在伊朗-图兰地区作为生物多样性和特有性重要中心的重要性。Nepeta glomerulosa属于唇形科（Lamiaceae）中较大的属之一，包含超过200个物种。该属的多样性中心分布在西喜马拉雅、欧亚大陆和西南亚（Budantsev 1993；Harley等人2004；Sharma等人2021）。由于其多种治疗特性，如利尿、祛痰、抗痉挛和抗菌作用，Nepeta属植物引起了制药行业的兴趣（Amirmohammadi等人2019；Talebi等人2022）。Nepeta glomerulosa是一种半特有物种，传统上用于治疗各种皮肤和胃肠道疾病、肺炎及瘙痒症状（Amirmohammadi等人2019；Talebi等人2018）。该植物为多年生草本植物，基部木质化，表面覆盖着白色绒毛状鳞片，并具有无柄的腺毛。叶片呈卵形，叶脉网状明显，边缘有锯齿。花为紫色，排列成穗状花序；每朵花都由白色至紫色的苞片衬托。花萼与苞片长度相等，有5条脉，下部锯齿通常比上部长。果实为长方形坚果状（Jamzad 2012；Rechinger 1982）。Nepeta glomerulosa分布于伊朗的Alborz、Zagros和Khorassan-Kopet Dagh山脉，以及阿富汗的Paropamisus山脉（Bakhshi等人2021；Rechinger 1982；图1）。它生长在岩石和砾石坡上，常与Astragalus L.、Artemisia L.、Amygdalus L.和Pistacia L.等植物群落共存。此外，它也出现在伊朗-图兰生物地理区域的干燥河流和泉水中（Karami等人2022）。该物种表现出显著的形态多样性，包括营养性状（如植株高度）和生殖性状（如苞片颜色和大小、花萼尺寸、花序形态）的变异。此外，花序形态在植物生命周期中会发生显著变化，尤其是在开花开始和结束的过渡阶段（Bakhshi等人2021；图2）。Nepeta glomerulosa包含四个亚种，其中三个亚种的分布范围有重叠，因此难以区分。这些亚种均为伊朗特有：（1）N. stapfiana，分布于伊朗西部和中部；（2）N. glomerulosa，分布于中部地区；（3）N. carmanica，分布于伊朗中部、西部、东南部和东部（Bakhshi等人2021；Rechinger 1982）；第四个亚种N. ghorana发现于阿富汗西南部的Ghorat（Govaerts 2009）。

**2 材料与方法**

**2.1 样本信息**

本研究从Nepeta glomerulosa的已知分布范围内收集了188个个体用于形态分析，以及134个个体用于分子分析。所选种群涵盖了该物种的全部地理、生态和环境变异，包括海拔、生境类型和气候条件的差异。在每个种群中，以足够的空间间隔采样，以减少采集到亲缘关系密切的植物的可能性，并仅选择健康成熟的个体。野外采集工作在主要生长季节（5月至6月）进行，此时诊断性形态特征已经完全发育。用于分子分析的叶片材料在野外采集后立即保存，同时制备凭证标本并存放在马什哈德费尔多西大学（FUMH）的标本馆中。遗憾的是，我们未能在阿富汗西南部的Parapamis山脉进行野外工作，因为那里有N. glomerulosa subsp. ghorana亚种。采样的分类单元、凭证信息和种群ID见表1。

**2.2 分子研究**

共使用了来自18个种群的188株植物（每个种群5-23株）来分析遗传多样性和种群结构，详见表1。从每个种群采集新鲜叶片并在硅胶中干燥。

**2.2.1 基因组DNA提取和SRAP扩增**

我们采用了改良的CTAB协议从唇形科植物中提取DNA（Talebi等人2022；Zagorcheva等人2020），然后使用NanoDrop 2000c分光光度计和1.5%凝胶来评估提取DNA的数量和质量。评估了36对不同的SRAP引物对，其中12对引物根据PCR反应效率和条带数量被选中（表2）。SRAP-PCR扩增程序借鉴了之前用于其他唇形科物种的方法，并进行了修改（Aghaei等人2017；Hmissi等人2024）。扩增循环参数如下：初始变性94°C 3分钟，随后是5个循环的94°C变性60秒、35°C退火60秒、72°C延伸60秒；接下来的35个循环中，初始变性94°C 60秒、50°C退火60秒、72°C延伸60秒，最后是一个72°C的延伸步骤10分钟。扩增产物在4°C下保存。扩增后，通过1.5%琼脂糖凝胶进行电泳分离。

**2.2.2 数据统计与分析**

扩增片段根据同源条带的存在（1）或缺失（0）进行评分，只有可重复的条带被纳入分析。然后将条带数据转换为每个引物组合的二元矩阵（Robarts和Wolfe 2014）。我们使用POPGENE版本1.32（Yeh等人1999）比较不同引物组合的扩增结果，包括位点数量、片段大小和多态位点百分比。多态信息含量（PIC）通过公式PIC = 1 ? p2 ? q2计算（p表示条带比例，q表示无条带比例）（Sharma等人2016）。计算了遗传指数，包括样本大小、等位基因数量（Na）、有效等位基因数量（Ne）（Kimura和Crow 1964）、Nei基因多样性（He）（Nei 1973）和Shannon信息指数（I）（Lewontin 1972），以及不同Nepeta glomerulosa种群内的遗传多样性。还确定了总遗传多样性（Ht）、种群内多样性（Hs）和遗传分化系数（Gst）。基因流（Nm）使用公式Nm = 0.5(1—Gst)/Gst计算。18个种群的主坐标分析（PCoA）使用GenAlex版本6.5（Peakall和Smouse 2006）进行。分子方差分析（AMOVA）使用Arlequin版本3.5.2.2（Excoffier和Lischer 2010）进行，以评估种群内部和种群间的遗传变异。Phylip版本3.6981软件包（Felsenstein 1993）被用来计算遗传距离矩阵。基于遗传距离矩阵，使用Phylip软件包构建了未加权对组平均法（UPGMA）（Sneath和Sokal 1963）树。为了调查N. glomerulosa样本系的种群遗传结构，使用了STRUCTURE 2.3的common-line版本进行贝叶斯种群分配分析（Pritchard等人2000）。分别使用K=1–12的值来运行STRUCTURE。每次运行的初始燃烧期为100,000次，伴随着100,000次蒙特卡洛马尔可夫链迭代。使用R语言实现的pophelper包（Francis 2017）来找到最优的K值。Ubuntu 22.0.4中的LibreOffice用于定制样本系的图形表示。使用ΔK方法（Evanno等人2005）计算了亚种群的数量。根据成员概率（Q值）大于70%或小于70%，N. glomerulosa样本系被分别分类为亚种群或混合亚种群（Liu等人2013）。ClustVis版本2.0是一个用于多变量数据聚类可视化的Web工具（BETA [https://biit.cs.ut.ee/clustvis]），用于生成热图，展示不同种群间SRAP位点的相关性和差异，并进行主成分分析（PCA）。

2.3 形态学研究

在形态学研究中，共采集了18个种群中的134株植物（每种群3-8株），并使用《伊朗植物志》（Jamzad 2012）和《伊朗植物志》（Rechinger 1982）进行鉴定。在第一阶段，评估了40个形态特征（数据S1）。我们使用IBM SPSS Statistics软件（版本26）来获取相关性。对于定量数据，我们使用了皮尔逊相关性；对于定性数据，我们使用了斯皮尔曼相关性。然后，对于不相关且不呈正态分布的特征，我们使用了方差分析（ANOVA）；对于不呈正态分布或定性的特征，我们使用了克鲁斯卡尔-瓦利斯检验（Kruskal-Wallis test）；具有显著差异的特征被选为下一阶段的分析对象。最终，选择了五个主要特征（植株高度（ph）、苞片长度（bl）、萼片长度（cl）和花冠长度（ctl）以及苞片颜色（bc）进行最终分析（数据S2）。测量使用了尺子、数字卡尺、清洁器和光学显微镜。

2.3.1 形态学分析

在R版本4.4.1中，对来自18个种群的134个个体的形态测量数据进行了主成分分析（PCA），分析了四个定量特征（植株高度（ph）、苞片长度（bl）、萼片长度（cl）和花冠长度（ctl）以及作为定性特征的苞片颜色（bc）。分析使用了PAST版本4.17（Hammer等人2001）来进行冗余分析（RDA），以模拟环境预测因子、基因组变异和形态测量变异之间的线性关系。

2.4 环境数据

在采样过程中，记录了伊朗不同地区每个种群的地理和土壤特征。地理变量包括经度、纬度、海拔、坡度和朝向，使用手持GPS设备（精度±3-5米）、测斜仪和指南针进行现场测量。选择这些变量是因为它们能够捕捉空间异质性和潜在的基因流障碍。从根区（深度5-30厘米）收集土壤样本，以评估土壤特性，包括氮（N）、磷（P）、钾（K）、碳（C）、钙（Ca2+）、pH值和电导率（EC），这些因素直接影响植物营养、生长和局部适应。土壤分析遵循标准实验室协议：总氮使用凯氏定氮法测量，钾使用火焰光度法测量，土壤有机碳使用Walkley-Black方法测量，电导率（EC）使用EC计测量，磷使用Olsen方法测量，pH值使用Jenway pH计测量。环境数据的冗余分析（RDA）和Mantel检验分别使用Past软件（版本3.14）和R（版本4.4.1）进行。

3 结果

3.1 分子结果

3.1.1 SRAP扩增产物的多态性分析

共筛选了36对SRAP引物组合。12对的扩增条带清晰且分布均匀，用于估计N. glomerulosa的遗传多样性。这些SRAP引物组合扩增了251个位点（249个多态性位点，占99.20%），大小范围为50-1700 bp。引物组合F9/R10（94.44%）和F9R20（95.45%）显示出最低的多态性，而引物组合F8/R8、F8R9、F8R10、F8R11、F9R8、F9/R9、F9R11、F9R12、F913和F9R21的多态性最高（100%）。多态性信息含量（PIC）的范围是从0.20到0.33（表2）。

3.1.2 遗传多样性和遗传结构分析

使用POPGENE程序得到的遗传多样性参数结果显示在表3中。物种水平的Na、Ne、H和I值分别为1.992、1.4286、0.2605和0.4048。在种群水平上，假设哈迪-温伯格平衡（Hardy–Weinberg equilibrium），种群DF（Na=1.7211, Ne=1.3608, H=0.2177, I=0.335）具有最高的遗传多样性，其次是种群SA（Na=1.1594, Ne=1.1275, H=0.0708, I=0.1014），其遗传多样性最低。平均Ht、Hs和Gst值分别为0.2663、0.1545和0.4199。Nm值为0.6906，小于1，表明种群间的基因交流受到限制。AMOVA的结果显示，30.1698%的变异存在于种群间，69.8301%的变异存在于种群内（表4）。

3.1.3 种群聚类分析

对18个N. glomerulosa种群的134个个体进行了主坐标分析（PCoA），分析了四个定量特征（植株高度（ph）、苞片长度（bl）、萼片长度（cl）和花冠长度（ctl）以及作为定性特征的苞片颜色（bc）。分析使用了PAST版本4.17（Hammer等人2001）进行冗余分析（RDA），以建模环境预测因子、基因组变异和形态测量变异之间的线性关系。

3.2 环境数据

在采样过程中，记录了伊朗不同地区每个种群的地理和土壤特征。地理变量包括经度、纬度、海拔、坡度和朝向，使用手持GPS设备（精度±3-5米）、测斜仪和指南针进行现场测量。选择这些变量是因为它们能够捕捉空间异质性和潜在的基因流障碍。从根区（深度5-30厘米）收集土壤样本，以评估土壤特性，包括氮（N）、磷（P）、钾（K）、碳（C）、钙（Ca2+）、pH值和电导率（EC），这些因素直接影响植物营养、生长和局部适应。土壤分析遵循标准实验室协议：总氮使用凯氏定氮法测量，钾使用火焰光度法测量，土壤有机碳使用Walkley-Black方法测量，电导率（EC）使用EC计测量，磷使用Olsen方法测量，pH值使用Jenway pH计测量。环境数据的冗余分析（RDA）和Mantel检验分别使用Past软件（版本3.14）和R（版本4.4.1）进行。

3.1.3 种群聚类分析

18个N. glomerulosa种群的主坐标分析（PCoA）显示，这些种群被分为两个不同的亚组（图3）。亚组I包括种群MB、ZGD、SA和QA，亚组II包括种群RLH、BJ、SHQ、SBR、SHS、SMO、YSD、SHF和BDR。种群CHD和DF分布在两个亚组之间。

3.1.4 热图和PCA

图S2显示了一个热图，展示了不同种群间SRAP位点的相关性和差异。顶部的分支图聚集了相似的位点，左侧的分支图聚集了相似的种群。沿垂直轴可以识别出3个主要簇，水平轴也显示了3个不同的簇。蓝色区域代表低值（低相关性），红色区域表示高值（高相关性），白色空间对应于缺失或中性的值。图S3展示了主成分分析（PCA）的结果。x轴代表主成分1（PC1），解释了总变异的21.1%，而y轴代表主成分2（PC2），解释了9.7%的变异。这两个成分共同提供了关于种群间遗传变异主要轴的见解。

3.1.5 冗余分析（RDA）

基于环境和遗传数据的冗余分析（RDA）结果显示在图6中。土壤因素，包括磷（P）、碳（C）、钙（Ca2+）、电导率（EC）和pH值，促进了种群SHQ、SHS、YSD和SHF的隔离。氮（N）因素与种群DF的区分有关，坡度影响了种群CHD的隔离。朝向在分离种群QS、SBR、SMO、BDR和CHM方面起到了作用。海拔与种群SA、QA、BLH和BJ的隔离有关。地理因素如纬度和经度进一步区分了种群TMK和MB/ZGD。

3.1.6 Mantel检验

Mantel检验的结果（数据S2）基于斯皮尔曼等级相关性，显示了遗传变量和地形变量之间的显著相关性（R=0.07279, p<0.05），以及遗传变量和地理变量之间的相关性（R=0.1872, p<0.01）。

3.2 形态学结果

3.2.1 主成分分析

对18个N. glomerulosa种群的主成分分析（PCA）显示，这些种群被分为三个不同的亚组（图7），其中两个亚组（亚组2；N. glomerulosa和亚组3；N. stapfiana）完全分离。第一个轴解释了52.4%的变异，第二个轴解释了19%的变异。在第一个轴中，属性cl、bc、ctl和bl起着最重要的作用；而在第二个轴中，属性bc和ph起着最重要的作用（图7和S4）。图7：在PowerPoint图查看器中打开

基于形态测量数据的18个N. glomerulosa种群的主成分分析（PCA）（1 = N. carmanica亚种；2 = N. glomerulosa亚种；3 = N. stapfiana亚种）。

3.2.2 冗余分析（RDA）
图8显示了基于形态测量和环境数据的冗余分析（RDA）结果。土壤因素（K、C、Ca2+、N、EC和PH）影响了SHS、SMO、BDR和CHM种群，使它们与其他种群分离。P因素导致DF种群与其他种群隔离。坡向因素导致SHQ、SBR和YSD种群与其他种群隔离。海拔和坡度因素导致TMK、RLH和SHF种群隔离，纬度因素导致MB和CHD种群隔离，而纬度与经度因素的组合导致ZGD、SA、QA和BJ种群与其他种群隔离。图8：在PowerPoint图查看器中打开

基于形态测量和环境数据的18个N. glomerulosa种群的冗余分析（RDA）；a = 基于种群的RDA；b = 基于个体的RDA。

3.2.3 曼特尔检验
根据Spearman等级相关性，曼特尔检验结果显示形态测量数据与地理数据之间存在显著相关性（R = 0.0821，p < 0.01）。然而，形态测量数据与地形数据之间未发现相关性（R2 = 0.0001621，p > 0.05）。

4 讨论
本研究首次对Nepeta glomerulosa的遗传和表型变异进行了分子和形态学特征描述，该物种是分布于伊朗Alborz、Zagros和Khorassan-Kopet Dagh山脉的特有唇形科植物。

4.1 N. glomerulosa的遗传多样性
新栖息地的殖民机会和有限的基因流使得植物种群变得孤立，从而降低了遗传多样性（Zhang等人，2019年）。迄今为止，N. glomerulosa的种群结构和遗传多样性尚未被报道。共扩增了251个位点，其中249个（99.20%）是多态的。Na、Ne、He和I等指数在样本集合间表现出变异性，并揭示了样本种群内部和之间的遗传多样性。从物种角度来看，N. glomerulosa的遗传多样性较低，He值为0.2605，低于先前遗传研究中报道的短命多年生植物的平均He值（0.55，Nybom，2004年）。从种群水平来看，Alborz山脉中的N. glomerulosa遗传多样性显著高于Zagros和Khorassan–Kopet Dagh山脉。Alborz山脉中的DF种群具有最高的遗传多样性，为72.11%，而Khorassan–Kopet Dagh山脉中的SA种群遗传多样性最低，为15.94%。东北部和南部种群的遗传多样性降低可能是由于北部种群与地理和生态隔离导致的基因流受限（Fetter和Weakley，2019年；McGreevy Jr.等人，2024年）。相比之下，伊朗北部没有显著的地理障碍，允许种群集中分布并增加基因流，从而增强了遗传多样性并减少了这些种群之间的分化（Hosseinian Yousefkhani和Nabizadeh，2025年）。此外，Alborz山脉可能是该物种的起源中心，这可能有助于其高遗传多样性。多种因素，包括基因流、选择、突变和遗传漂变，决定了种群间的遗传分化（Gst）。基因流可以抵消遗传分化（Oh和Oh，2022年）。Gst表示种群间变异与物种内部总遗传变异的比率，是评估种群遗传分化的重要指标。Gst值大于0.25表示种群间有显著的分化（Wright，1984年）。Nybom（2004年）和Hamrick（1987年）的研究表明，Gst接近0.5的物种表现出一定程度的自花授粉（Vilperte等人，2021年；Wang、Dai等人，2023年；Wang等人，2014年）。这一观察结果支持N. glomerulosa存在自花授粉的现象，其Gst值为0.4199，这归因于相邻种群之间的较大距离。在没有传粉者的情况下，自花授粉可以增加种群增长并确保物种的繁殖。这种策略有助于在该生态位中获得竞争优势（Sun等人，2005年）。然而，没有传粉者的自花授粉产生的种子质量较差，会降低N. glomerulosa的种群增长率（Law等人，2010年），这可能解释了其种群规模较小的原因（Karami等人，2022年）。

4.2 N. glomerulosa的遗传结构
尽管N. glomerulosa分布在不同的山脉系统中，但其遗传结构显示出相对较弱的分化。缺乏强烈的区域遗传分离表明种群间有广泛的基因流，尤其是在大多数种群集中的Zagros山脉。这种模式与大部分种群中存在混合现象一致，表明历史或持续的连通性限制了遗传分化。种群间观察到的广泛混合可能反映了先前分化谱系之间的反复杂交，这在分布范围重叠且连通性间歇性的山地植物物种中很常见（Abbott等人，2013年；Goulet等人，2016年；Mallet，2005年；Shriner，2023年）。在N. glomerulosa中，这种混合可能是由于过去的气候波动周期性地扩大了适宜的栖息地，使得目前地理上分离的种群之间能够接触。Zagros种群在多个遗传簇中的主导地位支持了Zagros山脉作为连接Alborz和Yazd–Kerman地区的主要遗传交换走廊的作用。Zagros山脉复杂的地形和异质性栖息地可能促进了局部分化和反复的基因流，导致了观察到的镶嵌遗传结构，而不是明确的系统地理学分离。

4.3 环境因素的影响
环境异质性在塑造N. glomerulosa种群间的遗传分化中起着重要作用。海拔影响温度和降水模式，低海拔地区降水的影响更大，而高海拔地区温度变得更为关键（Karami等人，2022年）。这些变化可能驱动了1500–4000米海拔范围内种群间的适应性分化。经度和纬度也促进了遗传分化。位于较高纬度的MB、ZGD和TMK种群与其他种群隔离（图6）。这表明地理距离可能成为基因流的障碍，导致随着经度或纬度的增加而产生更大的遗传距离。海拔、坡度和坡向的变化可以创建微生境不连续性，限制花粉和种子的传播，即使在短地理距离内也会导致部分隔离（Degen等人，2006年；Rousset，1999年；Torroba-Balmori等人，2017年）。这种地形障碍在山区尤其有效，可能解释了尽管整体区域结构较弱但仍存在局部分化的现象（Han等人，2025年）。在其他山地物种中也报告了类似的模式，其中环境梯度与地理相互作用产生了精细的遗传结构，而没有完全隔离（McGreevy Jr.等人，2024年；Morente-López等人，2018年；Wang和Bradburd，2014年）。在N. glomerulosa中，这表明遗传结构不是由严格的地理分离驱动的，而是由环境过滤和间歇性基因流之间的相互作用决定的。

4.4 N. glomerulosa的保护
在野外采样过程中，我们反复观察到几个N. glomerulosa种群受到严重的人为压力。这些观察结果虽然没有作为分析数据集的一部分进行量化，但一致显示了过度放牧、与城市扩张相关的栖息地丧失以及无控制的药物采集的影响。该物种的脆弱性因其种群规模小和分布范围有限而加剧。特别是山区似乎受到不成比例的影响，因为这些地区的快速城市化和人类活动的增加与栖息地退化和采集压力的加剧相吻合。基于这些野外观察，保护工作应优先考虑保护N. glomerulosa种群最易受干扰的山区栖息地。推荐的策略包括规范采集行为、建立保护区以及开发离体保护措施，如种子库。长期保护需要当地社区、政策制定者和研究人员的协调行动，以保护这一具有生态和药用价值的物种。

5 结论
本研究对伊朗-图兰地区的半特有药用植物Nepeta glomerulosa的遗传多样性和种群结构进行了全面评估。我们使用SRAP标记对来自伊朗主要山脉（Alborz、Zagros和Khorassan-Kopet Dagh）的18个种群的188个样本进行了分析。该物种表现出中等程度的遗传多样性（He = 0.2605），其中Alborz山脉的多样性最高，表明它可能是起源中心。尽管如此，这种多样性低于短命多年生植物的平均水平，这引发了关于保护的担忧。识别出三个不同的遗传簇，对应于不同的地理区域，几乎一半的样本中观察到显著的混合现象，表明基因流和种群分化复杂。地理和环境因素，如海拔、土壤条件和地理距离，在塑造遗传结构中起着关键作用。鉴于该物种的分布范围有限、种群数量少以及过度放牧和栖息地退化带来的风险，立即采取保护措施至关重要，特别是在敏感的山区。对历史联系的分析表明其起源于Alborz地区，为了解其进化历史提供了宝贵的见解。这些结果为N. glomerulosa的成功保护方法奠定了基础，并强调了将遗传和环境因素纳入战略规划的重要性。未来需要更广泛的采样研究，特别是来自阿富汗的研究，以及下一代测序技术的应用，以深入理解这一具有药用价值的物种。及时的保护措施，包括现场保护、可持续采集实践和离体保护措施，对于确保其持续生存至关重要。地方社区、政策制定者和专家之间的合作对于成功实施保护措施至关重要。

作者贡献
Sahar Karami：概念化（同等贡献）、数据管理（同等贡献）、正式分析（同等贡献）、方法论（同等贡献）、撰写初稿（同等贡献）。
Hamid Ejtehadi：概念化（同等贡献）、数据管理（同等贡献）、资金获取（同等贡献）、项目管理（同等贡献）、撰写-审阅和编辑（同等贡献）。
Jamil Vaezi：数据管理（同等贡献）、撰写-审阅和编辑（同等贡献）。
Hamid Moazzeni：概念化（同等贡献）、数据管理（同等贡献）、资金获取（同等贡献）、项目管理（同等贡献）、监督（同等贡献）、撰写-审阅和编辑（同等贡献）。

致谢
我们衷心感谢Mohammad Reza Amiri（FUMH）、Shahram Bahadori（TUH）和Parisa Khoshniat在野外采样期间的宝贵帮助。我们还要感谢Mohammad Reza Joharchi（FUMH）在形态测量分析方面的指导，以及Ramiar Majidi慷慨分享的摄影材料。第一作者（S.K.）得到了马什哈德Ferdowsi大学（资助编号3/53411）的财务支持。

利益冲突
作者声明没有利益冲突。

数据可用性声明
作者确认支持本研究发现的数据可以在文章[和/或其支持信息]中找到。