不结球白菜(
Brassica campestris L.)属于十字花科芸薹属,是一年生或二年生草本植物,原产于中国,栽培历史悠久。在长江中下游地区以其鲜嫩、风味独特、营养丰富等特点深受人们喜爱。不结球白菜全国栽培面积约为140万hm
2(2100万亩),已成为我国第三大蔬菜作物,年产值1600亿元(人民币,下同)以上
[1],在蔬菜的生产和供应方面占有重要的市场地位
[2-3]。随着人们生活水平的不断提高,对不结球白菜的产量和品质也提出了更高的要求。因此,培育高产优质的不结球白菜品种已成为当前育种工作的核心目标。通过系统地收集和保存不结球白菜的种质资源,进行遗传多样性研究和鉴定,对促进育种进展和品种创新具有至关重要的意义。
目前表型性状是作物育种和生物学研究的起点和最终落脚点,开展表型性状调查并利用表型数据进行遗传多样性分析是大白菜种质资源高效利用的有效途径,该方法在玉米
[4]、萝卜
[5]、辣椒
[6]、甜瓜
[7]等多种作物上均有报道。国内外研究主要集中在品种的评价与鉴定方面,涵盖了遗传多样性、耐热性、表型特征以及分子标记等多个方面的分析。在不结球白菜遗传多样性分析中韩建明
[8]最早对125份不结球白菜进行遗传多样性分析,通过变异性分析、主成分分析、聚类分析发现普通白菜和其它种类之间存在亲缘关系。何晓丽等
[9]则通过多次田间调查和遗传多样性分析,对54份不结球白菜种质资源进行了全面评估,发现所收集的种质资源在质量性状和数量性状方面具有丰富的多样性和差异性。范伟强等
[10]对天津地区的33个大白菜品种进行的表型多样性分析,结果显示大白菜品种间的变异性在短缩茎高方面最为显著,而叶球横径差异较小。并且通过主成分分析和聚类分析将样本分为6个组群,为大白菜新品种的推广和更新提供了参考依据。胡海娇等
[11]通过对48份不结球白菜种质资源分别进行了田间表型鉴定和遗传多样性分析,研究发现多个数量性状存在显著的变异,不同种质资源可以通过系统聚类分为不同的类群,这些研究为不结球白菜的品种改良和资源利用提供了重要的数据支持。在耐热材料分析中余长春等
[12]通过遗传多样性分析了自交不亲和系材料,并通过高温环境测定热害试验,发现不结球白菜(编号C72)热害指数最低,表明该材料具有较强的高温耐受性。另外Miao等
[13]利用单核苷酸多态性分子标记技术,分析不同中国白菜品种的遗传多样性。Chand等
[14]研究印度芥菜的遗传成分和多样性,基于不同的形态生理性状进行分析,利用统计学方法和遗传标记技术,评估不同品种之间的遗传变异和多样性,探讨其与形态和生理性状的关系。
本试验收集了47个地区的共196份不结球白菜种质资源,对196份不结球白菜进行表型遗传多样性分析、相关性分析、主成分分析、聚类分析和熵权-TOPSIS综合评价,旨在为不结球白菜育种质资源优异表型性状的发掘及其在育种上的高效应用提供参考。
1 材料与方法
1.1 材料
试验材料为实验室前期从47个地区收集的196份种质,由山西农业大学园艺学院蔬菜所提供,信息见
附表1。主要类型有:菘菜、快菜、长梗白菜、青梗菜、奶白菜(
图1)。
1.2 试验方法
试验地点位于山西省晋中市榆次区东阳镇山西农业大学试验基地(37.82°N,112.87°E)。2023年9月1日,于温于温室大棚进行穴盘育苗,每个材料育苗60株。田间试验设计3次重复,每份材料种植15株,2023年9月20日定植到田间。采用完全随机区组设计,小区种植面积为3 m2。每份材料种植15株;种植模式为平畦覆膜栽培,膜上定植,每畦定植双行,株行距15 cm×15 cm,小区的四周设置保护行。在生长期统一进行栽培管理,防治病虫害。
2023年秋季进行田间鉴定,统一调查中心柱纵径、株型、束腰、株高、菜心充实度、叶片形状、叶片尖端形状、叶片叶缘波状、叶面泡状程度、叶片边缘状态、叶片颜色、叶面光泽度、叶脉鲜明度、叶数量、叶长、叶柄长、叶宽、叶柄宽、叶柄厚、叶柄颜色和叶柄横切面形状共21个表型性状(
表1和
表2)。根据GB/T19557.5-2017《植物品种特异性、一致性和稳定性测试指南大白菜》和《不结球白菜种质资源描述规范和数据标准》
[15]采集表型数据,每个材料随机调查5株,取平均值进行数据统计分析。
1.3 数据统计与分析
应用Excel 2016对源数据进行统计与整理。应用SPSS 26.0软件对21个性状进行香农多样性指数(h=-∑[(pi)×In(pi)])、变异系数计算,相关性分析、主成分分析和系统聚类分析;利用Origin 2021绘图软件绘制了相关性和系统聚类分析热图。最后利用熵权-TOPSIS综合评价模型对196份不结球白菜种质进行综合评分排序。
2 结果与分析
2.1 196份不结球白菜种质资源主要表型性状多样性分析
在196份不结球白菜种质资源中,21个表型性状呈现出明显的差异,变异系数介于13.96%至50.78%,表明变异性丰富(
表3)。中心柱纵径的变异系数最大,为50.78%,变异为2.00~30.00 cm;其余性状的变异系数依次降低,中柱纵径50.78%>菜心充实度46.20%>叶柄长44.83%>叶柄厚40.04%>叶面泡状程度36.13%>叶柄颜色34.26%>株高34.01%>束腰31.15%>叶宽30.86%>叶长30.24%>叶片形状29.53%>叶柄宽27.67%>叶片颜色26.73%>叶片叶缘波状21.99%>叶片尖端形状21.21%>叶面光泽度20.88%>叶数量19.99%>叶片边缘状态18.49%>叶柄横切面形状17.29%>叶脉鲜明度13.96%>株型11.96%。结果表明:21个表型性状变异系数均大于10%,意味着这196份不结球白菜在中柱纵径,叶片数,株高等方面差异性较大,具有丰富的遗传多样性,有利于种质资源的筛选和比较。
2.2 相关性分析
通过不结球白菜种质资源相关性分析(
图2)可知,在21个表型性状指标中,不结球白菜的中柱纵径与株型、叶数量、叶柄颜色、叶柄横切面形状呈正相关,与株高、叶面光泽度、叶长、叶柄长和叶宽呈负相关;束腰与菜心充实度、叶数量、叶柄颜色呈显著正相关,与株高、叶片性状、叶柄长和叶宽呈负相关;株高与叶长、叶宽、叶柄长和叶柄厚显著正相关,与叶片边缘状态、叶片颜色和叶柄颜色呈负相关;叶片形状与叶宽、叶面光泽、叶面颜色和叶片尖端形状呈正相关,与叶数量、叶长和叶柄长呈负相关;叶数量与株型、叶柄颜色、中柱纵径和菜心充实度呈正相关,与叶片形状和叶片尖端形状呈负相关;叶柄颜色与中柱,株型,束腰呈显著正相关,与叶面泡状程度、叶宽、叶片叶缘波状、叶片形状呈负相关;叶长与株高、叶柄长、叶宽呈正相关,与叶柄颜色、叶片颜色,束腰呈负相关。分析后发现表型性状指标间相关性显著,关联系数较强,指标信息间存在重叠,更有利于优良不结球白菜种质筛选。
2.3 196份不结球白菜种质资源主要表型性状指标的主成分分析
如
表4所示,前5个主成分的贡献率为:5.71%、8.48%、9.19%、14.54%和21.29%对应的特征值为1.26、1.87、2.02、3.20和4.68。这些值均大于1,说明它们涵盖了主要的表型性状信息。因此,选取前5个主成分进行分析是合理且可靠的。第1主成分的特征值为4.68,贡献率为21.29%,以株高、叶长、叶柄长和叶宽,为主要性状。第2主成分的特征值为3.20,贡献率为14.54,以菜心充实度、株形和束腰为主要性状,说明第2主成分反映的是不结球白菜的株型和商品性信息。第3主成分的特征值为2.02,贡献率为9.19%,以束腰、叶片尖端形状、叶柄宽和叶柄厚为主要性状。第4主成分的特征值为1.87,贡献率为8.48%,以中柱纵径、叶数量和叶柄横切面形状为主要性状。说明第4主成分反映的是植株的耐抽薹信息。第5主成分的特征值为1.26,贡献率为5.71%,以株型、叶脉鲜明度和叶面光泽度为主要性状。
2.4 不结球白菜主要性状聚类分析
通过观察表型,发现中柱纵径值越高说明材料越容易抽薹,属于不耐寒材料。中柱值越低说明材料不易抽薹,属于耐寒材料。另外我们发现叶片数量和植株高度和产量是相关的,叶片数量越大和株高越高,说明产量可能越大。
对196份不结球白菜种质资源的21个性状进行聚类分析,在欧氏距离0.98处将其分为5个组群(
图3),计算各类群性状的平均值(
表5)。第1组群包括13份不结球白菜种质,该类种质平均中柱纵径9.2 cm,为各类群最高,株高19.7 cm、叶数量16.84、叶长20.84 cm、叶柄长8.76 cm、叶宽11.21 cm、叶柄宽3.8 cm、叶柄厚0.79 cm。特征为具有较高的中柱纵径9.2 cm,耐抽薹性差,不适合耐抽薹材料的选育。第2组群包括24份不结球白菜种质,该类种质平均中柱纵径6.77 cm、株高14.08 cm、叶数量11.91、叶长15.5 cm、叶柄长7.58 cm、叶宽12.15 cm、叶柄宽3.69 cm、叶柄厚0.65 cm。中柱纵径为6.77 cm,特征为株高较小,菜心充实度低,适用于耐热矮棵青梗菜的选育。第3组群包括30份不结球白菜种质,该类种质平均中柱纵径5.6 cm,株高29.7 cm、叶数量14.9、叶长29.41 cm、叶柄长15.90 cm、叶宽17.48 cm、叶柄宽3.88 cm、叶柄厚0.77 cm。特征为中柱纵径最小,株型大,植株高,叶面积大,耐抽薹性较好,产量高,适合用于大棵耐抽薹材料的选育。第4组群包括37份不结球白菜种质,该类种质平均中柱纵径6.19 cm、株高20.85 cm、叶数量16、叶长20.79 cm、叶柄长9.26 cm、叶宽11.47 cm、叶柄宽4.83 cm、叶柄厚0.92 cm。特征为中柱纵径6.19 cm,植株高度为20.85 cm,束腰性强、菜心充实度高、商品性好,可用于耐热型青梗菜的选育。第5组群包括92份不结球白菜种质,该类种质平均中心柱纵径7.5 cm、株高16.5 cm、叶数量15.1、叶长16.96 cm、叶柄长6.8 cm、叶宽9.41 cm、叶柄宽3.3 cm、叶柄厚0.55 cm。特征为叶宽和叶柄宽为最小,其它性状处于中等。
2.5 196份不结球白菜基于熵权-TOPSIS法的不结球白菜综合评价
通过对不结球白菜21个表型性状指标的筛选,发现15个性状指标可以从多方面进行熵权-TOPSIS综合评价,这15个性状指标均与选育目标相关,能够较好反应该资源的综合性质,如
表6所示。不结球白菜各表型性状信息熵和权重成反比,在15个性状中,菜心充实度的信息熵最小(0.93),权重最大(0.15),说明该指标离散程度较大,对材料综合评价影响较大,叶柄长信息熵最大(0.99),权重最小(0.01),对材料综合评价影响较小。
通过对196不结球白菜份材料的中柱纵径、株型、束腰,株高等15个主要表型性状进行综合评估,筛选出了排名前20的材料数据见
表7和附
图1,其中代号189、184、66为综合评价前三名(
表8)。
3 讨论
种质资源是自然界中非常宝贵的材料和基因库,在优质种质发掘、优异性状和基因利用方面发挥重要作用,是新品种选育的基础
[16]。
变异系数和遗传多样性指数是种质资源多样性分析的常用参数和指标。变异系数越高,说明性状的离散程度越大,个体之间的差异越明显
[17]。研究表明变异系数大于20%,可说明不同样本的多样性丰富
[18]。本研究分析发现,中柱纵径、菜心充实度和叶柄长度变异较大,变异数较高。根据中柱总径的数值变异分析,表明在耐抽薹方面不同材料差距较大。菜心充实度和叶柄长度也有较大的变异系数,表明在产量和商品性方面也具有较大的变化。结合中柱纵径等变异数据,推测可能是植物抽薹以后,叶柄长度也有变化。中柱纵径变异性分析与肖艳等
[19]、韩睿等
[20]研究结果一致。本研究中遗传指数最大的是株型,与韩睿等
[19]研究结果不一致,原因可能为材料和气候不同导致植物表型差异。
相关性分析和主成分分析是表型性状分析中常用的2种方法,用于评估2个或多个变量之间的关系强度和方向,能够将多个复杂的原始指标转化为几个互不相关且具有代表性的综合指标,通过提取主要因子,从而更加客观地筛选出优良种质
[21]。本研究从21个表型性状中选取5个主成分,贡献率最高为21.29%,累计贡献率为59.20%,株高、叶长、叶柄长、束腰、中柱纵径、株型、叶数量等指标占据主要优势。这5个主成分基本能够说明原始指标的大部分信息,可作为不结球白菜的表型综合指标。研究发现不结球白菜的中柱纵径与株型、叶数量、叶柄颜色、横切面形状和菜心充实度呈极显著正相关,与株高、叶面光泽度、叶长、叶柄长和叶宽呈负相关。这与余阳俊等
[22]研究结果一致。
熵权⁃TOPSIS综合评价和聚类分析是一种客观赋权方法,常用于多指标评价问题中,以确定各指标的权重。该方法基于信息熵的原理,利用数据的离散程度来确定指标的权重
[23-24]。TOPSIS法的核心思想是通过比较各个备选方案与理想解的接近程度来进行排序,其中理想解是所有指标中的最佳值集合,负理想解是所有指标中的最差值集合
[25-26]。余忠浩等
[27]采用TOPSIS法对109个大豆品种进行了分析,从中筛选出7个表现优异的品种,为大豆品种选育提供了科学依据。王文迪等
[28]采用TOPSIS法对9份大豆进行了分析,对大豆进行了综合评价和排序,发现了通辽地区大豆选育的优良亲本。本研究采用熵权-TOPSIS综合评价法,对196份不结球白菜进行综合排序,筛选出20个优异材料。结合聚类分析在欧氏距离0.98时将196份材料分为了5个组群,第1组群叶片数最多,中柱纵径最大,说明此类材料不耐抽薹,不耐寒。第2组群类株高和叶柄较小,适用于耐热材料型矮棵青梗菜的选育。第3组群特征为中柱纵径短,但株型大,叶片大,耐抽薹性较好,产量高,适合用于大棵耐抽薹材料的选育。第4组群特征为中柱纵径适中,此类材料抽薹后,植株高度为20.85 cm,束腰性强,商品性好,可用于耐热型青梗菜的选育。第5组群叶宽和叶柄宽为最小,其他性状处于中等。
不结球白菜的表型鉴定在农业育种、品质提升、抗性研究等多个领域具有重要的指导意义,为改良不结球白菜的栽培和生产提供了科学依据。目前,本研究只针对表型性状进行分析,下一步可结合分子标记和基因编辑技术进行更深层次的分析。
4 结论
本研究通过遗传多样性分析、相关性分析、主成分分析、聚类分析和利用15个农艺性状进行了熵权-TOPSIS综合评价,对196份不结球白菜的21个表型性状进行多角度分析。由结果可知,所选不结球白菜材料遗传差异大,变异丰富。通过熵权-TOPSIS筛选出20个优质材料。
附加材料:附表1(试验材料编号及来源)、附
图1(排序前20材料图片)见中国知网增强出版内容。
山西农业大学科技创新提升工程项目(CXGC2023031)