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. 引言
【研究意义】咖啡与可可、茶并称为三大饮料作物,主要分布在拉丁美洲、亚洲及非洲等热带地区,2017 年咖啡收获面积 1.09×103万hm2、产量921万t[1]。咖啡属Coffea约有100多个种,生产栽培的主要有小粒种(Coffea arabica L.)和中粒种(Coffea canephora Pierre),小粒种约占总产量 的60%、中粒种约占总产量 的40%[2]。我国引种咖啡已有100多年历史,主要分布在云南、海南、四川、台湾等地区,经过十多年的发展,目前已形成云南、四川小粒种咖啡优势产区和海南中粒种咖啡优势产区,2018年全国种植面积12万hm2,年产量17万t、农业产值20多亿元[3],是产区农民脱贫致富的主要经济来源。然而,我国咖啡产业存在种植品种单一、品质较差和比较效益偏低等突出问题,加强优良咖啡品种选育势在必行。因此,开展咖啡种质资源鉴定评价及遗传多样性研究,对咖啡优良品种的选育及咖啡产业发展均具有重要意义。【前人研究进展】我国收集保存了大量的咖啡种质资源并进行分析评价。周华等[4] 通过对其根、茎、叶、花、果等进行评价,明确了不同性状的表型分类。白学慧等[5]对96份咖啡种质抗锈病性状进行了评价。一些学者采用基于DNA水平的各种分子标记技术对我国咖啡种质资源进行了遗传多样性分析,明确其亲缘关系[6-10]。植物表型性状标记具有简单、明显和易于识别的特征,是进行种质资源鉴定和描述等研究最直接和最基本的方法和手段[11-13]。聚类分析可揭示种质资源个体间的遗传相似性,还能表明类群间的遗传差异和相互间的亲缘关系。主成分分析可揭示种质资源的生物学特点和生物学意义,能从多个因素指标中找出数量较少并能控制所有变量因子的主成分,其中每个主成分都能反映原始变量的大部分信息,且所含信息互不重复[13]。聚类分析和主成分分析可为作物种质资源鉴定评价和分类提供参考依据,已在豌豆[14]、陆地棉[15]、甜高粱[16-17]等作物进行了相关研究。【本研究切入点】表型性状鉴定是咖啡种质资源研究最基本的途径,咖啡因、绿原酸、蛋白质、脂肪及蔗糖等主要生化成分是咖啡风味及营养价值的主要决定因素[18]。虽然我国收集保存了大量的咖啡种质资源,但许多种质资源尚未充分鉴定评价与开发利用,因而加强咖啡种质资源鉴定评价是咖啡新品种选育研究的重要基础工作。目前未见应用聚类分析及主成分分析方法对咖啡种质资源主要农艺性状和品质性状进行鉴定和分析评价的报道。 【拟解决的关键问题】本研究以60份小粒种咖啡种质为试验材料,测量其8个农艺性状和5个品质性状,采用相关性分析、主成分分析和聚类分析相结合的方法,综合分析咖啡种质农艺性状和品质性状的遗传多样性,以期筛选优良种质或发掘咖啡资源的优异性状,为有效保护咖啡种质资源及品种改良提供参考。
1. 材料与方法
1.1 供试材料
供试材料为60份小粒种咖啡种质(表1),取自德宏热带农业科学研究所(农业部瑞丽咖啡种质资源圃),主要为波邦、铁皮卡、卡杜艾、卡蒂姆、卡杜拉、萨其姆等品种类型。
表 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 1.2 试验方法
1.2.1 农艺性状测定
鲜果重(kg):在果实盛熟期,采摘单株鲜果,并称取其重量,待全株果实采收结束时统计单株成熟鲜果重量。
干豆重(g):在果实盛熟期,采摘鲜果。去中果皮,洗掉种子表面果胶,置于50℃烘箱烘72 h,至带种皮豆含水量约为12%时,测得的数据即为干豆重。
种子大小(长、宽、厚):随机抽取20粒种子。用游标卡尺测量每粒种子基部至顶部的长度,即为种子长;种子最大横切面的最大直径,即为种子宽;种子最大横切面的最小直径,即为种子厚。
鲜干比(%):带种皮豆样品,去种皮制成标准商品豆,计算商品豆重占鲜果重的比例,重复3次,计算平均值。
百粒重(g):随机取100粒带种皮豆称重,重复3次,计算平均值。
出米率(%):随机称取(100±1)g带种皮豆样品,去种皮制成标准商品豆并称重,计算商品豆质量占带种皮豆质量的比例,重复3次,计算平均值。
1.2.2 品质性状测定
咖啡因按GB/T 5009.139-2014的方法[18]测定,绿原酸按GB/T 22250-2008的方法[19]测定,蔗糖按GB 5009.8-2016第二法[20]测定,蛋白质用杜马斯燃烧法测定。
脂肪测定采用项目组先前研制的“咖啡脂肪含量测定的超声超离方法”,将样品与石油醚混合液超声提取(60℃水浴、时间3 ~7 min、功率250 W),通过冷冻离心分离固相与液相(温度4℃、转速5 000~8 000 r·min−1、时间5~9 min),烘干固相物质,计算脂肪含量。
1.2.3 数据分析
采用 Microsoft Excel 2007 软件进行数据处理,计算平均值、标准差和变异系数,并参照谢向誉等[21]的方法计算遗传多样性指数(H′),即H′= −ΣPilnPi,Pi为性状第i级别内材料份数占总份数的百分比。采用SPSS 19.0软件进行相关性分析和主成分分析,利用 k-均值分类法进行聚类分析[22]。
2. 结果与分析
2.1 农艺性状表型评价
不同咖啡种质农艺性状变化较大,统计分析结果见表2。鲜果重变幅为0.79~2.97 kg,平均1.87 kg;干豆重变幅为300.49~592.68 g,平均390.58 g;鲜干比变幅为16.55%~48.79%,平均21.50%;种子长变幅为10.74~12.69 mm,平均11.74 mm;种子宽变幅为7.36~9.05 mm,平均7.98 mm;种子厚变幅为4.40~5.70 mm,平均 4.88 mm;百粒重变幅为15.06~27.13 g,平均19.51 g;出米率变幅为78.63%~86.60%,平均84.01%。
表 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 由表2可知,60份咖啡种质农艺性状存在很大差异,8个农艺性状的变异系数变幅为1.97%~25.15%,差异非常明显。变异系数的大小顺序为鲜干比>鲜果重>干豆重>百粒重>种子厚>种子宽>种子长>出米率。8个农艺性状中,变异系数以鲜干比最高(变异系数为25.15%),鲜干比最大值为48.79%、最小值为16.55%;其次是鲜果重(变异系数为19.85%),鲜果重最大值为2.97 kg、最小值为0.79 kg;以出米率变异系数最低(变异系数为1.97%),出米率最大值为86.80%、最小值为78.63%。鲜干比、鲜果重、干豆重变异系数大,说明这些性状存在着丰富的变异,具有较好的遗传改良基础。出米率、种子长、种子厚的变异系数较小,说明他们一致性较强,具有较稳定的遗传特性。
8个农艺性状的遗传多样性指数在1.43~2.04,其遗传多样性指数由高到低依次为种子长>种子厚>种子宽、出米率>干豆重>鲜果重、百粒重>鲜干比。种子长的遗传多样性指数最高,其后依次是种子厚、种子宽和出米率,说明种子大小具有丰富的遗传多样性,但由于他们的变异系数较小,基于种子大小对咖啡种质资源改良筛选潜力较小。鲜干比的遗传多样性指数最低,遗传多样性较窄。综上所述,60份咖啡种质资源的8个农艺性状存在较明显的差异,表现出不同程度的遗传多样性。
2.2 品质性状评价
不同咖啡种质品质性状变化较大,统计分析结果见表3。蛋白质变幅为11.98%~15.43%,平均13.40%;脂肪变幅为9.03%~15.91%,平均13.02%;绿原酸变幅为8.10~24.56 mg·g−1,平均16.59 mg·g−1;咖啡因变幅为7.62~11.80 mg·g−1,平均9.89 mg·g−1;蔗糖变幅为7.10~10.00 g·hg−1,平均8.48 g·hg−1。
表 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 60份咖啡种质品质性状存在很大差异,蛋白质、绿原酸等5个性状的变异系数为5.2%~22.68%,差异非常明显。在5个品质性状表现中,变异系数的大小顺序为绿原酸>脂肪>咖啡因>蔗糖>蛋白质。绿原酸变异系数最大,最高含量为24.56 mg·g−1、最低含量仅为8.10 mg·g−1、平均值为16.59 mg·g−1,说明绿原酸离散程度较大,稳定性较差,具有较好的遗传改良基础;脂肪的变异系数次之,为12.40%;咖啡因和蔗糖的变异系数分别为9.84%和9.42%;蛋白质变异系数最小,为5.20%。蛋白质、蔗糖、咖啡因的变异系数均较小,说明他们一致性较强,具有较稳定的遗传特性。
5个品质性状遗传多样性指数在0.55~2.07,其遗传多样性指数由高到低依次为绿原酸>脂肪>蛋白质>咖啡因>蔗糖。绿原酸遗传多样性指数最高(2.07),说明具有丰富的遗传多样性,而且其变异系数也最大(22.68%),表明基于绿原酸对咖啡种质资源筛选改良潜力较大。脂肪(2.06)、蛋白质(2.05)和咖啡因(1.99)的遗传多样性指数也较高,说明这些性状具有丰富的遗传多样性,但其变异系数较小,基于这些性状对咖啡种质资源进行改良筛选的潜力较有限。蔗糖的遗传多样性指数最低,说明供试咖啡种质在蔗糖方面的遗传多样性最小。综上所述,供试咖啡种质的主要品质性状具有丰富的遗传多样性。
2.3 农艺性状、品质性状的相关性分析
对咖啡农艺性状和品质性状进行统计和分析(表4),结果表明,鲜果重与绿原酸呈极显著负相关;干豆重与鲜干比、种子长、蛋白质呈极显著正相关,与出米率呈显著正相关,其中与种子长和鲜干比的相关系数分别为0.755和0.714,明显大于与其他性状间的相关系数,说明干豆重与种子长和鲜干比存在紧密关联;鲜干比与种子长呈极显著正相关;种子长与百粒重、出米率、蛋白质呈极显著正相关,与蔗糖呈显著正相关;种子厚与百粒重、出米率、蛋白质呈极显著正相关,与脂肪呈显著正相关;百粒重与出米率、蛋白质、脂肪呈极显著正相关;出米率与蛋白质、脂肪呈极显著正相关;蛋白质与脂肪呈极显著正相关;咖啡因与蔗糖呈极显著负相关。其余性状间相关不显著。
表 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.综上所述,鲜果重越大,绿原酸含量越低;干豆重越大,鲜干比和种子长越大,出米率和蛋白质含量也越高;种子越长越厚,百粒重、出米率也越大,蛋白质含量也越高;出米率越高,蛋白质、脂肪含量越高;蛋白质含量越高,脂肪含量越高;咖啡因含量越高,蔗糖含量越低。
2.4 主成分分析
8个农艺性状和5个品质性状的主成分分析结果(表5)显示,前4个主成分(PC1+ PC2+ PC3+ PC4)累计贡献率为69.17%。PC1贡献率为29.17%,种子长、种子宽、种子厚、百粒重和出米率在第1主成分有较高的载荷,特征值中贡献最大的是百粒重,故第1主成分为百粒重的构成因子;PC2贡献率为16.25%,鲜果重、干豆重、脂肪、绿原酸在第2主成分有较高的载荷,特征值中贡献最大的是鲜果重,故第2主成分为鲜果重的构成因子;PC3贡献率为12.35%,蛋白质、咖啡因在第3主成分有较高的载荷,其特征向量值中蛋白质贡献最大,故第3主成分为蛋白质构成因子;PC4贡献率为10.86%,鲜干比、绿原酸在第4主成分有较高的载荷,特征向量值中鲜干比贡献最大,故第4主成分为鲜干比构成因子。可见,前4个主成分的12个性状是构成小粒种咖啡种质多样性的主要因素,也是小粒种咖啡育种中应重点考虑的性状指标。
表 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 2.5 聚类分析
基于农艺性状、品质性状差异,对60份小粒种咖啡资源进行聚类分析,结果如图1所示。在遗传距离为25.0处,60份小粒种咖啡资源分为两个大类群,第I类有53份资源,第II类有7份资源,主要表现为鲜果重差异。第I类资源鲜果重均在2.20 kg以下,第II类资源鲜果重均在2.20 kg以上。可见,小粒种咖啡种质资源的鲜果重存在明显的差异。第I类又分为2个亚类,第I-1类包括49份资源,包含了卡蒂姆和萨其姆及未知品种类型的大多数资源,主要表现为比第I-2类有较大的鲜果重和种子宽,但绿原酸含量偏低,均在15.59 mg·g−1以下。第I-1类又分为a和b两个分支,a分支较b分支鲜果重大,绿原酸、咖啡因含量高。第I-2类有4份资源,分别是180、204、282和284,主要特点是鲜果重小,均在1.30 kg以下。第II类又分为两个亚类,第II-1类有3份资源,分别是 10、51和88;第II-2类有4份资源,分别是55、78、93和140。第II-2类较第II-1类具有更高的鲜果重和干豆重,鲜果重均在2.60 kg以上,干豆重均在538.00 g以上。综上所述,第II类鲜果重大,尤其是II-2 类,可作为选育和改良品种的基本材料。
3. 讨论与结论
遗传多样性是生物多样性的基础,生物遗传变异均发生在分子水平,在不断进化过程中形成了丰富的遗传多样性[23-24]。种质资源遗传差异大小直接影响性状遗传改良和研究利用的效果,优异基因能否被发掘与有效利用影响着新品种选育或新材料创制[25-26]。目前,我国云南和海南等地已开展对咖啡种质资源的收集保存,周华等[4]的研究结果表明了不同咖啡种质资源之间存在丰富的遗传多样性,闫林等[6-8]、黄丽芳等[9-10]、王晓阳等[27]的研究结果也表明小粒种咖啡遗传关系较紧密,但从农艺性状和品质性状方面进行种质资源遗传多样性分析还未见报道。变异系数反映了某性状数据的离散程度,其大小与性状的变异范围呈正相关,变异系数越大,表明性状的变异程度越大[28]。本研究对60份小粒种咖啡种质资源的8个农艺性状和5个品质性状进行分析比较,鲜干比、鲜果重、干豆重、绿原酸的离散程度大,稳定性差,具有较好的遗传改良基础;出米率和蛋白质的离散程度较小,性状表现较为稳定。遗传多样性信息指数是评价种质资源的主要指标之一,值越高,表明某性状的多样性越丰富[29]。本研究中平均遗传多样性指数1.82,其中种子长(2.04)、绿原酸(2.07)、脂肪(2.06)、蛋白质(2.05)遗传多样性指数较大,鲜干比(1.43)和蔗糖(0.55)较小。不同咖啡种质间的性状差异较大,表现出丰富的遗传多样性,且改良潜力大。
主成分分析是将多个变量转化为少数几个综合因子,每个主成分基本上又能反映原始变量的大部分信息,其所含信息互不重复。本研究的前4个主成分提供了原性状69.17%的信息,反映了所有表型性状的大部分相关信息。种子厚、百粒重、鲜果重、鲜干比、蛋白质对咖啡性状多样性构成起主要作用,可见咖啡种质资源的分类受到多个性状的共同影响。这些性状不仅可作为咖啡种质鉴定评价的重要性状,也可作为区分咖啡种质资源的主要性状。前4个主成分间包含的性状信息具有一定的相关性,且各性状类型差异较大,各主成分性状载荷值反映了主要性状的育种选择潜力。综上所述,在育种工作中,抓住种子厚、百粒重、鲜果重、鲜干比、蛋白质等主要因子,可望缩短选育新品种或新品系的年限。
农艺性状大多受多基因控制,不同性状间存在不同程度的相关性,某一性状的改变会导致其他相关性状发生变化,因此在育种过程中改良重要性状的同时可同步改良次要性状。本研究发现,种子长和种子厚与百粒重、出米率、蛋白质均呈极显著正相关,分别与蔗糖和脂肪呈显著正相关;出米率与蛋白质、脂肪呈极显著正相关;咖啡因与蔗糖呈极显著负相关。说明在咖啡品种选育过程中,通过筛选种子长和厚的种质作为亲本材料,可增加百粒重、出米率、蛋白质、蔗糖和脂肪等产量性状和品质性状;筛选咖啡因含量低的种质作为亲本材料,可增加咖啡蔗糖含量,提升品质。
本研究利用聚类分析法对60份咖啡种质的13个性状进行了聚类分析,让性状相近的聚为一类,使其同质化程度最高。每个组群具有特定的形态特征,以便更系统地分析咖啡种质资源的特点。第I类为单株鲜果重小的资源材料,这类材料具有较小的鲜果重和种子宽,绿原酸含量偏低等特点。第II类的鲜果重大,具有较大的鲜果重和种子宽,但绿原酸含量偏低。可见,小粒种咖啡种质资源的鲜果重存在明显的差异,尤其是II-2 类的4份种质,鲜果重和干豆重最大,可作为新品种选育和改良的基本材料。
60份咖啡种质资源存在丰富的遗传多样性,可通过引种驯化、杂交、嫁接和分子育种等技术把其优异基因延续下去。针对现有咖啡种质资源具有丰富的农艺性状和品质性状多样性等特征,今后应在现有研究的基础上,结合全基因组测序和代谢组学等方法,对各表型性状进行精准鉴定、关联分析,挖掘重要性状优异基因,并与常规育种方法相结合,加快我国咖啡种质资源的创新利用和新品种选育进程。
-
表 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)
-
期刊类型引用(4)
1. 何安乐,黄丽芳,胡丽松,庞永青,王晓阳,伍宝朵,董云萍,闫林,李学俊. 42份海南咖啡种质资源品质性状鉴定分析. 分子植物育种. 2024(03): 871-880 . 百度学术
2. 黄丽芳,龙宇宙,李金芹,董云萍,王晓阳,陈鹏,王宪文,闫林. 低温胁迫对小粒种咖啡幼苗生理特性的影响. 中国农业科技导报. 2023(02): 60-67 . 百度学术
3. 王灿,方先媛,张大英,李学俊,郑维全,李志刚,杨建峰,祖超,鱼欢. 柬埔寨主要胡椒产区胡椒品质比较分析. 热带农业科学. 2021(07): 84-88 . 百度学术
4. 石玉涛,郑淑琳,吴伟伟,罗盛财,陈夏婷,王飞权,李远华,张渤. 武夷名丛茶树种质资源茶多糖抗氧化活性分析. 福建农业学报. 2020(07): 801-810 . 本站查看
其他类型引用(0)