盐胁迫下黄河流域棉花品种主要性状的通径分析及优异种质筛选.pdf

櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄櫄 马保国酵素菌肥对大蒜的增产效应及培肥土壤效果初探 中国农学通报，（）：韩剑宏，刘泽霞，张连科，等环保酵素对盐碱土关键化学性质的影响 水土保持通报，（）：郜胜涛，黄荣松，张涛，等微生物菌肥替代部分化肥施用对土壤肥力及水稻产量的影响 云南农业科技，（）：，：，（）：高亮，谭德星酵素菌生物有机肥在潍坊滨海盐土上的应用效果研究 现代农业科技，（）：，翟敬华，蒋细旺酵素对红菜薹生长发育及土壤特性的影响 长江大学学报（自然科学版），（）：，杜磊，雄东毕，李成奇，等盐胁迫下黄河流域棉花品种主要性状的通径分析及优异种质筛选 江苏农业科学，（）：盐胁迫下黄河流域棉花品种主要性状的通径分析及优异种质筛选杜磊，雄东毕，李成奇，赵海红（运城学院生命科学系，山西运城 ；新疆阿克苏地区农业农村局，新疆阿克苏 ）摘要：面对黄河流域盐碱地开发利用重大政策机遇期，选育耐盐碱新品种有利于促进盐碱地棉花产业发展。在 溶液模拟盐胁迫和大田自然盐碱胁迫条件下，对 份棉花种质资源进行萌发期和苗期盐胁迫试验。综合考虑萌发期和苗期的耐盐情况，通过因子分析和聚类分析，筛选出 份耐盐材料和 份盐敏材料。通过建立逐步回归方程得到耐盐综合得分最优预估方程 预估 ，筛选出相对地下干物质量、相对发芽率、相对出苗率和相对 酶活性等 个能显著体现棉花耐盐表现的指标。通径分析结果显示，相对地下干物质量的直接通径系数最大，表明相对地下干物质量对耐盐综合得分的直接作用最大；相对地下干物质量的决策系数最大，为 且大于 ，表明其对耐盐综合得分的影响贡献最大且起促进作用。因此，选育棉花耐盐品种，要注重对相对地下干物质量的选择，同时兼顾对相对发芽率、相对出苗率和相对 酶活性的选择。关键词：棉花；盐胁迫；因子分析；聚类分析；通径分析中图分类号：文献标志码：文章编号：（）收稿日期：基金项目：新疆生产建设兵团重点研发计划（编号：）；山西省高等学校科技创新项目（编号：）；运城学院博士科研启动项目（编号：）；运城学院特色农产品发展学科群科研项目（编号：）。作者简介：杜磊（），男，山西运城人，博士，副教授，主要从事棉花抗逆育种研究。：。黄河流域棉区作为我国三大主要棉区之一，近年来随着粮棉冲突加剧，棉花种植面积逐年减少，种植区域不断向盐碱地、旱地转移 。黄河流域盐碱地面积约 万 ，约占我国盐碱地面积的五分之一。党的十八大以来，在习近平总书记的亲自谋划和部署下，黄河流域生态保护和高质量发展上升为重大国家战略，至此黄河流域盐碱地开发利用迎来重大政策机遇期。棉花作为中等耐盐作物，选育耐盐碱新品种有利于促进盐碱地棉花产业发展。棉花耐盐性在不同生长期和品种间表现不同，其中萌发期和苗期的耐盐能力相对较弱，是鉴定棉花品种耐盐性的关键时期。同一棉花品种的耐盐性在萌发期和苗期存在差异，只有对 个时期的耐盐性进行综合分析，才能更准确地评价其耐盐性。黄河流域棉花种质资源极其丰富，很多传统品种性状优良但以早熟、抗病、高产最为突出，育种工作中对棉花的耐盐碱性相对重视不够 。统计方法中的因子分析和通径分析已应用于棉花、玉米、葡萄、水稻等作物的综合评价 。基于因子分析和通径分析，本研究以黄河流域培育推广的棉花品种群体为试验对象，通过综合分析萌发期和苗期的耐盐性，比较不同棉花品种间的耐盐性差异，并对耐盐性状进行鉴定筛选，以期为棉花耐盐育种工作提供数据参考。江苏农业科学 年第 卷第 期 材料与方法 试验材料棉花材料由山西省农业科学院棉花研究所和南京农业大学棉花研究所提供，共计 份来自黄河流域的棉花品种。所有材料均经过多代自交。试验方法 萌发期耐盐鉴定采用双层滤纸法进行萌发期耐盐鉴定。空白对照采用蒸馏水溶液处理，盐胁迫采用 的 溶液处理。每处理 粒种子，重复 次。分别测定第 天、第 天的发芽数目，第 天的胚根长度 。苗期耐盐鉴定试验于 连续 年在山西省运城市盐湖区解州镇许贾村盐碱地进行，通过 年检测棉花生长情况和土壤盐含量，确定了 块棉田，在试验当年雨水偏少的 月其土壤相关指标基本稳定，年平均值如 下：总 盐 含 量 为 ，速 效 磷 为 ，有机质含量为 ，碱解氮含量为 ，速效钾含量为 ，电导率为 ，值为 。选取正常农业生产地块作为对照，年连续 年平均值如下：土壤总盐含量为 ，速效磷含量为 ，有机质含量为 ，碱解氮含量为 ，速效钾含量为 ，电导率为 ，值为 。每处理重复 次，单个棉花品种为 个处理。田间试验采用随机区组试验，田间一膜 行种植，每份材料种植 行，行长为 ，株距为 ，平均行距为 。田间管理同大田常规管理，播种 之后，测量相关指标并取样进行有关分析。指标测定分别测量空白对照和盐胁迫下各性状的数值，再计算各性状的相对值。本试验选取了相对发芽势（）、相对发芽率（）、相对胚根长（）等 个萌发期耐盐性状相对值和相对叶绿素（）、相 对 苗 高（）、相 对 地 下 干 物 质 量（）、相对地上干物质量（）、相对 酶活性（）、相对 酶活性（）、相对 酶活性（）、相对 含量（）、相对根冠比（）、相对出苗率（）和相对立苗率（）等 个苗期的耐盐性状相对值。发芽情况：以胚根长 种子长作为发芽标准，统计发芽势、发芽率；发芽势：发芽数 供试种子数 ；发芽率：发芽数 供试种子数 ；胚根长：以所有已发芽种子第 天的平均胚根长为指标，用 根系分析软件获取根长；叶绿素含量：采用 型测定仪，测定棉株倒 叶中部；苗高：采用直尺测量棉株从子叶节到主茎生长点的长度；干质量：将棉株轻柔的冲洗干净，分离棉株地上部和地下部，杀青 ，烘干至恒质量称质量；、活性及 含量测定：取棉株倒 叶测量；根冠比：地下部鲜质量 地上部鲜质量；出苗率：出苗种子数 供试种子数 ；立苗率：幼苗 叶 心时测量，立苗种子数 供试种子数 ；相对值：盐胁迫处理数值 对照处理数值；变异系数：标准差 平均值 。数据处理采用 和 对试验数据进行相关性分析、因子分析、聚类分析、正态性检验、逐步回归分析、间接通径系数和决策系数分析。间接通径系数、决策系数计算方法参考王巧妹等的方法 。结果与分析 耐盐性状的表型多样性分析对 份黄河流域棉花品种的 个耐盐性状进行研究，相关描述统计结果见表 。相关研究显示，含量和根冠比同棉花耐盐性负相关，其他 个性状同棉花耐盐性正相关 。由表 可知，含量和根冠比在盐胁迫下所得数据的均值同对照组相比，依次增加了 和 ，胚根长、酶活性、地上干物质量、苗高、酶活性、酶活性、发芽势、叶绿素含量、地下干物质量、发芽率、立苗率和出苗率在盐胁迫下所得数据的均值同对照组相比，依次降低了 、和 。盐胁迫下胚根长降幅最大，预示着盐胁迫对胚根生长存在明显抑制作用。发芽率、立苗率和出苗率在盐胁迫下所得数据的均值同对照组相比变化较小，表江苏农业科学 年第 卷第 期书书书表 耐盐性状的描述统计处理项目发芽势（）发芽率（）胚根长（）叶绿素含量（值）苗高（）地下干物质量（）地上干物质量（）活性（，）活性 （）活性 （）含量（）根冠比出苗率（）立苗率（）对照组最小值 最大值 均值 标准差 变异系数 盐胁迫组最小值 最大值 均值 标准差 变异系数 明作为耐盐先锋植物，棉花的耐盐性相对较强，在盐胁迫下其发芽和出苗立苗所受影响较小。根冠比在盐胁迫下的均值比对照组大，地上干物质量和地下干物质量在盐胁迫下所得数据的均值同对照组相比分别降低了 和 ，上述结果均表明盐胁迫对地上部分的抑制作用大于地下部分。对于变异系数，由表 可知，对照组和盐胁迫时大部分性状的变异系数绝对值较大，表明本研究所涉及的 份棉花品种的耐盐性存在较大差异，预示着可以筛选到适宜的耐盐品种和盐敏品种。耐盐性状的因子分析因子分析需要各性状间具有较强的相关性。各耐盐性状的 （）检测和 球度检测分析结果见表 。值为 且大于 ；球度检测值为 且大于 ，对应的伴随概率为 且 ，表明变量间存在相关性，可以进行因子分析和公共因子提取 。表 和 的检验检验类型指标数值 检验 值 球形度检验近似卡方 通过因子分析从 个耐盐性状中提取出 个公共因子，结果见表 。特征值分别为 、和 ，均大于 。个公共因 子 的 贡 献 率 分 别 为 、和 ，累计贡献率为 ，放弃的其他 个耐盐性状的累计贡献率不到 ，信息丢失较少。表 特征向量分析显示：公共因子 与相对叶绿素、相对苗高、相对地下干物质量和相对地上干物质量相关，公共因子 与相对发芽势、相对发芽率相关，公共因子 与相对出苗率、相对立苗率相关，公共因子 与相对 含量、相对根冠比相关，公共因子 与相对 酶活性、相对 酶活性、相对 酶活性和相对胚根长相关。棉花品种耐盐性综合评价个公共因子的得分函数计算公式，依据得分系数矩阵（表 ）整理如下：公共因子得分函数 江苏农业科学 年第 卷第 期表 各公共因子的贡献率和特征向量项目特征向量因子 因子 因子 因子 因子 特征值 贡献率 累计贡献率 注：表示耐盐性状在各因子中的最大绝对值。表 公共因子得分系数矩阵性状得分系数因子 因子 因子 因子 因子 （）（）（）（）（）（）（）（）（）（）（）（）（）（）；公共 因 子 得 分 函 数 ；公共因子 得分函数 ；公共 因 子 得 分 函 数 ；公共因子 得分函数 。将 个公共因子的贡献率和累计贡献率作为权数进行加权分析，可得各棉花品种的耐盐综合得分计算公式：（）。基于因子分析获得各棉花品种的耐盐综合得分，分值越大，品种的耐盐性越高。基于系统分析法的组间联接平均法进行聚类分析，将 份棉花品种 分 为 类（图 ），类 品 种 嘉 星 号（）、鲁棉 号（）、国欣棉 号（）、邯 （）、冀 邯 号（）和中棉所 号（）等 个为盐敏材料，耐盐综合得分范围为 ，可知棉花品种的耐盐综合得分小于 ，可认为是盐敏材料。类品种豫棉 号（）、陕棉 号（）、邯郸 （）、中棉所江苏农业科学 年第 卷第 期 号（）、冀棉 （）、中棉所 号（）、中 棉 所 号（）、冀棉 号（）、中棉所 号（）、豫棉 号（）、豫早 （）和仁和 号（）等 个为耐盐材料，耐盐综合得分范围为 ，可知棉花品种的耐盐综合得分大于 ，可认为是耐盐材料。类品种为较耐盐材料，类品种为不耐盐材料。耐盐综合得分与耐盐性状的逐步回归分析通过 的单样本 检验对 个棉花品种的耐盐综合得分进行正态性检验，结果见表 ，正态性检验显著性水平是 ，耐盐综合得分的数值分布属于正态分布。对盐胁迫下 个棉花品种的耐盐综合得分与 个耐盐性状的关系进行逐步回归分析 。以耐盐综合得分为因变量，个耐盐性状等为自变量，通过 进行逐步回归分析，输出结果见表 、表 。逐步回归分析结果（表 ）显示，相关系数 和决定系数 随着自变量依次加入方程而逐步增大，说明引入的自变量对因变量的作用在增加。模型 表 正态性检验正态性检验耐盐综合得分 正态参数 ，均值标准差 最极端差别绝对值 正 负 渐近显著性（双侧）的相关系数 和决定系数 分别为 和 ，剩余因子 槡的数值较小，说明模型较全面地考虑了对耐盐综合得分有影响的自变江苏农业科学 年第 卷第 期表 模型汇总模型调整的 标准估计的误差 表 各模型变量的系数模型自变量非标准化系数（）标准误差标准系数（试用版）值值（常量）（常量）（常量）（常量）（）（）（）（）注：为 的缩写，代表回归系数；为对回归系数的显著性检验。量 。依据表 模型 的回归系数，建立耐盐综合得分与各耐盐性状的线性回归方程，得到耐盐综合得分最优预估方程 预估 。显著性分析结果表明，相对地下干物质量、相对发芽率、相对出苗率和相对 酶活性的 值均小于 ，表明自变量与因变量间差异显著，上述线性回归方程有统计学意义。相对地下干物质量、相对发芽率、相对出苗率和相对 酶活性可以作为关键指标，在棉花耐盐鉴定时优先测定并代入上述线性回归方程，对各棉花品种耐盐性进行快速预估鉴定。耐盐综合得分与各耐盐性状的间接通径系数和决策系数计算耐盐综合得分与各耐盐性状的相关性分析结果见表 。可知 与 、的相关系数分别为 、，与 、的相关系数分别为 、，与 的相关系数为 。、的直接通径系数分别为 、（表 ）。、分别对 的间接通径系数具体计算如下：书书书表 相关性及检验输出结果模型相关系数 注：、分别表示在 和 水平上差异显著。江苏农业科学 年第 卷第 期 ，；，（）；，。其他间接通径系数计算结果见表 。结果表明，的直接通径系数最大，表明相对地下干物质量对耐盐综合得分的直接作用最大。的间接通径系数最大，表明相对 酶活性的间接作用最大，即其通过相对地下干物质量、相对发芽率、相对出苗率作用于耐盐综合得分的间接作用最大。决定系数作为通径分析中的决策指标，可以把自变量对响应变量的作用进行综合排序 。其计算结果见表 ，可知相对地下干物质量的决定系数最大，为 且大于 ，表明其对耐盐综合得分的影响贡献最大且起促进作用。相对 酶活性决定系数最小，为 且大于 ，表明其对耐盐综合得分的影响贡献最小且起促进作用。表 耐盐综合得分与各耐盐性状的通径分析自变量与 值的简单相关系数直接通径系数间接通径系数 合计决定系数 讨论与结论 棉花的耐盐性存在器官差异棉花的各个器官中，根部最直接、最先感受盐胁迫。龚江等研究发现，棉花可通过促进根系的生长应对盐胁迫 。黄雅婕等在生产实践中发现耐盐棉花品种最直观的表现是地下根系相对发达 。本研究通径分析结果显示，相对地下干物质量的直接通径系数最大，表明其对棉花品种耐盐性的直接作用最大；相对地下干物质量的决定系数最大且为正数，表明其对棉花品种耐盐性的影响贡献最大且起促进作用。前人研究结果也表明，棉花的根系同其耐盐性具有很高相关性 。但这些研究缺乏针对根系的从形态特征到生理生化特征的全面系统阐述，这也是后续研究的重点。超氧化物歧化酶（）增强棉花耐盐能力的途径研究表明，过量表达 基因，可提高植物的耐盐性。等研究发现，随着盐浓度的升高，耐盐棉花品种 酶活性迅速上升，酶活性变化不大；非耐盐棉花品种 酶活性变化不大，酶活性呈下降趋势 。赵远伟等研究发现棉花可通过提高 、和 活性来提升耐盐性。植物耐盐调控机制研究表明，盐胁迫会导致植株体内活性氧的过量生成，而抗氧化酶在活性氧清除过程中发挥重要作用 。本研究结果为 增强棉花耐盐能力增加了一种新的解释。棉花的耐盐性属于数量性状，除了受耐盐基因直接作用外，还受基因互作、多基因累加效应等的间接作用 。本研究通径分析结果显示，相对 酶活性的间接通径系数最大，表明其对棉花品种耐盐能力的间接作用最大，即 基因通过作用于其他耐盐性状基因，间接影响棉花品种耐盐能力的作用最大；决策系数分析结果显示，相对 酶活性决策系数大于 ，表明其对棉花品种耐盐性起促进作用，即升高 酶活性可增强棉花耐盐能力。耐盐综合得分最优预估方程在实践应用中的便捷性、可靠性和局限性通过逐步回归分析，得到耐盐综合得分最优预估方程 预估 。生产实践中仅需测定相对地下干物质量、相对发芽率、相对出苗率和相对 酶活性等 个性状就可预估棉花品种耐盐性，可大大提高工作效率。棉花的耐盐性属于复杂数量性状，受多基因调控多因素影响，且不同的发育时期耐盐性有差异，因此采用多个时期多个指标进行耐盐性评价更具有可靠性。上述预估方程中的 个性状，包含了萌发期和苗期的性状，包含了生理和生化性状，具备多时期多指标的特点，因此其预估结果是可靠的。本研究得到的 个耐盐材料中，上述 个性状的数值均明显高于盐敏材料。使用该方程预估棉花品种耐盐性时，可通过比较预估得分和因子分析江苏农业科学 年第 卷第 期得分获得预估精度，发现耐盐材料平均预估精度为 ，盐敏材料平均预估精度为 。表明在生产实践中，使用该方程对棉花耐盐性进行预估时，预估的分值越高，预估精度也越高；反之，预估的分值越低，预估精度也越低。即棉花品种自身的耐盐性越好，使用该方程进行耐盐性预估的准确度越高，反之则越低。通过因子分析和聚类分析，获得各棉花品种的耐盐综合得分，从 份黄河流域棉花品种中筛选到 份耐盐材料和 份盐敏材料，耐盐综合得分大于 ，认为是耐盐材料；小于 认为是盐敏材料。通过通径分析，获得可依据分值高低预估棉花品种耐盐性的方程：预估 ，依据方程可知棉花品种耐盐性同相对地下干物质量、相对发芽率、相对出苗率和相对 酶活性等个性状关系最密切。其中相对地下干物质量对棉花品种耐盐性的直接影响力最大，表明棉花根系同其耐盐性具有很高相关性；相对 酶活性通过影响其他性状来间接影响棉花品种耐盐性的促进作用最大，进一步证实了高 酶活性对棉花品种耐盐性起促进作用。参考文献：董合忠中国棉花种业和原棉品质的国际竞争力分析 中国棉花，（）：，牛铭铭，王欢，王春琳，等沿黄盐碱地美国红鱼池塘养殖技术 河南水产，（）：，朱建峰，杨秀艳，武海雯，等植物种子萌发期耐盐碱性提高技术研究进展 生物技术通报，（）：胡根海，付远志，张晓红耐旱耐盐棉花种子萌发的影响因子分析初报 种子，（）：张莹莹，卢道文，孙海潮，等黄淮海区域玉米收获期籽粒含水量与主要农艺性状的相关分析和通径分析 江苏农业科学，（）：和雅妮，奚晓军，查倩，等葡萄果实矿质元素与品质的相关性及通径分析 江苏农业科学，（）：张进帅，黄翠红，王慧，等水稻单株产量相关性状的相关性及其通径分析 江苏农业科学，（）：费克伟，罗晓丽，司怀军，等五个棉花品种抗旱性与 活性相关性分析 作物杂志，（）：杜磊中国陆地棉品种耐盐性评价及耐盐优异等位变异的挖掘 南京：南京农业大学，：李志博，曹连莆，魏亦农，等荒漠绿洲区棉花幼苗耐寒生理特性及其评价指标分析 棉花学报，（）：黄雅婕陆地棉耐盐碱评价体系的建立 乌鲁木齐：新疆农业大学，：王巧妹，王继武，管志勇，等茶用菊七月白 苏菊 号 代群体与产量相关性状的变异分析及高产植株筛选 江苏农业学报，（）：宋小园，朱仲元，刘艳伟，等通径分析在 逐步线性回归中的实现 干旱区研究，（）：余如刚，王雪茄，王国良，等紫花苜蓿品种苗期耐盐性分析及评价指标筛选 草地学报，（）：张国伟，路海玲，张雷，等棉花萌发期和苗期耐盐性评价及耐盐指标筛选 应用生态学报，（）：龚江，谢海霞盐分胁迫下滴灌时间对棉花生长和产量的影响 中国棉花，（）：，马瑛瑛钾钙调控对盐、旱胁迫条件下棉株生长、生理及产量的影响 杨凌：西北农林科技大学，：王庆惠，韩伟，侯银莹，等不同耐盐品种棉花根系主要指标对盐分胁迫的响应 应用生态学报，（）：刘，钱玉源，张海娜，等 个陆地棉品种（系）苗期耐盐性鉴定评价 作物研究，（）：张超，马晓丽，卢晓峰，等盐分胁迫下土施甲哌对棉苗叶片生理和根系形态的影响 江苏农业科学，（）：李函利，王省芬，王国宁，等低酚棉种质资源耐盐性鉴定 植物遗传资源学报，（）：，（）：，（），：，（）：赵远伟棉花幼苗根系对盐胁迫的响应及机制 保定：河北农业大学，：，（），（）：，：，（）：张晓佩陆地棉长链非编码 和 的功能和作用机制研究 泰安：山东农业大学，：江苏农业科学 年第 卷第 期