基于隶属函数法对甘蓝型油菜种质苗期耐低氮胁迫能力的综合评价.pdf

陈敬东，余忆，鲁金春子，等基于隶属函数法对甘蓝型油菜种质苗期耐低氮胁迫能力的综合评价 江苏农业科学，（）：基于隶属函数法对甘蓝型油菜种质苗期耐低氮胁迫能力的综合评价陈敬东，余忆，鲁金春子，戴希刚，沈金雄，傅廷栋，伍晓明，曾长立，万何平（江汉大学生命科学学院 湖北省汉江流域特色生物资源保护开发与利用工程技术研究中心，湖北武汉 ；华中农业大学，湖北武汉 ；中国农业科学院油料作物研究所 农业农村部油料作物生物学与遗传育种重点实验室，湖北武汉 ）摘要：油菜作为我国重要的油料作物，其生长发育过程中对氮素营养需求量较大，但氮肥利用率却较低，在很大程度上制约了油菜的正常生长。探究苗期甘蓝型油菜耐低氮胁迫能力的评价方法，筛选耐低氮胁迫能力强的种质，可以减少氮肥施用量，保护生态环境。以 份甘蓝型油菜种质材料通过水培试验，考察低氮水平（）和正常氮水平（）的地上部鲜质量、地下部鲜质量、地上部干质量、地下部干质量、叶片数、株高、总根长、根平均直径、侧根数等性状，以各性状的耐低氮系数为衡量油菜耐低氮胁迫能力的指标，利用隶属函数法对不同基因型油菜进行耐低氮胁迫能力的综合评价，筛选出低氮胁迫耐受型材料个，低氮胁迫敏感型材料 个，低氮胁迫中间型材料 个。值与地上部鲜质量、地下部鲜质量、地上部干质量、地下部干质量、叶片数、株高、总根长和侧根数呈显著正相关，可作为油菜苗期耐低氮胁迫能力鉴定的辅助指标，对探究油菜耐低氮胁迫遗传机制，选育耐低氮胁迫油菜新种质具有重要意义。关键词：甘蓝型油菜；耐低氮胁迫能力；隶属函数法；苗期；氮肥利用率中图分类号：文献标志码：文章编号：（）收稿日期：基金项目：国家重点研发计划（编号：）；湖北省汉江流域特色生物资源保护开发与利用工程技术研究中心开放基金（编号：）。作者简介：陈敬东（），男，湖北红安人，硕士，从事作物遗传育种研究。：。通信作者：曾长立，博士，教授，从事油菜种质资源研究，：；万何平，博士，讲师，从事油菜遗传研究，：。甘蓝型油菜（）因产量较高、适应能力较强，在我国油菜种植面积中占比达到 左右 。年我国自产植物油占植物油消费总量的 ，国产油料难以满足国民生产消费需求，供需不平衡的现象与日俱增 。因此，选育优质油菜种质，提高油料自给率，对保证我国油料安全具有重要意义。在油菜生长发育过程中，氮（）作为大量元素，对于其产量、品质具有重要影响 。油菜处于苗期时，充足的氮素营养保证了其正常的生长发育，但在生产实践过程中，成熟期却常常因过量施用氮肥，导致油菜大量吸收氮素营养，造成其体内营养元素平衡性被破坏，反而抑制了油菜的正常生长 。与其他作物相比，油菜对氮素营养需求量相对较大，然而其氮肥利用效率却相对较低，生产者需要大量施用氮肥以保证其产量 。同时，过量施用氮肥会使盈余的氮素营养通过多种途径流失，从而可能导致一系列的环境问题 。除作为油料作物外，油菜已开发出饲用、菜用、观赏等多种用途。江西婺源、山西汉中等地以油菜花海为依托，打造了一批享誉全国的农业观光旅游景区，为推动当地经济社会发展提供了新动力 。因此，须要筛选出一批耐低氮胁迫能力强、适应性广的油菜品种（品系），以节约成本、保护生态环境。作物的耐低氮胁迫能力是一个复杂的数量遗传性状，可能受到多种功能基因的共同调控，且对环境因子的变化十分敏感，难以直接评价 。目前，关于油菜低氮胁迫耐受种质和耐低氮胁迫机理研究已经取得了一些进展。研究者们通过选择低氮耐受的种质进一步改良油菜氮效率（）。等研究结果表明，半矮秆杂交种在多种条件下都表现出更高的氮收获指数和氮利用效率，且通过试验确定其在低氮条件下为氮高效基因型 。邹小云等通过对抽薹期不同氮效率油菜的各性状进行分析，发现油菜抽薹期氮素吸收效率由根平均直径、最大叶长和叶宽的乘积、叶片的蒸腾速率以江苏农业科学 年第 卷第 期及油菜根系硝酸还原酶的活性共同决定 。等发现在高氮和低氮 种条件下，根系形态（）良好或氮吸收效率（）较高的种质，其 也相对较高，且 对 的影响比 对 的影响更加明显 。上述研究通常比较多个农艺性状，并从中筛选出某项效应最大的单一指标从而评价油菜种质耐低氮胁迫能力，然而油菜的耐低氮胁迫能力是一个复杂的性状，不同低氮耐受型材料的耐性机制也不尽相同，利用单一指标筛选耐受材料的方法具有一定的局限性，其鉴定结果的有效性尚需进一步试验验证。而主成分分析法可以准确地计算多个指标所占权重，并从多个具有一定相关性的指标中计算出新的综合变量，具有较好的准确性。本研究在前人基础上，利用正常氮（）和低氮（）种氮素营养水平，通过对 份油菜材料 个指标的耐低氮胁迫能力的综合评价，分析比较了综合指标与油菜苗期生物量之间的关系，以验证低氮耐受指标的准确性，并对各指标进行评价，确定了油菜耐低氮能力评价的适宜指标，筛选获得了一批低氮耐受材料，以期为油菜耐低氮品种选育提供理论支持。材料与方法 材料供试甘蓝型油菜自交系材料共 份（表 ），均由华中农业大学国家油菜工程技术研究中心提供，份油菜资源主要来自于湖北、湖南、四川、重庆等省（市），前期研究结果表明，该群体材料具有广泛的遗传变异 。供试材料在 年 月 日播种于湖北省武汉市蔡甸区江汉大学油菜种质资源试验区内。试验设计与处理为了探索能维持油菜最低生长要求的氮素营养水平，随机挑选 个油菜材料，利用油菜水培方法培养至 叶期进行不同浓度氮营养处理（浓度：、）。根据结果，选取 浓度为低氮胁迫浓度，浓度为正常氮营养水平。为精确考察 份油菜资源的耐低氮能力，等前期建立了一套大规模油菜水培系统 ，可以同时培养超过 株油菜幼苗。在 份油菜自交系中每个基因型选出 粒种子进行发芽处理，再从中选出 株长势相似的健康幼苗，移入水培系统中培养 周 。每周更换 次培养液，周培养液分别为 霍格兰溶液、霍格兰溶液、霍格兰溶液和 霍格兰溶液。从第 周开始，将同种基因型 株幼苗均分为 组，分别以正常氮（）水平（）和低氮（）水平（）对油菜进行处理，营养液配方见表 。配方中，调节 和 （）浓度控制正常氮营养液和低氮营养液的氮含量，因 和 （）浓度导致的 和 含量的改变则利用 和 补足。以后每 为种条件更换次培养液，培养周后收获并进行表型精准鉴定。本次利用无重复随机区组设计试验。表型鉴定材料收获后分成地上部（）、地下部（），考察其叶片数（）、株高（）、鲜质量（、）。取地下部根系，去除表面杂质后用去离子水多次冲洗，用爱普生扫描仪（）对甘蓝型油菜 根 系 进 行 扫 描，通 过 根 系 形 态 分 析 软 件（）分析获得总根长（）、根平均直径（）、侧根数（）等油菜根系形态指标。地上部干质量（）和地下部干质量（）采用 杀青 后，烘干至恒质量测定其数值。指标计算及数据分析低氮胁迫耐受指数：为了排除因基因型背景差异而对试验结果可能造成的影响，选择计算油菜 个表型性状的低氮胁迫耐受指数（、和 ，为低氮条件下的表型性状数值与正常氮条件下表型性状数值之间的比值）。隶属函数值公式如下：（），。式中：表示第 个综合指标；表示第 个综合指标中的最大值；表示第 个综合指标中的最小值。权重公式如下：，。式中：表示第 个综合指标在全部综合指标中的重要程度，即权重；为各基因型油菜第 个综合指江苏农业科学 年第 卷第 期表 份甘蓝型油菜自交系材料构成编号材料名称来源类型编号材料名称来源类型 湖北半冬性 湖北半冬性 贵州半冬性 湖北半冬性 重庆半冬性 湖北半冬性 四川半冬性 湖北半冬性 四川半冬性 湖北半冬性 湖南半冬性 湖北半冬性 湖南半冬性 安徽半冬性 湖南半冬性 青海春性 湖南半冬性 青海春性 湖南半冬性 陕西冬性 湖南半冬性 陕西半冬性 湖南半冬性 陕西冬性 湖北半冬性 陕西半冬性 湖北半冬性 湖北半冬性 湖北半冬性 湖北半冬性 湖北半冬性 重庆半冬性 湖北半冬性 重庆半冬性 湖北半冬性 重庆半冬性 未知半冬性 重庆半冬性 湖北半冬性 重庆半冬性 湖北半冬性 重庆半冬性 湖北半冬性 湖北半冬性 湖北半冬性 湖北半冬性 江苏半冬性 湖北半冬性 湖北半冬性表 不同氮素水平营养液配方化学试剂名称正常氮处理（）低氮处理（）（）标的贡献率。综合评价指标公式如下：（），。式中：为各基因型在低氮胁迫条件下的耐低氮胁迫能力综合评价指标。上述公式参考文献 。采用 进行数据处理，采用 绘图。结果与分析 油菜幼苗在不同浓度氮营养水平下的表现由图 可知，不同氮营养水平下油菜幼苗的长势呈现明显的变化。其中高氮水平下，油菜生长正常，叶片颜色呈绿色，而随着氮浓度的降低，油菜生长会因为氮素营养的逐渐缺乏，而受到不同程度的抑制，具体表现为叶片变小、颜色变黄，出现紫色叶片。在低氮浓度下，油菜幼苗生长极其缓慢。低氮水平对不同油菜品种的抑制效果存在一定的差异，说明不同基因型油菜的耐低氮能力不同。综合考虑，选取 浓度为低氮胁迫浓度，浓江苏农业科学 年第 卷第 期度为正常氮营养水平。不同氮水平条件下 份甘蓝型油菜材料表型分析及相关性分析由表 可知，在正常氮条件下，甘蓝型油菜植株的 平均值为 ，变动范围为 ，变异系数（）为 ；的平均值为 ，变 动 范 围 为 ，为 ；的平均值为 ，变动范围为 ，为 ；的平均值为 ，变 动 范 围 为 ，为 ；的平均值为 张，变动范围为 张，为 ；的平均值为 ，变动范围为 ，为 ；的平均值为 ，变动范围为 ，为 ；的平均值为 ，变动范围为 ，为 ；的平均值为 条，变动范围为 条，为 。在低氮条件下，甘蓝型油菜植株的 平均值为 ，变动范围为 ，为 ；的平均值为 ，变动范围为 ，为 ；的平均值为 ，变 动 范 围 为 ，为 ；的平均值为 ，变动范围为 ，为 ；的平均值为 张，变动范围为 张，为 ；的平均值为 ，变动范围为 ，为 ；的平均值为 ，变动范围为 ，为 ；的平均值为 ，变动范围为 ，为 ；的平均值为 条，变动范围为 条，为 。、和 在低氮条件下的表型数值都极显著低于正常氮条件下的表型数值，而 和 在低氮条件下的表型数值都极显著高于正常氮条件下的表型数值。在 和 种氮浓度下，个农艺性状都呈现出正态分布的趋势（图 ），表明这些性状可能均为数量性状。个低氮胁迫耐受指数（）均表现出较大的江苏农业科学 年第 卷第 期表 份甘蓝型油菜种质在 种氮水平下的性状统计性状（）（）（）（）（张）正常低 正常低 正常低 正常低 正常低 均值 最小值 最大值 标准差 变异系数（）性状 （）（）（）（条）正常低 正常低 正常低 正常低 均值 最小值 最大值 标准差 变异系数（）注：表示在 水平上差异显著。遗传变异，的平均值为 ，变动范围为 ，为 ；的平均值为 ，变动范围为 ，为 ；的平均值为 ，变动范围为 ，为 ；的平均值为 ，变动范围为 ，为 ；的平均值为 ，变动范围为 ，为 ；的平均值为 ，变 动 范 围 为 ，为 ；的平均值为 ，变动范围为 ，为 ；的平均值为 ，变 动 范 围 为 ，为 ；的平均值为 ，变动范围为 ，为 。结果表明，低氮胁迫条件对不同基因型油菜的 、和 存在显著的抑制效果，其中 受影响最大。同时，低氮胁迫条件也会显著促进 和 。此外，个低氮胁迫耐受指数都有正态分布的特点（图 ）。相关性分析结果（图 ）显示，各性状的低氮胁迫耐受指数之间还存在着一定的关联性，这导致各指标所蕴含的信息可能存在交叉重复，且各性状在油菜低氮胁迫耐受能力中发挥的作用也不完全一致，因此，直接以低氮胁迫耐受指数评价油菜品种耐低氮胁迫能力的方法是不准确的。主成分分析利用 软件对 个低氮胁迫耐受指数进行主成分分析（表），前个主成分（综合指标 、）的 贡 献 率 分 别 达 到 了 、，累计贡献率合计为 ，其他可忽略。这样就将原本 个低氮胁迫耐受指数转化为了 个相互独立的综合指标。决定第 综合指标大小的低氮胁迫耐受指数主要是、和 ，可描述为生物量因子、根系形态因子；决定第 综合指标大小的主要是 和 ，可描述为根系形态因子；决定第 综合指标大小的主要是 和 ，可描述为植株形态因子、根系形态因子；决定第 综合指标大小的主要是 ，可描述为植株形态因子。油菜耐低氮胁迫能力的综合评价和聚类分析利用公式计算不同基因型油菜所有综合指标的隶属函数、权重和耐低氮胁迫能力综合评价 值（表）。以 值作为衡量不同基因型油菜耐低氮胁迫能 力 的 强 弱，其 中 材 料 （）和 （）的 值较大，具有较强的耐低氮胁迫能力。采用欧氏距离（）和最大距离法（）对 值进行聚类分析（图 ），将 个基因型根据耐低氮胁迫能力分为 类，以综合衡量不同甘蓝型油菜种质的耐低氮胁迫能力。其中 份材料被划分为低氮胁迫耐受型种质（表 ），其 值均大于 ；份材料被划分为低氮胁迫敏感型种质，其 值均小于 ；其他 份材料则被划分为低氮胁迫中间型种质，其 值处于 之间。江苏农业科学 年第 卷第 期 线性回归预测油菜耐低氮胁迫能力主成分分析结果表明，生物量因子和根系形态因子中 的贡献率相对较高。因此以 、和 作为自变量，以 值作为因变量，建立回归方程：（，）。该方程展示出 个低氮胁迫耐受指数和 值之间的显著线性关系，简化了性状鉴定工作。与基于单一性状的评价方法相比，该方法准确性高、筛选效率好。讨论与结论 讨论作物的耐低氮胁迫能力被广泛认为是受多基因调控的数量性状，每个相关的功能基因所发挥的江苏农业科学 年第 卷第 期表 各综合指标的系数及贡献率综合指标 贡献率（）作用可能相对微小；且功能基因与环境因子的互作会直接影响作物耐低氮胁迫能力的表现，这使得耐低氮胁迫能力性状不稳定，可能受到各种因素的影响 。因此，选择正确的耐低氮胁迫能力评价方法是最终准确衡量作物耐低氮胁迫能力的关键。对于评价油菜耐低氮胁迫能力的指标，目前较为常用的有产量、生物量和根系形态等 ，。或多或少存在着相关性的多个性状所提供的低氮耐受信息相互交叉重叠，展现出数量性状的特点。因此，利用常规鉴定方法通过单一指标来评价油菜的耐低氮胁迫能力，可能因基因型差异而造成植株对低氮胁迫的反应不尽相同，使得结果较为片面；且某个单一指标也可能与其他指标存在相互作用，使得结果准确性差，难以准确而全面地反映不同基因型耐低氮胁迫能力的强弱。需要通过测定多种指标去综合评价油菜的耐低氮胁迫能力，使得试验结果更加真实可靠。利用多元统计分析方法，即主成分分析法和隶属函数法，能够将不同的性状指标，转化成新的综合指标，其数目更少、相对独立性更强，可以较为全面地表现作物不同性状指标的全部信息，利用统计学方法计算出材料间的低氮耐受综合指标值和 值，从而准确且全面地反映出作物的耐低氮胁迫能力。随着人们对农业生产实际中过量施用氮肥可能造成的环境问题的逐渐关注，如何在保证作物产量的前提下减少氮肥的使用越来越成为科研人员的研究热点。在低氮胁迫条件下选育氮高效即耐低氮胁迫能力强的品种（系），是解决这一问题的关键因素之一。目前，对水稻、小麦、玉米、大豆、高粱等作物在低氮胁迫条件下的氮高效遗传机制及种质筛选方法已经有了大量的研究 。甘蓝型油菜不同基因型间耐低氮胁迫能力差异显著，具有广泛的遗传变异 。因此，选育耐低氮胁迫能力强的油菜种质以减少氮肥施用具有广阔的研究前景。邹小云等利用 差异较大的甘蓝型油菜种质湘油 和 为亲本，以杂交和回交方法构建了 个世代，通过主基因配合多基因混合遗传模型，发现不同的氮素营养效率鉴定指标对应着不同的遗传江苏农业科学 年第 卷第 期表 油菜各基因型耐低氮胁迫能力的综合指标值（）、权重（）、（）值、值材料编号 （）（）（）（）权重（）江苏农业科学 年第 卷第 期表 低氮胁迫耐受型甘蓝型油菜种质材料编号材料名称 模型，因而对种质氮素营养效率的改良也应该按照领域和目的的不同，选择最佳的氮素营养效率类别；他们也证明了通过不同氮效率油菜种质间的杂交和回交可以选育出氮高效的基因型 。张浩等在高氮和低氮 种营养条件下，以不同基因型油菜种质为材料，通过 年 点的田间试验，发现油菜具有较大的减氮增产增效潜力，且不同的基因型间潜力差异较大，他们还以产量均值为依据，将供试材料划分为不同氮效率的 种类型 。利用主成分分析法和隶属函数法综合评价作物耐低氮胁迫能力的研究已经在谷子、玉米和小麦等作物中得到了初步应用 。秦璐等也利用该方法对甘蓝型油菜盛花期耐低氮胁迫能力进行了比较分析，但在苗期对甘蓝型油菜耐低氮胁迫能力鉴定的研究尚未见报道 。本研究以各性状的低氮胁迫耐受指数为标准评价油菜单一指标的耐低氮胁迫能力，利用主成分分析法将 个油菜农艺性状简化为了 个相互独立的综合指标。这种方法既能较为全面地反映原性状指标的信息，且更突显出了不同基因型油菜的性状差异。根据贡献率计算出各综合指标相应的隶属函数值，加权后即可得到不同基因型油菜耐低氮胁迫能力的综合评价值，即值。进一步通过聚类分析，将 份油菜种质划分为 个类型，其中 、等 个材料被划分为低氮胁迫耐受型种质，、等 个材料被划分为低氮胁迫敏感型种质，、等 个材料则被划分为低氮胁迫中间型种质。作物抗逆性种质筛选是一项工作量巨大、时间成本较高的工作，如何在种质资源鉴定过程中，筛选出操作简单、结果可靠的鉴定指标对于作物抗逆育种具有重要意义。目前已经有许多与作物耐低氮胁迫能力相关的农艺和生理指标被应用于研究，生物量、根系形态、株高、叶片数等指标均对低氮胁迫有显著响应，可作为鉴定油菜耐低氮胁迫能力的指标 ，。本研究根据主成分分析结果，推测 、和 这 个低氮胁迫耐受指数对油菜耐低氮胁迫能力可能具有较大的影响，因此以 个低氮胁迫耐受指数为自变量，以 值作为因变量进行了多元回归分析，发现这 个指标确实与 值之间存在显著线性关系，在一定程度上可以作为油菜苗期耐江苏农业科学 年第 卷第 期低氮胁迫能力的鉴定指标。结论本研究通过对 份油菜种质 个农艺性状的低氮胁迫耐受指数进行主成分分析，简化为了 个综合指标。利用隶属函数法计算出各材料的 值，并依据 值进行聚类，将 份油菜种质划分为个类型，低氮胁迫耐受型种质 个，低氮胁迫敏感型种质 个，低氮胁迫中间型种质 个。通过多元线性回归，建立了以 、和 计算不同基因型 值的回归方程，简化了油菜耐低氮胁迫能力鉴定的评价指标，使得鉴定过程更为迅速、简单，为进一步改良现有优质油菜种质的耐低氮胁迫能力提供了理论基础。本研究重点关注了低氮胁迫对油菜苗期农艺性状的影响，而尚未对生理性状进行研究。因此，后续有必要继续开展相关研究，筛选到可靠的耐低氮胁迫能力较强的油菜种质资源。参考文献：沈金雄，傅廷栋我国油菜生产、改良与食用油供给安全 中国农业科技导报，（）：刘成，冯中朝，肖唐华，等我国油菜产业发展现状、潜力及对策 中国油料作物学报，（）：任涛，鲁剑巍中国冬油菜氮素养分管理策略 中国农业科学，（）：，（）：，（），（）：张福锁，申建波，危常州，等绿色智能肥料：从原理创新到产业化实现 土壤学报，（）：陈淑婷，李华，陈飞平油菜花海在乡村旅游中的应用 农村经济与科技，（）：万何平，何冰冰，陈敬东，等一种高通量鉴定甘蓝型油菜耐低氮种质的方法及其应用 江苏农业科学，（）：，（）：邹小云，官梅，官春云甘蓝型油菜氮素高效吸收的植株形态与生理特性研究 作物杂志，（）：，：，（），：，：熊洁，邹小云，陈伦林，等油菜苗期耐铝基因型筛选和鉴定指标的研究 中国农业科学，（）：，（）：王建强，韩配配，李银水，等不同氮效率油菜苗期根系形态及养分累积差异研究 中国油料作物学报，（）：巢成生，王玉乾，沈欣杰，等甘蓝型油菜苗期氮高效吸收转运特征研究 中国农业科学，（）：秦俭，蒋开锋，罗婧，等不同氮肥水平下香型杂交水稻产量特点与品种筛选 西南农业学报，（）：张恒，陈艳琦，任杰莹，等西南麦区小麦苗期氮高效品种筛选及指标体系构建 四川农业大学学报，（）：，杨豫龙，赵霞，王帅丽，等黄淮海中南部玉米氮高效品种筛选及产量性状分析 玉米科学，（）：远月丽大豆苗期氮高效利用种质筛选及遗传分析 北京：中国农业科学院，张飞，王佳旭，张旷野，等低氮逆境下氮高效高粱群体微环境及光合荧光应答效应研究 山西农业大学学报（自然科学版），（）：杨睿，伍晓明，安蓉，等不同基因型油菜氮素利用效率的差异及其与农艺性状和氮营养性状的关系 植物营养与肥料学报，（）：邹小云，官春云 甘蓝型油菜氮素营养效率的遗传效应分析 分子植物育种，（）：张浩，李双，叶祥盛，等甘蓝型油菜减氮增效潜力评价及种质资源筛选 中国油料作物学报，（）：连盈，卢娟，胡成梅，等低氮胁迫对谷子苗期性状的影响和耐低氮品种的筛选 中国生态农业学报，（）：师赵康低氮胁迫下玉米幼苗性状鉴定及叶片中低氮响应基因的转录组分析 太谷：山西农业大学，胡成梅小麦苗期耐低氮胁迫相关性状的全基因组关联分析 太谷：山西农业大学，秦璐，韩配配，常海滨，等甘蓝型油菜耐低氮种质筛选及绿肥应用潜力评价 作物学报，（）：，（），：江苏农业科学 年第 卷第 期