利用新型STI×YSI复合指数筛选耐盐伊朗香蜂花地方品种：一种识别有潜力候选者的方法学方法

《Scientia Horticulturae》：Screening salt-tolerant Iranian Moldavian balm landraces using a novel STI?×?YSI composite index: A methodological approach with identification of promising candidates

　　

背景：香蜂花（*Dracocephalum moldavica* L.）是一种有价值的药用植物，但其伊朗基因型对盐胁迫的响应仍缺乏系统研究。本研究引入了一种基于新型复合指数进行初步筛选耐盐基因型的方法学方法。方法：研究人员将10个伊朗基因型置于三个盐度水平（0、75和150 mM NaCl）下，采用完全随机设计，设置三个重复（n=90）。测量了农艺性状、抗氧化酶活性、丙二醛（MDA）和精油谱。研究人员开发了一种新型复合指数（STI×YSI），并通过Bootstrap分析（1000次重抽样）进行了验证。采用主成分分析（PCA）、聚类分析和多变量方差分析（MANOVA）进行基因型区分。结果：盐度显著降低了地上部干重（降低52%，从15.21 g降至7.29 g，p<0.01），同时增加了抗坏血酸过氧化物酶（APX）活性（从0.386增至0.517 μmol mg^-1 protein min^-1，p<0.001）。STI×YSI指数解决了高基线产量与毒物兴奋效应之间的悖论。PCA（解释83.1%方差，KMO=0.76）和聚类分析（共表相关系数=0.82）确定Azarshahr（0.70）、Maragheh（0.64）和Ajabshir（0.57）为得分最高的基因型，且排名稳定（>90% Bootstrap稳定性）。香叶醇的双相响应（75 mM NaCl时增加，150 mM NaCl时降低）支持了生长-防御权衡理论。SOD和脯氨酸的响应无法可靠评估，因为本研究的统计功效对这些性状不足（0.38–0.45）。结论：Azarshahr、Maragheh和Ajabshir是本次初步筛选中的有前景候选者。然而，这些发现需要在田间验证中补充离子测量（Na⁺/K⁺比值）并增加重复次数，才能用于育种推荐。该方法学框架（STI×YSI）为初步基因型筛选提供了一种稳健的方法。

**研究背景与问题提出**

土壤盐碱化是全球干旱和半干旱地区作物生产力的主要限制因素，全球超过8.33亿公顷土地受盐碱影响，且预计到2050年将大幅增加。伊朗约50%的耕地受盐碱化威胁，但本土药用植物基因型对盐胁迫响应的信息仍然匮乏，尤其对于具有药用和经济价值的物种。香蜂花（*Dracocephalum moldavica* L.）为唇形科一年生草本植物，其精油富含香叶醇、香叶醛和乙酸香叶酯等活性化合物，具有抗菌和镇静作用。然而，伊朗本土香蜂花基因型的耐盐性遗传多样性尚未系统研究，存在以下关键空白：（1）缺乏用于香蜂花基因型排序的标准耐盐指数计算；（2）受基线产量影响的STI（耐盐指数）与代表相对降低的YSI（产量稳定指数）之间存在概念悖论；（3）尚未开展结合PCA验证（KMO和Bartlett检验）及共表相关系数聚类分析的综合分析。此外，先前研究缺乏同时测量农艺、生化、精油性状与离子积累（Na⁺、K⁺）的数据，而离子稳态是耐盐机制的核心。本研究故意将离子测量排除在初步筛选阶段之外，以优先进行成本效益评估，所有关于耐受机制的结论均为暂定结论。毒物兴奋效应（低剂量刺激、高剂量抑制）给耐盐筛选带来挑战：表现出毒物兴奋效应的基因型（YSI>1）可能看似耐盐，但其低基线产量使其商业价值低。因此，本研究旨在：（1）同时评估盐度对10个伊朗基因型农艺性状、抗氧化酶和精油成分的影响；（2）提出新型复合指数（STI×YSI）以解决高基线产量与相对产量降低之间的悖论；（3）通过PCA和聚类分析初步筛选耐盐基因型。研究假设包括：香蜂花基因型间存在显著耐盐遗传变异；复合指数优于单一指数；耐盐基因型表现出独特的次生代谢物积累和抗氧化酶活性模式。该研究发表在《Scientia Horticulturae》。

**主要关键技术与方法**

研究人员采用完全随机设计，将10个伊朗香蜂花地方品种（采集自东阿塞拜疆和西阿塞拜疆省，凭证标本存放于Maragheh大学标本馆，馆藏号5121）置于3个NaCl浓度（0、75、150 mM）下，3次重复（n=90）。测量农艺性状（株高、叶数、地上部干重等）、生化指标（APX、GPX、SOD活性、MDA、脯氨酸含量）及精油成分（GC-MS分析）。计算传统耐盐指数（STI、SSI、GMP、MP、TOL、YSI），并开发新型复合指数STI×YSI，通过Bootstrap重抽样（1000次）验证稳定性。采用PCA（KMO和Bartlett检验）、聚类分析（Ward法，共表相关系数）及MANOVA进行基因型区分。

**研究结果**

**3.1 统计假设验证与事后功效分析**
通过Kolmogorov-Smirnov检验和Levene检验验证数据正态性和方差齐性。事后功效分析显示，地上部干重（0.92）、APX（0.87）和MDA（0.85）功效充足（≥0.8），但SOD（0.45）和脯氨酸（0.38）功效低，导致这些结果不具有结论性。

**3.2 农艺性状与基因型×盐度互作**
方差分析显示，基因型显著影响株高和花序长度（p<0.01），盐度显著影响大多数农艺性状（除花序数）。基因型×盐度互作仅对地上部干重显著（p<0.01，η2_p=0.86）。Malekan在非胁迫下地上部干重最高（22.14 g），而Hashtrud在150 mM NaCl下干重（15.08 g）显著高于75 mM（7.86 g），表现出毒物兴奋效应。Miandoab、Marand等基因型在严重盐胁迫下干重严重降低。

**3.3 耐盐指数与新型复合指数**
基于STI，Malekan（1.24）、Azarshahr（1.18）和Ajabshir（1.09）最高。但Hashtrud因低基线产量（11.98 g）而STI低（0.49），尽管其YSI=1.26。新型复合指数STI×YSI解决了该悖论：Azarshahr（0.70）、Maragheh（0.64）和Hashtrud（0.62）数值排名最高，但Hashtrud因低基线产量和不稳定排名（48.7% Bootstrap稳定性）被排除在育种推荐之外。Bootstrap分析显示Azarshahr、Maragheh和Ajabshir的排名在94.2%重抽样中稳定。

**3.4 抗氧化酶与丙二醛（MDA）**
盐度显著增加APX活性（从0.386增至0.517 μmol mg^-1 protein min^-1，p<0.01）和MDA含量（从0.619增至1.110 nmol g^-1 FW，p<0.01）。基因型间APX和MDA差异显著。Azarshahr和Maragheh具有最高APX活性和最低MDA，显示高效膜保护。GPX活性仅Ajabshir和Hashtrud显著增加。SOD和脯氨酸因统计功效低（0.45和0.38）结果不具结论性，可能为假阴性，非生物学无效。

**3.5 精油成分**
共鉴定31种化合物，主要成分为乙酸香叶酯（38.26%）、香叶醛（28.31%）、香叶醇（15.53%）和Z-柠檬醛（6.42%），合计占88.52%。香叶醇和香叶醛存在显著基因型×盐度互作（p<0.01和p<0.05）。香叶醇呈现双相响应：75 mM NaCl时在多个基因型中增加（如Miandoab从13.92%升至22.70%），150 mM NaCl时降低。精油产量随盐度显著降低，Malekan产量最高（0.310 ml/pot）。

**3.6 相关性矩阵**
STI与地上部干重（r=0.91）和香叶醇（r=0.68）强正相关，与SSI强负相关（r=-0.87）。MDA与GPX呈弱正相关（r=0.32，不显著）。

**3.7 主成分分析（PCA）**
KMO=0.76，Bartlett检验显著（χ2=268.4，p<0.001），前4个主成分解释83.1%方差。旋转后：PC1（耐性与品质维持，25.7%）与STI（0.912）、地上部干重（0.889）、香叶醇（0.862）正相关；PC2（生长与生物量，23.8%）；PC3（精油品质，17.3%）；PC4（抗氧化防御，16.3%）。双标图显示Cluster 2基因型（Azarshahr、Maragheh、Ajabshir、Malekan）位于PC1正侧，与耐性相关向量一致。

**3.8 多变量分析与优劣基因型鉴定**
聚类分析（共表相关系数=0.82）和PCA一致识别Azarshahr、Maragheh、Ajabshir和Malekan为优良基因型。MANOVA确认两簇显著分离（Wilks' Λ=0.183，F(18,42)=5.27，p<0.001，η2_p=0.69）。Cluster 2基因型地上部干重、香叶醇和GPX活性分别比Cluster 1高13.08%、15.50%和65.01%。

**讨论与结论**

讨论部分首先指出盐度显著降低农艺性状，与先前研究一致，但缺乏离子测量限制了离子毒性贡献的确定。基因型×盐度互作表明存在可资利用的遗传多样性。STI×YSI复合指数成功解决了基线产量与毒物兴奋效应间的悖论，推荐优先考虑具有高STI和高基线产量的基因型（如Azarshahr、Maragheh、Ajabshir，Malekan为次选）。APX活性增加是关键的抗氧化防御机制，而SOD和脯氨酸的无效结果因低统计功效而不具结论性，需增加重复数和多次采样。香叶醇的双相响应支持生长-防御权衡理论，适度盐胁迫可提高香叶醇含量，但严重胁迫降低产量。与已有研究相比，优良香蜂花基因型的耐盐指数与商业薄荷、罗勒相当甚至更优。研究局限性包括：缺乏Na⁺/K⁺测量、SOD和脯氨酸统计功效不足、仅温室单季实验、基因型数量有限。未来需在田间条件下进行多地点多年验证，增加重复数，并整合转录组/蛋白质组分析。

**研究结论**
本研究首次系统筛选了10个伊朗香蜂花基因型对盐胁迫的响应，引入了新型STI×YSI复合指数框架以解决基线产量-毒物兴奋效应悖论。盐度显著降低农艺性状（株高降低45%，叶数减少61%，地上部干重减少52%），APX活性显著增加（从0.386增至0.517单位，p<0.001），表明关键抗氧化防御机制。STI×YSI复合指数成功区分了高基线产量与毒物兴奋效应，Bootstrap分析确认优良基因型排名高度稳定（94.2%）。基于该指数，经PCA（KMO=0.76，解释83.1%方差）和聚类分析（共表相关系数=0.82）验证，Azarshahr（0.70）、Maragheh（0.64）和Ajabshir（0.57）被识别为有前景候选者，Malekan（0.52）因高基线产量（22.14 g）作为次选。这些基因型相对于敏感基因型地上部干重高13.08%，香叶醇高15.50%，GPX活性高65.01%。然而，本研究结果仅为初步筛选，存在关键局限：（1）缺乏Na⁺和K⁺测量，无法解释离子耐受机制；（2）SOD（0.45）和脯氨酸（0.38）统计功效低，结果不具结论性；（3）仅在温室单季实验，需田间验证。最终结论：Azarshahr、Maragheh和Ajabshir为有前景候选者，Malekan为次选，为后续研究提供了明确方向。但在进行育种推荐之前，必须通过田间实验（测量Na⁺/K⁺比值，增加重复数n≥5）进行确认。