Genetic Diversity on Selected Agronomic Traits of Arabica Coffee Germplasm
-
摘要:目的 分析小粒种咖啡种质资源农艺性状和品质性状的多样性,为小粒种咖啡品种选育提供科学依据。方法 以60份小粒种咖啡种质资源为材料,进行8个农艺性状和5个品质性状测定,并采用相关性分析、主成分分析和聚类分析等方法进行多样性分析。结果 在8个农艺性状中,变异系数以鲜干比最高(25.15%),其次为鲜果重(19.85%),以出米率最低(1.97%);遗传多样性指数以种子长为最高(2.04)。在5个品质性状中,变异系数以绿原酸最高(22.68%),其次为脂肪(12.40%),蛋白质为最低(5.20%);遗传多样性指数以绿原酸最高(2.07),其次为脂肪(2.06)和蛋白质(2.05)。性状相关性分析结果表明,干豆重与出米率呈显著正相关,种子大小与百粒重、出米率、蛋白质呈极显著正相关,百粒重与出米率、蛋白质、脂肪呈极显著正相关,出米率与蛋白质、脂肪呈极显著正相关,蛋白质与脂肪呈极显著正相关,咖啡因与蔗糖呈极显著负相关。主成分分析结果表明,前4个主成分因子(PC1、PC2、PC3、PC4)包含了12个农艺性状,累计贡献率达69.17%,表明这些性状是造成小粒种咖啡种质多样性的主要因素。聚类分析结果表明,60份小粒种咖啡种质可分为两个大类群,第I类群有53份资源,第II类群有7份资源。第II类群的鲜果重大,尤其是II-2类群的鲜果重和干豆重。结论 60份咖啡种质资源存在丰富的遗传多样性,可通过引种驯化、杂交、嫁接和分子育种等技术将其优异基因延续下去。第II-2 类种质的鲜果重和干豆重大,可作为选育和改良品种的基本材料。Abstract:Objective Diversity on the agronomic traits relating to breeding purpose of arabica coffee germplasms in collection was studied.Method From 60 arabica coffee germplasms, 8 agronomic and 5 quality traits were targeted and subjected them to the correlation, principal component, and cluster analyses for the genetic diversity study.Result Among the 8 agronomic traits, the fresh/dry ratio had the highest coefficient of variation of 25.15% followed by fresh fruit weight of 19.85%, while green bean rate of 1.97% being the lowest. Seed length ranked the highest on the genetic diversity indices. On the 5 quality traits, the chlorogenic acid content showed the highest coefficient of variation of 22.68% followed by fat of 12.40%, and protein of 5.20% being the lowest of all. Chlorogenic acid also had the greatest genetic diversity index of 2.07 followed by fat of 2.06 and protein of 2.05. The dry bean weight of the germplasms significantly correlated with their green bean rate; the seed size with 100-seed weight, green bean rate and protein; the 100-seed weight with green bean rate, protein and fat; and, the green bean rate with protein and fat. The contents of protein and fat showed an extremely significant correlation; whereas, those of caffeine and sucrose an extremely significant inverse correlation. The principal component analysis indicated the top 4 main components (i.e., PC1-PC4) included the 12 agronomic traits and had a combined contribution of 69.17% of the total constituting the major factors that affected the diversity. A cluster analysis classified the 60 germplasms into two groups, 53 in Group I and 7 in Group II. Group II was high on the fresh fruit weight, especially, Subgroup II-2, which exhibited high fresh and dry bean weights.Conclusion The arabica coffee germplasms under study were rich in diversity. The desirable genes could be preserved by domestication, hybridization, grafting and/or molecular breeding. Subgroup II-2, being high on fresh and dry bean weights, could conceivably be used as a parent for new variety breeding.
-
0. 引言
【研究意义】玉米(Zea mays L.)是我国的主要粮食作物之一[1],粮饲兼用,也是重要的工业原料,在国家粮食安全保障体系中占有举足轻重的地位[2]。随着国家对知识产权保护范围的拓宽和对科技人才的不断重视,育种人对品种权的保护意识也在不断增强[3]。再者,玉米种质资源的有效挖掘利用、分子辅助选育、基因工程等生物新技术在玉米育种过程中的应用,使得玉米育种年限不断缩短,近年申请品种审定和品种权保护的玉米新品种数量也在急剧上升[4]。植物新品种审定、登记和保护制度是农业知识产权中最重要的一种制度,在鼓励育种创新、丰富我国种质资源、保障粮食安全生产等方面有重要的推动作用[5−6]。然而,无论是申请品种审定还是品种权保护,在其实质性审查时都要进行特异性(distinctness)、一致性(uniformity)和稳定性(stability)测试(简称DUS 测试)[7−8]。当前,玉米DUS测试是参照《植物品种特异性、一致性和稳定性测试指南 玉米》(GB/T 19557.24—2018)[7](以下简称《指南》)的要求开展。《指南》中要求测试的数量性状,一般受环境条件影响较大,在田间种植试验中通常需要设置一些标准品种,以矫正环境变化对数量性状表达的影响[8]。《指南》属于国家标准,在全国范围内流通使用。但经过标准品种的种植,发现《指南》中提供的标准品种在云南生态区域生长不太适宜,导致标准品种示例的性状表达状态在云南低纬高原生态区域存在不同程度变化,出现标准品种不标准等问题,限制了玉米DUS测试工作有序开展,影响测试结果可靠性和科学性[9−10]。【前人研究进展】李娟等[11]基于DUS测试,对玉米标准品种在贵州地区的表达差异及多样性进行分析,建议可将用于研究的36个标准品种调减至31个。刘迪发等[12]对玉米DUS测试标准品种在海南性状表达差异性进行分析,结果表明不同标准品种、不同性状在海南的表现与《指南》赋值的差异程度不同。徐振江等[13]对玉米DUS测试标准品种在华南亚热带生态区的表现与优化进行研究,研究结果认为存在一些标准品种的性状表达与《指南》赋值的符合系数不高。任丽等[14]用符合系数对标准品种在上海地区的性状表现与水稻测试指南中的代码进行差异性分析,并进一步对部分性状的标准品种进行了优化。徐振江等[15]对水稻新品种DUS测试目测性状标准品种进行筛选优化,筛选出16个标准品种就已涵盖原有的23个标准品种代表的性状,大大减少了标准品种数量。赵连佳等[16]对红花(Carthamus tinctorius L.)测试指南中的32份标准品种42个相关性状进行调查分析,结果表明大部分红花标准品种及性状在新疆的表现与测试指南相符。Pérez等[17]在不同环境下对13种板栗(Castanea mollissima BL.)性状进行测试筛,选出较稳定性显著性强的5个性状以确定DUS试验中使用的不同描述符的性质。数量性状的表达存在不同程度的差异,LIU等[18]以亚洲毛茛(Ranunculus japonicus Thunb.)DUS测试指南的制定为例,应用统计方法研究DUS测试数量性状筛选、分级以及标准品种选择。【本研究切入点】已报道的相关作物研究大多集中在对标准品种在不同区域性状表达的稳定性和多样性进行研究分析[19],对标准品种进行优化或筛选出更合适标准品种的报道很少,特别是对《指南》中玉米标准品种筛选优化以及有关低纬高原生态区域玉米DUS测试标准品种筛选的研究几乎没有报道。因此,筛选适合在低纬高原生态区域种植的玉米DUS测试标准品种迫在眉睫。【拟解决的关键问题】本研究以80个杂交玉米品种和20个玉米自交系为研究材料,参照《指南》[7]中的要求和《植物品种特异性 一致性 稳定性测试总论》[8](以下简称《总论》)中标准品种筛选原则,筛选出一套适合于云南等低纬高原区域种植的玉米DUS测试标准品种,并对筛选的标准品种进行聚类分析,为筛选适合在低纬高原区域种植的玉米标准品种等相关研究奠定理论基础。
1. 材料与方法
1.1 试验材料
本试验种植了100个适合在低纬高原生态区域生长的玉米品种,包括80个杂交玉米品种和20个玉米自交系。考虑到品种的代表性和尽可能选择遗传背景差异大的品种,从30家不同育种公司和单位收集品种,每家收集了2~5个品种。试验材料具体信息如表1所示。
表 1 100个试验材料品种名称及品种类型Table 1. Names and types of 100 maize varieties tested品种编号
Varieties
number品种名称
Varieties
name品种类型
Varieties
types品种编号
Varieties
number品种名称
Varieties
name品种类型
Varieties
types品种编号
Varieties
number品种名称
Varieties
name品种类型
Varieties
types品种编号
Varieties
number品种名称
Varieties
name品种类型
Varieties
types1 AA323 CS 26 恩白玉9号 CH 51 洱之玉918 CH 76 川单99 CH 2 AA128 CS 27 龙白玉2号 CH 52 富滇8号 CH 77 荣玉99 CH 3 AA235 CS 28 龙白玉999 CH 53 澄白2号 CH 78 弥玉369 CH 4 A001 CS 29 隆白玉3号 CH 54 秋庆402 CH 79 万丰玉 3612 CH 5 S01 CS 30 龙白玉3号 CH 55 砚白1号 CH 80 奎单828 CH 6 S02 CS 31 贵白玉2号 CH 56 凉都玉1号 CH 81 锋玉217 CH 7 S03 CS 32 贵白玉1号 CH 57 林新13号 CH 82 万丰玉939 CH 8 DT006 CS 33 弥玉129 CH 58 宣农玉7号 CH 83 云谷99 CH 9 DT12832 CS 34 腾丰玉198 CH 59 靖黄7号 CH 84 瑞喜731 CH 10 DT33631 CS 35 澄白2号 CH 60 云瑞 8311 CH 85 顶立313 CH 11 DT3363 CS 36 兆白99 CH 61 弥玉002 CH 86 滇喜518 CH 12 BC57 CS 37 隆白玉109 CH 62 云瑞909 CH 87 金黄10号 CH 13 TRWX-9 WS 38 古白玉1号 CH 63 云瑞208 CH 88 瑞黑糯1号 WH 14 DT5133 CS 39 靖青贮1号 CH 64 田单60 CH 89 云黄13号 CH 15 DT37465 CS 40 靖白1号 CH 65 洱之玉666 CH 90 亲明玉1号 CH 16 DT1461 CS 41 瑞佳甜4号 SH 66 珍谷 68 CH 91 富好2号 CH 17 DT37466 CS 42 云甜玉12号 SH 67 云瑞202 CH 92 水金玉211 CH 18 DT1581 CS 43 花超99 SH 68 金滇 1588 CH 93 足玉817 CH 19 DT6413212 CS 44 双色年华 SH 69 宣晟5号 CH 94 胜玉209 CH 20 DT43774 CS 45 毫发520 SH 70 禾圣707 CH 95 迈单520 CH 21 宣成白1号 CH 46 云甜玉11号 SH 71 禾圣668 CH 96 红玉618 CH 22 神禾788 CH 47 瑞佳甜3号 SH 72 易金273 CH 97 红玉711 CH 23 澄白3号 CH 48 玉研212 CH 73 云良10号 CH 98 高农玉203 CH 24 宣白玉10号 CH 49 子玉68 CH 74 渝单958 CH 99 雅玉2号 CH 25 澄白1号 CH 50 科元玉518 CH 75 西南玉99 CH 100 瑞甜7号 SH 表中字母“CS、WS、CH、SH、WH”分别表示“普通玉米自交系、糯质玉米自交系、普通玉米杂交种、甜玉米杂交种、糯质玉米杂交种”。
CS: common maize inbred lines; WS: waxy maize inbred line; CH: common maize hybrid; SH: sweet maize hybrid; and WH: waxy maize hybrid.1.2 试验方法
1.2.1 田间形态性状测试
田间试验设计与管理:田间测试于2022年4月至2023年11月在农业农村部植物新品种DUS测试(昆明)分中心嵩明基地进行两个周期的形态性状测试。地块平整,肥力中等偏上。试验地地理坐标为25.18° N,103.25° E,海拔为
1900 m,年均气温15 ℃左右,年降雨量在1200 ~1400 mm,年平均无霜期约235 d。试验采用随机区组设计,2次重复,塘播,行距75 cm,塘距35 cm,定苗时每穴留2株,4行区。四周设保护行,常规管理同大田一致。本试验材料于2022年4月29日播种,5月22日进行间苗、定苗。整个生育期除草2次,防病虫害喷洒药3次,中耕培土、起垄各1次。1.2.2 性状观测与判别
性状调查方法、时期、部位及判别严格按照《指南》[7]中的具体要求进行,共测试了100个试验材料的37个数量性状(表2)。
表 2 37个数量性状编号及名称Table 2. Codes and names of 37 quantitative characteristics性状编号
Characters
number性状名称
Characters name性状编号
Character
number性状名称
Characters name性状编号
Characters
number性状名称
Characters name2 幼苗:第一叶鞘花青甙显色强度 16 雄穗:侧枝弯曲程度 29 果穗:穗柄长度 3 叶片:绿色程度 17 雄穗:最低位侧枝以上主轴长度 30 果穗:长度 4 散粉期 18 *雄穗:最高位侧枝以上主轴长度 31 果穗:直径 5 抽丝期 19 雄穗:一级侧枝数目 32 果穗:穗行数 6 植株:上部叶片与茎秆夹角 20 雄穗:侧枝长度 33 果穗:形状 7 叶片:弯曲程度 21 茎秆:“之”字型程度 35 仅适用于甜玉米:籽粒黄色程度 9 雌穗:花丝花青甙显色强度 22 茎秆:支持根花青甙显色强度 36 仅适用于甜玉米:籽粒长度 10 雄穗:花药花青甙显色强度 23 叶片:长度 37 仅适用于甜玉米:籽粒宽度 11 雄穗:小穗密度 24 叶片:宽度 38 仅适用于甜玉米:籽粒皱缩程度 12 雄穗:颖片基部花青甙显色强度 25 叶:叶鞘花青甙显色强度 39 籽粒:类型 13 雄穗:颖尖花青甙显色强度 26 植株:穗位高度 43 穗轴:颖片花青甙显色强度 14 雄穗:颖片中部花青甙显色强度 27 植株:高度 15 雄穗:侧枝与主轴夹角 28 植株:穗位高与株高比率 1.2.3 玉米标准品种筛选
玉米标准品种筛选遵循《总论》中标准品种选择的原则进行[8],主要要求标准品种样品要易获得,性状表达的波动性不大,性状表达的变异范围尽可能覆盖,标准品种的数量尽可能少,必要时需要征求相关专家的意见。获得100个试验材料的37个形态性状数据后,进行标准品种筛选时对《指南》[7]中有9个表达状态的性状,遵循筛选的标准品种尽可能优先覆盖代码为1、3、5、7、9对应的状态;对有5个表达状态的性状,标准品种尽可能优先覆盖代码为2、4对应的状态;对有4个表达状态的性状,标准品种尽可能优先覆盖代码为2、3对应的状态。在这样的原则下,筛选的标准品种在满足能充分发挥示例作用下,可以最大限度减少所需要的标准品种数量,从而减少测试成本。在标准品种确定时,通常将符合某个性状某种表达状态的少数品种优先作为候选标准品种。再根据需要的目标性状表达状态逐一增加候选品种数量,最终再优化品种数量,删除冗余品种。
1.2.4 玉米标准品种聚类分析
综合两年的测试数据,以筛选的标准品种形态性状数据为基础,目测性状直接使用观测值,对于测量数据需先转化成代码值再进行分析。利用NTSYSpc2.10e统计软件中similarity程序计算相似系数,以Clustering程序中SHAN进行非加权组平均法(unweighted pair group method with arithmetic mean, UPGMA)聚类分析[20],绘制玉米品种的遗传聚类系统树。
2. 结果与分析
2.1 37个性状描述结果
由表3、4中性状描述的最大值、最小值及覆盖范围可知,研究材料在37个性状的绝大多数表达状态均有分布,但在极少数性状的个别表达状态上没有分布。结果表明,研究材料的来源还不够丰富,一定程度上也说明个别性状的表达状态设置存在不合理性,实际育种和生产中并不存在那么多的表达状态。由表3、4可知,37个性状的变异系数在5.70%~110.89%。其中,变异系数在15%~20%的性状有5个,小于15%的性状有10个,其他性状的变异系数均大于20%,说明大多数性状上获得的数据离散程度较高,100个研究材料仍然具有代表性。
表 3 基于研究材料的目测性状描述结果Table 3. Descriptions on appearances of maize under study性状编号
Characters
number最大值
Max最小值
Min覆盖范围
Coverage
area变异系数
Coefficient of variation/
%2 8 1 1~8 27.69 3 3 2 2~3 20.41 4 9 1 1~9 24.98 5 9 1 1~9 25.08 6 6 2 2~6 23.89 7 6 1 1~6 33.73 9 7 1 1~7 56.57 10 7 1 1~7 42.86 11 5 2 2~5 18.88 12 8 1 1~8 71.48 13 8 1 1~8 35.34 14 8 1 1~8 48.33 15 7 2 2~7 29.08 16 6 1 1~6 42.93 21 3 1 1~3 30.93 22 8 1 1~8 51.13 25 7 1 1~7 76.17 29 4 1 1~4 38.24 33 3 1 1~3 15.07 35 6 4 4~6 10.00 36 7 6 6~7 7.19 37 7 6 6~7 5.70 38 4 3 3~4 14.34 39 8 1 1~8 45.39 43 9 1 1~9 110.89 表 4 基于研究材料的测量性状描述结果Table 4. Descriptions on measured qualities of maize under study性状编号
Characters
number最大值
Max最小值
Min极差
Range覆盖范围
Coverage
area平均值
Average
value变异系数
Coefficient of
variation/%17 47.6 21.7 25.9 21.7~47.6 37.3 14.02 18 37.5 15.4 22.1 15.4~37.5 27.6 15.11 19 24.5 1.8 22.7 1.8~24.5 9.1 53.42 20 30.5 9.7 20.8 9.7~30.5 21.3 20.67 23 112.0 42.0 70.0 42.0~112.0 90.1 14.21 24 12.1 5.9 6.2 5.9~12.1 10.0 11.65 26 208.6 24.4 184.2 24.4~208.6 115.6 31.42 27 433.7 137.4 296.3 137.4~433.7 290.9 19.55 28 0.6 0.2 0.4 0.2~0.6 0.4 18.93 30 21.7 11.1 10.6 11.1~21.7 17.3 11.37 31 5.6 3.5 2.1 3.5~5.6 4.8 8.46 32 18.5 11.8 6.7 11.8~18.5 14.8 9.63 2.2 玉米标准品种筛选结果
通过对100个试验材料进行田间种植和数据采集,并遵循《总论》中标准品种选择的原则进行筛选。标准品种筛选结果如表5所示,共筛选到玉米标准品种29个,包括1个糯质玉米自交系、8个普通玉米自交系、1个糯质玉米杂交种、6个甜玉米杂交种和13个普通玉米杂交种。标准品种可覆盖37个数量性状,示例159种表达状态,平均每个品种可示例5.5种性状表达状态。玉米品种AA128、毫发520可分别示例19个和18个性状的表达状态,S02可示例12个性状的表达状态,S01可示例11个性状的表达状态,DT1581可示例10个性状的表达状态,A001可示例8个性状的表达状态,DT12832、龙白玉999可示例7个性状的表达状态,龙白玉3号、瑞黑糯1号可示例6个性状的表达状态,瑞佳甜3号、砚白1号、宣农玉7号可示例5个性状的表达状态,DT5133、隆白玉109、万丰玉939可示例4个性状的表达状态,S03、TRWX-9、云甜玉12号、富好2号、瑞甜7号可示例3个性状的表达状态,腾丰玉198、田单60、西南玉99可示例2个性状的表达状态,双色年华、云甜玉11号、瑞喜731、云黄13号、雅玉2号可示例1种性状的表达状态。平均每个性状约有4.3个性状表达状态可被标准品种示例,其中,最多的性状有5个性状表达状态可被标准品种示例,如散粉期;最少的性状只有1个性状表达状态可被标准品种示例,如叶片绿色程度。
表 5 玉米标准品种筛选结果Table 5. Selected standard varieties of maize品种编号Varieties number 品种名称
Varieties name性状编号及表达状态代码
Characters number and code of expression status2 AA128 2(7),3(2),4(7),5(7),7(1),11(5),14(5),16(1),17(5),18(5),20(3),22(1),23(3),25(1),27(6),29(4),31(5),32(5),43(1) 4 A001 10(3),11(3),12(6),13(3),22(7),23(4),25(7),29(2) 5 S01 2(5),9(1),10(1),14(1),15(3),21(3),22(3),24(3),26(1),27(3),28(1) 6 S02 4(5),5(5),9(3),13(7),19(1),20(1),21(2),22(5),24(1),28(3),30(5),39(2) 7 S03 26(5),27(9),30(7) 9 DT12832 6(2),20(5),23(5),24(9),25(5),28(5),32(6) 13 TRWX-9 23(2),26(2),30(3) 14 DT5133 26(5),27(7),31(7),33(1) 18 DT1581 4(9),5(9),7(5),9(7),16(3),24(5),26(7),28(7),31(3),32(4) 28 龙白玉999 6(4),12(8),14(3),17(7),19(7),27(9),30(5) 30 龙白玉3号 19(5),24(9),26(9),31(7),32(7),33(2) 34 腾丰玉198 9(5),24(3) 37 隆白玉109 19(9),23(7),24(5),33(3) 42 云甜玉12号 23(3),28(1),36(7) 44 双色年华 6(6) 45 毫发520 4(3),5(3),13(1),17(3),18(3),19(3),23(1),24(1),26(1),27(1),28(3),30(3),31(5),32(5),36(6),37(6),38(3),39(6) 46 云甜玉11号 35(6) 47 瑞佳甜3号 26(3),27(3),35(4),37(7),38(4) 55 砚白1号 7(3),10(5),16(5),23(5),28(9) 58 宣农玉7号 13(5),26(7),27(7),32(3),39(4) 64 田单60 12(2),30(7) 75 西南玉99 26(5),43(5) 82 万丰玉939 15(7),18(7),25(3),28(5) 84 瑞喜731 43(7) 88 瑞黑糯1号 14(7),15(5),27(5),31(3),39(8),43(9) 89 云黄13号 12(4) 91 富好2号 10(7),20(7),28(7) 99 雅玉2号 43(3) 100 瑞甜7号 2(1),4(1),5(1) 第三列中括号前面的数字表示该标准品种可示例的性状编号,括号里的数字为该标准品种可示例的表达状态代码。
Datum in front of parentheses on 3rd column is code of agronomic trait; that within parentheses, code of expressive status.2.3 标准品种的聚类分析
以筛选得到的29个玉米标准品种形态性状数据为基础,利用NTSYSpc 2.10e统计软件对筛选的玉米标准品种进行聚类分析,分析结果如图1所示。由图可知,本试验筛选的29个玉米标准品种间的相似系数为
0.7841 ~0.9098 ,品种编号分别为28和34、44和45两对品种间的相似系数最高,均为0.9098 ,其余大多数品种间的相似系数在0.8051 ~0.8889 。以0.7841 为阈值时,可以将29个标准品种分为两个类群,第Ⅰ类群包含6个品种,约占21%,主要为甜质、糯质型杂交玉米品种;第Ⅱ类群包含23个品种,约占79%,主要包括1个糯质玉米自交系、8个普通玉米自交系和14个普通玉米杂交种。3. 讨论
3.1 玉米标准品种的筛选
《指南》[7]中提供的标准品种在某些性状覆盖的表达状态较多,在某些性状可示例的标准品种较少,示例分布不均衡,不能起到很好的示例矫正作用。品种的一致性和稳定性较好,就能很好地发挥参考作用[8]。当前《指南》中提供的标准品种主要以自交系为主,杂交中占很小的比例,而且在低纬高原生态区域生长适应性不好,所以性状表达差异比较明显,标准品种的作用未能得到很好的发挥。由于玉米自交系区域适应没有杂交种强,其性状表达较杂交种更容易受环境因素影响,所以本研究设计的材料中杂交种材料较自交系要多。最终筛选的29个标准品种中杂交种有20个品种,约占总数的70%,自交系有9个品种。从筛选结果来看,本试验筛选的29个标准品种在有些性状的两端表达状态未能发挥示例作用,在一定程度也说明本研究涉及材料的多样性不够。
3.2 玉米标准品种的示例作用分析
与《指南》[7]中提供的玉米标准品种相比,《指南》中涉及59个玉米标准品种,覆盖37个数量性状,示例164个表达状态,平均每个标准品种示例2.8个性状表达状态。2012年版本[21]中的标准品种共计55个,覆盖37个数量性状,示例159个表达状态,平均每个标准品种示例4.3个性状表达状态。本试验最终筛选了29个玉米品种作为标准品种,覆盖37个数量性状,示例159个表达状态,平均每个品种可示例5.5个性状表达状态。由此可以看出,29个标准品种尽管示例的性状表达状态数量没有比之前两个版本《指南》中标准品种示例的多,但平均每个品种可示例的性状表达状态数量较多,对整个测试工作降低成本、减少工作量的效果明显。通过比较两个测试指南版本中[7,21]的玉米标准品种,最新指南版本有59个玉米标准品种,之前版本中有55个玉米标准品种,但两者之间的差异品种只有8个。说明最新《指南》中的玉米标准品种大多数仍用十多年前筛选的标准品种,但由于全球气候环境的变化及品种自身对生态条件适应能力的减弱,我们有必要加强玉米标准品种筛选方面的研究,及时更换标准品种,真正使标准品种发挥作用。
3.3 玉米标准品种聚类分析
根据聚类结果可以知道筛选的玉米标准品种间的相似系数为
0.7841 ~0.9098 ,较前人[11−13]对标准品种研究结果要高,虽然筛选的标准品种性状覆盖范围大、示例性强,但品种间的相似性也比较高。以0.7841 为阈值时,可以将29个标准品种分为两个类群,第Ⅰ类群包含6个品种,主要为甜质、糯质型杂交玉米品种;第Ⅱ类群包含23个品种,主要为普通型玉米品种。从聚类结果来看,本试验筛选的29个标准品种间相似度较高,说明此次筛选的标准品种间的遗传背景差异不大。从另外一方面也说明本试验中使用的37个数量性状的品种区分力一般。但本试验筛选的29个玉米标准品种可示例159个表达状态,覆盖的性状表达状态还是相对比较高,标准品种筛选效果比较理想。4. 结论
通过本试验研究,从100个玉米试验材料中筛选到29个玉米品种可作为标准品种,品种间相似系数为
0.7841 ~0.9098 ,覆盖37个数量性状,示例159个表达状态,平均每个品种可示例约5.5个性状表达状态,品种的示例能力较《指南》中标准品种还强。该29个玉米品种可作为云南等低纬高原区域进行DUS测试的标准品种,为今后进一步使用和标准品种相关研究奠定基础和提供参考。 -
表 1 供试咖啡种质资源
Table 1 Coffee germplasms used in this study
编号
No.种质代码
Germplasm code品种类型
Varietal type编号
No.种质代码
Germplasm code品种类型
Varietal type1 5 波邦 Bourbon 31 143 未知 Unknown 2 10 卡杜艾 Catuai 32 145 未知 Unknown 3 35 卡蒂姆 Catimor 33 147 卡蒂姆 Catimor 4 36 卡蒂姆 Catimor 34 148 卡蒂姆 Catimor 5 37 卡蒂姆 Catimor 35 149 卡蒂姆 Catimor 6 40 未知 Unknown 36 150 卡蒂姆 Catimor 7 41 波邦 Bourbon 37 151 未知 Unknown 8 42 卡蒂姆 Catimor 38 154 未知 Unknown 9 45 卡杜拉 Caturra 39 159 未知 Unknown 10 49 卡蒂姆 Catimor 40 162 卡蒂姆 Catimor 11 50 卡蒂姆 Catimor 41 163 卡蒂姆 Catimor 12 51 波邦 Bourbon 42 166 卡蒂姆 Catimor 13 55 未知 Unknown 43 169 卡蒂姆 Catimor 14 61 未知 Unknown 44 170 卡蒂姆 Catimor 15 78 卡蒂姆 Catimor 45 173 卡蒂姆 Catimor 16 86 波邦 Bourbon 46 174 卡蒂姆 Catimor 17 87 铁皮卡 Typica 47 175 卡蒂姆 Catimor 18 88 波邦 Bourbon 48 178 卡蒂姆 Catimor 19 91 卡杜艾 Catuai 49 180 卡蒂姆 Catimor 20 92 卡杜艾 Catuai 50 182 未知 Unknown 21 93 未知 Unknown 51 204 未知 Unknown 22 118 卡杜拉 Caturra 52 282 未知 Unknown 23 119 未知 Unknown 53 284 未知 Unknown 24 120 未知 Unknown 54 296 未知 Unknown 25 122 卡蒂姆 Catimor 55 390 萨其姆 Sachimor 26 129 未知 Unknown 56 391 萨其姆 Sachimor 27 132 未知 Unknown 57 393 萨其姆 Sachimorr 28 138 未知 Unknown 58 394 萨其姆 Sachimor 29 139 未知 Unknown 59 395 萨其姆 Sachimor 30 140 卡杜拉 Caturra 60 435 未知 Unknown 表 2 农艺性状变异系数及遗传多样性分析结果
Table 2 Coefficient of variation and genetic diversity on agronomic traits of coffee germplasms
性状
Characteristics平均值
Average最小值
Min.最大值
Max.标准差
S变异系数
CV/%遗传多样性指数
H′鲜果重 Fresh weight/kg 1.87 0.79 2.97 0.37 19.85 1.88 干豆重 Dry weight/g 390.58 300.49 592.68 71.19 18.23 1.91 鲜干比 Fresh/dry ratio/% 21.50 16.55 48.79 0.11 25.15 1.43 种子长 Seed length/mm 11.74 10.74 12.69 0.44 3.71 2.04 种子宽 Seed width/mm 7.98 7.36 9.05 0.32 4.07 1.92 种子厚 Seed thickness/mm 4.88 4.40 5.70 0.26 5.38 1.99 百粒重100-seed weight/g 19.51 15.06 27.13 19.51 10.81 1.88 出米率 Green bean rate/% 84.01 78.63 86.80 1.65 1.97 1.92 表 3 品质性状变异系数及遗传多样性分析结果
Table 3 Coefficient of variation and genetic diversity on quality traits of coffee germplasms
性状
Characteristics平均值
Average最小值
Min.最大值
Max.标准差
S变异系数
CV/%遗传多样性指数
H′蛋白质 Protein/% 13.40 11.98 15.43 3.45 5.20 2.05 脂肪 Fat/% 13.02 9.03 15.91 6.88 12.40 2.06 绿原酸 Chlorogenic acid/(mg·g−1) 16.59 8.10 24.56 16.46 22.68 2.07 咖啡因 Caffeine/(mg·g−1) 9.89 7.62 11.80 4.18 9.84 1.99 蔗糖 Sucrose/(g·hg−1) 8.48 7.10 10.00 2.9 9.42 0.55 表 4 农艺性状、品质性状相关性分析
Table 4 Correlation between agronomic and quality traits of coffee germplasms
变量
Variable鲜果重
Fresh weight干豆重
Dry weight鲜干比
Fresh/dry ratio种子长
Seed length种子宽
Seed width种子厚
Seed thickness百粒重
100-seed weight出米率
Green bean rate蛋白质
Protein脂肪
Fat绿原酸
Chlorogenic acid咖啡因
Caffeine蔗糖
Sucrose鲜果重
Fresh weight1 干豆重
Dry weight0.057 1 鲜干比
Fresh/dry ratio0.053 0.714** 1 种子长
Seed length0.012 0.755** 0.965** 1 种子宽
Seed width0.217 0.125 −0.025 −0.104 1 种子厚
Seed thickness0.045 0.186 −0.057 −0.033 0.061 1 百粒重
100-seed weight0.205 0.164 0.363 0.368** 0.047 0.418** 1 出米率
Green bean rate0.021 0.308* 0.491 0.435** 0.185 0.374** 0.742** 1 蛋白质
Protein0.07 0.382** 0.526 0.477** 0.168 0.468** 0.757** 0.803** 1 脂肪
Fat0.012 0.141 0.183 0.177 −0.026 0.276* 0.467** 0.477** 0.499** 1 绿原酸
Chlorogenic acid−0.706** −0.009 −0.093 −0.115 0.398 0.114 −0.087 −0.01 0.072 0.097 1 咖啡因
Caffeine0.017 −0.051 −0.18 −0.187 0.115 −0.08 −0.109 −0.305 −0.133 −0.188 0.028 1 蔗糖
Sucrose−0.048 0.167 0.178 0.271* −0.316 0.099 0.003 0.063 0.112 0.066 −0.175 −0.426** 1 注:* 表示在 0.05 水平上显著相关;**表示在 0.01 水平上显著相关。
Note:* and ** represented significant correlation (P<0.05) and extremely significant correlation (P<0.01) respecively.表 5 农艺性状、品质性状主成分分析
Table 5 Principal components of agronomic and quality traits of coffee germplasms
变量
Variable主成分
Principal componentPC1 PC2 PC3 PC4 鲜果重 Fresh weight −0.359 0.801 −0.078 −0.313 干豆重 Dry weight −0.150 0.669 −0.471 0.365 鲜干比 Fresh / dry ratio 0.173 −0.300 −0.326 0.837 种子长 Seed length 0.653 0.090 0.205 −0.209 种子宽 Seed width 0.751 0.439 0.184 0.028 种子厚 Seed thickness 0.904 0.019 0.070 0.031 百粒重100-seed weight 0.919 0.164 0.146 0.007 出米率 Green bean rate 0.697 0.163 0.053 0.221 蛋白质 Protein −0.196 −0.089 0.768 0.204 脂肪 Fat −0.330 0.558 −0.060 0.012 绿原酸 Chlorogenic acid −0.106 0.527 0.280 0.502 咖啡因 Caffeine −0.398 0.230 0.648 0.103 蔗糖 Sucrose 0.476 0.174 −0.274 −0.287 特征值 Eigen value 3.86 2.11 1.61 1.41 贡献率 Contribution rate/% 29.71 16.25 12.35 10.86 累计贡献率
Cumulative contribution rate/%29.71 45.96 58.31 69.17 -
[1] FAO Statistics[EB/OL]. (2019-01-18)[2019-08-12]. http://faostat.fao.org.
[2] VAN DER VOSSEN H, BERTRAND B, CHARRIER A. Next generation variety development for sustainable production of Arabica coffee (Coffea arabica L.): a review [J]. Euphytica, 2015, 204(2): 243−256. DOI: 10.1007/s10681-015-1398-z
[3] 赵青云, 普浩杰, 王秋晶, 等.咖啡果皮不同堆沤处理养分含量及对咖啡植株生长的影响[J/OL]. 热带作物学报, 2019:1-8[2019-09-21]. http://kns.cnki.net/kcms/detail/46.1019.s.20190920.1323.006.html. ZHAO Q Y, PU H J, WANG Q J, et al.Nutrient contents of coffee peel with different composting treatments and its effects on coffee plant growth[J/OL]. Chinese Journal of Tropical Crops, 2019:1-8[2019-09-21].http://kns.cnki.net/kcms/detail/46.1019.s.20190920.1323.006.html. (in Chinese)
[4] 周华, 张洪波, 夏红云, 等. 咖啡种质资源多样性研究 [J]. 中国热带农业, 2015(5):23−27. DOI: 10.3969/j.issn.1673-0658.2015.05.007 ZHOU H, ZHANG H B, XIA H Y, et al. Genetic diversity analysis of coffee germlasms [J]. China Tropical Agriculture, 2015(5): 23−27.(in Chinese) DOI: 10.3969/j.issn.1673-0658.2015.05.007
[5] 白学慧, 夏红云, 李锦红, 等. 咖啡种质资源抗锈性初步鉴定 [J]. 热带农业科学, 2014, 34(7):60−64. DOI: 10.3969/j.issn.1009-2196.2014.07.012 BAI X H, XIA H Y, LI J H, et al. Resistance identification of coffee germplasm resources to rust [J]. Chinese Journal of Tropical Agriculture, 2014, 34(7): 60−64.(in Chinese) DOI: 10.3969/j.issn.1009-2196.2014.07.012
[6] YAN L, OGUTU C, HUANG L F, et al. Genetic diversity and population structure of coffee germplasm collections in China revealed by ISSR markers [J]. Plant Molecular Biology Reporter, 2019, 37(3): 204−213. DOI: 10.1007/s11105-019-01148-3
[7] 闫林, 黄丽芳, 王晓阳, 等. 咖啡种质资源遗传多样性的ISSR分析 [J]. 南方农业学报, 2019, 50(3):491−499. YAN L, HUANG L F, WANG X Y, et al. Genetic diversity of coffee gemeplasms by ISSR analysis [J]. Journal of Southern Agriculture, 2019, 50(3): 491−499.(in Chinese)
[8] 闫林, 黄丽芳, 王晓阳, 等. 基于ISSR标记的咖啡资源遗传多样性分析 [J]. 热带作物学报, 2019, 40(2):300−307. YAN L, HUANG L F, WANG X Y, et al. Genetic diversity of coffee germplasms by ISSR markers [J]. Chinese Journal of Tropical Crops, 2019, 40(2): 300−307.(in Chinese)
[9] 黄丽芳, 董云萍, 王晓阳, 等. 利用RAPD标记分析咖啡种质资源的遗传多样性 [J]. 热带作物学报, 2014, 35(12):2313−2319. DOI: 10.3969/j.issn.1000-2561.2014.12.001 HUANG L F, DONG Y P, WANG X Y, et al. Genetic diversity analysis of coffee germplasms by RAPD markers [J]. Chinese Journal of Tropical Crops, 2014, 35(12): 2313−2319.(in Chinese) DOI: 10.3969/j.issn.1000-2561.2014.12.001
[10] 黄丽芳, 董云萍, 王晓阳, 等. 云南咖啡资源遗传多样性的RAPD分析 [J]. 中国热带农业, 2017(5):48−52. DOI: 10.3969/j.issn.1673-0658.2017.05.012 HUANG L F, DONG Y P, WANG X Y, et al. Genetic diversity analysis of coffee germlasms from Yunnan by RAPD Markers [J]. China Tropical Agriculture, 2017(5): 48−52.(in Chinese) DOI: 10.3969/j.issn.1673-0658.2017.05.012
[11] 陶爱芬, 祁建民, 林培青, 等. 红麻优异种质产量和品质性状主成分聚类分析与综合评价 [J]. 中国农业科学, 2008, 41(9):2859−2867. DOI: 10.3864/j.issn.0578-1752.2008.09.039 TAO A F, QI J M, LIN P Q, et al. Cluster analysis and evaluation of elite kanaf germplasm based on principal components [J]. Scientia Agricultura Sinica, 2008, 41(9): 2859−2867.(in Chinese) DOI: 10.3864/j.issn.0578-1752.2008.09.039
[12] 廖丽, 陈玉华, 赵亚荣, 等. 地毯草种质资源形态多样性 [J]. 草业科学, 2015, 32(2):248−257. DOI: 10.11829/j.issn.1001-0629.2014-0223 LIAO L, CHEN Y H, ZHAO Y R, et al. Morphology diversity of Axonopus compressus germplasm [J]. Pratacultural Science, 2015, 32(2): 248−257.(in Chinese) DOI: 10.11829/j.issn.1001-0629.2014-0223
[13] 魏忠芬, 李慧琳, 奉斌, 等. 贵州紫苏种质资源表型性状的遗传多样性 [J]. 西南农业学报, 2017, 30(1):45−52. WEI Z F, LI H L, FENG B, et al. Genetic Diversity of Phenotype Characters of Perilla frutescens Germplasm Resources in Guizhou [J]. Southwest China Journal of Agricultural Sciences, 2017, 30(1): 45−52.(in Chinese)
[14] 万述伟, 宋风景, 郝俊杰, 等. 271份豌豆种质资源农艺性状遗传多样性分析 [J]. 植物遗传资源学报, 2017, 18(1):10−18. WANG S W, SONG F J, HAO J J, et al. Genetic Diversity of Agronomic Traits in 271 Pea Germplasm Resources [J]. Journal of Plant Genetic Resources, 2017, 18(1): 10−18.(in Chinese)
[15] 代攀虹, 孙君灵, 何守朴, 等. 陆地棉核心种质表型性状遗传多样性分析及综合评价 [J]. 中国农业科学, 2016, 49(19):3694−3708. DOI: 10.3864/j.issn.0578-1752.2016.19.003 DAI P H, SUN J L, HE S P, et al. Comprehensive evaluation and genetic diversity analysis of phenotypic traits of core collection in upland cotton [J]. Scientia Agricultura Sinica, 2016, 49(19): 3694−3708.(in Chinese) DOI: 10.3864/j.issn.0578-1752.2016.19.003
[16] KISUA J, MWIKAMBA K, MAKOBE M, et al. Genetic diversity of sweet and grain sorghum populations using phenotypic markers [J]. International Journal of Biosciences, 2015, 6(9): 34−46. DOI: 10.12692/ijb/6.9.34-46
[17] 王黎明, 焦少杰, 姜艳喜, 等. 不同来源甜高粱种质资源的表型遗传多样性分析 [J]. 植物遗传资源学报, 2014, 15(2):411−416. WANG L M, JIAO S J, JIANG Y X, et al. Genetic diversity analysis on sweet Sorghum germplasm resources of different origins based on agronomical traits Genetic diversity analysis of sweet Sorghum germplasm resources from different origins using agronomical traits Genetic diversity analysis of sweet Sorghum germplasm resources from different origins using agronomical traits [J]. Journal of Plant Genetic Resources, 2014, 15(2): 411−416.(in Chinese)
[18] 黄家雄, 吕玉兰, 程金焕, 等. 不同海拔对小粒种咖啡品质影响的研究 [J]. 热带农业科学, 2012, 32(8):4−7. DOI: 10.3969/j.issn.1009-2196.2012.08.002 HUANG J X, LV Y L, CHENG J H, et al. Preliminary study on the influence of different altitudes on the quality of coffee Arabica [J]. Chinese Journal of Tropical Agriculture, 2012, 32(8): 4−7.(in Chinese) DOI: 10.3969/j.issn.1009-2196.2012.08.002
[19] 王光路, 李长文, 刘志达, 等. HPLC法同时测定绿咖啡豆中绿原酸及咖啡碱含量 [J]. 食品研究与开发, 2017, 38(23):155−159. DOI: 10.3969/j.issn.1005-6521.2017.23.030 WANG G L, LI C W, LIU Z D, et al. Simultaneous quantitative analysis of chlorogenic acid and caffeine in raw coffee beans by HPLC [J]. Food Research and Development, 2017, 38(23): 155−159.(in Chinese) DOI: 10.3969/j.issn.1005-6521.2017.23.030
[20] 高贤玉, 张发明, 柏天琦, 等. 莲雾果实糖酸含量分析 [J]. 热带农业科学, 2019, 39(3):75−79, 90. GAO X Y, ZHANG F M, BAI T Q, et al. Sugar and Organic Acid Contents in the Fruit of Wax Apple [J]. Chinese Journal of tropical agriculture, 2019, 39(3): 75−79, 90.(in Chinese)
[21] 谢向誉, 陆柳英, 曾文丹, 等. 31份木薯种质资源的鉴定评价及遗传多样性分析 [J]. 南方农业学报, 2017, 48(3):393−400. XIE X Y, LU L Y, ZENG W D, et al. Identification, evaluation and genetic diversity analysis of 31 cassava germplasm resources [J]. Journal of Southern Agriculture, 2017, 48(3): 393−400.(in Chinese)
[22] 朱慧珺, 张耀文, 赵雪英, 等. 山西省绿豆种质资源的遗传多样性分析 [J]. 山西农业科学, 2019, 47(9):1540−1543, 1602. DOI: 10.3969/j.issn.1002-2481.2019.09.11 ZHU H J, ZHANG Y W, ZHAO X Y, et al. Genetic diversity analysis of mung bean germplasm resources in Shanxi Province [J]. Journal of Shanxi Agricultural Sciences, 2019, 47(9): 1540−1543, 1602.(in Chinese) DOI: 10.3969/j.issn.1002-2481.2019.09.11
[23] 马克平. 试论生物多样性的概念 [J]. 生物多样性, 1993, 1(1):20−22. DOI: 10.3321/j.issn:1005-0094.1993.01.006 MA K P. The concept of biological diversity [J]. Chinese Biodiversity, 1993, 1(1): 20−22.(in Chinese) DOI: 10.3321/j.issn:1005-0094.1993.01.006
[24] 时圣明, 潘明佳, 王洁, 等. 分子鉴定技术在中药中的应用 [J]. 中草药, 2016, 47(17):3121−3126. DOI: 10.7501/j.issn.0253-2670.2016.17.027 SHI S M, PAN M J, WANG J, et al. Application of molecular identification techniques in Chinese materia Medica [J]. Chinese Traditional and Herbal Drugs, 2016, 47(17): 3121−3126.(in Chinese) DOI: 10.7501/j.issn.0253-2670.2016.17.027
[25] 刘同金, 张晓辉, 沈镝, 等. 欧洲山芥种质资源的表型遗传多样性分析 [J]. 植物遗传资源学报, 2015, 16(3):528−534. LIU T J, ZHANG X H, SHEN D, et al. Analysis on genetic diversity of Barbarea vulgaris germplasm resources based on phenotypic traits [J]. Journal of Plant Genetic Resources, 2015, 16(3): 528−534.(in Chinese)
[26] 聂石辉, 彭琳, 王仙, 等. 鹰嘴豆种质资源农艺性状遗传多样性分析 [J]. 植物遗传资源学报, 2015, 16(1):64−70. NIE S H, PENG L, WANG X, et al. Genetic diversity of agronomic traits in chickpea(Cicer arietinum L.) germplasm resources [J]. Journal of Plant Genetic Resources, 2015, 16(1): 64−70.(in Chinese)
[27] 王晓阳, 黄丽芳, 闫林, 等. 基于SSR标记的我国咖啡种质资源遗传多样性分析及指纹图谱构建 [J]. 热带农业科学, 2018, 38(12):30−36, 42. WANG X Y, HUANG L F, YAN L, et al. Genetic diversity analysis and fingerprinting of coffee germplasm in China by SSR markers [J]. Chinese Journal of Tropical Agriculture, 2018, 38(12): 30−36, 42.(in Chinese)
[28] 曾潮武, 梁晓东, 李建疆. 新疆春小麦种质资源主要农艺性状的遗传多样性分析 [J]. 分子植物育种, 2017, 15(9):3740−3750. ZENG C W, LIANG X D, LI J J. Genetic diversity analysis in main characters of spring wheat germplasm in Xinjiang [J]. Molecular Plant Breeding, 2017, 15(9): 3740−3750.(in Chinese)
[29] 潘存祥, 许勇, 纪海波, 等. 西瓜种质资源表型多样性及聚类分析 [J]. 植物遗传资源学报, 2015, 16(1):59−63. PAN C X, XU Y, JI H B, et al. Phenotypic diversity and clustering analysis of watermelon germplasm [J]. Journal of Plant Genetic Resources, 2015, 16(1): 59−63.(in Chinese)
-
期刊类型引用(12)
1. 潘浩男,潘贞珍,环秀菊,覃潇敏,郑树芳,陆宇明,康专苗,陶亮,王文林,覃振师. 澳洲坚果桂热1号人工授粉有效性及花粉直感研究. 果树学报. 2025(01): 151-161 . 百度学术
2. 孔广红,陶丽,贺熙勇,倪书邦,陈丽兰,肖晓明. 不同授粉组合对澳洲坚果良种HAES863坐果及果实大小的影响. 中国果树. 2024(03): 93-97 . 百度学术
3. 孔广红,余文才,李玲,陶丽,贺熙勇. 云南西双版纳澳洲坚果开花与结实特性. 南方农业学报. 2024(04): 1089-1098 . 百度学术
4. 马静,李扬,陶亮,贺熙勇,龙青姨,宫丽丹. 澳洲坚果种质资源花序表型性状多样性分析. 云南农业大学学报(自然科学). 2024(04): 131-140 . 百度学术
5. 刘姚,张金云,吴世军. 广东省澳洲坚果种业发展现状与展望. 果树资源学报. 2023(01): 1-6 . 百度学术
6. 万继锋,卢智强,邹明宏,陈菁,宋喜梅,杨倩,董丽红,罗炼芳,曾辉. 澳洲坚果成熟期间果实内果皮颜色、果仁品质变化及最适采收期预测. 中国南方果树. 2023(06): 72-75+84 . 百度学术
7. 林玉虹,杨祥燕,蔡元保,潘如军,崔明勇,叶维雁,曾黎明,黄思婕,李冰. 68份澳洲坚果种质资源的表型分类研究. 南方农业学报. 2023(08): 2379-2387 . 百度学术
8. 万继锋,曾辉,邹明宏,陈菁,宋喜梅,杨倩,罗炼芳,陆超忠. 澳洲坚果种质资源果实表型多样性分析. 中国南方果树. 2022(06): 134-140 . 百度学术
9. 万继锋,曾辉,邹明宏,陈菁,罗炼芳,陆超忠. 八个澳洲坚果品种成年树抽梢、开花结果特性及栽培要点. 热带农业科学. 2021(03): 43-47 . 百度学术
10. 谭秋锦,韦媛荣,黄锡云,张涛,许鹏,宋海云,王文林,郑树芳. 10份澳洲坚果种质果实性状与营养成分分析. 果树学报. 2021(05): 672-680 . 百度学术
11. 万继锋,曾辉,邹明宏,陈菁,罗炼芳,陆超忠. 湛江地区澳洲坚果成年树抽梢生物学特性的观测分析. 中国南方果树. 2021(05): 65-70+75 . 百度学术
12. 项亚男. 澳洲坚果形态解剖学及产业研究进展. 南方农业. 2021(26): 151-152 . 百度学术
其他类型引用(0)